Method for measuring frequency of phasor of power system

Abstract

The present invention relates to a method for measuring the frequency of phasor of a power system, belonging to the technical field of power systems and automation thereof. A grid-connected photovoltaic generation system is provided with a control module, a data buffer module, an inverter, a sampling module and an analog-to-digital conversion module; a voltage transformer is connected to the grid side of the grid-connected photovoltaic generation system; the voltage transformer is electrically connected to the analog-to-digital conversion module; the analog-to-digital conversion module converts an analog signal sent by the voltage transformer into a digital signal; the sampling module samples the digital signal and transmits the digital signal into the data buffer module; the data buffer module is provided with a plurality of data windows; and the control module receives signals from the data windows.

Claims

1. A method for measuring a frequency of phasor of a power system, comprising the following steps: (1) providing a grid-connected photovoltaic generation system with a control module (10), a data buffer module (20), an inverter (30), a sampling module (50) and an analog-to-digital conversion module (60), a voltage transformer (70) being connected to a grid side of the grid-connected photovoltaic generation system, the voltage transformer (70) being electrically connected to the analog-to-digital conversion module (60), the analog-to-digital conversion module (60) converting an analog signal sent by the voltage transformer (70) into a digital signal, the sampling module (50) sampling the digital signal and transmitting the digital signal into the data buffer module (20), the data buffer module (20) being provided with a plurality of data windows (21), the control module (10) receiving signals from the data windows (21); (2) the control module (10) performing Fourier transform calculation to the signals of the data windows (21) in the data butter module (20): $\begin{array}{c}X\left(M,l\right)=\underset{n=-\left(N-1\right)/2}{\overset{\left(N-1\right)/2}{.Math.}}{x}_a\left(n+l\right)h\left(n\right){}^{-j\left(n+l\right){}_0}\\ =\underset{n=-\left(N-1\right)/2}{\overset{\left(N-1\right)/2}{.Math.}}\left[\underset{k=0}{\overset{K}{.Math.}}\frac{1}{k!}{a^{\left(k\right)}\left({}_{n+l}\right)}^k\right].Math.\\ {}^{j\left(n+l\right){}_0}.Math.h\left(n\right).Math.{}^{-j\left(n+l\right){}_0}+\\ \underset{n=-\left(N-1\right)/2}{\overset{\left(N-1\right)/2}{.Math.}}{\left[\underset{k=0}{\overset{K}{.Math.}}\frac{1}{k!}{a^{\left(k\right)}\left({}_{n+l}\right)}^k\right]}^{*}.Math.\\ {}^{-j\left(n+l\right){}_0}.Math.h\left(n\right).Math.{}^{-j\left(n+l\right){}_0}\\ \underset{k=0}{\overset{K}{.Math.}}\left[a^{\left(k\right)}.Math.\frac{1}{k!}\underset{n=-\left(N-1\right)/2}{\overset{\left(N-1\right)/2}{.Math.}}{\left({}_{n+l}\right)}^k.Math.h\left(n\right)\right]+\\ \underset{k=0}{\overset{K}{.Math.}}\left\{{\left[a^{\left(k\right)}\right]}^{*}.Math.\frac{1}{k!}\underset{n=-\left(N-1\right)/2}{\overset{\left(N-1\right)/2}{.Math.}}{\left({}_{n+l}\right)}^k.Math.h\left(n\right).Math.{}^{-j2\left(n+l\right){}_0}\right\}\\ =\underset{k=0}{\overset{K}{.Math.}}\left[a^{\left(k\right)}.Math.H\left(k,l,0\right)\right]+\underset{k=0}{\overset{K}{.Math.}}\left\{{\left[a^{\left(k\right)}\right]}^{*}.Math.H\left(k,l,-2{}_0\right)\right\}\\ =D_{K,l}.Math.A+C_{K,l}.Math.A^{*},\end{array}$ wherein, X(M,l) is a transformed value of a measured result of phasor, M is a phasor reference point in the data windows (21), l is a point distance from each sampling point to the reference point, N is a time span of the data windows (21), h(n) is a window function of the data windows (21), .sub.0 is a time factor of the window function of the data windows (21), x.sub.a (n+l) is a value of a n+1.sup.th sampling point in the data windows (21), a.sup.(k) is a k.sup.th derivative when the signals of the data windows (21) are at t=0, a factor of the k.sup.th derivative a.sup.(k) is $H\left(k,l,\right)=\frac{1}{k!}\underset{n=-\left(N-1\right)/2}{\overset{\left(N-1\right)/2}{.Math.}}{\left({}_{n+l}\right)}^k.Math.h\left(n\right).Math.{}^{j\left(n+l\right)}\},$ and, D.sub.K,l and C.sub.K,l are a column of the factor of the derivative a.sup.(k) and a column matrix a factor of a conjugate transpose derivative [a.sup.(k)]*, respectively; (3) the phasor reference point M being unified into one point, and the transformed value X(M,l) of the measured result of phasor being phase-shifted and calculated to obtain a phasor estimated value {circumflex over (X)}(M,l):
{circumflex over (X)}(M,l)=e.sup.jl.sup.0.Math.{circumflex over (X)}(M+l,0); (4) the control module (10) combining the transformed value X(M,l) of the measured result of phasor with the phasor estimated value {circumflex over (X)}(M,l) and performs recombination:
E.Math.X=D.Math.A+C.Math.A*, wherein, e.sup.jl.sup.0 is a vector corresponding to the sampling point of the data windows (21), a diagonal matrix E=diag(e.sup.j.sup.0.sup.l.sup.ML, . . . , e.sup.j.sup.0.sup.l.sup.M+L) is constructed relative to e.sup.jl.sup.0; a diagonal matrix X=[{circumflex over (X)}(M,L), . . . , {circumflex over (X)}(M,L)].sup.T is constructed relative to the matrix of {circumflex over (X)}(M,l); D=[(D.sub.K,l.sub.ML).sup.T, . . . , (D.sub.K,l.sub.m+L).sup.T].sup.T; C=[(C.sub.K,l.sub.ML).sup.T, . . . , (C.sub.K,l.sub.M+L).sup.T].sup.T; A=)[a.sup.(0), . . . , a.sup.(K)].sup.T; meanwhile, the corresponding phasor estimated value {circumflex over (X)}(M,l) at the phasor reference point M in the data windows (21) is represented by Taylor's series, and a modified phasor value X(M) corresponding to the phasor reference point M in the data windows (21) is obtained by a least square method; and (5) the control module (10) performing data processing to the modified phasor value X(M) to obtain a modified frequency value F(M).

