Next Article in Journal
Experimental Investigation of Micro Cooling Units on Impingement Jet Array Flow Pressure Loss and Heat Transfer Characteristics
Next Article in Special Issue
Spatio-Temporal and Power–Energy Scheduling of Mobile Battery Storage for Mitigating Wind and Solar Energy Curtailment in Distribution Networks
Previous Article in Journal
Comparative Study of Stepwise Optimization and Global Optimization on a Nine-Phase Flux-Switching PM Generator
Previous Article in Special Issue
Techno-Economic Planning and Operation of the Microgrid Considering Real-Time Pricing Demand Response Program
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Augmented Power Dispatch for Resilient Operation through Controllable Series Compensation and N-1-1 Contingency Assessment

1
Department of Electrical Engineering, School of Automation, Guangdong University of Technology, Guangzhou 510006, China
2
Brunel Interdisciplinary Power Systems Research Centre, Department of Electronic and Electrical Engineering, Brunel University London, London UB8 3PH, UK
3
Center for Electric Power and Energy, Department of Electrical Engineering, Technical University of Denmark, 2800 Kongens Lyngby, Denmark
*
Authors to whom correspondence should be addressed.
Energies 2021, 14(16), 4756; https://doi.org/10.3390/en14164756
Submission received: 24 June 2021 / Revised: 28 July 2021 / Accepted: 3 August 2021 / Published: 5 August 2021
(This article belongs to the Collection Feature Papers in Smart Grids and Microgrids)

Abstract

:
Research on enhancing power system resilience against extreme events is attracting significant attention and becoming a top global agenda. In this paper, a preventive augmented power dispatch model is proposed to provide a resilient operation. In the proposed model, a new N-1-1 security criterion is proposed to select disruptive N-1-1 contingency cases that might trigger cascading blackouts, and an iterative contingency assessment process based on the line outage distribution factor is proposed to deal with security constraints. In terms of optimization objectives, two objectives related to power flow on the transmission line are considered to reduce the possibility of overload outages. Controllable series compensation devices are also considered in the model to improve the power flow distribution. Case studies conducted on the modified IEEE 30-bus, 118-bus and Polish 2382-bus systems show that the power flow solution of the proposed power dispatch model can avoid some branches from undertaking excessively heavy loads, especially lines forecasted to be affected by extreme events. The results of blackout simulations through a hidden failure cascading outage simulation model show that the average power losses of the proposed model are reduced by around 40% in some cases as compared to the classical economic dispatch model.

1. Introduction

1.1. Background and Motivation

Blackouts caused by extreme weather events have been occurring with increasing frequency in the past few decades and have caused considerable economic losses to electric power utilities and society in general [1,2,3]. Research on enhancing power system resilience against these extreme events is attracting significant attention and becoming a top global agenda [4,5]. Resilience is defined as the ability of a power system to withstand extreme events, adapt its operation and structure for preventing or mitigating their impact, and recover quickly after its exposition to such events [6]. Generally, approaches to enhance power system resilience can be summarized into two categories, namely infrastructure hardening measures and operational resilience enhancement strategies [7,8]. Infrastructure hardening measures might be more effective than the operational ones but are less economical and hard to implement [9]. Operational strategies could take advantage of incoming information about extreme events, applying specific actions that are adaptive to the prevailing conditions [10]. This paper focuses on enhancing the operational resilience of a power grid and proposes an augmented power dispatch model to offer a preventive resilient generation dispatch solution for power system operators.

1.2. Related Works and Research Gap

The contribution of preventive scheduling of power generation resources to system resilience has been well recognized [11]. Several studies have been carried out based on the classical unit commitment and economic dispatch problems. In [12], a proactive resource allocation model based on a stochastic unit commitment problem was proposed to minimize the restoration cost and improve the operational resilience of electrics power grids against hurricanes. Reference [13] proposed a resilience-constrained unit commitment model considering the impact of microgrids in which contingency scenarios with higher weather-related failure probability were considered. In general, microgrids can keep the system stable operation and minimize blackout. Reference [14] applied the mean-variance Markowitz theory and value-at-risk to assess the risk of system cost variability due to the presence of high penetration of renewables such as PV. The risk of system cost can further be reduced through energy storage participation. Reference [15] proposed a day-ahead two-stage stochastic unit commitment optimization framework to offer a preventive operation strategy to power system operators. Reference [16] presented a day-ahead and intraday multi-objective dispatch model for an integrated biomass concentrated solar system considering the conditional value-at-risk. Reference [17] developed a sequentially proactive operation strategy in which a Markov process was used to model the uncertain sequential transition of system states due to damages caused by extreme events. The transition probabilities were evaluated according to weather-dependent failure rates. Our previous work [18] proposed a resilience-constrained unit commitment framework in which a proportional hazard model was used to calculate the forced outage rate of transmission lines affected by extreme weather events. Reference [19] proposed a two-stage robust unit commitment taking into account the weather-related failure probability of transmission lines.
The main difference between the above resilience-based generation dispatch optimization models and the traditional ones is the considered contingency set. Contingency cases with higher failure probability related to extreme events are considered in the security constraints in resilience-based models. Fragility curves of system components in different extreme events are considered to model the impact of the forecasted extreme event on component outages in the above papers. However, fragility curves presented in the literature are difficult to implement in practice since it is hard to obtain accurate parameters of the curve. Moreover, it is also difficult to forecast or measure the specific impact of extreme events on every electrical component. For example, it is not easy to get the wind speed and angle acting on a transmission line during a hurricane. What we can know from current extreme weather events prediction and waring service system in power grids are the areas or components that will be affected and the general hazard level for the event. This gives challenges in real operation to access the efficacy of the above methods. Another problem of the existing works is that most of the studies only focus on changing the security constraint to provide a resilient operation strategy. The security constraints only ensure that the operation is safe for the defined contingency set. The system might be violated if other high-order contingency cases occur.

1.3. Contributions of This Paper

In this paper, both security constraints and the objective function of the optimization model are modified to offer a resilient generation dispatch solution. As for the security constraints, two types of contingency sets are considered. The first type is the N-1 contingency of all transmission lines. The second type is the proposed N-1-1 contingency set. Given that heavily loading lines have a high probability of overload outages or hidden failures [20], transmission lines with heavy post-contingency power flow of N-1 outage are selected as N-1-1 contingency sets in this work. By considering the security constraint of these N-1-1 contingency cases, the further propagation of N-1 contingency cases caused by the extreme event might be reduced. It can be imagined that the number of the defined N-1-1 contingency cases is affected by the generation dispatch solution. In turn, the dispatch solution may be changed after considering the contingency constraints. Therefore, an iterative process is needed to address this issue. The iterative contingency filtering process proposed in [21] is extended in this paper to generate the defined N-1-1 contingency cases and incorporate the security constraints of the N-1 and N-1-1 contingency sets into the optimization model.
The security constraints mentioned above only ensure that the operation is safe for the defined contingency set. The system might be violated if other high-order contingency cases occur (for example, the N-3 contingency case). Of course, we can avoid this situation by considering higher-order security constraints in the optimization model, but considering higher-order security constraints may lead to unsolvable problems.
In this paper, we try to cope with this situation by adjusting the optimization criteria of the objective function. A key reason why initial failures of a few components caused by extreme events develop into cascading failures that eventually lead to blackouts is the massive power flow transfer in the grid after the initial failure [22]. If we can improve the power flow distribution in the network to prevent massive power flow transfer after disturbances, the probability of blackouts caused by cascading failures could be reduced [23]. Referring to our previous study [24], two objectives related to power flow are considered in this paper to improve the power flow distribution in the transmission system. However, improving the distribution of power flow by simply adjusting the generation of generators is limited in some cases because the line flow is limited by the impedance and capacity of the transmission line. In this paper, the adjustment of transmission line reactance by controllable series compensation devices is considered in the optimization model. Controlled series compensation devices are already part of the transmission network in many countries, and it is expected that the number of controlled series compensation devices installed in the grid will increase significantly in the future with advances in control methods [25]. This paper focuses on devices that compensate for the reactance of the line they are installed in [26], such as static synchronous series compensators (SSSC) and thyristor-controlled series capacitors (TCSC). The main difficulty in incorporating series compensation devices into the DC power flow-based security-constrained economic dispatch problem is the computational complexity caused by the flexibility of line reactance. Reference [27] reformulated this nonlinear programming problem as a mixed-integer linear programming problem by using the big M relaxation technique. A more computationally efficient method was proposed in [28] based on the PTDF-based power flow model. Referring to the above study, this paper jointly optimizes the generator output and the reactance adjustment value of the transmission line by a controllable series compensation device to improve the power flow distribution in the transmission system. The main contributions of this paper are summarized below:
  • A new N-1-1 security criterion is defined to select the disruptive contingency cases that might trigger cascading failures. It is defined from the perspective of preventing further propagation of N-1 contingency. The security constraints of the defined contingency set are formulated by using the line outage distribution factors and are taken into account in the proposed power dispatch model through an iterative contingency filtering process.
  • Two objectives related to active power flow on transmission lines are considered in the dispatch model to avoid full or heavy loads on lines in the transmission system, thus reducing the probability of massive power flow transfer and overload cascading outages after the initial outage. The adjustment of transmission line reactance by controllable series compensation devices is considered in the optimization model.
  • The proposed augmented power dispatch model is nonlinear due to the line flow-related objective function and the consideration of controllable series compensation devices. Linear relaxation techniques are introduced to convert the model into a mixed-integer linear program. Although MILP is less challenging compared to the original nonlinear programming, it is not preferred, especially for large-scale real power systems. A computationally efficient two-stage solution is proposed to further reduce computational complexity.
