Parameter Identification of Power Grid Subsynchronous Oscillations Based on Eigensystem Realization Algorithm

: The subsynchronous oscillation caused by the resonance between power electronic devices and series compensation devices or weak power grids introduced by large-scale renewable energy generation greatly reduces the transmission capacity of the system and may endanger the safe operation of the power system. It even leads to system oscillation instability. In this paper, based on the advantages of a simple solution, a small amount of calculation and anti-noise of ERA, a method of subsynchronous oscillation parameter identification based on the eigensystem realization algorithm (ERA) is proposed. The Hankel matrix in the improved ERA is obtained by splicing the real part matrix and the imaginary part matrix of the synchrophasor, thus solving the problem of angular frequency conjugate constraints of two fundamental components and two oscillatory components which are not considered in the existing ERA. The solution to this problem is helpful to improve the accurate parameter identification results of ERA under the data window of 200 ms and weaken the limitation caused by the assumption that the synchrophasor model is fixed. The practicability of the improved method based on PMU is verified by the synthesis of ERA and the actual measurement data. Compared with the existing ERA, the improved ERA can accurately identify the parameters of each component under the ultra-short data window and realize the dynamic monitoring of power system subsynchronous oscillation.


Introduction
Resonance between power electronic devices and series compensators or weak grids frequently causes subsynchronous oscillations (SSOs) in modern power systems with renewable generations.SSO occurs when voltage and current in a power system have a frequency lower than synchronization frequency.In some cases, there may also be a supersynchronous component with a frequency higher than the synchronization frequency.SSOs can reduce transmission capacity, threaten power system safety, and cause instability.Therefore, it is essential to closely monitor the dynamic process of SSOs in the power system and implement control measures in time to restrain SSOs effectively [1][2][3][4][5].
There are several methods for detecting subsynchronous oscillatory dynamic processes.In the literature [6], digital fault recorders are utilized for online monitoring of subsynchronous oscillations.Electrical and mechanical signals are both effective for SSO monitoring.This expands our monitoring capabilities and helps us identify potential issues proactively, saving costs and time.In references [7][8][9], Prony analysis, a model-based parameter identification method, is introduced.This method has high accuracy when the signal-noise ratio is high.However, the algorithm contains a large amount of computation and takes a long time.Moreover, all the methods mentioned above require the use of instantaneous data from the fault recorders (FRs).Collecting FR data on time is challenging due to the installation of FRs on various buses.Additionally, observing SSO globally is challenging due to the lack of a unified timing standard for FRs.
At present, the Wide Area Measurement System (WAMS) and Phasor Measurement Unit (PMU) are widely used in power systems.The most significant advantage of subsynchronous oscillation analysis based on synchrophasor data is that the data measured by different power stations have a unified time scale and comparability.SSO parameter identification methods based on synchrophasors can be classified as discrete Fourier transform (DFT) algorithm-based or not.
In DFT-based algorithms, literature [10] utilizes WAMS/PMU to analyze the impact of the DFT algorithm and a synchrophasor sampling rate on the amplitude spectrum.It provides a correction ratio for the accurate identification of SSO parameters.However, to accurately extract sub-synchronous oscillation component parameters, this method requires a large data window of 10 s.This is less reliable and real-time in identifying rapidly changing SSOs.Therefore, some scholars have proposed an interpolation DFT method combined with a Hann window [11,12].This approach not only reduces the data window to 2 s, but also effectively eliminates spectrum aliasing influence.Reference [13] is built upon by literature [11] through an interpolation method utilizing the damped Rife-Vincent window function to enhance SSO identification accuracy.This method shortens the data window to 2 s.
Among non-DFT algorithms, there are mainly two modal parameter extraction algorithms based on synchrophasors [14], namely the matrix beam method [15][16][17] and the feature system implementation algorithm [18].The non-DFT algorithms are less limited by the resolution of the spectrum.Therefore, a time window of 2 s can be used instead.Although the method above can ensure the accuracy of identifying sub-synchronous oscillation parameters to a certain extent, 2 s is still a relatively long time window.In addition, the above method does not consider the constraint of angular frequency conjugation, which reduces the accuracy of identifying the SSO parameters.
The eigensystem realization algorithm (ERA) is a classical modal parameter extraction algorithm.To extract the modal parameters, we can construct the Hankel matrix and system matrix.This transformation will convert the problem into a matrix decomposition problem.By doing so, we can obtain the angular frequency and coefficient of each mode directly.This is a simple and low-computation approach to solving the problem.In theory, there is an angular frequency conjugate constraint between the two modes of fundamental components.There is also an angular frequency conjugate constraint between the two modes of the oscillatory component.However, the ERA method cannot consider the angular frequency conjugate constraints as synchrophasors are in the complex number field.Therefore, the identification error of subsynchronous oscillation parameters is significant in the actual situation.
During the research process, we conducted an analysis based on fixed models, but in actual situations, these models are often affected by various factors and changes.A longer data window will be more limited to the assumption of fixed models, because as the data window becomes longer, the data within it may exhibit more complexity as time changes.The non-conjugation of angular frequency reduces the accuracy of subsequent parameter identification.Therefore, based on practical applications, it is necessary to shorten the data window.However, it was found in research that the identification accuracy of ERA is greatly reduced after shortening the data window.Therefore, we improved the algorithm to maintain high accuracy even in a shorter data window, which is our fundamental starting point.
This paper improves ERA's Hankel matrix construction, solving the issue of obtaining angular frequencies without conjugations in the existing ERA.The Hankel matrix is constructed by combining the real and imaginary parts of the synchrophasor, and it can retain all the information of synchrophasors.The accuracy of identifying fundamental wave parameters and subsynchronous and supersynchronous oscillation parameters is significantly improved.At the same time, the improved ERA shortens the data window to 200 ms, which can better mitigate the limitations of assuming the model to be constant and is more conducive to real-time monitoring of the dynamic process of SSO in power systems.Finally, the improved ERA methodology is presented and validated using synthetic, simulated, and measured PMU data.It is then compared with the existing ERA.The ERA can accurately identify the frequency, amplitude, and phase of the fundamental, subsynchronous, and supersynchronous components of a single synchronization signal using a 200 ms data window.
The structure of this article is as follows.In the second part, the synchrophasors model and the extraction of modal parameters are studied.The third section introduces SSO parameter identification based on improved ERA.In section four, the performance of the proposed method is evaluated using synthetic, simulated, and measured PMU data.

