1. Introduction
The reliable and safe operation of power transmission grids is of paramount importance for the prosperity of modern society. Power outages and interruptions in the U.S. have been estimated to cost around 150 billion dollars per year. The major blackout of Northeast America in 2003 resulted, alone, in a USD 6 billion economic loss for the region [
1,
2]. Moreover, the social consequences of power interruptions, such as those related to transportation, food storage and credit card operations, just to mention a few, are no less serious than economic consequences [
3].
Electrical power blackouts are caused by cascading failures, usually initiated by the failure of a limited set of components, often caused, in turn, by external events such as lightning, ice, and other extreme weather conditions. Other components failures and disconnections may then propagate across the transmission network, due to the following power load re-distribution among the still-functioning components, which may pass their design load capacities and fail or be disconnected to avoid further severe damage.
Various research efforts have focused on the risk of cascading failures in power transmission networks, analyzing temporal series of blackout data [
4] to build probabilistic distributions by developing simulation models to describe cascade dynamics [
5,
6,
7,
8] and by defining and calculating quantitative indicators for identifying system vulnerabilities to cascading failures [
9,
10].
These efforts are also of relevance for the successful development of the new concept of the “smart grid” [
3], which involves the evolution of the current centralized power generation structures into distributed ones that are interconnected via a properly designed and operated information and communication network (ICT) for improving its observability and controllability [
11]. Indeed, the interconnection of the power network to a fast and reliable ICT network allows, in principle, thorough online monitoring of the system with improved coordination of the control and protection systems; automated control strategies may then promptly detect and isolate small contingencies before they give rise to catastrophic cascading failures, or, at the very least, they may contribute to the effective mitigation of their damage [
12,
13,
14]. Understanding of the dynamics of power system cascading failures is essential for devising more reliable, safe, and economically sound infrastructure designs to improve mitigation and control strategies and to provide more accurate risk assessments. The complexity of the cascading failure process in an electrical power grid is mainly due to the concurrence of the physics of power flows and discrete stochastic network topology reconfigurations, as components become disconnected from the grid as a result of the failure cascade.
For cascading failure analysis, the power transmission network is typically modeled as a graph
, where
is the set of vertices (nodes) representing generators, transmission, and load buses and
is the set of edges (arcs) representing the power lines linking the network components. Usually, it is assumed that only elements belonging to one of the sets
or
are subjected to failure; consequently, the analyses are carried out under the assumptions of node removals [
7,
15,
16,
17] or edge removals [
5,
6,
10,
18]. In this work, we restrict our attention to only edge removals, since power transmission line failures are more common than bus failures [
19].
Under this graph-based framework, cascading failure models can be relatively simple, conceptually. First, a capacity is a assigned to each edge , to represent the maximum amount of power flow that can safely flow through line . Then, a single edge is assumed to fail and it is removed from the network, giving rise to a power redistribution transient according to given rules. During the transient, whenever in a line, such a line can fail with probability . When a line fails, it is removed from the set , thus further modifying the topology of the network; then, the power flowing in the network is redistributed again, possibly giving rise to further failures or disconnections.
Typically, two approaches are considered for cascading failure modeling: (i) complex network (abstract) models and (ii) realistic power flow models considering alternate current (AC) or linearized direct current (DC) schemes.
In the first kind of model, such as the Motter–Lai model [
7], a generic flow unit is assumed to travel along the shortest paths joining pairs of nodes, without specification and consideration of electrical parameters; extensions of this model can be found, for example, in Ref. [
6], where the properties of network connection efficiency are introduced for evaluating the flows; in Ref. [
10], where the exceeding loads are propagated on the neighboring components by a simple topological re-dispatch rule; and in Ref. [
20], where a random flow model is proposed. Although simple and computationally not expensive, these models could be too far from the real physical behavior of power systems [
21]. Nevertheless, they have been applied to analyze a broad range of networks ranging from randomly generated networks to some real power transmission networks [
15,
17].
The second class of models aims to overcome the limitations of complex network models by considering a more realistic electrical network model whose power flows are governed by Kirchoff laws. Then, a system of nonlinear equations describing the physics of the alternate current (AC) power flow can be derived [
22]. Under normal operating conditions, a solution of the system can be obtained by means of the Newton–Raphson method. However, under extremely dynamic operating conditions, such as those arising in cascading failures, the Newton–Rapshon approach may be too slow or may even fail to converge, thus hampering any analysis requiring repeated network model evaluations (e.g., for probabilistic risk assessment, uncertainty propagation, or sensitivity analysis). For these reasons, often, the linearized version of the AC power flow is preferred, which assumes a direct current (DC) model and other simplifying assumptions [
22]. The major problem related to the use of this kind of model is that the electrical data of actual power networks, such as generator and load locations, transformer voltage magnitudes, and line impedances and capacities, are available to the scientific community only for a handful of systems, e.g., IEEE test systems (the Power System Test Archive, UWEE). This fact has limited the application of cascading failure models to a restricted set of case studies with very specific topologies and electrical parameters, significantly hindering the generality of the results obtained and of the conclusions drawn. These limitations should not be underestimated, since many researchers have already pointed out the importance of the topological configuration and electrical properties in the assessment of the network behavior with respect to cascading failure [
9,
15,
23] and voltage stability [
24]. Moreover, many of the control and mitigation algorithms presented in the literature [
12,
13,
14] are developed on the basis of the power flow models discussed above; thus, in general, their performances have been tested on the same restricted set of case studies, leaving many unanswered questions about their robustness and efficacy in different operating contexts.
The first objective of this work is that of overcoming some of the limitations of the models described above by exploiting the recently proposed algorithm for generating random artificial power networks characterized by topological and electrical properties in statistical agreement with those of real power transmission systems [
25]. To this end, we first extend the algorithm of Ref. [
25] in order to be able to also sample the power supply and demand locations and magnitudes. The extension allows for further enlargement of the space of grid configurations available for statistical analysis and provides a means for possibly capturing the role played by the uncertain distributions of the loads and generators, which are mainly due to the continuous shift towards the increasingly distributed generation schemes of the new-generation power grids [
3]. Note that the random networks generated with the algorithm in Ref. [
25] do not in general fulfill the
network design basic criterion, which requires that at least two line failures are necessary for a cascading failure event to be triggered [
26]. Then, the algorithm is coupled to a model of cascading failures based on a DC power flow approximation and relying on a minimal, proportional re-dispatch control scheme, which maintains the power balance in each network island formed during the cascades. A dynamic power inertia model taken from Ref. [
12] is also introduced to realistically account for the temporal evolution of the cascading failures.
The proposed computational model allows us to perform statistical analyses on different families of power grids with respect to their cascading failure behaviors.
An additional original contribution of this work, derived from the computer simulations generated with the tool described above, is the definition of new metrics for ranking the control and mitigation effort requirement of individual accidental scenarios and/or of the power grid configurations with respect to cascading failures. In order to do so, the major cascade consequences—i.e., the final load shedding, the number of lines disconnected due to overloads, and the cascade durations—are properly accounted for in the proposed definitions.
Finally, a genetic algorithm-based procedure aiming to optimize the proposed metrics is proposed in order to identify possible strategies for improving the robustness of an existing power network with respect to cascading failures.
The proposed computational approach is demonstrated with regard to random power transmission networks derived from the IEEE14 and the IEEE118 reference grids (Power System Test Archive, UWEE).
The paper is structured as follows.
Section 2 reviews some of the topological properties of the real data available on power transmission networks and briefly explains the RT-nestedSW algorithm of Ref. [
25];
Section 3 recalls the DC power flow model, illustrates the details of the cascading failure model employed throughout the whole work, and introduces the standard measures usually employed to quantify cascading failure damage;
Section 4 illustrates the simulation results and introduces the new metrics for the cascade scenarios and the sampled networks; and finally,
Section 5 concludes the work and outlines possible future research issues.
2. Random Generation of Power Grid Configurations
Previous works on modeling power transmission networks have pointed out the need to generate random power grid test cases of scalable sizes. A few solutions have been proposed for this purpose. For example, Ref. [
27] proposed ring-like structures to study the pattern and the velocity of contingency propagation, and Ref. [
18] used a tree structure network to study critical points and detect transition points in power system blackouts. However, both types of random network generators neither accurately reflect the topologies of real power transmission grids nor account for the electrical properties of the networks (e.g., loads, generation, and impedances), which together determine the behavior of a power transmission system including the dynamic propagation of oscillations and disturbances in the grid. The recent RT-nested Small World (RT-nestedSW) algorithm introduced by Ref. [
25] significantly overcomes both the topological and the electrical deficiencies of previous approaches.
In this section, we recall the main features of the RT-nestedSW algorithm introduced by Ref. [
25] for randomly generating plausible realistic power grid topologies, and we also explain the extension we made in order to be able to generate load and generator schemes with random positions and magnitudes. For further details, the interested reader is referred to Ref. [
25].
From a purely topological point of view, a power network can be represented by an undirected graph
, with
nodes representing the substations in the system (i.e., generators, loads and transmission nodes) and
links representing the transmission lines interconnecting the substations. Often, the topologies of power grids are considered small-world networks [
28]. However, although real power network systems bear some similarities with small-world networks, in Ref. [
25], it was shown that they have significantly better connectivity scaling laws. In fact, the average nodal degree
is basically constant and does not scale with the network size, as in the case of small-world grids [
28]. A detailed analysis of the topological properties of real power transmission networks is out of the scope of this work; for more complete studies of these systems and their relationships with standard random graph models, the interested reader is referred to [
25,
29,
30].
In order to be able to reproduce the random wiring of a power transmission network, Ref. [
25] proposed a new approach based on nesting several small-world sub-networks into a regular lattice. By doing so, it is possible to generate networks characterized by a connectivity lying between those of one-dimensional and two-dimensional lattices with better scaling properties than the standard small-world model. The RT-nested SW algorithm proceeds in a hierarchical way for generating random power grids: first, it produces connected sub-networks, and then it interconnects the sub-networks by means of lattice connections. Finally, it generates line impedances by sampling from proper probability distributions estimated on the basis of the available data.
The generation of the sub-networks stems from an algorithm different from the classical one proposed by Ref. [
28], with differences lying mainly on the link selection and rewiring procedures. With regard to the link selection, instead of generating links by connecting the most immediate
neighboring nodes to form a regular lattice, the RT-nested SW algorithm selects a number
(sampled from a geometric distribution with an expected value
) of links at random from a local neighborhood of
, with
being a properly defined distance threshold. The local neighborhood for node
is defined as the group of nodes with a mutual node index difference less than
:
. With regard to link rewiring, in Ref. [
25], the authors exploited the apparent correlation between the rewires for building a Markov chain with transition probabilities
in order to select clusters of nodes and, therefore, groups of links to be rewired. After the Markov chain is run
times, clusters of nodes are obtained, which are labeled “1” if they have to be rewired or “0” if not. Then, by a specific probability
, some links are selected to rewire from all the links originating from each “1” cluster of nodes, and the corresponding local links are rewired to outside “1” clusters.
In the second step of the algorithm, lattice connections are sampled from neighboring sub-networks to form the main connected network.
Finally, line impedances are randomly generated from a heavy-tailed distribution properly estimated from the data available. The
realizations are then sorted by magnitude in ascending order and grouped into local links, rewire links, and lattice connections, according to the range of values they belong to, as shown in
Figure 1. The ranges are properly defined on the basis of physical considerations. The line impedances in each interval are then assigned randomly to the corresponding group of links in the topology.
It is important to underline that the data used to estimate the distributions adopted in the above procedure are those related to a single, specific real power transmission network. The RT-nested SW algorithm, then, allows us to generate plausible networks that are statistically similar to the grid chosen as the reference. The work of Ref. [
25] also provides the statistical tools needed to extract the required distributions from the reference network.
In particular, in this work, we aim to obtain transmission grids similar to the historical IEEE14 and IEEE118 test cases (from the UWEE Power System Test Case Archive). In order to do so, we fix the number of generators (
) equal to the number of buses in the IEEE system under consideration, which have non-zero power generation (
. Moreover, the nominal power production vector
, where each component
is the power generated by node
in our network, is a simple random permutation of
, the nominal power production vector for the IEEE system under consideration. The total power demand
is, then, randomly allocated among the nodes which have
, so that the power balance constraint within the network is satisfied, in agreement with the DC power flow modeling of the system [
22]:
The portion of
to be assigned to node
with
is sampled from a symmetric Dirichlet distribution in the following way:
where
is a
dimensional vector of all ones.
The maximum power that each generator can produce is assumed to be a fraction
larger than nominal:
where
is the power produced by generator
. The line capacities
,
are assigned by the following procedure:
The power produced by each generator is assumed to be at its maximum: for each ; accordingly, is increased by the same fraction for each in order for (1) to be satisfied.
The corresponding power
flowing in each line
of the transmission network is computed by means of a DC power flow model with no losses, solved by the Matlab function MATPOWER [
31].
The line capacities are then assigned as follows:
where
is a parameter playing the same role of
for the power generation. Alternatively, different strategies for more realistically assigning both the maximum generator power and the line capacities could in principle be followed, based, for example, on sampling values of both
and
from properly devised probability distributions or on calculating them on the basis of physical laws and/or engineering practices. However, we believe this is outside the scope of the present work, which is mainly methodological, and that the assumptions made are sufficient to capture the average properties of cascading failures in power transmission networks. At the end of this procedure, we obtain a “topologically and electrically” characterized power grid, which can be used for simulating cascading failures within a DC power flow approximation of the power network behavior.
In what follows, the symbol will be used to denote the complete power transmission network configuration, i.e., its network topology, the loads, the generators, and the link impedances and capacities.
5. Conclusions
In this work, we have investigated some statistical properties of power transmission networks with the general objective of identifying common strengths and weaknesses with respect to their cascading failure behavior. In order to do so, we have integrated a random power grid generator with a cascading failure simulator based on a DC approximation of the power flows. The integrated algorithm has allowed us to perform a systematic analysis of the cascading failure dynamics within a broad set of networks that are individually different in terms of power grid topologies, loads configurations, generators configurations and line impedances but bear the same statistical properties typical of actual power transmission grids. In particular, in this work, we have referred to both IEEE14- and IEEE118-like networks.
The analysis has led to the identification of an unexpected bi-modal behavior of the cascading failures, leading to load shedding, with respect to the final number of lines disconnected due to overloading. It has been shown that this bi-modality is strictly correlated to the formation of large islands during the cascade propagation, which, under the re-dispatch strategy adopted in this work, amplifies the network damage in terms of load losses.
Then, we developed a new metric for quantifying the control and mitigation requirements of individual scenarios with respect to cascading failures; this metric accounts for the duration of a cascading failure event initiated by a single-line failure and its consequences in terms of final load shedding. By averaging the values of the metric over all possible line failures in a single-grid model, we obtained a new metric for the robustness of a whole network with respect to cascading failures; the standard deviation of the metric can also be interpreted as a metric itself, since it captures the presence of weak lines whose failures lead to anomalous propagating behaviors in power networks, which otherwise appear to be robust. The metric has been shown to be able to correctly identify as critical those randomly generated scenarios and/or configurations that present evident design flaws; however, at the same time, it also identifies those scenarios/configurations with more subtle and unexpected deficiencies, which are otherwise very difficult to capture. The calculation of the proposed metric has also allowed us to perform a rough sensitivity analysis of the features most influencing the cascading failure behavior of power transmission networks, which confirmed the expected importance of the topology and of the correlated impedances’ distribution.
Finally, we have shown that the proposed metrics can be effectively exploited within a genetic algorithm search scheme for the identification of optimal improvements to an existing power grid in terms of both line impedances and loads at the nodes. The new metrics can then be effectively included in practical optimizations, where the necessary modifications to national and over-national power grids must be chosen by taking into account several other possibly conflicting objectives such as economics, congestion issues, political considerations, etc.
There are many other potential uses of the computational approach developed in this work. One of the most promising seems to be that of automatically generating and selecting critical scenarios/configurations to be used as worst-case scenarios for testing, validating, and improving the robustness of the control or mitigation strategies to be adopted during cascading failure events; control and mitigation effort metrics, such as those introduced in this work, may serve to compare the performances of different controllers, possibly also accounting for associated economic requirements.
Current and future research efforts are focused on strengthening the applicability of the proposed approach, which is mainly methodological, to practical problems. This entails, first of all, modifying the random network sampling algorithm in order to be able to generate only secure configurations and, secondly, accessing more powerful and parallelized computational resources in order to be able to study the dynamics of more realistic large power transmission networks.