2. The method for measuring the frequency of phasor of a power system according to claim 1, wherein the inverter (30) combines the modified phasor value X(M) with the modified frequency value F(M) to output the phasor and frequency of a drive signal.

3. The method for measuring the frequency of phasor of a power system according to claim 2, wherein the control module (10) is electrically connected to a communication unit (40), the inverter (30) outputting and transmitting the phasor and frequency of the drive signal to a network via the communication unit (40), the communication unit (40) performing communication through IEEE standard C37.118.

4. The method for measuring the frequency of phasor of a power system according to claim 1, wherein the calculation formula between the modified phasor value X(M) and the modified frequency value is F(M) is: $F\left(M\right)=\frac{\mathrm{angle}\left[X\left(M_1\right).Math.{X\left(M_2\right)}^{*}\right]}{2\left(T_2-T_1\right)},$ wherein, M.sub.1 and M.sub.2 are window center sampling points of two different data windows 21, T.sub.1 and T.sub.2 are times corresponding to M.sub.1 and M.sub.2, respectively, and angle[X(M.sub.i).Math.X(M.sub.2)*] is an included angle between the phasor values X(M.sub.1) and X(M.sub.2).

5. The method for measuring the frequency of phasor of a power system according to claim 1, wherein the sampling module (50) is electrically connected to a GPS receiver (80), the GPS receiver (80) marking the digital signals collected by the sampling module by an absolute time mark (50), the data buffer module (20) receiving the signals of the data windows (21) marked by the absolute time mark, the modified phasor value X(M) and the modified frequency value F(M) being synchronous with the absolute time mark of the signals of the data windows (21).

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a structural diagram of a grid-connected photovoltaic generation system according to the present invention;