SSO Model with Supersynchronous Component
The power system's instantaneous signal model, denoted as x(t), is established.This model consists of a fundamental component with a possible frequency deviation, a pair of frequency-coupled subsynchronous components, and supersynchronous components.Formula (1) demonstrates these components.
The synchrophasor sequence Ẋ(k) corresponding to the instantaneous signal x (t) can be obtained by DFT of Formula (1), as shown in (2).
where Ẋ0 (k), Ẋsub (k), and Ẋsup (k) are the fundamental components, subsynchronous components, and supersynchronous components of the synchrophasor, respectively.The subsynchronous component and the supersynchronous component are a pair of coupling components.Their frequencies satisfy the relationship of f sub + f sup = 2 f N , where f N is the nominal frequency.
According to the theory of DFT, each component contains a positive frequency component and a negative frequency component [12,19].The subsynchronous component and the super-synchronous component are a pair of coupling components, so the oscillatory component ( Ẋ+ s (k), Ẋ− s (k)) can be used to replace the sum of the subsynchronous component and the supersynchronous component, as shown in (3) and (4).
where "+, −" denotes positive and negative frequencies and "*" indicates conjugation; besides, "α" and "β" are the angular frequencies of fundamental and oscillatory components, respectively, and are expressed in the following formula together with Q( f , l): where "N" is the number of instantaneous data in DFT synchrophasors calculation, and " f s " is the synchrophasor data reporting frequency, and The modal model of the synchrophasor under subsynchronous oscillation is expressed as ( 6)-( 8) x sub e jϕ sub + Q * ( f sup , +1)x sup e −jϕ sup ( 7) where "R 1 , R 2 , R 3 , R 4 " are four constant parameters, "ω 1 , ω 2 " are the positive and negative angular frequencies of the fundamental components, and "ω 3 , ω 4 " are the positive and negative angular frequencies of the oscillating components.Note that in (6), only "k" is a variable that represents the characteristics of the angular frequency over time.

