Next Article in Journal
Remote Monitoring of Joints Status on In-Service High-Voltage Overhead Lines
Previous Article in Journal
Effect of LaNi3 Amorphous Alloy Nanopowders on the Performance and Hydrogen Storage Properties of MgH2
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Different Interval Lengths in a Rolling Horizon MILP Unit Commitment with Non-Linear Control Model for a Small Energy System

Institute of Energy Systems, Hamburg University of Technology, Denickestr. 15, 21073 Hamburg, Germany
*
Author to whom correspondence should be addressed.
Energies 2019, 12(6), 1003; https://doi.org/10.3390/en12061003
Submission received: 15 February 2019 / Revised: 8 March 2019 / Accepted: 11 March 2019 / Published: 14 March 2019
(This article belongs to the Section F: Electrical Engineering)

Abstract

:
In this paper, a fixed electricity producer park of both a short- and long-term renewable energy storage (e.g., battery, power to gas to power) and a conventional power plant is combined with an increasing amount of installed volatile renewable power. For the sake of simplicity, the grid is designed as a single copper plate with island restrictions and constant demand of 1000 MW; the volatile input is deducted from scaled 15-min input data of German grid operators. A mixed integer linear programming model is implemented to generate an optimised unit commitment (UCO) for various scenarios and configurations using CPLEX® as the problem solver. The resulting unit commitment is input into a non-linear control model (NLC), which tries to match the plan of the UCO as closely as possible. Using the approach of a rolling horizon the result of the NLC is fed back to the interval of the next optimisation run. The problem’s objective is set to minimise CO2 emissions of the whole electricity producer park. Different interval lengths are tested with perfect foresight. The results gained with different interval lengths are compared to each other and to a simple heuristic approach. As non-linear control model a characteristic line model is used. The results show that the influence of the interval length is rather small, which leads to the conclusion that realistic forecast lengths of two days can be used to achieve not only a sufficient quality of solutions, but shorter computational times as well.

1. Introduction

The expansion of volatile renewable energy (VRE) is a major project, that many countries decided to undergo as part of their commitments under the UN Paris Climate Change Conference of 2016 [1]. Various energy system configurations (i.e., distribution of technologies, their installed rated power and capacities) are suggested to meet the goals of e.g., the German Energiewende [2]. Energy systems can consist of VRE plants (e.g., wind, photovoltaic, concentrating solar power), of renewable storage facilities (e.g., batteries, power to gas to power, pumped hydro) and of conventional storage facilities or “primary energy storages” [3] (e.g., combined cycle gas turbine (CCGT), coal-fired steam power plant). To avoid confusion, in this paper the term “energy storage” is used solely for renewable storage facilities or “secondary energy storages” [3]. “Primary energy storages” are referred to as conventional power plants throughout this paper.
For this work a small energy system without transmission losses (often referred to as “copper plate”) and no possible external source or sink of power (often referred to as “island restriction”), in the following referred to as “copper plate island”, is the test system. The copper plate island consists of a CCGT power plant, a short-term and a long-term energy storage and varying instalments of solar and wind power. Installed solar and wind power are equally distributed across all configurations and their installed power is given as VRE power in MWeach, meaning the installed power each generation type has in a given configuration. As demand for the energy system, a constant power of 1000 MW is defined.
To utilise given energy system configurations to their maximum potential as well as test their feasibility, time series simulation is a commonly used tool. These simulations need some kind of unit commitment (UC) for dispatching the plants in the given energy system.
The methods for UC can vary from heuristic to non-linear optimisation problems, from deterministic to stochastic with all their methodical benefits and drawbacks [4,5]. A popular approach is the formulation of the energy system as a mixed integer linear problem (MILP) for small [6] to larger systems [7,8]. Most often the optimisation problem is formulated to find the best solution for the entire time series (e.g., one year) at once, which can become very complex and time consuming and might make local computing resources incompatible [9]. To solve this problem, the use of a rolling horizon became more and more popular in recent years [5,10,11,12].
In [13] a comprehensive review of energy system modelling tools can be found. Seventy-five models were categorised in purpose, approach, methodology, temporal resolution, modelling horizon, spatial resolution and coverage and availability. These models were then discussed and their features listed. More than half of the presented models use a linear programming or MILP approach. The use of rolling horizons is however rarely mentioned.
Recently developed frameworks such as the one of [11] demonstrate the significance of the rolling horizon approach. This framework was developed to show the effect of different interval lengths (e.g., five hours or four days) on a one-day unit commitment problem. It would allow for modification of the MILP results in order to adjust for uncertainties or new information, but was not used in [10]. The framework was updated in [14] to model start-up times and tested for different rolling horizon interval lengths between two and 24 h. [15] presented a different framework, where the rolling horizon was implemented for a similar reason, but added an optional economic dispatch algorithm to adjust the suggested UC. This is to model market behaviour more properly.
This use of cascading UC to mimic the behaviour of a day-ahead and intra-day market was implemented by [16]. The market behaviour was modelled by a 24 h MILP unit commitment with a stochastically generated VRE forecast for the day ahead and a 4 h non-linear point algorithm to minimise the cost associated with diverging from the day-ahead UC.
A similar approach of stochastic programming for the UC was investigated and tested against a deterministic approach by [17]. The stochastic problem was tested in various configurations of rolling horizon lengths for the problem and solution aggregation on a larger energy system with electric grid constraints. In a similar fashion this research topic was investigated by [18], with the difference that in [18] various lengths for the final unit commitment—minimising the divergence from the stochastic forecast—were tested. Another difference lies in the configuration of the underlying energy system, whereby both [17,18] use different kinds of short-term energy storage devices and no long-term energy storage devices.
Investigation of lengths for the rolling horizon was done by [19] to find the optimal use of power plant flexibility considering the start-up and maintenance related problems of CCGT power plants.
Focusing on a thermal energy storage device for the medium term, [20] investigated a rolling horizon with a week-long forecast against a previous research with a day-ahead forecast and found the results to be better for the week-long unit commitment.
To adjust for non-linear behaviour of the system [20,21] did a piecewise linearisation of the price-effects that influenced the objective function. The focus of [21] was the usage of upper-, lower-bound or mean value of the linearised piece from the price function.
This paper focusses on energy system feasibility for future energy systems, for which [22] reviewed various energy system simulation models. The difference to the previously discussed systems is that future energy systems contain long-term energy storage devices. In this case, the UC might need longer forecasting lengths, as their name implies. It should be noted however, that [16,20,21] still adjust their models for non-linearities despite their rather short interval length.
For a better integration of long-term storages [23] integrates the rolling horizon approach to apply different levels of detail. For the exact UC, a complex model is used with short horizons, which is thought to not encompass the characteristics of long-term storages. A less complex model is used for the entire time frame, to set the characteristics for the exact planning.
Similarly [24] applied the use of a less complex model, but for a whole year simulation. The research investigated the use of an evolutionary algorithm for the energy system configuration with an underlying two level MILP UC. For the lower level a one-day rolling horizon MILP UC was implemented. The parameters of the MILP were generated from a non-linear model. Unfortunately the model of [24] does not ensure, that the long-term storage is filled to its initial state of charge (SOC) at the end of the simulation.
The rolling horizon is done here primarily to save simulation time, as [25] proposed, where the effect of different lengths of the rolling horizon and their effect on simulation time und quality of the solution was investigated. The quality of the solution was judged against the solution of a MILP UC solution for the entire time frame at once.
This paper introduces a different modelling approach for the solution of UC problems, in which the solution of a MILP UC is handed to a non-linear control model (NLC). That way the feasibility of the proposed UC can be tested. Further the MILP formulation uses the approach of [26] to model part-load behaviour more correctly, where the parameters of the MILP and initial values in a rolling horizon operation are generated from the NLC. With this approach, the effect of the interval length (see Section 2.1) in a rolling horizon problem becomes the research focus of this paper. To our knowledge, this has not been addressed with an NLC before.
To determine the effect of the interval lengths the results of MILP unit commitment optimisations (UCO) with varying temporal settings are compared to a unit commitment with a heuristic approach (UCH). The CO2 emissions of the electricity producer park are the minimisation objective of the MILP, which is why the specific CO2 emissions of the electricity producer park are used as the primary quality measure of the UC results. Additionally, the share of storages, meaning the amount of energy delivered by the energy storage devices in respect to total demand, acts as an indicator of difference between UCO and UCH.
A characteristic line model, namely CharL (see Supplementary Materials), is inherent to the UCH and used as an NLC for the UCO, instead of simply aggregating the MILP results (see Section 2.1).

2. Model Structure and Materials

