Dynamic State Estimation for Synchronous Machines Based on Adaptive Ensemble Square Root Kalman Filter

Dynamic state estimation (DSE) for generators plays an important role in power system monitoring and control. Phasor measurement unit (PMU) has been widely utilized in DSE since it can acquire real-time synchronous data with high sampling frequency. However, random noise is unavoidable in PMU data, which cannot be directly used as the reference data for power grid dispatching and control. Therefore, the data measured by PMU need to be processed. In this paper, an adaptive ensemble square root Kalman filter (AEnSRF) is proposed, in which the ensemble square root filter (EnSRF) and Sage–Husa algorithm are utilized to estimate measurement noise online. Simulation results obtained by applying the proposed method show that the estimation accuracy of AEnSRF is better than that of ensemble Kalman filter (EnKF), and AEnSRF can track the measurement noise when the measurement noise changes.


Introduction
In order to obtain the optimal control strategy of generator in power system, phasor measurement unit (PMU) is required to observe the generator data (power angle, electric angular velocity, etc.) to determine its state.However, there are random noises and measurement errors in PMU measurements, which cannot be directly used as the true state value [1][2][3].Statistical method is used to calculate the estimated value of generator state variables and predict the future state.This method is called dynamic state estimation (DSE).Therefore, the prediction of generator's future state by DSE is helpful for power system regulation and has important significance to the stability of power network.
In recent years, the problem of electromechanical transient state estimation of generators has attracted wide attentions.Generator is a complex nonlinear system, and its DSE needs corresponding nonlinear filtering algorithm.In the field of generator DSE algorithms, many studies mainly focus on extended Kalman filter (EKF), unscented Kalman filter (UKF) and cubature Kalman filter (CKF) [4,5].In [6], EKF was implemented to estimate the dynamic state of generator, where the fourth-order model of generator with unknown input was utilized.However, the linearization process of the generator model by EKF will lead to large truncation errors, which cannot be applied to the strongly nonlinear systems.In [7,8], the deterministic sampling filter-UKF algorithm was applied to estimate the dynamic state of generator.Due to the poor flexibility of UKF in parameter selection, the CKF method was proposed by S. Haykin [9,10].In [11], the CKF was used for the DSE of generator, which could achieve much better performance than both EKF and UKF.
On the other hand, considering the system uncertainties, in [12], based on the second-order model of generator, adaptive interpolation and H ∞ theory was introduced into EKF to estimate the dynamic state of generator, by which the robustness was improved.In [13], an online state estimation method for synchronous generators based on UKF was proposed, and the fourth-order nonlinear model of generators was used to verify that the algorithm has not only high filtering speed and accuracy but also strong robustness.In [14], square root filtering was introduced into CKF to ensure the nonnegative qualitative and numerical stability of covariance matrix, by which, the asymmetric or nonnegative covariance of CKF in the iteration of DSE could be avoided and the estimation accuracy of CKF could also be improved.To the best of the authors' knowledge, there have been few results concerning DSE for generator by using ensemble Kalman filtering algorithm [15], which has been widely utilized in meteorology.
In this paper, based on the improved EnKF-EnSRF, which can approximate the nonlinear system by random sampling method, an adaptive ensemble root mean square Kalman filter algorithm is proposed.The proposed method does not need to interfere with the measured values, avoiding the problem of underestimating the analysis error covariance in EnKF and improving the filtering precision.Meanwhile, the Sage-Husa algorithm is introduced to estimate measurement noise online.When there is a deviation in measurement noise, the method can revise the covariance of measurement noise dynamically and filter it, which reduces the influence of noise error on filtering.Extensive simulation results show that the proposed method can not only achieve a higher filtering precision than the conventional EnKF but also with a stronger robustness to noise.Hence, the main contributions of this paper are:

•
Proposing a novel, robust AEnSRF method applicable for measurement noise estimation.

•
The proposed robust AEnSRF method is based on the combination between the EnSRF and the Sage-Husa algorithm.

•
The proposed AEnSRF does not need to interfere with the measured values, avoiding the problem of underestimating the analysis error covariance in EnKF and improving the filtering precision.

•
The suggested technique utilizes Sage-Husa algorithm to adjust the covariance matrices of measurement noise dynamically and filter it when there is a deviation in measurement noise.

