Next Article in Journal
A Lightweight Secure and Energy-Efficient Fog-Based Routing Protocol for Constraint Sensors Network
Next Article in Special Issue
Modelling Bidding Behaviour on German Photovoltaic Auctions
Previous Article in Journal
Core Temperature Estimation for a Lithium ion 18650 Cell
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Day-Ahead Market Modelling of Large-Scale Highly-Renewable Multi-Energy Systems: Analysis of the North Sea Region towards 2050

1
Department of Management, Technical University of Denmark (DTU), 2800 Lyngby, Denmark
2
Department of Wind Energy, Technical University of Denmark (DTU), 4000 Roskilde, Denmark
3
Department of Electrical Power Engineering and Mechatronics, Tallin University of Technology (TalTech), 12616 Tallinn, Estonia
*
Author to whom correspondence should be addressed.
Energies 2021, 14(1), 88; https://doi.org/10.3390/en14010088
Submission received: 9 November 2020 / Revised: 15 December 2020 / Accepted: 16 December 2020 / Published: 25 December 2020

Abstract

:
This paper proposes a mathematical model in order to simulate Day-ahead markets of large-scale multi-energy systems with a high share of renewable energy. Furthermore, it analyses the importance of including unit commitment when performing such analysis. The results of the case study, which is performed for the North Sea region, show the influence of massive renewable penetration in the energy sector and increasing electrification of the district heating sector towards 2050, and how this impacts the role of other energy sources, such as thermal and hydro. The penetration of wind and solar is likely to challenge the need for balancing in the system as well as the profitability of thermal units. The degree of influence of the unit commitment approach is found to be dependent on the configuration of the energy system. Overall, including unit commitment constraints with integer variables leads to more realistic behaviour of the units, at the cost of considerably increasing the computational time. Relaxing integer variables significantly reduces the computational time, without highly compromising the accuracy of the results. The proposed model, together with the insights from the study case, can be especially useful for system operators for optimal operational planning.

Graphical Abstract

1. Introduction

Energy systems, including all energy vectors, such as heating, transportation, and agriculture, are converting to electricity-based energy usage, due to climate change and environmental concerns. The European Commission has the vision of decarbonising the whole energy system by 2050 [1]. Denmark, in this direction, has the ambition to completely phase-out coal by 2030 [2]. Renewable based electricity generators, such as hydro, wind, solar photovoltaic (PV), or biomass, are replacing the carbon-based generators. Many of these renewable energy sources are inherently variable in nature, such as wind, solar, or micro-hydro. Consequently, increasing the share of such variable renewable energy (VRE) sources in electricity systems increases the variability and uncertainty in the full energy system.
Maintaining a stable and secure operation in the electricity system with large share of VRE can be very challenging for the power system operators. Major challenges that are involved in operational planning are the estimation of operational reserves [3], or determining the ramp requirements and flexibility for the generators. This information is then used in order to mitigate the impact of variability in the electric power system.
Operational planning should include the co-optimisation of all the sectors to avoid infeasibilities and sub-optimal solutions. For example, the estimation of reserves from combined heat and power (CHP) units while only performing the optimisation of electrical power systems can create infeasibilities for the units in real-time due to heating constraints and, thereby, challenging the security of the system. Exploiting the synergies of multi-energy analysis can also include, for instance, planning maintenance of the units when considering the needs of the different parts of the energy sector.
The operation of the system generally takes place in real life through different energy markets. Examples of these market for electricity are day-ahead (DA), intra-day, or balancing markets. This paper focuses on simulating the operation of the DA market.
It is generally relevant to consider the unit commitment (UC) problem when simulating the operation of the DA market [4]. Different UC modelling approaches have been researched for many years [4,5,6]. Mixed integer programming (MIP) is a widely used methodology for UC. However, MIP based methods are computationally expensive. Lagrangian relaxation offers saving in computational time without compromising the accuracy [7].
The simulation of DA operation of large-scale multi-energy systems when considering UC constraints is computationally challenging. In the literature, either such a problem has been handled for a small system (for example, a six-bus system in [8]) or a compromise has been made in terms of simulation horizon (24 h in [9]), technologies involved (such as only for electrical systems), or scenario years analysed. Recently, UC based studies have also been applied for multi-carrier energy systems, but, again, the application has been limited to small test systems [10]. Performing these operational planning studies for short term can lead to infeasible solutions and it is otherwise sub-optimal for reasons, such as being unable to consider long-term constraints, for example, yearly schedule of hydro reservoirs. The integration of VRE has also been limited in the sense of modelling the details of weather dependencies.
Table 1 shows an analysis of the available literature concerning relevant features for the scope of this paper. It can be observed that there is no available literature that can handle all of the features.
The novelty of this paper involves developing a methodology for integrating all of the features that are mentioned in Table 1 to simulate the DA market and using the developed methodology to provide analysis and recommendations for future scenarios of very large energy systems. The efficacy of the developed methodology is applied for case studies of large-scale energy systems of North Sea countries for 2020, 2030, and 2050 energy scenarios, including a detailed representation of renewable variations. The computational and practical challenges in modelling and implementation for such a large system are discussed. The computational cost and accuracy of the results for different UC modelling approaches are also compared. Even though the existing literature has deeply investigated the UC problem, introducing such analysis for the context of this study (large scale, different scenario years, multi-energy system) can be relevant in understanding the advantages and disadvantages of the UC approach when simulating the DA market.
The paper is structured, as follows. Section 2 explains the mathematical model-based methodology that was applied to model the DA market operation. Because the scale of the problem is large in terms of technologies, geography, and time period, special considerations need to be taken in order to reduce computational complexity, as described in Section 3. Section 4 presents the case study. Section 5 shows the results and discusses the limitations of the study, while Section 6 summarises the conclusions. The symbols used are shown in the nomenclature section at the end of the article.

2. Mathematical Modelling