In this work, the used model’s base functionality is a model of characteristic lines, such as the efficiency curve of a thermal power plant. Figure 1 exemplarily shows characteristic lines of the relative efficiency over the relative power for combined and open cycle gas turbines according to [27]. For application in the model, those curves are scaled with net efficiencies and used to compute the use of fuel and emissions of CO2 at any given point of time in the simulation.

2.1. Program Flow of CharL

The CharL library can be used to calculate time series of given energy system configurations. It is possible to simulate on a UCH or by incorporation of a UCO generated operation plan as the simplified program flow in Figure 2 shows. The displayed scheme starts after the initial data reading and initialisation and begins with assessing whether the unit commitment is done heuristically or via MILP optimisation. In the latter case necessary parameters for the MILP model (see Section 2.6.3) are optimised. Both paths then start the modelling process for all time steps t . Every time step starts with the assessment of VRE power. If the time step is at the start of a new optimisation period and the unit commitment is not done heuristically, a MILP is formulated after reading the VRE power forecast for the current interval and technological constraints (e.g., storage level) from the previous time step. Afterwards the MILP is handed to the commercial solver CPLEX [28], whose final solution is saved (see Figure 2) for all timesteps of the period k (see Figure 3). After the optimisation is finished or when the unit commitment is done heuristically, the program continues at the time step t and operates first the storage units and then the conventional power plant units as close as possible to the proposed unit commitment. The proposed unit commitment is followed for next timesteps until the beginning of a new period starts, where a new MILP is formulated and solved. The level for the storage units and the emissions are calculated after the proposed powers of the respective units were deemed valid. The results are evaluated after finishing the last time step. The program is held entirely in C++ with the external libraries FLOPC++ [29] and OSI for the formulation of the MILP and the Qt library solely for the GUI (see Supplementary Materials).
As previously stated, the formalisation of the linear problem is repeated at every new beginning of a period but formulated for the length of an interval. The term “interval” hereby refers to the time period that the UCO problem is formulated for, while the term “period” means the time period for which the result of the UCO is saved and input into the NLC. The interval length therefore is always equal or greater than the period length. Greater interval lengths lead to an overhang that is (re-)solved in the following period(s). The overhang prevents consecutive operation plans from being independent of each other. This process can be seen in Figure 3 and is called a rolling horizon.

2.2. Renewable Energy Modelling

For complexity and computational resources the renewable energy profile was modelled using scaled solar and wind production of 2015 as made publicly available by the German transmission system operators (TSO) [30,31,32,33,34,35,36,37]. The scaling was done using the monthly installed capacities of solar and wind, which are provided by the TSOs [38]. The resulting production profiles for wind and solar of 2015 can be seen in Figure 4. Lastly, for the greater interval length and therefore bigger optimisation problems, these curves are averaged to result in hourly temporal resolution.

2.3. Modelling Residual Power Generators

The residual power generators are the (renewable) energy storage devices’ discharge units and the conventional power plants, which could be considered a discharge unit of a conventional energy storage device (see Section 1). This results in a similar modelling approach, namely using the efficiency curve to calculate the needed power of the input form of energy.
Since energy storage devices can only discharge, if there is energy already stored (i.e., after a charging process), the charging process has to be modelled using a second efficiency curve. Additionally, the energy storage devices have a self-discharge rate when not in use, as [39] defined self-discharge only under open-circuit conditions. Accordingly, a fixed round-trip efficiency, meaning the efficiency of power to storage to power, is not used in this work.
Furthermore, maximum load gradients are considered for all units. The emissions of power plants are calculated using their fuel’s respective specific emission per MWhth as an additional parameter. The NLC considers start-up times as the corresponding functions are implemented. For this paper, the residual power generators are parameterised so that the start-up times have no influence on the UC. That way the MILP model (see Section 2.5) can resemble the NLC more closely.

2.4. Heuristic Unit Commitment

The UCH follows a very simple approach to integrate any excess energy of VRE plants into energy storage devices and meet any resulting residual load from the energy storage devices, if possible. Any further resulting residual load will be matched by the conventional thermal power plants.
The UCH iteratively solves the entire time series until the storage levels of all storages at the initial and final timestep do not change any further from one run to the next within a tolerance of 1 MWh.

2.5. MILP Formulation

The MILP formulation is designed to mirror the modelling approach of CharL as close as possible. The following equations are written as they are programmed using the FLOPC++-library [29] to address the solver CPLEX [28]. Using Gurobi [40] as the solver would have been possible as well, but was disregarded due to practical reasons.

2.5.1. Constraints to Model Thermal Power Plants

Thermal power plants are modelled with a semi-continuous operating range. A binary variable ψ states, whether the unit is online or offline, while P min and P max are minimum and maximum limits for the power generation P , which is fed into the grid. For every period k K the constraints need to be formulated for all generating plants g G and time steps t T k :
P ( g , t ) P min ( g ) ψ ( g , t ) ; g G t T k . P ( g , t ) P max ( g ) ψ ( g , t ) ; g G t T k .
In addition to the operating range of the units, ramping limits R are applied using the following equations considering the length of a time step τ . It is possible to differentiate between the maximum upward and downward gradient:
P ( g , t 1 ) P ( g , t ) R down ( g ) τ ; g G t T k . P ( g , t 1 ) P ( g , t ) R up ( g ) τ ; g G t T k .
For t = T k , 0 the variable for t 1 is replaced by data from the non-linear control of the previous period.
The efficiency of a power generation unit usually shows non-linear characteristics, which cannot be integrated directly into a MILP formulation. In accordance with [26], the efficiency characteristic is implemented using the load-unspecific basic demand. Hereby, the power on the unit side (and therefore the fuel demand) is described with P u :
P ( g , t ) = a ( g ) P u ( g , t ) a ( g ) b ( g ) P nom ( g ) ψ ( g , t ) ; g G t T k .
Hereby a constant loss factor or basic demand, which is calculated in relation to the nominal power of the plant P nom , is subtracted in addition to using a constant load-unspecific efficiency a . The parameterisation of a and b allows for modelling the non-linear efficiency characteristic with sufficient accuracy (see Section 2.6.3). In comparison to a piecewise linearisation, less variables and especially no additional binary variables ought to be used, which leads to lower computational effort [41].

2.5.2. Constraints to Model Energy Storage Devices

For units of the generic storage class s S , power and ramping limits are formulated according to Section 2.5.1. It is distinguished between the operating conditions “charging” and “discharging”:
P ch ( s , t ) P ch , min ( s ) ψ ch ( s , t ) ; s S t T k . P charge ( s , t ) P ch , max ( s ) ψ charge ( s , t ) ; s S t T k . P discharge ( s , t ) P dis , min ( s ) ψ dis ( s , t ) ; s S t T k . P discharge ( s , t ) P dis , max ( s ) ψ dis ( s , t ) ; s S t T k .
P ch ( s , t 1 ) P ch ( s , t ) R ch , down ( s ) ; s S t T k . P ch ( s , t 1 ) P ch ( s , t ) R ch , up ( s ) ; s S t T k . P dis ( s , t 1 ) P dis ( s , t ) R dis , down ( s ) ; s S t T k . P dis ( s , t 1 ) P dis ( s , t ) R dis , up ( s ) ; s S t T k .
For t = T k , 0 the variable for t 1 is replaced by data from the non-linear control of the previous period. Because the operating conditions are exclusive, simultaneous charging and discharging of the storage is inhibited by Equation (6):
ψ ch ( s , t ) + ψ dis ( s , t ) 1 ; s S t T k .
The efficiency characteristics for both operating conditions are modelled also using the load-unspecific basic demand as described in Section 2.5.1 (see Equation (5)):
P ch , u ( s , t ) = a ch ( s ) P ch ( s , t ) a ch ( s ) b ch ( s ) P ch , nom ( s ) ψ ch ( s , t ) ; s S t T k . P dis , u ( s , t ) = a dis ( s ) P dis ( s , t ) a dis ( s ) b dis ( s ) P dis , nom ( s ) ψ dis ( s , t ) ; s S t T k .
To apply the efficiency losses for a storage unit, a distinction is made between the power on the grid side P ch and the power as seen by the storage P ch , u .
Since Equation (7) yields the possibility for no power on the storage but power on the grid side, which in return is beneficial for the objective function (14), the following two equations are necessary and state that the efficiency of each operating mode ought to be always below one:
P ch , u ( s , t ) P ch ( s , t ) ; s S t T k . P dis , u ( s , t ) P dis ( s , t ) ; s S t T k .
The model only considers the net capacity of a storage C net . The upper and lower bounds for the storage level E can therefore be set as follows:
E _ ( s , t ) = 0 ; s S t T k . E ¯ ( s , t ) = C net ( s ) ; s S t T k .
The storage level of an energy storage device is calculated based on the storage level of the previous time step, the self-discharge rate ξ and the charging/discharging power on the storage side P ch , u and P dis , u , respectively, in the current time step:
E ( s , t ) = E ( s , t 1 ) [ 1 ξ ( s ) τ ] + τ [ P ch , u ( s , t ) P dis , u   ( s , t ) ] ; s S   t T k .
For t = T k , 0 the variable for t 1 is replaced by data from the non-linear control of the previous period. To ensure a closed energy balance over the whole simulation time frame, boundary conditions need to guarantee the equality of the initial storage level in the first interval of the simulation and the final storage level of the last interval. One problem regarding the rolling horizon approach is the missing communication between the first and the last loop as to what initial, and therefore final, storage level might be suboptimal.
To surpass that problem, the initial storage level E init is passed to the model as data point. The value for each storage unit is received via the UCH approach described in Section 2.4. The following constraints ensure a closed energy balance for each storage within a tolerance of one percent net capacity:
E ( s , 0 ) = E init ( s ) ;   s S ;   if   first   period
and:
E _ ( s , t period ) = E init ( s ) ± 0.01 C net ( s ) ;   s S E ¯ ( s , t period ) = E init ( s ) ± 0.01 C net ( s ) ;   s S }   if   t period = T total

