Power system studies require the topological structures of real-world powernetworks; however, such data is confidential due to important securityconcerns. Thus, power grid synthesis (PGS), i.e., creating realistic powergrids that imitate actual power networks, has gained significant attention. Inthis letter, we cast PGS into a graph distribution learning (GDL) problem wherethe probability distribution functions (PDFs) of the nodes (buses) and edges(lines) are captured. A novel deep GDL (DeepGDL) model is proposed to learn thetopological patterns of buses/lines with their physical features (e.g., powerinjection and line impedance). Having a deep nonlinear recurrent structure,DeepGDL understands complex nonlinear topological properties and captures thegraph PDF. Sampling from the obtained PDF, we are able to create a large set ofrealistic networks that all resemble the original power grid. Simulationresults show the significant accuracy of our created synthetic power grids interms of various topological metrics and power flow measurements.
Quick Read (beta)
Deep Generative Graph Distribution Learning
for Synthetic Power Grids
Power system studies require the topological structures of real-world power networks; however, such data is confidential due to important security concerns. Thus, power grid synthesis (PGS), i.e., creating realistic power grids that imitate actual power networks, has gained significant attention. In this letter, we cast PGS into a graph distribution learning (GDL) problem where the probability distribution functions (PDFs) of the nodes (buses) and edges (lines) are captured. A novel deep GDL (DeepGDL) model is proposed to learn the topological patterns of buses/lines with their physical features (e.g., power injection and line impedance). Having a deep nonlinear recurrent structure, DeepGDL understands complex nonlinear topological properties and captures the graph PDF. Sampling from the obtained PDF, we are able to create a large set of realistic networks that all resemble the original power grid. Simulation results show the significant accuracy of our created synthetic power grids in terms of various topological metrics and power flow measurements.
.tifpng.pngconvert #1 \OutputFile \AppendGraphicsExtensions.tif \bstctlciteIEEEexample:BSTcontrol
MANY power system studies require the actual topologies of real-world power networks with real bus and line coordinates as well as real physical charactristics such as power injections and line impedance values. However, due to the confidentiality concerns, there are very few datasets of real power grids such as the IEEE test cases  and the Polish grid  that are publicly available for such studies. Even these datasets do not contain important information regarding the geographical locations of buses, and lack the diverse characteristics of real topologies. In recent studies, several works are focused on the spatial models for power networks that are realistic but synthesized without revealing any confidential information. [3, 4] create power networks based on the locations of cities and power plants in Texas; however, no geographical or performance comparison with the actual power grid is provided to justify the approach. Although [3, 4] contain useful engineering details, they do not devise a general algorithm that can work with any realistic power network; thus, in practice, these approaches are limited. In this line of research,  is the only learning-based power grid synthesis (PGS) model presented in recent literature that applies a Gaussian Mixture Model (GMM) to compute the distributions of bus locations in the Western Interconnect network. However, due to lack of generalization, GMM is theoretically incapable of learning the physical properties of power grid components. Moreover, the data required by  exceeds the power grid data as this model represents the demand based on the average residential power usage and city populations. Also, similar to  and , the GMM approach  puts strict assumptions on the node degree and line distributions which do not necessarily hold.
This letter presents a novel deep graph distribution learning (DeepGDL) algorithm to create synthetic power grids. We present a recurrent model that efficiently captures node/edge distributions of arbitrary complex networks. Leveraging a deep architecture, DeepGDL captures complex nonlinear manifolds of nodes/edges in the actual power network, and generates synthetic power grids that imitate the original network. To the best of our knowledge, DeepGDL is the only deep learning-based PGS solution. Also, in contrast to all previous works, the presented algorithm needs no additional information other than the network topology and the bus/line physical properties; hence, it can be applied as a general framework to any PGS problem with minimum amounts of data.
II DeepGDL for Power Grid Synthesis
Let us consider the actual power grid obtained from the Columbia University Synthetic Power Grid (CUSPG) dataset  that contains 14430 buses and 18884 lines located in Western US. We represent this network as an undirected weighted graph where denotes the -th node while an edge is a transmission line connecting two nodes and . For each node , contains four features including the latitude , longitude , power supply (i.e., power generated at the bus ), and power demand . A weight is defined to represent the reactance of the line corresponding to each . Our objective is to learn the PDFs of nodes as well as edges , and sample from them to generate synthetic power grids.
Algorithm 1 is the pseudocode of our presented model. First, using the modularity optimization in , is decomposed into communities (i.e., subgraphs) where the nodes inside each are densely connected while the nodes in distinct communities have sparse connections. For each , the adjacency matrix is modeled based on the node ordering that maps any to a row/column index in . Having the set of all node permutations denoted by , one can write a new adjacency matrix defined by ; hence, we model as an arbitrarily-sized sequence of node orderings each with a unique nodal and edge feature sequence. Considering the ordering , a nonlinear function maps to a 2-dimensional node feature tensor using the following formulation:
where is the nodal feature of the -th node w.r.t the ordering . Similarly, we define that maps to a 2-dimensional edge feature tensor computed by:
where is the weighted adjacency vector of node ; thus, for an undirected weighted , the tensor uniquely defines the graph . As a result, one can recover by the mapping ; hence, in order to learn the distribution of the nodes and edges in communities denoted by , we maximize the likelihood of their observation in actual network by:
where if is true and zero otherwise. Using this optimization, we maximize the observation of ’s communities, hence finding a model that can imitate .
Let us write a decomposition for in (3) using conditional probabilities of the observed sequence of nodal and edge features implied by for any community :
The recursive characteristics of the graph features in (4) motivates us to model nodal and edge features by a recurrent neural network. Therefore, we model using a deep Gated Recurrent Unit (GRU) network. For each community , at each iteration of our deep recurrent model, the GRU computes a graph state using the previous node/edge observations by the following recurrent structure:
Here, and in are the intermediate nodal and edge features of observed at the -th round of the recursive formulation. At each round , the temporal features, and , are computed for the graph observation using a sigmoid function with weights and , respectively. Then, the GRU’s updating vector is computed using the intermediate graph features, and , as well as the GRU state at the previous round . Here, denotes a tangent hyperbolic function with weight applied to compute . The GRU state is finally computed as a linear regression of and weighted by the temporal feature . After running the recursion (5) for rounds, the final GRU state is our extracted graph state .
At each iteration , we model the node/edge distribution of each community , and generate a new synthetic community corresponding to using the following formulations:
|Topological Properties||Power Flow Statistics (MW)|
|171.18 (11.34)||30.42 (1.20)||307.18 (15.59)||4,590.25 (217.01)|
where and are two sigmoidal neural networks (SNNs) with the input and output vectors and , respectively. The vector encodes the nodal feature distribution of the -th node in . Similarly, encodes the edge weight distribution corresponding to this node. Sampling from these two distributions in (6), we create a new node (i.e., the i-th synthesized node) of . We run (6) to create new nodes and edges until the GRU reaches a predefined state, which terminates the graph generation process and outputs for every community . Each of the generated communities is a component of our final synthetic power network .
As shown in , in dense regions of real-world power grids, the low-degree buses are generally connected to high-degree buses. Thus, in order to merge the generated components and form a single-component robust network , we connect the low-degree nodes in dense regions, to high-degree nodes until the becomes single-component. The resulting is finally reported as the output.
III Simulation Results
First, we set several hyperparameters defined in the DeepGDL algorithm. The GRU state variables , and , are considered to have dimensions. Both and are modeled as a SNN with neurons in its first layer and neurons in the second layer. The gradient descent is considered to have a learning rate of while its batch size is . In this study, graph communities are considered for . We generate synthetic power grids via DeepGDL. Fig. 1(a) and (b) depict the actual network and one of our synthetic networks , respectively. Each community is shown by a distinct color. Moreover, Fig. 1(c) and (d) show the amount of power supply at each bus (node) in the two networks. To compare the performance of the generated grids with the actual grid, we conduct cascading failure experiments using DC power flow. Fig. 1(e) and (f) demonstrate the Yield (i.e. ratio of demand supplied after cascades to the actual value of demand) in the two power grids for cascades initiatd in regions of radius km. Each point shows the computed Yield at the center of its corresponding region. As shown in Fig. 1, accurately resembles the actual grid regarding the topological characteristics and power flow measurements. Table I reports detailed topological measurements and power flow statistics obtained by , DeepGDL algorithm, and the state-of-the-art GMM  model. It is shown that DeepGDL has significant improvements over GMM in terms of average node degree , average path length , network Diameter , Density , Modularity , and Average Clustering Coefficient . Also, Table I shows the substantial accuracy improvements of DeepGDL over the state-of-the-art in terms of power flow measurements.
This letter presents a deep graph distribution learning algorithm for the problem of power grid synthesis. A novel recurrent model is proposed to capture deep nodal and edge features from a realistic power network. To the best of our knowledge, DeepGDL is the only deep learning model that can synthesize large-scale networks. Simulation results show significant accuracy of DeepGDL in terms of topological measurements and power flow analysis metrics.
-  “IEEE benchmark systems,” www.labs.ece.uw.edu/pstca/, accessed: 2019-01-20.
-  “Polish grid,” www.pserc.cornell.edu/matpower/, accessed: 2019-01-21.
-  K. M. Gegner, A. B. Birchfield, T. Xu, K. S. Shetye, and T. J. Overbye, “A methodology for the creation of geographically realistic synthetic power flow models,” in Power and Energy Conference at Illinois (PECI), 2016 IEEE, 2016, pp. 1–6.
-  A. B. Birchfield, T. Xu, K. M. Gegner, K. S. Shetye, and T. J. Overbye, “Grid structural characteristics as validation criteria for synthetic networks,” IEEE Transactions on power systems, vol. 32, no. 4, pp. 3258–3265, 2017.
-  S. Soltan, A. Loh, and G. Zussman, “A learning-based method for generating synthetic power grids,” IEEE Systems Journal, no. 99, pp. 1–10, 2018.
-  S. Soltan, A. Loh, and g. Zussman, “Columbia university synthetic power grid with geographical coordinates,” Tech. Rep., 1 2018.
-  X. Lu, K. Kuzmin, M. Chen, and B. K. Szymanski, “Adaptive modularity maximization via edge weighting scheme,” Information Sciences, vol. 424, pp. 55–68, 2018.