The methodology used in this paper to simulate the DA market can be divided into four stages: DA optimisation (Section 2.1), VRE simulations (Section 2.2), storage and planned maintenance optimisation (Section 2.3), and stochastic outage simulations (Section 2.4). The stages are linked, as shown in the flow chart of Figure 1. Section 2.5 presents the sensitivity cases studied in this paper, which focus on the UC modelling approach.
The optimisations and simulations, except for the VRE simulations, are performed with the energy system model Balmorel [24], an energy system tool, deterministic, open source [25], with a bottom-up approach. It has been traditionally used in order to model the electricity and district heating sectors, although it is being developed to increase its capabilities and include more sectors [26].
The temporal representation used is composed of years y Y , which are composed of seasons s S (in this paper days), which are disaggregated into time steps t T (in this paper hours).
The geographical representation used is composed of countries, which are composed of regions r R , which are composed of areas a A . The regions in the model represent copperplates for the transmission of electricity. The areas in the model represent copperplates for the transmission of heat.

2.1. Day-Ahead Optimisation

The DA optimisation of the energy markets is performed on a daily basis in order to model the behaviour of the spot market of electricity. This means that consecutive dispatch optimisations of 24 h are performed. The results are linked from day-to-day, meaning that operational decisions made in previous days can have a limiting effect on the operation on the following days, depending on the flexibility of the units g GGG .
The storage content at the beginning of each day for the DA optimisation, as well as the planned maintenance, is fixed from the storage and planned maintenance optimisation. Additionally, the availability of the units within the day is also affected by stochastic outage, which, together with the planned maintenance, determine the final availability of the units to participate in the markets. A detailed mathematical formulation of the problem can be found in Appendix A of [27]. For the sake of space limitation, only a few equations are shown here.

2.1.1. Objective Function

The objective function includes all the operational costs of the studied system in the studied time steps and during the solving phase the value of the objective function is minimised (Equation (1)). Furthermore, the time steps considered correspond to one day. In this paper, the costs have been aggregated into variable operational and maintenance costs ( c y v o m ), and emission tax costs ( c y e m i ). The fixed operational and maintenance costs ( C y f o m ), which depend on the installed capacity, are also included in the objective function as a parameter. The disaggregation of these costs can be found in Equations (A2) to (A6) of [27].
min y Y c y v o m + c y e m i + C y f o m

2.1.2. System Constraints

There are two main system constraints that need to be fulfilled in every time step. One of them is the electricity balance in each region. The electricity production ( p g , a , y , s , t e l ) of electricity generation units ( EL ) minus the electricity storage loading ( s t o l g , a , y , s , t ) of storage units ( STO ) and minus the electricity demand resulting from heat production ( p g , a , y , s , t h ) of electricity-to-heat units ( PTOH ) in all areas that belong to the region ( AR ) must equal the inelastic electricity demand of the region ( D r , y , s , t e l ) plus the export flows ( x r , r , y , s , t ) minus import flows ( x r , r , y , s , t ) with other regions. The transmission losses ( x l o s s ) are allocated to electricity import:
a AR g EL p g , a , y , s , t e l g STO , g EL s t o l g , a , y , s , t g PTOH p g , a , y , s , t h η g , s , t = D r , y , s , t e l + r R x r , r , y , s , t x r , r , y , s , t · ( 1 x l o s s ) r R , y Y , s S , t T
The other system constraint is the heat balance in each area, where the heat production ( p g , a , y , s , t h ) of the units generating heat ( HEAT ) minus the endogenous loading of heat storage ( s t o l g , a , y , s , t ) must equal the inelastic heat demand ( D a , y , s , t h ):
g HEAT p g , a , y , s , t h g STO , g HEAT s t o l g , a , y , s , t = D a , y , s , t h a A , y Y , s S , t T
Ancillary services in the electricity sector were not included to simplify the problem.

2.1.3. Technological Constraints

The generation technologies in the model are split into three categories: dispatchable generation units, pure storage units, and VRE units. Dispatchable technologies include electricity-only, combined heat and power, heat pumps, boilers, geothermal, and hydroelectric power with reservoirs and seasonal inflow. These units are subject to several technical constraints. The operation of these technologies depends on the available units in each time step and the specific technical characteristics of each type. The equations that are relevant for each of these technologies are described from (A29) to (A48) in the report [27]. The storage units modelled without seasonal inflow that can only modify their storage level based on their active loading and unloading are defined as pure storage units. These technologies are hydro pumping, electric batteries, pit heat storage, or hot water tanks. The relevant equations for these units are Equations (A49) to (A71) in [27]. The non-dispatchable technologies included are solar heating, solar PV, wind offshore and onshore, and hydro-run-of-river. The equations limiting their production are (A72) and (A73) in [27]. Furthermore, electricity trade is allowed between regions, which are assumed to be copper plates, and it is limited to the available transmission capacity between the regions in each time step. The relevant equation limiting their operation is (A74) in [27]. Heat trade between areas is not allowed.

2.2. Modelling of Renewable Generation Including Fluctuations

The CorRES model [28] is used for the simulation of the VRE generation time series being used as inputs for Balmorel. The CorRES model is based on data from the weather research and forecasting model [29]. In addition to modelling the spatiotemporal dependencies in solar PV and wind generation, CorRES allows for the modelling of VRE technology development impacts on the VRE time series. For the presented case study, the assumed VRE technology developments towards 2050 are described in [30]. The resolution of the simulated VRE generation time series is hourly; the data are aggregated to the regions that were used in Balmorel.

2.3. Storage and Planned Maintenance Optimisation

When simulating the DA market, it is important to capture that some of the decisions that market participants take into account future market price expectations, rather than just planning for the next 24 h. The planned maintenance and use of storage are part of these long-term decisions. In this paper, these two decisions are obtained by performing a full year dispatch optimisation. Planned maintenance and storage content at the beginning of each season are saved and forced in the DA optimisation. Planned maintenance is also used in the stochastic outage simulation in order to calculate the available capacity for production that can suffer an unexpected outage.
The formulation is similar to the one shown in Section 2.1, with a few exceptions. In this optimisation, the time steps considered correspond to one year. All of the equations can be found in [27]. In this paper, only some equations are shown.