2.5.3. Constraints to Model the Energy System

The electric demand P demand and the renewable power generation P volatile are given as input data for the system balance, which is formulated in Equation (13), to match demand and generation in the grid for each time step:
P demand ( t ) = P volatile ( t ) + G P ( g , t ) S P ch ( s , t ) + S P dis ( s , t ) + P slack + ( t ) P slack ( t ) ; t T k .
Slack variables P slack for a positive and negative deviation of load and generation are established to ensure solvability of the optimisation problem.

2.5.4. Objective Function

As a model of the electricity market is not yet integrated in the optimisation problem, the objective function (14) is set to minimise the emissions throughout the regarding time frame:
min { T τ [ G P u ( g , t ) e ( g ) S P ch , u ( s , t ) e stor + S P dis , u ( s , t ) e stor + P slack + ( t ) e pen + + P slack ( t ) e pen ] }
The usage of the slack variables should be avoided to achieve equality of demand and supply. However, in time steps where storages are full and there is high renewable generation, those slack variables are needed to find valid solutions. Therefore, the specific emissions for the slack variables e pen are set to high values as a penalty, so that those variables are only contributing to the solution when no other can be found.
The penalty for power surpluses serves two functions in that it not only serves the solution process, but also incentivises the use of storage devices. This becomes necessary when the excess of energy cannot be used within the current interval. The virtual storage emissions are set to incentivise the usage of more efficient storages. When not in place an overuse of the less efficient long-term storage was observed as the long-term storage device may reduce conventional and penalty usage by the same amount as the short-term storage within a single interval.

2.6. Temporal and Technical Scenario Configuration

Throughout this paper, a scenario is defined as a set of period and interval length and configuration is defined as the installed powers, capacities and set of units. Due to the complexity of the constraints and the high temporal resolution of the data, a copper plate island with a rather small number of actors is chosen. This makes the result not applicable to larger systems, which is acceptable for this paper, since the effect of different interval lengths in UCO is the main research target of this paper rather than an accurate representation of a real energy system.
The interval lengths were chosen to represent possible ranges from forecast lengths of high resolution, local weather models [42], like the two and three days intervals, up to forecast lengths of low resolution, global weather models [43] with the 5 and 14 days intervals. The 30 days interval was set to test whether there is need for even longer forecast lengths and the 365 days interval was set as it is commonly used in UC models [6,7,8,41,44]. To achieve or rather allow a similar computation time the MILP-solver settings have to be adjusted for the varying interval lengths. Those settings can be seen in Table 1, where interval and period lengths are given as time in days.
For the 365 days scenario the maximum time abort criteria was set comparatively low, since the program tended to crash otherwise.
Abort criteria for the MILP solver to stop its solution process can be found in Table 1 as well. Abort criteria are the time for the solution process and the relative gap. The values were chosen to limit the overall solution time to a similar value, but do not necessarily lead to equal overall solution times (see Table A5, Table A6, Table A7, Table A8, Table A9 and Table A10). Other than that, the default settings of FLOP++ and CPLEX for a branch-and-bound solution were used.

2.6.1. The Actors

For the sake of simplicity and reproducibility, a constant demand of 1000 MW all throughout the year was chosen. This is deemed acceptable for the research target at hand as the UC problem resulting from the VRE generation is seen as sufficiently complex. As conventional power plant, using the primary storage natural gas, a combined cycle gas turbine (CCGT) power plant was chosen with a rated power output of 1000 MW, so that it is always able to cover any remaining residual load. As a short-term energy storage device, a battery with 5000 MWhnet is chosen. A long-term energy storage device consists of an electrolyser for conversion of power to hydrogen and its own CCGT plant for reconversion of hydrogen to power. As storage capacity 240,000 MWhnet is chosen. The storage capacities translate to their respective storage type. Both storages are rated at 1000 MW, which for the electrolyser results in a higher maximum load of around 1600 MW due to its overload capacities [45]. The installed rated powers of the VRE types solar and wind are set between 1000 and 5000 MW each, referred to as MWeach in the following (see Table 2), since both technologies being equally sized in all configurations. The installed VRE power in the following graphs is to be understood as the value both units each are rated at.

2.6.2. Parameterisation of Initial Storage Levels

As previously stated, the initial storage cannot be optimised in a rolling horizon approach and need to be entered into the model as input data. The following initial states of charge are used in the simulations for E init (see Equations (11) and (12)) of the long-term storage.
The values of Table 2 were determined from the storage level calculated using the heuristic approach for the respective configurations, so that the results are comparable. The 365 days interval (scenario f) is the exception to this rule, since it had access to both initial and end storage level within a single interval. Therefore in the respective scenarios E init is implemented as a variable rather than input data, so that the UCO was allowed to optimise those storage levels for the case of single interval to solve.

2.6.3. Parameterisation of Efficiency Characteristics

As described in Section 2.5.1 and Section 2.5.2, the efficiency of all plants is modelled using the load-unspecific base demand. To obtain the necessary parameters a and b , the optimisation problem stated in Equation (15) is solved for each individual plant. Therefore, (non-linear) characteristics stating efficiencies η at corresponding relative loads p are used as input data. The number of elements contained in the input characteristics is described by set I . Efficiencies lower than 10% will be ignored, due to the efficiency being the denominator, which gives low values too much weight:
min ( i = 1 I δ + ( i ) + δ ( i ) )
subject to:
1 η ( i ) = 1 a + b 1 p ( i ) + δ + ( i ) δ ( i ) ;   i I
1 η max 1 a + b 1 p ( i ) ;   i I
1 / a is the algebraic representation. Programatically this is solved by using a new variable representing the inverse of a to achieve linearity.
In the optimisation problem δ + and δ represent the positive and negative deviations between the parameterised efficiencies and the input data, which are necessary as slack variables and ought to be minimised.
In Figure 5 the non-linear efficiency curves used in the UCH (dashed) are shown. The curve for the CCGT is determined scaling the relative efficiency presented in [27] and shown in Figure 1 with a maximum efficiency of 61%, while the electrolyser curve stems from [45]. Those curves are also given as input data for the parameter optimisation of Equation (15). The resulting curves (dotted) are also visualised in Figure 5.
The short term storage is modelled using a constant efficiency, resulting from the square root of round-trip efficiency stated in [46], for both charging and discharging operation. This in return results in parameter b set to zero by the optimisation and a constant efficiency for the optimiser. The resulting parameters a and b used in the UCO are stated in Table 3.
Minimum and maximum discharge power for the battery are 0 and 100% respectively as it is determined by the battery management system [46], all other minimum and maximum power operation points can be inferred from Figure 5. Additional parameters used in the simulation are listed in Table A2 in Appendix A.

3. Results

