CN1635680A - Quick speed vari-data window phasor solving method - Google Patents

Abstract

This invention refers to a method for obtaining quickly the phasor of meta-data window, which contains the computer sampling the voltage and current in high voltage line through mutual inductor, filtering attenuation component, partially summing sampling value by Fourier series algorithm, calculating phasor, transmitting phasor to opposite side of transmission line. Said invention reduces CPU calculation with short data window and raises the action speed of current differential protection, adaptive to the differential protection of high voltage power transmission line.

Description

Quick Variable Data Window Phasor Extraction Method

Technical field

The invention belongs to the electroporation field, especially a kind of quick Variable Data Window Phasor Extraction Method is asked for electric current and voltage phasor in the high-voltage fence current differential protection.

Background technology

Along with high performance 32 bit DSP digital signal processors gradually adopting in protective relaying device, the real time data sampling disposal ability can be brought up to every power frequency period 96 point samplings.Yet the traffic rate that offers current differential protection at present is lower.If the transmission sampled value will be wasted lot of data information.And adopt the way that transmits phasor, and the abundant information composite phasor of high-speed sampling can be sent to offside, can additionally not increase the burden of communication again, also alleviated the CPU amount of calculation simultaneously.Yet traditional phasor acquiring method (as the half cycle fourier algorithm of widely using at present, full-wave fourier algorithm) data window is longer, influences the responsiveness of current differential protection.

Half cycle fourier algorithm (earlier through difference, the filtering DC component):

i(k)＝i(k)-i(k-n)

$I_{\mathrm{real}}=\frac{4}{N}\underset{k=0}{\overset{N/2-1}{\mathrm{\Σ}}}\left[i\left(k\right)\cos \left(\frac{2\mathrm{k\π}}{N}\right)\right]$

$I_{\mathrm{imag}}=-\frac{4}{N}\underset{k=0}{\overset{N/2-1}{\mathrm{\Σ}}}\left[i\left(k\right)\sin \left(\frac{2\mathrm{k\π}}{N}\right)\right]$

Full-wave fourier algorithm (earlier through difference, the filtering DC component):

i(k)＝i(k)-i(k-n)

$I_{\mathrm{real}}=\frac{2}{N}\underset{k=0}{\overset{N-1}{\mathrm{\Σ}}}\left[i\left(k\right)\cos \left(\frac{2\mathrm{k\π}}{N}\right)\right]$

$I_{\mathrm{imag}}=-\frac{2}{N}\underset{k=0}{\overset{N-1}{\mathrm{\Σ}}}\left[i\left(k\right)\sin \left(\frac{2\mathrm{k\π}}{N}\right)\right]$

I in the above-mentioned expression formula (k), i (k-n) is a sample sequence, I _Real, I _ImagBe respectively the real part and the imaginary part of phasor.

The basic thought of fourier algorithm originates from Fourier series, and principle is simple, and the computational accuracy height has obtained extensive use in Microcomputer Protection.But because the data window width of this algorithm is longer, full-wave fourier algorithm needs a power frequency period (20ms) data window, and the half cycle fourier algorithm needs half power frequency period (10ms) data window, has reduced the protection responsiveness.Adopt above-mentioned phasor transfer approach, when system breaks down, could calculate phasor when data window reaches 1/2nd power frequency periods after fault, add data transmission period, make the differential protection responsiveness impossible fast.

Summary of the invention

At above-mentioned the deficiencies in the prior art, the object of the present invention is to provide a kind of quick phasor acquiring method, under the prerequisite that guarantees reliability, improve the responsiveness of current differential protection effectively.

To achieve the above object of the invention, technical scheme of the present invention is: quick Variable Data Window Phasor Extraction Method may further comprise the steps:

A. computer is sampled to the voltage and current of high-tension line by instrument transformer, and the filter attenuation components;

B. use the Fourier series algorithm ask for a certain component of sampled value part and;

C. to the part of sampled value component with carrying out phasor calculation;

D. phasor is transmitted the offside of index electric line.

Above-mentioned said filter attenuation components is meant uses m the output that filtering algorithm is tried to achieve, and it equals (m current sampling data-m-6 current sampling data) * 2.563, that is:

out _m＝(i _m-i _m-6)·2.563

The said Fourier of above-mentioned steps b part and computational methods be:

${\mathrm{Real}}_p=\underset{k=p\·24-24+1}{\overset{p\·24}{\mathrm{\Σ}}}\sin \left(\frac{2\·k\·\mathrm{\π}}{9}\right)\·\mathrm{ou}{t}_k$