•
The proposed method mitigates the adverse influence of noise error on the filtering result.

•
The proposed method has a higher filtering precision than the conventional EnKF and a stronger robustness to noise.
The rest of the paper is organized as follows: In Section 2, the second-order generator DSE model is presented.Section 3 describes the main algorithm in detail.IEEE 9-bus system and the real power grid system are used to test and illustrate the usefulness of the proposed method, and some necessary comparison results with ENKF and ENSRF are also performed in Section 4. Finally, the conclusions are provided in Section 5.

Mathematical Model of Dynamic State Estimation
State space model can be used to describe the dynamic system as where k indicates the time instant, W k represents system noise, V k represents measurement noise, Φ denotes the state transition matrix, Γ is the input matrix, and C is the observation matrix.The statistical characteristics of noise W k and V k are given as follows where Q is the variance matrix of system process noise, R is the variance matrix of system measurement noise, δ k j = 1, k = j 0, k j .

Equation of Measurement and State Estimation
In the electromechanical transient process of power system, the system topology and bus voltage cannot be acquired in real time.Thus, the rotor power angle and angular velocity, which do not abruptly change and are constrained by the rotor differential equation, are selected as the state variables of the generator DSE.
The second-order differential equations of synchronous generators can be expressed as follows where δ indicates the power angle, ω represents the electric angular velocity; T m denotes the mechanical torque, T e is the electromagnetic torque; P m is the mechanical power, P e is the electromagnetic power, T J represents the generator's rotor's inertia time constant, and D is the damping coefficient.
The differential equation of generator model in (3) can be rewritten as where the unit of δ is degree.For convenience, one has In general, the first three terms of Taylor expansion of general state transition matrix can meet the precision requirements where T represents the sampling period.
The power angle and angular velocity of generator can be measured directly by PMU, which are usually selected as the measurements According to (4) to (7), the discrete form of generator DSE model can be obtained

Error Analysis
Equations ( 4) to (9) show that the imprecision of the model parameters T J , D and the errors in the measurement of P e and P m will lead to the imprecision of the model.The obtained value of T J is generally accurate.In general, the value of D is very small, indicating the influence of mechanical friction and wind resistance on the operation of generators.Therefore, the measurement errors of P e and P m are the main components of generator process noise.Actually, it is difficult to obtain precise P m value, assuming that P m value is constant when there is only governor but no shut-off and fast shut-off valve; the change of P m caused by governor action is treated as the process noise of the system.
According to [16,17], the measurement error is generally between 1% and 2%.Considering the influence, the process noise variance matrix is set as where P e0 is the standard unit of output electromagnetic power in steady state.
The measurement errors of δ and ω are mainly caused by PMU measurement errors.Ideally, the measurement error variance of δ is 2 2 , ω is 10 −6 , and the measurement noise variance matrix of the system is set as R = diag(2 2 , 10 −6 ) (11)

