Multi-Time Period Optimal Dispatch Strategy for Integrated Energy System Considering Renewable Energy Generation Accommodation

: With the target of carbon peaking and carbon neutrality, renewable energy generation (REG) develops rapidly. The increasing penetration of REG brings along the problems of ﬂuctuation in power ﬂow and the possible abandonment of wind and photovoltaics (PV) generation. In this context, the so-called integrated energy system (IES) becomes a promising solution to the accommodation of REG thanks to energy storage systems and coupling devices inside. In this paper, the optimal operation model of an IES is ﬁrst presented, with the schemes of green certiﬁcate trading and carbon emission right trading included to provide economic incentives for accommodating REG. Next, in order to address the problem of uncertainty in REG, the devices in the IES are divided into three types based on regulation ﬂexibility, and a multi-time period optimal dispatching scheme is proposed, including day-ahead optimal scheduling, rolling optimal dispatching, and real-time control strategy. Finally, it is demonstrated by simulation results of a numerical example that the proposed method not only promotes the accommodation capability for REG but can also cope well with contingencies.


Introduction
Under the target of carbon peaking and carbon neutrality, the development of renewable energy generation (REG) advances rapidly.It is anticipated that renewable capacity will expand by over 1800 GW, or over 60%, by 2026, among which China ranks the first, accounting for 43% of global growth [1].However, the output of REG units, such as photovoltaics and wind farms, depends heavily on the meteorological condition, which is hard to predict accurately.The intermittency in REG and the difference between the generation prediction and the actual output may disrupt the supply-demand balance of a power system, bringing along fluctuations in power flow.Then, if the operation plan of the power system is not adjusted accordingly and timely, the REG may not be accommodated by the power system, and hence some wind and solar generation has to be curtailed.On the other hand, if the power output from REG is insufficient, power shortage may occur.Both cases have to be properly addressed so as to promote the development of REG [2,3].In this context, the integrated energy system (IES), with the energy storage systems and coupling devices inside, becomes a promising solution to the increasing penetration of REG [4,5].When there is excessive REG, the IES can store the electricity or convert it into other kinds of energy such as heat and cold; when there is insufficient REG output, the IES can still meet the users' load demand with other energy sources, such as natural gas.
Energies 2022, 15, 4329 2 of 18 Many publications are available regarding the modeling, planning, operation, and assessment of an IES.In [6][7][8], the notions of the so-called Energy Hub and Energy Internet are introduced, and the general modeling and simulation approaches of the IES are reviewed.Some research is devoted to extending the scope of the IES, incorporating Power-to-Gas (P2G) devices and dynamic pipeline networks [9], energy storage systems in buildings [10], demand response of users [11], or distributed generators [12] into the optimal dispatching model, and proves that the extended IES is more robust and economic under multi-energy flow scenarios.Some research focuses on promoting the accommodation of REG.In [13], a weighting factor is introduced to balance the minimization of grid operation cost and REG curtailment.In [14], the hybrid multi-objective particle swarm optimization (HMOPSO) method is adopted to minimize economic cost, maximize REG utilization, and minimize expected energy not supplied (EENS) at the same time.In [15], the cost of REG curtailment caused by flexibility shortage of conventional generators is included in the objective function.However, the above optimal dispatching models cannot provide direct economic incentives for IES operators to maximally accommodate REG and reduce carbon emission.The market mechanisms of tradable green certificates (TGCs) [16][17][18] and carbon emission right trading [19][20][21] are designed to address this problem, but little research on the optimal operation of the IES has been conducted with both trading schemes considered.
To cope with uncertainties in IES operation, there are mainly two approaches.One is to mathematically characterize the uncertainty in the optimization model with chanceconstrained programming [22], stochastic differential equations [23], stochastic programming [24], and so on.Another branch is the multi-time scale control scheme, because the smaller the time scale is, the more accurate the prediction could be.A two-stage scheduling strategy is proposed in [25], where the rolling optimization aims at minimizing the operation cost and startup/shutdown penalty under the time of use tariff mechanism to formulate the large time-scale schedule, and the dynamic adjustment stage deals with the small time-scale uncertainties of REG and loads.In [26], a day-ahead scheduling model and an intraday scheduling model are first presented, and the intraday model is then divided into three sub-layers of slow control, intermediate control, and fast control to perform rolling optimization.
In general, the existing research on the optimal dispatch of an IES mainly focuses on the economy and reliability, and there are relatively few studies on promoting REG accommodation with market mechanisms.Additionally, the current studies on multi-time scale control schemes mainly present a framework with control strategies of different time scales for the whole IES, but few of them study the flexibility characteristics of devices and formulate different control strategies for different devices.Therefore, the major contributions of our work are as follows: (1) The TGC trading scheme and carbon trading scheme are incorporated into the modeling of the IES operation so as to provide economic incentives for promoting REG accommodation as well as to quantify the REG accommodation.(2) Based on the ramping characteristics, the devices in the IES are classified into three types to participate in different stages of the operation optimization with different margins, allowing full play to the IES's flexibility.(3) A multi-time period optimal dispatch scheme of an IES is presented to address the correction in REG output prediction, dividing the operation planning of the IES into day-ahead optimal scheduling, rolling optimal dispatching, and real-time control.(4) The effects of each stage in the proposed dispatch strategy are compared and analyzed under normal operation and emergency, respectively, providing useful guidance for further study in the IES control strategy formulation.
The rest of this paper is organized as follows.In Section 2, the optimal operation model of an IES with schemes of TGCs and carbon trading is established.In Section 3, devices in the IES are categorized into three types based on flexibility, and a multi-time period optimal dispatching scheme is presented.In Section 4, a numerical example is employed to demonstrate the presented method and the results are analyzed.Finally, in Section 5, concluding remarks are given.