In a first step, the usage of the UCH is compared to the UCO with 2 days interval rolling horizon (scenario a) for all VRE power configurations. Figure 6 shows the specific CO2 emissions (dots) of the plants of the copper plate island and the storage share (square) in dependence on installed VRE power. The results show that the UCO model (dotted lines) leads to slightly lower CO2 emissions than the UCH (solid lines) for all configurations.
In terms of the storage share, the UCO delivers a widely different result than the UCH, with higher values in all configurations but especially for small installed VRE power. This is further discussed in Section 3.1. Both UCH and UCO show a maximum share of storages below the maximum installed VRE power. This is due to the ever-rising share of direct VRE supply, leaving less demand to be covered by the energy storage devices.
Lastly, Figure 6 explains why the 1000 MWeach configuration shows a lower absolute improvement than the 2000 MWeach configuration, with the storage share of the UCH. Since the UCH only integrates surplus of VRE power supply, it is obvious that the 1000 MWeach configuration has negligible surplus of VRE to be integrated, thus the UCO has only the CCGT left to optimise. In practice, this is done using the short-term storage.
The reduction of CO2 emissions due to optimisation in comparison to the UCH is shown in absolute and relative terms in Figure 7 also for the interval length of two days. The figure shows that the relative CO2 reduction through UCO for installed VRE power below 3000 MWeach is generally below 10%. Between 3000 and 3750 MWeach the absolute reduction is rather constant across VRE configurations while the relative reduction shows increasing values. For 4000 MWeach both absolute and relative reduction show lower values, which can be explained with the maximum share of storage to be integrated as seen in Figure 6 and discussed above. For higher installed capacities, the denominator becomes so insignificant that the relative reduction is amplified while the absolute CO2 reduction becomes remarkably small.
The main result can be seen in Figure 8, where the reduction of CO2 emissions through UCO relative to the heuristic approach (UCH) is shown in relation to the interval length, which in this context implies the corresponding period lengths of Table 1. The configurations with 1000 MWeach (circles) and 2000 MWeach (squares) VRE power show no significant influence of the interval length across the scenarios (interval lengths). The configuration with 3000 MWeach shows a slight advantage for a whole year optimisation run, which is explicable with the maximum possible integration of storage power (see Figure 6) for the assumed residual power units’ configuration. For the 4000 and 5000 MWeach configurations the whole year optimisation run shows, against expectations, negative emission reduction, which can be explained by the linearisation effects, that are discussed in Section 3.2 and visible in Figure 10.
All interval lengths between 2 and 30 days with their corresponding period lengths do not seem to have a major influence on the quality of the optimisation result, but show significant changes in computational time (see Table A3 in Appendix A). Further, the graph shows that the absolute reduction of emissions through UC decreases with increased VRE power installation, the exception being the 2000 MWeach configuration, as discussed above.

3.1. Result of Unit Commitment

The goal of optimisation is to find non-trivial solutions [44]. Figure 9, depicting a time series excerpt, clearly shows a non-trivial result from the implemented UCO model (dashed and dotted lines). For the UCH the long-term storage device’s discharge unit (solid line) follows the residual load as the UCH emptied the short-term storage device already at that time frame (therefore not shown), which results from the pre-set order of operation. For the hours 11:30 to 14:15 however the long-term storage device’s discharge unit cannot go below its minimum load and therefore feeds to the grid with minimum load, creating the need for curtailment. The UCO makes a better use of storage devices in that the long-term storage (dotted line) rarely follows the residual load, but rather is running mostly at full power or not at all. The short-term storage device (dashed line) is charged with the surpluses created by the long-term storage device and discharges when the long-term storage device is not utilised. This happens predominantly at low residual loads, where the long-term storage device’s discharge unit would have a relative efficiency lower than the round-trip efficiency of the short-term storage device. While this seems an unrealistic unit dispatch, it is the most efficient plan, in energy terms, for the formulated mathematical problem. Hence, this can be considered an optimal strategy—especially as the curtailment due to minimum load of the long-term storage device’s discharge unit is avoided. In Figure 9, a scattering of the UC can be observed and is deemed unrealistic. It can be prevented by including intertemporal constraints such as minimum downtimes and start-up times in the mathematical problem formulation, which has not been implemented in the UCO version used in this paper. The resulting scattering of the tactic is further discussed in Section 4.
As previously discussed, the biggest difference between the UCO and the UCH can be seen in the share of storage in the supplied energy. This results from the different use of the energy storage devices. While the UCH only stores surpluses of renewables, the optimisation finds non-trivial solutions where the CCGT of the long-term storage is used at its best efficiency, charging the battery to avoid a less efficient operating point later and instead discharging the battery. This is the expected behaviour of the optimisation routine, since its target is to minimise the CO2 emissions and the efficiency spread of the CCGT between minimum and maximum operating point is 20%-rel. and the losses due to battery usage lie around 10%.

3.2. Results for Different Interval Lengths

As described in Section 3.1, the characteristic line model acts as an NLC for the optimised operation plan. While the algorithm tries to set the power of each unit according to plan, it does not use the linearised efficiency curves shown in Figure 5 as dotted lines, but the non-linear characteristics (dashed lines). This results in differences between the direct MILP solution and the final UCO solution after non-linear control by algorithm.
For storage units that means that the storage level might be over- or underestimated and the unit is not running with the planned power or, due to those discrepancies in storage level, at all. For the conventional residual unit, the difference results in only negligibly different emissions when running as planned. Due to the storage deviations unplanned run-times may occur for the conventional units. These run-times result e.g., from timesteps, where the MILP can use the energy storage devices because of non-zero SOC, but the NLC calculated zero SOC (see Figure 10) and needs to use the conventional power plants to cover the residual load. The timesteps and duration of the different run-times vary with interval and period length of the rolling horizon.
For a closer look, time series of the storage level are shown in Figure 10 for the long-term storage with an installed VRE power of 3000 MWeach. The dashed lines represent the MILP results, while the solid lines show the storage level after the NLC. The time series are presented for a 5 days interval (scenario c, blue lines) and a 365 days optimization (scenario f, green lines).
It can clearly be seen, that the storage level behaviours for both scenarios differ widely, especially from the end of January to November. In case of the 5 days interval, the storage is discharged completely at the end of January. Afterwards, the energy in the storage is charged only up to 80 GWh of energy once, with various other minor peaks, but no significant accumulation of energy. Lastly, starting in November, enough renewable surplus leads to an accumulation of energy, until the initial storage level is either reached or exceeded again.
In the 365 days interval, with perfect foresight for the entire time frame, the long-term storage is not depleted at the beginning of the year, but operated between 50 and 90% SOC. The peaks throughout the year, however, qualitatively correspond to the 5 days interval time series. After a complete discharge at the end of November, both scenarios progress with an almost equal storage level.
Furthermore, it can clearly be seen that a significant difference in storage level is mounting up as long as the zero SOC point is avoided by the MILP solution. The reason lies in the parameterisation of the long-term charge unit efficiency curve, as displayed in Figure 5. The optimiser calculates a 13 rel.-% higher efficiency for maximum power than the NLC, which in return accumulates a higher storage level for long operating hours at higher charging power (see Figure 10). As Figure 5 suggests we observe an underestimation for operating hours at low charging power.
When looking at the same VRE configuration, but shorter interval length of 5 days (blue lines) the differences in storage level between the MILP solution and the solution after NLC become significantly smaller. This is due to the rolling horizon, which is aligning the optimisation storage level more frequently (at the end of each optimisation period for a new interval) with the NLC results. For the 365 days interval and period length, there is no feedback loop, so that differences can accumulate (see Figure 10).
This is also the main cause why seemingly optimal solutions for long optimisation intervals might be found, but are not necessarily feasible by the non-linear control calculation and therefore lead to higher emissions than even the simplest heuristic approach.
While the linearisation effects mainly influence scenarios with long optimisation intervals (in this paper especially the 365 days optimisation in scenario f), the VRE configuration influences the magnitude of those effects. With higher installed VRE capacity, the usage of the long-term storage with high loads increases, so that higher deviations due to linearisation effects can accumulate. This can be observed in Figure 8.

4. Discussion

Several points need to be discussed in order to correctly categorise the findings of this paper. As previously stated, the constant power demand does not represent a real grid situation. The complexity still was deemed sufficient for the test of interval lengths, which is supported by Figure 9, where the UCO comes to a different solution than the UCH in almost every time step. The quantitative values of specific emissions and storage share therefore do not represent a prognosis for future energy systems. Additionally, it is emphasized, that the configuration of residual power supply was not optimised. Lastly, with the chosen boundary conditions it could be argued, that either the efficiency curves or the number of units should have been adjusted to correctly represent the installed powers. This was deemed too complex within the scope of this research.
The linearisation effects presented in Section 3.2 can mainly be attributed to the efficiency curve parameterisation in Section 2.6.3. A different formulation of the optimisation target might resolve in a closer fit around maximum power. The results for high installed VRE power might be improved by that. An adjustment of the parameter optimisation’s objective (see Equation (15)) was not done however, since the 365 days interval is considered unrealistic and other interval lengths were matching the NLC’s result close enough (see Figure 10). This rather shows the missing robustness of a 365 days interval length for the chosen problem.
A shorter period length for the 365 days interval might fix the previously mentioned problems, but that was considered too distant from the scope of this paper due to large computational effort. In general, the results for the different interval lengths might be tied closely to the corresponding period lengths. This difference is only addressed by the 2 and 3 days interval lengths, where no significant effect occurred.
Another point that needs discussion is the scattering observed in Figure 9. The CCGT of the long-term storage is working for several short time periods following each other rather closely. This is a behaviour that would not be expected in reality since there are temporal constraints on down-time and maximum possible ramps to consider. This was deemed beyond of scope for this paper, but will be worked on to better approximate real behaviour. As it is expected to change the results only in quantity and appearance it is seen as a minor problem for the findings of this paper.
As a last shortcoming of the model it should be emphasized, that the implemented heuristic is very simple. It could be argued that a more complex heuristic of conditions and loops could meet the results of the optimisation more closely. A rather simple heuristic was chosen, since the results were meant to cover the baseline expectation.

