Login| Sign Up| Help| Contact|

Patent Searching and Data


Title:
DYNAMIC STATE ESTIMATION OF AN OPERATIONAL STATE OF A GENERATOR IN A POWER SYSTEM
Document Type and Number:
WIPO Patent Application WO/2019/016547
Kind Code:
A1
Abstract:
Example embodiments described herein are directed towards dynamic state estimation of an operating state of a generator in a power system. Such estimation is performed for an individual generator in real time with improved accuracy and without the use of Global Position System (GPS) synchronization.

Inventors:
SINGH ABHINAV (GB)
PAL BIKASH (GB)
Application Number:
PCT/GB2018/052032
Publication Date:
January 24, 2019
Filing Date:
July 18, 2018
Export Citation:
Click for automatic bibliography generation   Help
Assignee:
IMPERIAL INNOVATIONS LTD (GB)
International Classes:
H02J3/24; G01R31/34
Other References:
ESMAEIL GHAHREMANI ET AL: "Dynamic State Estimation in Power System by Applying the Extended Kalman Filter With Unknown Inputs to Phasor Measurements", 1 November 2011, IEEE TRANSACTIONS ON POWER SYSTEMS, IEEE SERVICE CENTER, PISCATAWAY, NJ, US, PAGE(S) 2556 - 2566, ISSN: 0885-8950, XP011387422
KAMWA I ET AL: "Performance of Demodulation-Based Frequency Measurement Algorithms Used in Typical PMUs", 1 April 2004, IEEE TRANSACTIONS ON POWER DELIVERY, IEEE SERVICE CENTER, NEW YORK, NY, US, PAGE(S) 505 - 514, ISSN: 0885-8977, XP011109778
ABHINAV KUMAR SINGH: "Decentralized estimation and control for power systems", 1 January 2014 (2014-01-01), XP055509622, Retrieved from the Internet [retrieved on 20180925]
AMINIFAR FARROKH ET AL: "Power System Dynamic State Estimation With Synchronized Phasor Measurements", 1 February 2014, IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT, IEEE SERVICE CENTER, PISCATAWAY, NJ, US, PAGE(S) 352 - 363, ISSN: 0018-9456, XP011535853
AMINIFAR FARROKH ET AL: "Synchrophasor Measurement Technology in Power Systems: Panorama and State-of-the-Art", IEEE ACCESS, vol. 2, 9 January 2015 (2015-01-09), pages 1607 - 1628, XP011570876, DOI: 10.1109/ACCESS.2015.2389659
Attorney, Agent or Firm:
HGF LIMITED (LONDON OFFICE) (IM)
Download PDF:
Claims:
CLAIMS

A method for Dynamic State Estimation, DSE, of an operational state of a generator in a power system, the method comprising:

receiving measured voltage and current analog signals associated with the generator;

sampling the received measured signals to voltage and current discrete signals;

estimating magnitude, phase and frequency variables of the voltage and current, as well as associated voltage and current variances, respectively, for each variable, using the discrete voltage and current signals as an input in a Discrete Fourier Transform, DFT;

calculating a power variable and associated power variance based on the estimated magnitude, phase and frequency variables of the voltage and current, as well as associated voltage and current variances for each variable; and

estimating the dynamic state of the generator using at least a subset of the estimated magnitude, phase and frequency variables of the voltage and current, as well as associated voltage and current variances for each variable, and the calculated power variable and the associated power variance using a state estimator.

The method of claim 1 , wherein the estimating the dynamic state of the generator is based on, at least in part, estimating a relative angle as a difference between a rotor angle and the estimated voltage phase.

The method of claim 2, wherein the estimating the relative angle further comprises estimating the relative angle according to Αάι = (ωι - f ), where a is the difference between the rotor angle and the estimated voltage phase, ω is rotor speed and fy is the frequency variable of the voltage of the fh generator. 4. The method of any of claims 1-3, wherein the sampling further comprises multiplying the voltage and current discrete signals with a window function.

5. The method of claim 4, wherein the multiplying with the window function further comprises multiplying according to Z(2) =∑¾=ο Y^ e N~ =— eieW [A -—) -

2] v fs '

— e~J9W (λ +— ), where Z is a DFT of a product of Y and h, N denotes total number of samples for fining the DFT, /( denotes a discrete sample, Y denotes a sinusoidal signal with harmonics and noise, h is a window function, W is a DFT of the window function, f is the frequency of Vs fundamental component in Hz, fs is a sampling frequency for DFT in Hz, and λ e {0, 1, (N - 1)}.

6. The method according to any of claims 1-5, wherein the subset of the estimated magnitude, phase and frequency variables of the voltage and current is the estimated magnitude of the voltage, the estimated frequency of the voltage and the estimated magnitude of the current.

7. The method according to claim 6, wherein the estimated magnitude of the voltage and the estimated frequency of the voltage and associated voltage variances are inputs to the state estimator, and are utilized as pseudo-inputs, and are provided in differential equations of the state estimator, and the estimated magnitude of the current and the calculated power and associated current and power variances are provided as an algebraic equation and are given as measurement inputs to the state estimator.

The method according to claim 7, wherein the differential equations are represented as xlk = xlk + T0gl(xlk , iilk , wlk) + vlk = xlk = gl(xlk, ulk, wlk) + vlk, where x is a column vector of a state, k = (k - 1), To is a sampling period in seconds, g is a discrete form of a differential function, u is a column vector of pseudo-inputs, w is a column vector of noise, /' is the ith generation unit, and k is the kth sample and the algebraic equation is represented as ik yik _

hl{xlk, ulk, wlk) + wlk, where y is a column vector of a measurement, Pe denotes active electrical-power output of a machine, / and \m represents an analogue stator current and its magnitude, respectively, ¼ and Vq represents d-axis and q-axis stator voltages, respectively, Id and lq represents d-axis and q-axis stator currents, respectively.

The method according to any of claims 1-5, wherein the subset of the estimated magnitude, phase and frequency variables of the voltage and current is the estimated magnitude of the current, the estimated frequency of the current and the estimated magnitude of the voltage.

10. The method according to claim 9, wherein the estimated magnitude of the current and the estimated frequency of the current and associated current variances are inputs to the state estimator and are used in differential equations of the state estimator, and the estimated magnitude of the voltage and the calculated power and associated voltage and power variances are represented as an algebraic equation and is given as a measurement input to the state estimator.

1 1. The method according to any of claims 1-5, wherein the subset of the estimated magnitude, phase and frequency variables of the voltage and current is the estimated magnitude of the voltage, the estimated magnitude of the current and the estimated frequency of the voltage.

12. The method according to claim 11 , wherein the estimated magnitude of the voltage and the estimated magnitude of the current and associated voltage and current variances are inputs to the state estimator and are used in differential equations of the state estimator, and the estimated frequency of the voltage and the calculated power and associated voltage and power variances are represented as an algebraic equation and is given as a measurement input to the state estimator, wherein the calculated power is reactive power.

13. The method according to any of claims 1-5, wherein the subset of the estimated magnitude, phase and frequency variables of the voltage and current is the estimated magnitude of the current, the estimated frequency of the current and the estimated magnitude of the voltage.

14. The method according to claim 13, wherein the estimated magnitude of the current and the estimated frequency of the current and associated current variances are inputs to the state estimator and are used in differential equations of the state estimator, and the estimated magnitude of the voltage and the calculated power and associated voltage and power variances are represented as an algebraic equation and is given as a measurement input to the state estimator.

15. The method according to any of claims 1-14, wherein the calculated power is real power or reactive power of the generator.

16. The method according to any of claims 1-15, wherein the estimation of the

magnitude, phase and frequency variables of the voltage and current, as well as associated variance for each variable further comprises estimating according to / =

jin N jin N jin N

fs lz0+2z1+9z2. z0 _ eieB+e~iec _ l-e fs _ l-e fs _ l-e fs

Z0-2Z1 + Z2 ' Z ~ eJdE + e-JdF' f f ]3 ' ~ f f ]3 ' ~ ίΕ- _\ίΕ_ Ϋ'