Components
An IES, also known as an Energy Internet, is a comprehensive energy production, supply, and marketing unit.Various devices for the production, conversion, storage, and transmission of electricity, gas, heat, and cold energy coordinates with each other to provide integrated energy services.Different structures of IESs are designed for different regions with different load conditions, and the structure of the IES studied in this paper is shown in Figure 1, in which the photovoltaics are chosen as representative REG units.
The rest of this paper is organized as follows.In Section 2, the optimal operation model of an IES with schemes of TGCs and carbon trading is established.In Section 3 devices in the IES are categorized into three types based on flexibility, and a multi-time period optimal dispatching scheme is presented.In Section 4, a numerical example is employed to demonstrate the presented method and the results are analyzed.Finally, in Section 5, concluding remarks are given.

Components
An IES, also known as an Energy Internet, is a comprehensive energy production supply, and marketing unit.Various devices for the production, conversion, storage, and transmission of electricity, gas, heat, and cold energy coordinates with each other to provide integrated energy services.Different structures of IESs are designed for different regions with different load conditions, and the structure of the IES studied in this paper is shown in Figure 1, in which the photovoltaics are chosen as representative REG units.As shown in Figure 1, the IES interacts with the power grid and natural gas network and supplies electric load, gas load, heating load, and cooling load to users.There are mainly two kinds of devices: standalone devices for the production and storage of a single type of energy and coupling devices that involve the conversion of different types of energy.In the IES to be studied below, standalone devices include photovoltaic generators (pv), electric storage (es), heat storage (hs), and cold storage (cs); coupling devices include gas turbine (gas-electricity, ge), gas boiler (gas-heat, gh), electric boiler (electricity-heat eh), electric chiller (electricity-cold, ec), absorption chiller (heat-cold, hc), and waste heat boiler (heat-heat, whb).It is worth noting that the waste heat boiler is used to collect the heat in flue gases from the combustion of the gas turbine to generate hot water or steam allowing the subsequent use of heat energy.The waste heat boiler, together with the gas turbine and absorption chiller, comprises the Combined Cooling, Heating, and Power (CCHP) unit, which can supply electricity, heat, and cold energy at the same time.As shown in Figure 1, the IES interacts with the power grid and natural gas network, and supplies electric load, gas load, heating load, and cooling load to users.There are mainly two kinds of devices: standalone devices for the production and storage of a single type of energy and coupling devices that involve the conversion of different types of energy.In the IES to be studied below, standalone devices include photovoltaic generators (pv), electric storage (es), heat storage (hs), and cold storage (cs); coupling devices include gas turbine (gas-electricity, ge), gas boiler (gas-heat, gh), electric boiler (electricity-heat, eh), electric chiller (electricity-cold, ec), absorption chiller (heat-cold, hc), and waste heat boiler (heat-heat, whb).It is worth noting that the waste heat boiler is used to collect the heat in flue gases from the combustion of the gas turbine to generate hot water or steam, allowing the subsequent use of heat energy.The waste heat boiler, together with the gas turbine and absorption chiller, comprises the Combined Cooling, Heating, and Power (CCHP) unit, which can supply electricity, heat, and cold energy at the same time.