2.3.1. Available Units

The number of available units available for generation ( n g , a , y , s , t a v , o n ) is constrained by the total number of units (installed fuel capacity of a unit type ( F C g , a , y ) divided by the size of each unit type ( U S g g e n )), and the amount of units on planned maintenance ( n g , a , y , s n a v , p m ). The number of units on planned maintenance influence the availability factor of the units when simulating the DA market.
F C g , a , y U S g g e n n g , a , y , s n a v , p m n g , a , y , s , t a v , o n g GGG , a A , y Y , s S , t T
The maximum number of units on maintenance is limited by the total installed number of units:
F C g , a , y U S g g e n n g , a , y , s n a v , p m g GGG , a A , y Y , s S

2.3.2. Yearly Maintenance Requirement

The minimum maintenance time of each technology ( M M T g ) needs to be respected during the year:
s S n g , a , y , s n a v , p m = M M T g · F C g , a , y U S g g e n g GGG , a A , y Y

2.3.3. Uninterrupted Maintenance

Units starting up maintenance ( n g , a , y , s n a v , p m , s u ) must be on maintenance during the minimum maintenance time without interruption, i.e., the maintenance must take place in consecutive seasons:
s = 1 M M T g n g , a , y , s s n a v , p m , s u n g , a , y , s n a v , p m g GGG , a A , y Y , s S

2.3.4. Logical Conditions

The number of units on maintenance depends on the units starting or stopping maintenance ( n n a v , p m , s d ):
n g , a , y , s n a v , p m n g , a , y , s 1 n a v , p m = n g , a , y , s n a v , p m , s u n g , a , y , s n a v , p m , s d g GGG , a A , y Y , s S
In order to reflect the discrete nature of the generation units that are part of the energy system, the following variables are restricted to be positive integer variables Z + . However, this constraint is relaxed due to the computational complexity in the optimisation, as mentioned in Section 3.2.
n g , a , y , s , t a v , o n , n g , a , y , s n a v , p m , n g , a , y , s n a v , p m , s u , n g , a , y , s n a v , p m , s d Z + g GGG , a A , y Y , s S , t T

2.4. Stochastic Outage Simulations

Unexpected operational problems can lead to making units being unavailable until the problem is fixed, which can influence market prices. Hence, it is relevant to capture these occurrences. Using, as input parameters, the size of a single unit, the total capacity of a unit type in each area, the planned maintenance (if previously calculated), and the probability of an unexpected outage, for each time step and unit in the system, Monte carlo simulations are performed in order to simulate these outages. The outcome from these simulations is then fed to the DA optimisation, and the relevant variables fixed. This approach is applied to all units except the VRE ones, since their availability is part of the time series used. The formulation can be found in Equations (A85) and (A86) of [27].

2.5. Sensitivity Cases: Unit Commitment Modelling Approaches

In order to analyse the importance of the UC modelling approach when modelling the DA market, three different sensitivity cases of UC modelling approaches are studied in the DA optimisations: (1) adding UC constraints with integer commitment variables (UC-MIP), (2) adding UC constraints with relaxed commitment variables (UC-RMIP), and (3) not adding constraints nor corresponding commitment variables (NO-UC).

3. Special Considerations for Long-Term Operational Planning of Large-Scale Energy System

3.1. Unit Commitment Assumptions

Introducing UC in the optimisation allows for an improved representation of conventional generation, at the cost of increasing considerably computational complexity due to the use integer variables. Solving a large-scale MIP problem can be intractable. In order to deal with this problem, one can either relax the integer variables, or limit the technologies modelled with integer variables. In this paper, the second approach is considered so the impact of the different optimisation approaches can be evaluated. The technologies that are modelled with UC integer variables are almost all type of fuel-based thermal plants, i.e., gas turbines, steam turbines, combined cycle turbines, and boilers. Engines were not included, since they are very fast and their size is generally much smaller than other generators, making their impact negligible. The rest of the technologies, i.e., hydro reservoirs, other storage, electricity-to-heat (P2H), and VRE, were not modelled with UC variables in order to reduce the complexity of the problem.

3.2. Simplifications in Storage and Planned Maintenance Optimisation

Optimal planning of maintenance is solved as a relaxed mixed integer problem (RMIP), including all days of the years, but with 1 every 3 h to reduce the complexity. Planned maintenance is only computed for the units that were modelled with UC (see Section 3.1).

4. Case Study: The North Sea Offshore Grid

The study case used in this paper corresponds to the offshore grid scenario presented in [30], developed as part of the project North Sea Offshore Network—Denmark. The scenario focuses on the following countries: Germany, Great Britain, Norway, Denmark, Belgium, and Netherlands. The sectors included are the electricicity and district heating sectors. The study case shows, towards 2050, a high share of VRE, transmission interconnection, and partial electrification of the district heating sector in the countries in focus. The capacity development (Figure 2) was highly influenced by the assumptions on increasing CO2 EU ETS price: 6.06, 76.70, and 130.38 €2012/ton in 2020, 2030, and 2050, respectively. More details about the scenario can be read in [30,31,32].

5. Results and Discussion

First, Section 5.1 focuses on key results that were obtained from the storage and planned maintenance optimisations. Section 5.2 presents the results from the Day-Ahead optimisations, focusing on the influence of VRE penetration and the UC modelling approach. Section 5.3 discusses the limitations of the study. Costs and prices are in €2012, which is the base year used for economic data in the model.

5.1. Storage and Planned Maintenance Optimisation

5.1.1. Planned Maintenance