Ensemble Kalman Filter
Regarding the nonlinear systems and in order to track their state trajectories, the extended Kalman filter and the nonlinear filtering method which combines the techniques of unscented transformation and numerical difference without computing Jacobian matrix are widely used.However, when the system's nonlinearity is strong and the noise is non-Gaussian, the performance of these algorithms will be degraded significantly, resulting in unreliable state estimation.Therefore, ensemble Kalman filter, a sub-optimal filtering method, has attracted more and more attention.Based on Kalman filtering, ensemble prediction is introduced into ensemble Kalman filtering.Sequential Monte Carlo simulation method is employed to create the initial set of samples representing the state statistics.The state transition function is used to calculate the new set of states at the current time for each sample in the set, then the samples of the new set of states are calculated.Mean and covariance are used to get the optimal value of the current time estimate [18].
In the traditional Kalman filter, the calculation formulas of state prediction error covariance and analysis error covariance are presented as where X denotes the state vector, P represents the state covariance, T is the matrix transposition.The superscript f indicates the prediction state, the superscript a is the analysis state, and the superscript t is the real state.
In the EnKF, the posterior distribution of the state function is approximated by the ensemble elements, and the degree of approximation depends on the number of elements in the ensemble.Assuming that the number of set elements is N, XE f represents the set storing all predicted states, and XE a represents the set storing all analyzed states.When N tends to infinity, the real value of the state can be replaced by the mean value of the set of state estimators.It follows from [19] that the prediction error covariance and the analysis error covariance can be calculated by The EnKF uses the finite element in the ensemble for the estimation of the state error covariance, thus avoiding the prediction of the covariance matrix.
In the process of EnKF, firstly, according to the known state priori information, an initial state set of N elements is obtained by SMC method, then a prediction set XE f is obtained by predicting each element in the initial state set through Kalman filter, in which the elements are expressed by X f i .Secondly, the observation set of state by SMC method as well, and the prediction state set is modified by covariance of XE f to get the analysis set XE a of state.The elements are expressed by X a i , and the average value of the state analysis set XE a is the optimal state estimation value.The analysis set of states is used for prediction [20,21].
The iteration formula of ensemble Kalman filter algorithm is given as follows (1) State prediction: (2) State update: where X are the deviation matrices of elements and average values in predicted state set and measured state set, respectively, P f XZ,k+1 is the covariance of predictive state error and observation state error at k + 1 time, P f ZZ,k+1 is the covariance of observation state error at k + 1 time.
The EnKF algorithm not only overcomes the problem that Kalman filter is limited to linear systems, but also avoids the problem of large amount of computing resources when calculating prediction error covariance.Tangent linear model and adjoint model are no longer needed.The method also gives the optimal value of the estimated result and the confidence interval of the estimated value at the same time.
The EnKF is a sub-optimal estimation method which uses the ensemble to describe the traditional Kalman filter algorithm.A prediction set estimation is used to calculate the state prediction error covariance matrix of the Kalman gain matrix.If each element of the set is updated with the same observation value and the same gain, the set will systematically underestimate the state update error covariance matrix.It even can lead to subsequent analysis degradation and filtering divergence.This problem can be alleviated by adding random perturbations to the measurements with correct statistical data.However, the introduction of random perturbation will inevitably increase the additional error of sampling error related to measurement error covariance.This additional error will reduce the covariance accuracy of the observation state error as well as the estimation accuracy.

Ensemble Square Root Kalman Filter
In EnSRF, the calculation of the analytic state variables is divided into two parts [22], one is the update of the mean value of the analytic state set, the other is the update of the mean deviation of the analytic state set sample where X a represents the mean of analysis state set, X a indicates the mean deviation of analysis state set samples, X f denotes the mean deviation of prediction state set, X f is the mean deviation of prediction state set samples, Z is the observation value obtained by PMU, Z is the random observation deviation obeying the probability distribution of observation error, K is the gain of traditional Kalman filter, K is the mean deviation gain of the update set.
In EnKF, when all sets of samples are updated with the same observation value and the same gain [23], i.e., Z = 0 and K = K, the covariance of the analysis state set can be expressed as (25).In the traditional Kalman filter, the covariance expression of state analysis error is shown in (26).As can be seen from ( 25) and (26), the value of P a in EnKF lacks the term KRK T , so EnKF has the problem of systematic underestimation EnSRF eliminates the problem of systematic underestimation of EnKF without disturbing the observed values [22][23][24].In EnSRF, (24) can be reformulated as K is substituted for (25) and the analysis error covariance of the set K needs to satisfy (26).Furthermore, if K = αK and α is a constant, one can further derive The traditional Kalman filtering algorithm is a linear unbiased minimum variance estimation algorithm under standard conditions.Kalman filter can obtain better estimation results when the dynamic system model and noise statistics are known.However, in the practical application, it is difficult to obtain accurate system model as well as accurate noise statistical characteristics, which will affect the accuracy of Kalman filter estimation results.