5. Conclusions

In this paper, a different modelling approach for rolling horizon UC was presented. For that approach, a UCO model was implemented and its positive effect in contrast to a simple UCH could be clearly shown. As the main result it is observed, that the interval length for a rolling horizon approach has an almost negligible effect on the objective value, namely the CO2 emissions of the energy park, with this approach, even though, the detailed unit commitment can vary significantly (see Figure 10)). However, the difference is mainly to be found in computational effort.
A 2 days optimisation interval is therefore suggested as sufficient considering quantitative reduction of the optimisation target, computation time (see solver statistics in the Appendix A) and realistic lengths of weather forecasts [42,43].
The variation of the configuration (installed volatile power capacity in this study) however has a bigger effect on the optimisation objective (CO2 emissions in this study) than the UCO. Considering the massive computation time of a UCO, it is suggested to explore configuration variations using a simple heuristic and to optimise a small set of resulting configurations with a UCO only afterwards.
The effects of the linearisation of a given system as discussed in Section 4 make a strong case for rolling horizon and against whole-year optimisation runs. Additionally, a system configuration optimised according to a whole year optimisation run might lead to different results than a system configuration optimised according to a rolling horizon run. Since the forecast length that is inherent to a whole year run is not achievable, it is suggested to optimise an energy system configuration with a rolling horizon approach as the final configuration of any energy system should work with realistically achievable foresight.

Outlook

While the proposed optimisation model already allows a detailed technical description of the energy system, there are several possibilities for model improvement in upcoming versions. As discussed in Section 4, the optimiser’s results tend to include scattering of generation and storage units. Intertemporal constraints could resolve that problem. Such intertemporal constraints will be included and looked upon in future versions. It is of great interest how that effects both the solution and the computation time since the problem is more complex but also tighter, which might make the decision trees of the branch-and-bound considerably smaller.
In this paper, the optimisation goal is set to minimise the CO2 emissions of the underlying energy system. Another common approach is the minimisation of (operational) costs, which is currently being implemented allowing multi-objective optimisations (costs and CO2 emissions).
A major problem with using a rolling horizon for UC problems is the control of storage levels for long-term storages, especially with short interval length, because the optimiser is not aware and therefore does not plan for eventually upcoming periods of continuously low supply of VRE power (“dark doldrums”). The use of long-term storage strategies e.g., by the use of minimum support points is an interesting topic for investigation.
Lastly, it would be interesting to see whether the findings of this paper are applicable to different configurations and more complex systems.

Supplementary Materials

Software repository is under https://collaborating.tuhh.de/iet/CharL. Commit “e451c36f“ of 17 September 2018 was used for this research. Due to licensing issues the access is currently granted on an individual basis. Additionally, we used (all links last accessed on 13 March 2019): FLOPC++ found under https://projects.coin-or.org/FlopC++ (The FLOPC++ implementation of revision 362 was altered so that a result reaching the time limit is still considered “optimal” if the relative gap is below 10%.); OSI found under https://projects.coin-or.org/Osi (OSI stands for Open Solver Interface); Qt Library found under https://www.qt.io/.

Author Contributions

Conceptualization, G.E. and T.Z.; methodology, G.E. and T.Z.; software, G.E. and T.Z.; validation, none; formal analysis, G.E., T.Z. and A.K.; investigation, G.E. and T.Z.; resources, G.E. and T.Z.; data curation, G.E. and T.Z.; writing—original draft preparation, G.E., T.Z. and A.K.; writing—review and editing, G.E., T.Z. and A.K.; visualization, G.E. and T.Z.; supervision, A.K.; project administration, A.K.; funding acquisition, G.E.

Funding

Funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)–Projektnummer 392323616 and the Hamburg University of Technology (TUHH) in the funding programme Open Access Publishing.

Acknowledgments

We would like to thank Christopher Ball for language support. Funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)–Projektnummer 392323616 and the Hamburg University of Technology (TUHH) in the funding programme Open Access Publishing.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Nomenclature

SymbolUnitDescription
a -Efficiency parameter
b -Efficiency parameter
C net MWhnetNet capacity of storage
E MWhStorage level
e kg/MWhthSpecific CO2 emissions
G -Set of conventional generation units
g -Index of conventional generation unit in set G
I -Set of points in efficiency curve
i -Index of points in set I
k -Index of optimisation period
K -Set of optimisation periods
P MWPower of unit on grid side
P u MWPower of unit on unit side
p -Relative load
R MW/minMaximum power gradient of unit
S -Set of storage units
s -Index of storage unit in set S
T -Set of timesteps
T k -Set of timesteps of interval for period k
T k , 0 -First timestep of interval for period k
t -Timestep in Set T
t period -Number of timesteps in a period
δ -Deviation in parameterisation problem
η -Efficiency
ξ %/hSelf-discharge rate of storage
τ hLength of timestep
ψ -Status variable for unit (0 = off, 1 = on)

Terms and Abbreviations

Term/AbbreviationDescription
ConfigurationTypes, powers and capacities of installed units
ScenarioSet of interval and period length
MWeachInstalled power for each solar and wind generation in MW
CCGTCombined Cycle Gas Turbine
CharLCharacteristic Line Model
FlopC++Formulation of Linear Optimization Problems in C++
MILPMixed Integer Linear Problem
OSIOpen Solver Interface
SOCState of Charge
UCHUnit Commitment Heuristic
UCOUnit Commitment Optimisation
NLCNon-linear control (model)
TSOTransmission System Operator
VREVolatile Renewable Energy

Appendix A