Figure 3 shows the share of installed capacity under planned maintenance of district heating units burning waste in Denmark. Units burning waste tend to show high capacity factors and, hence, the scheduling of their maintenance is of relevance. The results show that most of the planned maintenance takes place during summer, which is when the district heat demand is lowest. For CHP units, by 2020, most of the maintenance takes place in July (89%), and towards 2050, a larger share of the maintenance takes place in earlier months. The fact that, in 2050, maintenance in May is 34% of the total could be linked to solar PV generation, since, in Denmark, the production of this technology is the highest in this month, favoring the use of P2H and, hence, leading to less need for CHP units to be operative. For district heating boilers, no significant difference towards 2050 is observed.

5.1.2. Planned Storage Use

The aggregated planned energy content along the year of hydro reservoirs with seasonal inflow, one of the types of storage considered in the model, is shown in Figure 4 for Norway. The purpose of this figure is to analyse the influence of VRE penetration towards 2050 on the use of this type storage during the year. The minimum value of each profile has been subtracted. The results show that the reservoirs are mainly filled during the summer and discharged during winter. The maximum energy content is higher in 2020 than in 2030 and 2050. This is a result of storing less energy in the reservoir during the year, which is linked to using hydro energy for balancing VRE. Furthermore, the evolution of the content in 2020 is smoother than in 2030 and 2050, which is influenced by the volatility and seasonality of VRE generation. These results strengthen the importance of performing full year optimisations towards 2050.

5.2. Day-Ahead Optimisation

5.2.1. Annual Production

Figure 5 depicts the aggregated generation of electricity and heat per year, scenario, and technology in the countries in focus for different UC approaches. The penetration of VRE in the electricity sector towards 2050 is remarkable, at the expense of decreasing the use of thermal technologies. The results show that the share of CO2 free generation increases from 64% in 2020 to 91%, which is linked to the assumed VRE penetration. On the heating side, the generation of P2H units increases towards 2050 at the expense of thermal units. Introducing UC constraints slightly increases the aggregated production of thermal power units and reduces VRE generation.
The average capacity factor of thermal power units for electricity with/without heat generation is shown in Figure 6 for the UC-MIP approach. The difference with NO-UC and UC-RMIP is negligible. The electricity capacity factor is calculated with rated electricity generation capacity, whereas the heat one uses the rated heat generation capacity. The results show that, towards 2050, the average electricity capacity factor decreases for CHP back pressure units, CHP extraction units, and electricity-only units, especially for the last two types. The heat capacity factor also tends to decrease. These results suggest that the massive penetration of VRE might challenge the profitability of thermal units towards 2050.

5.2.2. System Costs

Table 2 shows the disaggregated operational costs development in the countries in focus per UC approach. The total operational costs increase in 2030 with respect to 2020, and then decrease in 2050 when compared to 2030. This development is linked to VRE penetration and CO2 tax assumptions. By 2030, there is still considerable fossil generation, leading to high CO2 tax costs. By 2050, even though the CO2 tax is higher than in 2030, there is much less fossil generation, which, together with further VRE penetration, leads to considerably less variable costs. Using UC-RMIP leads to much closer results to UC-MIP than using NO-UC, although the difference decreases towards 2050, again due to less fossil generation. Not including UC constraints, i.e., using NO-UC, underestimates both variable and CO2 costs, since it overestimates the flexibility of the units when ignoring relevant costs, such as start-up.

5.2.3. Hourly Electricity Balance

Figure 7 depicts the hourly electricity balance for four representative days for the years 2020, 2030, and 2050, for Great Britain and with the UC-MIP approach. The graph includes aggregated generation and demand per type, as well as electricity prices. The prices are derived from the dual variable of the electricity balance equation (Equation (2)). The penetration of VRE replaces most thermal generation towards 2050. Such a large penetration of VRE energy highlights the need for proper planning of balancing resources towards 2050, since forecast errors could challenge the correct operation of the system.

5.2.4. Electricity Prices

Figure 8 shows the influence of the UC approach in the cumulative probability curves for the electricity prices in each year, for one of the modelled regions of Denmark, i.e. DK1. The influence of the optimisation approach decreases towards 2050, which is explained with the decrease of thermal power capacity and use towards 2050, since they are the ones that are affected by the different UC approaches. UC-MIP leads to overall higher prices in all price-range due to forcing discrete block sizes to be on/off for thermal units. The number of hours with very low prices (where VRE curtailment sets the price) are higher with UC-RMIP, and especially with UC-MIP. These low prices generally correspond to those hours where it is cheaper (or the only feasible way) in order to increase the use of more expensive generation units, rather than starting/shutting them. In the mid-range prices, the prices for UC-RMIP tend to be the lowest, which can be explained in a similar way to VRE curtailment, but, instead, what it is being “curtailed” is the next available cheaper generator. On the high-range prices, NO-UC underestimates high prices, which can be explained with the non-consideration of restrictive constraints, like minimum-on/off time, minimum production, and ramping. The order of magnitude of the prices is highly influenced by the CO2 tax assumption development (Section 4).

5.2.5. Curtailment

The influence of the UC approach on curtailment for different technology types in the countries in focus per year are shown in Table 3. Not considering UC costs leads to less curtailment. Wind offshore is curtailed more often than onshore due to operation costs assumptions, i.e., variable costs for offshore are more expensive than for onshore. The impact of not considering UC costs on curtailment increases towards 2050, when there is higher VRE penetration.

5.2.6. Average Revenue of Wind and Solar Pv Units

Figure 9 shows the influence of the UC approach on the variability of the average revenue per energy unit sold in each region and year from the operation in the DA market for wind and solar PV units. The results show a decrease in the variability of the average revenue across countries towards 2050, mainly due to grid expansion. This highlights the importance of large-scale energy system analysis. Solar PV’s average revenue is higher in 2020 and 2030 than wind unit’s, but it is slightly smaller in 2050. With a UC-MIP approach, the average revenues are slightly higher for both technologies. The yearly levels of average revenue are directly linked to the DA price development.

5.2.7. Electricity-Only Thermal Plant Operation