Sage-Husa Algorithm
In 1969, Sage and Husa presented an improved noise filtering algorithm.The algorithm can compute both matrices of process noise variance Q and measurement noise variance R [25,26] online and in real time when the system noise mean and covariance are unknown.Comparing with the adaptive filtering algorithm based on maximum likelihood criterion, this algorithm adds a forgetting factor, which strengthens the influence of newer observations on filtering and weakens the influence of older observations on the results.This algorithm effectively reduces the influence of model error on filtering, restrains filter divergence, and improves estimation accuracy [27].
The steps of Sage-Husa algorithm are given as: (1) Filtering equation: (2) Time-varying noise estimator: where qk+1 represents the mean of process noise at k + 1; rk+1 is the mean of k + 1, e k+1 is the residual ), b is the forgetting factor and 0 < b < 1.Generally, Sage-Husa algorithm can simultaneously estimate process noise Qk+1 and measurement noise Rk+1 .In fact, the algorithm cannot estimate both Qk+1 and Rk+1 at the same time.It can be seen from ( 37) and (39) that the calculation of Qk+1 and Rk+1 depends on e k+1 .The change of e k+1 will affect the calculation of Qk+1 and Rk+1 at the same time, which cannot guarantee the accuracy of the estimation results.In addition, there are minus signs in the calculation of Qk+1 and Rk+1 , which may affect the semi-positive and positive definiteness of Qk+1 and Rk+1 , leading to filter divergence.
In the electromechanical transient process of generator, the influence of measurement noise Rk+1 on the estimation is more important.Generally speaking, process noise Q is a constant.Assuming that both process noise as well as measurement noise are Gauss, white noise with mean value equals to zero, i.e., r = 0 and q = 0.
Sage-Husa algorithm can be simplified as follows.
(1) Filtering equation: Xk+1|k = F k Xk|k (41) (2) Measurement noise estimator: In order to ensure the positive covariance of measurement noise, (47) is replaced by (48).It can be seen from the formula that a certain filtering-precision is sacrificed to ensure the stability of the filter.

Adaptive Ensemble Square Root Kalman Filter
The proposed AEnSRF combines ensemble root mean square Kalman filter with Sage-Husa algorithm.Comparing with ensemble Kalman filter, the proposed AEnSRF not only has a higher filtering-precision but also can be effectively corrected while measurement noise deviates, so as to suppress filter divergence and improve estimation accuracy.
The specific steps are as follows: (1) Initialization: The initial value ω 0 of generator rotor angular velocity is 1.The initial value δ of generator power angle, electromagnetic power P e and mechanical power P m is its steady-state operation value.(2) The initial value of the state set with the number of set elements is generated by Monte Carlo method, and the initial estimation error covariance matrix is taken as the unit matrix.

Basic Test Analysis
The IEEE-9 bus system is firstly taken as the test system in this paper.The single-line diagram of it is given in Figure 2. The classical second-order model of generator considering governor is adopted in simulation.The inertia time constants of G1, G2, and G3 are respectively 47.28, 12.8, 6.02 s; and the damping coefficient is set to 2. The initial values of generator state variables and control variables are set to the values of the system's stable operation before fault occurred.Suppose that in

