Review on Energy Resilience
Selection and peer-review under responsibility of the scientific committee of CUE2019 Copyright © 2019 CUE
Applied Energy Symposium 2019: Low carbon cities and urban energy systems October 16-18, 2019, Xiamen, China
Paper ID: 0088
QUANTIFYING THE IMPACTS OF WINDSTORMS ON RESILIENCE OF URBAN ENERGY SYSTEMS
Minglei Bao1, Yang Yang1, Kexin Wang2, Yi Ding 1*, Xiaoming Zhou1
1 College of Electrical Engineering, Zhejiang University, Hangzhou, China
2 Economic & Technology Research Institute of State Grid Shandong Electric Power Company, Jinan, China
ABSTRACT Faced with ever-increasing environmental pollution
and public concerns about energy security, urban energy systems (UESs) need to be constructed to improve the efficiency and reliability of energy utilization. However, due to significant impacts of extreme weather on the operation of UESs, it is important to develop a model to evaluate the resilience of UESs. In this paper, a synthetic model is proposed to quantify the impacts of windstorms on the resilience of UESs. Firstly, the optimization model of UESs under contingency states is developed to determine the generation re-dispatch and load shedding. Moreover, considering the effects of high winds, the line fragility model is utilized to calculate the unavailability of power lines according to the surrounding wind speeds. On the basis, the simulation framework for the resilience evaluation of UESs is developed utilizing Monte Carlo Simulation (MCS) technique to quantify the impacts of high winds on the UESs. Finally, the proposed methods are validated using the urban energy test systems.
Keywords: resilience; urban energy systems, High Winds, line fragility, Monte Carlo Simulation
NONMENCLATURE
Abbreviations
UESs urban energy systems MCS Monte Carlo Simulation CHP combined heat and power plants
Symbols
,i j subscript of nodes in electric system ,m n subscript of nodes in gas system
l superscript of contingency state
ig C generation cost of generator g at
node i l
ig P generation output of generator g
at node i
iL C compensation cost of load
curtailments
, l l
iL iL P Q active and reactive load curtailments
at node
ms C production cost of gas source at
node m l
ms W gas production of gas source s at
node m 0 0 ,
iL iL P Q active and reactive load demand at
normal state 0 l
i phase angle of voltage at node i
l
ij Y element of admittance matrix
l
i V voltage of electric node i
min max ,
i i V V min/max voltage at node i
min max ,
ig ig P P min/max generator output at node
i max max
, iL iL
P Q max active/reactive load shedding l
ij S power flow on line ij
max
ij S limits of power flow on line ij
l
m gas pressure at node m
l
mn f gas flow on pipeline mn
mn M constant pipeline flow coefficient
min max ,
m m min/max gas pressure at node m
min max ,
mn mn f f min/max gas flow on pipeline mn
L
l
m W gas load at node m
min max ,
ms ms W W min/max gas production of gas
source
2 Copyright © 2019 CUE
M number of nodes in gas system N number of nodes in electricity system
1. INTRODUCTION The increasing environmental pollution and energy
crisis entail the need for the construction of future low- carbon cities. As important “life floods” to low-carbon cities, urban energy systems (UESs) could provide clean, efficient and sustainable energy to guarantee the rapid development of society [1]. Through operational optimization of UESs, different energy vectors (heat, gas and electricity) could be converted efficiently to supply multi-energy to consumers.
Recently, affected by climate change (high temperatures, frequent windstorms and heavy rains), the resilient operation of UESs under severe weather can be threatened [2]. It is reported that the total annual economic loss of distribution electricity systems caused by severe weather ranges from 20 to 50 billion dollars [3]. Due to the dramatic consequences of severe weather, it is urgent to evaluate the resilience of UESs under extreme weather to guarantee the reliable provision of energy to consumers.
Resilience is defined as the ability of systems to resist the possible disturbances and recover to normal operation, which has been widely investigated in ecological, social and engineering systems [4]. In terms of power systems, numerous efforts have been paid to quantifying and modelling the impacts of severe weather on the system operation. Reference [2] proposes a conceptual framework to analyze the resilience of power systems caused by windstorms. In reference [3], the operational resilience and infrastructure resilience of power systems are investigated using time-dependent resilience metrics. In reference [4], a three-stage analysis framework is developed to investigate the change of system performance during resistive, absorptive and restorative stages. The previous researches mainly focus on the resilience evaluation of power systems, without considering the integration between different energy systems [5]. Faced with the rapid development of UESs, additional efforts need to be devoted to the resilience evaluation of UESs.
In this paper, a synthetic model is developed to evaluate the resilience of UESs induced by extreme weather. Firstly, the energy hub is modelled to describe the coupling relationship between its input energy and output energy. On the basis, the operation model of UESs under contingency states caused by component failures is proposed to determine the generation re-dispatch and
load curtailments. Moreover, the fragility model of power lines under high winds is developed to characterize the relationship between line unavailability and wind speeds. Based on the fragility model of power lines, the simulation framework for resilience analysis of UESs is developed using Monte Carlo Simulation (MCS) technique. Case studies illustrate the validity of the proposed model.
2. MODEL OF URBAN ENERGY SYSTEMS The traditional UESs consist of distribution electric
system, natural gas system and energy hub, as shown in Fig.1 [1]. The UESs obtain electricity and gas from common-utilities to satisfy the multi-energy demands of consumers. The energy conversion between different energy vectors could be realized through the energy hub.
2.1 Energy hub model
In this paper, the energy hub depicted in Fig.2 consists of a power transformer, a heat pump and combined heat and power plants (CHP) [1]. The input energy consists of electricity and gas, while the output energy consists of electricity heat. The energy hub could optimize the energy conversion between heat pump and CHP to improve system efficiency.
The coupling relationship between input energy and
output energy could be expressed as:
Fig 1. Description of urban energy systems.
Distribution electricity system
Natural gas system
Energy hub Energy
consumers Electric
load
Heat
load
Gas
load
H L
Power transformer
Heat pump
CHP
Electricity
Heat
Electricity
Gas
PEH
FEH
aPEH
(1-a)PEH
Le
Lh
Fig 2. Structure of energy hub.
3 Copyright © 2019 CUE
1
ee EH
h h EH
L P
L F
(1)
where represents the dispatch factor of electricity flow;
e and
h represent the electrical efficiency and
thermal efficiency of CHP, respectively; represents
heat pump efficiency; e
L and h
L represent the
electricity demand and heat demand of energy hub, respectively;
EH P and
EH F represent the input
electricity flow and gas flow of energy hub.
2.2 UES model under contingency states
Considering the impacts of high winds, the component failures can make the UESs deviate from its normal operation state. In order to maintain the reliable operation of UESs, measures such as load shedding or generation re-dispatch will be adopted [6, 9]. For contingency state l , the load shed and generation re- dispatched could be determined with optimization model which is described in (2)-(16). The objective function is to minimize total system cost including energy cost and load shedding cost.
1 1
min ( ) N M
l l l
ig ig iL iL ms ms
i m
C P C P C W
(2)
The objective function (2) is subject to the following constraints:
0 , 1
cos( ) N
l l l l l l l l
ig iL iL EH i i j ij i j
j
P P P P V V Y
(3)
0 1
sin( ) N
l l l l l l l
ig iL iL i j ij i j
j
Q Q Q V V Y
(4)
min maxl
i i i V V V (5)
min maxl
ig ig ig P P P (6)
min maxl
ig ig ig Q Q Q (7)
max
0 l
iL iL P P (8)
max
0 l
iL iL Q Q (9)
max
0 l
ij ij S S (10)
,
, 1
l
EH iee
l h h EH i
PL
L F
(11)
2 2 2
sgn( , ) l l l l l
m n mn mn m n f M
(12)
min maxl
m m m (13)
min maxl
mn mn mn f f f (14)
L
1
0 M
l l l
ms m mn
m
W W f
(15)
min maxl
ms ms ms W W W (16)
where Eq.(3) is active power balance equation; Eq.(4) is reactive power balance equation; Eq.(5) is bus voltage constraint of electric distribution system; Eq.(6) is active power output constraints of generators; Eq.(7) is reactive power output constraints of generators; Eq.(8) and Eq.(9) are electric load curtailment limits; Eq.(10) is power flow constraints of power lines; Eq.(11) is coupling relationship between input and output energy of energy hub; Eq.(12) is pipeline flow equation. Eq.(13) is nodal pressure constraints. Eq.(14) is gas pipeline flow constraints. Eq.(15) is nodal gas flow balance equation. Eq.(16) is gas source production constraints.
3. RESILIENCE EVALUATION OF UESS A synthetic model can be developed to evaluate the
impacts of high winds on the resilience of UESs. Considering gas pipelines are usually buried underground, only the impacts of high winds on distribution electricity system are modelled in this paper. Due to the integration between electricity distribution system and natural gas system, the failures in electricity system due to high winds can also have high impacts on the gas system. If the line failures lead to the power imbalance in the electricity system, the production of gas will be re-dispatched to adjust the power output of CHP. Therefore, although there are no component failures in natural gas system, the operation condition of gas system can also be adjusted considering the integration between different energy vectors. In this section, the modelling of line fragility and resilience evaluation of UESs to windstorms are illustrated.
3.1 Line fragility model
Considering the impacts of high winds, the fragility of power lines is closely related to their surrounding wind speeds. According to [2], the unavailability of a line is a piecewise function of wind speeds, which can be expressed as:
0
_
, if
( ) if
1, if
ij critical
ij ij v critical collasape
critical
p v v
p v p v v v v
v v
, (17)
where ijp v denotes the unavailability of power line
ij as the function of wind speed of v ; 0 ij
p denotes the
unavailability of line for good weather state; _ij vp v
4 Copyright © 2019 CUE
denotes the linear relation between unavailability and wind speeds between
critical v and
collasape v .
3.2 Simulation framework for resilience analysis of UESs
As illustrated in Fig.3, the process for the resilience evaluation of UESs can be divided into three steps. The first step is to initialize system parameters, including the data of distribution electricity system, gas system and energy hub.
The second step is the resilience analysis of UESs
using MCS technique [10, 11]. Firstly, the parameters about wind speeds in different regions are set, based on which the unavailability of power lines can be determined. According to the unavailability induced by high winds, the failure lines could be determined through random sampling. For certain failures, the optimization model in (2)-(15) can be applied to determine the load shedding and re-dispatch of generation. Repeat the previous procedures until the stopping criteria of MCS technique are satisfied.
Based on the load curtailments for each simulation time, the resilience indices of UESs can be calculated. The resilience indices of UESs are defined as the amount of electric load demand which are still connected after high winds, which can be expressed as:
min maxl
ms ms ms W W W (18)
where t and T represent the simulation time and the total sampling number of MCS technique.
4. CASE STUDIES AND DISCUSSIONS The test system in this paper is composed of IEEE 33-
node distribution system and the modified 7-node gas system. The network topology and the parameters of test system can be found in [1]. The distribution energy resources in electricity system are assumed to be storage devices. The electrical efficiency and thermal efficiency of CHP are assumed to be 0.4 and 0.45 [7]. The heat pump efficiency is assumed to be 3 [8]. The line fragility curve is determined according to [3], which is shown in Fig.4. In order to evaluate the impacts of windstorms on the resilience of UESs, three cases with different winds speeds are analyzed, including 31m/s, 33m/s, 36m/s and 40m/s wind speeds. Moreover, the resilience of UESs under different wind speeds are compared with that of IEEE 33-node distribution system. The simulation results are given in Fig.5.
It can be noted from Fig.5 that the increase of wind
speeds can reduce the resilience of UESs. This is mainly because that the more power lines will be attacked with the increase of wind speeds and more electric load will be cut for the reliable operation of UESs. The resilience of UESs is for 31m/s, 33m/s, 36m/s and 40m/s wind speeds, respectively. Besides, it can be noted that UESs are more resilient than distribution electricity system for the same wind speeds. This is mainly because that the gas system can provide electricity to distribution electricity system through CHP for contingency states.
Initialize system
parameters
Set parameters, t = 1
Random sampling
Update system
parameters
Run optimization
model
Record system
information
Meeting stopping
criteria of MCS
Calculate resilience
indices
t = t+1
Y
N
Fig 3. Simulation framework for resilience evaluation of urban energy systems
Fig 4. Relationship between unavailability of lines and winds speeds
5 Copyright © 2019 CUE
5. CONCLUSIONS This paper proposes a synthetic model to quantify
the impacts of windstorms on the resilience of UESs. Firstly, the operation model of UESs under contingency states is proposed to determine the generation re- dispatch and load curtailments. Moreover, the simulation framework for resilience analysis of UESs is developed using MCS technique. The simulation results show that the increase of wind speeds can reduce the resilience of UESs. Besides, the UESs can be more resilient to windstorms than distribution electricity system due to the support of gas system for contingency states.
The studies in this paper mainly focus on the impacts of windstorms on the resilience of UESs composed of power, gas and heat systems. In the future work, the resilience of different energy systems (e.g. water system) under different extreme weather (e.g. earthquake) can be evaluated.
ACKNOWLEDGEMENT The research is supported by the China NSFC under
Grant 71871200 and National Natural Science Foundation China and Joint Programming Initiative Urban Europe Call (NSFC– JPI UE) under grant 71961137004.
REFERENCE [1] X. Jin, Y. Mu, H. Jia, J. Wu, X. Xu, and X. Yu, "Optimal day-ahead scheduling of integrated urban energy systems," Applied energy, vol. 180, pp. 1-13, 2016. [2] M. Panteli and P. Mancarella, "Modeling and evaluating the resilience of critical electrical power infrastructure to extreme weather events," IEEE Systems Journal, vol. 11, pp. 1733-1742, 2015. [3] M. Panteli, P. Mancarella, D. N. Trakas, E. Kyriakides, and N. D. Hatziargyriou, "Metrics and quantification of
operational and infrastructure resilience in power systems," IEEE Transactions on Power Systems, vol. 32, pp. 4732-4742, 2017. [4] M. Ouyang, L. Dueñas-Osorio, and X. Min, "A three- stage resilience analysis framework for urban infrastructure systems," Structural safety, vol. 36, pp. 23- 31, 2012. [5] M. Bao, Y. Ding, C. Singh, and C. Shao, "A Multi-State Model for Reliability Assessment of Integrated Gas and Power Systems Utilizing Universal Generating Function Techniques," IEEE Transactions on Smart Grid, 2019. [6] H. Zhou, M. Bao, B. Zhu, Z. Xiang, L. Xu, Z. Gao, et al., "A strategy for contingency management with aggregated thermostatically controlled loads in china’s electricity market," in IOP Conference Series: Earth and Environmental Science, 2019, p. 032010. [7] G. Chicco and P. Mancarella, "Matrix modelling of small-scale trigeneration systems and application to operational optimization," Energy, vol. 34, pp. 261-273, 2009. [8] C. Shao, Y. Ding, J. Wang, and Y. Song, "Modeling and Integration of Flexible Demand in Heat and Electricity Integrated Energy System," IEEE Transactions on Sustainable Energy, vol. 9, pp. 361-370, 2018. [9] P. Wang, Y. Ding, and Y. Xiao, "Technique to evaluate nodal reliability indices and nodal prices of restructured power systems," IEE Proceedings Generation, Transmission and Distribution, vol. 152, pp. 390-396, 2005. [10] Y. Ding, C. Singh, L. Goel, J. Ø stergaard, and P. Wang, "Short-Term and Medium-Term Reliability Evaluation for Power Systems With High Penetration of Wind Power," IEEE Transactions on Sustainable Energy, vol. 5, pp. 896- 906, 2014. [11] M. Chaudry, J. Wu, and N. Jenkins, "A sequential Monte Carlo model of the combined GB gas and electricity network," Energy Policy, vol. 62, pp. 473-483, 2013.
Fig 5. Simulation results of UESs and distribution electricity system for different wind speeds