$\mathrm{Imaginar}{y}_p=\underset{k=p\·24-24+1}{\overset{p\·24}{\mathrm{\Σ}}}\cos \left(\frac{2\·k\·\mathrm{\π}}{96}\right)\·\mathrm{ou}{t}_k$

Wherein: p Fu Shi of Real p=part and real part, p Fu Shi of Imaginary p=partly and imaginary part.

The computational methods of the said phasor of above-mentioned steps c are:

$\mathrm{Phasor}\mathrm{Rea}{l}_p=a_p\·\underset{k=1}{\overset{p}{\mathrm{\Σ}}}\mathrm{Rea}{l}_k+c_p\·\underset{k=1}{\overset{p}{\mathrm{\Σ}}}\mathrm{Imaginar}{y}_k$

$\mathrm{Phasor}\mathrm{Imaginar}{y}_p=b_p\·\underset{k=1}{\overset{p}{\mathrm{\Σ}}}\mathrm{Rea}{l}_k+d_p\·\underset{k=1}{\overset{p}{\mathrm{\Σ}}}\mathrm{Imaginar}{y}_k$

Wherein, p=1,2,3,4

a ₁＝6.451 b ₁＝-4.279 c ₁＝-4.279 d ₁＝7.012

a ₂＝2 b ₂＝0 c ₂＝0 d ₂＝2

a ₃＝1.376 b ₃＝-0.295 c ₃＝-0.295 d ₃＝1.415

a ₄＝1 b ₄＝0 c ₄＝0 d ₄＝1

The real part of p phasor of Phasor Real p=, the imaginary part of p phasor of Phasor Imaginary p=.

Above-mentioned steps c further comprises: calculate phasor when data window after the fault reaches 1/4th power frequency periods first, when reaching 1/2nd power frequency periods, calculates by data window after the fault phasor for the second time, when data window after the fault reaches 3/4ths power frequency periods, calculate phasor for the third time, after data window after the fault reaches a power frequency period, adopt complete all Fourier filtering algorithms to ask for phasor; Said minimum data window is 1/4th cycles.

Advantage of the present invention and beneficial effect are: it can be sent to offside with the abundant information composite phasor of high-speed sampling, can additionally not increase the burden of communication, has alleviated the CPU amount of calculation, and data window is shorter, has improved the responsiveness of current differential protection.Data window is short; not only improved the protection responsiveness; be beneficial to the stable operation of electric power system; also make protection have the saturated ability of anti-TA; this is because TA has the one linear transform period about 5ms before saturated; the short data window can make full use of this period and whether differentiate troubles inside the sample space, thereby greatly improves the performance of current differential protection.

Description of drawings

Fig. 1 is a main flow chart of the present invention;

Fig. 2 is the current waveform figure of explanation Fu Shi part of the present invention and calculating.

Embodiment

The present invention is further described in detail below in conjunction with the drawings and specific embodiments.

Consult Fig. 1 main flow chart of the present invention, quick Variable Data Window Phasor Extraction Method, be the application sample value do the Fourier series algorithm part and, ask for phasor then, the minimum data window that calculates phasor according to quick Variable Data Window Phasor Extraction Method is 1/4th cycles.The performing step of the method is: at first the voltage and current of high-tension line is sampled by instrument transformer by computer, and filter attenuation components, the filter attenuation components is meant uses m the output that filtering algorithm is tried to achieve, it equals (m current sampling data-m-6 current sampling data) * 2.563, that is:

out _m＝(i _m-i _m-6)·2.563

Wherein: out _mM output of=filtering algorithm, i _m=the m current sampling data.

Use then the Fourier series algorithm ask for sampled value part and, p Fu Shi partly and real part be:

$\mathrm{Rea}{l}_p=\underset{k=p\·24-24+1}{\overset{p\·24}{\mathrm{\Σ}}}\sin \left(\frac{2\·k\·\mathrm{\π}}{96}\right)\·\mathrm{ou}{t}_k$

P Fu Shi part and imaginary part be:

$\mathrm{Imaginar}{y}_p=\underset{k=p\·24-24+1}{\overset{p\·24}{\mathrm{\Σ}}}\cos \left(\frac{2\·k\·\mathrm{\π}}{96}\right)\·\mathrm{ou}{t}_k$

Wherein, p=1,2,3,4.

When reaching 1/4th power frequency periods, calculates by data window after the fault phasor first, when reaching 1/2nd power frequency periods, calculates by data window after the fault phasor for the second time, when data window after the fault reaches 3/4ths power frequency periods, calculate phasor for the third time, after data window after the fault reaches a power frequency period, adopt complete all Fourier filtering algorithms to ask for phasor.Minimum data window in the embodiment of the invention is 1/4th cycles.

Calculate phasor afterwards again, the real part of p phasor is:

