openSDUC 1.3.31

Open Stochastic Daily Unit Commitment of Thermal and ESS Units (openSDUC) - Version 1.3.31 - September 21, 2024

Open Stochastic Daily Unit Commitment of Thermal and ESS Units (openSDUC)

It determines the system operation for supplying the demand at minimum cost.

The openSDUC model presents a decision support system for defining the generation operation of a large-scale electric system.

The scope of the model corresponds to one day or one week divided in load levels: 01-01 00:00:00+01:00 to 01-21 23:00:00+01:00. The time division allows a flexible representation of the periods for evaluating the system operation. For example, by 84 periods of two hours or by the 168 hours of the week.

It considers stochastic short-term yearly uncertainties (scenarios) related to the system operation. The operation scenarios are associated with renewable energy sources and electricity demand.

The model formulates an optimization problem including binary and continuous operation decisions.

The unit commitment (UC) model is a based on a tight and compact formulation including operating reserves. It considers different energy storage systems (ESS), e.g., pumped-storage hydro, battery, etc.

The main results of the model can be structured in these topics:

• Operation: the output of different units and technologies (thermal, storage hydro, pumped-storage hydro, RES), RES curtailment

• Emissions: CO2

• Marginal: Short-Run Marginal Costs (SRMC)

A careful implementation has been done to avoid numerical problems by scaling parameters, variables and equations of the optimization problem allowing the model to be used for large-scale cases.

Installation

There are 2 ways to get all required packages under Windows. We recommend using the Python distribution Anaconda. If you don’t want to use it or already have an existing Python (version 3.8 | 3.9 recommended, 2.7 is supported as well) installation, you can also download the required packages by yourself.

Miniconda (recommended)

1. Miniconda. Choose the 64-bit installer if possible.

   1. During the installation procedure, keep both checkboxes “modify the PATH” and “register Python” selected! If only higher Python versions are available, you can switch to a specific Python Version by typing conda install python=<version>

   2. Remark: if Anaconda or Miniconda was installed previously, please check that python is registered in the environment variables.

2. Packages and Solver:

   1. Launch a new command prompt (Windows: Win+R, type “cmd”, Enter)

   2. Install CBC solver via Conda by conda install -c conda-forge coincbc. If you have any problem about the installation, you can also follow the steps that are shown in this link.

   3. Install openSDUC via pip by pip install openSDUC

Solvers

HiGHS

The HiGHS solver can also be used. It can be installed using: pip install highspy. This solver is activated by calling the openTEPES model with the solver name ‘appsi_highs’.

Gurobi

Another recommendation is the use of Gurobi solver. However, it is commercial solver but most powerful than open-source solvers for large-scale problems. As a commercial solver it needs a license that is free of charge for academic usage by signing up in Gurobi webpage. It can be installed using: conda install -c gurobi gurobi and then ask for an academic or commercial license. Activate the license in your computer using the grbgetkey command (you need to be in a university domain if you are installing an academic license).

GLPK

As an easy option for installation, we have the free and open-source GLPK solver. However, it takes too much time for large-scale problems. It can be installed using: conda install -c conda-forge glpk.

CBC

The CBC solver is our recommendation if you want a free and open-source solver. For Windows users, the effective way to install the CBC solver is downloading the binaries from this link, copy and paste the cbc.exe file to the PATH that is the “bin” directory of the Anaconda or Miniconda environment. It can be installed using: conda install -c conda-forge coincbc.

Mosek

Another alternative is the Mosek solver. Note that it is a commercial solver and you need a license for it. Mosek is a good alternative to deal with QPs, SOCPs, and SDPs problems. You only need to use conda install -c mosek mosek for installation and request a license (academic or commercial). To request the academic one, you can request here. Moreover, Mosek brings a license guide. But if you are request an academic license, you will receive the license by email, and you only need to locate it in the following path C:\Users\(your user)\mosek in your computer.

Users

If you are not planning on developing, please follows the instructions of the Installation.

Once installation is complete, openSDUC can be executed in a test mode by using a command prompt. In the directory of your choice, open and execute the openSDUC_run.py script by using the following on the command prompt (Windows) or Terminal (Linux). (Depending on what your standard python version is, you might need to call python3 instead of python.):

openSDUC_main

Then, four parameters (case, dir, solver, and console log) will be asked for.

Remark: at this step only press enter for each input and openSDUC will be executed with the default parameters.

After this in a directory of your choice, make a copy of the 9n or sSEP case to create a new case of your choice but using the current format of the CSV files. A proper execution by openSDUC_Main can be made by introducing the new case and the directory of your choice. Note that the solver is glpk by default, but it can be changed by other solvers that pyomo supports (e.g., gurobi, mosek).

Then, the results should be written in the folder who is called with the case name. The results contain plots and summary spreadsheets for multiple optimised energy scenarios, periods and load levels as well as the investment decisions.

Note that there is an alternative way to run the model by creating a new script script.py, and write the following:

from openSDUC.openSDUC import openSDUC_run

openSDUC_run(<case>, <dir>, <solver>)