Table 4 shows the influence of the UC approach in the utilization factor of electricity-only thermal generation units using natural gas or nuclear fuels. The installed capacity of these units is also shown. The results show a decrease in utilization for all of the considered units towards 2050. Units burning natural gas half their capacity factor with respect to 2020 in 2030 and stay almost constant by 2050. Nuclear power capacity factor also decreases considerably towards 2050. The UC approach impacts the utilization of nuclear power by 2050, overestimating its flexibility when using NO-UC. These results are influenced by the decrease of nuclear power in the system towards 2050, and the increasing penetration of VRE. This result challenges the profitability of these thermal units towards 2050.
The impact of the UC modelling approach on hourly nuclear operation, which is highly affected by UC constraints, is shown in Figure 10, which shows the hourly operation of nuclear power plants for four consecutive days in Great Britain. The figure depicts that, towards 2050, with higher VRE penetration, unless the UC costs and constraints are considered, nuclear units start up and shut down with high frequency. This operation of nuclear plants might be extremely challenging, and perhaps unrealistic. These results strengthen the importance of multiple scenario year analysis as well as the UC modelling approach used.

5.2.8. Computational Time

Table 5 shows the influence of the UC modelling approach on the computational time of the DA optimisation of each full year. The computational time that is required to solve the UC-MIP case is much higher than the others. The computational time decreases significantly towards 2050, since the number of thermal power plants is considerably reduced in this scenario. Using UC-RMIP leads to closer computational times to NO-UC than to UC-MIP.

5.3. Limitations of the Study

The simplifications undertaken to reduce the complexity of the problem limit the findings of this study.
The flexibility of technologies for which UC constraints were not applied could have been overestimated, especially for hydro power units.
Further work should include adequacy analysis and/or ancillary services requirements in the simulations. Both of these aspects can be highly relevant for the correct operation of the system, especially towards 2050, with more VRE penetration. Furthermore, future research should also include stronger sector coupling, since it can considerably influence the generation of VRE technologies and the need for flexibility ([23,26,33]).
The full-year foresight assumed for the planned maintenance and storage optimisation can be unrealistic, due to the high uncertainty on, e.g., weather. Performing the analysis with several weather years and/or performing stochastic optimisation would help to understand the role of uncertainty in the results. In this paper, for simplification, a unique weather year was used (2012), and stochastic features are mainly modelled through unplanned outages.

6. Conclusions

This paper proposes a mathematical model to simulate Day-ahead markets of large-scale multi-energy systems with a high penetration of renewable energy towards 2050. Furthermore, it analyses the influence of UC modelling on the results. The results highlight the importance of analysing multiple scenario years with long time series, including several sectors, as well as the value of not restricting the analysis to the small scale.
The results for the studied case of the North Sea region show how the penetration of VRE towards 2050 challenges thermal units’ traditional operation in the electricity and district heating sectors towards 2050. Furthermore, the high penetration of wind and solar is likely to challenge the need for balancing in the system and, hence, flexibility will be very important towards 2050.
The influence of the UC approach is found to be dependent on the scenario year. Generally, including UC constraints with integer variables leads to more realistic behaviour of the units, at the cost of considerably increasing the computational time. Relaxing integer variables significantly reduces the computational time, but medium-level prices are underestimated. Not including UC constraints leads to an underestimation of costs, VRE curtailment, VRE’s average revenue per energy unit sold, as well as price volatility. It also overestimates the flexibility of the thermal units. Hence, depending on the purpose of the analysis, it is recommended to think carefully on which UC modelling approach to use and acknowledge the limitations. When the focus is on prices and revenues, using UC constraints with integer variables is preferable; otherwise, relaxing the integer variables is encouraged.
Even though the study case is performed on the North Sea region, the findings of this paper are valid to other regions of the world, where high VRE penetration is expected.
The proposed model, together with the insights that were obtained from the study case, can be specially useful for system operators, who can use this model to perform operational planning studies towards 2050.

Author Contributions

Conceptualization, J.G.-B., K.D., H.K. and M.J.K.; Data curation, J.G.-B.; Formal analysis, J.G.-B.; Investigation, J.G.-B., K.D., H.K. and M.J.K.; Methodology, J.G.-B. and H.K.; Resources, J.G.-B.; Software, J.G.-B. and H.K.; Supervision, J.G.-B., K.D. and M.J.K. Validation, J.G.-B.; Visualization, J.G.-B.; Writing—original draft, J.G.-B.; Writing—review & editing, J.G.-B., K.D., H.K. and M.J.K. All authors have read and agreed to the published version of the manuscript.

Funding

This paper is a part of Flex4RES and NSON-DK projects, which are funded by the Nordic Energy Research (grant 76084), and Danish Energy Agency, EUDP (grant 64018-0032; previously ForskEL), respectively. Matti Koivisto acknowledges support from the PSfuture project (La Cour Fellowship, DTU Wind Energy).

Acknowledgments

The authors thank Hans Ravn for his contributions.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Acronyms
V R E Variable Renewable Energy
V R E Variable Renewable Energy
C H P Combined Heat and Power
P 2 H Power-to-Heat
D A Day Ahead
U C Unit Commitment
M I P Mixed Integer Linear Programming
R M I P Relaxed Mixed Integer Linear Programming
Sets
g GGG Generation and pure storage units
y Y Years
s S Seasons
t T Time steps
a A Areas
r R Regions
a r AR Areas in regions
Subsets
STO GGG Pure storage units
PTOH GGG Electricity to heat generation units
EL GGG Technologies delivering electricity to consumers
HEAT GGG Technologies delivering heat to consumers
Parameters
F C g , a , y Installed input fuel consumption capacity [MW]
x l o s s Transmission loss [-]
U S g , a , y , s , t g e n Unit size of input fuel capacity of a generation unit [MW/unit]
C f o m Fixed annual cost [€]
M M T g Minimum maintenance time [days]
D y , r , s , t e l Exogenous gross electricity consumption rate [MW]
D a , y , s , t h Exogenous gross heat consumption rate [MW]
η g , s , t Unit efficiency [-]
Positive decision variables
p g , a , y , s , t e l Net delivered electricity [MW]
p g , a , y , s , t h Net delivered heat [MW]
s t o l g , a , y , s , t Storage loading rate [MW]
x r , r , y , s , t Transmission flow [MW]
n g , a , y , s , t a v , o n Units available for generation on [-]
n g , a , y , s , t n a v , p m Units not available for generation on planned maintenance [-]
n g , a , y , s , t n a v , p m , s u Units not available for generation starting up planned maintenance [-]
n g , a , y , s , t n a v , p m , s d Units not available for generation shutting down planned maintenance [-]
c g , a , y v o m Variable operational and maintenance annual cost [€]
c g , a , y C O 2 Carbon dioxide tax annual cost [€]