(2) FIG. 2 is a flowchart of a method for measuring the frequency of phasor of a power system according to the present invention; and

(3) FIG. 3 is a schematic diagram of a data window of a conventional sampling module according to the present invention.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

(4) In order to further understand the features, technical means, specific objects and functions of the present invention, the present invention will be further described as below in details with reference to the accompanying drawings by specific embodiments.

(5) As shown in FIG. 1 and FIG. 2, the present invention provides a method for measuring the frequency of phasor of a power system, implemented through the following steps:

(6) (1) A grid-connected photovoltaic generation system is provided with a control module 10, a data buffer module 20, an inverter 30, a communication unit 40, a sampling module 50 and an analog-to-digital conversion module 60. A voltage transformer 70 is connected to the grid side of the grid-connected photovoltaic generation system in order to reduce the working voltage of the grid side. The voltage transformer 70 is electrically connected to the analog-to-digital conversion module 60. The analog-to-digital conversion module 60 converts an analog signal sent by the voltage transformer 70 into a digital signal. The sampling module 50 is electrically connected to a GPS (Global Positioning System) receiver 80. The sampling module 50 samples the digital signal and transmits the digital signal into the data buffer module 20. The GPS receiver 80 marks the digital signals collected by the sampling module 50 by an absolute time mark. The data buffer module 20 is provided with a plurality of data windows (21). The control module (10) receives the signals marked by the absolute time mark from the data windows (21), and carries out calculation to modify the output phasor and frequency of the inventor 30.

(7) (2) The control module 10 performs Fourier transform to signals of the data windows 21 in the data buffer module 20:

(8) $\begin{array}{c}X\left(M,l\right)=\underset{n=-\left(N-1\right)/2}{\overset{\left(N-1\right)/2}{.Math.}}{x}_a\left(n+l\right)h\left(n\right){}^{-j\left(n+l\right){}_0}\\ =\underset{n=-\left(N-1\right)/2}{\overset{\left(N-1\right)/2}{.Math.}}\left[\underset{k=0}{\overset{K}{.Math.}}\frac{1}{k!}{a^{\left(k\right)}\left({}_{n+l}\right)}^k\right].Math.\\ {}^{j\left(n+l\right){}_0}.Math.h\left(n\right).Math.{}^{-j\left(n+l\right){}_0}+\\ \underset{n=-\left(N-1\right)/2}{\overset{\left(N-1\right)/2}{.Math.}}{\left[\underset{k=0}{\overset{K}{.Math.}}\frac{1}{k!}{a^{\left(k\right)}\left({}_{n+l}\right)}^k\right]}^{*}.Math.\\ {}^{-j\left(n+l\right){}_0}.Math.h\left(n\right).Math.{}^{-j\left(n+l\right){}_0}\\ \underset{k=0}{\overset{K}{.Math.}}\left[a^{\left(k\right)}.Math.\frac{1}{k!}\underset{n=-\left(N-1\right)/2}{\overset{\left(N-1\right)/2}{.Math.}}{\left({}_{n+l}\right)}^k.Math.h\left(n\right)\right]+\\ \underset{k=0}{\overset{K}{.Math.}}\left\{{\left[a^{\left(k\right)}\right]}^{*}.Math.\frac{1}{k!}\underset{n=-\left(N-1\right)/2}{\overset{\left(N-1\right)/2}{.Math.}}{\left({}_{n+l}\right)}^k.Math.h\left(n\right).Math.{}^{-j2\left(n+l\right){}_0}\right\}\\ =\underset{k=0}{\overset{K}{.Math.}}\left[a^{\left(k\right)}.Math.H\left(k,l,0\right)\right]+\underset{k=0}{\overset{K}{.Math.}}\left\{{\left[a^{\left(k\right)}\right]}^{*}.Math.H\left(k,l,-2{}_0\right)\right\}\\ =D_{K,l}.Math.A+C_{K,l}.Math.A^{*},\end{array}$

(9) where, X(M,l) is a transformed value of the measured result of phasor, M refers to phasor reference points in the data windows, l is a point distance from each sampling point to the reference point in the data windows, N is the time span of the data windows, h(n) is a window function of the data windows, .sub.0, is the time coefficient of the window function of the data windows, x.sub.a(n+l) is the value of the (n+1)th sampling point in the data windows, a.sup.(k) is the kth order derivative of the signals of the data windows at t=0, the coefficient of the kth order derivative a.sup.(k) is