Basic Test Analysis
The IEEE-9 bus system is firstly taken as the test system in this paper.The single-line diagram of it is given in Figure 2. The classical second-order model of generator considering governor is adopted in simulation.The inertia time constants of G1, G2, and G3 are respectively 47.28, 12.8, 6.02 s; and the damping coefficient is set to 2. The initial values of generator state variables and control variables are set to the values of the system's stable operation before fault occurred.Suppose that in the IEEE 9 bus system, three-phase short-circuit fault occurs at the beginning of branch 4-8 at 0.8 s.After 0.36 s, the switch on both sides of branch jumps off and the fault is removed.The whole simulation time is 6 s.The real values of power angle and angular velocity of generator are obtained by Bonneville Power Administration (BPA) simulation software, and the measured values are obtained by adding the real values to the random noise where its mean value obeys normal distribution with zero standard deviation.In this paper, the DSE algorithm is implemented in Matlab.The initial values of the estimated generator state variables are set to stable operation value, and the initial state error covariance matrix is set to unit matrix.
In order to further compare the difference of filtering results of each algorithm, the root mean square error (RMSE) is used where n is the sampling period.The process noise variance matrix and the measurement noise variance matrix of the generator are selected as the corresponding values in [10] and [11].Figure 3 is a comparison of the filtering results of G1 by utilizing the EnKF and EnSRF algorithms under normal noise.Table 1 gives the RMSE results of EnKF and EnSRF algorithms.The real values of power angle and angular velocity of generator are obtained by Bonneville Power Administration (BPA) simulation software, and the measured values are obtained by adding the real values to the random noise where its mean value obeys normal distribution with zero standard deviation.In this paper, the DSE algorithm is implemented in Matlab.The initial values of the estimated generator state variables are set to stable operation value, and the initial state error covariance matrix is set to unit matrix.
In order to further compare the difference of filtering results of each algorithm, the root mean square error (RMSE) is used where n is the sampling period.The process noise variance matrix and the measurement noise variance matrix of the generator are selected as the corresponding values in [10,11].Figure 3 is a comparison of the filtering results of G1 by utilizing the EnKF and EnSRF algorithms under normal noise.Table 1 gives the RMSE results of EnKF and EnSRF algorithms.The real values of power angle and angular velocity of generator are obtained by Bonneville Power Administration (BPA) simulation software, and the measured values are obtained by adding the real values to the random noise where its mean value obeys normal distribution with zero standard deviation.In this paper, the DSE algorithm is implemented in Matlab.The initial values of the estimated generator state variables are set to stable operation value, and the initial state error covariance matrix is set to unit matrix.
In order to further compare the difference of filtering results of each algorithm, the root mean square error (RMSE) is used where n is the sampling period.
The process noise variance matrix and the measurement noise variance matrix of the generator are selected as the corresponding values in [10] and [11].Figure 3 is a comparison of the filtering results of G1 by utilizing the EnKF and EnSRF algorithms under normal noise.Table 1 gives the RMSE results of EnKF and EnSRF algorithms.As it might be noted from Figure 3 and Table 1, EnKF and AEnSRF both show good filtering effect without deviation of measurement noise.Because the observation value in AEnSRF does not introduce disturbance, the problem of underestimation in EnKF filtering is avoided.The filteringprecision of AEnSRF is improved comparing with EnKF.
Figure 4 is a sketch of the noise estimator tracking the measured noise in AEnSRF without deviation from the measured noise.As it might be noted from Figure 3 and Table 1, EnKF and AEnSRF both show good filtering effect without deviation of measurement noise.Because the observation value in AEnSRF does not introduce disturbance, the problem of underestimation in EnKF filtering is avoided.The filtering-precision of AEnSRF is improved comparing with EnKF.
Figure 4 is a sketch of the noise estimator tracking the measured noise in AEnSRF without deviation from the measured noise.As it might be noted from Figure 3 and Table 1, EnKF and AEnSRF both show good filtering effect without deviation of measurement noise.Because the observation value in AEnSRF does not introduce disturbance, the problem of underestimation in EnKF filtering is avoided.The filteringprecision of AEnSRF is improved comparing with EnKF.
Figure 4 is a sketch of the noise estimator tracking the measured noise in AEnSRF without deviation from the measured noise.From the simulation, we can see that in the process of tracking power angle measurement noise by AEnSRF, the standard deviation of 0.14 s is stable at 1.9 degrees, and the angular velocity is too small (almost 0).Combining with Figure 4, it can be remarked that AEnSRF can track the measured noise more accurately.

Noise Measurement and Analysis on the IEEE 9-Bus System
In this subsection, G2 is taken as an example to test the estimation effect of AEnSRF on generator state as well as the tracking of measurement noise when it is biased.Assuming that the actual value of the standard deviation of power angle measurement noise is 3 degrees, the value is set as 2 degrees in the filtering program to detect the tracking effect of AEnSRF on power angle measurement noise.
Figure 5 depicts the comparison results of G2 by utilizing EnKF and AEnSRF algorithms, respectively.It can be obviously noted from Figure 5 that both of them can show good estimation results under the consideration of measurement noise deviation condition.
From the simulation, we can see that in the process of tracking power angle measurement noise by AEnSRF, the standard deviation of 0.14 s is stable at 1.9 degrees, and the angular velocity is too small (almost 0).Combining with Figure 4, it can be remarked that AEnSRF can track the measured noise more accurately.