Capacity
First of all, each device in the IES is bounded by the capacity constraint in (1), where P t i represents the output of each device as well as the power transmitted through the power grid and natural gas network at time t, S t i represents the state of charge of each storage device at time t, P i , S i , P i , and S i are the corresponding lower bounds and upper bounds.Note that for conventional devices, the upper bound is determined by the installed capacity, while for the photovoltaic generators, the output is also bounded by the solar condition.P i ≤ P i t ≤ P i , ∀i ∈ {pv, es, hs, cs, ge, gh, eh, ec, hc, whb, grid, gas}

Conversion
The coupling device takes one type of energy as input and produces another type of energy, following the power conversion constraints defined in (2), where η i is the conversion efficiency of each coupling device, P t i,out and P t i,in are the corresponding output and input at time t.
There is also auxiliary power consumption in the IES.For example, some devices may use a water pump during operation.The auxiliary power P t w is calculated in (3), where P t i is the output and k i is the corresponding coefficient.

Storage
The energy storing and releasing process of the storage system is described by (4), where S t i , P t i,s , and P t i,r are respectively the state of charge, charging power, and discharging power at time t; σ i , η i,s , and η i,r are respectively the loss coefficient, charging efficiency, and discharging efficiency of the device.
To ensure the operation of the IES on the next day and assuming that each day is identical, the storage system needs to restore the initial capacity at the end of each day, as shown in (5).
Considering that an energy storage system cannot be charging and discharging at the same time, constraint in (6) needs to be enforced.
However, ( 6) is a quadratic constraint which introduces nonlinearity into the optimization model.By introducing two binary indicators of the storage status, B t i,s and B t i,r , (6) can be converted into a set of linear constraints as in (7).When the storage device is storing energy (charging), B t i,s = 1 and B t i,r = 0; when the storage device is releasing energy (discharging), B t i,s = 0 and B t i,r = 1.

