1. Introduction
One of the essential elements in 4th generation district heating is a transition to smart energy systems with a better integration across sectors allowing to react to fluctuating renewable energy sources [
1]. Besides work looking into providing flexibility for renewable electricity integration within the power sector, work on a cross-sectoral integration including all sectors and/or energy carriers has been published [
2]. Hence, tools and approaches integrating different energy systems are needed [
3]. For providing flexibility to the power system, energy storage systems are crucial. As thermal energy storage tanks are several orders of magnitude cheaper than pumped hydro storage, an integration of the heating sector with the power market offers economically interesting flexibility options [
4]. Indeed, the thermal inertia of a heating grid can be used as heat storage as well, if it is considered in heat generation planning. Not requiring any investment into dedicated storage tanks, this approach allows to leverage additional flexibility at minimum cost.
Optimal scheduling of heat and electricity production units is traditionally done using unit commitment models with one energy balance for heat and one energy balance for power, neglecting grid dynamics [
5]. Integration of heating grid dynamics into scheduling of combined heat and power plants (CHPs) has been studied in several past publications. Most start from detailed physical models of the heating grid trying to achieve a model suitable for optimization by simplifying the grid topology [
6] or assuming simplifications like fixed mass flows or fixed time delays in an iterative optimization scheme [
7] or introducing a multi-step approach combining different simulation and optimization steps like static hydraulic, static thermal and dynamic thermal model [
8] or a unit commitment and a dispatch model [
9]. Li et al considers the building thermal inertia in addition to pipeline dynamics with fixed mass flows [
10]. Dai et al solve the same problem as Li considering pipeline dynamics with fixed mass flows and building thermal inertia in an iterative optimization scheme [
11]. As well considering thermal inertia of buildings and heating grid dynamics [
12] proposes an optimization model with fixed transport delays calculated based on design flows which allows to use off-the-shelf solvers. A very simple model of the heating grid as one dynamic energy mass is proposed and compared with other methods in a simulation in [
13]. Groß starts from a detailed physical model as well, but in contrast to others trains a linear regression model using multiple simulations. This regression function represents a direct link between energy stored within the heating grid and supply temperatures at the heat producer allowing a fast and efficient optimization [
14]. Giraud uses a similar function linking supply temperature changes at the heat producer with changes of the heat energy stored in the grid for optimization of supply temperature reducing heat and pumping losses. Its parameters are iteratively updated by a detailed grid simulation [
15].
All approaches offer feasible optimization solutions, however as first a detailed grid model is needed, engineering efforts to setup such a solution are high. In this paper we are introducing a new method to consider the thermal inertia of a heating grid as storage in CHP scheduling. Based on a classic unit commitment model for CHPs not considering grid dynamics presented at the beginning of
Section 3, we will first introduce an outer approximation of the thermal inertia of a heating grid in
Section 3.1. This formulation is inspired by models for dedicated energy storage (like batteries or heat storage tanks). In
Section 3.2 the heat energy stored in the grid will be modeled based on past and present supply temperatures. In contrast to [
14,
15] this function is created using time delays to the consumers as well as their share of consumption, allowing a more fast modeling as the creation of a detailed physical grid model is avoided. The fundamental dynamics of heating grids are explained in
Section 2.
Section 4 presents a motivating case study using the formulations of
Section 3.2 and
Section 5 discusses the results and shows an outlook on future developments and applications.
2. District Heating Grid Dynamics
In district heating grids heat energy is transported with water at different temperatures via a pipe network. Hence, hydraulic and temperature dynamics might be of interest when the thermal inertia of a heating grid should be used as energy storage. As hydraulic dynamics are very fast compared to the temperature dynamics in a heating grid, they can be neglected when looking at the energy dynamics of a heating grid [
8]. The temperature dynamics are however rather slow and have an important effect on the energy balances of a heating grid.
The water transporting the thermal energy is heated at production sites to a controllable supply temperature and transported to consumers via supply pipes. The consumers are controlling the mass flow of heated water depending on their heat demand and the supply temperature. This consumed water is cooled down and returned to the production sites via separate return pipes. The change of inner energy of water at producer and consumer in a time slot
can be described with:
where
is the mass flow of water between supply and return in time step
,
is the supply temperature,
is the return temperature and
is the specific heat capacity of water (assumed to be constant).
As the hydraulic dynamics can be neglected, it is assumed that mass flow changes induced by the consumers instantly influence the mass flow and hence the heat production at the producers. Changes in supply temperature by the producers however, are not directly seen at consumer stations, as the changed temperature first needs to propagate to the consumers with the water in the supply pipes. As the maximum velocity of water in a supply pipe is limited, there is a delay between an increase of supply temperature at the producer and the related increase of supply temperature at the consumer. This delay results in an inherent thermal energy storage, charging when the supply temperature is increased and discharging when the supply temperature is decreased.
This storage using the thermal inertia of a heating grid is explained with the simple example with one producer and one consumer having constant heat demand shown in
Figure 1a. In t2 the supply temperature at the generation is increased and following Equation (1) the heat output of the generation is increasing proportionally. As production in the following time slots is bigger than consumption, the heating grid shows a charging behavior. The increased supply temperature takes until t4 to reach the consumer. Now the consumer having a constant demand reduces the mass flow to keep Equation (1) in balance. Due to the fast hydraulics this mass flow change directly propagates to the producer and hence production is reduced to the original level and charging of the grid ends. When the supply temperature at the consumer is reduced in t7 we can see a similar inversed behavior discharging the grid until t10. The surface between produced energy and consumed energy is the energy charged or discharged in the heating grid.
For two consecutive consumers we get a behavior as shown in
Figure 1b. A change of supply temperature at the consumer in time t2 is having a similar effect of increased generation as before. When the supply temperature reaches the first consumer in t3 the production is decreased only by the share of the load/mass flow of this costumer (in our case 50%). When the second consumer is reached in t6 the remaining overproduction is reduced. A decrease in supply temperature at the producer in t8 shows as before an inverted storage profile discharging the heating grid.
3. Optimal Unit Commitment Models
Unit commitment problems are a very well-studied field for electric power generation [
16]. Mixed-integer linear programs have shown to be an efficient way to formulate and solve this problem with off-the-shelf open source or commercial solvers. Most unit commitment models are time discrete models. Hence, the time horizon is discretized into several time slots of equal size and for each time slot
the variables of the model are optimized. Key goal is to find optimal setpoints for the electric power output
and running status
(binary variable: 0:off/1:on) of each generator
in each time slot
. Electric power output and running status are linked with:
For each generator
in each time slot
with
and
being the minimum and maximum possible electric output of the generator. It is assumed that energy demand
for every time slot is known. If it is possible to sell and buy electricity from the electricity market, we introduce positive variables
and
being the electricity bought and sold to the energy markets. If no electric grid constraints are considered, we get the resulting electric energy balance:
Classic goals of unit commitment are either minimization of cost or maximization of profit. As objective function minimizing the overall cost we get:
With and being the time varying prices for selling or buying electricity at the energy markets and being the cost for running generation in a time slot and being the variable cost of generation of generator .
If heat generation of combined heat and power plants (CHPs) is considered, we get an additional variable
representing the heat output of the CHP. Like Equation (2) we get:
With
being the maximum heat output of CHP
. Electric output and heat output of a generator are linked. For gas turbines the power-to-heat ratio is usually constant. For CHPs with extraction condensing turbines both are linked as shown in
Figure 2 [
17].
If heating grid dynamics are not considered, we can formulate a heat energy balance like the electric energy balance Equation (3):
with
being the known overall heat demand and
being the thermal losses per time slot
, which are neglected in the following optimization. With
being the variable cost for heat generation the objective function becomes:
3.1. Outer Approximation of Heating Grid Dynamics Using State-of-Charge Formulation
The integration of dedicated thermal storage with storage tanks into a unit commitment problem can be achieved by introducing a state-of-charge variable
at each time slot
as well as a variable for the charging and discharging heat power
per time slot
. The state-of-charge can implicitly be formulated as:
where
is limited by the maximum charging and discharging power and
is limited by the maximum and minimum energy storage capacity of the thermal storage tank.
Introducing the state-of-charge formulation the heat energy balance Equation (6) without consideration of heat losses becomes:
With adapted limits of
and
this formulation can approximate the heat storage capability of a heating grid as shown in
Figure 1. For a given mass flow
the limits of the grid charging and discharging power can be calculated based on Equation (1):
with
being the maximum possible supply temperature and
being the minimum possible supply temperature in time slot
. The minimum possible supply temperature can be calculated upfront based on expected heat demand or outdoor temperature.
If the maximum time delay to the last consumer
is known, an upper limit of energy stored in the grid at time
can be formulated as a summation of the last upper bounds of charging power:
However, this formulation is only an upper bound of the maximum of the heat energy stored within the grid as only the last and not all consumers are feed with delay
(
Figure 1b). Hence, the presented Equations (8)–(11) leads to an outer approximation of the storage capabilities of a heating grid. It might be used for an approximation of economic potential but is not suitable for operations planning as results might be outside the feasible operation region.
3.2. Modeling of Heating Grid Dynamics with Estimation of Stored Energy
Looking at
Figure 1 the heat stored within the heating grid seems to be linked to the supply temperature. Hence, a function directly linking
to
is possible for heating grids with one major generation site if mass flows vary little. Groß proposes to train such a model based on simulations using linear regression [
14]. Here we are deriving such a function based on transport delays to consumers and consumers’ share in overall consumption. It is assumed that demand forecasts for all consumers are available. Additionally, it is assumed that initial supply temperatures as well as respective mass flows/transport delays are given. They could be calculated by a supply temperature optimization like [
15].
As in the previous section we can formulate the heat energy balance not considering heat losses:
However, to get a more exact representation of the thermal inertia of a heating grid, we need to calculate
differently. Assuming constant mass flow and hence constant transport delays, the supply temperature at the producer at time t2–t3 in
Figure 1a is increased by:
Leading to an instantaneous increase in CHP heat output by:
As soon as the supply temperature increase reaches the consumer after delay
in time t4–t5 the CHP output is decreased by:
Not considering heat losses the energy stored in the heating grid follows:
As the heat load is constant in
Figure 1a, the changes in heat output of the CHP directly are equal to changes in the energy charged to or discharged from the grid. Hence, combining Equation (14) and Equation (15) we get:
with
being the supply temperature at the producer,
being the time delay between producer and consumer and
being a constant mass flow.
If a time varying mass flow
is introduced as a parameter, Equation (17) becomes:
To consider multiple consumers, as in the example in
Figure 1b, we introduce a matrix M based on time delay to and share of consumption of a consumer. In this matrix row 1 has non-zero entries whenever a water volume starting at time one at the producer reaches a consumer; row 2 is having non-zero entries when a water volume starting at time two at the producer reaches a consumer et cetera. The non-zero values in M matrix correspond the share of consumption of the consumer being reached. Those shares are assumed to be constant. Hence every row sums up to a value of 1 if all consumers are reached within the time horizon.
Looking at the example given in
Figure 1b an increase in supply temperature at the producer in t2 reaches the first consumer in t4 and the second consumer in t6. Having equal share in overall consumption non-zero entries in row t2 are 0.5 in column t4 and 0.5 in column t6 (
Table 1). Having a row for each time instant allows to account for variable time delays if mass flows are not constant.
Extending the calculation of
Equation (18) with matrix M (
Table 1) to allow multiple consumers we get:
With Equation (12), Equation (19) and
Table 1 we know have a model which allows to directly calculate the currently charged or discharged thermal energy in the heating grid based on a supply temperature profile, assuming demands of and transport delays to consumers as well as mass flows at the producer are known.
4. Motivating Case Study
As a motivating case study, we have a closer look at a heating grid with one CHP and two consumers (comparable to
Figure 1b). The CHP is having a maximum electric output of 500 MW and a maximum heat output of 400 MW. Electricity is produced using an extraction condensing turbine and hence power-to-heat ratio is variable within the limits shown in
Figure 2. Heat consumer 1 having a load share of about 45% and a peak load of 120 MW is connected to the CHP with a 20 km pipe with 0.7 m diameter. Heat consumer 2 having a load share of about 55% and a peak load of 145 MW is following consumer 1 with a 10 km pipe with 0.7 m diameter. Electricity prices are taken from EPEX Day-Ahead Spot Market for Germany. The optimization horizon is one day with 24 time slots (one per hour). Matrix M as well as mass flows are calculated assuming that water in the pipes is flowing at maximum velocity of 1.5 m/s.
If supply temperature can be increased for heat storage by a maximum of 30 K, the production cost can be reduced by 2.4% (several k€/day) compared to a scenario without heat storage. Corresponding results with supply temperature and heat power profiles are shown in
Figure 3.
Table 2 shows achievable savings for different maximum increases of supply temperature. As expected the savings increase, with increasing flexibility of supply temperature and hence increasing storage capacity of the heating grid. As shown in
Table 3 variations of the demand share have a rather small influence on the results. Variations of velocity lead to a bit larger changes, however remain about one order of magnitude smaller than the achieved savings. In almost all cases the supply temperature profile does not change. Only for a reduction of velocity by 10% the supply temperature profile shows minor changes like a lower temperature increase of 22.6 K instead of 30 K in hour 6. The optimization problem is modeled in Julia/JuMP [
18] and solved with Gurobi 8.0 [
19] on an Intel® Xeon® v2 with 2.6 GHz (2 processors) and 4 GB RAM. Solution times for all problems shown in
Table 2 and
Table 3 are around 16 ms never exceeding 20 ms, hence an application for online-optimization seems in reach even for larger heating grids and real-world cases.
5. Discussion/Comparison of Modeling Approaches
In this paper we introduced two different ways to model the thermal inertia of heating grids without consideration of thermal losses. First, an outer approximation was derived using a simple state-of-charge formulation. As this representation did not lead to accurate results, second, an estimation of the stored energy using past and present supply temperatures was introduced. This formulation is reducing the modeling effort compared to many other approaches, as only maximum supply temperature deviations, transport delays to consumers and their respective load shares need to be parametrized. Hence, an extension to real world district heating grids is possible without major efforts, if consumers are grouped into consumption areas having similar transport delay from the producer. This consumption areas and respective transport delays can be used to create the M matrix.
As some but not all approaches in literature, this formulation only allows us to model dynamics of a heating grid feed by one generation site. Another limitation is that constant transport times from heat source to consumers are assumed. This leads to an error as when an increased temperature reaches a consumer the corresponding mass flow and hence velocity and transport time are changed to keep Equation (1) in balance. If temperature increases are rather short and small, this is neglectable, but with long and large temperature increases this can decrease velocity importantly. One way to overcome this limitation could be to use iterative updates with a simulation model like in [
7].
The results of the motivating case study show that the presented formulation is computationally very efficient and outperforms many existing algorithms. The model is quite robust against small parameter variations and it allows interesting cost savings of several thousand Euros per day as shown with calculations for different maximum supply temperature increases. As a high temperature variation leads to increased component fatigue (e.g., for pipes), in real world applications the maximum increase of supply temperature will most likely be limited to 10K or 20 K being intra-day temperature variations which can be seen in today’s operations as well.
In the case study only one pipe with a large diameter is considered and hence heat losses can be neglected. However, if a distribution grid is included, increase of thermal losses is becoming more important. Hence, future work will try to overcome the inaccuracies of the fixed mass flows, address the integration of heat losses and apply the second approach to real-world heating grids.
Author Contributions
Conceptualization, L.M. and A.A.H.; methodology, L.M. and A.A.H.; software, A.A.H.; validation, L.M. and A.A.H.; formal analysis, L.M.; investigation, A.A.H.; resources, L.M.; data curation, L.M.; writing—original draft preparation, L.M.; writing—review and editing, L.M.; visualization, L.M.; supervision, S.H.; project administration, L.M.; funding acquisition, L.M. and S.H.
Funding
The authors gratefully acknowledge funding by the German Federal Ministry of Education and Research (BMBF) within the Kopernikus Project ENSURE ‘New ENergy grid StructURes for the German Energiewende’.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Lund, H.; Werner, S.; Wiltshire, R.; Svendsen, S.; Thorsen, J.E.; Hvelplund, F.; Mathiesen, B.V. 4th Generation District Heating (4GDH): Integrating smart thermal grids into future sustainable energy systems. Energy 2014, 68, 1–11. [Google Scholar] [CrossRef]
- Lund, H.; Østergaard, P.A.; Connolly, D.; Mathiesen, B.V. Smart energy and smart energy systems. Energy 2017, 137, 556–565. [Google Scholar] [CrossRef]
- Lund, H.; Østergaard, P.A.; Chang, M.; Werner, S.; Svendsen, S.; Sorknæs, P.; Thorsen, J.E.; Hvelplund, F.; Mortensen, B.O.G.; Mathiesen, B.V.; et al. The status of 4th generation district heating: Research and results. Energy 2018, 164, 147–159. [Google Scholar] [CrossRef]
- Lund, H.; Østergaard, P.A.; Connolly, D.; Ridjan, I.; Mathiesen, B.V.; Hvelplund, F.; Thellufsen, J.Z.; Sorknæs, P. Energy storage and smart energy systems. Int. J. Sustain. Energy Plan. Manag. 2016, 11, 3–14. [Google Scholar] [CrossRef]
- Atanasoae, P. The operating strategies of small-scale combined heat and power plants in liberalized power markets. Energies 2018, 11, 3110. [Google Scholar] [CrossRef]
- Wigbels, M.; Bohm, B.; Sipilae, K. Dynamic Heat Storage Optimisation and Demand Side Management; 2005:8DHC-05.06; IEA Annex VII.: Sittard, Netherlands, 2005. [Google Scholar]
- Li, Z.; Wu, W.; Shahidehpour, M.; Wang, J.; Zhang, B. Combined heat and power dispatch considering pipeline energy storage of district heating network. IEEE Trans. Sustain. Energy 2016, 7, 12–22. [Google Scholar] [CrossRef]
- Tröster, S. Zur Betriebsoptimierung in Kraft-Wärme-Kopplungssystemen unter Berücksichtigung der Speicherfähigkeit des Fernwärmenetzes. Ph.D. Thesis, Universität Dortmund, Dortmund, Germany, 1999. [Google Scholar]
- Schweiger, G.; Larsson, P.-O.; Magnusson, F.; Lauenburg, P.; Velut, S. District heating and cooling systems —Framework for Modelica-based simulation and dynamic optimization. Energy 2017. [Google Scholar] [CrossRef]
- Li, P.; Wang, H.; Lv, Q.; Li, W. Combined heat and power dispatch considering heat storage of both buildings and pipelines in district heating system for wind power integration. Energies 2018, 10, 893. [Google Scholar] [CrossRef]
- Dai, Y.; Chen, L.; Min, Y.; Chen, Q.; Hao, J.; Hu, K.; Xu, F. Dispatch model for CHP with pipeline and building thermal energy storage considering heat transfer process. IEEE Trans. Sustain. Energy 2019, 10, 192–203. [Google Scholar] [CrossRef]
- Gu, W.; Wang, J.; Lu, S.; Luo, Z.; Wu, C. Optimal operation for integrated energy system considering thermal inertia of district heating network and buildings. Appl. Energy 2017, 199, 234–246. [Google Scholar] [CrossRef]
- Lesko, M.; Bujalski, W. Modeling of district heating networks for the purpose of operational optimization with thermal energy storage. Arch. Thermodyn. 2017, 38, 139–163. [Google Scholar] [CrossRef]
- Groß, S. Untersuchung der Speicherfähigkeit von Fernwärmenetzen und deren Auswirkungen auf die Einsatzplanung von Wärmeerzeugern. Ph.D. Thesis, TU Dresden, Dresden, Germany, 2012. [Google Scholar]
- Giraud, L. Modélisation dynamique et gestion avancée de réseaux de chaleur. Ph.D. Thesis, Université Grenoble Alpes, Grenoble, France, 2016. [Google Scholar]
- Tahanan, M.; van Ackooij, W.; Frangioni, A.; Lacalandra, F. Large-scale Unit Commitment under uncertainty. 4OR-Q. J. Oper. Res. 2015, 13, 115–171. [Google Scholar] [CrossRef]
- Cerbe, A. Prozessnahe Einsatzoptimierung mit BoFiT unter Berücksichtigung der Netzrestriktionen. Verbundprojekt Bessere Ausnutzung von Fernwärmeanlagen; Pilotprojekt HEW Hamburg: Hamburg, Germany, 2002. [Google Scholar]
- Dunning, I.; Huchette, J.; Lubin, M. JuMP: A modeling language for mathematical optimization. SIAM Rev. 2017, 59, 295–320. [Google Scholar] [CrossRef]
- Gurobi Optimization, I. Gurobi Optimizer Reference Manual. Available online: http://www.gurobi.com (accessed on 12 December 2018).
© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).