The remainder of the paper is organized as follows: section two introduces the proposed augmented power dispatch optimization model and its solution methodology. Section three presents the simulation results of the case studies conducted on the IEEE 30-bus, 118-bus, and Polish 2383-bus systems. Section four presents the discussions. Conclusions are given in section five.

2. The Proposed Augmented Power Dispatch Model and Solution Methodology

2.1. Definition of the Defined N-1-1 Contingency Set

The classical security-constrained economic dispatch mainly considers the N-1 security constraints and is not resilient enough to ensure the safe operation of the grid. However, since there is a large number of transmission lines in the grid, it is impractical to have a solution that satisfies the security constraints for all N-k contingencies [29]. Contingency selection is often employed to reduce computational complexity.
As mentioned in the introduction, two types of contingency sets are considered in this paper. The first type is the N-1 contingency of all lines. The second type is the proposed N-1-1 contingency set. It is defined from the point of preventing the further propagation of N-1 contingency cases. Although the N-1 security constraints have already ensured that the loss of any single line does not result in power flow violation of the remaining online lines, there exist some lines whose post-contingency power flows exceed steady-state operating limits and are close to the emergency operating limits. These lines have a high probability of hidden failure and are also susceptible to cascading outages if disturbances occur. Therefore, in this work, these lines are considered as potentially disruptive N-1-1 contingency cases. By considering these contingency cases in the optimal power dispatch problem, the further propagation of N-1 contingency cases caused by an extreme event might be reduced.
Since the selection of the defined N-1-1 contingency cases is based on the analysis of the post-contingency power flow for each N-1 contingency case, the number of N-1-1 contingency cases is affected by the power dispatch solution. In turn, the solution may be changed after considering the security constraints of these contingency cases. We will first introduce the power dispatch optimization model in the next section and then the method for handling the security constraints of the defined contingency set.

2.2. Mathematical Formulation of the Proposed Power Dispatch Optimization Model

The security constraints of the defined contingency set only ensure that the power flow solution is secure for specific contingency cases. If other high-order contingency cases occur, such as the N-3 contingency case, the system might be violated. As mentioned in the introduction, the optimal criteria of the generation dispatch optimization model are modified in this paper to improve the resilience of the solution. Two objectives related to power flow are considered in this paper to reduce the probability of overload outage of N-k contingency cases. Firstly, the disaster forecast information is used to select the lines that will be affected by the upcoming extreme events. It can be seen that only the forecasted information about whether a component is affected by an extreme weather event is needed, which makes it more applicable for system operators. The first objective is to reduce the loading rate of lines that are forecast to be affected by extreme weather events, as these lines have a higher probability of outages. The second one minimizes the difference in the loading rate of all transmission lines. These two objectives interact with each other. By optimizing both objectives in the optimal power flow problem [30], not only the loading rate of the transmission line can be reduced, but also the overall utilization of the grid can be improved. By defining the loading rate r l as the ratio of the active power flow on a transmission line to its transmission capacity, the detailed mathematical representation for each objective is given below.
r l = | f l f l max |
F 1 = l S AL r l
F 2 = 1 N l S L | r l l S L r l N |
  F 3 = g S G ( a g p g 2 + b g p g + c g )