fs [fsl fs [fsl fs [fs 1 β _ lz0F-Z C _ 1 (ZpF-ZiC 7 _ 8πΖο . ΓΩΩίΐλ - fN •6 - ^ZlB-Z0E 2jLnlz1B-Z0E) ' Ym - [N{Befl + Ce-fl

+ 1 }]'LKUV) ~ 2nJ Y N(.N2-1) MJ N J jV0V-l) > fa 771

2CRB(?m);d0 = 6CRB(§) , where Y denotes a sinusoidal signal with harmonics and noise, h is a window function, Z is a DFT of the product of Y and h, f is a frequency of Ys fundamental component in Hz, Ym is a magnitude of Y's fundamental component, fs is a sampling frequency for DFT in Hz, N is a total number of samples for finding DFT, CRB is a Cramer-Rao bound, σ denotes standard deviation with o2 as its variance, Θ denotes a phase in Ys fundamental component, and fo denotes a voltage base value in Hz.

The method according to any of claims 1-16, wherein the calculation of the power variable and the associated variance further comprises calculating according to

Peik = cos(¾fe - § k); σρ = [σ ( *)2 + cos2fe - §jk) +

{Vm )2{lm)2 [ ik + where Pe is active electrical power output of

the generator, Vm is voltage magnitude, \m is current magnitude, θν, and θι, are phases of stator voltage and current, respectively, σ denotes standard deviation and o2 its variance.

The method according to any of claims 1-17, wherein the estimating of the DSE further comprises estimating according to

e

electrical power output of the generator, Vm is voltage magnitude, \m is current magnitude, θν, and θι, are phases of stator voltage and current, respectively, σ denotes standard deviation and o2 denotes variance, ω and ω' represent column vectors of noise in y and u, respectively, P denotes a covariance matrix, f denotes frequency of Y's fundamental component in Hz, u and v represent column vectors of pseudo-inputs and process noise, respectively, represents a composite state vector, g denotes a discrete form of differential functions, / denotes a column vector of measurements, /' denotes the ith generator unit, and k denotes a discrete sample.

19. The method according to any of claims 1-18, wherein the DFT is an interpolated DFT.

20. The method according to any of claims 1-19, wherein the state estimator is an

Unscented Kalman filter, UKF.

An apparatus for Dynamic State Estimation, DSE, of an operational state of a generator in a power system, the apparatus comprising:

a transceiver to receive measured voltage and current analog signals associated with the generator;

a processor to sample the received measured signals to voltage and current discrete signals;

the processor to estimate magnitude, phase and frequency variables of the voltage and current, as well as associated voltage and current variances, respectively, for each variable, using the discrete voltage and current signals as an input in a Discrete Fourier Transform, DFT;

the processor to calculate a power variable and associated power variance based on the estimated magnitude, phase and frequency variables of the voltage and current, as well as associated voltage and current variances for each variable; and

the processor to estimate the dynamic state of the generator using at least a subset of the estimated magnitude, phase and frequency variables of the voltage and current, as well as associated voltage and current variances for each variable, and the calculated power variable and the associate power variance using a state estimator.

22. The apparatus of claim 21 , wherein the processor is to estimate the dynamic state of the generator based on, at least in part, an estimation of a relative angle as a difference between a rotor angle and the estimated voltage phase.

23. The apparatus of claim 22, wherein the processor is to estimate the relative angle the relative angle according to Αάι = (ωι - f ), where a is the difference between the rotor angle and the estimated voltage phase, ω is rotor angle and fy is the frequency variable of the voltage of the fh generator.

24. The apparatus of any of claims 21-23, wherein the processor is to sample via a multiplication of the voltage and current discrete signals with a window function.

25. The apparatus of claim 24, wherein the processor is to multiply the voltage and

current discrete signals with the window function according to Z(X) =

(λ + ), where Z is a DFT of a product of Y and h, N denotes total number of samples for fining the DFT, k denotes a discrete sample, Y denotes a sinusoidal signal with harmonics and noise, h is a window function, W is a DFT of the window function, f is a frequency of Vs fundamental component in Hz, fs is a sampling frequency for DFT in Hz, and λ e {0, 1 (N - 1)}.

26. The apparatus according to any of claims 21-25, wherein the subset of the estimated magnitude, phase and frequency variables of the voltage and current is the estimated magnitude of the voltage, the estimated frequency of the voltage and the estimated magnitude of the current.

27. The apparatus according to claim 26, wherein the estimated magnitude of the

voltage and the estimated frequency of the voltage and associated voltage variances are inputs to the state estimator, and are utilized as pseudo-inputs, and are provided in differential equations of the state estimator, and the estimated magnitude of the current and the calculated power and associated current and power variances are provided as an algebraic equation and are given as measurement inputs to the state estimator.

28. The apparatus according to claim 24, wherein the differential equations are

represented as xlk = xlJc + T0gl(xlli, ulli, wlli) + vlJc = xlk = gl(xlJi , ulJi , wlJc) + vlJc, where x is a column vector of a state, k = (k - 1), To is a sampling period in seconds, g is a discrete form of a differential function, u is a column vector of pseudo-inputs, w is a column vector of noise, /' is the ith generation unit, and k is the Qik

e

kth sample and the algebraic equation is represented as y ik _ r

fik

, where y is a column vector of

a measurement, Pe denotes active electrical-power output of a machine, / and \m represents an analogue stator current and its magnitude, respectively, ¼ and Vq represents d-axis and q-axis stator voltages, respectively, Id and lq represents d-axis and q-axis stator currents, respectively.

The apparatus according to any of claims 21 -25, wherein the subset of the estimated magnitude, phase and frequency variables of the voltage and current is the estimated magnitude of the current, the estimated frequency of the current and the estimated magnitude of the voltage.

The apparatus according to claim 29, wherein the estimated magnitude of the current and the estimated frequency of the current and associated current variances are inputs to the state estimator and are used in differential equations of the state estimator, and the estimated magnitude of the voltage and the calculated power and associated voltage and power variances are represented as an algebraic equation and is given as a measurement input to the state estimator.

The apparatus according to any of claims 21 -25, wherein the subset of the estimated magnitude, phase and frequency variables of the voltage and current is the estimated magnitude of the voltage, the estimated magnitude of the current and the estimated frequency of the voltage. 32. The apparatus according to claim 31 , wherein the estimated magnitude of the

voltage and the estimated magnitude of the current and associated voltage and current variances are inputs to the state estimator and are used in differential equations of the state estimator, and the estimated frequency of the voltage and the calculated power and associated voltage and power variances are represented as an algebraic equation and is given as a measurement input to the state estimator, wherein the calculated power is reactive power.

The apparatus according to any of claims 21 -25, wherein the subset of the estimated magnitude, phase and frequency variables of the voltage and current is the estimated magnitude of the current, the estimated frequency of the current and the estimated magnitude of the voltage.

The apparatus according to claim 33, wherein the estimated magnitude of the current and the estimated frequency of the current and associated current variances are inputs to the state estimator and are used in differential equations of the state estimator, and the estimated magnitude of the voltage and the calculated power and associated voltage and power variances are represented as an algebraic equation and is given as a measurement input to the state estimator.

The apparatus according to any of claims 21 -34, wherein the calculated power is real power or reactive power of the generator.

The apparatus according to any of claims 21 -35, wherein the processor is to estimate the magnitude, phase and frequency variables of the voltage and current, as well as associated variance for each variable according to f =— iz°+2Zi +9¾

° ' N Z0 -2Z1 + Z2

f h \ _ _ CRB(? ) = 2*1 = 2 = _ 2CR^) 2 =

2nJ Y N(.N2 -1) M J N J Y£N(.N- 1) J o2 Ym

2CRB(?m); 8g = 6CRB (§) , where Y denotes a sinusoidal signal with harmonics and noise, h is a window function, Z is a DFT of the product of Y and h, f is a frequency of Vs fundamental component in Hz, Ym is a magnitude of Y's fundamental component, fs is a sampling frequency for DFT in Hz, N is a total number of samples for finding DFT, CRB is a Cramer-Rao bound, σ denotes standard deviation with o2 as its variance, Θ denotes a phase in Vs fundamental component, and fo denotes a voltage base value in Hz.

37. The apparatus according to any of claims 21 -36, wherein the processor is to

calculate the power variable and the associated power variance according to eife = t¾fe fe cos(¾fe - §jk ; dH2 k = [σ¾( *)2 + cos2fe - §jk) +

{Vm )2 {lm )2 [ ik + where Pe is active electrical power output of

the generator, Vm is voltage magnitude, \m is current magnitude, θν, and θι are phases of stator voltage and current, respectively, σ denotes standard deviation and o2 its variance.

The apparatus according to any of claims 21-37, wherein the processor is to estimate the DSE according to w ik _ pik _

cos{Bik - §jk); <¾k = [¾fc( fe)2 + cos2fe - §jk)

. pik _

(WW [<¾* + ¾*] sin W - §n ; X* = [^] ; k =

; Xlk = gl(xlJc, ulJc) + vlJc ylk = hl(xlk, ulk) + ωιΙζ, where Pe is active

electrical power output of the generator, Vm is voltage magnitude, \m is current magnitude, θν, and θι are phases of stator voltage and current, respectively, σ denotes standard deviation and o2 its variance, ω and ω' represent column vectors of noise in y and u, respectively, P denotes a covariance matrix, f denotes frequency of Y's fundamental component in Hz, u and v represent column vectors of pseudo- inputs and process noise, respectively, X represents a composite state vector, g denotes a discrete form of differential functions, y denotes a column vector of measurements, /' denotes the ith generator unit, and k denotes a discrete sample.

39. The apparatus according to any of claims 21-38, wherein the DFT is an interpolated DFT.

40. The apparatus according to any of claims 21-39, wherein the state estimator is an Unscented Kalman filter, UKF. 41. A computer readable medium having executable instructions stored thereon which, when executed by an apparatus for Dynamic State Estimation, DSE, of an operational state of a generator in a power system, cause the apparatus to:

receive measured voltage and current analog signals associated with the generator;

sample the received measured signals to voltage and current discrete signals; estimate magnitude, phase and frequency variables of the voltage and current, as well as associated voltage and current variances, respectively, for each variable, using the discrete voltage and current signals as an input in a Discrete Fourier Transform, DFT; calculate a power variable and associated power variance based on the estimated magnitude, phase and frequency variables of the voltage and current, as well as associated voltage and current variances for each variable; and

estimate the dynamic state of the generator using at least a subset of the estimated magnitude, phase and frequency variables of the voltage and current, as well as associated voltage and current variances for each variable, and the calculated power variable and the associate power variance using a state estimator.

The computer readable medium of claim 41 , wherein the instructions further cause the apparatus to estimate the dynamic state of the generator is based on, at least in part, an estimation of a relative angle as a difference between a rotor angle and the estimated voltage phase.

The computer readable medium of claim 36, wherein the instructions further cause the apparatus to estimate the relative angle according to Αάι = (ωι - f ), where a is the difference between the rotor angle and the estimated voltage phase, ω is rotor speed and fy is the frequency variable of the voltage of the fh generator.

The computer readable medium of any of claims 41-43, wherein the instructions further cause the apparatus to sample by a multiplication of the voltage and current discrete signals with a window function.

The computer readable medium of claim 44, wherein the instructions cause the apparatus to multiply with the window function according to Z(X) =

(λ + ), where Z is a DFT of a product of Y and h, N denotes total number of samples for fining the DFT, k denotes a discrete sample, Y denotes a sinusoidal signal with harmonics and noise, h is a window function, W is a DFT of the window function, f is the frequency of Y s fundamental component in Hz, fs is a sampling frequency for DFT in Hz, and λ e {0, 1 (N - 1)}.

The computer readable medium according to any of claims 41-45, wherein the subset of the estimated magnitude, phase and frequency variables of the voltage and current is the estimated magnitude of the voltage, the estimated frequency of the voltage and the estimated magnitude of the current.

47. The computer readable medium according to claim 46, wherein the estimated magnitude of the voltage and the estimated frequency of the voltage and associated voltage variances are inputs to the state estimator, and are utilized as pseudo-inputs, and are provided in differential equations of the state estimator, and the estimated magnitude of the current and the calculated power and associated current and power variances are provided as an algebraic equation and are given as measurement inputs to the state estimator.

48. The computer readable medium according to claim 47, wherein the differential

equations are represented as xlk = xlJc + T0gl(xlk , ύι1 , wlJc) + vlJc = xlk = gl {^xlk , ulk , wlk) + vlJc, where x is a column vector of a state, k = (k - 1) , To is a sampling period in seconds, g is a discrete form of a differential function, u is a column vector of pseudo-inputs, w is a column vector of noise, /' is the ith generation unit, and k is the kth sam le and the al ebraic e uation is represented as ylk =

, where y is a column

vector of a measurement, Pe denotes active electrical-power output of a machine, / and lm represents an analogue stator current and its magnitude, respectively, Vd and Vq represents d-axis and q-axis stator voltages, respectively, Id and lq represents d- axis and q-axis stator currents, respectively.

49. The computer readable medium according to any of claims 41 -45, wherein the

subset of the estimated magnitude, phase and frequency variables of the voltage and current is the estimated magnitude of the current, the estimated frequency of the current and the estimated magnitude of the voltage.

50. The computer readable medium according to claim 49, wherein the estimated

magnitude of the current and the estimated frequency of the current and associated current variances are inputs to the state estimator and are used in differential equations of the state estimator, and the estimated magnitude of the voltage and the calculated power and associated voltage and power variances are represented as an algebraic equation and is given as a measurement input to the state estimator.

The computer readable medium according to any of claims 41 -45, wherein the subset of the estimated magnitude, phase and frequency variables of the voltage and current is the estimated magnitude of the voltage, the estimated magnitude of the current and the estimated frequency of the voltage.

52. The computer readable medium according to claim 51 , wherein the estimated

magnitude of the voltage and the estimated magnitude of the current and associated voltage and current variances are inputs to the state estimator and are used in differential equations of the state estimator, and the estimated frequency of the voltage and the calculated power and associated voltage and power variances are represented as an algebraic equation and is given as a measurement input to the state estimator, wherein the calculated power is reactive power.

53. The computer readable medium according to any of claims 41 -45, wherein the

subset of the estimated magnitude, phase and frequency variables of the voltage and current is the estimated magnitude of the current, the estimated frequency of the current and the estimated magnitude of the voltage.

54. The computer readable medium according to claim 53, wherein the estimated

magnitude of the current and the estimated frequency of the current and associated current variances are inputs to the state estimator and are used in differential equations of the state estimator, and the estimated magnitude of the voltage and the calculated power and associated voltage and power variances are represented as an algebraic equation and is given as a measurement input to the state estimator.

55. The computer readable medium according to any of claims 41 -54, wherein the

calculated power is real power or reactive power of the generator.

56. The computer readable medium according to any of claims 41 -55, wherein the

instructions further cause the apparatus to estimate the magnitude, phase and frequency variables of the voltage and current, as well as associated variance for lz0 + 2Z1+9Z2. Z0 _ e&B+e~fic β _ 1- each variable according to / ; C =

^ Z0 -2Z1 +Z2 ' Z ~ eJdE+e-JdF ' ~ JN

here Y denotes a sinusoidal signal with harmonics and noise, h is a window function, Z is a DFT of the product of Y and h, f is a frequency of Vs fundamental component in Hz, Ym is a magnitude of Y's fundamental component, fs is a sampling frequency for DFT in Hz, N is a total number of samples for finding DFT, CRB is a Cramer-Rao bound, σ denotes standard deviation with o2 as its variance, Θ denotes a phase in

Vs fundamental component, and fo denotes a voltage base value in Hz.

57. The computer readable medium according to any of claims 41-56, wherein the

instructions further cause the apparatus to calculate the power variable and the associated variance according to eife = 9^ 1 cos(§yk - 9\k , <¾fc = [σ2ί¾2 +

(i¾fe)2<¾] cos2ik - 9 k) + (i¾fe)¾)2 [<¾* + sin2fe - 9 k), where Pe is active electrical power output of the generator, Vm is voltage magnitude, \m is current magnitude, θν, and 0/ are phases of stator voltage and current, respectively, σ denotes standard deviation and o2 its variance.

58. The computer readable medium according to any of claims 41-57, wherein the

instructions further cause the apparatus to estimate the DSE according to

electrical power output of the generator, Vm is voltage magnitude, \m is current magnitude, θν, and θι are phases of stator voltage and current, respectively, σ denotes standard deviation and o2 denotes variance, ω and ω' represent column vectors of noise in y and u, respectively, P denotes a covariance matrix, f denotes frequency of Y's fundamental component in Hz, u and v represent column vectors of pseudo-inputs and process noise, respectively, represents a composite state vector, g denotes a discrete form of differential functions, / denotes a column vector of measurements, /' denotes the ith generator unit, and k denotes a discrete sample.

The computer readable medium according to any of claims 41-58, wherein the DFT is an interpolated DFT.

60. The computer readable medium according to any of claims 41-59, wherein the state estimator is an Unscented Kalman filter, UKF. 61. A power generation system comprising a plurality of power generators and an

apparatus for Dynamic State Estimation, DSE, of an operational state of a generator in the power system, wherein in the apparatus is an apparatus according to any one of claims 21-40.

Description:
DYNAMIC STATE ESTIMATION OF AN OPERATIONAL STATE OF A GENERATOR IN A POWER SYSTEM

[0001] Example embodiments described herein are directed towards dynamic state estimation of an operating state of a generator in a power system. Such estimation is performed for an individual generator in real time with the use of calculated variances, thereby providing improved estimation accuracy. According to some embodiments, the estimation may utilize the calculation of a relative angle of an individual generator thereby providing estimation without the use of Global Position System (GPS) synchronization. BACKGROUND

[0002] A disturbance in a power system (such as a fault) can initiate spontaneous oscillations in the power-flows in transmission lines. These oscillations grow in magnitude within few seconds if they are undamped or poorly damped. This may lead to loss in synchronism of generators or voltage collapse, ultimately resulting in wide-scale blackouts. The power blackout of August 10, 1996 in the Western Electricity Co-ordination Council region is a famous example of blackouts caused by such oscillations.

[0003] A generator's voltage, current and power are sinusoidal quantities, and since each sinusoid has a magnitude and a phase (which are together known as a phasor), these quantities can either be represented as sine waves, or as phasors. The conversion of sine waves to phasors is done by phasor measurement units (PMUs). In order to estimate the rotor angle any Dynamic State Estimation (DSE) algorithm available in modern power systems typically utilize synchronized measurements obtained using PMUs.

SUM MARY

[0004] One problem with current state estimation methods is the level of error in the estimations. A further problem with current state estimation is time synchronization is that it has associated noise and synchronization-errors. Synchronization errors increase the total vector error (TVE) of PMU measurements. As synchronized measurements are used for DSE, these errors can get propagated to the estimated states and deteriorate the overall accuracy and robustness of estimation. It is also not possible to completely eliminate time synchronization as it is inherently required for estimation of rotor angles.

[0005] Example embodiments directed herein provide for more accurate estimations as a greater number of inputs and measurements, as well as variances for such inputs and measures, are used in the estimation. Furthermore, some of the example embodiments presented herein provide for a means of decentralized DSE where the need for time synchronization may be eliminated. [0006] Accordingly, the example embodiments presented herein are directed towards an apparatus, computer readable medium and corresponding method for Dynamic State Estimation (DSE) of an operational state of a generator in a power system. The apparatus comprises a transceiver to receive measured voltage and current analog signals associated with the generator. The apparatus further comprises a processor to sample the received measured signals to voltage and current discrete signals. The processor is further to estimate magnitude, phase and frequency variables of the voltage and current, as well as associated voltage and current variances, respectively, for each variable, using the discrete voltage and current signals as an input in a Discrete Fourier Transform (DFT). The processor is also to calculate a power variable and associated power variance based on the estimated magnitude, phase and frequency variables for the voltage and current, as well as associated voltage and current variances for each variable. The processor is further to estimate the dynamic state of the generator using at least a subset of the estimated magnitude, phase and frequency variables of the voltage and current, as well as associated voltage and current variances for each variable, and the calculated power variable and the associate power variance using a state estimator.

[0007] The example embodiments described herein provide a means of DSE with improved accuracy as a greater number of measurements and inputs are used in conjunction with the DFT and state estimator. Furthermore, the example embodiments presented herein improve the accuracy of such estimates with the use of calculated variances.

[0008] According to some of the example embodiments, the estimation of the dynamic state of the generator may further comprise estimating a relative angle as a difference between a rotor angle and an estimated voltage phase. The processor may provide such estimation.

[0009] According to some of the example embodiments, with the use of the estimation of a relative angle, time synchronization is not needed. Therefore, DSE may be provided in a decentralized manner where the dynamic states may be utilized for decentralized control of a particular generator.

BRIEF DESCRIPTION OF THE DRAWINGS

[0010] The foregoing will be apparent from the following more particular description of the examples provided herein, as illustrated in the accompanying drawings in which like reference characters refer to the same parts throughout the different views. The drawings are not necessarily to scale, emphasis instead being placed upon illustrating the examples provided herein: Figure 1 is an illustration of a working example featuring a 16-machine, 68-bus power system model, according to some of the example embodiments described herein;

Figure 2 is a graphical comparison of DSE for α 13 , ω 13 , and E q ' 13 for a base case, according to some of the example embodiments presented herein;

Figure 3 is a graphical comparison of estimation errors for α 13 , ω 13 , and E q ' 13 for the base case, according to some of the example embodiments presented herein;

Figure 4 is a graphical comparison of DSE for ¾ 1 , ψ and V r 13 for the base case, according to some of the example embodiments presented herein;

Figure 5 is a graphical comparison of estimation errors for

¾ 1 , ψ and V r 13 for the base case, according to some of the example embodiments presented herein;

Figure 6 is a graphical comparison of DSE for ω 13 for varying noise levels, according to some of the example embodiments presented herein;

Figure 7 are tables featuring root mean square errors for DSE and DSE with PMU, according to some of the example embodiments;

Figure 8 is an operational node diagram of an apparatus for DSE of a generator, according to some of the example embodiments presented herein;

Figure 9 is an example node configuration of the apparatus described in Figure 8, according to some of the example embodiments described herein; and

Figure 10 is a flow diagram depicting example operations which may be taken by the apparatus of Figures 8 and 9, according to some of the example embodiments described herein.

NOMENCLATURE

[0011] The following sections provides a listing of the mathematical nomenclature, and their corresponding meaning, as used throughout this document:

a difference of rotor angle and stator voltage phase in rad

0 denotes a zero matrix (or vector) of appropriate size

χ denotes a state sigma point

y denotes a measurement sigma point

g, g discrete and continuous forms of differential functions, respectively

h a column vector of the system algebraic functions denotes an identity matrix of appropriate size

the Kalman gain matrix

denotes a covariance matrix or cross-covariance matrix column vectors of pseudo-inputs and process noise, respectively column vectors of noise in y and u , respectively

column vectors of states and measurements, respectively composite state vector

rotor angle in rad

denote estimated and predicted values, respectively

denotes the A th component of a DFT

rotor-speed and its synchronous value in rad/s, respectively sub-transient emfs due to d axis damper coil in p.u.

sub-transient emfs due to q axis damper coil in p.u. denotes standard deviation, with o 2 as variance

phase of Vs fundamental component in rad

phases of stator voltage and current, respectively, in rad rotor damping constant in p.u.

transient emf due to flux in q-ax\s damper coil in p.u.

transient emf due to field flux linkages in p.u.

field excitation voltage in p.u.

frequency of Vs fundamental component in Hz

sampling frequency for interpolated-DFT method in Hz frequency of V n p.u. and its base value in Hz, respectively generator inertia constant in s

Hann window function

analogue stator current and its magnitude, respectively, in p.u. denote the f h generation unit and V-ϊ, respectively

d-axis and q-ax\s stator currents, respectively, in p.u.

A*" and ( -1)th samples and the I th sigma-point, respectively

AVR gain in p.u.

the ratio (X" d - X ) / (X d ' - X ) K d2 the ratio (¾ - X" d )/(X d ' - X t )

K q1 the ratio

K q2 the ratio {x q ' - X" q )/{x q ' - m, n number of states in x and X, respectively, (n=m+2)

N, M total samples for finding DFT and total generation units, respectively

P e active electrical-power output of a machine in p.u.

Rs armature resistance in p.u.

t system time in s

T, To denote transpose and UKF's sampling period (in s), respectively Te, TmO electrical and mechanical torques, respectively, in p.u.

T r time constant for the AVR filter in s

T d ' 0 , T q ' Q d-axis and q-ax\s transient time constants, respectively, in s

T d ' 0 , T q ' 0 d-axis and q-ax\s sub-transient time constants, respectively, in s

V, Vm analogue stator voltage and its magnitude, respectively, in p.u. Vd, V q d-axis and q-ax s stator voltages, respectively, in p.u.

V r , Vre f AVR-filter voltage and AVR-reference voltage, respectively, in p.u.

W DFT of Hann window function

Xd, X q d-axis and q-ax s synchronous reactances, respectively, in p.u.

X d ' 0 , X q ' 0 d-axis and q-ax s transient reactances, respectively, in p.u.

X d o, X q o d-axis and q-ax\s sub-transient reactances, respectively, in p.u.

X armature leakage reactance in p.u.

Y denotes a sinusoidal signal with harmonics and noise

Ym magnitude of Vs fundamental component in p.u.

Z DFT of the product of Y and h

DETAILED DESCRIPTION

[0012] In the following description, for the purposes of explanation and not limitation, specific details are set forth, such as particular components, elements, techniques, etc. in order to provide a thorough understanding of the examples provided herein. However, the examples may be practiced in other manners that depart from these specific details. In other instances, detailed descriptions of well-known methods and elements are omitted so as not to obscure the description of the examples provided herein.

[0013] Example embodiments described herein are directed towards dynamic state estimation of an operating state of a generator in a power system. A disturbance in a power system (such as a fault) can initiate spontaneous oscillations in the power-flows in transmission lines. In order to monitor and control such oscillations and related dynamics which cause instability, the operating state of the system needs to be estimated in realtime, with update rates which are in time scales of ten milliseconds or less (as the time constants associated with such oscillations are not more than ten milliseconds), and this real-time estimation of operating state is known as dynamic state estimation (DSE).

[0014] The dynamic states which are estimated and obtained as outputs from DSE algorithms are angles, speeds, voltages and fluxes of the rotors of all the generators in the power system. The inputs which are given to DSE algorithms are some measurable time- varying quantities such as voltage and current of the stator, and some measurable time- invariant quantities such as resistances, reactances, inertia and other constants for the generator. The constant quantities are measured beforehand, and are used as parameters in DSE algorithms.

[0015] A generator's voltage, current and power are sinusoidal quantities, and since each sinusoid has a magnitude and a phase (which are together known as a phasor), these quantities can either be represented as sine waves, or as phasors. The conversion of sine waves to phasors is done by phasor measurement units (PMUs). In order to estimate the rotor angle any Dynamic State Estimation (DSE) algorithm available in power system current systems typically utilize synchronized measurements obtained using PMUs.

[0016] Example embodiments presented herein provide a means of improved DSE where greater number of inputs and measurements, as well as variances for such inputs and measurements, are utilized in the estimation. Thus, providing improved accuracy for the measurements. According to some of the example embodiments, a means for DSE is provided without the use of time synchronization.

[0017] One problem with time synchronization is that it has associated noise and synchronization-errors. Synchronization errors increase the total vector error (TVE) of PMU measurements. As synchronized measurements are used for DSE, these errors can get propagated to the estimated states and deteriorate the overall accuracy and robustness of estimation. It is also not possible to completely eliminate time

synchronization as it is inherently required for estimation of rotor angles.

[0018] According to some of the example embodiments, it has been appreciated that although time synchronization is typically used for the estimation of the rotor angle, it is not needed for estimation of other dynamic states, such as rotor speed, rotor voltages and fluxes, as these states are not defined with respect to a common reference angle. Thus, if the dynamic model which is used for estimation can be modified in such a way that rotor angle is replaced with another angle which does not require time-synchronization, then this can minimize the effects of synchronization on accuracy and robustness of estimation.

[0019] Some of the example embodiments provide an algorithm for DSE which realizes the above concept. This is done by modifying the estimation model to estimate a relative angle (which does not require synchronization) instead of rotor angle. One such angle is the difference between the rotor angle and the generator terminal voltage phase, also known as the internal angle of the generator. As the rotor angle and the voltage phase have a common reference angle, this reference angle gets cancelled in the difference of the two quantities. Thus, the internal angle, rotor speed, voltages and fluxes can be estimated using the modified estimation model without requiring any synchronized measurements. These dynamic states can then be utilized for decentralized control of the generator. It should be noted that, according to some of the example embodiments, if the estimation of rotor angle is specifically required then it can be indirectly estimated as the sum of the estimated internal angle and the measured terminal voltage phase obtained using PMU.

[0020] An example advantage of some of the example embodiments provided herein is to provide a system and method of DSE in which dynamic states are estimated without any time synchronization by incorporating internal angle in estimation model, which in turn ensures robustness of the method to synchronization errors.

[0021] The error in phasor measurements considered in several existing methods of DSE is much less than 1 % TVE. Such methods of DSE do not consider realistic errors in measurements. The example embodiments presented herein considers and remains accurate for varying levels of errors in measurements - from 0.1 % to 10%. Furthermore, none of the currently available methods take into account GPS synchronization errors.

[0022] As synchronization is not required for estimation of the states, according to some of the example embodiments, DSE for these states can be performed using the analogue measurements directly acquired from current transformers (CTs) and voltage transformers (VTs). This is particularly beneficial for decentralized control purposes.

[0023] According to some of the example embodiments, a dual-stage estimation process has been proposed in which a Discrete Fourier transform (DFT) and a state estimator. According to some of the example embodiments, the DFT may be an interpolated DFT and the state estimator may be an unscented Kalman filtering (UKF). According to some of the example embodiments the DFT and the state estimator have been combined as two stages of estimation. The DFT stage dynamically provides estimates of means and variances of the inputs required by the state estimation stage, and this continuous updating of variances may provide noise-robustness of the proposed example

embodiments. In existing methods of DSE for power systems, only static estimates of measurement variances are provided to the estimator.

[0024] According to some of the example embodiments, analytical expressions have been obtained for the means and variances of the parameter estimates of a sinusoidal signal (which are given as input to the state estimation stage from the DFT stage). Most of these expressions are currently not available in literature.

[0025] Rest of the description is organized as follows. First, decoupled equations which may be utilized in conjunction with the example embodiments presented herein are discussed under the subheading 'Power System Dynamics in a Decoupled Form'.

Thereafter, the process for estimation of magnitude, phase and frequency of the analogue signals of terminal voltage and current is discussed under the subheading 'DFT based Estimation'. An explanation as to how these estimates may be further used for DSE using a state estimator is provided under the subheading 'State Estimation'. A working example featuring simulations to demonstrate the development of the example embodiments presented herein is provided under the subheading 'Case Study. An example node configuration of an apparatus that may be utilized in providing DSE is discussed under the subheading 'Example Node Configuration' . Finally, an operational flow is provided under the 'Example Operations'.

[0026] Power System Dynamics in a Decoupled Form

[0027] A power system comprises a wide variety of elements, including generators, their controllers, transmission lines, transformers, relays and loads. All these elements are electrically coupled to each other, and, therefore, in order to define a power system using dynamic mathematical equations, knowledge of the models, states and parameters of all these constituent elements is useful. Acquiring this knowledge in real-time is not feasible as power systems span wide geographic regions, which are as large as a country, or even a continent. Therefore, according to some of the example embodiments, dynamic equations of the power system in a decoupled form is utilized, so that the real-time estimation of dynamic states can be conducted. According to some of the example embodiments, the estimation may be performed in a decentralized manner.

[0028] Such a decoupling of system equations may be achieved if a generator and its controller(s) is considered as a decentralized unit, and the stator terminal voltage magnitude, Vm, and its phase, 0V , are treated as 'inputs' in the dynamic equations, instead of considering them as algebraic quantities or measurements. This concept is referred to herein as 'pseudo-inputs' for decoupling the equations.

[0029] In order to estimate the internal angle, which is the difference between the rotor angle and the voltage phase, instead of estimating the rotor angle, the decoupled equations and the pseudo-inputs for a generator are altered. The altered decoupled equations are given by equations (1)-(1 1), derived using the sub-transient model of machines with four rotor coils in each machine, known as IEEE Model 2.2 as provided in "IEEE Guide for Synchronous Generator Modeling Practices and Applications in Power System Stability Analyses," IEEE Std 1 110-2002 (Revision of IEEE Std 11 10-1991), pp. 1- 72, 2003. In these equations, the altered pseudo-inputs are Vm and voltage frequency, fV , and i refers to the system's ith machine or generator, 1≤ i≤ M. Slow dynamics of the speed-governor have been ignored in this model (although they can also be added, if required). Also, model of a static automatic voltage regulator (AVR) is included with the model of each machine.

Αά ι = ( - f ) (1) cd - E d - i X q - Xq ) K ql I q + K q2 j j (3)

E d ~ E q + ( X d ~X d )

E d n - (4)

1 qo

¾ = ψη [Ε + {X a d - XiYd - d] (5)

'do

= ψ 1 qo η [-Ej' + ( - xiYd - ¾] (6) where, (7)

'r

Efd = Ka [Vref ~ Vr ]> Efdmin— Efd— Efdmax (8)

& = ^d 2 + Pq 2 . K = - sin a 1 , V q l = l¾ cos a 1 (1 1)

[0030] The above equations may be written in the following composite state-space form which may be used for DSE (here pseudo-inputs are denoted by u ', and the process noise and the noise in pseudo inputs have been included, and denoted by v' and ω',

respectively). x l = g n (x u a , w n ) + v i ; u' i - w n = [l¾ f ] T (12)

[0031] DFT based Estimation

[0032] The example embodiments described herein make use of an interpolated DFT for estimating the parameters of a sinusoidal signal. It should be appreciated that other forms of DFTs may be utilized in the estimation of the parameters. The use of an interpolated DFT based estimation has the example advantage of being both fast and accurate enough for real-time control applications in power systems. According to some of the example embodiments, the DFT may be used for finding the estimates of frequency, magnitude and phase of the fundamental components of measurements obtained from CTs and PTs.

[0033] The fundamental component of a sinusoidal signal can be extracted by multiplying the signal with a suitable window function which eliminates other harmonics and higher frequency components in the signal, followed by finding its DFT. One such function is, for example, a Hanning window function given by hk = sin2( nk/N ), and if this function is multiplied with N samples of an analogue signal Y(t) sampled at a frequency fs, then DFT of the product is given by Ζ(λ) as follows. ZW =∑ Y k k e-^ = ^ Y e^w ( - ) - ^ Y e-^w ( + ) (13)

[0034] where, Ym, Θ and f are magnitude, phase and frequency of Y's fundamental component, respectively; λ ε {0, 1 , ...,Ν - 1}; and W(A) is the following DFT of Hanning window function.

WW (14) [0035] A concept in interpolated-DFT based estimation is to approximate W(A) with the following expression, provided that N » 1 and λ « N.

[0036] By substituting equation (15) in equation (13), Ζ(λ) can be expressed as follows for N » 1 and λ « N.

[0037] where Y m , Θ and / denote the estimates of Y m , Θ and f, respectively. As equation (16) has three unknowns (which are Y m , Θ and /), three distinct equations are required to estimate these unknowns. This may be done by choosing any three distinct values of λ in equation (16) (say λ = 1 , λ = 2 and λ = 3). The obtained values of Y m , Θ and / will have associated estimation errors which will depend on N and on the values of λ which are used for generating the three distinct equations. More precisely, these estimation errors are inversely proportional to N 4 , and, hence, N should be as large as practically feasible. In the example embodiments described herein N is taken to be in the order of 10 3 , as this is the highest order for N for which interpolated-DFT can run on a state-of-the-art DSP processor without overloading it, however it should be appreciated such values are presented herein merely as an example. Overloading refers to overall processor usage of above 95%. Also, for a given N, the estimation errors are minimized if the choices for λ are taken as λ = 0, λ = 1 and λ = 2, provided that— < 2.1 ; otherwise, for 2.1 <— < 3, the errors are minimized if fs fs

the choices are λ = 1 , λ = 2 and λ = 3. The value of— should not be greater than 3 as then fs

the delay in obtaining the estimated values becomes too large (that is, more than two cycles, or more than 0.04 s for a 50Hz power system), and at the same time it should not be too small as then the accuracy of estimation is diminished. According to some of the example embodiments, an intermediate value of— ~ 1 .5 has been taken and, hence, the fs

former choices of λ = 0, λ = 1 and λ = 2 are applicable.

[0038] — is an unknown quantity as f needs to be estimated. But because of power fs

system operational requirements, f should remain within 5% of the base system frequency, fo (which is usually 50 Hz or 60 Hz), and, hence, if N and f s are chosen such that— = 1.5, fs then ^ « 1.5.

fs

[0039] The 3 equations which are obtained by putting λ = 0, λ = 1 and λ = 2 in equation (16) can be written in matrix form as follows.

[0040] Equation (17) implies that the product of a square-matrix and a column vector is equal to a zero vector, when both the matrix and the vector have non-zero elements. This occurs if the columns of the matrix are linearly dependent, that is, the determinant of the matrix is zero, given as follows.

[0041] Simplification of the above determinant gives / as follows. f fs IZ 0 + 2Z 1 + 9Z 2

(19) ' N J Z 0 -2Z 1 + Z 2

[0042] Θ may be obtained by substituting the above value of / back into equation (16) and eliminating Y m . To do this, the equation which is obtained by putting λ = 0 in equation (16) is divided by the equation obtained by putting λ = 1 in equation (16), which comes as follows.

[0043] Solving for e ie using equation (20) gives the following expression.

[0044] Using equations (21) and (16) (with λ = 0), Y m comes as follows.

[0045] where B, C, and e ie are given by equations (20)-(21).

[0046] It should be noted that Y m , Θ and / are real quantities, but they are obtained as functions of complex quantities (given in the right hand sides (RHSs) of equations (19), (21) and (22), respectively). Hence, these quantities will have negligible but finite imaginary parts associated with them because of finite computational accuracy of any computational device. Thus, during implementation, the imaginary parts should be ignored and only the real parts of RHSs should be assigned to Y m , Θ and /. Also, as Y m and / are strictly positive, absolute values of real parts of respective RHSs should be assigned to them.

[0047] The variance of the above estimate of / in equation (19) is approximately twice the minimum possible variance which is theoretically achievable using any unbiased estimator (known as Cramer-Rao bound (CRB)). CRB for frequency estimation of a sinusoidal signal and is given by CRB(f) (in Hz 2 ) as follows.

[0048] where σγ is the variance of noise in Y (in p.u.). CRBs for Y m and Θ are given by CRB(? m ) (in p.u.) and CRB(0) (in rad 2 ), respectively, as follows.

[0049] The variances of Y m and Θ are found to be approximately two and six times the above CRBs in equation (24), respectively; and hence, the estimated variances of /, Θ and Y m are given by dj (in p.u.), dy m (in p.u.) and 8g (in rad 2 ), respectively, as follows. d2 = 2€ ) . ^ = 2Cff g( ?m ). d 2 = 6CffS (£) (25)

[0050] where CRB(/), CRB(? m ) and CRB(0) are given by equations (23)-(24). Estimates of means and variances obtained above are given as inputs to the state estimation stage (e.g., UKF), as detailed in the next section.

[0051] An example advantage of obtaining the analytical expressions for dj , 8y m , and 8g in equations (23)-(25) is that these variances may be continuously updated and provided to the dynamic estimator (which is the state estimation stage, for example, the UKF stage) along with /, Θ and Y m , thereby improving the accuracy of dynamic state estimation. [0052] State Estimation

[0053] State estimation will be described herein using an UKF as an example. However, it should be appreciated that other forms of state estimation, or types of Kalman Filters, may be employed. UKF is a nonlinear method for obtaining dynamic state estimates of a system. It employs the idea that performing DSE is easier if the distribution of state estimates is transformed, than if the system model itself is transformed through

linearization. System linearization requires computation of Jacobian matrices and is a mathematically challenging task for a high order power system model, especially if it needs to be done at every iteration. Since linearization is not required in UKF, and, moreover, it has higher accuracy and similar computational speeds as that of linear methods of DSE, UKF has been used for performing DSE according to some of the example embodiments, however, other Kalman Filters or forms of state estimation may also be employed. UKF is a discrete method and, hence, the system given by (12) may be discretized before UKF may be applied to it. Discretizing (12) at a sampling period To, by approximating x l with [x lk - x lk )/T 0 , gives the following equation (where k and k represent the kth and (k - 1) th samples, respectively). xik = x Ck + T 0 g i (x ik , u ik , w ik ) + v ik => x ik = g^x^. u^. w ) + v ik (26)

[0054] In the above model, i¾ fe and f v k (found using the DFT method) are used in the pseudo-input vector ii lk as follows. u ik = [V fr k Y = [V fj k ] T + w ik (27)

[0055] UKF also utilizes a measurement model besides the above process model. The estimates of active power, e ife (defined by equation (10)), and stator current magnitude, /, (defined by equation (1 1)), which are obtained using the DFT method are used as measurements for UKF. After incorporating the measurement noise, w lk , the

measurement model is given as follows. _ h l {x lk , U lk , W lk ) + W (28)

[0056] u lk and y lk are estimated quantities and have finite variances which may be included in the process and measurement models, respectively. This may be done by including w lk and w lk in the models as the following zero-mean noises.

[0057] where ¾ and P„ denote the covariance matrices of w lk and w , respectively. In order to find the estimates and the variances in equations (27)-(30), the stator voltage, V'(t), and stator current, l'(t), measured using VT and CT, respectively, are processed using the DFT method. Thus, V m k , f v k , θ , d k , σ , and <¾ fc are obtained by putting

vm fv "v

Y(t)=V'(t) in equations (13)-(22) and updating these estimates and variances for every k th sample. Similarly, 1%, f, lk , 6\ k , σ , σ , and <¾ fc are obtained by putting Y(t)=l'(t). As

Pe k = V m k Im cos(§y k - 6\ k ) (from equation (10)) and the means values and variances of Vm , Im , 6y k , and 6\ k are known, the mean value of e ife (denoted as e ife ) and its estimated variance (denoted as <¾ fc can be represented in terms of these known quantities, and have been obtained as follows (here it should be noted that by definition θ and 6\ k lie in the interval (-ττ/2, ττ/2], hence, they should be 'unwrapped' by adding or subtracting suitable multiples of π to them, in order to find cos{o v k - 9\ k )).

yik jik

COS (31)

+(t¾ fe ) 2 ( fe ) 2 [¾ + ¾ i ^in 2 (^

[0058] Thus, the four quantities which are utilized by the state estimation stage, for example the UKF stage, from the DFT stage are it 1 , y l , P , and Pj , given by equations (27)-(30). These quantities should be updated ever To s, as this is the sampling period of the UKF stage. Also, in equation (26), both x lk and w lk are unknown quantities and may be combined together as a composite state vector X lk with a composite covariance matrix Ρψ defined as follows.

[0059] Here P lk is the covariance matrix of x lk , and ¾ is the cross-covariance matrix of x lk and w lk . With the above definition, the model in equations (26)-(28) is redefined as follows.

X ik = g i ( x Ck > i- + v Ck . yik = h i ( X ik > 6 ifc) + w ik (33)

[0060] With equation (33) as a model and χ ί0 as a steady state estimate of x lk and with the knowledge of g l , h l , u l , y l , P , Pj and the process noise covariance matrix, v ife , the filtering equations of the state estimator, for example, the UKF, for the k th iteration of the i th unit are given as follows.

Operation 1 : Initialize

if (k==1) then initialize x ik = x i0 , w' i = 0 2xl , P ik = p{° P = 0 m½2 , P ik = pj° in

equation (32) to get Ρψ & X ik .

else reinitialize w' ik and in equation (32) according to equation (29), leaving the rest of the elements in X ik and pf unchanged.

Operation 2: Generate sigma points

Operation 3: Predict states

Operation 4: Predict measurements

Operation 5: State estimation Update (e.g., Kalman Update)

K lk = X ik _ X ik- + K ik(yik _

Operation 6: Output and time update

output X ik and Ρ , k ^ (k + 1), go to Operation 1. [0061] Case Study

[0062] Figure 1 illustrates a model 16-machine, 68-bus benchmark test system that has been used for the case study presented herein as an example. The test system has been utilized with MATLAB-Simulink running on Windows 7 has been used for its modelling and simulation. A detailed description of the system (including various parameters) is given in B. Pal and B. Chaudhuri, Robust Control in Power Systems. New York, U.S.A.: Springer, 2005; and A.K. Singh, B.C. Pal, "Report on the 68-bus, 16-machine, 5-area system," IEEE PES Task Force on Benchmark Systems for Stability Controls, version 3.3, pp. 1-41 , Dec. 2013. Static AVRs are used in all the machines, and their parameters are given in A.K. Singh, B.C. Pal, "Decentralized Control of Oscillatory Dynamics in Power Systems Using an Extended LQR," IEEE Trans. Power Syst, vol. 31 , no. 3, pp. 1715-1728, May 2016.

[0063] According to some of the example embodiments, the robust dynamic state estimator (as discussed under the subheadings 'Power System Dynamics in a Decoupled Form' and 'State Estimation') runs at the location of each generation unit, and provides dynamic state estimates for the unit. The measurements which are provided to the estimator are V (t) and l(t), and are generated by adding noise to the simulated analogue values of terminal voltage and current of the unit. As explained under the subheading 'DFT based Estimation', N, fo and f s are taken as 1200, 50 Hz and 40000 Hz, respectively. The sampling period of state estimation stage (e.g., UKF stage), To, is taken as 0.01s, and thus, the estimates obtained from the DFT stage are also updated every 0.01 s. Also, Pv l is found. For comparison with the proposed estimator, another UKF based dynamic state estimator which uses PMU measurements also runs at each unit's location and is termed as DSE-with-PMU. Estimate of the internal angle in case of DSE-with-PMU is obtained by subtracting the measurement of terminal voltage phase from the estimate of rotor angle.

[0064] The measurement error for the robust DSE method is the percentage error in the analogue signals of V(t) and l(t), while the measurement error for DSE-with-PMU method is the total vector error (TVE) in the phasor measurements of terminal voltage and current. As the measurement errors for the two estimators are of two different kinds, these methods may not be directly compared for the same noise levels. Nevertheless, performance of the two methods for standard measurement errors can be compared, as specified by IEEE. As mentioned in these standards, the measurement error in CTs/VTs should be less than 3%, while the standard error for PMUs is 1 % TVE. Hence, in the base case for comparison, the measurement error for robust DSE is taken as 3%, while for the DSE-with-PMU method, it is taken as 1 % TVE. The system starts from a steady state in the simulation. Then at t = 1s, a disturbance is created by a three-phase fault at bus 54 and is cleared after 0.18s by opening of one of the tie-lines between buses 53-54. The simulated states, along with their estimated values for the base case for one of the units (the 13 th unit), have been plotted in Figure 2 and Figure 4. Corresponding estimation errors, which is the difference of estimated and simulated values, have also been plotted in Figure 3 and Figure 5.

[0065] It can be seen in Figures 2-5 that for robust DSE the plots of estimated values almost coincide with those of the simulated values and the estimation errors are low, but for the DSE-with- PMU method the difference between the simulated and estimated values is apparent and the estimation errors are much higher. This shows that the example embodiments perform accurately with standard measurement errors in CTs/VTs, while DSE-with-PMU fails to do so with standard errors in PMU measurements.

[0066] Robustness of the example embodiments has been tested against varying noise levels in measurements. Figure 6 shows the estimation results for ω 13 for two more cases: in the first case the noise levels are one-third the base case, while in the second case the noise levels are thrice the base case. Also, root mean squared errors (RMSEs) for varying error levels have been calculated and tabulated in Table l-Table II, which are illustrated in Figure 7, for the two methods. It can be observed that the performance of the proposed method remains robust to errors up to 3%, and even for 10% error, its performance deteriorates only to a small extent. On the other hand, DSE-with-PMU method does not perform accurately for error levels above 0.3% TVE, that is, it is not accurate for 1 % TVE and 3% TVE. It should be noted that various TVE levels have been generated by varying the synchronization errors in the simulated PMU measurements.

[0067] Computational feasibility of the example embodiments may be inferred from the fact that the entire simulation, including simulation of the power system, with two estimators at each machine, runs in real-time on MATLAB-Simulink running on Windows 7 on a personal computer with Intel Core 2 Duo, 2.0 GHz CPU and 2 GB RAM. The expression 'real-time' here means that 1 second of the simulation takes less than 1 second of processing time. Also, the total execution time for all the operations for the proposed method for one time step (that is for one iteration) is 0.44 millisecond. Thus, the method can be easily implemented using current technologies as the update rate required by the proposed method is 10 milliseconds.

[0068] Example Node Configuration [0069] Figure 8 illustrates a functional example node configuration of an apparatus 10 for DSE of an operational state of a generator in a power system, as described herein. The apparatus 10 may be configured to receive measured voltage (V) and current (I) signals associated with the generator. Upon undergoing signal processing, the signals are input in a Discrete Fourier Transform (DFT). The DFT is used to estimate magnitude, phase and frequency variables of the voltage (V m , θ ν , and fv) and current (l m , θ,, and fi), as well as associated voltage and current variances, respectively, for each variable.

[0070] The apparatus may be further configured to calculate a power variable (P e ), and corresponding power variance, based on the estimated magnitude, phase and frequency variables of the voltage (V m , θ ν , and fv) and current (l m , θ ν , and fi), as well as associated voltage and current variances, respectively, for each variable.

[0071] Thereafter, the dynamic state (DSE) of the generator may be estimated using at least a subset of the estimated magnitude, phase and frequency variables of the voltage (Vm, θ ν , and fv) and current (l m , θ ν , and fi), as well as associated voltage and current variances, respectively, for each variable, as well as the calculated power variable (P e ), and corresponding power variance.

[0072] In Figure 8, two examples of such subsets are provided. The first example includes a subset comprising the estimated magnitude of the voltage (V m ), the estimated frequency of the voltage (fv), the estimated magnitude of the current (l m ) and the calculated power (P e ) variable as well as all associated variances. The second example includes a subset comprising the estimated magnitude of the current (l m ), the estimated frequency of the current (If), the estimated magnitude of the voltage (V m ) and the calculated power (P e ) variable as well as all associated variances.

[0073] The chosen subset may thereafter be input in to the state estimator, for example a UKF. The portions of the subset may be provided to the state estimator as differential equations and portions of the subset may be provided to the state estimator as an algebraic equation. The output of the state estimator is the DSE of the operational state of a generator in the power system. According to some of the example embodiments, as the DFT and the state estimator make use of estimated and calculated variances, a more accurate DSE may be provided.

[0074] Figure 9 illustrates an example hardware node configuration of the apparatus 10 of Figure 8. The apparatus 10 may comprise any number of transceiver 12 which may be configured to receive and transmit any form of measurement, instructions, processing or estimation related information as described herein. According to some of the example embodiments, the transceiver may also comprise a single transceiving interface or any number of receiving and/or transmitting interfaces.

[0075] The apparatus 10 of Figure 9 may also comprise at least one processor 14 which may be configured to process received measurement signals, estimate various variables and variances, calculate a power variable and variance and estimate the DSE of the generator. The processor 14 is further configured to provide an operation discussed herein. The processor 14 may be any suitable computation logic, for example, a microprocessor, digital signal processor (DSP), field programmable gate array (FPGA), or application specific integrated circuity (ASIC) or any other form of circuitry.

[0076] The apparatus 10 may further comprise at least one memory 16 that may be in communication with the transceiver and the processor. The memory 16 may store received or transmitted data, processed data, and/or executable program instructions. The memory may be any suitable type of machine readable medium and may be of a volatile and/or non-volatile type.

[0077] Example Node Operations

[0078] Figure 10 illustrates a flow diagram depicting example operations which may be taken by the apparatus 10 in Dynamic State Estimation (DSE) of an operational state of a generator in a power system as described herein. It should be appreciated that Figure 10 comprises some operations which are illustrated in a solid border and some operations which are illustrated with a dashed border. The operations which are comprised in a solid border are operations which are comprised in the broadest aspect. The operations which are comprised in a dashed boarder are example aspects which may be comprised in, or a part of, or are further operations which may be taken in addition to the operations of the broader example aspects. It should be appreciated that these operations need not be performed in order. Furthermore, it should be appreciated that not all the operations need to be performed. The example operations may be performed in any order and in any combination.

[0079] Operation 20

[0080] The apparatus 10 is to receive 20 measured voltage and current analog signals associated with the generator. The transceiver 12 is to receive measured voltage and current analog signals associated with the generator. [0081] Operation 22

[0082] The apparatus 10 is to sample 22 the received measurement signals to voltage and current discrete signals. The processor 14 is to sample the received measurement signals to voltage and current discrete signals.

[0083] Example operation 24

[0084] According to some of the example embodiments, the sampling 22 further comprises multiplying 24 the voltage and current discrete signals with a window function. The processor may multiple the voltage and current discrete signals with the window function.

[0085] According to some of the example embodiments, the multiplication with the window function may further comprise a multiplication according to Ζ(λ) =∑¾=ο Y k h k e N ~ =- Y ^ e' e W (λ - ψ) ~ Y je-' e W (λ + ). Example operation 24 is further described under at least the subheading 'DFT based Estimation'.

[0086] Operation 26

[0087] The apparatus is further to estimate 26 magnitude, phase and frequency variables of the voltage and current, as well as associated voltage and current variances, respectively, for each variable, using the discrete voltage and current signals as an input in a DFT. The processor 14 is to estimate the magnitude, phase and frequency variables of the voltage and current, as well as associated voltage and current variances, respectively, for each variable, using the discrete voltage and current signals as an input in a DFT.

[0088] According to some of the example embodiments, the estimation 26 may be

2CRB ^ = 2CRB(? m ); dg = 6CRB(§) An example advantage of operation 26 is that with the calculation of the variances for respective variables, DSE with greater accuracy may be achieved. Operation 26 is further described under at least the subheading FT based Estimation'.

[0089] Operation 28

[0090] The apparatus is further to calculate a power variable and associated power variance based on the estimated magnitude, phase and frequency variables of the voltage and current, as well as associated voltage and current variances for each variable. The processor 14 is to calculate the power variable and associated power variance based on the estimated magnitude, phase and frequency variables of the voltage and current, as well as associated voltage and current variances for each variable.

[0091] According to some of the example embodiments, the calculated power may be a real power or a reactive power. According to some of the example embodiments, the calculation of the power variable and the associate variance further comprises a calculation according to e ife = +

(i¾ fe ) 2 <¾] cos 2 fe - 9 k ) + (i¾ fe )¾) 2 [<¾* + sin 2 fe - 9 k ). Operation 28 is further described under at least the subheading 'State Estimation'.

[0092] Example operation 30

[0093] According to some of the example embodiments, the DFT is an interpolated DFT. Example operation 30 is further described under at least the subheading 'DFT based Estimation'.

[0094] Operation 32

[0095] The apparatus 10 is further to estimate the dynamic state of the generator using at least a subset of the estimated magnitude, phase and frequency variables of the voltage and current, as well as associated voltage and current variances for each variable, and the calculated power variable and the associated power variance using a state estimator. The processor 14 is to estimate the dynamic state of the generator using at least the subset of the estimated magnitude, phase and frequency variables of the voltage and current, as well as associated voltage and current variances for each variable, and the calculated power variable and the associated power variance using the state estimator. [0096] According to some of the example embodiments, the estimating 32 may further comprise estimating the DSE according to

= fe^) 2 + ^ fe ) 2 ¾] COs2 ^

[0097] Example operation 34

[0098] According to some of the example embodiments, the estimating 32 of the dynamic state is based on, at least in part, estimating 34 a relative angle as a difference between a rotor angle and the estimated voltage phase. The processor 14 may estimate the relative angle as the difference between the rotor angle and the estimated voltage phase.

[0099] According to some of the example embodiments, the estimating of the relative angle further comprises estimating the relative angle according to Αά ι = (ω ι - fy). An example advantage of example operation 34 is that with the use of a relative angle, time synchronization may be avoided. Example operation 34 is further described under at least the subheading 'Power System Dynamics in a Decoupled Form'.

[00100] Example operation 36

[00101] According to some of the example embodiments, the estimating 32 is based on a subset of the estimated magnitude of the voltage, the estimated frequency of the voltage and the estimated magnitude of the current. The processor 14 may base the estimation of the dynamic states on a subset comprising the estimated magnitude of the voltage, the estimated frequency of the voltage and the estimated magnitude of the current.

[00102] Example operation 38

[00103] According to the embodiments in which the subset is as described in example operation 36, the estimated magnitude of the voltage and the estimated frequency of the voltage and associated voltage variances are inputs to the state estimator, and are utilized as pseudo-inputs, and are provided in differential equations of the state estimator, and the estimated magnitude of the current and the calculated power and associated current and power variances are provided as an algebraic equation and are given as measurement inputs to the state estimator.

[00104] According to such example embodiments, the differential equations are represented as x lk = x lk + T 0 g l (x lk , ii lk , w lk ) + v lk = x lk = g l {x lk , u lk , w lk ) + v lk , where x is a column vector of a state, k = (k - 1) , To is a sampling period in seconds, g is a discrete form of a differential function, u is a column vector of pseudo-inputs, w is a column vector of noise, i is the i th generation unit, and k is the k th sample and the algebraic equation is represented as y lk =

h l (x lk , ii lk , w lk ) + w lk . Example operation 38 is further described under at least subheading 'State Estimation'.

[00105] Example operation 40

[00106] According to some of the example embodiments, the estimating 32 is based on a subset of the estimated magnitude of the current, the estimated frequency of the current and the estimated magnitude of the voltage. The processor 14 is to estimate the dynamic state based on the subset of the estimated magnitude of the current, the estimated frequency of the current and the estimated magnitude of the voltage.

[00107] Example operation 42

[00108] According to the embodiments in which the subset is as described in example operation 40, the estimated magnitude of the current and the estimated frequency of the current and associated current variances are inputs to the state estimator and are used in differential equations of the state estimator, and the estimated magnitude of the voltage and the calculated power and associated voltage and power variances are represented as an algebraic equation and is given as a measurement input to the state estimator.

[00109] It should further be appreciated that other such subsets may be utilized. For example, according to some of the example embodiments, the subset of the estimated magnitude, phase and frequency variables of the voltage and current is the estimated magnitude of the voltage, the estimated magnitude of the current and the estimated frequency of the voltage. According to such embodiments, the estimated magnitude of the voltage and the estimated magnitude of the current and associated voltage and current variances are inputs to the state estimator and are used in differential equations of the state estimator, and the estimated frequency of the voltage and the calculated power and associated voltage and power variances are represented as an algebraic equation and is given as a measurement input to the state estimator, wherein the calculated power is reactive power.

[00110] According to some of the example embodiments, the subset of the estimated magnitude, phase and frequency variables of the voltage and current is the estimated magnitude of the current, the estimated frequency of the current and the estimated magnitude of the voltage. According to such example embodiments, the estimated magnitude of the current and the estimated frequency of the current and associated current variances are inputs to the state estimator and are used in differential equations of the state estimator, and the estimated magnitude of the voltage and the calculated power and associated voltage and power variances are represented as an algebraic equation and is given as a measurement input to the state estimator. Example operation 42 is further described under at least subheading 'State Estimation'.

[00111] Example operation 44

[00112] According to some of the example embodiments, the estimating 32 further comprises estimating the DSE with the use of a UKF as the state estimator. The processor 14 may estimate the DSE with the use of the UKF as the state estimator.

[00113] The various example embodiments described herein are described in the general context of method steps or processes, which may be implemented in one aspect by a computer program product, embodied in a computer-readable medium, comprising computer-executable instructions, such as program code, executed by computers in networked environments. A computer-readable medium may comprise removable and non-removable storage devices comprising, but not limited to, Read Only Memory (ROM), Random Access Memory (RAM), compact discs (CDs), digital versatile discs (DVD), etc. Generally, program modules may comprise routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. Computer-executable instructions, associated data structures, and program modules represent examples of corresponding acts for implementing the functions described in such steps or processes. [00114] Throughout the description and claims of this specification, the words "comprise" and "contain" and variations of them mean "including but not limited to", and they are not intended to (and do not) exclude other moieties, additives, components, integers or steps. Throughout the description and claims of this specification, the singular encompasses the plural unless the context otherwise requires. In particular, where the indefinite article is used, the specification is to be understood as contemplating plurality as well as singularity, unless the context requires otherwise.

[00115] Features, integers, characteristics or groups described in conjunction with a particular aspect, embodiment or example of the invention are to be understood to be applicable to any other aspect, embodiment or example described herein unless incompatible therewith. All of the features disclosed in this specification (including any accompanying claims, abstract and drawings), and/or all of the steps of any method or process so disclosed, may be combined in any combination, except combinations where at least some of such features and/or steps are mutually exclusive. The example

embodiments not restricted to the details of any foregoing embodiments. The example embodiments extend to any novel one, or any novel combination, of the features disclosed in this specification (including any accompanying claims, abstract and drawings), or to any novel one, or any novel combination, of the steps of any method or process so disclosed.

[00116] The reader's attention is directed to all papers and documents which are filed concurrently with or previous to this specification in connection with this application and which are open to public inspection with this specification, and the contents of all such papers and documents are incorporated herein by reference.

[00117] In the drawings and specification, there have been disclosed exemplary embodiments. However, many variations and modifications may be made to these embodiments. Accordingly, although specific terms are employed, they are used in a generic and descriptive sense only and not for purposes of limitation, the scope of the embodiments being defined by the following claims.