Modal Parameter Extraction Problem
The extraction of frequency, amplitude, and phase for each component of subsynchronous oscillation from a vast amount of synchrophasor data is a typical problem in modal parameter extraction.ERA is a classical method used for this purpose.The ERA uses the Hankel matrix and the system matrix to transform the problem of extracting modal parameters into a matrix decomposition problem.This allows for the direct determination of the angular frequencies (ω 1 , ω 2 , ω 3 , ω 4 ) and coefficients (R 1 , R 2 , R 3 , R 4 ) for each component.This approach offers simplicity and computational efficiency.
The central idea of ERA is to use the time continuity of synchrophasor sequence Ẋ(k) to construct two Hankel matrices, and the two elements of two Hankel matrices in the same position have the same time difference.Through the relationship between the two Hankel matrices, a system characteristic matrix can be obtained, that is, diag[e ω 1 , e ω 2 , e ω 3 , e ω 4 ].Then, ( 8) is used to obtain coefficients (α, β), and ( 6) is used to obtain coefficients (R 1 , R 2 , R 3 , R 4 ).Finally, phases (ϕ 0 , ϕ sub , ϕ sup ) of the fundamental component, the subsynchronous component, and the supersynchronous component can be further obtained based on R m , (5), and (7).

Parameter Identification of SSO Based on
where K is the number of synchrophasors and the Hankel matrixes are matrices of (K − L) rows and (L + 1) columns.
The Y matrix can be decomposed in the following form: Each synchronphasor piece of data in the Y matrix is replaced with the synchronphasor data of the next measurement, and the new matrix is marked Y ′ .Then, the relationship between Y ′ and Y is shown in (13).
where Z 0 = diag[e w 1 , e w 2 , e w 3 , e w 4 ].Each row and column in Y is characterized by temporal continuity; the same goes for Y ′ .Therefore, the values of four angular frequencies can be obtained only by solving the Z 0 matrix.According to Formula (13), the system matrix Z 0 can be solved by the following formula: From matrix Z 0 , "ω 1 , ω 2 , ω 3 , ω 4 " can be obtained, and R m can be further calculated based on ω m and (6).

Defects of Existing ERA
The existing ERA directly constructs the Hankel matrix based on the synchrophasor in the complex number field and then constructs the complex number field system matrix.By solving the eigenvalues of the system matrix, we can obtain ω 1 , ω 2 , ω 3 , ω 4 , and R 1 , R 2 , R 3 , R 4 are solved by the least square method.This method can be solved accurately under ideal conditions.
However, in practical applications, when directly solving the eigenvalues of the system matrix Z 0 , the obtained ω 1 and ω 2 do not satisfy the conjugate relations of jα and −jα.ω 3 and ω 4 also do not satisfy the conjugate relation of jβ and −jβ.Thus, it is impossible to solve α and β accurately.As a result, the parameter identification error of each component is significant.
The complex number field matrix constructed by ERA makes it challenging to guarantee pairwise conjugates of the four eigenvalues obtained by directly solving the system matrix.If the Hankel matrix can be constructed in the real number domain, the angular frequency conjugate constraint problem can be solved in theory.

