The present disclosure generally relates to utility distribution networks, and in particular, to a system and associated method for time-synchronized topology and state estimation in real-time unobservable utility distribution systems.
State and topology estimation in utility distribution systems is imperative to ensure safe and effective operation. One challenge with this task is that a utility distribution network may not be equipped with a sufficient quantity of synchrophasor measurement devices and/or other sensors required for classical state estimation.
It is with these observations in mind, among others, that various aspects of the present disclosure were conceived and developed.
The present patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee.
Corresponding reference characters indicate corresponding elements among the view of the drawings. The headings used in the figures do not limit the scope of the claims.
A system and associated methods for deep neural network (DNN)-based distribution system state estimation (DSSE) of real-time unobservable systems are disclosed herein. The system aims to perform DNN-based DSSE under varying topologies using transfer learning.
Real-time monitoring and control of distribution networks was traditionally deemed unnecessary because it had radial configuration, unidirectional power flows, and predictable load patterns. However, the fast growth of behind-the-meter (BTM) generation, particularly solar photovoltaic (PV), electric vehicles, and storage, is transitioning the distribution system from a passive load-serving entity to an active market-ready entity, whose reliable and secure operation necessitates real-time situational awareness. Phasor measurement units (PMUs), distribution-PMUs (D-PMUs), and/or micro-PMUs, collectively referred to as synchrophasor measurement devices (SMDs) in this disclosure, have been introduced into the distribution system to provide fast (sub-second) situational awareness by enabling time-synchronized estimation. However, due to the high cost of installation, the numbers of SMDs in a typical distribution network are not large enough to provide an independent assessment of the system state. The assumption of Gaussian noise in synchrophasor measurements has also been disproved recently.
At the same time, modern distribution systems are being equipped with advanced metering infrastructure (AMI) in the form of smart meters. Hence, a combination of smart meter data with SMD data for facilitating distribution system state estimation (DSSE), particularly in three-phase unbalanced distribution systems, has also been tried out. However, smart meters typically measure energy consumption from 15 minute to hourly time intervals and report their readings after a few hours or even a few days. These two aspects make smart meter data unsuitable for real-time DSSE. Moreover, smart meter data are not time-synchronized by default, which makes their direct integration with SMD data a statistical challenge. Real-time DSSE is possible using SMDs alone. However, as mentioned above, placing large numbers of SMDs in a distribution network is cost prohibitive.
Lastly, topology of a distribution network changes with time. This implies that impacts of changes in the configuration of the distribution system must be considered when performing DSSE. Although approaches have been proposed in the past to estimate topology and state estimates independently, a joint framework for identifying the topology and estimating the states at high-speed in distribution systems that are unobservable by SMDs has not been proposed before.
One implementation of the system is applied to estimate voltages in distribution systems under extremely unobserved conditions (e.g., with few synchrophasor measurement devices available). The system can represent voltage phasors in an unbalanced distribution system at node g as xg=[xg1, xg2, xg3]T where the superscripts are phase indices, and xgk=Vgk<θgk, where Vgk and θgk are voltage phase magnitude and voltage phase angle, respectively, of phase k at node g. The overall system state x=[x1, . . . , xn]T is the column vector that includes voltage phasors of all the nodes, where the unbalanced distribution system includes a plurality of nodes (e.g., n total nodes where g∈[1, n]). Similarly, all available measurements can be stacked to create the column vector z=[z1, . . . , zm]T (e.g., for m total sensors). In classical state estimation, the relationship between the measurements and the states is given by:
where h(⋅) is the measurement function, whose structure depends on the type of sensor(s) providing the measurements as well as the observability of the unbalanced distribution system by the sensor(s), and e represents the noise in the measurements. Time-synchronized state estimation in distribution networks using classical approaches, such as least-squares, can be formulated as shown in (2) below.
Classical state estimation requires the distribution system to be completely observed by synchrophasor measurement devices (SMDs). In other words, m>>n. However, it is highly unlikely that, at least in the near future, a distribution system will be equipped with as many SMDs as is required for complete observation of the distribution system by them alone. To circumvent the problem of scarcity of SMDs and other sensors for performing time-synchronized DSSE, the system disclosed herein applies a Bayesian approach in which a state, x, and a measurement, z, are treated as random variables. A minimum mean squared error (MMSE) estimator is applied to minimize the estimation error as shown in (3) below.
The MMSE estimator directly minimizes the estimation error while classical estimators, such as least-squares, minimize the modeling error embedded via the measurement function h(⋅) that relates the measurements with the states. By circumventing the need for the measurement function using the Bayesian MMSE state estimator, the system disclosed herein bypasses real-time observability requirements. Note that in (3), there are two underlying challenges to computing the conditional mean. First, the conditional expectation, which is defined by:
requires the knowledge of p(x|z), which is the joint probability density function (PDF) between x and z. When the number of SMDs is scarce, the PDF between SMD data and all the voltage phasors (i.e., the states of the distribution system) is unknown or impossible to specify. This makes direct computation of {circumflex over (x)}*(z) intractable. Second, even if the underlying joint PDF is known, finding a closed-form solution for (4) can be difficult. A DNN can be used to approximate the MMSE state estimator as a DNN has excellent approximation capabilities; e.g., the DNN for DSSE developed here finds a mapping (denoted by (⋅)) that relates x and z.
Even if a DNN can successfully approximate (x|z) for a given topology, once the topology changes, the distribution of the inputs (i.e., SMD measurements) for which the DNN had been trained for, change. This is best realized from the fact that the direction of the currents in a distribution feeder can reverse when topology change occurs due to status changes of the switches. Thus, there is a need to update the trained DNN after topology changes. One way to do this is to train the DNN for DSSE afresh for every new topology. However, doing so may take a very long time. An alternate (better) solution is to employ transfer learning to transfer the knowledge gained from the old topology to the new topology.
Transfer learning improves the learning of the target prediction function in the target domain using the knowledge available in the source domain and task. A domain comprises two parts: a feature space, , and a marginal probability distribution, P(z). Given , a task includes two parts: a label space, , and an objective prediction (mapping) function, (⋅). In DNN-based DSSE under varying topologies, does not change as the same SMD measurements will be used for different topologies. However, P(z) changes because loads are served by different paths when topology changes, i.e., when S≠T. Similarly, does not change because the number of states (i.e., voltage phasor at each node) and their nature are the same. However, (⋅), must be retrained for the target domain, i.e., ≠. Now, it is clear from this problem set-up that the objective is to induce transfer of knowledge gained from S and (old topology) to T and (new topology). Thus, the system can employ inductive transfer learning to attain the desired objective. Four approaches have been proposed for implementing inductive transfer learning, namely, feature-representation transfer, instance transfer, relational-knowledge transfer, and parameter transfer. Here, the system can apply parameter transfer to update the DNN for DSSE as the DNN's parameters can be used for multiple domains.
Parameter-based transfer learning methods include parameter-sharing and fine-tuning. Parameter-sharing assumes that the parameters are highly transferable due to which the parameters in the source domain can be directly copied to the target domain, where they are kept “frozen”. Fine-tuning assumes that the parameters in the source domain are useful, but they must be trained with limited target domain data to better adapt to the target domain. Since there is no guarantee that the parameters of the DNN-based DSSE will be highly transferable for different topologies, fine-tuning is used here to update the weights of the DNN for DSSE when topology changes. Essentially, fine-tuning provides a more effective initialization (than random initialization) by using the weights from the previously well-trained DNN. By doing this, the system bypasses the need for large amounts of data (and time) for DNN re-training.
The basic structure of the second DNN 120 for DSSE is shown in
The second DNN 120 for DSSE shown in
As opposed to the second DNN 120 that was built for DSSE and apples regression-based operations to estimate states, the first DNN 110 (shown in
It is possible that a distribution system may not have any SMD to begin with or the number of SMDs may not be sufficient to perform DNN-based TI and/or DSSE. For such situations, the following algorithm is proposed to place SMDs in the system. In this algorithm, Budget refers to the budget allocated for SMD placement, TIaccuracy and DSSEaccuracy are the minimum desired accuracy for TI and DSSE, respectively, DSSEcorr is the minimum correlation between each pair of nodes, and M is the number of nodes in the system. For placing SMDs for DNN-based TI, a greedy search algorithm called sequential forward selection was used, while the correlation between the nodes, which was needed for placing SMDs for DNN-based DSSE, was determined using Spearman correlation coefficient (SCC).
An overview of a process 200 applied by the system 100 disclosed herein is shown in
The process 200 includes: receiving a set of real-time SMD measurements from a distribution system; applying the set of real-time SMD measurements as input to a first deep neural network for topology identification of the distribution system to determine a topology of the distribution system; and applying the set of real-time SMD measurements as input to a second deep neural network for distribution system state estimation of the distribution system, where an output of the second deep neural network includes a set of estimated states of the distribution system based on the set of real-time SMD measurements. The first deep neural network is a classification-based deep neural network configured for determining the topology of the distribution system based on the set of real-time SMD measurements. The second deep neural network is a regression-based deep neural network configured for estimating a state and a measurement of a node of the distribution system based on the set of real-time SMD measurements.
During real-time operation, the process 200 can include: determining, based on the output of the first deep neural network, if a topology of the distribution system is consistent with a training topology used to train the second deep neural network for distribution system state estimation; and updating the second deep neural network using a transfer learning methodology prior to application of the set of real-time SMD measurements as input to the second deep neural network.
The process 200 can also include: training the first deep neural network for topology identification and the second deep neural network for distribution system state estimation through an offline learning process that iteratively models a set of modeled SMD measurements for one or more possible topologies of the distribution network and trains the first deep neural network and the second deep neural network based on the set of modeled SMD measurements for one or more possible topologies of the distribution network. This step can further include: receiving a set of historical smart meter data for the distribution network; modeling unbalanced phase power flow for an i-th topology of the distribution network; modeling, for the i-th topology, placement of one or more modeled SMDs at one or more nodes of the distribution network; and obtaining, for the i-th topology, a set of modeled SMD measurements including modeled voltage phasor measurements and modeled current phasor measurements representing SMD measurements captured by one or more modeled SMDs at the one or more nodes of the distribution network. Based on the topology being modeled and the set of modeled SMD measurements, the process 200 can further include at least one of the following: training the second deep neural network based on the set of modeled SMD measurements for the i-th topology, the i-th topology being a base topology of the distribution network; training the first deep neural network based on the set of modeled SMD measurements for the i-th topology, the i-th topology being a last possible topology of the distribution network, and storing the set of modeled SMD measurements in association with the i-th topology.
Modeling placement of one or more modeled SMDs at one or more nodes of the distribution network includes: applying (for TI-based node selection) a sequential forward selection methodology to select the one or more nodes of the distribution network for SMD placement. Alternatively (or jointly), modeling placement of one or more modeled SMDs at one or more nodes of the distribution network includes: receiving a minimum correlation value representative of target minimum correlation between one or more nodes of the distribution network; determining a Spearman Correlation Coefficient value between each respective modeled voltage phasor measurement; constructing one or more Spearman Correlation Matrices based on the Spearman Correlation Coefficient value between each respective modeled voltage phasor measurement; applying a hierarchical clustering methodology to the one or more Spearman Correlation Matrices to obtain one or more clusters; and selecting a common node of the one or more nodes of the distribution network in each cluster of the one or more clusters for SMD placement.
Iteratively obtaining modeled SMD measurements can include: obtaining a set of modeled error-free SMD measurements including modeled error-free voltage phasor measurements and modeled error-free current phasor measurements; and adding measurement noise to the set of modeled error-free SMD measurements to obtain the set of modeled SMD measurements including the modeled voltage phasor measurements and the modeled current phasor measurements.
In this study, the DNN-based DSSE was tested for a real-world distribution system for which historical smart meter data was available at all the nodes and SMD measurements were only possible at the feeder-head. Note that most power utilities have real-time measurements only at the feeder-head of their distribution systems. Therefore, it was of interest to evaluate the performance of the proposed method for DSSE when SMD placement cannot be done due to budget constraints and only the existing measurements can be used. There are 648, 665, and 637 nodes in phase A, phase B, and phase C, respectively, of this distribution system, whose voltage magnitudes and angles must be estimated for different operating conditions. Additionally, this feeder is a renewable-rich distribution system with 766 household/commercial rooftop solar PV units. Therefore, this is an ideal distribution system to evaluate the performance of the proposed method for performing DSSE for a large-scale real-world distribution system with high penetration of renewable resources.
The implementation procedure shown in
The reasonably good results reported in Table I demonstrate the applicability of the process 200 for large real-world systems. Hence, power utilities, which do not own enough SMDs to fully observe their distribution systems, can use the process 200 to perform time-synchronized DSSE.
2. Implementation of Transfer Learning and DNN-Based DSSE when Topology Changes for a Real-World Distribution System Located in the U.S. Midwest (See
This distribution system did not have any SMDs placed apriori. However, it has nine switches, which resulted in 84 feasible topologies. Therefore, an integrated algorithm for SMD placement was used to determine the suitable locations for placing SMDs in this distribution system for performing DNN-based TI and DSSE. Six locations were identified as shown with green ovals and red arrows in
Next. It was assumed that the topologies of this distribution system are changed in accordance with Table III.
When topology changes occur, after correctly identifying the new topology using DNN-based TI (note that the DNN-based TI had 99.19% accuracy for this distribution system), the DNN trained for doing DSSE for the old topology, must be updated. As described in the process 200 (see
Lastly, the angle MAE results are now compared with and without fine-tuning of the DNN for DSSE. It is observed from
Device 300 comprises one or more network interfaces 310 (e.g., wired, wireless, PLC, etc.), at least one processor 320, and a memory 340 interconnected by a system bus 350, as well as a power supply 360 (e.g., battery, plug-in, etc.).
Network interface(s) 310 include the mechanical, electrical, and signaling circuitry for communicating data over the communication links coupled to a communication network. Network interfaces 310 are configured to transmit and/or receive data using a variety of different communication protocols. As illustrated, the box representing network interfaces 310 is shown for simplicity, and it is appreciated that such interfaces may represent different types of network connections such as wireless and wired (physical) connections. Network interfaces 310 are shown separately from power supply 360, however it is appreciated that the interfaces that support PLC protocols may communicate through power supply 360 and/or may be an integral component coupled to power supply 360.
Memory 340 includes a plurality of storage locations that are addressable by processor 320 and network interfaces 310 for storing software programs and data structures associated with the embodiments described herein. In some embodiments, device 300 may have limited memory or no memory (e.g., no memory for storage other than for programs/processes operating on the device and associated caches). Memory 340 can include instructions executable by the processor 320 that, when executed by the processor 320, cause the processor 320 to implement aspects of the system and the methods outlined herein.
Processor 320 comprises hardware elements or logic adapted to execute the software programs (e.g., instructions) and manipulate data structures 345. An operating system 342, portions of which are typically resident in memory 340 and executed by the processor, functionally organizes device 300 by, inter alia, invoking operations in support of software processes and/or services executing on the device. These software processes and/or services may include DSSE processes/services 390, which can include aspects of system 100 and process 200 described herein and/or implementations of various modules described herein. Note that while DSSE processes/services 390 is illustrated in centralized memory 340, alternative embodiments provide for the process to be operated within the network interfaces 310, such as a component of a MAC layer, and/or as part of a distributed computing network environment.
It will be apparent to those skilled in the art that other processor and memory types, including various computer-readable media, may be used to store and execute program instructions pertaining to the techniques described herein. Also, while the description illustrates various processes, it is expressly contemplated that various processes may be embodied as modules or engines configured to operate in accordance with the techniques herein (e.g., according to the functionality of a similar process). In this context, the term module and engine may be interchangeable. In general, the term module or engine refers to model or an organization of interrelated software components/functions. Further, while the DSSE processes/services 390 is shown as a standalone process, those skilled in the art will appreciate that this process may be executed as a routine or module within other processes.
It should be understood from the foregoing that, while particular embodiments have been illustrated and described, various modifications can be made thereto without departing from the spirit and scope of the invention as will be apparent to those skilled in the art. Such changes and modifications are within the scope and teachings of this invention as defined in the claims appended hereto.
This is a U.S. non-provisional patent application that claims benefit to U.S. Provisional Application Ser. No. 63/459,791, filed on Apr. 17, 2023, which is herein incorporated by reference in its entirety.
This invention was made with government support under DE-AR0001001 awarded by the Department of Energy. The government has certain rights in the invention.
Number | Date | Country | |
---|---|---|---|
63459791 | Apr 2023 | US |