Robert Qiu1,2,3 Lei Chu2,3 Xing He2,3 Zenan Ling2,3 and Haichun Liu2,3
1 Tennessee Technological University, Cookeville, TN, 38505, USA
2 Department of Electrical Engineering, Shanghai Jiaotong University, Shanghai, 200240, China
3 Research Center for Big Data Engineering and Technology, State Energy Smart Grid Research and Development Center, Shanghai, China
A cornerstone of the smart grid is the advanced monitorability on its assets and operations. Increasingly pervasive installation of phasor measurement units (PMUs) allows the so‐called synchrophasor measurements to be taken roughly 100 times faster than the legacy supervisory control and data acquisition (SCADA) measurements, time‐stamped using the Global Positioning System (GPS) signals to capture the grid dynamics. On the other hand, the availability of low‐latency two‐way communication networks will pave the way to high‐precision real‐time grid state estimation and detection, remedial actions upon network instability, and accurate risk analysis and post‐event assessment for failure prevention.
In this chapter, we firstly model spatiotemporal PMU data in large‐scale grids as random matrix sequences. Secondly, some basic principles of random matrix theory (RMT), such as asymptotic spectrum laws, transforms, convergence rate, and free probability, are introduced briefly in order to improve the understanding and application of RMT technologies. Lastly, the case studies based on synthetic data and real data are developed to evaluate the performance of the RMT‐based schemes in different application scenarios (i.e., state evaluation and situation awareness).
The modern power grid is one of the most complex engineering systems in existence; the North American power grid is recognized as the supreme engineering achievement in the 20th century [1]. The complexity of the future's electrical grid is ever increasing: (1) the evolution of the grid network, especially the expansion in size; (2) the penetration of renewable/distributed resources, flexible/controllable electronic units, or even prosumers with dual load‐generator behavior [2]; and (3) the revolution of the operation mechanism, e.g., demand‐side management. Also, financial, environmental, and regulatory constraints are pushing the electrical grid toward its stability limit.
Generally, power grids have experienced three ages—G1, G2, and G3 [3]. The network structures are depicted in Figure 23.1 [4]. Their data flows and energy flows, as well as corresponding data management systems and work modes, are quite different [5], which are shown in Figure 23.2 and Figure 23.3, respectively.
G1 was developed from the power system around 1900 to 1950, featured by small‐scale isolated grids. For G1, units interchange energy and data within the isolated grid to keep generation‐consumption balance. The units are mostly controlled by themselves, i.e., operating under individual‐work mode. As shown in Figure 23.1(a), each apparatus collects designated data and makes corresponding decisions only with its own application. The individual‐work mode works with an easy logic and little information communication. However, it means few advanced functions and inefficient utilization of resources. It is only suitable for small grids or isolated islands.
G2 was developed from power grids about 1960 to 2000, featured by zone‐dividing large‐scale interconnected grids. For G2, units interchange energy and data with adjacent units. The units are dispatched by a control center, i.e., they are operating under team‐work mode. The regional team leaders, such as local dispatching centers, substations, and microgrid control centers, aggregate their own team members (i.e., units in the region) into a standard black‐box model. These standard models will be further aggregated by the global control center for control or prediction purposes. The two aggregations above are achieved by four steps: data monitoring, data preprocessing, data storage, and data processing. The description above can be summarized by by dotted blue lines in Figure 23.3. In general, the team‐work mode conducts model‐based analysis and mainly concerns system stability rather than individual benefit; it does not work well for smart grids with 4Vs data.
The development of G3 was launched at the beginning of the 21st century, and for China, it is expected to be completed around 2050 [ 3]. Figure 23.1(c) shows that the clear‐cut partitioning is no longer suitable for G3, as well as the team‐work mode, which is based on the regional leader. For G3, the individual units, rather than the regional center (if it still exists), play a dominant role. They are appropriately self‐controlled with high intelligence, resulting in much more flexible flows for both energy exchange and data communication [6]. Accordingly, the group‐work mode is proposed. Under this mode, the individuals freely operate under the supervision of the global control centers [ 5]. VPPs [7], MMGs [8], for instance, are typically G3 utilities. These group‐work mode utilities provide a relaxed environment to benefit both individuals and the grids: the former, driven by their own interests and characteristics, are able to create or join a relatively free group to benefit mutually from sharing their own superior resources; meanwhile, these utilities are often big and controllable enough to be good customers or managers to the grids.
Data are more and more easily accessible in smart grids. Figure 23.4 shows numerous data sources: information communication technology (ICT), advanced metering infrastructure (AMI), supervisory control and data acquisition (SCADA), sensor technology (ST), phasor measurement units (PMUs), and intelligent electronic devices (IEDs) [9]. Hence, data with features of volume, velocity, variety, and veracity (i.e., 4Vs data) [10] are inevitably generated and daily aggregated. Particularly, the “4Vs” are elaborated as follows:
As mentioned above, smart grids are always huge in size and complex in topology; big data analytics and a data‐driven approach become natural solutions for the future grid [13–16]. Driven by data analysis in high dimension, big data technology works out data correlations (indicated by statistical parameters) to gain insight to the inherent mechanisms. Actually, big data technology has already been successfully applied as a powerful data‐driven tool for numerous phenomena, such as quantum systems [17], financial systems [18, 19], biological systems [20], as well as wireless communication networks [21–23]. For smart grids, the data‐driven approach and data utilization are currently important topics, as evidenced in the special issue of “Big Data Analytics for Grid Modernization” [24]. This special issue is most relevant to our book in spirit. Several SA topics are discussed as well. We highlight anomaly detection and classification [25, 26], social media such as Twitter in [27], the estimation of active ingredients such as PV installations [28, 29], and finally the real‐time data for online transient stability evaluation [30]. In addition, we point out research about the improvement in wide‐area monitoring, protection and control (WAMPAC) and the utilization of PMU data [31–34], together with the fault detection and location [35–37]. Xie et al., based on principal component analysis (PCA), proposes an online application for early event detection by introducing a reduced dimensionality [38]. Lim et al. studies the quasi‐steady‐state operational problems relevant to the voltage instability phenomena [39]. These works provide primary exploration of the big data analysis in the smart grid. Furthermore, a brief account for random matrix theory (RMT), which can be seen as basic analysis tools for spatial‐temporal grid data processing, is elaborated in the following subsection.
The last two decades have seen the rapid growth of RMT in many science fields. Brilliant mathematical work in RMT shed light on the challenges of classical statistics. In this subsection, we present a brief introduction to the main development of RMT. The application‐related account, with particular attention paid to recently rising RMT‐based technology that are relevant for the smart grid, is elaborated in Section 23.3.
The research of random matrices began with the work of Wishart in 1928, which focused on the distribution of the sample covariance matrices. The first asymptotic results on the limit spectrum of large random matrices (energy levels of nuclei) were obtained by Wigner in 1950s in a series of works [40–43] which ultimately lead to the well‐known semicircle law [44]. Another breakthrough was presented in [45], which studied the distribution of eigenvalues for empirical covariance matrices. Based on these excellent works, RMT became a vibrant research direction of its own. Plenty of brilliant works that branched off the early physical and statistical applications were put forward in the last decades. For the sake of brevity, here we only show two remarkable results that turned out to be related to a large number of research hotspots in economics, communications, and the smart grid. One of the most striking advances is the discovery of the Tracy Widom distribution of extreme eigenvalues and another one is the single ring law, which described the limit spectrum of eigenvalues of non‐normal square matrices [46]. Interested readers are referred to monographs [47–49] for more details.
We will end this section by providing the structure of the remainder of this chapter.
Firstly, Section 23.2 gives a tutorial account of existing mathematical works that are relevant to the statistical analysis of random matrices arising in smart grids. Specially, Section 23.2.1 introduces data collected from the widely applied phasor measurement unit and data modeling using linear and nonlinear combinations of random matrices. Section 23.2.2 focuses on asymptotic spectrum laws of the major types of random matrices. Section 23.2.3 presents three dominant transforms that play key roles in describing the limit spectra of random matrices. Recent results on the convergence rate to the asymptotic limits are contained in Section 23.2.4. Section 23.2.5 is dedicated to free probability theory, which is demonstrated as a practical tool for smart grids.
Secondly, we begin with some representative problems arisen from the wide deployment of synchronous phasor measurement units that capture various features of interest in smart grids. We then show how random matrix theory has been used to characterize the data collected from synchronous phasor measurement and tackle the problems in the era of big data. In particular, Section 23.3.1 provides some basis hypothesis tests that remain fundamental to research into the behavior of the data in smart grids. Section 23.3.2 concerns stability assessment from some recently developed data‐driven methods based on RMT. Section 23.3.3 focuses on situation awareness for smart grids from linear eigenvalue statistics. The early event detection problem is studied in detail using free probability in Section 23.3.4.
In this section, we provide comprehensive existing mathematical results that are associated with the analysis of statistics of random matrices arising in smart grids. We also describe some new results on random matrices and other data‐driven methods that were inspired by problems of engineering interest.
Before the comprehensive utilization of the RMT framework, we try to build a model for spatiotemporal PMU data using large dimensional random matrices.
It is well accepted that the transient behavior of a large electric power system can be illustrated by a set of differential and algebraic equations (DAEs) as follows [50]:
where are the power state variables, e.g., rotor speeds and the dynamic states of loads, represent the system input parameters, define algebraic variables, e.g., bus voltage magnitudes, denote the time‐invariant system parameters. , and are the sample time, number of system variables, and bus, respectively. The model‐based stability estimators [51, 52] focus on linearization of nonlinear DAEs in (23.1) and (23.2) which gives
where , are the Jacobian matrices of with respect to and . is a diagonal matrix whose diagonal entries equal and is the correction time of the load fluctuations. denotes a diagonal matrix whose diagonal entries are nominal values of the corresponding active or reactive of loads; is assumed to be a vector of independent Gaussian random variables.
It is noted that estimating the system stability by solving the equation (23.3) is becoming increasingly more challenging [ 49] as a consequence of the steady growth of the parameters, say, , , and . Besides, the assumption that follows a Gaussian distribution would restrict the practical applications.
On the other hand, as a novel alternative, the lately advanced data‐driven estimators [38, 39, 52, 53] can assess stability without knowledge of the power network parameters or topology. However, these estimators are based on the analysis of individual window‐truncated PMU data. In this chapter, we seek to provide a method with the ability of continuous learning of power system for spatiotemporal PMU data.
Firstly, we provide a novel method for modeling spatiotemporal PMU data. Figure 23.5 illustrates the conceptual representation of the structure of the spatiotemporal PMU data. More specifically, let denote the number of the available PMUs across the whole power network, each providing measurements. At the th time sample, a total of measurements, say , are collected. With respect to each PMU, the measurements could contain many categories of variables, such as voltage magnitude, power flow and frequency, etc. In this chapter, we develop PMU data analysis assuming each type of measurements is independent. That is, we assume that at each round of analysis, . Given time periods of seconds with Hz sampling frequency in the th data collection. Let and , a sequence of large random matrix
is obtained to represent the collected voltage PMU measurements.
As illustrated in Figure 23.5, is a large random matrix with independent identical distributed entries. Here we also include other forms of basic random matrices that are relevant to the applications in smart grids as follows.
Gaussian unitary ensemble (GUE): Let be Wigner matrix, also known as Gaussian unitary ensemble GUE, and . satisfies the following conditions:
For the convenience of analysis, we can denote GUE as . Besides, the joint of ordered eigenvalues of GUE is [54, 55]
Laguerre unitary ensemble (LUE): Let be Gaussian random variables with and . The so called Wishart matrix or Laguerre unitary ensemble LUE can be expressed as . The of for is [47, 55]
Large random matrix polynomials:
where are analytical functions, is a GUE, and is a LUE. See more details in Section 23.2.5.
In this subsection, we provide a brief introduction to the asymptotic spectrum laws of the large basic random matrices, as shown in Section 23.2.1. There are remarkable results describing the asymptotic spectrum laws. Here special attention is paid to the limit behavior of marginal eigenvalues as the data dimensions tend to infinity.
We start with the GUE matrix whose entries are independent identical distributed zero‐mean (real or complex) Gaussian ensembles. As shown in [47], as , the empirical distribution of eigenvalues of converges to the well‐known semicircle law, whose density can be represented as
Also shown in [40], the same result could be obtained for a symmetric whose diagonal entries are 0 and whose lower‐triangle entries are independent and take the values with equal probability.
If no attempt is made to symmetrize the square matrix , then the eigenvalues of are asymptotically uniformly distributed on the unit circle of the complex plane. This is referred to as the well‐known Griko circle law, which is elaborated in the following theorem.
The semicircle law and circular law explain the asymptotic property of large random matrices with independent entries. However, as illustrated in Section 23.1.2, the key issues in smart grids involve the singular values of rectangular large random matrices . The LUE matrices have dependent eigenvalues of interest even if has independent entries. Let the matrix aspect ratio ; the asymptotic theory of singular values of was presented by the landmark work [ 45] as follows.
As and , the limit distribution of the eigenvalues of converges to the so‐called Marcenko‐Pastur law, whose density function is
where and
Analogously, when , the limit distribution of the eigenvalues of converges to
In addition to Wigner's semicircle law above and Marchenko‐Pastur law, we are also interested in the single ring law developed by Guionnet, Krishnapur, and Zeitouni (2011) [ 46]. It describes the empirical distribution of the eigenvalues of a large generic matrix with prescribed singular values, i.e., an matrix of the form with some independent Haar‐distributed unitary matrices and a deterministic matrix whose singular values are the ones prescribed. More precisely, under some technical hypotheses, as the dimension tends to infinity, if the empirical distribution of the singular values of converges to a compactly supported limit measure on the real line, then the empirical eigenvalues distribution of converges to a limit measure on the complex plane that depends only on The limit measure is rotationally invariant in , and its support is the annulus with such that
The transforms of large random matrices are especially useful to study the limit spectral properties and to tackle the problems of polynomial calculation of random matrices. In this subsection, we will review the useful transforms including Stieltjes transform, R transform, and S transform suggested by problems of interest in power grids [49, 56].
We begin with the Stieltjes transform of , which is defined as follows.
An important application of the Stieltjes transform is its relationship with the limit spectrum density of .
It is noted that the signs of the and coincide. This property should be emphasized in the following examples where the sign of the square root should be chosen.
For GUE and LUE matrices, the corresponding Stieltjes transforms are shown in the following examples.
Another two important transforms, which we elaborate in the following, are the R transform and the S transform. The key point of these two transforms is that the R/S transform enables the characterization of the limiting spectrum of a sum/product of random matrices from their individual limiting spectra. These properties turn out to be extremely useful in the following subsection. We start with the blue function, that is, the functional inverse of the Stieltjes transform , which is defined as
and then the R transform is simply defined by
Two important properties of the R transform are shown in the following.
let , and be the R transforms of matrices , and , respectively. We have
For any ,
This additivity law can be easily understood in terms of Feynman diagrams; we refer interested readers to references [ 49] for details. The above properties of the R transform enable us to do the linear calculation of the asymptotic spectrum of random matrices.
Another important transform of engineering significance in RMT is the S transform. The S transform is related to the R transform and is defined by
An interesting property of the S transform is that the S transform of the product of two independent random matrices equals the product of the S transforms:
Note that (23.16) is known as multiplication law of the S transform. For the sake of brevity, see Section 23.2.5 for more details.
In this section, we investigate the spectral asymptotics for GUE and LUE matrices. We are motivated by the practical problems introduced in [ 49]. Let be the empirical spectral distribution function of GUE or LUE matrices and be the distribution function of the limit law (semicircle law for GUE matrices and Marchenko‐Pastur law for LUE matrices). Here, we study the convergence rate of the expected empirical distribution function to . Especially, the bound
is mainly concerned in the following.
The rate of convergence for the expected spectral distribution of GUE matrices has attracted much attention due to its increasingly appreciated importance in applied mathematics and statistical physics. Wigner initially looked into the convergence of the spectral distribution of GUE matrices [57]. Bai [58] conjectured that the optimal bound for in the GUE case should be of order . Bai and coauthors in [59] proved that . Gotze and Tikhomirov in [60] improved the result in [ 59] and proved that . Bai et al. in [61] also showed that on the condition that the 8th moment of satisfies . Girko in [62] stated as well that assuming uniform bounded 4th moment of . Recently, Gotze and Tikhomirov proved an optimal bound as follows.
The convergence of the density (denoted by ) of the standard semicircle law to the expected spectral density is proved by Gotze and Tikhomirov in the following theorem.
For LUE matrix with spectral distribution function , let as ; it is well known that convergences to the Marchenko‐Pastur law with density
where . The bound
for the convergence rate is shown in the following theorems.
Considering the case , a similar result is shown in Theorem 23.2.9.
Interested readers are referred to [63] for technical details and Section 23.3.3 for applications in the smart grid.
Free probability theory, initiated in 1983 by Voiculescu in [64], together with the results published in [65] regarding asymptotic freeness of random matrices, has established a new branch of theories and tools in random matrix theory. Here, we provide some of the basic principles and then examples to enhance the understanding and application of free probability theory.
Let be self‐adjoint elements which are freely independent. Consider a self‐adjoint polynomial in n non‐commuting variables and let be the element . Now we introduce the method [66, 67] to obtain the distribution of out of the distributions of .
Let be a unital algebra and be a subalgebra containing the unit. A linear map
is a conditional expectation if
and
An operator‐valued probability space consists of and a conditional expectation . Then, random variables are free with respect to (or free with amalgamation over ) if whenever are polynomials in some with coefficients from and for all and . For a random variable , we denote the operator‐valued Cauchy transform:
whenever is invertible in . In order to have some nice analytic behavior, we assume that both and are ‐algebras in the following; will usually be of the form , the ‐matrices. In such a setting and for , this is well defined and a nice analytic map on the operator‐valued upper half‐plane:
and it allows to give a nice description for the sum of two free self‐adjoint elements. In the following we will use the notation
Let be a complex and unital ‐algebra and let self‐adjoint elements . is given. Then, for any non‐commutative polynomial , we get an operator by evaluating at . In this situation, knowing a linearization trick [68] means to have a procedure that leads finally to an operator
for some matrices of dimension , such that is invertible in if and only if is invertible in . Hereby, we put
Let be given. A matrix
where
is called a linearization of , if the following conditions are satisfied:
To introduce the following corollary, which will enable us to shift for to a point
lying inside the domain in order to get access to all analytic tools that are available there.
Let be a non‐commutative ‐probability space, self‐adjiont elements which are freely independent, and a self‐adjoint polynomial in non‐commuting variables . We put . The following procedure leads to the distribution of .
by using the fixed point iteration for the operator‐valued free additive convolution. The Cauchy transform of is then given by
Finally, we obtain the desired distribution of by applying the Stieltjes inversion formula.
For readers' convenience, experimental results obtained in various conditions are also presented. Specially, Figures 23.6, 23.7 and 23.8 illustrate the theoretical limit spectra and empirical one in the case of the polynomial of random matrices introduced in Example . Figures 23.9, 23.10 and 23.11 present simulation results for the Example . We see that the theoretical results agree remarkably with the numerical simulations in various conditions.
In this section, we elaborate some of the more representative problems described in Section 23.1 that capture various features of interest in smart grids and we show how random matrix results have been used to tackle the problems that arise in the large power grid with wide deployment of PMU equipments. Besides, we also conclude with some state‐of‐art data driven methods for comparison.
Considering the data model introduced in Section 23.2.1, the problem of testing hypotheses on means of populations and covariance matrices is addressed. We start with a review of traditional multivariate procedures for these tests. Then we develop adjustments of these procedures to handle high‐dimensional data in smart grids.
As depicted in Section 23.2.1, a large random matrix flow is adopted to represent the massive streaming PMU data in one sample period. Instead of analyzing the raw individual window‐truncated PMU data [38, 39] or the statistic of [52, 53], a comprehensive analysis of the statistic of is conducted in the following. More specifically, denote as the covariance matrix of th collected PMU measurements; we want to test the hypothesis:
It is worthy noting that the hypothesis (23.24) is a famous testing hypothesis in multivariate statistical analysis that aims to study samples share or approximately share some distribution and consider using a set of samples (data streams denoted in equation (23.4) in this paper), one from each population, to test the hypothesis that the covariance matrices of these populations are equal.
The LR test [69] and CLR test [70] as introduced in Section 23.3.1 are most commonly test statistics for the hypothesis in ( 23.24). These tests can be understood by replacing the population covariance matrix by its sample covariance matrix . While direct substitution of by brings invariance and good testing properties as shown in [ 69] for normally distributed data. The test statistic may not work for high‐dimensional data as demonstrated in [71, 72]. Besides, the estimator has unnecessary terms, which slow down the convergence considerably when the dimension of PMU data is high [72, 73]. In such situations, to overcome the drawbacks, a trace criterion [ 72] is more suitable to the test problem. Specially, instead of estimating the population covariance matrix directly, a well‐defined distance measure exploiting the difference among data flow is conducted, that is, the trace‐based distance measure between and is
where is the trace operator. Instead of estimating , and by sample covariance matrix‐based estimators, we adopt the merits of the U‐statistics [74]. Especially, for ,
is proposed to estimate . It is noted that represents summation over mutually distinct indices. For example, says summation over the set . Similarly, the estimator for can be expressed as
The test statistic that measures the distance between and is
Then the proposed test statistic can be expressed as:
As , the asymptotic normality [ 73] of the test statistic (23.28) is presented in the following:
Let , the false alarm probability (FAP) for the proposed test statistic can be represented as
where . For a desired FAP , the associated threshold should be chosen such that
Otherwise, the detection rate (DR) can be denoted as
It is noted that the computation complexity of the proposed test statistic in (23.30) is which limits its practical application. Here, we proposed an effective approach to reduce the complexity of the proposed test statistic from to by principal component calculation and redundant computation elimination. For simplicity, we briefly explained the technical details in our recent work, which is available at .
In this section, we evaluate the efficacy of the proposed test statistic for power system stability. For the experiments shown in the following, the real power flow data were of a chain‐reaction fault happened in the China power grids in 2013. The PMU number, the sample rate, and the total sample time are , , and , respectively. The chain‐reaction fault happened from to . Let . Figure 23.12 shows that the mean and variance of agree well with theoretical ones. Based on the results in Figure 23.12 and event indicators (23.29), the occurrence time and the actual duration of the event can be identified as and , respectively. The location of the most sensitive bus can also be identified using the data analysis above. The result shown in Figure 23.13 illustrates that 17 and 18 PMU are the most sensitive PMUs, which is in accordance with the actual accident situation.
Situation awareness (SA) is of great significance in power system operation, and a reconsideration of SA is essential for future grids [ 24]. These future grids are always huge in size and complex in topology. Operating under a novel regulation, their management mode is much different from previous one.
All these driving forces demand a new prominence for the term “situation awareness” (SA). The SA is essential for power grid security; inadequate SA is identified as one of the root causes for the largest blackout in history—the August 14, 2003, blackout in the United States and Canada [75].
In [76], SA is defined as the perception of the elements in an environment, the comprehension of their meaning, and the projection of their status in the near future. This chapter is aimed at the use of model‐free and data‐driven methodology for the comprehension of the power grid.
The massive data compose the profile of the actual grid—present state; SA aims to translate the present state into perceived state for decision making [77].
The proposed methodology consists of three essential procedures as illustrated in Figure 23.14(b): (1) big data model—to model the system using experimental data for the RMM; (2) big data analysis—to conduct high‐dimensional analyses for the indicator system as the statistical solutions; (3) engineering interpretation—to visualize and interpret the statistical results for human beings for decision making.
Power grids operate in a balance situation obeying
where and are the power injections on node , while and are the injections of the network satisfying
For simplicity, combining (23.33) and (23.34), we obtain
where is the vector of power injections on nodes depending on , . is the system status variables depending on , , while is the network topology parameters depending on , .
For a system state with certain fluctuations, and thus randomness in data sets, we formulate the system as
With a Taylor expansion, (23.36) is rewitten as
Equ. ( 23.34) shows that is linear with ; it means that . On the other hand, the values of system status variables are relatively stable, and we can ignore the second‐order term and higher‐order terms. In this way, we turn (23.37) into
Suppose the network topology is unchanged, i.e., . From (23.38), we deduce that
On the other hand, suppose the power demand is unchanged, i.e., . From ( 23.38), we obtain that
where
Note that , i.e., the inversion of the Jacobian matrix , expressed as
Thus, we describe the power system operation using a random matrix—if there is an unexpected active power change or short circuit, the corresponding change of system status variables , i.e. , , will obey (23.39) or (23.40) respectively.
For a practical system, we can always build a relationship in the form of with a similar procedure as (23.35) to ( 23.40); it is linear in high dimensions. For an equilibrium operation system in which the reactive power is almost constant or changes much more slowly than the active one, the relationship model between voltage magnitude and active power is just like the multiple input multiple output (MIMO) model in wireless communication [49, 78]. Note that most variables of vector are random due to the ubiquitous noises, e.g., small random fluctuations in . In addition, we can add very small artificial fluctuations to make them random or replace the missing/bad data with random Gaussian variables. Furthermore, with the normalization, we can build the standard random matrix model (RMM) in the form of , where is a standard Gaussian random matrix.
The data‐driven approach conducts analysis requiring no prior knowledge of system topologies, unit operation/control mechanism, causal relationship, etc. It is able to handle massive data all at once; the large size of the data, indeed, enhances the robustness of the final decision against the bad data (errors, losses, or asynchronization). Compared with classical data‐driven methodologies (e.g., PCA), the RMT‐based counterpart has some unique characteristics:
We adopt a standard IEEE 118‐node system as the grid network (Figure 23.15) and the events is shown in Table 23.1.
Table 23.1 Series of Events.
is the power demand of node 52.
Stage | E1 | E2 | E3 | E4 |
Time (s) | ||||
(MW) | 0 |
The power demand of nodes are assigned as
where and are the element of a standard Gaussian random matrix; , . Thus, the power demand on each node is obtained as the system injections (Figure 23.16(a)); the voltage can also be obtained (Figure 23.16(b)). Suppose we sample the voltage data at 1 Hz, and the data source is denoted as . The number of dimensions is and the sampling time span is .
Suppose that the power demand data (Fig. 23.16(a)) are unknown or unqualified for SA due to low sampling frequency or bad quality. For further analysis, we just start with data source (Figure 23.16(b)) and assign the analysis matrix as (4 minutes' time span). First, we conduct category for the system operation status; the results are shown as Figure 23.16(c). In general, according to the raw data source and the analysis matrix size, we divide our system into 8 stages. Note that it is a statistical division— , and are transition stages, and their time span is right equal to the length of the analysis matrix minus one, i.e, . These stages are described as follows:
We also select two typical data cross‐sections for stage and : during period at the sampling time , and 2) during period at the sampling time .
Besides, as discussed in 23.2.1, we build up the RMM from the raw voltage data. Then, is employed as a statistical indicator to conduct anomaly detection. For the selected data cross section and , their M‐P law and ring law analysis are shown as Figure 23.17(a), 23.17(b), 23.17(c) and 23.17(d). With moving slide window (MSW), the curve is obtained as Figure 23.17(e).
Fig 23.17 shows that when there is no signal in the system, the experimental RMM well matches the ring law and M‐P law, and the experimental value of LES is approximately equal to the theoretical value. This validates the theoretical justification for modeling rapid fluctuation at each node with additive white Gaussian noise, as shown in Section 23.2.1. On the other hand, the ring law and M‐P law are violated at the very beginning ( ) of the step signal. Besides, the proposed high‐dimensional indicator , is extremely sensitive to the anomaly. At , the starts the dramatic change as shown in the curve as Figure 23.17(e), while the raw voltage magnitudes are still in the normal range as shown in Figure 23.16(c). Moreover, we design numerous kinds of LES and define The results are shown in Figure 23.18 and prove that different indicators have different characteristics and effectiveness; this suggests another topic to explore in the future.
Furthermore, we investigate the SA based on the high‐dimensional spectrum test. The sampling time is set as and . The following Lemma 23.2.7 and Lemma 23.2.9,
(span and ), and (span and ) are selected. The results are shown in Figure 23.19 and Fig. 23.20. These results validate that empirical spectral density test is competent to conduct anomaly detection—when the power grid is under a normal condition, the empirical spectral density and the ESD function are almost strictly bounded between the upper bound and the lower bound of their asymptotic limits. On the other hand, these results also validate that GUE and LUE are proper mathematical tools to model the power grid operation.
The curve (also called nose curve) and the smallest eigenvalue of the Jacobian matrix [ 39] are two clues for steady stability evaluation. In this case, we focus on the E4 part during which P Node‐52 keep increasing to break down the steady stability. The related curve and curve, respectively, are given in Figure 23.21(a) and Figure 23.21(b). Only using the data source , we choose some data cross section, as shown in Figure 23.21(a). The RMT‐based results are shown as Figure 23.22. The outliers become more evident as the stability degree decreases. The statistics of the outliers are similar to the smallest eigenvalue of the Jacobian matrix, Lyapunov exponent, or the entropy in some sense.
For further analysis, we take the signal and stage division into account. Generally speaking, sorted by the stability degree, the stages are ordered as . According to Figure 23.18, we make Table 23.2. The high‐dimensional indicators and have the same trend as the stability degree order. These statistics have the potential for data‐driven stability evaluation.
Table 23.2 Indicator of Various LESs at Each Stage.
MSR | T2 | T3 | T4 | DET | LRF | |
: Theoretical Value | ||||||
0.8645 | 1338.3 | 10069 | 8.35E4 | 48.322 | 73.678 | |
— | 665.26 | 93468 | 1.30E7 | 1.3532 | 1.4210 | |
[0240:0500, 261]: Small fluctuations around 0 MW | ||||||
0.995 | 1.010 | 1.040 | 1.080 | 0.959 | 1.014 | |
6E –6 | 78.38 | 3.03E4 | 7.14E6 | 0.4169 | 0.3908 | |
1 | 1 | 1 | 1 | 1 | 1 | |
[0501:0739, 239]: A step signal (0 MW 30 MW) is included | ||||||
0.9331 | 1.280 | 2.565 | 7.661 | 0.5453 | 1.284 | |
1.49E1 | 1.64E2 | 1.16E3 | 8.63E3 | 3.43E1 | 3.97E1 | |
[0740:0900, 161]: Small fluctuations around 30 MW | ||||||
0.9943 | 1.010 | 1.039 | 1.084 | 0.9568 | 1.015 | |
0.8608 | 0.9121 | 0.9476 | 1.234 | 0.8972 | 1.101 | |
[0901:1139, 239]: A step signal (30 MW 120 MW) is included | ||||||
0.8742 | 2.054 | 1.06E1 | 7.22E1 | 7E ‐2 | 1.597 | |
5.49E1 | 2.06E3 | 3.87E4 | 8.54E5 | 1.52E2 | 1.62E2 | |
[1140:1300, 161]: Small fluctuations around 120 MW | ||||||
0.9930 | 1.019 | 1.067 | 1.135 | 0.9488 | 1.021 | |
0.7823 | 1.053 | 1.189 | 1.135 | 0.7310 | 0.9255 | |
[1301:1539, 239]: A ramp signal (119.7 MW ) is included | ||||||
0.9337 | 1.295 | 2.787 | 9.615 | 0.5316 | 1.294 | |
8.50E1 | 7.41E2 | 5.63E3 | 5.17E4 | 2.14E2 | 2.30E2 | |
[1540:2253, 714]: Steady increase ( 358.1 MW) | ||||||
0.8906 | 1.717 | 6.530 | 3.48E1 | 0.1483 | 1.545 | |
1.35E1 | 3.28E2 | 5.33E3 | 1.10E5 | 6.11E1 | 6.85E1 | |
[2254:2500, 247]: Static voltage collapse (361.9 MW ) | ||||||
0.4259 | 1.02E1 | 2.11E2 | 4.65E3 | –1.4E1 | 1.08E1 | |
1.94E3 | 5.81E5 | 1.20E8 | 3.2E10 | 9.02E4 | 9.62E4 |
a) ; .
The key for correlation analysis is the concatenated matrix , which consists of two parts—the basic matrix and a certain factor matrix , i.e., . For more details, see our previous work [ 53]. The LES of each is computed in parallel, and Figure 23.23 shows the results.
In Figure 23.23, the blue dotted line (marked with “None”) shows the LES of basic matrix , and the orange line (marked with “Random”) shows the LES of the concatenated matrix ( is the standard Gaussian random matrix). Figure 23.23 demonstrates that: (1) node 52 is the causing factor of the anomaly; (2) sensitive nodes are 51, 53, and 58; and (3) nodes 11, 45, 46, etc., are not affected by the anomaly. Based on this algorithm, we can continue to conduct behavior analysis, e.g., detection and estimation of residential PV installations [84]. Behavior analysis is a big topic. Because of space limitations, we will not expand it here.
Following [85], we build the statistic model for power grids. Considering random vectors observed at time instants we form a random matrix as follows
In an equilibrium operating system, the voltage magnitude vector injections with entries and the phase angle vector injections with entries experience slight changes. Without dramatic topology changes, rich statistical empirical evidence indicates that the Jacobian matrix keeps nearly constant, and so does . Also, we can estimate the changes of and only with the classical approach. Thus, we rewrite (23.43) as:
where , and Here and are random matrices. In particular, is a random matrix with Gaussian entries.
Multivariate linear or nonlinear polynomials perform a significant role in problem modeling, so we build our models on the basis of random matrix polynomials. Here, we study two typical random matrix polynomial models.
The first case is the multivariate linear polynomial:
The second one is the self‐adjoint multivariate nonlinear polynomial:
Here, both and are the sample covariance matrices. The asymptotic eigenvalue distributions of and can be obtained via basic principles of free probability theory, as introduced above. The asymptotic eigenvalue distributions of are regarded as the theoretical bounds.
We formulate our problem of anomaly detection in terms of the same hypothesis testing as [ 85]: no outlier exists , and outlier exists .
where is the standard Gaussian random matrix.
Generate , from the sample data through the preprocess in 23.3.4. Compare the theoretical bound with the spectral distribution of raw data polynomials. If an outlier exists, will be rejected, i.e., signals exist in the system.
The data sampled from the power grid is always non‐Gaussian, so we adopt a normalization procedure in [ 80] to conduct data preprocessing. Meanwhile, we employ the Monte Carlo method to compute the spectral distribution of the raw data polynomial according to asymptotic property theory. See details in Algorithm 1.
Our data fusion method is tested with simulated data in the standard IEEE 118‐bus system. Detailed information of the system refers to the case118.m in Matpower package and Matpower 4.1 User's Manual [86]. For all cases, let the sample dimension . In our simulations, we set the sample length equal to , i.e., , and select six sample voltage matrices presented in Table 23.3, as shown in Figure 23.24. The results of our simulations are presented in Figure 23.25 and Figure 23.26. The outliers existed when the system was abnormal, and its sizes become large when the anomaly becomes serious.
Table 23.3 System status and sampling data.
Cross Section (s) | Sampling (s) | Descripiton |
Reference, no signal | ||
Existence of a step signal | ||
Steady load growth for Bus 22 | ||
Steady load growth for Bus 52 | ||
Chaos due to voltage collapse | ||
No signal |
a) We choose the temporal end edge of the sampling matrix as the marked time for the cross section. E.g., for , the temporal label is 217, which belongs to . Thus, this method is able to be applied to conduct real‐time analysis.
Motivated by the immediate demands of tackling the tricky problems raised from large‐scale smart grids, this chapter introduced RMT‐based schemes for spatiotemporal big data analysis. Firstly, we represent the spatiotemporal PMU data as a sequence of large random matrices. This is a crucial part for power state evaluation, as it turns the big PMU data into tiny data for practical uses. Rather than employing the raw PMU data, a comprehensive analysis of PMU data flow, namely, RMT‐based techniques, is then proposed to indicate the state evaluation. The core techniques include streaming PMU data modeling, asymptotic properties analysis, and data fusion methods (based on free probability). Besides, the case studies based on synthetic data and real data are also included with the aim to bridge the technology gap between RMT and spatiotemporal data analysis in smart grids.
The current work based on RMT provides a fundamental exploration of data analysis for spatiotemporal PMU data. Much more attention is to be paid to this research direction, such as classification of power events and load forecasting. It is also noted that this work provides data‐driven methods that are new substitutes for power system state estimation. The combination of power system scenario analysis, spectrum sensing mechanisms, networking protocols, and big data techniques [15, 34, 49, 87] is encouraged to be investigated for better understanding of the power system state.
3.129.87.138