References

  1. European Commission. 2050 Long-Term Strategy. Available online: https://ec.europa.eu/clima/policies/strategies/2050_en#tab-0-0 (accessed on 22 December 2020).
  2. Energi-, Forsynings-og Klimaministeriet. Denmark’s Draft Integrated National Energy and Climate Plan. 2018. Available online: https://ens.dk/sites/ens.dk/files/EnergiKlimapolitik/denmarks_draft_integrated_energy_and_climate_plan.pdf (accessed on 22 December 2020).
  3. Das, K.; Litong-Palima, M.; Maule, P.; Altin, M.; Hansen, A.D.; Sørensen, P.E.; Abildgaard, H. Adequacy of frequency reserves for high wind power generation. IET Renew. Power Gener. 2016, 11, 1286–1294. [Google Scholar] [CrossRef] [Green Version]
  4. Padhy, N.P. Unit commitment-a bibliographical survey. IEEE Trans. Power Syst. 2004, 19, 1196–1205. [Google Scholar] [CrossRef]
  5. Zhuang, F.; Galiana, F.D. Towards a more rigorous and practical unit commitment by Lagrangian relaxation. IEEE Trans. Power Syst. 1988, 3, 763–773. [Google Scholar] [CrossRef]
  6. Frangioni, A.; Gentile, C.; Lacalandra, F. Tighter approximated MILP formulations for unit commitment problems. IEEE Trans. Power Syst. 2008, 24, 105–113. [Google Scholar] [CrossRef] [Green Version]
  7. Li, T.; Shahidehpour, M. Price-based unit commitment: A case of Lagrangian relaxation versus mixed integer programming. IEEE Trans. Power Syst. 2005, 20, 2015–2025. [Google Scholar] [CrossRef] [Green Version]
  8. Zhao, C.; Jiang, R. Distributionally robust contingency-constrained unit commitment. IEEE Trans. Power Syst. 2017, 33, 94–102. [Google Scholar] [CrossRef]
  9. Shi, J.; Oren, S.S. Stochastic unit commitment with topology control recourse for power systems with large-scale renewable integration. IEEE Trans. Power Syst. 2017, 33, 3315–3324. [Google Scholar] [CrossRef]
  10. Yan, M.; Zhang, N.; Ai, X.; Shahidehpour, M.; Kang, C.; Wen, J. Robust two-stage regional-district scheduling of multi-carrier energy systems with a large penetration of wind power. IEEE Trans. Sustain. Energy 2018, 10, 1227–1239. [Google Scholar] [CrossRef]
  11. Poncelet, K.; Delarue, E.; D’haeseleer, W. Unit commitment constraints in long-term planning models: Relevance, pitfalls and the role of assumptions on flexibility. Appl. Energy 2020, 258, 113843. [Google Scholar] [CrossRef]
  12. Palmintier, B. Flexibility in generation planning: Identifying key operating constraints. In Proceedings of the 2014 Power Systems Computation Conference, Wrocław, Poland, 18–22 August 2014; pp. 1–7. [Google Scholar]
  13. Palmintier, B.; Webster, M. Impact of unit commitment constraints on generation expansion planning with renewables. In Proceedings of the 2011 IEEE Power and Energy Society General Meeting, Detroit, MI, USA, 24–28 July 2011; pp. 1–7. [Google Scholar]
  14. Pavičević, M.; Mangipinto, A.; Nijs, W.; Lombardi, F.; Kavvadias, K.; Jiménez Navarro, J.P.; Colombo, E.; Quoilin, S. The potential of sector coupling in future European energy systems: Soft linking between the Dispa-SET and JRC-EU-TIMES models. Appl. Energy 2020, 267, 115100. [Google Scholar] [CrossRef]
  15. Belderbos, A.; Valkaert, T.; Bruninx, K.; Delarue, E.; D’haeseleer, W. Facilitating renewables and power-to-gas via integrated electrical power-gas system scheduling. Appl. Energy 2020, 275, 115082. [Google Scholar] [CrossRef]
  16. Hua, B.; Baldick, R.; Wang, J. Representing Operational Flexibility in Generation Expansion Planning Through Convex Relaxation of Unit Commitment. IEEE Trans. Power Syst. 2018, 33, 2272–2281. [Google Scholar] [CrossRef]
  17. Palmintier, B.S.; Webster, M.D. Impact of Operational Flexibility on Electricity Generation Planning With Renewable and Carbon Targets. IEEE Trans. Sustain. Energy 2016, 7, 672–684. [Google Scholar] [CrossRef]
  18. Aunedi, M.; Pantaleo, A.M.; Kuriyan, K.; Strbac, G.; Shah, N. Modelling of national and local interactions between heat and electricity networks in low-carbon energy systems. Appl. Energy 2020, 276, 115522. [Google Scholar] [CrossRef]
  19. Koltsaklis, N.E.; Dagoumas, A.S. Incorporating unit commitment aspects to the European electricity markets algorithm: An optimization model for the joint clearing of energy and reserve markets. Appl. Energy 2018, 231, 235–258. [Google Scholar] [CrossRef]
  20. Koltsaklis, N.E.; Nazos, K. A stochastic MILP energy planning model incorporating power market dynamics. Appl. Energy 2017, 205, 1364–1383. [Google Scholar] [CrossRef]
  21. Clegg, S.; Mancarella, P. Integrated Electrical and Gas Network Flexibility Assessment in Low-Carbon Multi-Energy Systems. IEEE Trans. Sustain. Energy 2016, 7, 718–731. [Google Scholar] [CrossRef]
  22. Schwele, A.; Kazempour, J.; Pinson, P. Do unit commitment constraints affect generation expansion planning? A scalable stochastic model. Energy Syst. 2019, 11, 247–282. [Google Scholar] [CrossRef]
  23. Brown, T.; Schlachtberger, D.; Kies, A.; Schramm, S.; Greiner, M. Synergies of sector coupling and transmission reinforcement in a cost-optimised, highly renewable European energy system. Energy 2018, 160, 720–739. [Google Scholar] [CrossRef] [Green Version]
  24. Wiese, F.; Bramstoft, R.; Koduvere, H.; Pizarro Alonso, A.; Balyk, O.; Kirkerud, J.G.; Tveten, Å.G.; Bolkesjø, T.F.; Münster, M.; Ravn, H. Balmorel open source energy system model. Energy Strategy Rev. 2018, 20, 26–34. [Google Scholar] [CrossRef]
  25. Github. Balmorel Model. Available online: https://github.com/balmorelcommunity/Balmorel (accessed on 22 December 2020).
  26. Bermúdez, J.G.; Münster, M.; Jensen, I.G.; Koivisto, M.J.; Kirkerud, J.; Chen, Y.K. The role of sector coupling in the green transition: A least-cost energy system development in North Europe towards 2050. TechRxiv 2020. [Google Scholar] [CrossRef]
  27. Kanellas, P.; Das, K.; Gea-Bermúdez, J.; Sørensen, P. Balancing Tool Chain: Balancing and Automatic Control in North Sea Countries in 2020, 2030 and 2050. 2020. Available online: https://orbit.dtu.dk/en/publications/balancing-tool-chain-balancing-and-automatic-control-in-north-sea (accessed on 22 December 2020).
  28. Koivisto, M.; Das, K.; Guo, F.; Sørensen, P.; Nuño, E.; Cutululis, N.; Maule, P. Using time series simulation tool for assessing the effects of variable renewable energy generation on power and energy systems. WIREs Energy Environ. 2019, 8, e329. [Google Scholar] [CrossRef]
  29. Skamarock, W.; Klemp, J.; Dudhia, J.; Gill, D.; Barker, D.; Duda, M.; Huang, X.; Wang, W.; Powers, J. Description of the advanced research WRF version 3. 2008. Available online: https://opensky.ucar.edu/islandora/object/technotes:500 (accessed on 22 December 2020).
  30. Koivisto, M.; Gea-Bermúdez, J.; Sørensen, P. North Sea offshore Grid development: Combined optimization of grid and generation investments towards 2050. IET Renew. Power Gener. 2019, 14, 1259–1267. [Google Scholar] [CrossRef] [Green Version]
  31. Gea-Bermúdez, J.; Das, K.; Pade, L.; Koivisto, M.; Kanellas, P. NSON-DK Day-Ahead Market Operation Analysis in the North Sea Region towards 2050. 2019. Available online: https://orbit.dtu.dk/en/publications/nson-dk-day-ahead-market-operation-analysis-in-the-north-sea-regi (accessed on 22 December 2020).
  32. Gea-Bermúdez, J.; Pade, L.; Koivisto, M.; Ravn, H. Optimal generation and transmission development of the North Sea region: Impact of grid architecture and planning horizon. Energy 2020, 191, 116512. [Google Scholar] [CrossRef]
  33. Koivisto, M.; Gea-Bermúdez, J.; Kanellas, P.; Das, K.; Sørensen, P. North Sea region energy system towards 2050: Integrated offshore grid and sector coupling drive offshore wind power installations. Wind. Energy Sci. 2020, 5, 1705–1712. [Google Scholar] [CrossRef]