Conservation
The energy conservation constraints of an IES include electricity, gas, heat, and cold, as shown in ( 8 where P t hc,h is the heat input of the absorption chiller at time t; P t whb,i is the heat input of the waste heat boiler at time t; k ge,h is the waste heat coefficient of the gas turbine.

Objective Function
To promote the economy and environmental friendliness of the IES, in addition to the external energy cost C e and the operation cost C op , the cost of tradable green certificates C tgc and carbon emission trading C em are also included in the objective function of optimal dispatching:

External Energy
As shown in Figure 1, the expenses of external energy include the costs of purchasing natural gas from the gas network and purchasing electricity from the grid as in (11), where ξ t gas and ξ t grid are respectively the price of buying natural gas and electricity at time t.
C e = T ∑ t=0 P t gas ξ t gas + P t grid ξ t grid ∆t (11)

Operation and Maintenance
The operation and maintenance cost C op is defined in (12), where ξ i,op is the unit output energy cost of device i.

Tradable Green Certificate
The TGC is a market-based cost-efficient tool to regulate the deployment of renewable energy.The TGC cost model can be expressed as a piecewise function with three stages as shown in (13) [27].Firstly, if the actual accommodated REG ratio is greater than the specified quota ratio, the IES is profitable by selling the REG power, which is reflected in the model as selling TGCs.Secondly, if the actual accommodated REG ratio is smaller than the specified quota ratio, but still within the penalty margin, the gap needs to be filled by purchasing TGCs.Thirdly, if the actual REG is woefully inadequate and exceeds the penalty margin, the IES not only needs to purchase TGCs but will also be penalized heavily.
where k, u, and s are respectively the specified REG ratio, actual REG ratio, and the penalty margin; ξ tgc and ξ p are respectively the price of the TGC and the penalty; P t a , P t pv , and P t grid are respectively the total generation output, accommodated photovoltaic output, and electricity purchased from the grid at time t.The specified ratio k is usually no less than 15% [28].

Carbon Emission Right
Similar to the mechanism of the TGC but targeting the carbon emission, the IES can gain profits if the carbon emission is smaller than the allocated quota, and may need to purchase carbon emission rights if the actual emission is larger than the quota [29], as shown in (14), where ξ em is the price of carbon emission rights, E t D is the carbon emission quota of the system at time t, and M t D is the actual emission.
The allocated emission quota is calculated in (15), which is proportional to the electric load.The actual carbon emission of the system is calculated in (16), i.e., each carbon emission source contributes proportionally to its output.In our optimization model, the carbon emission mainly comes from the purchased electricity M t grid , the CCHP unit (the gas turbine (ge), waste heat boiler (whb), and the absorption chiller (hc)) M t cchp , and the gas boiler M t gb .For the purchased electricity from the grid, the carbon emission quota is proportional to its quantity; for the CCHP unit, the emission quota is calculated based on the equivalent electricity generation; for the gas boiler (gh), the quota is proportional to its heat supply.In view of the intermittency and uncertainty of REG such as photovoltaics, its output prediction is of high difficulty, making the single day-ahead scheduling method far from meeting the economic, technical, and environmental needs of the IES.However, the shorter the prediction period, the more accurate the prediction result.Therefore, in the following section, a multi-time period optimal dispatching scheme is proposed, dividing the operation of the IES into three stages, as shown in Figure 2. Firstly, the day-ahead scheduling plan of the IES is determined based on the day-ahead prediction of REG output and load demands, which is an hourly plan.Secondly, with more information during the operation day, a rolling prediction of REG output and load demands with higher precision can be obtained, and the dispatching of the IES is updated, which is at the hour or minute level.Finally, in case of an emergency, a real-time control strategy will respond in minutes or seconds to ensure the secure and economic operation of the IES.

Multi-Time Period Optimal Dispatching Scheme
In view of the intermittency and uncertainty of REG such as photovoltaics, its outpu prediction is of high difficulty, making the single day-ahead scheduling method far from meeting the economic, technical, and environmental needs of the IES.However, the shorter the prediction period, the more accurate the prediction result.Therefore, in the following section, a multi-time period optimal dispatching scheme is proposed, dividing the operation of the IES into three stages, as shown in Figure 2. Firstly, the day-ahead scheduling plan of the IES is determined based on the day-ahead prediction of REG output and load demands, which is an hourly plan.Secondly, with more information during the operation day, a rolling prediction of REG output and load demands with higher precision can be obtained, and the dispatching of the IES is updated, which is a the hour or minute level.Finally, in case of an emergency, a real-time control strategy wil respond in minutes or seconds to ensure the secure and economic operation of the IES.Based on the response time and ramping characteristics, the devices in the IES are classified into three types to participate in the dispatching schemes of different time scales, as shown in Table 1.Since the operation of the IES is based on the day-ahead scheduling, the type A device mainly includes the base-load units that are not suitable for frequent start-stop or power adjustment, and will operate with a low margin during day-ahead scheduling.The type B device includes some flexible units such as the gas turbine.It participates in day-ahead scheduling with a high margin and participates in rolling optimization with a low margin.Furthermore, some devices with strong flexibility such as the electricity storage units are selected as type C devices to deal with emergencies, participating in day-ahead scheduling and rolling dispatching with a high margin, and participating in real-time control with a low margin.The external energy network is adjustable in all scales of dispatching schemes to ensure energy balance.

Day-Ahead Optimal Scheduling
The day-ahead scheduling plan of the IES is based on the day-ahead prediction of energy supply and load demands.The steps are shown in Figure 3 and are listed as follows.
The day-ahead scheduling plan of the IES is based on the day-ahead pred energy supply and load demands.The steps are shown in Figure 3 and are follows.
1. Collect historical load data and predict the demand curve of electricity, gas, h cold of the next day.2. Collect meteorological information and predict the REG output.3. Add margin constraints for type A, B, and C device into the optimizatio described in Section 2 and feed in the source and load data above.4. Obtain a day-ahead scheduling plan of the IES.

1.
Collect historical load data and predict the demand curve of electricity, gas, heat, and cold of the next day.

2.
Collect meteorological information and predict the REG output.

3.
Add margin constraints for type A, B, and C device into the optimization model described in Section 2 and feed in the source and load data above.4.
Obtain a day-ahead scheduling plan of the IES.

Rolling Optimal Dispatching
Due to the uncertainty of REG output, the scheduling plan formulated above may not satisfy the power balance.When the actual generation output is less than prediction and there is a power shortage, it is necessary to purchase electricity from the external grid with a high price; when the actual generation output is greater than prediction and there is a power redundancy, there will be "abandonment of wind and light".For the prediction of REG output, the closer the time is, the higher the accuracy.Therefore, based on the rolling prediction of REG output, the dispatching plan of the IES can be revised on the time scale of hours or minutes.The steps are shown in Figure 4 and are listed as follows.

1.
At time t update the REG output prediction and load prediction of [t, T] based on the data of [0, t].

2.
Type A devices do not participate in the rolling dispatching, so fix its output to the day-ahead scheduling plan.

3.
Add margin constraints for type B and C devices into the optimization model and feed in the updated source and load data.

4.
Obtain an updated operation plan of the IES and dispatch accordingly in [t, t + 1].

5.
Increment t and repeat steps 1 to 4.

Real-Time Control Strategy
The IES should also be capable of dealing the emergencies such as a sudden energy supply, sudden increase in load demands, or failures of devices.The int of the IES with renewable generators, energy conversion devices, energy storage and energy networks can strengthen the coupling relationship between the power natural gas system, and thermal system through the internal complementar energy flow, providing robust energy services and improving energy ut efficiency.In the case when the change is too sudden to be predicted, if the sys operates under the day-ahead and rolling plan, the economic and environmental will be affected.Therefore, taking the sudden change in REG output as an exampl time control strategy to respond to emergencies is proposed, as shown in

Real-Time Control Strategy
The IES should also be capable of dealing the emergencies such as a sudden drop in energy supply, sudden increase in load demands, or failures of devices.The integration of the IES with renewable generators, energy conversion devices, energy storage systems, and energy networks can strengthen the coupling relationship between the power system, natural gas system, and thermal system through the internal complementary multi-energy flow, providing robust energy services and improving energy utilization efficiency.In the case when the change is too sudden to be predicted, if the system still operates under the day-ahead and rolling plan, the economic and environmental benefits will be affected.Therefore, taking the sudden change in REG output as an example, a real-time control strategy to respond to emergencies is proposed, as shown in Figure 5.

1.
At time t, determine whether to respond to the emergency by comparing the output change ∆Pδt to the threshold S th .If the output change is greater than the threshold, a real-time control scheme will be generated, correcting the operation curve of the IES on the time scale of minutes or seconds.

2.
Update the prediction of emergency during [t, t + t 0 ].

3.
Type A and B devices do not participate in the real-time control, so fix its output to the rolling dispatching plan.

4.
Add margin constraints for type C devices into the optimization model and feed in the updated source and load data.

5.
Obtain an optimal operation plan and control the IES accordingly in [t, t + 1]. 6.
Increment t and repeat steps 1 to 5.

7.
When the emergency is over, the IES returns to the original dispatching plan.
3. Type A and B devices do not participate in the real-time control, so fix its output to the rolling dispatching plan.4. Add margin constraints for type C devices into the optimization model and feed in the updated source and load data. 5. Obtain an optimal operation plan and control the IES accordingly in [ , 1] t t + .
6. Increment t and repeat steps 1 to 5.
7. When the emergency is over, the IES returns to the original dispatching plan.

Settings of the Example
The energy consumption data in a representative summer day of a residential area and the operation parameters of an IES are employed to illustrate the proposed multi-time period optimal dispatching scheme.The electricity load demand and the cooling load demand are shown in Figure 6, the time-of-day tariff is shown in Figure 7, the device parameters of the IES are shown in Table 2, the economic parameters are shown in Table 3, and the margin parameters of the proposed control strategy are shown in Table 4.

Settings of the Example
The energy consumption data in a representative summer day of a residential area and the operation parameters of an IES are employed to illustrate the proposed multi-time period optimal dispatching scheme.The electricity load demand and the cooling load demand are shown in Figure 6, the time-of-day tariff is shown in Figure 7, the device parameters of the IES are shown in Table 2, the economic parameters are shown in Table 3, and the margin parameters of the proposed control strategy are shown in Table 4.

Dispatching Results in Normal Operation
The day-ahead prediction, rolling prediction, and actual photovoltaic output of a representative summer day are shown in Figure 8.The day-ahead prediction is a smooth curve with a peak at noon, but due to the meteorological uncertainties during the day, the actual output varies a lot from the day-ahead prediction.With the intraday information, the rolling prediction is much closer to the actual output.
Economic results of the IES are listed in Table 5, from which we can see that, in terms of the total cost, the dispatching strategy of only day-ahead scheduling leads to the most expensive outcome, the dispatching strategy of day-ahead scheduling together with rolling optimization results in the outcome with much less costs, and the proposed three-stage control scheme produces the most economic outcome.The reduction in total cost mainly comes from the reduction in purchasing electricity from the power grid.Too much cost spent on purchasing electricity indicates that there is a large gap between the supply and demand of the IES, which is not only an economic burden but may also hinder the security and stability of power system operation.When the IES operates according to the day-ahead scheduling plan only, any variation of the PV generation output from the day-ahead prediction has to be compensated by the external power grid, hence there is a large expense in buying electricity from the power grid.When the IES operates with rolling optimization, the output prediction is updated throughout the day and the dispatching of the IES can be adjusted accordingly and continually, hence the purchased electricity from the external grid is much less, and the grid cost is dropped by 82%.When real-time control is further incorporated, the IES can respond much more timely and effectively to the output fluctuation with the help of the energy coupling and storage devices inside, further reducing the grid cost by 33%.The costs of TGC trading and carbon trading also vary significantly between "dayahead" and "day-ahead + rolling", which can reflect the relationship between the generation outputs of REG units and conventional generators.The more the REG is dispatched, the smaller the TGC cost (larger income) would be; the more the generation outputs from conventional generators such as gas turbines are, the larger the carbon trading cost (less income) would be.For further examination, the detailed dispatching results of the IES for "day-ahead" and for "day-ahead + rolling" are demonstrated in Figure 9, and the corresponding PV generation dispatching results are shown in Figure 10.When the IES operates according to the day-ahead scheduling plan only, all types (A, B, and C) of devices follow the day-ahead hourly schedule, so the PV generation output cannot be fully dispatched even if there is more power output available, resulting in abandonment of PV generation.When the IES operates with rolling optimization, the operation of type B and C devices can be adjusted with rolling updates of PV generation output prediction.As a result, the generation outputs of gas turbine and electric storage systems are reduced when there is more PV generation output available, and are increased when the PV generation output drops, hence reducing the abandonment of PV generation as well as the purchasing cost from the power grid.

Dispatching Results in Emergency
The photovoltaic outputs in normal operation and emergency are shown in Figure 11, from which we can see that due to some situations, the photovoltaic output suddenly drops at around 9:00.Economic results of the IES are shown in Table 6, and the dispatching results for "dayahead + rolling" and for "day-ahead + rolling + real-time" are demonstrated in Figure 12.It can be seen from Table 6 that the change in grid cost dominates.Comparing the dayahead scheduling and rolling optimization, though in both schemes a certain amount of electricity has to be bought from the power grid in response to the emergency, rolling optimization needs less external electricity during normal hours, so its grid cost decreases by 54% compared to day-ahead scheduling.In the real-time control scheme, thanks to the flexibility of the type C devices, the power gap is quickly filled by discharging electric storage units, hence the grid cost is further reduced by 64%.To see the advantage of the proposed multi-time period dispatching strategy, it can be seen from Figure 12 that the sudden drop of PV generation output results in a huge power gap.In the "day-ahead + rolling" scheme, the gap is filled by purchasing a large amount of external electricity from the power grid, which not only results in a large cost but may also lead to the violation of the transmission capacity.In the "day-ahead + rolling + real-time" scheme, the electric storage unit releases power to respond to the accident, and the IES buys electricity from the power grid later at noon when the electricity price is low to charge the storage unit, hence reducing grid cost and minimizing the risk of cascading faults.

Conclusions
As a comprehensive energy service provider with various energy coupling and storage devices, the IES is a feasible solution to accommodate the fluctuant REG output.In this paper, the operation model of the IES is first presented, and the schemes of TGC trading and carbon trading are incorporated.Secondly, a multi-time period optimal dispatching scheme is proposed, and the devices in the IES are classified into three types according to their flexibility to participate in different scales of the scheme.Finally, the performance of the presented dispatching scheme is demonstrated with a numerical example, and simulations of the IES operation are carried out under normal operation and emergency respectively.The results show that, firstly, the day-ahead scheduling drafts a plan for the IES operation, then the rolling optimization refines the plan with the updated prediction during the day, and finally the real-time control responds to the unpredicted conditions with the flexible electric storage devices.The results demonstrate that the proposed IES model, together with the multi-time period optimal dispatching scheme, can not only promote the accommodation of renewable energy but can also cope with contingencies smoothly.

Figure 1 .
Figure 1.Structure of a sample IES.

2 .
Type A devices do not participate in the rolling dispatching, so fix its outp day-ahead scheduling plan.3. Add margin constraints for type B and C devices into the optimization mo feed in the updated source and load data.4. Obtain an updated operation plan of the IES and dispatch accordingly in [ , t5.Increment t and repeat steps 1 to 4.

Figure 1 .
At time t , determine whether to respond to the emergency by comparing th change t Pδ Δ to the threshold th S .If the output change is greater t

Figure 5 .
Figure 5. Illustration of real-time control strategy.

Figure 6 .
Figure 6.Load demand profiles of the example.

Figure 5 .
Figure 5. Illustration of real-time control strategy.

Figure 6 .
Figure 6.Load demand profiles of the example.

Figure 7 .
Figure 7. Electricity price of the example.

Figure
Figure Detailed dispatching results of the IES in normal operation.

Figure 10 .
Figure 10.PV dispatching results in normal operation.
): Cload are respectively the electricity, gas, heating, and cooling load at time t; P t pv is the output of photovoltaics at time t; P t es,r , P t hs,r , and P t cs,r are respectively the discharging power of electric storage, heat storage, and cold storage at time t; P t es,s , P t hs,s , and P t cs,s are respectively the charging power of electric storage, heat storage, and cold storage at time t; P t ge,g and P t ge,e are respectively the gas consumed and electricity produced by the gas turbine at time t; P t eh,e and P t eh,h are respectively the electricity consumed and heat produced by the electric boiler at time t; P t Specifically, the input of the absorption chiller and waste heat boiler comes from the waste heat generated by the gas turbine: and P t ec,c are respectively the electricity consumed and cold produced by the electric chiller at time t; P t w is the auxiliary power consumption at time t; P t gh,g and P t gh,h are respectively the gas consumed and heat produced by the gas boiler at time t; P t hc,c is the cold output of the absorption chiller at time t; P t whb,o is the heat output of the waste heat boiler at time t.h ≤ k ge,h P t ge,e , ∀t ∈ [0, T]

Table 1 .
Different types of devices.

Table 5 .
Operation Results in Normal Operation.

Table 6 .
Operation Results in Emergency.
Figure 12.Detailed dispatching results of the IES in emergency.
Lower/upper bound of the power output of device i S i /S i Lower/upper bound of the state of charge of storage device i η i Conversion efficiency of energy coupling device i η i,s /η i,r Charging/discharging efficiency of energy storage device i i State of charge of storage device i at time t