Table A1. Hardware setup used for the calculations.
Table A1. Hardware setup used for the calculations.
Setup NameProcessorRAM
PC 1Intel Xeon E5-1620 @3.6 GHz64 GB RAM UDIMM
PC 2Intel i7-3770 @3.4 GHz16 GB RAM UDIMM
PC 3Intel i7-4790 @3.6 GHz32 GB RAM UDIMM
PC 4Intel i7-6700 @3.4 GHz32 GB RAM UDIMM
Table A2. Additional parameters used in the simulation.
Table A2. Additional parameters used in the simulation.
ParameterUnitValueExplanation/Source
CCGT Ramp Rate%/min100Beyond of scope
CCGT Start-up Timesh0Beyond of scope
Specific Fuel EmissiontCO2/MWhth0.202[47]
Ramp Rate Storages (all)%/min100Beyond of scope
Start-up Times Storages (all)h0Beyond of scope
Self-discharge Rate Short-Term%/h0.01[48] (Middle of spread)
Self-discharge Rate Long-Term%/h0.0006875[46] (Middle of spread)
Virtual positive penalty emissiontCO2/MWh106-
Virtual negative penalty emissiontCO2/MWh102-
Virtual storage emissiontCO2/MWh10−3-
Table A3. Specific CO2 emissions for simulated scenarios.
Table A3. Specific CO2 emissions for simulated scenarios.
Installed Power Wind and Solar in MW EachSpecific CO2 Emissions in g/kWh
Heuristic2 Days3 Days5 Days14 Days30 DaysYear
1000236.77230.06230.06230.07230.10230.10230.11
1250214.13202.76-----
1500187.69175.55175.53----
1750161.06149.72-----
2000136.29126.16126.13126.13126.17126.19126.39
2250113.34104.43-----
250092.2684.5184.48----
275072.7165.95-----
300054.6449.5249.4849.4749.4649.4648.88
325040.7235.74-----
350028.8723.7023.65----
375018.0812.77-----
400010.318.358.348.398.508.7112.48
42507.745.67-----
45005.243.223.20----
47502.921.00-----
50000.740.000.000.000.020.038.59
Table A4. Share of storages for simulated scenarios.
Table A4. Share of storages for simulated scenarios.
Installed Power Wind and Solar in MW EachShare of Storages in %
Heuristic2 Days3 Days5 Days14 Days30 DaysYear
10000.0112.5012.5612.6212.5912.6212.53
12500.2612.37-----
15001.1712.1112.11----
17502.7212.16-----
20004.5512.6112.5912.6312.6112.6112.50
22506.6113.40-----
25008.8814.4014.35----
275011.1315.67-----
300013.4216.8916.8616.8516.7616.6916.92
325015.0317.75-----
350016.2918.5018.50----
375017.4319.26-----
400017.9118.6718.6618.6818.5718.6217.96
425017.1817.88-----
450016.5717.1817.19----
475016.0416.55-----
500015.5215.7215.7215.7215.6315.6313.85
Table A5. Statistics for 2 days optimisation of MILP solver.
Table A5. Statistics for 2 days optimisation of MILP solver.
Installed Power Each Wind and Solar in MWStatistics of Rolling Horizon Regarding the Time of the SolutionStatistics of Rolling Horizon Regarding the Relative Gap of the Solutions All in -
Average Time in s# of Times Time Limit ReachedMin.Max.AverageStandard Deviation# of Times No Solution Found
1000255.226802.81 × 10−58.30 × 10−57.54 × 10−40
1250296.508107.53 × 10−53.35 × 10−44.14 × 10−30
1500232.866106.01 × 10−52.28 × 10−42.74 × 10−30
1750249.276501.13 × 10−47.34 × 10−41.26 × 10−20
2000298.448002.08 × 10−41.47 × 10−32.61 × 10−21
2250350.709601.36 × 10−46.60 × 10−41.00 × 10−20
2500373.6710202.79 × 10−41.56 × 10−32.25 × 10−20
2750331.379101.39 × 10−44.88 × 10−44.60 × 10−34
3000357.329803.48 × 10−43.67 × 10−36.99 × 10−24
3250361.4999−4.00 × 10−63.47 × 10−41.65 × 10−32.23 × 10−24
3500306.087905.56 × 10−43.56 × 10−34.40 × 10−22
3750303.4183−3.31 × 10−34.83 × 10−42.32 × 10−32.76 × 10−22
4000372.8310304.84 × 10−41.52 × 10−31.21 × 10−22
4250385.11109−5.83 × 10−47.64 × 10−43.08 × 10−34.95 × 10−21
4500416.41116−1.78 × 10−47.68 × 10−42.41 × 10−31.99 × 10−23
4750396.9711408.16 × 10−43.28 × 10−34.81 × 10−21
5000474.4013501.06 × 10−33.79 × 10−35.75 × 10−21
Table A6. Statistics for 3 days optimisation of MILP solver.
Table A6. Statistics for 3 days optimisation of MILP solver.
Installed Power Each Wind and Solar in MWStatistics of Rolling Horizon Regarding the Time of the SolutionStatistics of Rolling Horizon Regarding the Relative Gap of the Solutions all in -
Average Time in s# of Times Time Limit ReachedMin.Max.AverageStandard Deviation# of Times No Solution Found
1000625.848902.07 × 10−55.58 × 10−55.55 × 10−40
1500696.9610105.91 × 10−52.19 × 10−43.03 × 10−30
2000936.5113801.79 × 10−41.16 × 10−32.07 × 10−20
25001125.3315802.78 × 10−41.97 × 10−33.29 × 10−20
30001064.41151−1.00 × 10−62.55 × 10−41.58 × 10−32.91 × 10−21
3500900.4712501.01 × 10−37.11 × 10−38.81 × 10−22
4000931.9713107.44 × 10−43.31 × 10−35.32 × 10−21
45001126.8016501.18 × 10−33.60 × 10−33.54 × 10−20
50001136.87164−2.73 × 10−31.63 × 10−35.42 × 10−34.91 × 10−24
Table A7. Statistics for 5 days optimisation of MILP solver.
Table A7. Statistics for 5 days optimisation of MILP solver.
Installed Power Each Wind and Solar in MWStatistics of Rolling Horizon Regarding the Time of the SolutionStatistics of Rolling Horizon Regarding the Relative Gap of the Solutions All in -
Average Time in s# of Times Time Limit ReachedMin.Max.AverageStandard Deviation# of Times No Solution Found
10001105.311041.00 × 10−62.37 × 10−53.78 × 10−52.55 × 10−40
20002217.89107.501.41 × 10−43.98 × 10−42.87 × 10−30
30002445.49120.502.03 × 10−44.67 × 10−45.54 × 10−30
40002205.14106.509.62 × 10−42.82 × 10−32.26 × 10−20
50002229.26109.501.26 × 10−33.21 × 10−32.66 × 10−20
Table A8. Statistics for 14 days optimisation of MILP solver.
Table A8. Statistics for 14 days optimisation of MILP solver.
Installed Power Each Wind and Solar in MWStatistics of Rolling Horizon Regarding the Time of the SolutionStatistics of Rolling Horizon Regarding the Relative Gap of the Solutions All in -
Average Time in s# of Times Time Limit ReachedMin.Max.AverageStandard Deviation# of Times No Solution Found
10001975.74137.00 × 10−63.16 × 10−55.46 × 10−52.41 × 10−40
20006005.35381.00 × 10−51.16 × 10−42.14 × 10−41.51 × 10−30
30006487.444302.51 × 10−43.97 × 10−42.29 × 10−30
40003013.661807.99 × 10−42.80 × 10−32.00 × 10−20
50003940.3825−1.30 × 10−52.40 × 10−36.98 × 10−34.33 × 10−20
Table A9. Statistics for 30 days optimisation of MILP solver.
Table A9. Statistics for 30 days optimisation of MILP solver.
Installed Power Each Wind and Solar in MWStatistics of Rolling Horizon Regarding the Time of the SolutionStatistics of Rolling Horizon Regarding The Relative Gap of the Solutions All in -
Average Time in s# of Times Time Limit ReachedMin.Max.AverageStandard Deviation# of Times No Solution Found
100011,190.5816.51.30 × 10−55.79 × 10−57.75 × 10−53.42 × 10−40
200021,609.5436.54.50 × 10−52.08 × 10−41.54 × 10−46.99 × 10−40
300020,163.36341.50 × 10−57.14 × 10−46.50 × 10−42.52 × 10−30
40007155.091202.77 × 10−46.53 × 10−43.39 × 10−30
50009712.9915.502.56 × 10−39.59 × 10−35.88 × 10−20
Table A10. Statistics for 365 days optimisation of MILP solver.
Table A10. Statistics for 365 days optimisation of MILP solver.
Installed Power Each Wind and Solar in MWStatistics of Rolling Horizon Regarding the Time of the SolutionStatistics of Rolling Horizon Regarding the Relative Gap of the Solutions All in -
Average Time in s# of Times Time Limit ReachedRelative GapStandard Deviation# of Times No Solution Found
1000100,036.5411.77 × 10−400
2000100,047.7414.43 × 10−400
3000100,076.5611.04 × 10−300
4000100,092.6513.19 × 10−200
5000100,053.8616.59 × 10−300