Noise Measurement and Analysis on the IEEE 9-Bus System
In this subsection, G2 is taken as an example to test the estimation effect of AEnSRF on generator state as well as the tracking of measurement noise when it is biased.Assuming that the actual value of the standard deviation of power angle measurement noise is 3 degrees, the value is set as 2 degrees in the filtering program to detect the tracking effect of AEnSRF on power angle measurement noise.
Figure 5 depicts the comparison results of G2 by utilizing EnKF and AEnSRF algorithms, respectively.It can be obviously noted from Figure 5 that both of them can show good estimation results under the consideration of measurement noise deviation condition.The simulation results of G2 show that under AEnSRF algorithm, the RMSE of power angle estimation is 0.0516, the RMSE of electric angular velocity estimation is 0.0011.Under EnKF algorithm, the RMSE of power angle estimation is 0.0536 and the RMSE of electric angular velocity estimation is 0.0012.By comparing both algorithms, AEnSRF has higher filtering-precision.
Figure 6 is the simulation diagram of AEnSRF tracking curve for measurement noise when there is a deviation in measurement noise.As it can be noticed from this figure, when the power angle measurement noise is biased, the Sage-Husa noise estimator with on-line real-time estimation noise  The simulation results of G2 show that under AEnSRF algorithm, the RMSE of power angle estimation is 0.0516, the RMSE of electric angular velocity estimation is 0.0011.Under EnKF algorithm, the RMSE of power angle estimation is 0.0536 and the RMSE of electric angular velocity estimation is 0.0012.By comparing both algorithms, AEnSRF has higher filtering-precision.
Figure 6 is the simulation diagram of AEnSRF tracking curve for measurement noise when there is a deviation in measurement noise.As it can be noticed from this figure, when the power angle measurement noise is biased, the Sage-Husa noise estimator with on-line real-time estimation noise is introduced by AEnSRF, which makes the estimated standard deviation of power angle measurement noise stabilize around 3 degrees at about 1.1 s. is introduced by AEnSRF, which makes the estimated standard deviation of power angle measurement noise stabilize around 3 degrees at about 1.1 s.

Noise Measurement and Analysis on a Real Power System
For further verification of the proposed algorithm's effectiveness, simulation tests are also carried out on a real power grid system.One generator outgoing circuit has a short-circuit fault, and the generator parameters are based on the actual data.The initial value of the state variable is the steady value, and the initial error covariance matrix is the unit matrix.Assuming that the actual value of the standard deviation of power angle measurement noise is 2.5 degrees, the value is set to 2 degrees in the filtering program.
Comparing EnKF with the proposed method, Figures 7 and 8 show the simulation results of DSE of generator in the presence of measurement noise deviation.The simulation results of G2 illustrate that under the AEnSRF, the RMSE of power angle estimation is 0.0435, the RMSE of electric angular velocity estimation is 0.0010.Under the EnKF, the RMSE of power angle estimation is 0.0438 and the RMSE of electric angular velocity estimation is 0.0013.

Noise Measurement and Analysis on a Real Power System
For further verification of the proposed algorithm's effectiveness, simulation tests are also carried out on a real power grid system.One generator outgoing circuit has a short-circuit fault, and the generator parameters are based on the actual data.The initial value of the state variable is the steady value, and the initial error covariance matrix is the unit matrix.Assuming that the actual value of the standard deviation of power angle measurement noise is 2.5 degrees, the value is set to 2 degrees in the filtering program.
Comparing EnKF with the proposed method, Figures 7 and 8 show the simulation results of DSE of generator in the presence of measurement noise deviation.The simulation results of G2 illustrate that under the AEnSRF, the RMSE of power angle estimation is 0.0435, the RMSE of electric angular velocity estimation is 0.0010.Under the EnKF, the RMSE of power angle estimation is 0.0438 and the RMSE of electric angular velocity estimation is 0.0013.Comparing with EnKF, AEnSRF has higher filtering-precision.In AEnSRF, the noise estimator tracks the power angle measurement noise at 0.3 s and stabilizes at about 2.5 degrees.This shows that AEnSRF has a good tracking effect on measurement noise.

Conclusions
In this paper, a novel DSE method was proposed based on the EnSRF; a simplified Sage-Husa adaptive Kalman filter with relatively simple principle and good practicability was introduced.In the proposed AEnSRF method, set members were utilized to approximate the posterior distribution of the real state without random disturbance to the measured values.Comparing with EnKF, AEnSRF could enhance the accuracy of filtering.At the same time, Sage-Husa noise estimator was added between the prediction and correction steps to estimate the measurement noise online and in real time, which effectively avoided the problem that the filtering-precision was reduced or even divergent due to the deviation of the measurement noise.Comparing with EnKF, AEnSRF has higher filtering-precision.In AEnSRF, the noise estimator tracks the power angle measurement noise at 0.3 s and stabilizes at about 2.5 degrees.This shows that AEnSRF has a good tracking effect on measurement noise.