(10) $H\left(k,l,\right)=\frac{1}{k!}\underset{n=-\left(N-1\right)/2}{\overset{\left(N-1\right)/2}{.Math.}}{\left({}_{n+l}\right)}^k.Math.h\left(n\right).Math.{}^{j\left(n+l\right)}\},$
and, D.sub.K,l and C.sub.K,l are a coefficient column matrix of the derivative a.sup.(k) and a coefficient column matrix of a conjugate transposed derivative [a.sup.(k)]*, respectively.

(11) (3) The phasor reference points M are unified into one point, and the transformed value X(M,l) of the measured result of phasor is phase-shifted and calculated to obtain an estimated phasor value (M,l):
(M,l)=e.sup.jl.sub.0.Math.(M+l,0)

(12) (4) The control module 10 combines the transformed value X(M,l) of the measured result of phasor with the estimated phasor value (M,l) and performs recombination:
E.Math.X=D.Math.A+C.Math.A*

(13) where, e.sup.jl.sub.0 is a vector corresponding to the sampling point of the data windows; a diagonal matrix E=diag(e.sup.j.sup.0.sup.l.sup.ML, . . . , e.sup.j.sup.0.sup.l.sub.M+L) is constructed relative to e.sup.jl.sub.0; a diagonal matrix X=[(M,L), . . . , (M,L)].sup.T is obtained from the matrix of (M,l); D=[(D.sub.K,l.sub.ML).sup.T, . . . , (D.sub.K,l.sub.m+L).sup.T].sup.T; C=[(C.sub.K,l.sub.ML).sup.T, . . . , (C.sub.K,l.sub.M+L).sup.T].sup.T; A=[a.sup.(0), . . . , a.sup.(K)].sup.T.

(14) Meanwhile, the corresponding estimated phasor value (M,l) at the phasor reference points M in the data windows is represented by Taylor series, and a modified phasor value X(M) corresponding to the phasor reference point M in the data windows is obtained by the least square method; and

(15) (5) The control module 10 performs data processing to the modified phasor value X(M) to obtain a modified frequency value F(M):

(16) $F\left(M\right)=\frac{\mathrm{angle}\left[X\left(M_1\right).Math.{X\left(M_2\right)}^{*}\right]}{2\left(T_2-T_1\right)},$

(17) where, M.sub.1 and M.sub.2 are central sampling points of two different data windows 21, T.sub.1 and T.sub.2 are the time corresponding to M.sub.1 and M.sub.2, respectively, and angle[X(M.sub.1).Math.X(M.sub.2)*] is an included angle between the phasor values X(M.sub.1) and X(M.sub.2).

(18) (6) The modified phasor value X(M) and the modified frequency value F(M) are synchronous with the absolute time mark of the signals of the data windows. The inverter 30 combines the modified phasor value X(M) with the modified frequency value F(M) to output the phasor and frequency of a driving signal, and transmitting the phasor and frequency of the driving signal to a network via the communication unit 40 after marking the absolute time mark. The communication unit 40 performs communication over IEEE standard C37.118.

(19) In conclusion, by using Taylor series to represent the corresponding estimated phasor value X(M,l) at the phasor reference points M in the data) windows for re-modification, and obtaining the precise modified phasor value X(M) corresponding to the phasor reference points M in the data windows by the least square method, the method for measuring the frequency of phasor of a power system provided by the present invention makes the automatic safety device of the generation system free from the interference of the reference level and is flexible, avoids the phenomenon of frequency spectrum leakage caused by asynchronous sampling, and effectively eliminates or weakens an oscillation effect caused by dynamic characteristics.

(20) The foregoing embodiment merely shows one implementation way of the present invention. Although this embodiment has been described specifically and in details, it is not intended to limit the scope of the present invention. It is to be pointed out that, a person of ordinary skill in the art may make various changes and improvements without departing from the idea of the present invention, and these changes and improvements shall fall into the protection scope of the present invention. Therefore, the protection scope shall be subject to the appended claims.