CN106645919B - The full spectral power of power grid based on three-phase instantaneous power vibrates phasor method for synchronously measuring - Google Patents
The full spectral power of power grid based on three-phase instantaneous power vibrates phasor method for synchronously measuring Download PDFInfo
- Publication number
- CN106645919B CN106645919B CN201611249360.0A CN201611249360A CN106645919B CN 106645919 B CN106645919 B CN 106645919B CN 201611249360 A CN201611249360 A CN 201611249360A CN 106645919 B CN106645919 B CN 106645919B
- Authority
- CN
- China
- Prior art keywords
- frequency range
- power
- frequency
- phase
- oscillation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R21/00—Arrangements for measuring electric power or power factor
- G01R21/001—Measuring real or reactive component; Measuring apparent energy
- G01R21/002—Measuring real component
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
A kind of device carrying out electric system entire spectrum oscillation of power phasor synchro measure and analysis based on the total instantaneous power of three-phase.This method is measured to obtain the total instantaneous active power of three-phase by three-phase instantaneous voltage and transient current recording in grid branch measurement point, and is that each measuring point obtains high-precision markers by time synchronism apparatus.Then frequency-division section processing technique is used, seeks the parameter for the oscillation mode that each frequency range includes using specific filtering parameter, down-sampled rate parameter, the calculating parameters such as data window and material calculation for each frequency range；And the phase angle with the initial phase at whole scale moment (10 seconds, the whole second, 0.1 second etc.) as oscillation phasor, to realize the expression that can carry out the oscillation phasor of wide area analysis.In order to further increase the precision of phasor calculation, the selection at whole scale moment avoids the boundary point of computation window；Frequency, amplitude, initial phase, damped coefficient and the markers of the oscillation of power for meeting threshold value detected in each frequency range are finally exported in real time.
Description
Technical field
The application belongs to power system measuring analysis field, in particular to various oscillation of powers carry out in a kind of pair of power grid
Line measuring technique.
Background technique
The oscillation of power of various frequencies can all bring different harm to operation of power networks in electric system, run in control centre
Personnel wish can know that the oscillation of power which frequency mode is occurring in power grid, and how is the coverage area of oscillation, vibrate function
Rate is in the flowing relation of equipment room, and each equipment is to information such as the percentage contributions of hunting power.
Recently as the new energy power plant such as high voltage direct current, soft straight, wind-powered electricity generation, photovoltaic be widely used inverter, frequency converter, can
The FACTS technologies such as string benefit are controlled, are even more a large amount of non-power frequency component occur in power grid.These non-power frequency ingredients include subsynchronous vibration
Swing, electrical resonance, control system resonance, harmonic wave, fractional harmoni, m-Acetyl chlorophosphonazo and these waveforms and power frequency and they between
The waveforms complexity such as phase inter-modulation frequency conversion electricity.The presence of these ingredients can cause following harm, such as sub-synchronous oscillation meeting
Endanger the service life of the big axis of generating set, m-Acetyl chlorophosphonazo can cause voltage fluctuation and flicker, the overload (m-Acetyl chlorophosphonazo of passive filter
Resonance amplification caused by), interference power line on control, protection and communication signal, cause Mechatronic Systems low-frequency oscillation, influence with electricity
Press through control equipment and the certain household electrical appliance normal work etc. that zero point is synchronization signal.Therefore the non-power frequency ingredient of power grid
Must control in certain level hereinafter, this just need to generation, development, sphere of action and the oscillation source of these non-power frequency ingredients into
Row is paid close attention to.
However, the existing grid monitoring system based on PMU and RTU in power-management centre be for power frequency and its near
The frequency of range is measured, and is the measuring tool based on voltage phasor, electric current phasor.Such as although PMU to dispatch
Center personnel is provided with the ability that Fast synchronization observation is carried out to electrical power system wide-area dynamic event, but PMU is designed from principle
On be only applicable to that the electromechanical dynamic process near power frequency is measured and monitored, be not particularly suited for sub-synchronous oscillation, harmonic wave, point
The measurement and monitoring of the electromagnetic transients such as harmonic wave, m-Acetyl chlorophosphonazo for several times.This is because for example for sub-synchronous oscillation, it is typical
Mechanical oscillation frequency or power sub-synchronous oscillation frequency are 10~45Hz, and corresponding voltage sub-synchronous oscillation frequency is 5~40Hz,
After it is to power-frequency voltage amplitude modulation, the sum of two sinusoidal signals can be converted into using the product to sum formula of trigonometric function, point
It Wei not the subsynchronous sinusoidal signal of 10~45Hz and the supersynchronous sinusoidal signal of 55~95Hz.And international and national standard provides
The effective range of PMU is 45~55Hz, even if the frequency more than the range is not filtered, amplitude also will significantly disappear
Subtract, is also lost or has distorted the letter of the higher-orders of oscillation such as sub-synchronous oscillation according to this calculated power of voltage and current information
Breath.Therefore, PMU is not suitable for the monitoring of sub-synchronous oscillation.Similarly PMU is not suitable for harmonic wave, fractional harmoni and m-Acetyl chlorophosphonazo yet
Measurement and monitoring.
Although the recorder data of failure plant stand can also be called in real time in failure, although this original transient state recorder data
Frequency synthesis ingredient comprising power frequency and other frequency contents and power frequency, but operations staff cannot be bent from this original recording
Line intuitively observes the change procedure of non-power frequency ingredient.The system having in recent years, including the monitoring of power quality device, sub-synchronous oscillation
Device, although control can be sent by non-power frequency ingredient, such as in the amplitude of sub-synchronous oscillation ingredient, frequency, decay factor parameter
Central master station, but oscillation phase and synchronous time mark due to lacking non-power frequency ingredient cause operation, analysis personnel cannot be intuitive
Analyze each energy exchange relationship of the grid equipment in oscillatory process for participating in oscillation, it is difficult to realize determining for disturbing source and oscillation source
Position.In addition, the wave-record chart of typically at least thousand Hz or more, due to the limit value of existing communication resource and storage resource, it is impossible to real
When a wide range of long-time on send scheduling station, therefore operations staff cannot dynamically monitor in real time for a long time these non-power frequencies at
Generation, the development process divided, even if using the method to recorder data down-sampled rate sampling Shang Song control centre, but can draw again
Enter aliasing effect, distorts the essential information of oscillation.
It should further be noted that the oscillation that operation of power networks is concerned about, such as sub-synchronous oscillation is function relevant to generator's motor shaft
Rate oscillation, the existing sub-synchronous oscillation detection device based on transient state recording are subsynchronous to monitor by analysis ac flow
(wave recording device used in electric system at present does not calculate instantaneous power in real time) of oscillation, in most cases electric current shakes
Swinging with the sum of the frequency of oscillation of power is power frequency 50Hz, that is, there is complementary relationship, but same in the multi-mode of a generator time
In step oscillation, complementary relationship is simultaneously not present in the frequency of some frequency of power oscillation and current oscillation, this is existing based on electric current record
One limitation of wave analysis oscillation.
The measurement of the generators direct current such as excitation voltage, exciting current, revolving speed is observed by PMU in addition, also having at present
And analysis, the monitoring to sub-synchronous oscillation is realized, but other frequencies is dry in the direct current signals such as excitation voltage, exciting current
It is more to disturb ingredient, i.e., noise is relatively low, therefore monitoring accuracy is poor.
Finally of particular note is that can not directly copy existing PMU directly to set common non-oscillation frequency
Count out phasor measuring set.Because different according to power-frequency voltage, current curve calculating power-frequency voltage, electric current phasor from routine,
The usually not pre-determined and invariable standard frequency of oscillation of power (such as calculate used in voltage, electric current phasor
Power frequency 50Hz standard frequency), this is just to calculate instantaneous phase angle according to discrete Fourier transform DFT to bring difficulty, because of discrete Fu
In leaf transformation theoretical basis be have a standard constant rate of speed rotation vector.Fortunately sometime power flow is flowed to
Analysis, take decisive action is the relative phase angle between two oscillating vectors at this moment.Therefore it can pass through selection time-varying
Some common reference phasor, the comparison of angle between the vector of the different location Lai Shixian two.As long as each device all follows same rule
Then look for the common reference phasor of time-varying.
Summary of the invention
For the above problem present in oscillation of power monitoring in existing power grid, the invention proposes one kind to be based on three-phase wink
When power power grid entire spectrum oscillation of power phasor method for synchronously measuring.Emphasis solves following problems:
(1) it realizes and the direct convenience of various known and unknown power oscillation of power system is observed.It solves existing subsynchronous
Device can not achieve to there is no the power sub-synchronous oscillations of complementary relationship to realize the problem of observing with electric current sub-synchronous oscillation.
(2) the solution ununified standard reference frequency of oscillation frequency amount or common reference phasor by no means, cannot be based on conventional
The problem of discrete fourier algorithm (DFT) in PMU device carries out phasor calculation.
(3) solve while monitoring the oscillation phasor of different frequency range, the especially higher-order of oscillation, caused computationally intensive, communication
Measure big problem.
(4) accuracy of phasor calculation is improved as far as possible.
The present invention specifically uses following technical scheme:
A kind of power grid entire spectrum oscillation of power phasor method for synchronously measuring based on three-phase instantaneous power；It is characterized by:
Measuring to obtain power grid three-phase always by three-phase instantaneous voltage and transient current recording in grid branch measurement point instantaneously has
Function power handles the parameter for seeking the oscillation mode that each frequency range includes using frequency-division section, and with the whole of non-computational window border point
Phase angle of the initial phase at scale moment as oscillation phasor, will finally meet the phasor of the oscillation mode of oscillation amplitude threshold requirement
Expression of results output and Shang Song main website.
The power grid entire spectrum oscillation of power phasor method for synchronously measuring based on three-phase instantaneous power, including following step
It is rapid:
Step 1: measuring to obtain power grid by the three-phase instantaneous voltage of grid branch measurement point and the recording of transient current corresponding
The total instantaneous active power curve of three-phase of tested branch, and whole moment second markers is marked using clock synchronization apparatus；
pr=uaia+ubib+ucic
Wherein, ua、ub、ucThe three-phase voltage of the tested branch of certain respectively directly measured from voltage transformer is instantaneous
Value；Ia, ib, ic are respectively the three-phase current instantaneous value of certain the tested branch directly measured from current transformer, prIt is three
Mutually total instantaneous active power；
Step 2: by the total instantaneous active power data p of above-mentioned three-phaserIt is distributed to the interior of 5 different frequency range processing modules in real time
It deposits in buffer area, uses fx1Indicate the initial frequency of frequency range X, the period of the corresponding initial frequency of frequency range X is Tx1=1/fx1；fx2Table
Show the termination frequency of frequency range x, the corresponding period for terminating frequency is Tx2=1/fx2；Each frequency range buffer area can at least store Tx1×4
The data of second；Wherein, x value 1-5 respectively indicates 1-5 different frequency range；The frequency range of frequency range 1 to frequency range 5 is to be sequentially increased
's；
Step 3: using recursion bandpass filter in real time to the three-phase for entering core buffer in each frequency range processing module
Total instantaneous active power data are filtered；
Step 4: the total instantaneous power data of filtered three-phase being handled using down-sampled rate in each frequency range processing module；
Step 5: with the week of initial frequency corresponding to this frequency range in the processing module of remaining each frequency range in addition to frequency range 5
2 times of values of phase are 2Tx1To calculate data window duration, with the cycle T of termination frequency corresponding to this frequency rangex2For calculating cycle
The moving step length for calculating data window, grid branch three-phase total wink of data window is calculated this section using frequency spectrum analysis method
When power data carry out spectrum analysis；For the processing module of frequency range 5 then with 0.05 second to calculate data window, 0.01 second all to calculate
Phase carries out spectrum analysis to the total instantaneous power data of grid branch three-phase in window at this section using frequency spectrum analysis method；
In spectrum analysis, hunting power amplitude is taken to be greater than given threshold, and can preceding 5 frequency of the accounting greater than 5%
Mode is as calculated result, that is, result of spectrum analysis；
Step 6: each frequency range chooses t respectivelyx2-Tx1The nearest following moment before/2 moment is the calculating for needing to export
As a result t at the time ofx-out, wherein tx2The finish time of data window is calculated for current frequency range X, and by the calculating of the calculating data window
As a result by when window start time tx1It is converted to tx-outMoment；
For frequency range 1, t is chosen12-T11At the time of whole 10 second moment is the calculated result for needing to export before/2 moment
t1-out；
For frequency range 2, t is chosen22-T21T at the time of whole moment second is the calculated result for needing to export before/2 moment2-out；
For frequency range 3, t is chosen32-T31At the time of whole 0.1 second moment is the calculated result for needing to export before/2 moment
t3-out；
For frequency range 4, t is chosen42-T41At the time of whole 0.01 second moment is the calculated result for needing to export before/2 moment
t4-out；
For frequency range 5, t is chosen52-T51At the time of whole 0.01 second moment is the calculated result for needing to export before/2 moment
t5-out；
Step 7: storing on the spot that each frequency range calculates in real time as a result, retaining at least 2 weeks；
Step 8: by the real-time Shang Song main website of the calculated result of each frequency range or according to main website to the interval for sending data in each frequency range
Sent on it is required that, frequency, amplitude, initial phase, damped coefficient including the oscillation of power for meeting threshold value detected in each frequency range with
And markers；
Step 9: each frequency range processing module waiting meets this frequency range calculating cycle Tx2After duration, return step 5.
The present invention still further comprises following preferred embodiment:
In step 1, markers precision is higher than 1 microsecond, original instantaneous recording curve sample rate at least 4800Hz or more.
In step 2,5 different frequency ranges are provided that
Frequency range 1 is load oscillation or power plant heat power engineering system oscillation frequency bands: f11-f12Range be (0.01Hz,
0.1Hz]；
Frequency range 2 is electric system section low-frequency oscillation frequency range: f21-f22Range be (0.1Hz, 1Hz]；
Frequency range 3 is oscillation or generator speed-adjusting and control system oscillation frequency bands between generator: f31-f32Range be (1Hz,
10Hz]；
Frequency range 4 is subsynchronous or supersynchronous oscillation frequency bands: f41-f42Range be (10Hz, 100Hz]；
Frequency range 5 is each secondary industrial frequency resonance frequency range: f51-f52Range be (100Hz, 1000Hz]；
In step 3, -3dB the free transmission range of recursion bandpass filter corresponding to each frequency range processing module is [fx1,
fx2]。
In step 4, in the processing module of frequency range 1, every 1 is retained to the total instantaneous power data of filtered power grid three-phase
Second whole moment second transient data；
In the processing module of frequency range 2, to the total instantaneous power data of filtered power grid three-phase, every N1 point retains newest
1 point, wherein N1 is 1/1000 smallest positive integral less than or equal to former sample rate；
In the processing module of frequency range 3, instantaneous power data total for filtered power grid three-phase, down-sampled sample rate
Become the 1/100 of former sample rate；
In the processing module of frequency range 4, instantaneous power data total for filtered power grid three-phase, down-sampled sample rate
Become the 1/10 of former sample rate；
In the processing module of frequency range 5, down-sampled sample rate is taken as the smallest crude sampling rate more than or equal to 4800
Approximate number.
In steps of 5, the frequency spectrum analysis method is Prony algorithm or EMD algorithm.
The given threshold is taken as in 30 second period before current calculation window starts total instantaneous active power curve and goes
Except the 5% of the average value of maximum point and smallest point absolute value after dc component.
In steps of 5, calculated result information includes: the frequency f for the hunting power that each frequency range includesx, amplitude Ax1, initial phase
φx1, decay factor σx, can accounting ex。
In step 6, each frequency mode for each frequency range being calculated and being selected by step 5, by tx1The power at moment shakes
The amplitude A swungx1With initial phase φx1, it is scaled tx-outThe amplitude A at momentx-outWith initial phase φx-out, initial phase be relative to
2 π f of moment sinxThe phase of t curve, specific formula is as follows:
φx-out=φx1+(tx-out-tx1)(2πfx)
Remaining information, i.e. the frequency f of hunting powerx, decay factor σx, can accounting exWhen still taking calculating data window starting point
Carve tx1When value.
According to the present invention the device of the method exploitation may be implemented to known to place plant stand various types, various frequencies or
The real-time monitoring and analysis of unknown oscillation of power form the function based on whole time scale initial phase for having synchronous time mark
Rate vibrates phasor, and in real time to send scheduling station in the acceptable traffic of existing power dispatch data network, is adjusting to realize
Spend wide area comparative analysis of the main website to the same oscillationg component of the Wide Area Power difference plant stand.
Method of the invention has technical effect beneficial below compared with the existing technology:
(1) it based on oscillatory occurences such as the total analysis of instantaneous power sub-synchronous oscillation of three-phase, resonance, realizes true to electric system
The direct observation for the oscillation of power just paid close attention to.It avoids based on the resulting frequency of current oscillation and the generator shafting be actually concerned about
The inconsistent problem of frequency of oscillation, and it is same to also achieve the power time that frequency complementary relationship is not present with electric current sub-synchronous oscillation
Walk the observation of oscillation.
(2) expression using 10 seconds, the initial phase at the whole scale moment such as whole second, 0.1 second as phasor phase angle, solution
The ununified standard reference frequency of non-oscillation frequency amount or common reference phasor, the measuring device of different location can not carry out
The problem of wide area Phasor Measurements compare.The common reference phasor of this time-varying is taken as the sinusoidal letter that whole scale moment initial phase angle is 0
Number.
(3) it uses and is calculated for different frequency range using the method reduction of different sample rates, different calculation windows, different step-lengths
Amount and the traffic.Following 5 frequency ranges are arranged in the characteristics of for power oscillation of power system, cover the complete of electric system concern
The frequency range of portion's oscillation of power: the load oscillation or the oscillation of power plant heat power engineering system of 0.01-0.1Hz；0.1Hz -1Hz electricity
The low-frequency oscillation of Force system section；Oscillation, generator speed regulation dispatch control system oscillation between 1Hz -10Hz generator；10Hz—100Hz
Subsynchronous and supersynchronous oscillation；The Harmonious Waves in Power Systems of 100Hz -1000Hz vibrates.
(4) by avoiding the calculated result using calculation window edge, the accuracy of phasor calculation is improved.
Detailed description of the invention
The full spectral power of power grid of the Fig. 1 based on three-phase instantaneous power vibrates phasor method for synchronously measuring flow diagram；
The parameter definition figure of Fig. 2 bandpass filter；
The full spectral power of power grid of the Fig. 3 based on three-phase instantaneous power vibrates phasor synchronous measuring apparatus physical function module map.
Specific embodiment
Technical solution of the present invention is described in further detail with specific embodiment with reference to the accompanying drawings of the specification.
Electric system entire spectrum oscillation of power phasor is carried out based on the total instantaneous power of three-phase the invention proposes one kind to synchronize
The method of measurement.This method is measured to obtain three-phase total wink by three-phase instantaneous voltage and transient current recording in grid branch measurement point
When active power, and obtain high-precision markers, then use frequency-division section processing technique, for each frequency range using specific filtering ginseng
Several, down-sampled rate parameter, calculating data window and material calculation parameter seek the parameter for the oscillation mode that each frequency range includes；With whole quarter
Phase angle of the initial phase at moment (10 seconds, the whole second, 0.1 second etc.) as oscillation phasor is spent, so that wide area analysis can be carried out by realizing
Oscillation phasor expression；In order to further increase the precision of phasor calculation, the selection at whole scale moment avoids computation window
Boundary point.Device based on the method for the invention exploitation is deployed in each power plant or substation, acquires generator in real time, becomes
Depressor or the instantaneous recorder data of the voltage and current of important line, and high-precision markers is obtained by global position system, it supervises accordingly
Survey power grid in various frequencies oscillation of power, and export the oscillation of power for meeting threshold value detected in each frequency range frequency,
Amplitude, initial phase, attenuation coefficient and markers, finally by dispatch data net by these information it is real-time on be sent to dispatching of power netwoks control
Central master station processed.The Comprehensive Analysis Software of scheduling station is based on these synchro measure information and carries out the whole network wide area simultaneous synthesis point
Analysis, so that it may tell the oscillation of power which frequency mode is occurring in operations staff's power grid, the coverage area of oscillation is how,
Hunting power is in the interactive relation of equipment room, situations such as each equipment is to the percentage contribution of hunting power, to help operations staff
The method found out inhibition and eliminate oscillation of power.
Realize the physical function module map of the device of the invention as shown in figure 3, original voltage and current analog signal enters
After device, voltage and current digital signal is converted by sampling hold circuit and A/D conversion circuit, sampling pulse used comes from
In high-precision clock and frequency divider, to provide high-precision markers for digital signal.Digital voltage current signal is through counting
The three-phase instantaneous active power signal of number is obtained after calculation, which is separately input to the corresponding physical treatment of each frequency range simultaneously
In module；In the physical treatment module of each frequency range, instantaneous power signal is after filtering and down-sampled rate is handled, with Tx2For meter
The period is calculated, with 2Tx1To calculate data window duration, the frequency spectrum vibrated and initial phase analysis, and drops analysis result reckoning and arrive
The non-boundary whole scale moment；The analysis result of last each frequency range is stored on the spot, and by sending scheduling station in communication system.This
The algorithm flow chart of inventive method as shown in Figure 1, the specific steps of which are as follows:
Step 1: measuring to obtain power grid by the three-phase instantaneous voltage of grid branch measurement point and the recording of transient current corresponding
The total instantaneous active power curve of three-phase of tested branch, and whole moment second markers, markers precision are marked using clock synchronization apparatus
Higher than 1 microsecond, original instantaneous recording curve sample rate at least 4800Hz or more.
pr=uaia+ubib+ucic
Wherein ua、ub、ucThe three-phase voltage of the tested branch of certain respectively directly measured from voltage transformer is instantaneous
Value；Ia, ib, ic are respectively the three-phase current instantaneous value of certain the tested branch directly measured from current transformer, prIt is three
Mutually total instantaneous active power；
Step 2: by the total instantaneous active power curve p of above-mentioned three-phaserIt is distributed to the memory of the processing module of 5 frequency ranges in real time
In buffer area, f is usedx1Indicate the initial frequency of frequency range X, the period of the corresponding initial frequency of frequency range X is Tx1=1/fx1；fx2It indicates
The termination frequency of frequency range X, the corresponding period for terminating frequency are T2X=1/fx2；Each frequency range buffer area can at least store TX1× 4 seconds
Data；Wherein, x value 1-5 respectively indicates 1-5 different frequency range；The frequency range of frequency range 1 to frequency range 5 is to be sequentially increased
's；This 5 different frequency ranges are provided that
Frequency range 1 is load oscillation or the oscillation of power plant heat power engineering system) frequency range: f21-f22Range be (0.01Hz-
0.1Hz]；
Frequency range 2 is electric system section low-frequency oscillation frequency range: f21-f22Range be (0.1Hz -1Hz]；
Frequency range 3 is oscillation or generator speed-adjusting and control system oscillation frequency bands between generator: f31-f32Range be (1Hz-
10H]；
Frequency range 4 is subsynchronous or supersynchronous oscillation frequency bands: f41-f42Range be (10Hz -100Hz]；
Frequency range 5 is each secondary industrial frequency resonance frequency range: f51-f52Range be (100Hz -1000Hz]；
Step 3: using recursion bandpass filter in real time to the three-phase for entering core buffer in each frequency range processing module
Total instantaneous power data are filtered；- 3dB the free transmission range of each bandpass filter is [fx1, fx2], see Fig. 2.
Step 4: in each frequency range processing module to the total instantaneous power data of filtered power grid three-phase using at down-sampled rate
It manages (sample rate of original transient data is at least more than 4800Hz).The specific method is as follows:
In the processing module of frequency range 1, are retained to the total instantaneous power data of filtered power grid three-phase every 1 second whole moment second
Transient data；
In the processing module of frequency range 2, to the total instantaneous power data of filtered power grid three-phase, every N1 point retains newest
1 point, wherein N1 is 1/1000 smallest positive integral less than or equal to former sample rate；
In the processing module of frequency range 3, instantaneous power data total for filtered power grid three-phase, down-sampled sample rate
Become the 1/100 of former sample rate；
In the processing module of frequency range 4, instantaneous power data total for filtered power grid three-phase, down-sampled sample rate
Become the 1/10 of former sample rate；
In the processing module of frequency range 5, down-sampled sample rate is taken as the smallest crude sampling rate more than or equal to 4800
Approximate number.
Step 5: with the week of initial frequency corresponding to this frequency range in the processing module of remaining each frequency range in addition to frequency range 5
2 times of values of phase are 2Tx1To calculate data window duration, with the cycle T of termination frequency corresponding to this frequency rangex2For calculating cycle
The moving step length for calculating data window calculates number to the section using frequency spectrum analysis methods such as Prony or Empirical mode decompositions
Spectrum analysis is carried out according to the total instantaneous power data of grid branch three-phase in window；The processing module of frequency range 5 was then with 0.05 second
Data window is calculated, 0.01 second is calculating cycle, carries out spectrum analysis to the total instantaneous power data of grid branch three-phase；It is analyzing
As a result take hunting power amplitude to be greater than given threshold in (always instantaneously has in 30 second period being taken as before current calculation window starts
The 5% of the average value of maximum point and smallest point absolute value after function power curve removal dc component), and energy accounting is greater than
5% preceding 5 frequency modes are as primary Calculation as a result, calculated result information includes: the frequency for the hunting power that each frequency range includes
Rate fx, amplitude Ax1, initial phase φx1, decay factor σx, can accounting ex；
According to the definition in the above method and step 2 to each band limits, it is known that the specific calculating parameter of each frequency range is such as
Under:
Frequency range 1: calculating data window duration 200 seconds, and calculating cycle is 10 seconds
Frequency range 2: calculating data window duration 20 seconds, and calculating cycle is 1 second
Frequency range 3: calculating data window duration 2 seconds, and calculating cycle is 0.1 second
Frequency range 4: calculating data window duration 0.2 second, and calculating cycle is 0.01 second
Frequency range 5: calculating data window duration 0.05 second, and calculating cycle is 0.01 second
Step 6: each frequency range chooses t respectivelyX2-Tx1The nearest following moment before/2 moment is the calculating for needing to export
As a result t at the time ofx-out, wherein tx2The finish time of data window is calculated for current frequency range X, and by the calculating of the calculating data window
As a result by when window start time tx1It is converted to tx-outMoment；
For frequency range 1, t is chosen12-T11At the time of whole 10 second moment is the calculated result for needing to export before/2 moment
t1-out；
For frequency range 2, t is chosen22-T21T at the time of whole moment second is the calculated result for needing to export before/2 moment2-out；
For frequency range 3, t is chosen32-T31At the time of whole 0.1 second moment is the calculated result for needing to export before/2 moment
t3-out；
For frequency range 4, t is chosen42-T41At the time of whole 0.01 second moment is the calculated result for needing to export before/2 moment
t4-out；
For frequency range 5, t is chosen52-T51At the time of whole 0.01 second moment is the calculated result for needing to export before/2 moment
t5-out；
For each frequency mode that above-mentioned each frequency range is calculated and selected by step 5, by calculated result by calculating data window
Start time tx1, it is converted to tx-outMoment needs tx1The amplitude A of the oscillation of power at momentx1With initial phase φx1, modification
For tx-outThe amplitude A at momentx-outWith initial phase φx-out, initial phase is relative to 2 π f of moment sinxThe phase of t curve, tool
Body calculation formula is as follows:
φx-out=φx1+(tx-out-tx1)(2πfx)
Remaining information, i.e. the frequency f of hunting powerx, decay factor σx, can accounting exWhen still taking calculating data window starting point
Carve tx1When value；
Step 7: storing on the spot that each frequency range calculates in real time as a result, retaining at least 2 weeks；
Step 8: by the real-time Shang Song main website of the calculated result of each frequency range or according to main website to the interval for sending data in each frequency range
Sent on it is required that, frequency, amplitude, initial phase, decay factor including the oscillation of power for meeting threshold value detected in each frequency range with
And markers；
Step 9: each frequency range processing module waiting meets this frequency range calculating cycle Tx2After duration, return step 5.
Applicant combines Figure of description to be described in detail and describe the embodiment of the present invention, but this field skill
Art personnel are it should be understood that above embodiments are only the preferred embodiments of the invention, and explanation is intended merely to help reader in detail
More fully understand spirit of that invention, and it is not intended to limit the protection scope of the present invention, on the contrary, any based on invention essence of the invention
Any improvement or modification made by mind should all be fallen within the scope and spirit of the invention.
Claims (9)
1. a kind of power grid entire spectrum oscillation of power phasor method for synchronously measuring based on three-phase instantaneous power, which is characterized in that institute
State method the following steps are included:
Step 1: measuring to obtain power grid by the three-phase instantaneous voltage of grid branch measurement point and the recording of transient current and be accordingly tested
The total instantaneous active power curve of the three-phase of branch, and whole moment second markers is marked using clock synchronization apparatus；
pr=uaia+ubib+ucic
Wherein, ua、ub、ucThe three-phase voltage instantaneous value of the tested branch of certain respectively directly measured from voltage transformer；
Ia, ib, ic are respectively the three-phase current instantaneous value of certain the tested branch directly measured from current transformer, PrIt is total for three-phase
Instantaneous active power；
Step 2: by the total instantaneous active power P of above-mentioned three-phaserIt is distributed to the core buffer of 5 different frequency range processing modules in real time
In, use fx1Indicate the initial frequency of frequency range X, the period of the corresponding initial frequency of frequency range X is Tx1=1/fx1；fx2Indicate frequency range X's
Frequency is terminated, the corresponding period for terminating frequency is Tx2=1/fx2；Each frequency range buffer area can at least store Tx1× 4 seconds data；
Wherein, X value 1-5 respectively indicates 1-5 different frequency range；The frequency range of frequency range 1 to frequency range 5 is sequentially increased；
Step 3: using recursion bandpass filter in real time to the three-phase total wink for entering core buffer in each frequency range processing module
When active power data be filtered；
Step 4: the total instantaneous power data of filtered three-phase being handled using down-sampled rate in each frequency range processing module；
Step 5: with the period of initial frequency corresponding to this frequency range in the processing module of remaining each frequency range in addition to frequency range 5
2 times of values are 2Tx1To calculate data window duration, with the cycle T of termination frequency corresponding to this frequency rangex2It is counted for calculating cycle
The moving step length for calculating data window calculates this section the total instantaneous function of grid branch three-phase of data window using frequency spectrum analysis method
Rate data carry out spectrum analysis；For the processing module of frequency range 5 then with 0.05 second to calculate data window, 0.01 second is calculating cycle,
Spectrum analysis is carried out to the total instantaneous power data of grid branch three-phase in window at this section using frequency spectrum analysis method；
Hunting power amplitude is taken to be greater than given threshold in spectrum analysis, and can preceding 5 frequency mode work of the accounting greater than 5%
For calculated result, that is, result of spectrum analysis；
Step 6: each frequency range chooses t respectivelyx2-Tx1The nearest following moment before/2 moment is the calculated result for needing to export
Moment tx-out, wherein tx2For current frequency range X calculate data window finish time, and by the calculated result of the calculating data window by
When window start time tx1It is converted to tx-outMoment；
For frequency range 1, t is chosen12-T11T at the time of whole 10 second moment is the calculated result for needing to export before/2 moment1-out；
For frequency range 2, t is chosen22-T21T at the time of whole moment second is the calculated result for needing to export before/2 moment2-out；
For frequency range 3, t is chosen32-T31T at the time of whole 0.1 second moment is the calculated result for needing to export before/2 moment3-out；
For frequency range 4, t is chosen42-T41T at the time of whole 0.01 second moment is the calculated result for needing to export before/2 moment4-out；
For frequency range 5, t is chosen52-T51T at the time of whole 0.01 second moment is the calculated result for needing to export before/2 moment5-out；
Step 7: storing on the spot that each frequency range calculates in real time as a result, retaining at least 2 weeks；
Step 8: by the real-time Shang Song main website of the calculated result of each frequency range or space requirement according to main website to data are sent in each frequency range
On give, frequency, amplitude, initial phase, decay factor including the oscillation of power for meeting threshold value detected in each frequency range are with timely
Mark；
Step 9: each frequency range processing module waiting meets this frequency range calculating cycle Tx2After duration, return step 5.
2. the power grid entire spectrum oscillation of power phasor synchro measure side according to claim 1 based on three-phase instantaneous power
Method, it is characterised in that:
In step 1, markers precision is higher than 1 microsecond, original instantaneous recording curve sample rate at least 4800Hz or more.
3. the power grid entire spectrum oscillation of power phasor synchro measure side according to claim 1 based on three-phase instantaneous power
Method, it is characterised in that:
In step 2,5 different frequency ranges are provided that
Frequency range 1 is load oscillation or power plant heat power engineering system oscillation frequency bands: f11-f12Range be (0.01Hz, 0.1Hz]；
Frequency range 2 is electric system section low-frequency oscillation frequency range: f21-f22Range be (0.1Hz, 1Hz]；
Frequency range 3 is oscillation or generator speed-adjusting and control system oscillation frequency bands between generator: f31-f32Range be (1Hz, 10Hz]；
Frequency range 4 is subsynchronous or supersynchronous oscillation frequency bands: f41-f42Range be (10Hz, 100Hz]；
Frequency range 5 is each secondary industrial frequency resonance frequency range: f51-f52Range be (100Hz, 1000Hz].
4. the power grid entire spectrum oscillation of power phasor synchro measure side according to claim 1 based on three-phase instantaneous power
Method, it is characterised in that:
In step 3, -3dB the free transmission range of recursion bandpass filter corresponding to each frequency range processing module is [fx1, fx2]。
5. the power grid entire spectrum oscillation of power phasor synchro measure according to claim 1 or 3 based on three-phase instantaneous power
Method, it is characterised in that:
In step 4, whole to the reservation of filtered power grid three-phase total instantaneous power data every 1 second in the processing module of frequency range 1
Second moment transient data；
In the processing module of frequency range 2, to the total instantaneous power data of filtered power grid three-phase, every N1 point retains newest 1
Point, wherein N1 is 1/1000 smallest positive integral less than or equal to former sample rate；
In the processing module of frequency range 3, instantaneous power data total for filtered power grid three-phase, down-sampled sample rate becomes
The 1/100 of former sample rate；
In the processing module of frequency range 4, instantaneous power data total for filtered power grid three-phase, down-sampled sample rate becomes
The 1/10 of former sample rate；
In the processing module of frequency range 5, down-sampled sample rate is taken as the smallest crude sampling rate approximate number more than or equal to 4800.
6. the power grid entire spectrum oscillation of power phasor synchro measure side according to claim 1 based on three-phase instantaneous power
Method, it is characterised in that:
In steps of 5, the frequency spectrum analysis method is Prony algorithm or EMD algorithm.
7. the power grid entire spectrum oscillation of power phasor synchro measure according to claim 1 or 6 based on three-phase instantaneous power
Method, it is characterised in that:
In steps of 5, the given threshold is taken as total instantaneous active power in 30 second period before current calculation window starts
The 5% of the average value of maximum point and smallest point absolute value after curve removal dc component.
8. the power grid entire spectrum oscillation of power phasor synchro measure side according to claim 1 based on three-phase instantaneous power
Method, it is characterised in that:
In steps of 5, calculated result information includes: the frequency f for the hunting power that each frequency range includesx, amplitude Ax1, initial phase φx1、
Decay factor σx, can accounting ex。
9. the power grid entire spectrum oscillation of power phasor synchro measure side according to claim 8 based on three-phase instantaneous power
Method, it is characterised in that:
In step 6, each frequency mode for each frequency range being calculated and being selected by step 5, by tx1The oscillation of power at moment
Amplitude Ax1With initial phase φx1, it is scaled tx-outThe amplitude A at momentx-outWith initial phase φx-out, when initial phase is relative to this
Carve sin2 π fxThe phase of t curve, specific formula is as follows:
φx-out=φx1+(tx-out-tx1)(2πfx)
Remaining information, i.e. the frequency f of hunting powerx, decay factor σx, can accounting exStill it takes and calculates data window start time
tx1When value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611249360.0A CN106645919B (en) | 2016-12-26 | 2016-12-26 | The full spectral power of power grid based on three-phase instantaneous power vibrates phasor method for synchronously measuring |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611249360.0A CN106645919B (en) | 2016-12-26 | 2016-12-26 | The full spectral power of power grid based on three-phase instantaneous power vibrates phasor method for synchronously measuring |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106645919A CN106645919A (en) | 2017-05-10 |
| CN106645919B true CN106645919B (en) | 2019-11-12 |
Family
ID=58836539
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611249360.0A Active CN106645919B (en) | 2016-12-26 | 2016-12-26 | The full spectral power of power grid based on three-phase instantaneous power vibrates phasor method for synchronously measuring |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106645919B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107179440B (en) * | 2017-05-15 | 2019-12-10 | 国网新疆电力公司 | Subsynchronous oscillation inter-harmonic extraction method for online self-adaption frequency variation |
| CN108362931A (en) * | 2017-12-25 | 2018-08-03 | 北京四方继保自动化股份有限公司 | A kind of synchronous data collection method applied to network stability control |
| CN108957129A (en) * | 2018-07-27 | 2018-12-07 | 中国南方电网有限责任公司 | A kind of new energy grid connection system time/supersynchronous oscillation self-adaptive monitoring method |
| CN110378020B (en) * | 2019-07-19 | 2020-04-10 | 华北电力大学 | Multi-frequency-band dynamic phasor electromagnetic transient simulation method and system for power grid commutation converter |
| CN111446890B (en) * | 2020-03-25 | 2022-02-11 | 武汉迈信电气技术有限公司 | Encoder with current detection function and current detection method |
| CN112600638B (en) * | 2020-12-15 | 2023-07-14 | 国网北京市电力公司 | Power grid clock synchronization method and device |
| CN113960363B (en) * | 2021-09-22 | 2023-09-15 | 国网新疆电力有限公司 | High-frequency oscillation phasor measurement method and device based on multicomponent Morlet wavelet |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1093171A (en) * | 1994-01-31 | 1994-10-05 | 清华大学 | The method of monitoring subsynchronous oscillation of electrical power system and monitor |
| JPH1052096A (en) * | 1996-07-31 | 1998-02-20 | Tokyo Electric Power Co Inc:The | Power system stabilizing device and excitation control system |
| CN101408567A (en) * | 2008-11-28 | 2009-04-15 | 北京四方继保自动化股份有限公司 | Large scale electric network low-frequency oscillation frequency division section detection method based on empirical mode decomposition |
| CN101924364A (en) * | 2010-07-23 | 2010-12-22 | 清华大学 | Method for estimating nonlinear state of substation-dispatching center two-stage distributed power grid |
-
2016
- 2016-12-26 CN CN201611249360.0A patent/CN106645919B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1093171A (en) * | 1994-01-31 | 1994-10-05 | 清华大学 | The method of monitoring subsynchronous oscillation of electrical power system and monitor |
| JPH1052096A (en) * | 1996-07-31 | 1998-02-20 | Tokyo Electric Power Co Inc:The | Power system stabilizing device and excitation control system |
| CN101408567A (en) * | 2008-11-28 | 2009-04-15 | 北京四方继保自动化股份有限公司 | Large scale electric network low-frequency oscillation frequency division section detection method based on empirical mode decomposition |
| CN101924364A (en) * | 2010-07-23 | 2010-12-22 | 清华大学 | Method for estimating nonlinear state of substation-dispatching center two-stage distributed power grid |
Non-Patent Citations (1)
| Title |
|---|
| 基于多个广域测量系统的联合低频振荡分析和监测;段刚等;《中国电力》;20150430;第48卷(第4期);正文第102页右栏第4段-第103页左栏第5段 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106645919A (en) | 2017-05-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106645919B (en) | The full spectral power of power grid based on three-phase instantaneous power vibrates phasor method for synchronously measuring | |
| Reza et al. | Accurate estimation of single-phase grid voltage parameters under distorted conditions | |
| CN105223418B (en) | The measurement method and measuring device of subsynchronous and supersynchronous harmonic phasor | |
| RU2437105C2 (en) | Sensor, method and system for monitoring power lines | |
| Roscoe | Exploring the relative performance of frequency-tracking and fixed-filter phasor measurement unit algorithms under C37. 118 test procedures, the effects of interharmonics, and initial attempts at merging P-class response with M-class filtering | |
| CN106841778B (en) | PMU-based subsynchronous and supersynchronous harmonic parameter processing method | |
| Dotta et al. | A matlab-based PMU simulator | |
| Zhan et al. | Universal grid analyzer design and development | |
| CN107851995A (en) | For detecting the method and apparatus for supplying the voltage in power network | |
| Jain et al. | An adaptive time-efficient technique for harmonic estimation of nonstationary signals | |
| JP2013044752A (en) | Phase identification system and method | |
| CN109061297B (en) | Sub/super synchronous oscillation on-line monitoring method and device | |
| CN112505407B (en) | Power grid broadband oscillation monitoring method, system and equipment and readable storage medium | |
| US10948527B2 (en) | Power system detection of sustained oscillations based on synchrophasor frequency measurements | |
| Zhang et al. | Shifting window average method for accurate frequency measurement in power systems | |
| Mokeev | Filter synthesis for PMU | |
| Vanfretti et al. | A PMU-based fast real-time sub-synchronous oscillation detection application | |
| Vanfretti et al. | Estimation of Eastern Denmark's electromechanical modes from ambient phasor measurement data | |
| CN109116109B (en) | Inter-harmonic online monitoring method and device | |
| CN104007408A (en) | Method and device for on-line detection of dynamic performance of PMU | |
| CN111257636B (en) | Subsynchronous and supersynchronous oscillation frequency identification method for power system | |
| CN117630541A (en) | Broadband oscillation monitoring method and device | |
| Wu et al. | Estimation of power system inertia using traditional swing equation, polynomial approximation and RV methods | |
| Bronzini et al. | Power system modal identification via wavelet analysis | |
| CN111273103B (en) | Electric power system oscillation identification method based on synchronous phasor complex field spectrum analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |