The resilience pre-disaster assessment stage consists of a “toughness” assessment, that is, the expected capability of power systems to keep stable operation and provide continued power supply for users against disasters. The core of toughness assessment is to analyze the possible damage and losses due to specific disasters in advance, which is the basis for disaster hazard prediction and loss mitigation.
Toughness, as an inherent performance of resilient power systems, is primarily determined by grid structure, equipment strength and smart operational measures. However, the duration time, characteristics of occurrence and damage mechanism to components of different high impact/low probability (HILP) events vary greatly, so even the toughness of the same power system will act differently under different circumstances. Concrete toughness analysis is needed under specific disasters.
2.2.1. The Metrics of Toughness
This paper proposes some metrics to reflect the toughness of power system from three aspects: the expected damage degree of system, the risk of system splitting, along with the generation and transmission margins.
● The expected damage degree of system
When the system is hit by disasters, the direct impact is the load loss due to the breakdown of components, which emphasizes the necessity of evaluating the expected load and component loss situation against disasters.
The loss of load probability (LOLP) and the expectation of demand not supplied (EDNS) are suitable for the assessment of toughness, and the definition of them can be found in [
22].
Average load loss ratio in case of load curtailment (ALRIL), represents the average proportion of load loss to the whole current load demand, which is:
where
Ns indicates the whole simulation times, and
Llc,i represents the amount of load shedding in the
ith simulation (the simulation method is presented in the next part).
Ld,i is the load demand of system in the
ith simulation.
The probability of components failed during disaster (PCFD) is:
Nc,i is the number of the failed components in the
ith simulation, and
Nc indicates the number of the whole components. The noteworthy thing is that PCFD should be calculated according to the types of vulnerable components separately.
● The risk of system splitting
When a large-scale disaster happens, the risk of system splitting or some nodes being isolated will increase, as components faults may create isolated parts in the system. The following two metrics are introduced to address these possibilities in the evaluation framework.
The possibility of system islanding (PSI) and isolated note (PIN), as in Equations (3) and (4) respectively:
Nsi is the number of simulations where system islanding occurs and
Nin is the number of simulations where isolated node exists among the whole simulation times
Ns.
● The generation and transmission margins
Power systems can be divided into several areas with close electrical connections, according to the impacted areas under disasters and the electrical characteristics of the system. Usually, the neighboring areas are linked by several transmission lines, which make up the transmission section. Currently, it is widely-used practice to divide the post-disaster system into several parts, and each part is restored in parallel. This parallel restoration adopts a simpler restoration strategy at each area while accelerating the restoration process, and can guarantee system stability.
Therefore, the generation margins of areas (GMA) and the transmission margin of transmission section (TMTS) between different areas are proposed, to reflect the expected maximum generation capacity within the area and the expected transmission capacity of the section against disasters:
Gg,max The maximum real power output of unit g;
Lb,real,u The real load amount on bus b under disasters;
Lb,real,be The real load demand on bus b before disasters;
Gnow The set of normal units within the area under disasters;
Gall The set of normal units within the area before disasters;
Ball The set of all buses within the area;
Ck,max The power transmission upper limit of line k;
Ck,u The power flowing on line k under disasters;
Ck,be The power flowing on line k before disasters;
lnow The set of normal lines among the transmission section under disasters;
lall The set of normal lines among the transmission section before disasters.
If the value of generation margin is small or even negative, it indicates that there will be a serious shortage of power generation within the area, and this will lead to higher internal load shedding risk once the transmission lines are damaged. Conversely, large value of generation margin means the sufficient regional power generation, which will help to avoid load loss.
If the value of transmission margin is small, it indicates that the transmission capacity of the section will be restricted dramatically, so the power demand of different areas is mainly supplied by their internal generation. By contrast, the large transmission margin means that the transmission capacity of the section is not badly disrupted.
2.2.2. Toughness Analysis Based on PMRM
Under disaster conditions, the probability of high-loss and medium-loss damage situations increases greatly, so the average-risk-based assessment metrics cannot fully reflect the severity of serious damage. PMRM, which was first proposed by Acbeck and Haimes in 1984 [
23], helps to show high and medium risks by analyzing the conditional probabilities and the conditional risk expected values (CREV). Hence, PMRM is also applicable to toughness analysis.
To illustrate the practical value of PMRM, it is used to analyze the toughness of power systems, by taking “EDNS” for example. The following two metrics are proposed according to the proportion of expected load loss (the value of EDNS) to the system load demand. We use the load loss ratios of 5% and 25% as the segmentation points. The load loss intervals of low loss, medium loss and high loss are designated as 0~5%, 5~25% and 25~100% separately.
Ri represents the expected load loss under low-loss, medium-loss and high-loss situations:
Pi represents the conditional risk probability under low-loss, medium-loss and high-loss situations:
SL,
SM, and
SH are the sets of low-loss, medium-loss and high-loss events.
P(
j) is the occurrence probability of the load shedding event
j, and the corresponding amount of load curtailed is
Loss(
j).
pL,
pM, and
pH are the occurrence probabilities of low-loss, medium-loss and high-loss events.
When the CREV or high-loss and medium-loss conditional probabilities exceed the warning boundaries set by decision makers, it is hard for power systems to resist a given disaster. Under this circumstance, it is urgent to take measures to reinforce and update the system, which effectively mitigate the serious harm and reduce the risk of cascading failures.
2.2.3. The Procedure of Toughness Assessment
Taking windstorms and lightning as examples, the toughness assessment procedure under different types of disasters is introduced. This paper refers to the idea of reliability evaluation and risk assessment, and divides the procedure of toughness assessment into several parts.
● The identification of disasters characteristics
To accurately analyze the toughness of power systems, we should make sure which type of disaster is to be considered. The probable affected region, and the duration and intensity of disasters should be identified at first. In this paper, windstorms and lighting are discussed to evaluate the resilience of the power system.
● The acquisition of the type and the failure/repair rates of vulnerable components
The types of vulnerable components vary greatly under different disasters. For example, cyber- attacks will impact communication and information storage devices (such as EMS and SCADA); earthquakes will presumably destroy substations; transmission towers and lines tend to collapse during high winds. By identifying the different disasters, we can determine the temporal and spatial damage mechanism on power system of disasters, identify the type of vulnerable components and then acquire their failure/repair rates.
It is widely accepted that the quadratic model is suitable to describe the pressure exerted on transmission lines by winds. The failure rates of transmission lines are formulated to be proportional to the square of the wind speed [
24,
25]:
when the wind speed is below
ωcrit,
λwind =
λnorm.
The average repair time is regarded to be constant for wind speeds below
ωcrit [
24] and increases linearly for higher wind speeds. The repair time for transmission lines under windstorms can be defined as:
where
ω(t) is the wind speed at time
t;
α is a scaling parameter;
ωcrit is the critical wind speed over which
λ is increased;
λnorm is the constant failure rate under normal condition;
rtnorm is the reference repair time under normal condition and
k is the repair time growth rate with
ω(t).
In order to compare the effects of different disasters on the toughness of power systems, the lightning is analyzed as well. The line failure rate under lighting strike conditions is obtained from [
26], as shown in Equation (11), where a linear relationship is used, considering the ground flash density:
Based on what is suggested by Cadini [
26], the repair rate model of transmission lines is built as:
where
rrnorm is the sonstant repair rate during normal conditions;
μ,
β are scaling parameters; and
Ng is the number of ground-flashes per unit of time and area (
occ/
h*
km2).
● State sampling of components
Power systems are traditionally designed to sustain single outage contingency (“N-1” security criterion). However, multiple and large-scale component faults may happen during HILP events.
Generally, to mitigate the outage, faulty component repair work should be carried out, especially for long-duration disasters, such as windstorms. Therefore, the procedure of failure-repair-failure needs to be taken into account. SMCS is suitable to sample the states of vulnerable components.
As for short-duration event, such as earthquake or lightning, the repair work is usually not assigned until the end of disasters. Sometimes disasters (like floods) may be so severe that the repair work cannot be implemented for the safety of repair crews. In this case, Nonsequential Monte Carlo Simulation is adopted to sample the state of components, allowing the failure possibility of vulnerable components to be determined
Due to the lack of related mathematical models, we combine Equations (11) and (12) to establish the failure possibility model of transmission lines from a statistical point:
● State analysis of power system
After obtaining the state of each component, the power system state should be analyzed. To reduce the losses of outage, this paper takes the minimum amount of load shedding as the objective and uses DC OPF to solve the problem if necessary. Then the metrics proposed in
Section 2.2.1 are calculated.