References

  1. United Nations. Framework Convention on Climate Change. In Adoption of the Paris Agreement; FCCC/CP/2015/L.9/Rev.1; United Nations: New York, NY, USA, 2015. [Google Scholar]
  2. Buttler, A.; Spliethoff, H. Bedarf und Auslastung konventioneller Kraftwerke im Zuge der Energiewende: Eine Metastudie. VGB PowerTech 2016, 3, 41–47. [Google Scholar]
  3. Sterner, M.; Stadler, I. Energiespeicher. Bedarf, Technologien, Integration; Springer: Berlin/Heidelberg, Gernamy, 2014; ISBN 978-3-642-37379-4. [Google Scholar]
  4. Abujarad, S.Y.; Mustafa, M.W.; Jamian, J.J. Recent approaches of unit commitment in the presence of intermittent renewable energy resources: A review. Renew. Sustain. Energy Rev. 2017, 70, 215–223. [Google Scholar] [CrossRef]
  5. Ommen, T.; Markussen, W.B.; Elmegaard, B. Comparison of linear, mixed integer and non-linear programming methods in energy system dispatch modelling. Energy 2014, 74, 109–118. [Google Scholar] [CrossRef]
  6. Mehleri, E.D.; Sarimveis, H.; Markatos, N.C.; Papageorgiou, L.G. Optimal design and operation of distributed energy systems: Application to Greek residential sector. Renew. Energy 2013, 51, 331–342. [Google Scholar] [CrossRef]
  7. Henning, H.-M.; Palzer, A. A comprehensive model for the German electricity and heat sector in a future energy system with a dominant contribution from renewable energy technologies: Part I: Methodology. Renew. Sustain. Energy Rev. 2014, 30, 1003–1018. [Google Scholar] [CrossRef]
  8. Scholz, Y. Renewable Energy Based Electricity Supply at Low Costs: Development of the REMix Model and Application for Europe. Ph.D. Dissertation, Universität Stuttgart, Stuttgart, Germany, January 2012. [Google Scholar]
  9. Sarker, M.R.; Wang, J. Unit Commitment on the Cloud. 2017. Available online: http://arxiv.org/pdf/1702.03886v1 (accessed on 13 March 2019).
  10. Morales-España, G.; Latorre, J.M.; Ramos, A. Tight and Compact MILP Formulation for the Thermal Unit Commitment Problem. IEEE Trans. Power Syst. 2013, 28, 4897–4908. [Google Scholar] [CrossRef]
  11. Silvente, J.; Kopanos, G.M.; Pistikopoulos, E.N.; Espuña, A. A rolling horizon optimization framework for the simultaneous energy supply and demand planning in microgrids. Appl. Energy 2015, 155, 485–501. [Google Scholar] [CrossRef]
  12. Stadler, P.; Ashouri, A.; Maréchal, F. Model-based optimization of distributed and renewable energy systems in buildings. Energy Build. 2016, 120, 103–113. [Google Scholar] [CrossRef]
  13. Ringkjøb, H.-K.; Haugan, P.M.; Solbrekke, I.M. A review of modelling tools for energy and electricity systems with large shares of variable renewables. Renew. Sustain. Energy Rev. 2018, 96, 440–459. [Google Scholar] [CrossRef]
  14. Silvente, J.; Kopanos, G.M.; Dua, V.; Papageorgiou, L.G. A rolling horizon approach for optimal management of microgrids under stochastic uncertainty. Chem. Eng. Res. Des. 2018, 131, 293–317. [Google Scholar] [CrossRef]
  15. Costley, M.; Feizollahi, M.J.; Ahmed, S.; Grijalva, S. A rolling-horizon unit commitment framework with flexible periodicity. Int. J. Electr. Power Energy Syst. 2017, 90, 280–291. [Google Scholar] [CrossRef]
  16. Chen, H.; Xuan, P.; Wang, Y.; Tan, K.; Jin, X. Key Technologies for Integration of Multitype Renewable Energy Sources—Research on Multi-Timeframe Robust Scheduling/Dispatch. IEEE Trans. Smart Grid 2016, 7, 471–480. [Google Scholar] [CrossRef]
  17. Schulze, T.; McKinnon, K. The value of stochastic programming in day-ahead and intra-day generation unit commitment. Energy 2016, 101, 592–605. [Google Scholar] [CrossRef]
  18. Saint-Pierre, A.; Mancarella, P. Active Distribution System Management: A Dual-Horizon Scheduling Framework for DSO/TSO Interface Under Uncertainty. IEEE Trans. Smart Grid 2017, 8, 2186–2197. [Google Scholar] [CrossRef]
  19. Hermans, M.; Bruninx, K.; Delarue, E. Impact of CCGT Start-Up Flexibility and Cycling Costs Toward Renewables Integration. IEEE Trans. Sustain. Energy 2018, 9, 1468–1476. [Google Scholar] [CrossRef]
  20. Bischi, A.; Taccari, L.; Martelli, E.; Amaldi, E.; Manzolini, G.; Silva, P.; Campanari, S.; Macchi, E. A rolling-horizon optimization algorithm for the long term operational scheduling of cogeneration systems. Energy 2017. [Google Scholar] [CrossRef]
  21. Brijs, T.; Geth, F.; Siddiqui, S.; Hobbs, B.F.; Belmans, R. Price-based unit commitment electricity storage arbitrage with piecewise linear price-effects. J. Energy Storage 2016, 7, 52–62. [Google Scholar] [CrossRef]
  22. Lopion, P.; Markewitz, P.; Robinius, M.; Stolten, D. A review of current challenges and trends in energy systems modeling. Renew. Sustain. Energy Rev. 2018, 96, 156–166. [Google Scholar] [CrossRef]
  23. Breithaupt, T.; Leveringhaus, T.; Hofmann, L. Heuristic solution of a joint electric control reserve and wholesale market model. Int. J. Model. Simul. Sci. Comput. 2018, 1940001. [Google Scholar] [CrossRef]
  24. Li, B.; Roche, R.; Miraoui, A. Microgrid sizing with combined evolutionary algorithm and MILP unit commitment. Appl. Energy 2017, 188, 547–562. [Google Scholar] [CrossRef]
  25. Marquant, J.F.; Evins, R.; Carmeliet, J. Reducing Computation Time with a Rolling Horizon Approach Applied to a MILP Formulation of Multiple Urban Energy Hub System. Procedia Comput. Sci. 2015, 51, 2137–2146. [Google Scholar] [CrossRef]
  26. Wolf, D.; Witt, M.; Bruckner, T. Auswirkung der fluktuierenden Stromeinspeisung aus Windenergie auf die CO2-Emissionen fossil befeuerter Kraftwerke. In Energiesysteme der Zukunft: Technologien und Investitionen Zwischen Markt und Regulierung 5; Internationale Energiewirtschaftstagung: Wien, Austria, 2007. [Google Scholar]
  27. Kehlhofer, R. Combined-Cycle Gas & Steam Turbine Power Plants, 3rd ed.; PennWell: Tulsa, OK, USA, 2009; ISBN 9781593701680. [Google Scholar]
  28. IBM ILOG CPLEX Optimization Studio; IBM: New York, NY, USA, 2017.
  29. Hultberg, T.H. FLOPC++ An Algebraic Modeling Language Embedded in C++. In Operations Research Proceedings 2006; Waldmann, K.-H., Stocker, U.M., Eds.; Springer: Berlin/Heidelberg, Germany, 2007; pp. 187–190. ISBN 978-3-540-69994-1. [Google Scholar]
  30. 50 Hertz. Archiv Photovoltaik. Available online: http://www.50hertz.com/de/Kennzahlen/Windenergie/Archiv-Photovoltaik (accessed on 12 July 2018).
  31. 50 Hertz. Archiv Windenergie. Available online: http://www.50hertz.com/de/Kennzahlen/Windenergie/Archiv-Windenergie (accessed on 12 July 2018).
  32. Amprion. Einspeisung von Photovoltaik. Available online: https://www.amprion.net/Netzkennzahlen/Photovoltaikeinspeisung/index-3.html (accessed on 12 July 2018).
  33. Amprion. Einspeisung von Windenergie. Available online: https://www.amprion.net/Netzkennzahlen/Windenergieeinspeisung/index-3.html (accessed on 12 July 2018).
  34. Tennet. Tatsächliche und Prognostizierte Windenergieeinspeisung. Available online: http://www.tennettso.de/site/de/Transparenz/veroeffentlichungen/netzkennzahlen/tatsaechliche-und-prognostizierte-windenergieeinspeisung (accessed on 12 July 2018).
  35. TransnetBW. Erneuerbare Energien: Einspeisung Fotovoltaik. Available online: https://www.transnetbw.de/de/transparenz/marktdaten/kennzahlen (accessed on 12 July 2018).
  36. TransnetBW. Erneuerbare Energien: Einspeisung Windenergie. Available online: https://www.transnetbw.de/de/transparenz/marktdaten/kennzahlen (accessed on 12 July 2018).
  37. Tennet. Tatsächliche und Prognostizierte Solarenergieeinspeisung. Available online: http://www.tennettso.de/site/de/Transparenz/veroeffentlichungen/netzkennzahlen/tatsaechliche-und-prognostizierte-solarenergieeinspeisung_land (accessed on 12 July 2018).
  38. Netztransparenz. EEG-Anlagenstammdaten 2016. Available online: https://www.netztransparenz.de/EEG/Anlagenstammdaten (accessed on 12 July 2018).
  39. Blyr, A.; Sigala, C.; Amatucci, G.; Guyomard, D.; Chabre, Y.; Tarascon, J.-M. Self-Discharge of LiMn2 O 4/C Li-Ion Cells in Their Discharged State Understanding by Means of Three-Electrode Measurements. J. Electrochem. Soc. 1998, 145, 194–209. [Google Scholar] [CrossRef]
  40. Gurobi Optimization, L.L.C. Gurobi Optimizer Reference Manual. Available online: http://www.gurobi.com (accessed on 24 August 2018).
  41. Kanngießer, A. Entwicklung Eines Generischen Modells zur Einsatzoptimierung von Energiespeichern für die Techno-Ökonomische Bewertung Stationärer Speicheranwendungen. Ph.D. Dissertation, Technische Univsersität Dortmund, Dortmund, Germany, 29 April 2013. [Google Scholar]
  42. Baldauf, M.; Förstner, J.; Klink, S.; Reinhardt, T.; Schraff, C.; Seifert, A.; Stephan, K. Kurze Beschreibung des Lokal-Modells Kürzestfrist COSMO-DE (LMK) und seiner Datenbanken auf dem Datenserver des DWD No. 2.3. 2014. Available online: https://www.dwd.de/SharedDocs/downloads/DE/modelldokumentationen/nwv/cosmo_de/cosmo_de_dbbeschr_version_2_3_201406.pdf?__blob=publicationFile&v=5 (accessed on 24 July 2018).
  43. Reinert, D.; Prill, F.; Frank, H.; Denhard, M.; Zängl, G. Database Reference Manual for ICON and ICON-EPS; No. 1.2.3; Deutscher Wetterdienst: Dortmund, Germany, 2018. [Google Scholar]
  44. Voll, P.; Klaffke, C.; Hennen, M.; Bardow, A. Automated superstructure-based synthesis and optimization of distributed energy supply systems. Energy 2013, 50, 374–388. [Google Scholar] [CrossRef]
  45. Kopp, M.; Coleman, D.; Stiller, C.; Scheffer, K.; Aichinger, J.; Scheppat, B. Energiepark Mainz: Technical and economic analysis of the worldwide largest Power-to-Gas plant with PEM electrolysis. Int. J. Hydrogen Energy 2017, 42, 13311–13320. [Google Scholar] [CrossRef]
  46. Fuchs, G.; Lunz, B.; Leuthold, M.; Sauer, D.U. Technology Overview on Electricity Storage. Overview on the Potential and on the Deployment Perspectives of Electricity Storage Technologies; ISEA: Arlington, VA, USA, 2012. [Google Scholar]
  47. Peniche Garcia, R. Analysis of Renewable Energy Integration Options in Urban Energy Systems with Centralized Energy Parks. Ph.D. Dissertation, Hamburg University of Technology, Hamburg, Germany, 2017. [Google Scholar]
  48. Elsner, P.; Sauer, D.U. Energiespeicher. Technologiesteckbriefe zur Analyse „Flexibilitätskonzepte für die Stromversorgung 2050. “Technologiesteckbriefe zur Analyse” Flexibilitätskonzepte für die Stromversorgung 2050. 2015. Available online: https://energiesysteme-zukunft.de/publikationen/materialien/technologiesteckbriefe-zur-analyse-flexibilitaetskonzepte-fuer-die-stromversorgung-2050/ (accessed on 23 July 2018).