$\mathrm{Phasor}\mathrm{Rea}{l}_p=a_p\·\underset{k=1}{\overset{p}{\mathrm{\Σ}}}\mathrm{Rea}{l}_k+c_p\·\underset{k=1}{\overset{p}{\mathrm{\Σ}}}\mathrm{Imaginar}{y}_k$

The imaginary part of p phasor is:

$\mathrm{Phasor}\mathrm{Imaginar}{y}_p=b_p\·\underset{k=1}{\overset{p}{\mathrm{\Σ}}}\mathrm{Rea}{l}_k+d_p\·\underset{k=1}{\overset{p}{\mathrm{\Σ}}}\mathrm{Imaginar}{y}_k$

Wherein, p=1,2,3,4

a ₁＝6.451 b ₁＝-4.279 c ₁＝-4.279 d ₁＝7.012

a ₂＝2 b ₂＝0 c ₂＝0 d ₂＝2

a ₃＝1.376 b ₃＝-0.295 c ₃＝-0.295 d ₃＝1.415

a ₄＝1 b ₄＝0 c ₄＝0 d ₄＝1

The phasor that calculating is finished is sent to the offside of ultra-high-tension power transmission line at last.

It should be noted last that: above embodiment is the unrestricted technical scheme of the present invention in order to explanation only, although the present invention is had been described in detail with reference to the foregoing description, those of ordinary skill in the art is to be understood that: still can make amendment or be equal to replacement the present invention, and not breaking away from any modification or partial replacement of the spirit and scope of the present invention, it all should be encompassed in the middle of the claim scope of the present invention.

Claims (5)

1, quick Variable Data Window Phasor Extraction Method, it is characterized in that: this method may further comprise the steps:

A. computer is sampled to the voltage and current of high-tension line by current transformer, voltage transformer, and the filter attenuation components;

B. use the Fourier series algorithm ask for a certain component of sampled value part and;

C. to the part of sampled value component with carrying out phasor calculation;

D. phasor is sent to the offside of transmission line.

2, quick Variable Data Window Phasor Extraction Method according to claim 1 is characterized in that: the said Fourier of step b part and computational methods be:

${\mathrm{Real}}_p=\underset{k=p\·24-24+1}{\overset{p\·24}{\mathrm{\Σ}}}\sin \left(\frac{2\·k\·\mathrm{\π}}{96}\right)\·{\mathrm{out}}_k$

$\mathrm{Imaginar}{y}_p=\underset{k=p\·24-24+1}{\overset{p\·24}{\mathrm{\Σ}}}\cos \left(\frac{2\·k\·\mathrm{\π}}{96}\right)\·{\mathrm{out}}_k$

Wherein: p Fu Shi of Real p=part and real part, p Fu Shi of Imaginary p=partly and imaginary part.

3, quick Variable Data Window Phasor Extraction Method according to claim 1 is characterized in that: the said computational methods that get phasor of step c are:

$\mathrm{PhasorRe}{\mathrm{al}}_p=a_p\·\underset{k=1}{\overset{p}{\mathrm{\Σ}}}\mathrm{Re}{\mathrm{al}}_k+c_p\·\underset{k=1}{\overset{p}{\mathrm{\Σ}}}\mathrm{Im}{\mathrm{aginary}}_k$

$\mathrm{PhasorImagina}{\mathrm{ry}}_p=b_p\·\underset{k=1}{\overset{p}{\mathrm{\Σ}}}{\mathrm{Real}}_k+d_p\·\underset{k=1}{\overset{p}{\mathrm{\Σ}}}\mathrm{Imaginar}{y}_k$

p＝1，2，3，4

a ₁＝6.451 b ₁＝-4.279 c ₁＝-4.279 d ₁＝7.012

a ₂＝2 b ₂＝0 c ₂＝0 d ₂＝2

a ₃＝1.376 b ₃＝-0.295 c ₃＝-0.295 d ₃＝1.415

a ₄＝1 b ₄＝0 c ₄＝0 d ₄＝1

The real part of p phasor of Phasor Real p=,

The imaginary part of p phasor of Phasor Imaginary p=.

4, quick Variable Data Window Phasor Extraction Method according to claim 1, it is characterized in that: step c further comprises: calculate phasor when data window after the fault reaches 1/4th power frequency periods first, when reaching 1/2nd power frequency periods, calculates by data window after the fault phasor for the second time, when data window after the fault reaches 3/4ths power frequency periods, calculate phasor for the third time, after data window after the fault reaches a power frequency period, adopt complete all Fourier filtering algorithms to ask for phasor.

5, quick Variable Data Window Phasor Extraction Method according to claim 4 is characterized in that: said minimum data window is 1/4th cycles.

Images