where f l is the active power flow of line l . f l max is the active power flow capacity of line l in the normal operating state. S AL , S L , and S G are set of transmission lines that are forecast to be affected by an extreme weather event, all lines, and generators. N is the number of the line. p g is the active power output of generator g . a g , b g , c g are quadratic, linear, and no-load cost coefficients of generator g , respectively.
The first objective F 1 is created based on the sum of the loading rate of the lines affected by an extreme event, and the second objective F 2 is created based on the mean absolute deviation of the loading rate of all transmission lines. The second objective reflects how close the loading rate of a transmission line is to the average loading rate. It provides a measure of the homogeneity of power flow distribution. The difference in the power flow on all lines is reduced when this measure becomes smaller. With the cooperation of these two objectives, the loading rate of every transmission line is optimized to be in equilibrium at a low level. The last objective F 3 is the total production cost of all generators, and the production cost of a generator is usually a quadratic function of generation. Considering the above three objectives and the adjustment of transmission line reactance through controllable series compensation devices, we propose an augmented power dispatch model to improve the power flow distribution in the transmission system and minimize generation costs, thus improving the resilience of the power system under extreme weather events, as follows:
min   α l S AL r l + β l S L | r l l S L r l N | + γ g S G ( a g p g 2 + b g p g + c g )
p g min p g p g max     g S G
g S G K b , g p g d b l S L A b , l f l = 0     b S B
f l = f l int + f l facts     l S L
f l int = b S B S F l , b p b i n j     l S L
p b inj = g S G K b , g p g d b l S L A b , l f l facts     b S B
δ l min f l int Mz l f l facts δ l max f l int + M z l     l S L
δ l max f l int M ( 1 z l ) f l facts δ l min f l int + M ( 1 z l )     l S L
M z l f l int M ( 1 z l )     l S L
z l { 0 , 1 }     l S L
δ l min = { Δ b l min / b l     l S facts 0   l S facts
δ l max = { Δ b l max / b l     l S facts 0   l S facts
f l max f l f l max     l S L
r l = | f l f l max |     l S L
where α , β , λ are coefficients of the considered three optimization objective terms. p g min , p g max are the minimum and maximum active power output of generator g . d b is the active power demand at the bus b . S B is the set of buses. K b , g is the element located in row b and column g of the bus-generator incidence matrix K . A b , l is the element located in row b and column l of bus-line incidence matrix A . f l int is the initial active power flow on the transmission line without considering the control of controllable series compensation. f l facts is the active power flow on the transmission line injected by the controllable series compensation. p b inj is the active power injected to bus b . δ l min , δ l max are the minimum and maximum relative adjustment of the susceptance of line l . b l is the initial susceptance of line l . Δ b l min , Δ b l max are the minimum and maximum adjustment of the susceptance of line l . S facts is the set of lines with controllable series compensation. z l represents the direction of power flow f l and takes a value of 0 for positive power flows and a value of 1 for negative power flows. S facts is the set of lines equipped with controllable series compensation devices. M is a sufficiently large positive number. S F l , b is the shift factor between line l and bus b , which denotes a change in the power flow on line l when one unit of power is injected to bus b . Since the classical DC power flow model is used in the proposed model, the following equation can calculate the shift factor matrix of a given transmission network topology.
SF = XA T B 1
B 1 = T T ( TBT T ) 1 T
B = AX 1 A T
where SF is the shift factor matrix with the size of ( N × N B ) . S F l , b is in row l and column b of matrix SF . X is a diagonal matrix with elements representing the reactance of all lines with the size of ( N × N ) . A is a bus-lines incidence matrix with the size of ( N B × N ) . T is a reduced identity matrix with the size of [ ( N B 1 ) × N B ] . N , and N B are the number of all lines and all buses.
Constraints in (6) limit the minimum and maximum output of all generators. Equations in (7) ensure power balance for each bus. Active power flow on a transmission line is calculated by constraints (8) to (16), and limited by (17). In this paper, reactance control of series compensation devices is modelled as power flow injection to keep the initial shift factor matrix of the transmission network without compensation devices unchanged, thus improving the computational efficiency. The detailed model is shown below. Assuming that the susceptance adjustment range of a line equipped with a series compensation device is [ Δ b l min , Δ b l max ] , the active power flow on the line can be calculated by separating the effect of the series compensation device and the initial power flow of the line as follows [28]:
f l = b l ( θ l , to θ l , from ) + Δ b l ( θ l , to θ l , flom ) = f l int + Δ b l b l f l int Δ b l min Δ b l Δ b l max δ l = Δ b l b l f l facts = δ l f l int
where θ l , to and θ l , from are the voltage angle of the from bus and to bus of line l , respectively. Δ b l are the susceptance adjustment of line l by the control of the series compensation in line l . δ l is the relative adjustment of the susceptance of line l .
The equation f l facts = δ l f l int in (22) includes the production of adjustment variable δ and line flow f l . It needs further reformulation to avoid computational complexity. By using the big M reformulation technique and variable z l representing the direction of line flow f l , this equation can be reformulated as a set of mixed-integer linear constraints shown, as in (11) to (16).
The modelling of an augmented power dispatch model with controllable series compensation devices has been completed. However, it is difficult to solve the optimization model directly due to the absolute value terms in the objective function in (5) and the constraints in (18).
According to optimization theory with absolute values, the objective and constraints with absolute values can be reformulated as linear or mixed-integer linear by introducing some auxiliary variables. First, to linearize (18), the binary variable z l used in (11) to (13) and two new continuous variables m l and n l are introduced to reformulate each absolute value constraint as a set of linear constraints (23) to (26). To linearize the second objective term in (5), only two additional constraints (27) and (28), and one auxiliary variable u l , are needed for each term. As for the quadratic cost functions, they can be approximated by a set of piecewise linear functions along with constraints (29) and (30). Thus, the proposed nonlinear problem can be reformulated as a mixed-integer linear programming problem (MILP) as follows:
  min   α l S AL r l + β l S L u l + γ g S G ( F 3 ( p g min ) + I F g , q p x g , q )
s.t. (6) to (17)
f l f l max = m l n l     l S L
r l = m l + n l     l S L
0 m l 1 z l     l S L
0 n l z l     l S L
r l l S l r l N l u l     l S L
l S l r l N l r l u l     l S L
p g = p g min + q S q p x g , q     g S G
0 p x g , q p i g     g S G , q S G
where,
p i g = p g max p g min N Q     g S G
I F g , q = F 3 ( p g , q ) F 3 ( p g , q 1 ) p i g     g S G , q S Q
{ p g , q = p g min + 1 q p i g     g S G , q S Q p g , 0 = p g min     g S G
where p x g , q is the active power output of generator g at segment q . p i g is the maximum active power output of generator g at segment q .
Although MILP is less challenging compared to the original nonlinear program problems, it is not preferred, especially for large-scale real power systems. As mentioned before, the binary variables in the proposed problem determine the flow direction of the lines. Assuming that the power flow direction of each line is known, the binary variables can be fixed and the proposed problem will become a linear problem. Based on the above analysis, a computationally efficient two-stage solution method is introduced for large-scale real power systems. The first stage solves the classical economic dispatch model without the two proposed objectives and controllable series compensation devices and initializes the direction of the power flows for the second stage. The second stage enforces the power flows in the initial direction of the first stage; that is, the binary variables z l in (11) to (13) and (25) to (26) are fixed. The proposed power dispatch model becomes a linear problem. It can be easily solved. Though the solution may or may not be globally optimal, the significance of this method does not vanish in the absence of optimality; the algorithm can very quickly find a solution with a significantly lower cost and line loading rate.

2.3. Iterative Contingency Filtering Process

As described in Section 2.1, the defined N-1-1 contingency set is created after getting the optimal power flow dispatch solution, and the security constraints in turn affect the solution of the optimization model. Therefore, the iterative contingency filtering process proposed in [21] is extended in this paper to generate the defined N-1-1 contingency cases and incorporate the security constraints of the N-1 and N-1-1 contingency sets into the optimization model. The line outage distribution factor is used to calculate the post-contingency power flow on each online line as well as to create security constraints for active N-1 contingency and N-1-1 contingencies cases in the process. The post-contingency power flow is calculated by using the line outage distribution factor (LODF) of the DC power flow. The LODF is defined as the incremental active power flow on the monitored transmission line caused by the outage of a contingent line with a pre-contingency active power flow of one unit [31]. LODF under a single-line outage or multi-line outage can be generalized by using the power transfer distribution factor (PTDF) of a pre-contingency network with the following equations [32,33].
LODF = PTDF M , O 0 ( E PTDF O , O 0 ) 1
PTDF M , O 0 = X M 1 Φ T   B 1   Ψ
PTDF O , O 0 = X O 1 Ψ T B 1 Ψ
where LODF is a line outage distribution factor matrix with the size of ( N M × N O ) . PTDF M , O 0 and PTDF O , O 0 are power transfer distribution factor matrices with the size of ( N M × N O ) and ( N O × N O ) , respectively. E is an identity matrix with the size of ( N O × N O ) . X M and X O are diagonal matrices with elements representing the reactance of monitored lines and outaged lines with the size of ( N M × N M ) and ( N O × N O ) , respectively. Φ and Ψ are bus to monitored lines and bus to outaged lines incidence matrices with the size of ( N B × N M ) and ( N B × N O ) , respectively. N M and N O the number of monitored lines and outaged lines in a contingency. N B are the number of all buses.
Figure 1 shows the flowchart of this iterative contingency filtering process. The detailed mathematical equations for each step are presented below:
Step 1: Solve the power dispatch problem without contingency constraints and obtain the power flow solution. The classical DC power flow model is used in this paper.
Step 2: Calculate and analyze the post-contingency power flow of all N-1 contingency cases to generate a set of N-1 contingency violation lines set and a set of the defined N-1-1 contingency cases. For example, for a given N-1 contingency case of the outage of line i , the post-contingency power flow f j , i c of any monitored line j calculated by using LODF can be represented as follows [32]:
f j , i c = f j o + LODF j , i × f i 0
where f i 0 and f j 0 are the active power flow on transmission lines i and j before contingency. LODF j , i is the LODF between lines j and i when line i is on an N-1 contingency.
After getting the post-contingency power flow, if | f j , i c | is bigger than the emergency operating limits f j emax , line j will be marked as a violation line of line i and saved to the N-1 contingency violation lines set S vio , i , which will be used to create N-1 security constraints later. If | f j , i c | is smaller than emergency operating limits f j emax but bigger than the steady-state operating limits f j max , line j is marked as a potential N-1-1 contingency line following the outage of line i . Lines i and j are saved to the disruptive N-1-1 contingency line set S N - 1 - 1 , i . The above N-1 post-contingency power flow analysis is conducted for all N-1 contingency cases in the set of S N - 1 . Note that N-1 and N-1-1 contingency cases that will split the system are not considered.
Step 3: Calculate the post-contingency power flow of all N-1-1 contingency cases selected in Step 2. For example, for an N-1-1 contingency case of lines i and j , the post-contingency power flow f k , i j c of any monitored line calculated by using LODF can be directly represented as follows:
f k , i j c = f k 0 + LODF k , i × f i 0 + LODF k , i j , j × f j , i c + LODF k , i j , i × f i , i c
where LODF k , i j , i and LODF k , i j , j are the LODF between line k and lines i and j when lines i and j are on outage together or in sequence.
Equation (38) can be further simplified. By substituting (37) into (38) and given that f i , i c is equal to zero after the outage of line i , there is:
f k , i j c = f k 0 + LODF k , i j , j × f j 0 + ( LODF k , i + LODF k , i j , i × LODF j , i ) × f i 0
According to the LODF calculation equation shown in (34)
LODF k , i = PTDF k , i 1 PTDF i , i
LODF j , i = PTDF j , i 1 PTDF i , i
LODF k , i j , i = PTDF k , i × ( 1 PTDF j , j ) + PTDF k , j × PTDF j , i ( 1 PTDF i , i ) ( 1 PTDF j , j ) PTDF j , i × PTDF i , j
LODF k , i j , j = PTDF k , i × PTDF i , j + PTDF k , j × ( 1 PTDF i , i ) ( 1 PTDF i , i ) ( 1 PTDF j , j ) PTDF j , i × PTDF i , j
It can be further proved that:
LODF k , i + LODF k , i j , j × LODF j , i = LODF k , i j , i
By substituting (44) into (39), and a more simplified equation could be obtained to calculate the post-contingency power flow of an N-1-1 contingency case as below, which is the same as an N-2 contingency case [32]:
f k , i j c = f k 0 + LODF k , i j , j × f j 0 + LODF k , i j , i × f i 0
After getting f k , i j c , if | f k , i j c | is bigger than f k emax , line k will be marked as a violation line of N-1-1 contingency of lines i and j , and it is added into the violation lines set S vio , i j . The above N-1-1 post-contingency power flow analysis for all N-1-1 contingencies cases in the set of S N - 1 - 1 , i is conducted.
Step 4: Check if both N-1 and N-1-1 violation line sets are empty. If they are empty, stop. Otherwise, continue.
Step 5: Create N-1 and N-1-1 security constraints of the active contingency scenarios for the optimization model. With the obtained violation line sets S vio , i and S vio , i j , the N-1 and N-1-1 security constraints for each contingency case in the above two sets can be formulated as follows:
f j , i s = f j + LODF j , i × f i     j S vio , i , i S N - 1
f j emax f j , i s f j emax     j S vio , i , i S N - 1
f k , i j s = f k + LODF k , i j , j × f j + LODF k , i j , i × f i     k S vio , i j , j S N - 1 - 1 , i , i S N - 1
f k emax f k , i j s f k emax     k S vio , i j , j S N - 1 - 1 , i , i S N - 1
where f l , f j , i s , and f k , i j s are optimization variables of the optimal power dispatch model which represent the power flow on a line in normal operating state, N-1, and N-1-1 contingency cases, respectively.
Step 6: Add the above constraints (46)–(49) to the proposed augmented power dispatch problem and re-optimize and solve the problem until all security checks are satisfied. If security constraints could not be ensured, a load shedding variable is needed in the model to guarantee a feasible solution [34].
For the final power dispatch solution derived through the contingency assessment process described above, any N-1 contingency on all affected transmission lines would not result in a violation. In addition, there will be no disruptive N-1-1 contingency cases, and if there are, these contingency cases will be inactive. This means that their outages will not cause other lines to be overloaded.
Overall, with the consideration of the proposed optimization objective and security constraints, a power dispatch solution with better power flow distribution can be obtained as compared with the classical economic dispatch optimization model. The possibility of cascading failures is reduced, and the power system resilience is enhanced.

3. Case Studies

To evaluate the effectiveness of the proposed model and solution techniques, simulation studies on the modified IEEE 30-bus, IEEE 118-bus, and Polish 2383-bus systems were carried out.

3.1. IEEE 30-Bus System

Three cases were studied to test the effectiveness of the power dispatch solution obtained from the proposed model in enhancing power system resilience to extreme weather events.
Case 1 studies generation costs and power flow solutions for different power dispatch models. In this case, three different power dispatch models without considering the security constraints of contingency cases are studied. Model 1 is a classical single-objective economic dispatch model in which only the generation cost is considered as the optimization objective. Model 2 is the resilience-constrained economic dispatch model proposed in our previous work [24], where objectives related to power flows are considered, but the adjustment of transmission line reactance by controllable series compensation devices is not considered. Model 3 is the proposed augmented power dispatch model considering the control of series compensation devices. For the sake of brevity, in the following discussion Models 1 through 3 are denoted as ED, RCED, and APD, respectively.
Case 2 studies the effectiveness of the proposed N-1-1 security criterion and the contingency assessment process. In this case, the pre-contingency and post-contingency power flows of different models are analyzed and compared.
Case 3 studies the reliability and resilience performances of the above models by using a modified cascading collapse assessment with the consideration of hidden failure [20]. The probability distribution curve of the blackout size will be used to show the resilience performance of a power system in this case.
Detailed data for the test system are in Appendix A. S AL is formed by assuming lines 10, 16, 22, 29, 30, 33, 35, 37, and 38, which are affected by an extreme weather event. Assuming that all lines are equipped with controllable series compensation devices, the maximum adjustment value of line susceptance is set to 90% of the initial value. In all cases, for power dispatch models considering line flow adjustment objectives, the coefficients of the different objective terms are set to α = 1000 , β = 1000 , and γ = 1 , which are determined according to the order of magnitude of the different objective terms. In Cases 2 and 3, the emergency limit is set to 120% of the steady-state operation limit. All cases were tested in MATLAB 2016a using the Gurobi solver on a personal computer with a 3.20 GHz i5 processor and 8 GB RAM.

3.1.1. Case 1

Since the first and second objectives are introduced to minimize the active power flow in the transmission system, taking into account the series compensation devices, first the active power distribution in the transmission system was studied to investigate the effectiveness of the proposed objectives and the joint control action. Figure 2a,b show the loading rate distribution of affected lines and all lines, respectively. As shown in Figure 2, the loading rate of the affected lines and all lines of the RCED are similar to those of the ED. This indicates that, in some cases, the power flow distribution cannot be improved by controlling outputs of generators alone, even when the objectives related to power flow adjustment are considered. One reason for this result is that as load demand increases, the number of units committed to generation increases and the dispatch capacity to adjust power flows decreases. With the adjustment of transmission line impedance by controllable series compensation devices, the proposed augmented power flow dispatch model APD has good performance in improving the power flow distribution on the line. The loading rates of the affected lines in the APD are much lower than the others. The loading rates of all lines in the proposed model are also smoother than those of the other two models. That is, the inhomogeneity of power flow on the branches is small and the power flow distribution in the transmission system is uniform. In addition, the relatively small number of heavy and light load lines in the APD means that not only can heavily loading lines be avoided, ensuring safe system operation during extreme weather events but also light loading lines can be reduced, improving the overall utilization of transmission capacity.
To quantify the power flow distribution in the transmission system, we calculated seven indices, including the average loading rate (Ave. r l ) and (Ave. r a l ), the variance of the loading rate (Var. r l ), the maximum loading rate (Max. r l ), and the number of heavily loading lines (Num. r l = 1 , Num. r l > 0.8 and Num. r l > 0.6 ). The obtained results are shown in Table 1 below.
As shown in Table 1, the classical economic dispatch model ED has the worst results. Four lines operate at their operating limits (Num. r l = 1 ). These lines are more likely to experience overload outages or hidden failure outages when power flows fluctuate. Outages on these fully loading lines will also result in large-scale power flow transfers, which may result in cascading failure outages. As for the proposed model APD, the power flow distribution is significantly improved compared to the other two models, thanks to the adjustment of the transmission line reactance by the series compensation device. The average loading rate and loading rate variance for all lines, as well as the average loading rate for the affected lines, were significantly reduced (by almost half compared to the other models). As for the heavily loaded lines, only one transmission line operates at its capacity in the proposed model, while in the ED and RCED models there are four and three, respectively. As the loading rate decreases, the probability of relay malfunction is also greatly reduced, which will make the power system more secure. Compared with RCED, APD has a lower generation cost (that is, the proposed model can achieve a better power flow distribution with less cost). In summary, the simulation results obtained in this case validate the effectiveness of the proposed augmented power dispatch model and the solution.

3.1.2. Case 2

Case 2 investigates the effectiveness of the proposed N-1-1 security criteria and the contingency assessment process. In this case, security constraints for the defined N-1 and N-1-1 contingency cases are considered in the proposed augmented power dispatch model. In the following discussion, the model with security constraints is denoted as SCAPD. In this case, the pre-contingency and post-contingency power flow distributions of the different models are compared. The analysis of Case 1 shows that the APD has a good power flow distribution state when the power system is in steady-state operation. However, since it does not take into account the security constraints of any contingency case, the system may experience load shedding when transmission line failure occurs. To study this problem, first N-1 and N-1-1 contingency analyses are performed on the generation scheduling schemes obtained from different models to see how many power flows exceed the operating limits. The analysis results are shown in Table 2. S1 denotes the number of scenarios in which the post-contingency power flow exceeds the emergency limit for N-1 contingency cases. S2 denotes the number of scenarios in which the post-contingency power flow exceeds the steady-state operating limit but is less than the emergency limit for N-1 contingency cases. These scenarios are used to generate the proposed N-1-1 contingency cases. S3 denotes the number of scenarios in which the post-contingency power flow exceeds the emergency limit for the proposed N-1-1 contingency cases.
As can be seen from Table 2, scenario S1 for APD is much smaller than that for ED. However, in the APD model, there are still five lines with post-contingency power flows exceeding the emergency limit under the N-1 contingency. For example, when line 10 is out of service, the post-contingency loading rates of lines 40 and 41 are 1.282 and 1.296, respectively, which means that lines 40 and 41 will be out of service immediately after Line 10 is out of service. N-1 security constraints should be added to lines 40 and 41 to ensure the security of system operation. Scenarios S2 for APD shows that there are 21 scenarios in which the post-contingency power flow exceeds the steady-state operating limit but is less than the emergency limit of the transmission line. For example, when line 15 fails, line 10 has a post-contingency loading rate of 1.016. When there is a disturbance in the system, line 10 may experience an overload outage. If line 10 fails, two lines will be overloaded and this will develop into a cascading failure outage. Therefore, the proposed N-1-1 security constraint for preventing cascading blackouts should be considered in the power dispatching problem, especially in extreme weather conditions.
The proposed N-1-1 security constraints are added in the APD by using the iterative contingency assessment process described in Section 2.3. The model converged after two iterations. This validates the computational efficiency of the LODF-based contingency assessment process. the power flow distribution indices of SCAPD are shown in row three of Table 2. The SCAPD model has no overloaded lines in all N-1 and N-1-1 contingency analyses due to the security constraints of the N-1 and N-1-1 contingency cases. However, if only the traditional N-1 security constraint is considered in the APD model, there are two lines with post-contingency line flows exceeding the operating limits but under the emergency limits in all N-1 contingency analyses, as shown in the fourth row of Table 2. For example, when line 10 fails, line 41 has a post-contingency loading ratio of 1.2. Lines 10 and 41 will become one scenario in S2. If this N-1-1 contingency case occurs, it will result in three lines being overloaded, which could develop into a cascading outage. Therefore, by considering the security constraints of the defined N-1-1 contingency in the power dispatch model, to a certain extent, the risk of cascading failures in the power system can be reduced, thus improving the operational recovery of the system.

3.1.3. Case 3

In this case, the reliability and resilience performances of different power dispatch models in an extreme weather event are studied by using a modified cascading collapse assessment considering hidden failure. The initial first step of the cascading simulation procedure presented in [35] was changed to consider an initial trigger event caused by an extreme weather event, which was generated by Monte Carlo sampling. It is assumed that each line has a probability of outage in the initial state, but the affected lines have a higher probability of failure. In this case, for lines affected by the extreme weather event, each line has a probability of failure between 0 and 0.05, which is randomly generated in each simulation. The rest of the lines have a smaller random failure probability between 0 and 0.005. To avoid initial triggering events leading to system disconnections, a maximum of three lines were allowed to fail in the initial event. We performed one thousand simulations and calculated the percentage of load shedding in each simulation to represent the cascading blackout size. The probability distribution curves of the blackout sizes for different models are given in Figure 3, and the detailed results are shown in Table 3.
As shown in Figure 3, the proposed augmented power dispatch model with controllable series compensation devices (APD) has a better blackout size distribution than the classical ED and RCED; this result is consistent with the discussion in Case 1 and Case 2. SCAPD is the most resilient strategy among all models due to the consideration of N-1 and the proposed N-1-1 security constraints. Table 3 shows the detailed results of all cascade failure simulations. It can be seen that the average outage lines caused by extreme weather events at the initial trigger (the fifth column) are similar for all four models. This indicates that the different blackout sizes between models are not caused by the initial outage, but by subsequent cascading outages. The sixth column of Table 3 shows the average overload failure lines per simulation for the different models. The ED and RCED models have relatively more overload failures due to weaker power flow distribution conditions. The SCAPD produces the least overload failures. The average load shedding for SCAPD is only 2.57%, which is about 40% less compared to the results for ED and RCED. The probability of outage size exceeding 15% in the SCAPD model is very rare, at 0.005. From these simulation results, it can be seen that the proposed model gives a more secure power dispatch solution that is more resilient in extreme events.
Overall, power flow distribution conditions have a significant impact on the operational security of the system. The fewer highly loaded lines in a system, the lower the probability of overload cascading failures, thus making the system more resilient to extreme events. Considering the security constraints of the disruptive N-1-1 contingency cases in the prevention strategy contributes to the resilience of the power system under extreme weather conditions.

3.2. IEEE 118-Bus and Polish 2383-Bus Systems

Detailed data for the IEEE 118-bus test system are in Appendix A. For the model considering power flow-related objectives, the coefficients of the different objective terms are set as, α = 1 , 000 , 000 , β = 1 , 000 , 000 , and γ = 1 , depending on the order of magnitude of the different objective terms. It is assumed that lines 1 to 90 are affected by extreme weather events. Assuming that all lines are equipped with controllable compensation devices, and the maximum adjusted value of line susceptibility is set to 90% of the initial value. The limit of the contingency is set to 120% of the steady-state operating limit. Data for the Polish 2383-bus test system are MATPOWER “case2383wp”. It is assumed that lines 300 to 500 are affected by extreme weather events. Lines 1 to 100 are equipped with controllable compensation devices and the maximum adjusted value of line susceptibility is set to 20% of the initial value. The coefficients of the different objective terms are set as, α = 1 , 000 , 000 , β = 1 , 000 , 000 and γ = 1 . All cases were tested in MATLAB 2016a, using the Gurobi solver, on a PC with a 3.20 GHz i5 processor and 8 GB of RAM. Table 4 shows the results of the iterative contingency filtering process of the SCAPD model tested in the IEEE 118-bus system. Table 5 shows the results of the iterative contingency filtering process of the SCAPD model tested in the Polish 2383-bus system. Total iterations and computation time for solving SCAPD in different test systems are demonstrated in Table 6.
As shown in Table 4 and Table 5, as the number of iterations increases, the number of N-1 and N-1-1 violations decreases. Case studies on the IEEE 118-bus system converge after three iterations. Case studies on the Polish 2383-bus system converge after four iterations. Table 6 summarizes the total iteration and computation time for the three test systems used in the paper. It can be seen that with the increase of the system size, the computation time increases. However, the largest computation time is still within the acceptable range. These results illustrate the computational efficiency of the proposed two-stage and iterative contingency filtering solution techniques.
As for the effectiveness of the proposed power dispatch model, we will discuss the results of case studies for different power dispatch models in the IEEE 118-bus system in detail. Figure 4a,b show the loading rates of the lines affected by extreme weather events and all lines in the four models mentioned above. Figure 5 shows the probability distribution curves of blackout sizes obtained from 1000 cascading blackout simulations. The detailed generation cost, power flow distribution index, and cascading outage simulation results are shown in Table 7.
As shown in Figure 4a,b, in the study conducted for the IEEE 118-bus system, as with the IEEE 30-bus system, the loading rates of the affected lines in APD and SCAPD are lower than in the other models because of the objective one. The loading rates of all lines in APD and SCAPD are also smoother than the other models because of objective two. As can be seen in Table 7, the average loading rate in SCAPD is very small; no lines operate on the operating limit. The number of heavily loaded lines is only one-third of that of ED; this indicates that the joint control of generators and controllable compensation devices in collaboration with objectives one and two can reduce the overall loading rate of the transmission system.
For security and resilience performance, it can be seen from Table 7 that the classical ED model has 44 violations in the N-1 contingency analysis, while the proposed APD model only has 10. The classical ED model has 535 N-1-1 contingency cases, and these contingency cases would result in 719 violations. However, APD has only 57 potential N-1-1 contingencies, which is much smaller than ED. SCAPD has no violations in both N-1 and N-1-1 contingency analysis. The power system would operate more safely with a smaller number of violations. Figure 5 shows the probability distribution curves of outage sizes obtained from the 1000 cascading outage simulations.
As shown in Figure 5, the probability distribution curves of APD and SCAPD are much better than that of ED because of the lower loading rate. SCAPD gives the most resilient strategy among all models due to the consideration of N-1 and the proposed N-1-1 security constraints. The maximum load shedding percentage is only 3.86%, which is about 60% lower compared to ED. The average number of overload and hidden failure lines in each simulation is only 1.57 compared to 9.35 in ED. This indicates that security constraints of the N-1 and N-1-1 contingency cases effectively prevent the propagation of initial outages caused by extreme weather events. In summary, the proposed objective function and the adjustment of the transmission line reactance by a controllable series compensation device help to improve the power flow distribution. By considering the proposed N-1-1 security criteria and contingency assessment process, the power system operation scheme is more resilient to extreme weather events. The proposed model and solution methodology can be effectively applied to large-scale power systems.

4. Discussions

This paper proposes an augmented power dispatch model that takes into account controllable series compensation devices and disruptive N-1-1 contingency security constraints to provide a resilient operating scheme for power system operators to mitigate damage from extreme events.
Compared with the conventional economic dispatch model that only considers generation cost as the optimization objective, the proposed augmented power dispatch model aims to optimize the power flow distribution in the transmission system to reduce the loading rate of the transmission lines. Compared with the resilience-constrained economic dispatch proposed by [24], the adjustment of transmission line reactance by controllable series compensation devices is considered in this paper. Figure 2 and Figure 4 demonstrate the loading rate distribution of the above models in the two test systems. Table 1 and Table 7 show the comparison of operation cost and power flow distribution indices of these models. It can be seen that the average loading rate of the proposed model is smaller than the other two models and the number of heavy loading lines is reduced in the proposed model. The power dispatch model considering the defined N-1-1 and N-1 security constraints is compared with the model considering only the conventional N-1 security constraints. Table 2 shows the comparison of the pre-contingency and post-contingency power flow distributions for these models. The results indicate that the proposed N-1-1 security constraints have an impact on mitigating the further propagation of N-1 contingency cases and reducing the risk of overload cascade failures. However, it can be imagined that the considered security constraints cannot ensure the safe operation of the system for high order contingency cases. This is one of the reasons that the objective for power flow improvement is considered in the optimization model, which makes the proposed model much different from the existing related studies. Figure 3 and Figure 5 show the probability distribution curves of the blackout sizes for different models, which reflect the reliability and resilience performances of the models. Table 3 shows simulation results of cascading failure blackouts for different models tested on the IEEE 30-bus system. From the simulation results, it can be seen that the proposed model gives a more secure power dispatch solution that is more resilient in extreme events. Table 4 shows the results of the iterative contingency filtering process of the SCAPD model tested in the IEEE 118-bus system. Table 5 shows the results of the iterative contingency filtering process of the SCAPD model tested in the Polish 2383-bus system. Table 6 shows the total iterations and computation time for solving SCAPD in different test systems. These results illustrate the computational efficiency of the proposed two-stage and iterative contingency filtering solution techniques.
In future work, AC power flow models are considered to obtain more practical solutions for power system operators. Renewable energy and energy storage systems will also be considered to increase the application of the model.

5. Conclusions

In this paper, an augmented power dispatch model is proposed to provide a resilient operating scheme for power system operators to mitigate damage from extreme events. The major findings are summarized below:
  • By controlling series compensation devices to adjust the impedance of transmission lines, the proposed augmented power dispatch model can avoid some lines from taking excessive loads, especially those lines affected by extreme weather events.
  • The reduction of power flow on lines can reduce the possibility of overload cascade faults, thus reducing load shedding under extreme weather events and improving the resilience of the power system.
  • The proposed N-1-1 security criteria have an impact on mitigating the further propagation of N-1 contingency cases and reducing the risk of overload cascade failures.
  • The proposed iterative contingency assessment process enables us to solve the security-constrained power dispatch problem iteratively, reducing the problem size and computation time.

Author Contributions

L.L.L., C.S.L. and G.Y. guided the framework of the paper and provided professional guidance. L.H. did the simulation and wrote the paper. Z.H., Z.Z., N.T. and X.W. provided academic assistance and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Education Department of Guangdong Province: New and Integrated Energy System Theory and Technology Research Group [Project Number 2016KCXTD022]; National Natural Science Foundation of China (51907031); Brunel University London BRIEF Funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the author: Liping Huang (e-mail address: [email protected]).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The system data of the modified IEEE 30-bus and 118-bus test systems used in the paper is shown in Table A1, Table A2, Table A3, Table A4, Table A5, Table A6.
Table A1. Generator data of the modified IEEE 30-bus system.
Table A1. Generator data of the modified IEEE 30-bus system.
UBus
No.
Unit Cost CoefficientsPmax
(MW)
Pmin
(MW)
a
(MBtu)
b
(MBtu/MW)
c
(MBtu/MW2)
G1102.0000 0.0200 1200
G2201.7500 0.0175 1200
G32201.0000 0.0625 750
G42703.2500 0.0083 82.50
G52303.0000 0.0250 450
G61303.0000 0.0250 600
Table A2. Bus data of the modified IEEE 30-bus system.
Table A2. Bus data of the modified IEEE 30-bus system.
Bus No.d(MW)Bus No.d(MW)Bus No.d(MW)Bus No.d(MW)
10901711.7250
228.21107.54184.16264.55
33.121101912.35270
49.881214.56202.86280
501302122.75293.12
60148.062203013.78
729.641510.66234.16
839164.552411.31
Table A3. Line data of the modified IEEE 30-bus system.
Table A3. Line data of the modified IEEE 30-bus system.
Line
No.
From
Bus
To
Bus
X (pu)Flow
Limit
(MW)
Line
No.
From
Bus
To
Bus
X (pu)Flow
Limit
(MW)
1120.06123.52215180.2215.2
2130.19123.52318190.1315.2
3240.1761.752419200.0730.4
4340.04123.52510200.2130.4
5250.2123.52610170.0830.4
6260.1861.752710210.0730.4
7460.0485.52810220.1530.4
8570.1266.52921220.0230.4
9670.08123.53015230.215.2
10680.0430.43122240.1815.2
11690.2161.753223240.2715.2
126100.5630.43324250.3315.2
139110.2161.753425260.3815.2
149100.1161.753525270.2115.2
154120.2661.753628270.461.75
1612130.1461.753727290.4215.2
1712140.2630.43827300.615.2
1812150.1330.43929300.4515.2
1912160.230.4408280.230.4
2014150.215.2416280.0630.4
2116170.1915.2
Table A4. Generator data of the modified IEEE 118-bus system.
Table A4. Generator data of the modified IEEE 118-bus system.
UnitBus No.Unit Cost CoefficientsPmax (MW)Pmin (MW)UBus No.Unit Cost CoefficientsPmax (MW)Pmin (MW)
a
(MBtu)
b
(MBtu/MW)
c
(MBtu/MW2)
a
(MBtu)
b
(MBtu/MW)
c
(MBtu/MW2)
110400.010300028650200.02614730
240400.010300029660200.02614760
360400.010300030690200.0192415.60
480400.010300031700400.0103000
5100200.0221650032720400.0103000
6120200.118555033730400.0103000
7150400.010300034740400.0103000
8180400.010300035760400.0103000
9190400.010300036770400.0103000
10240400.010300037800200.02117310
11250200.045960038850400.0103000
12260200.0321242039870202.5003120
13270400.010300040890200.01621210
14310201.429321041900400.0103000
15320400.010300042910400.0103000
16340400.010300043920400.0103000
17360400.010300044990400.0103000
18400400.0103000451000200.04010560
19420400.0103000461030200.2504200
20460200.5263570471040400.0103000
21490200.0499120481050400.0103000
22540200.2084440491070400.0103000
23550400.0103000501100400.0103000
24560400.0103000511110200.2784080
25590200.0657650521120400.0103000
26610200.0637800531130400.0103000
27620400.0103000541160400.0103000
Table A5. Bus data of the modified IEEE 118-bus system.
Table A5. Bus data of the modified IEEE 118-bus system.
Bus
No.
D
(MW)
Bus
No.
D
(MW)
Bus
No.
D
(MW)
Bus
No.
D
(MW)
1153311296109130
260321776223192195
311733696309336
4117341776409490
50359965095126
615636936611796114
75737067849745
88438068098102
90398169099126
1004019870198100111
112104111171010166
1214142288723610215
131024354731810369
1442444874204104114
15270451597514110593
1675468476204106129
17334710277183107150
181804860782131086
19135492617911710924
2054505180390110117
214251518101110
2230525482162112204
23215369836011318
243954339843311424
25055189857211566
260562528663116552
27213573687011760
285158368814411899
297259831890
3006023490489
Table A6. Line data of the modified IEEE 118-bus system.
Table A6. Line data of the modified IEEE 118-bus system.
Line
No.
From
Bus
To
Bus
X (pu)Flow
Limit
(MW)
Line
No.
From
Bus
To
Bus
X (pu)Flow
Limit
(MW)
1120.09991769455560.0151176
2130.04241769555590.2158176
32120.06161769656570.0966176
4350.1081769756580.0966176
53120.161769856590.251176
6450.007983529956590.239176
74110.068817610059600.145176
8560.05417610159610.15176
95110.068217610260610.0135352
10670.020817610360620.0561176
117120.03417610461620.0376176
12890.030588010562660.218176
13850.026770410662670.117176
148300.050417610763590.0386352
159100.032288010863640.02352
1611120.019617610964610.0268176
1711130.073117611064650.0302352
1812150.070717611165660.037176
1912170.083417611265680.016176
20121170.1417611366670.1015176
2113150.244417611468690.037352
2214150.19517611568810.0202176
2315170.0437352116681160.00405352
2415190.039417611769700.127352
2515330.124417611869750.122352
2616170.180117611969770.101176
2717190.050517612070710.0355176
2817310.156317612170740.1323176
29171130.030117612270750.141176
3018190.049317612371720.18176
3119200.11717612471730.0454176
3219340.24717612574750.0406176
3320210.084917612675770.1999176
3421220.097176127751180.0481176
3522230.15917612876770.148176
3623240.0492176129761180.0544176
3723250.0835213077780.0124176
3823320.115317613177800.0485352
3924700.411517613277800.105176
4024720.19617613377820.0853176
4125270.16335213478790.0244176
4226250.038217613579800.0704176
4326300.08652813680960.182176
4427280.085517613780970.0934176
4527320.075517613880980.108176
46271150.074117613980990.206176
4728310.094317614081800.037176
4829310.033117614182830.03665176
4930170.038852814282960.053176
5030380.05417614383840.132176
5131320.098517614483850.148176
52321130.20317614584850.0641176
53321140.061217614685860.123176
5433370.14217614785880.102176
5534360.026817614885890.173176
5634370.009435214986870.2074176
5734430.168117615088890.0712352
5835360.010217615189900.032528
5935370.049717615289910.032176
6037390.10617615389920.0505176
6137400.16817615490910.0505528
6238370.037552815591920.1272176
6338650.098635215692930.032176
6439400.060517615792940.158176
6540410.0487176158921000.295176
6640420.183176159921020.0559176
6741420.13517616093940.0732176
6842490.32317616194950.0434176
6942490.32317616294960.0869176
7043440.2454176163941000.058176
7144450.090117616495960.0547176
7245460.135617616596970.0885176
7345490.186176166981000.179176
7446470.127176167991000.0813176
7546480.1891761681001010.1262176
7647490.06251761691001030.0525352
7747690.27781761701001040.204176
7848490.05051761711001060.229176
7949500.07521761721011020.112176
8049510.1371761731031040.1584176
8149540.2891761741031050.1625176
8249540.2911761751031100.1813176
8349660.09193521761041050.0378176
8449660.09193521771051060.0547176
8549690.3241761781051070.183176
8650570.1341761791051080.0703176
8751520.05881761801061070.183176
8851580.07191761811081090.0288176
8952530.16351761821091100.0762176
9053540.1221761831101110.0755176
9154550.07071761841101120.064176
9254560.009551761851141150.0104176
9354590.2293176

References

  1. Jufri, F.; Widiputra, V.; Jung, J. State-of-the-art review on power grid resilience to extreme weather events: Definitions, frameworks, quantitative assessment methodologies, and enhancement strategies. Appl. Energy 2019, 239, 1049–1065. [Google Scholar] [CrossRef]
  2. Vaccaro, A.; Pisica, I.; Lai, L.L.; Zobaa, A.F. A review of enabling methodologies for information processing in smart grids. Int. J. Electr. Power Energy Syst. 2019, 107, 516–522. [Google Scholar] [CrossRef]
  3. Lai, C.S.; Jia, Y.; Dong, Z.; Wang, D.; Tao, Y.; Lai, Q.H.; Wong, R.T.K.; Zobaa, A.F.; Wu, R.; Lai, L.L. A review of technical standards for smart cities. Clean Technol. 2020, 2, 19. [Google Scholar] [CrossRef]
  4. Bie, Z.; Lin, Y.; Li, G.; Li, F. Battling the extreme: A study on the power system resilience. Proc. IEEE 2017, 105, 1253–1266. [Google Scholar] [CrossRef]
  5. Xu, Z.; Lai, L.L.; Wong K., P.; Pinson, P.; Li, F.X. Guest editorial-special section on emerging informatics for risk hedging and decision making in smart grids. IEEE Trans. Ind. Inform. 2017, 13, 2507–2510. [Google Scholar] [CrossRef]
  6. Panteli, M.; Mancarella, P. The grid: Stronger, bigger, and smarter? Presenting a conceptual framework of power system resilience. IEEE Power Energy Mag. 2015, 13, 58–66. [Google Scholar] [CrossRef]
  7. Zhang, G.; Zhang, F.; Zhang, X.; Wu, Q.; Ke, M. A multi-disaster-scenario distributionally robust planning model for enhancing the resilience of distribution systems. Int. J. Electr. Power Energy Syst. 2020, 122, 106161. [Google Scholar] [CrossRef]
  8. Pepiciello, A.; Vaccaro, A.; Lai, L.L. An interval mathematic-based methodology for reliable resilience analysis of power systems in the presence of data uncertainties. Energies 2020, 13, 6632. [Google Scholar] [CrossRef]
  9. Tari, A.N.; Sepasian, M.S.; Kenari, M.T. Resilience assessment and improvement of distribution networks against extreme weather events. Int. J. Electr. Power Energy Syst. 2021, 125, 106414. [Google Scholar] [CrossRef]
  10. Panteli, M.; Trakas, D.N.; Mancarella, P.; Hatziargyriou, N.D. Power systems resilience assessment: Hardening and smart operational enhancement strategies. Proc. IEEE 2017, 105, 1202–1213. [Google Scholar] [CrossRef] [Green Version]
  11. Mahzarnia, M.; Moghaddam, M.P.; Baboli, P.T.; Siano, P. A review of the measures to enhance power systems resilience. IEEE Syst. J. 2020, 14, 4059–4070. [Google Scholar] [CrossRef]
  12. Arab, A.; Khodaei, A.; Khator, S.K.; Ding, K.; Emesih, V.A.; Han, Z. Stochastic pre-hurricane restoration planning for electric power systems infrastructure. IEEE Trans. Smart Grid 2015, 6, 1046–1054. [Google Scholar] [CrossRef]
  13. Eskandarpour, R.; Edwards, G.; Khodaei, A. Resilience-constrained unit commitment considering the impact of microgrids. In Proceedings of the 2016 North American Power Symposium (NAPS), Denver, CO, USA, 18–20 September 2016; pp. 1–5. [Google Scholar]
  14. Xie, C.; Wang, D.; Lai, C.S.; Wu, R.; Wu, X.; Lai, L.L. Optimal sizing of battery energy storage system in smart microgrid considering virtual energy storage system and high photovoltaic penetration. J. Clean. Prod. 2021, 281, 125308. [Google Scholar] [CrossRef]
  15. Sahraei-Ardakani, M.; Ou, G. Day-ahead preventive scheduling of power systems during natural hazards via stochastic optimization. In Proceedings of the 2017 IEEE Power & Energy Society General Meeting, Chicago, IL, USA, 16–20 July 2017; p. 1. [Google Scholar]
  16. Khaloie, H.; Vallee, F.; Lai, C.S.; Toubeau, J.F.; Hatziargyriou, N.D. Day-ahead and intraday dispatch of an integrated biomass-concentrated solar system: A multi-objective risk-controlling approach. IEEE Trans. Power Syst. 2021. [Google Scholar] [CrossRef]
  17. Wang, C.; Hou, Y.; Qiu, F.; Lei, S.; Liu, K. Resilience enhancement with sequentially proactive operation strategies. IEEE Trans. Power Syst. 2017, 32, 2847–2857. [Google Scholar] [CrossRef] [Green Version]
  18. Wang, Y.; Huang, L.; Shahidehpour, M.; Lai, L.L.; Yuan, H.; Xu, F.Y. Resilience-constrained hourly unit commitment in electricity grids. IEEE Trans. Power Syst. 2018, 33, 5604–5614. [Google Scholar] [CrossRef]
  19. Trakas, D.N.; Hatziargyriou, N.D. Resilience constrained day-ahead unit commitment under extreme weather events. IEEE Trans. Power Syst. 2020, 35, 1242–1253. [Google Scholar] [CrossRef]
  20. Jie, C.; Thorp, J.S.; Dobson, I. Cascading dynamics and mitigation assessment in power system disturbances via a hidden failure model. Int. J. Electr. Power Energy Syst. 2005, 27, 318–326. [Google Scholar]
  21. Tejada-Arango, D.A.; Sánchez-Martın, P.; Ramos, A. Security constrained unit commitment using line outage distribution factors. IEEE Trans. Power Syst. 2018, 33, 329–337. [Google Scholar] [CrossRef]
  22. Bao, Z.J.; Cao, Y.J.; Wang, G.Z.; Ding, L.J. Analysis of cascading failure in electric grid based on power flow entropy. Phys. Lett. A 2009, 373, 3032–3040. [Google Scholar] [CrossRef]
  23. Zhang, S.; Liu, C.C.; Gu, X.; Wang, T. Optimal transmission line switching incorporating dynamic line ratings. In Proceedings of the 2017 IEEE PES Innovative Smart Grid Technologies Conference Europe (ISGT-Europe), Torino, Italy, 26–29 September 2017; pp. 1–5. [Google Scholar]
  24. Wang, Y.; Huang, L.; Shahidehpour, M.; Lai, L.L.; Zhou, Y. Impact of cascading and common-cause outages on resilience-constrained optimal economic operation of power systems. IEEE Trans. Smart Grid 2020, 11, 590–601. [Google Scholar] [CrossRef]
  25. Ding, T.; Bo, R.; Li, F.; Sun, H. Optimal power flow with the consideration of flexible transmission line impedance. IEEE Trans. Power Syst. 2016, 31, 1655–1656. [Google Scholar] [CrossRef]
  26. Singh, A.K.; Pal, B.C. Dynamic Estimation and Control of Power Systems; Academic Press: Cambridge, MA, USA, 2018; p. 5. [Google Scholar]
  27. Sahraei-Ardakani, M.; Hedman, K.W. Day-ahead corrective adjustment of FACTS reactance: A linear programming approach. IEEE Trans. Power Syst. 2016, 31, 2867–2875. [Google Scholar] [CrossRef]
  28. Sahraei-Ardakani, M.; Hedman, K.W. Computationally efficient adjustment of FACTS set points in DC optimal power flow with shift factor structure. IEEE Trans. Power Syst. 2017, 32, 1733–1740. [Google Scholar] [CrossRef]
  29. Jia, Y.; Xu, Z.; Lai, L.L.; Wong, K.P. Risk-based power system security analysis considering cascading outages. IEEE Trans. Ind. Informat. 2016, 12, 872–882. [Google Scholar] [CrossRef]
  30. Lai, L.L.; Ma, J.T.; Yokoyama, R.; Zhao, M. Improved genetic algorithms for optimal power flow under both normal & contingent operation states. Int. J. Electr. Power Energy Syst. 1997, 19, 287–292. [Google Scholar]
  31. Liu, M.; Gross, G. Role of distribution factors in congestion revenue rights applications. IEEE Trans. Power Syst. 2004, 19, 802–810. [Google Scholar] [CrossRef]
  32. Guo, J.; Fu, Y.; Li, Z.; Shahidehpour, M. Direct calculation of line outage distribution factors. IEEE Trans. Power Syst. 2009, 24, 1633–1634. [Google Scholar]
  33. Guler, T.; Gross, G.; Liu, M. Generalized line outage distribution factors. IEEE Trans. Power Syst. 2007, 22, 879–881. [Google Scholar] [CrossRef]
  34. Fu, Y.; Shahidehpour, M.; Li, Z. AC contingency dispatch based on security-constrained unit commitment. IEEE Trans. Power Syst. 2006, 21, 897–908. [Google Scholar] [CrossRef]
  35. Salim, N.A.; Othman, M.M.; Musirin, I.; Serwan, M.S. Cascading collapse assessment considering hidden failure. In Proceedings of the 1st Int Conf Informatics Comput Intell ICI 2011, Bandung, Indonesia, 12–14 December 2011; pp. 318–323. [Google Scholar]
Figure 1. Flowchart of the proposed iterative contingency filtering process.
Figure 1. Flowchart of the proposed iterative contingency filtering process.
Energies 14 04756 g001
Figure 2. Loading rates for (a) lines affected by extreme weather and (b) all lines ranked from highest to lowest.
Figure 2. Loading rates for (a) lines affected by extreme weather and (b) all lines ranked from highest to lowest.
Energies 14 04756 g002
Figure 3. Probability distribution curves of blackout sizes for different models tested on the IEEE 30-bus system.
Figure 3. Probability distribution curves of blackout sizes for different models tested on the IEEE 30-bus system.
Energies 14 04756 g003
Figure 4. Loading rates for (a) lines affected by extreme weather and (b) all lines ranked from highest to lowest of the IEEE 118-bus system.
Figure 4. Loading rates for (a) lines affected by extreme weather and (b) all lines ranked from highest to lowest of the IEEE 118-bus system.
Energies 14 04756 g004
Figure 5. Probability distribution curves of blackout sizes for different models tested on the IEEE 118-bus system.
Figure 5. Probability distribution curves of blackout sizes for different models tested on the IEEE 118-bus system.
Energies 14 04756 g005
Table 1. Comparison of power flow distribution indices of different power dispatch models tested on the IEEE 30-bus system.
Table 1. Comparison of power flow distribution indices of different power dispatch models tested on the IEEE 30-bus system.
ModelCost ($)Ave. r l Ave. r l Var. r l Max. r l Num.
r l = 1
Num.
r l > 0.8
Num.
r l > 0.6
ED801.570.3520.790.0721446
RCED826.620.350.780.0681347
APD819.30.3230.3820.0331113
Table 2. Comparison of the pre-contingency and post-contingency power flow distributions for different models tested on the IEEE 30-bus system.
Table 2. Comparison of the pre-contingency and post-contingency power flow distributions for different models tested on the IEEE 30-bus system.
ModelsCost ($)Ave. r l Var. r l Max. r l S1S2S3
ED801.50.3520.0721135968
APD819.30.3230.033152135
SCAPD777.10.3130.0290.77000
SCAPD(N-1)802.60.3200.0290.93023
Table 3. Simulation results of cascading failure blackouts for different models tested on the IEEE 30-bus system.
Table 3. Simulation results of cascading failure blackouts for different models tested on the IEEE 30-bus system.
ModelsMax. LS%Probability
(LS% > 15%)
Ave.
LS%
Outages Caused
by the Event
Outages Caused
by Overload
ED19.5%0.0854.21%1.032.22
RCED20.21%0.0504.24%1.022.43
APD19.23%0.0152.98%1.051.52
SCAPD18.07%0.0052.57%1.031.05
Table 4. Results of the iterative contingency filtering process of the SCAPD model tested on the IEEE 118-bus system.
Table 4. Results of the iterative contingency filtering process of the SCAPD model tested on the IEEE 118-bus system.
IterationsGeneration CostNumber of N-1 Violation ScenariosNumber of N-1-1 Violation ScenariosComputation Time (s)
0582,259.310154.9
1575,352.53194.4
2563,046.5132.1
3572,850.9003.1
Table 5. Results of the iterative contingency filtering process of the SCAPD model tested on the Polish 2383-bus system.
Table 5. Results of the iterative contingency filtering process of the SCAPD model tested on the Polish 2383-bus system.
IterationsGeneration CostNumber of N-1 Violation ScenariosNumber of N-1-1 Violation ScenariosComputation Time (s)
01,851,036.565711174.6
12,006,177.359170.0
21,981,281.765195.9
32,018,414.201185.6
42,004,039.800171.2
Table 6. Total iterations and computation time for solving SCAPD in different test systems.
Table 6. Total iterations and computation time for solving SCAPD in different test systems.
Test SystemsTotal IterationsTotal Time (s)
IEEE 30-bus37.8
IEEE 118-bus314.5
Polish 2383-bus4897.3
Table 7. Pre-contingency and post-contingency power flow distribution indices and cascade blackout simulation results for different models tested on the IEEE 118-bus system.
Table 7. Pre-contingency and post-contingency power flow distribution indices and cascade blackout simulation results for different models tested on the IEEE 118-bus system.
ModelsCost ($)Ave. r l Var. r l Max. r l Num. r l = 1 Num. r l > 0.8 Num. r l > 0.6 S1S2S3Max.
LS%
Ave.
LS%
Outages Caused by the EventOutages
Caused by Overload
ED489,152.50.3820.0741818384453571910.93%1.23%2.239.35
RCED556,491.50.3010.05112824151111364.57%0.36%2.132.46
APD582,259.30.2570.041114131057134.12%0.30%2.112.26
SCAPD572,850.90.2510.0320.8601110803.86%0.19%2.171.57
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Huang, L.; Huang, Z.; Lai, C.S.; Yang, G.; Zhao, Z.; Tong, N.; Wu, X.; Lai, L.L. Augmented Power Dispatch for Resilient Operation through Controllable Series Compensation and N-1-1 Contingency Assessment. Energies 2021, 14, 4756. https://doi.org/10.3390/en14164756

AMA Style

Huang L, Huang Z, Lai CS, Yang G, Zhao Z, Tong N, Wu X, Lai LL. Augmented Power Dispatch for Resilient Operation through Controllable Series Compensation and N-1-1 Contingency Assessment. Energies. 2021; 14(16):4756. https://doi.org/10.3390/en14164756

Chicago/Turabian Style

Huang, Liping, Zhaoxiong Huang, Chun Sing Lai, Guangya Yang, Zhuoli Zhao, Ning Tong, Xiaomei Wu, and Loi Lei Lai. 2021. "Augmented Power Dispatch for Resilient Operation through Controllable Series Compensation and N-1-1 Contingency Assessment" Energies 14, no. 16: 4756. https://doi.org/10.3390/en14164756

APA Style

Huang, L., Huang, Z., Lai, C. S., Yang, G., Zhao, Z., Tong, N., Wu, X., & Lai, L. L. (2021). Augmented Power Dispatch for Resilient Operation through Controllable Series Compensation and N-1-1 Contingency Assessment. Energies, 14(16), 4756. https://doi.org/10.3390/en14164756

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