Figure 1. Relative efficiency curves for combined cycle gas turbines (solid line) and open cycle gas turbines (dashed line) according to [27].
Figure 1. Relative efficiency curves for combined cycle gas turbines (solid line) and open cycle gas turbines (dashed line) according to [27].
Energies 12 01003 g001
Figure 2. Simplified schematic program flow of CharL.
Figure 2. Simplified schematic program flow of CharL.
Energies 12 01003 g002
Figure 3. Simplified diagram of a rolling horizon with a ratio of period length to interval length of 1/3.
Figure 3. Simplified diagram of a rolling horizon with a ratio of period length to interval length of 1/3.
Energies 12 01003 g003
Figure 4. Input data for the model: solar (yellow) and wind (blue) power generation. The capacity factor represents the ratio of generated power to installed capacity.
Figure 4. Input data for the model: solar (yellow) and wind (blue) power generation. The capacity factor represents the ratio of generated power to installed capacity.
Energies 12 01003 g004
Figure 5. Efficiency curves as used in algorithm (dashed lines) and their linearised counter parts (dotted lines) as used for MILP for CCGT (which is also long-term storage discharge) (circles), and long-term storage charge unit (squares). The linearised counterparts are generated using Equation (15) with the parameters of Table 3.
Figure 5. Efficiency curves as used in algorithm (dashed lines) and their linearised counter parts (dotted lines) as used for MILP for CCGT (which is also long-term storage discharge) (circles), and long-term storage charge unit (squares). The linearised counterparts are generated using Equation (15) with the parameters of Table 3.
Energies 12 01003 g005
Figure 6. Specific CO2 emissions (circles) and share of storages of supplied demand (squares) for the heuristic approach (solid lines) and a UCO with a 2 days interval (dotted lines).
Figure 6. Specific CO2 emissions (circles) and share of storages of supplied demand (squares) for the heuristic approach (solid lines) and a UCO with a 2 days interval (dotted lines).
Energies 12 01003 g006
Figure 7. Absolute (circles) and relative (squares) reductions in CO2 emissions as a result of UCO in comparison to UCH for scenario a (2 days interval).
Figure 7. Absolute (circles) and relative (squares) reductions in CO2 emissions as a result of UCO in comparison to UCH for scenario a (2 days interval).
Energies 12 01003 g007
Figure 8. Reduction of specific CO2 emissions through UCO compared to heuristic approach for different installed VRE powers (expressed in MWeach) and interval lengths.
Figure 8. Reduction of specific CO2 emissions through UCO compared to heuristic approach for different installed VRE powers (expressed in MWeach) and interval lengths.
Energies 12 01003 g008
Figure 9. Time series excerpt of power of storages according to heuristic approach and UCO with 5 days interval on 17 January. 3000 MWeach of VRE are installed. Short-term storage power of UCH is not shown, since it is zero.
Figure 9. Time series excerpt of power of storages according to heuristic approach and UCO with 5 days interval on 17 January. 3000 MWeach of VRE are installed. Short-term storage power of UCH is not shown, since it is zero.
Energies 12 01003 g009
Figure 10. Time series of the storage level of a long-term storage for MILP result (dashed) and after control calculation (solid) for 3000 MWeach VRE power, using 5 days interval (scenario c, blue) and 365 days interval (scenario f, green).
Figure 10. Time series of the storage level of a long-term storage for MILP result (dashed) and after control calculation (solid) for 3000 MWeach VRE power, using 5 days interval (scenario c, blue) and 365 days interval (scenario f, green).
Energies 12 01003 g010
Table 1. Temporal configurations of the scenarios and abort criteria for MILP.
Table 1. Temporal configurations of the scenarios and abort criteria for MILP.
ScenarioHardware UsedInterval LengthPeriod LengthAbort Criteria
(see Table A1)Max. Time in sRel. Gap
aPC 12 days1 day12001.0 × 10−6
bPC 23 days1 day24001.0 × 10−6
cPC 35 days2 days36001.0 × 10−5
dPC 414 days5 days10,8001.0 × 10−5
ePC 330 days10 days21,6001.5 × 10−5
fPC 1365 days365 days100,0001.0 × 10−4
Table 2. Initial storage SOC depending on rated power of installed renewables.
Table 2. Initial storage SOC depending on rated power of installed renewables.
VRE Configuration in MWeachInitial SOC of Long-Term Storage in -
1000 … 22500.000
25000.056
27500.182
30000.417
32500.911
3500 … 50001.000
Table 3. Parameters a and b of Equation (15) for MILP.
Table 3. Parameters a and b of Equation (15) for MILP.
Unitab
CCGT/Long-term discharge0.6966390.2044030
Long-term charge0.6702190.0283414
Short-term (dis)charge0.9205000

Share and Cite

MDPI and ACS Style

Erichsen, G.; Zimmermann, T.; Kather, A. Effect of Different Interval Lengths in a Rolling Horizon MILP Unit Commitment with Non-Linear Control Model for a Small Energy System. Energies 2019, 12, 1003. https://doi.org/10.3390/en12061003

AMA Style

Erichsen G, Zimmermann T, Kather A. Effect of Different Interval Lengths in a Rolling Horizon MILP Unit Commitment with Non-Linear Control Model for a Small Energy System. Energies. 2019; 12(6):1003. https://doi.org/10.3390/en12061003

Chicago/Turabian Style

Erichsen, Gerrit, Tobias Zimmermann, and Alfons Kather. 2019. "Effect of Different Interval Lengths in a Rolling Horizon MILP Unit Commitment with Non-Linear Control Model for a Small Energy System" Energies 12, no. 6: 1003. https://doi.org/10.3390/en12061003

APA Style

Erichsen, G., Zimmermann, T., & Kather, A. (2019). Effect of Different Interval Lengths in a Rolling Horizon MILP Unit Commitment with Non-Linear Control Model for a Small Energy System. Energies, 12(6), 1003. https://doi.org/10.3390/en12061003

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop