Ecosyste.ms: Awesome

An open API service indexing awesome lists of open source software.

Awesome Lists | Featured Topics | Projects

https://github.com/SPF-OST/optihood

The package provides a complete python-based framework to OPTImize the investment in alternative energy technologies (heat pump, solar thermal, PV, etc.) as well as the operation of available energy resources (natural gas, sun, grid electricity, etc.) at a neighbourHOOD-scale
https://github.com/SPF-OST/optihood

Last synced: 3 months ago
JSON representation

The package provides a complete python-based framework to OPTImize the investment in alternative energy technologies (heat pump, solar thermal, PV, etc.) as well as the operation of available energy resources (natural gas, sun, grid electricity, etc.) at a neighbourHOOD-scale

Awesome Lists containing this project

README 

        
# optihood: energy modelling and optimization framework



optihood provides a complete python-based framework to **OPTI**mize the investment in alternative energy

technologies (heat pump, solar thermal, PV, etc.) as well as the operation of available energy resources (natural gas,

sun, grid electricity, etc.) for decentralized energy networks on a neighbour**HOOD**-scale. It is designed to facilitate 

the researchers and energy planners in optimizing an energy neighbourhood in terms of cost and/or environmental criteria. 

It enables the users to perform both single-objective and multi-objective optimization and analysis. The energy model and 

it's associated parameters can be defined easily using an excel file. Additionally, a variety of plotting methods are 

defined for easy and fast result visualization.



## Documentation

The documentation can be found on https://optihood.readthedocs.io/.



## Prerequisites

In order to use optihood, the following prerequisites are needed on your machine:



- [Python >=3.9](https://www.python.org/downloads/) is installed. Python 3.12 is recommended.

- Git is installed

- An active [Github](https://github.com/) account to clone the repo.

- A solver is installed. [Gurobi solver](https://www.gurobi.com/resource/parallelism-linear-mixed-integer-programming/)

is recommended, although other solvers like CBC, GLPK, Cplex could also be used.



## Installation



As of now, optihood is available as an open source code and needs to be installed from source. Please follow the

instructions mentioned below to complete the installation. The commands given below are suited for the Windows platform

and should be run from within the optihood directory in a Command Prompt. For other platforms, similar alternative

commands could be used.



1. Clone the optihood repo to a folder called 'optihood' on your local machine:

    ```commandline

    git clone https://github.com/SPF-OST/OptiHood.git

    ```

2. All the next commands should be run from within the optihood folder. Create a virtual environment and activate it:

    ```commandline

    py -3.12 -m venv venv

    venv\Scripts\activate

    ```    

   If using python version other than 3.12, replace -3.12 with -3.X (3.X being the python version).



3. Install the requirements into the created virtual environment. The requirements.txt file needs Python 3.12. If using another version, use pip tools to compile the provided requirements.in file before installing the requirements:

    ```commandline

    pip install wheel

    pip install -r requirements.txt

    ```

    It might be required to install C++ build tools. To do that, click on the link that appears with the error message and follow the instructions (it is the lapack package that is missing). In order to be able to install the missing package, it is required to have a complete Visual Studio instance and installing it with the "Desktop development with C++" workload.



4. Verify the installation of the oemof package and the solver by running the installation test in your virtual environment:

    ```commandline

    oemof_installation_test

    ```



    

    If the installation is successful a message similar to the following would display (the installed solver would be marked as working):

    ```commandline

    *****************************

    Solver installed with oemof

    cbc: not working

    glpk: not working

    gurobi: working

    cplex: not working

    

    *****************************

    oemof successfully installed.

    *****************************

    ```

5. To test whether the installation worked well, you could run a [basic example](https://github.com/SPF-OST/OptiHood/tree/main/data/examples/).



6. Run the test suite using:

```commandline

    pytest

```



## Defining an energy network



Optihood offers a several functionalities to define an energy network, optimize it and visualize the results, which

provides the user with a complete framework for optimization without the need to code by hand. 



An energy network can be defined as an object of the EnergyNetworkIndiv class or the EnergyNetworkGroup class. This

object then acts as the primary container for the model. The EnergyNetworkClass is the main parent class, from which 

two child classes EnergyNetworkIndiv and EnergyNetworkGroup are inherited. Either EnergyNetworkIndiv class or 

EnergyNetworkGroup class could be used to define an energy network. The choice mainly depends on whether the buildings 

are linked together (electrically and/or thermally) or not.



If the buildings within an energy network do not share electricity and/or heat, EnergyNetworkIndiv class is used:



    import EnergyNetworkIndiv

    network = EnergyNetworkIndiv(dateTimeIndex, tSH=35, tDHW=55)

    

Otherwise, if the buildings are expected to share energy (electrical and/or heat), EnergyNetworkGroup class is used:



    import EnergyNetworkGroup

    network = EnergyNetworkGroup(dateTimeIndex, tSH=35, tDHW=55)

    

The first parameter to be passed in both the cases is a Datetime index. This parameter gives the time range for an

optimization model. The Datetime index could be defined using date_range() in pandas:



    import pandas as pd

    dateTimeIndex = pd.date_range('2021-01-01 00:00:00', '2021-12-31 23:00:00', freq="60min")

    

The second and the third parameters tSH and tDHW define the temperatures for space heating and domestic hot water

production, respectively.



Once the 'network' object has been created, the next step then is to build the model from an input excel file which

defines different components which constitute the model, how they are connected and their associated parameters:



    network.setFromExcel(inputExcelFilePath, numberOfBuildings, clusterSize, opt)

    

'inputExcelFilePath' gives the path of the excel input file. 'numberofBuildings' is an integer parameter specifying the

number of buildings defined in the excel file. The last two parameters clusterSize and opt are optional. The 'opt'

parameter could be either 'costs' (default value) or 'env' depending on which criteria should be optimized. The

'clusterSize' parameter is used to provide a selected number of days which could be assumed representative of the entire

time range. For example: two typical days could be selected to model the entire year, which could represent two clusters

summer and winter. This would improve the optimization speed. If not given during the function call, the default value

of the clusterSize parameter assumes no day clusters.



## Input Excel File



The input excel file is used to define an optimization model and set the model parameters. Each sheet of this excel file

is structured to defin different components, such as buses, storages and transformers, their respective parameters,

connections between these components and the building to which they belong.



The input excel file typically has 9-10 sheets, each defining a different component type of the model.



### buses

This excel sheet defines the buses used in the energy network. Buses define the connections between different

components. Each row of this excel sheet represents a bus node in the model.



- **label** (string):

    label name of the bus. Should be unique for each building i.e. different buildings could have the same label for two

    buses but for one building the buses should not have duplicate labels.



- **active** (0 or 1):

    If set to 1 then the corresponding bus is active, otherwise (if 0) then the bus is inactive in the model. Could be

    used to optimize different scenarios with different buses without the need to re-construct the base excel file.



- **excess** (0 or 1):

    If set to 1 then the corresponding bus has the possibility of having an excess flow unbalanced to the demand. An

    additional sink node is automatically added in this case to accept this excess flow. If set to 0 then the bus cannot

    have an excess flow.



- **excess costs** (float):

    The associated cost of excess flow from a bus. Relevant only if **excess** is set to 1. Could be negative to denote

    monetary gains such as in case of PV feed-in.



- **building** (integer):

    Building number to which the bus belongs. Should be unique for each building.



### commodity_sources

This sheet defines the different commodity sources which serve as an energy input to the model. The parameters **label**,

**active** and **building** are analogous to the parameters described earlier for buses.



- **to** (string):

    Label of bus to which the energy from the commodity source flows. The corresponding bus label should exist in

    the buses sheet.



- **variable costs** (float):

    Cost per kW of the commodity source.



- **CO2 impact** (float):

    CO2 impact per kW of the commodity source.



### demand

The nodes related to the energy demand i.e. sink are defined in this sheet. The parameters **label**, **active** and

**building** are analogous to the parameters described earlier for buses.



- **from** (string):

    Label of the bus from which the energy flows to the demand node. The corresponding bus label should exist in

    the buses sheet.



- **fixed** (0 or 1):

    If set to 1, the energy demand profiles are fixed and a path to the demand profiles should then be given in the

    csv_data sheet of the input excel file. If set to 0, the optimizer will select the demand profiles for optimum

    operation (a total annual demand profile could be given in this case). At the present development stage only fixed

    demands are accepted.



- **nominal value** (float):

    The demand profiles series should be normalized values. this parameter then defines the maximum demand with which

    the normalized series is multiplied to obtain the actual demand profile series. If set to 1, then the given demand

    profile series is taken as it is (not normalized).



### transformers

The nodes related to the energy conversion units (or transformers) such as CHP, heat pump, etc. are given in this excel

sheet. The parameters **label**, **active** and **building** are analogous to the parameters described earlier for buses.



- **from** (string):

    Label of bus from which the energy flows to the transformer node. The corresponding bus label should exist in

    the buses sheet.



- **to** (strings separated by comma):

    Label of the bus(es) to which the energy flows from the transformer node. Different bus labels should be separated

    by a comma (,), for example: Bus1, Bus2, Bus3. The bus label(s) should exist in the buses sheet.



- **efficiency** (float values separated by comma):

    Conversion efficiency from input to the output(s) of the transformer node. Efficiencies should be separated by

    a comma (,) in the case with more than one outputs (i.e. when more than one buses are listed under **to**).



- **capacity_DHW** (float):

    Maximum capacity limit for domestic hot water (DHW) production from the transformer unit in kW. Could be left blank

    if the parameter is not relevant to a transformer (for example if a transformer does not produce DHW)



- **capacity_SH** (float):

    Maximum capacity limit for space heat (SH) production from the transformer unit in kW. Could be left blank if the

    parameter is not relevant to a transformer (for example if a transformer does not produce SH)



- **capacity_el** (float):

    Maximum capacity limit for electricity production from the transformer unit in kW. Could be left blank if the

    parameter is not relevant to a transformer (for example if a transformer does not produce electricity). Note for CHP

    this parameter acts as the main capacity against which optimization is performed.



- **capacity_min** (float):

    Minimum capacity to be installed in kW for an investment in the transformer unit.



- **lifetime** (non-negative integer):

    Lifetime of the node in years.



- **invest_base** (float):

    Base investment cost of the node.



- **invest_cap** (float):

    Investment cost per unit installed capacity (i.e. per kW) of the node.



- **maintenance** (float):

    Maintenance cost of the node. Given as a percentage of the base investment cost **invest_base**. For example:

    if the investment cost is to be taken as 5% of **invest_base**, then the value of **maintenance** field should be 0.05.

    Set as 0 if this cost is to be ignored.



- **installation** (float):

    Installation cost of the node. Given as a percentage of the base investment cost **invest_base**. For example:

    if the Installation cost is to be taken as 15% of **invest_base**, then **installation** field should be 0.15. Set

    as 0 if this cost is to be ignored.



- **planification** (float):

    Cost associated with planning. Given as a percentage of the base investment cost **invest_base**. For example:

    if the planning cost is 5% of **invest_base**, then **planification** should be 0.05. Set as 0 if this cost is to be

    ignored.



- **heat_impact** (float):

    Environmental impact for heat production. Set as 0 if there is no impact due to heat production.



- **elec_impact** (float):

    Environmental impact for electricity production. Set as 0 if there is no impact due to electricity production.



- **impact_cap** (float):

    Environmental impact per unit installed capacity i.e. per kW of the node. Set as 0 if there is no impact per

    installed capacity of the node.



### solar

This excel sheet defines the parameters related to the solar components such as solar thermal collector and PV panels.

The parameters **label**, **active** and **building** are analogous to the parameters described earlier for buses.

**from** and **to** parameters have been previously defined for commodity sources and demand sheets, respectively, while

the cost and environmental impact paramaters are described under transformers sheet.



- **connect** (string):

    Label of the bus which connects a solar collector to the model. This bus allows excess heat production from the solar

    collector. A node for heat sink is created automatically. The given bus label should exist in the buses sheet. This

    parameter is irrelevant for PV.



- **electrical_consumption** (float):

    Electrical consumption of the solar component. Given as a percentage fraction of the energy produced. 0.02 means

    the electrical consumption is 2% of the energy is produced.



- **peripheral_losses** (float):

    Peripheral losses of the solar component. Given as a percentage fraction of the energy produced. 0.05 means 5% of

    the energy produced is lost to the surrounding environment.



- **latitude** (float):

    Latitude of the geographical location where the solar collector/panel is placed. Given in degrees.



- **longitude** (float):

    Longitude of the geographical location where the solar collector/panel is placed. Given in degrees.



- **tilt** (float):

    Tilt angle of the solar collector/panel. Given in degrees.



- **azimuth** (float):

    Azimuth angle of the solar collector/panel. Given in degrees.



- **eta_0**, **a_1** and **a_2** (float):

    Efficiency parameters of the solar thermal collector. Solar thermal collector is linearized using the pre-calculations

    given in [oemof-thermal](https://oemof-thermal.readthedocs.io/en/latest/solar_thermal_collector.html).



- **temp_collector_inlet** (float):

    Inlet fluid temperature of the solar thermal collector. Given in degree C.



- **delta_temp_n** (float):

    Temperature difference between the inlet fluid and the mean fluid temperature in case of solar collector. For PV,

    this parameter denotes the temperature difference between the solar cells and the ambient.



- **capacity_max** (float):

    Maximum capacity limit in kW.



- **capacity_min** (float):

    Minimum possible capacity in kW for the installation of solar collector/panel.



### storages

This excel sheet defines the parameters related to the energy storage units such as battery and hot water tank. **label**,

**active** and **building** have been defined previously for buses excel sheet. A description of **from** and **to** has

been given in commodity sources and demand sheets, respectively. The cost and environmental impact parameters are

described in the transformers sheet section. **capacity_min** and **capacity_max** are described in the solar excel sheet

section.



- **efficiency inflow** (float):

    Charging efficiency of battery. This parameter is not relevant for thermal storages.



- **efficiency outflow** (float):

    Discharging efficiency of battery. This parameter is not relevant for thermal storages.



- **initial capacity** (float):

    initial capacity of the storage. This parameter is expressed as a fraction of the total storage capacity. 0 means storage is initially 

    assumed to be empty, 1 denotes that the storage is 100% full initially, while 0.5 means the storage is at 50% capacity initially.



- **capacity loss** (float):

    Losses from battery storage. This parameter is not relevant for thermal storages.



### stratified_storage

This excel sheet defines the parameters relevant to stratified thermal storage. The pre-calculations given in [oemof-thermal](https://oemof-thermal.readthedocs.io/en/latest/stratified_thermal_storage.html)

have been used to linearize the thermal hot water storage. The parameter names used here are similar to the parameters

defined in oemof-thermal.



### links

This excel sheet defines the parameters for electricity and space heating links. The buildings could share electricity

production and/or space heat production. Links allow this sharing to be possible. **label** and **active** have been

defined already for buses excel sheet. **invest_base** and **invest_cap** parameters (defined in the transformers sheet

section) are only relevant for space heating links in the present stage of development.



- **buildingA** (integer):

    Building number of the first building of the link. This should match with the values typically given in the

    **building** parameter in the other excel sheets.



- **buildingB** (integer):

    Building number of the second building of the link. This should match with the values typically given in the

    **building** parameter in the other excel sheets.



- **efficiency from A to B** (integer):

    Efficiency of energy transfer over the link from **buildingA** to **buildingB**.



- **efficiency from B to A** (integer):

    Efficiency of energy transfer over the link from **buildingB** to **buildingA**.



### csv_data

The paths to CSV files containing demand profiles, weather data and electricity impact data are to be given in this

excel sheet. **INFO** gives further information about each row.



### grid_connection

This excel sheet should not be modified by the users. It defines the separation of the flows from electricity grid and

the produced electricity flows to make sure that the grid electricity is not stored in batteries.



## Collaborators



SPF/OST & HEIG-VD