Conclusions
In this paper, a novel DSE method was proposed based on the EnSRF; a simplified Sage-Husa adaptive Kalman filter with relatively simple principle and good practicability was introduced.In the proposed AEnSRF method, set members were utilized to approximate the posterior distribution of the real state without random disturbance to the measured values.Comparing with EnKF, AEnSRF could enhance the accuracy of filtering.At the same time, Sage-Husa noise estimator was added between the prediction and correction steps to estimate the measurement noise online and in real time, which effectively avoided the problem that the filtering-precision was reduced or even divergent due to the deviation of the measurement noise.
i elements of prediction state set and measurement state set, respectively; X a i k are the i elements of analysis state set, X f k+1 , X a k+1 and Z f k+1 are the average values of predicted state set, analyzed state set and measured state set at k + 1 time instant, respectively, E f X,k+1 and E f Z,k+1

Figure 2 .
Figure 2. The single-line diagram of IEEE-9 bus test system.

Figure 2 .
Figure 2. The single-line diagram of IEEE-9 bus test system.

Figure 2 .
Figure 2. The single-line diagram of IEEE-9 bus test system.

Figure 4 .
Figure 4. Measurement noise tracking curve of G1 under AEnSRF: (a) Power angle measurement noise tracking curve of G1; (b) angular velocity measurement noise tracking curve of G1.

Figure 3 .
Figure 3.Comparison of G1 state filtering results under ensemble Kalman filter (EnKF) and AEnSRF: (a) Power angle dynamic estimation curve of G1; (b) electric angular velocity dynamic estimation curve of G1.

Figure 3 .
Figure 3.Comparison of G1 state filtering results under ensemble Kalman filter (EnKF) and AEnSRF: (a) Power angle dynamic estimation curve of G1; (b) electric angular velocity dynamic estimation curve of G1.

Figure 4 .
Figure 4. Measurement noise tracking curve of G1 under AEnSRF: (a) Power angle measurement noise tracking curve of G1; (b) angular velocity measurement noise tracking curve of G1.

Figure 4 .
Figure 4. Measurement noise tracking curve of G1 under AEnSRF: (a) Power angle measurement noise tracking curve of G1; (b) angular velocity measurement noise tracking curve of G1.

Figure 5 .
Figure 5. Comparisons of G2 state filtering results between EnKF and AEnSRF in the presence of measurement noise deviation: (a) Power angle dynamic estimation curve of G2; (b) electric angular velocity dynamic estimation curve of G2.

Figure 5 .
Figure 5. Comparisons of G2 state filtering results between EnKF and AEnSRF in the presence of measurement noise deviation: (a) Power angle dynamic estimation curve of G2; (b) electric angular velocity dynamic estimation curve of G2.

Figure 6 .
Figure 6.Tracking curve of G2 measurement noise by AEnSRF in the case of measurement noise deviation: (a) Power angle measurement noise tracking curve of G2; (b) angular velocity measurement noise tracking curve of G2.

Figure 6 .
Figure 6.Tracking curve of G2 measurement noise by AEnSRF in the case of measurement noise deviation: (a) Power angle measurement noise tracking curve of G2; (b) angular velocity measurement noise tracking curve of G2.

Figure 7 .Figure 7 .
Figure 7. Contrast chart of EnKF and AEnSRF for generator state filtering with deviation of measurement noise: (a) Power angle dynamic estimation curve of generator; (b) dynamic angular velocity estimation curve of generator.

Figure 7 .Figure 8 .
Figure 7. Contrast chart of EnKF and AEnSRF for generator state filtering with deviation of measurement noise: (a) Power angle dynamic estimation curve of generator; (b) dynamic angular velocity estimation curve of generator.

Figure 8 .
Figure 8. Tracking curve of AEnSRF for generator measurement noise with deviation of measurement noise: (a) Power angle measurement noise tracking curve of generator; (b) angular velocity measurement noise tracking curve of generator.

Table 1 .
The RMSE results of EnKF and EnSRF.

Table 1 .
The RMSE results of EnKF and EnSRF.

Table 1 .
The RMSE results of EnKF and EnSRF.