The Idea of Improved ERA
The basic idea of the improved ERA is to construct the real number domain Hankel matrix, construct the system matrix based on the real number domain Hankel matrix, and solve the eigenvalues of the system matrix.The eigenvalues obtained in this way can satisfy the angular frequency conjugate constraints of the positive and negative spectrums of the fundamental component and the positive and negative spectrums of the oscillatory component.That is, e ω 1 and e ω 2 satisfy the conjugate relations of e jα and e −jα ; e ω 3 and e ω 4 also satisfy the conjugate relation of e jβ and e −jβ .This solves the problem of the existing ERA directly calculating the eigenvalues of the complex number field system matrix, which leads to a large error.With the solution of this problem, the improved ERA can still maintain high identification accuracy under the ultra-short data window of 200 ms, which overcomes the limitation of the fixed hypothetical model to some extent.
Next, this section introduces the construction of the Hankel matrix in the real number field and the solution of modal parameters.

Construction of Real Number Field Hankel Matrix
Part 1: We use the complex number field Hankel matrices Y and Y ′ in Section 3.1 to construct two real number field Hankel matrices H and H ′ by extracting the real and imaginary parts of each element in the matrix.The detailed process of constructing the realvalued Hankel matrix is as follows.First, we express the real part matrix and the imaginary part matrix of the synchronized phasor sequence, as shown in Formulas ( 15) and ( 16), respectively.By observing these two matrices, we can see that each element in the matrices is a real number, and the real part matrix and the imaginary part matrix are the basis for constructing the real-valued Hankel matrix.
where Re represents the real part and Im represents the imaginary part.Second, the Hankel matrix constructed in the improved ERA method is a real-valued matrix, which is obtained by directly concatenating the real part matrix and the imaginary part matrix.This process is quite simple.Therefore, the real-valued Hankel matrix in the improved ERA can be expressed in the following formula: Similarly, the expression for another real-valued shifted Hankel matrix can be represented by Formula (18).
where H and H ′ are matrices with 2(K − L) rows and (L + 1) columns.
On the whole, the real number field matrix is composed of the real part matrix and the imaginary part matrix of the complex number field Hankel matrix.Compared to using only the real part matrix or the imaginary part matrix, its advantage is that all the information of the synchrophasor sequence can be fully retained, and matrices H and H ′ are very easy to obtain.Why the real part matrix and the imaginary part matrix can be spliced to construct a new real number field Hankel matrix is the focus of the following sections.
Part 2: We prove why the real part matrix and the imaginary part matrix can be spliced to obtain the real number field Hankel matrix, which is the core of this article.The real part matrix and the imaginary part matrix are composed of the real part and the imaginary part of the synchronous phasor sequence, respectively.Therefore, it is necessary to analyze the composition of the real and imaginary parts of the synchronous phasor.The analysis process is as follows.
The real and imaginary parts of fundamental and oscillatory components can be expressed using (3), (4), and (7), and their expressions are as follows: The above formula is simplified as follows: we let Then, according to (23), we can find Therefore, combining Formulas ( 3) and ( 19)-( 24), the real part and the imaginary part of the synchrophasor can be expressed as follows: Through ( 6), (25), and (26), we can determine that both the real part and the imaginary part of the synchrophasor sequence contain four modes.That is four modes with angular frequencies ±jα and ±jβ, which are consistent with the complex synchrophasor sequence modes.
The difference between the real and imaginary parts of a synchronized phasor after being split and the complete synchronized phasor lies only in the different constant parameters preceding the complex exponential, while the modal patterns remain the same.Therefore, the real part matrix and the imaginary part matrix can be spliced to construct a new real number field Hankel matrix.
Part 3: the real number field Hankel matrix H is decomposed by singular value decomposition, imitating "Y = Z 1 RZ 2 " and written in the form of "H = URV".The specific process is as follows.
The coefficients of the positive and negative components of fundamental components of Re[ Ẋ(k)] and Im[ Ẋ(k)] are not equal, but the coefficients of Ẋ(k) can be expressed in a particular proportion.Therefore, we let n (n = 1, 2) are four constants.Using Formulas (7), ( 23) and ( 27), the numerical details can be obtained by calculation as follows: Similarly, the expressions of γ The coefficients of the positive and negative components of the oscillation components of Re[ Ẋ(k)] and Im[ Ẋ(k)] are also not equal.They are also proportional to the coefficient of Ẋ(k), as shown in (30).
n (n = 3, 4) are four constants.Using Formulas ( 7), (23), and (30), the numerical details can be obtained by calculation as follows and the detailed derivation process is in Appendix A.
x sup e jϕ sup (31) Similarly, the expressions of γ ( 4 are shown in (32).Besides, the detailed derivation process is in Appendix A. where matrix U is a matrix with 2(K − L) rows and four columns, and matrix V is a matrix with four rows and (L + 1) columns.
According to the four Formulas, (10), ( 13), ( 15) and (33), we can determine that the expression of the Hankel matrix H ′ is as follows: where A = diag[e ω 1 , e ω 2 , e ω 3 , e ω 4 ] is the system matrix.So far, we have constructed the Hankel matrix of the real number field.In fact, although it is proved that it is a difficult process to concatenate the real part matrix and the imaginary part matrix, the constructed real field matrix is easy to obtain.This further illustrates the feasibility of this method.

Solution of Modal Parameters
In the beginning, we need to solve the system matrix A to obtain the angular frequencies of various components.Using (34), the equation for solving A is shown in (35).
After we obtain the matrix of system A = diag[e ω 1 , e ω 2 , e ω 3 , e ω 4 ], we can obtain ω 1 , ω 2 , ω 3 , ω 4 .The obtained ω 1 and ω 2 satisfy the conjugate relation of jα and −jα; ω 3 and ω 4 also satisfy the conjugate relation of jβ and −jβ, so angular frequency α of the fundamental component and angular frequency β of the oscillating component can be accurately solved.Now that ω 1 , ω 2 , ω 3 , and ω 4 are known, α and β can be obtained through Formula ( 8).The detailed process of solving the frequency, amplitude, and phase of each component is as follows: First, four constants R 1 , R 2 , R 3 , and R 4 are solved according to (36).
Second, according to Formulas ( 4) and ( 7), the angular frequency of the fundamental component, the angular frequency of the oscillating component, and the four constant parameters in the modal model, the calculated values of fundamental positive frequency component Ẋ+ 0 (k), fundamental negative frequency component Ẋ− 0 (k), oscillatory positive frequency component Ẋ+ s (k), and oscillatory negative frequency component Ẋ− s (k) are determined as in (37).
Then, combined with Formulas (4) and ( 5), the frequency, amplitude, and phase of the fundamental component are determined as shown in (38).
Finally, combined with Formulas (4) and ( 5) and f sub + f sup = 2 f N , the frequency, amplitude, and phase of the subsynchronous component and the frequency, amplitude, and phase of the supersynchronous component are determined according to the method of (39).
So far, the frequency, amplitude, and phase of the fundamental component, subsynchronous component, and supersynchronous component of the synchrophasor have been obtained.The modal parameters of the synchrophasor are finished.
The conclusion is that this algorithm has two outstanding advantages.On the one hand, by constructing the real-domain Hankel matrix, the problem of the angular frequency obtained by the existing ERA directly calculating the eigenvalues of the complex number field system matrix not satisfying the ±jα and ±jβ conjugate constraints is solved.On the other hand, the Hankel matrix is constructed by combining the real part and the imaginary part of the synchrophasor, which can completely retain the information of the synchrophasor and effectively improve the identification accuracy of fundamental parameters, subsynchronous oscillation parameters, and supersynchronous oscillation parameters.

Verification
To verify the correctness, feasibility, and practicality of the proposed power grid subsynchronous oscillation parameter identification method based on the eigensystem realization algorithm, the algorithm was validated using both synthetic and actual measured PMU data in MATLAB R2021a software.
The decision to use synthetic PMU data as the testing foundation was primarily based on practical considerations.While simulation scenarios and actual measured PMU data are of great value for power system dynamic analysis, they often cannot comprehensively cover all possible frequency ranges and oscillation characteristics.The operating states and fault conditions of actual power grids are complex and variable, making it difficult to guarantee that all types of subsynchronous oscillation data can be captured in real-world environments.To thoroughly test the performance and reliability of the algorithm, we need a dataset that allows for us to control the characteristics of the input data.Therefore, we chose synthetic oscillation data which can precisely set oscillation parameters according to needs, simulating various complex subsynchronous oscillation scenarios.By using these synthetic data, we can more systematically evaluate the performance of the algorithm under different conditions, ensuring its correctness and effectiveness in practical applications.
The choice to use actual measured PMU data is because they directly reflect the dynamic characteristics and oscillation features of the power system in actual operation.By applying these data to the algorithm presented in this paper, we can more accurately assess the feasibility of the algorithm in real-world scenarios.Additionally, using actual measured PMU data also helps us gain a deeper understanding of the power system's dynamic characteristics.Through the analysis of real data, we can discover phenomena that are difficult to observe in simulations or models, thus gaining a better understanding of the operation patterns and potential risks of the power system.
In the process of verification, the sampling frequency of instantaneous data is 2.0 kHz.First, the sampled instantaneous data are operated by DFT, and the synthetic PMU data are obtained.The frequency of reporting PMU data is 100 Hz.In the SSO parameter identification of ERA, a data window of 200 ms, that is, the Hankel matrix composed of 21 synchrophasor data, is selected.
Previous research [12] indicated that identification accuracy is most affected by f 0 and f sub .Therefore, by using the method of controlling variables, the SSO parameters of the synthesized PMU data are identified by changing f sub in the range of [5,45] Hz and x sub in the range of [5,100].At the same time, to evaluate the accuracy of the improved algorithm in identifying parameters under the influence of noise, zero-mean white noise is added to the instantaneous signal (1).According to literature [11], the measured signal-to-noise ratio (SNR) of the PMU is usually around 45 dB, so 40 dB is selected as the noise condition for the simulation experiment.
The performance of the improved ERA (referred to as "IERA") is compared with that of the existing ERA.The test results are listed in Table 1, and the results are analyzed as follows below.
namely "ERA 1 s", "ERA 200 ms", and "IERA 200 ms", is estimated by repeated identification of 20,000 times under ideal and noise conditions.The computer processor used is a 13th Gen Intel(R) Core(TM)i7-13620H.The identification times of the three algorithms are shown in Table 2.The identification times of the three algorithms are analyzed and compared as follows below.
From Table 2, it can be seen that the calculation time of the improved ERA is close to or even slightly smaller than that of the existing ERA 200 ms, and it is shorter than that of the existing ERA of 1 s.Therefore, the algorithm in this paper does not increase the parameter identification time due to the separation of the real and imaginary parts of the synchrophasors, which further demonstrates the feasibility and practicality of the algorithm in this paper.To further verify the effectiveness of the improved ERA, the simulated PMU data and the actual measured PMU data are analyzed.The simulation PMU data are generated based on Formula (1).The actual measured PMU data are from [12], which are recorded in an SSO event in the North China Power Grid.The current instantaneous data and amplitudes of synchrophasors of actual measured PMU data are shown in Figure 1.As shown in Figure 1, a 10 s data window is selected as the simulation object for the actual measured PMU data, and the SSO is changing rapidly.In this paper, the data window of the improved ERA is 200 ms.In addition, the existing ERA is used for comparison, and the lengths of the data window are 1 s and 200 ms, respectively, and there is the interpolation DFT algorithm with a data window length of 2 s for reference and comparison.Since there is little difference in the identification results of the fundamental components obtained by the three methods, only the identification results of the oscillatory components are introduced in detail.In addition, the subsynchronous component and the supersynchronous component are a pair of coupling components, and the identification result of the supersynchronous component is similar to that of the subsynchronous component, so the identification result of the subsynchronous component is analyzed as shown in Figure 2.  2 shows the identification results of subsynchronous oscillation components of the three methods when the length of the data window selected by the existing ERA is 1 s.It can be seen from the figure that there is little difference in the identification results of the three methods.However, the data window lengths of existing ERA and interpolation DFT are 1 s and 2 s, respectively.Theoretically, the longer the data window, the more information can be extracted, and the more accurate the identification result, the worse the real-time performance.The improved ERA only needs the data window of 200 ms to obtain the identification results with little difference in accuracy, but at the same time, the improved 200 ms ultrashort data window of ERA realizes the dynamic monitoring of power system SSO.  2 shows the identification results of subsynchronous oscillation components of the three methods when the length of the data window selected by the existing ERA is 200 ms.It can be seen from the figure that the f sub identification results of the subsynchronous oscillation components of the existing ERA are quite different from those of the interpolated DFT and the improved ERA in the 200 ms data window.The f sub identification results of the improved ERA are concentrated between 8.18 Hz and 8.40 Hz, while among the existing ERA f sub identification results, the smallest is 8.11 Hz, and the largest is even close to 8.60 Hz.It can be seen from the figure that the f sub identification result of the existing ERA exceeds that of the 8.50 Hz by five times, which seriously deviates from the identification result of the interpolated DFT.There is little difference in x sub identification results of the three methods.Moreover, the data window of 200 ms is 10 times shorter than that of 2 s, so the improved ERA can also reflect the real-time characteristics of SSO.

Conclusions
Based on the advantages of a simple solution, a small amount of calculation and anti-noise of ERA, an improved SSO parameter identification method based on ERA is proposed, and the dynamic real-time monitoring of SSO is realized.The existing ERA does not consider the angular frequency conjugate constraints of two fundamental components and two oscillatory components, but the improved ERA solves this problem and makes the results of parameter identification more accurate.At the same time, the Hankel matrix needed by the improved ERA can be easily obtained, and only the real part matrix and the imaginary part matrix of the synchrophasor need to be spliced.The improved method based on PMU is verified by the synthesis of ERA and the actual measurement data.Compared with the existing ERA, the improved ERA can identify parameters more accurately under the data window of 200 ms, and the ultra-short data window of 200 ms also realizes the dynamic monitoring of power system SSO.
In addition, there is also a certain limitation in improving ERA.This limitation is a common challenge faced by all identification methods, which is the limited ability to process signals that are too divergent.This will be a key area for improvement in future research work.However, the algorithm proposed in this article minimizes the impact of signal divergence on parameter identification results by shortening the data window.The improved ERA selects a 200 ms data window, which can effectively solve the most rapidly changing oscillation problems in practical applications.

Data Availability Statement:
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors..

Conflicts of Interest:
The authors declare no conflicts of interest.

Figure 1 .
Figure 1.The instantaneous signal of current and the amplitudes of the corresponding synchrophasor sequence for the actual measured PMU data case.(a) Instantaneous data.(b) Amplitude of synchrophasor.

4. 2 . 2 .
Figure (b) of Figure2shows the identification results of subsynchronous oscillation components of the three methods when the length of the data window selected by the existing ERA is 200 ms.It can be seen from the figure that the f sub identification results of the subsynchronous oscillation components of the existing ERA are quite different from those of the interpolated DFT and the improved ERA in the 200 ms data window.The f sub identification results of the improved ERA are concentrated between 8.18 Hz and 8.40 Hz, while among the existing ERA f sub identification results, the smallest is 8.11 Hz, and the largest is even close to 8.60 Hz.It can be seen from the figure that the f sub identification result of the existing ERA exceeds that of the 8.50 Hz by five times, which seriously deviates from the identification result of the interpolated DFT.There is little difference in x sub identification results of the three methods.Moreover, the data window of 200 ms is 10 times shorter than that of 2 s, so the improved ERA can also reflect the real-time characteristics of SSO.

Table 2 .
The average single identification time of the three algorithms (ms).