Figure 1. Flow chart of the methodology.
Figure 1. Flow chart of the methodology.
Energies 14 00088 g001
Figure 2. Installed capacity in the countries in focus (GW). The types considered are heat boilers, combined heat and power back pressure (CHP-BP), and CHP extraction (CHP-EXT), electricity-only thermal units (condensing), hydro reservoirs, hydro run-of-river (ROR), P2H, other electricity storage (PURE-EL-STO), heat storage (PURE-HEAT-STO), solar PV, wind offshore (WIND-OFF), and wind onshore (WIND-ON).
Figure 2. Installed capacity in the countries in focus (GW). The types considered are heat boilers, combined heat and power back pressure (CHP-BP), and CHP extraction (CHP-EXT), electricity-only thermal units (condensing), hydro reservoirs, hydro run-of-river (ROR), P2H, other electricity storage (PURE-EL-STO), heat storage (PURE-HEAT-STO), solar PV, wind offshore (WIND-OFF), and wind onshore (WIND-ON).
Energies 14 00088 g002
Figure 3. The share of capacity on maintenance of district heating units burning waste in Denmark.
Figure 3. The share of capacity on maintenance of district heating units burning waste in Denmark.
Energies 14 00088 g003
Figure 4. Planned energy content (TWh) along the year of hydro reservoirs in Norway at the beginning of each season, i.e., day. The minimum value of each profile has been subtracted.
Figure 4. Planned energy content (TWh) along the year of hydro reservoirs in Norway at the beginning of each season, i.e., day. The minimum value of each profile has been subtracted.
Energies 14 00088 g004
Figure 5. Production development per commodity, year, technology type, and unit commitment (UC) approach (TWh).
Figure 5. Production development per commodity, year, technology type, and unit commitment (UC) approach (TWh).
Energies 14 00088 g005
Figure 6. Average capacity factor development of thermal power units for electricity (left) and heat (right) generation using UC-MIP. The difference with NO-UC and UC-RMIP is negligible.
Figure 6. Average capacity factor development of thermal power units for electricity (left) and heat (right) generation using UC-MIP. The difference with NO-UC and UC-RMIP is negligible.
Energies 14 00088 g006
Figure 7. Hourly dispatch and electricity prices for four representative days of each year in Great Britain using UC-MIP. The demand profile corresponds to the exogenous profile. Endogenous demand from electricity storage or use of electricity-to-heat units in GB is negligible in the studied scenario.
Figure 7. Hourly dispatch and electricity prices for four representative days of each year in Great Britain using UC-MIP. The demand profile corresponds to the exogenous profile. Endogenous demand from electricity storage or use of electricity-to-heat units in GB is negligible in the studied scenario.
Energies 14 00088 g007
Figure 8. Probability distribution function of the hourly electricity price in DK1(one of the modelled regions of Denmark), for each year and UC modelling approach.
Figure 8. Probability distribution function of the hourly electricity price in DK1(one of the modelled regions of Denmark), for each year and UC modelling approach.
Energies 14 00088 g008
Figure 9. Regional variability of the average revenue per energy unit sold of wind and solar PV units in DA market towards 2050. Influence of UC modelling approach.
Figure 9. Regional variability of the average revenue per energy unit sold of wind and solar PV units in DA market towards 2050. Influence of UC modelling approach.
Energies 14 00088 g009
Figure 10. Hourly capacity factor of nuclear power in four consecutive days in Great Britain in different scenario years. The installed capacity in Great Britain of nuclear power is 8.9, 10.2 and 3.2 GW in 2020, 2030, and 2050, respectively. The impact of the UC modelling approach.
Figure 10. Hourly capacity factor of nuclear power in four consecutive days in Great Britain in different scenario years. The installed capacity in Great Britain of nuclear power is 8.9, 10.2 and 3.2 GW in 2020, 2030, and 2050, respectively. The impact of the UC modelling approach.
Energies 14 00088 g010
Table 1. Coverage of existing literature on day-ahead market modelling.
Table 1. Coverage of existing literature on day-ahead market modelling.
The Analysis Included:[5][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23]
a large-scale system
several scenario years
annual storage scheduling
full year hourly results
multi-energy markets
influence of unit commitment
planned maintenance scheduling
renewable energy sources
Table 2. Disaggregated operational costs (billion €).
Table 2. Disaggregated operational costs (billion €).
YearModelling
Approach
Fixed
Costs
Variable
Costs
Emission
Costs
Total
2020NO-UC8.212.31.822.4
UC-RMIP8.213.61.923.7
UC-MIP8.214.11.924.2
2030NO-UC9.613.56.329.4
UC-RMIP9.614.46.530.5
UC-MIP9.614.76.630.9
2050NO-UC10.68.56.125.2
UC-RMIP10.69.06.225.8
UC-MIP10.69.16.326.0
Table 3. Influence of UC approach in aggregated yearly curtailment per technology type in the countries in focus. The difference with respect to UC-MIP approach (TWh).
Table 3. Influence of UC approach in aggregated yearly curtailment per technology type in the countries in focus. The difference with respect to UC-MIP approach (TWh).
Technology TypeNO-UCUC-RMIP
202020302050202020302050
Wind offshore−0.9−4.3−13.9000
Wind onshore−0.30−0.1000
Total−1.2−4.3−14000
Table 4. Influence of UC modelling approach and year in annual capacity factor of electricity-only thermal generation units for different fuels.
Table 4. Influence of UC modelling approach and year in annual capacity factor of electricity-only thermal generation units for different fuels.
NATGASNUCLEAR
Annual Capacity
Factor (%)
Installed
Capacity (GW)
Annual Capacity
Factor (%)
Installed
Capacity (GW)
Modelling Approach202020302050202020302050202020302050202020302050
NO-UC17.17.67.664.763.750.795.285.555.123.510.73.2
UC-RMIP18.08.67.793.587.061.5
UC-MIP19.79.28.094.887.462.6
Table 5. Influence of year and UC modelling approach on the sum of computational time of all day-ahead (DA) optimisations (seconds).
Table 5. Influence of year and UC modelling approach on the sum of computational time of all day-ahead (DA) optimisations (seconds).
Modelling Approach202020302050Average
NO-UC1263165611681361
UC-RMIP11,02618,57916,22015,273
UC-MIP106,97696,43555,51486,308
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gea-Bermúdez, J.; Das, K.; Koduvere, H.; Koivisto, M.J. Day-Ahead Market Modelling of Large-Scale Highly-Renewable Multi-Energy Systems: Analysis of the North Sea Region towards 2050. Energies 2021, 14, 88. https://doi.org/10.3390/en14010088

AMA Style

Gea-Bermúdez J, Das K, Koduvere H, Koivisto MJ. Day-Ahead Market Modelling of Large-Scale Highly-Renewable Multi-Energy Systems: Analysis of the North Sea Region towards 2050. Energies. 2021; 14(1):88. https://doi.org/10.3390/en14010088

Chicago/Turabian Style

Gea-Bermúdez, Juan, Kaushik Das, Hardi Koduvere, and Matti Juhani Koivisto. 2021. "Day-Ahead Market Modelling of Large-Scale Highly-Renewable Multi-Energy Systems: Analysis of the North Sea Region towards 2050" Energies 14, no. 1: 88. https://doi.org/10.3390/en14010088

APA Style

Gea-Bermúdez, J., Das, K., Koduvere, H., & Koivisto, M. J. (2021). Day-Ahead Market Modelling of Large-Scale Highly-Renewable Multi-Energy Systems: Analysis of the North Sea Region towards 2050. Energies, 14(1), 88. https://doi.org/10.3390/en14010088

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