CN109884390A - A kind of subsynchronous/supersynchronous phasor measurement method of power system voltage harmonic wave - Google Patents
A kind of subsynchronous/supersynchronous phasor measurement method of power system voltage harmonic wave Download PDFInfo
- Publication number
- CN109884390A CN109884390A CN201910054291.5A CN201910054291A CN109884390A CN 109884390 A CN109884390 A CN 109884390A CN 201910054291 A CN201910054291 A CN 201910054291A CN 109884390 A CN109884390 A CN 109884390A
- Authority
- CN
- China
- Prior art keywords
- frequency
- amplitude
- subsynchronous
- supersynchronous
- spectral line
- 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.)
- Withdrawn
Links
Landscapes
- Measurement Of Resistance Or Impedance (AREA)
Abstract
A kind of subsynchronous/supersynchronous phasor measurement method of power system voltage harmonic wave carries out the processing of time domain Kaiser window to obtained voltage sample data comprising steps of 1) being acquired line number Data preprocess of going forward side by side to mains voltage signal；2) the local Discrete Fourier transformation that the data after adding window are successively carried out with frequency-division section, calculates the spectral line amplitude in each frequency range；3) it counts and arranges the spectral line amplitude in addition to fundamental frequency in each frequency range, obtain the spectral line of amplitude and specific；4) to obtained spectral line amplitude, double spectral line interpolation is carried out between two neighboring maximum amplitude, frequency correction coefficients are calculated, so that it is determined that the frequency of subsynchronous and supersynchronous component；5) frequency correction coefficients are utilized, the correction formula for acquiring amplitude and phase is calculated, obtain the amplitude and phase of subsynchronous and supersynchronous component.The present invention carries out Kaiser window function processing to sampled data in frequency range at equal intervals, is corrected processing to the particular spectral lines of acquisition, obtains subsynchronous and hypersynchronous frequency, amplitude and phase.
Description
Technical field
The present invention relates to a kind of subsynchronous/supersynchronous phasor measurement methods of power system voltage harmonic wave.
Background technique
In recent years, as large-scale wind power sends the continuous development with DC transmission engineering, the electric power that may induce outside
The problems such as system is subsynchronous and supersynchronous oscillation gradually highlights.Be mainly shown as in power grid the subsynchronous of diverging or permanent width and/
Or supersynchronous electric current, voltage and power harmonics, and frequency of oscillation can change with the change of grid operation mode.This is by shadow
The safe operation for ringing power grid and equipment results even in serious stability accident or equipment damage.Therefore it needs to secondary (super) same
It walks harmonic wave and carries out dynamic measurement.
Existing phasor measurement unit (Phasor Measurement Unit, PMU) is that do not have using power frequency phasor as target
There is (super) the synchronous policy requirement of special consideration time, and usually filters out the harmonic component in signal to improve power frequency phasor
Tactful precision.This not can accurately reflect the dynamic of time (super) synchronized harmonics necessarily, is not able to satisfy measurement request.PMU algorithm institute base
Fractional harmoni is not considered generally in signal model, it, not only cannot be quasi- once actual signal includes time (super) synchronization signal components
Time (super) synchronized phasor is really detected, and fundamental phasors measurement result can be seriously affected.And the accuracy of secondary (super) synchronized phasor
The positioning of harmonic wave will be directly affected, brings great difficulty to harmonics restraint.Kaiser window function is the one kind with adjustability
Harmonic amplitude and frequency can be improved by adjusting the main lobe width of window spectrum and the ratio of side lobe attenuation in zeroth order prolate spheroidal wave functions
Detection accuracy.But adjustment parameter selection is excessive, and frequency resolution can reduce, detection accuracy decline.
Summary of the invention
For above-mentioned technological deficiency present in the prior art, time (super) synchronized phasor can not be accurately detected in PMU, this
Invention provides a kind of subsynchronous/supersynchronous phasor measurement method of power system voltage harmonic wave, can accurately detect secondary (super) synchronous
Phasor, and adjustment parameter selection is less, will not reduce frequency resolution, detection accuracy is met the requirements.
Above-mentioned technical problem is solved, the following technical solution is employed by the present invention:
A kind of subsynchronous/supersynchronous phasor measurement method of power system voltage harmonic wave, it is characterized in that comprising the steps of:
(1) acquisition obtains the voltage sample data y (n) of 50 continuous power frequency periods of PMU installation point；N=0 in formula,
1 ... ..., N-1；Wherein, n is integer, indicates current sampling point；N is positive integer, indicates the total length of institute's sampled data；
(2) time-domain windowed processing is carried out using Kaiser window function to the voltage sample data y (n) that step (1) obtains；
The time-domain expression of Kaiser window function w (n) are as follows:
Wherein, I0(β) is Bessel function；N=0,1,2 ..., N/2；β is adjustable parameter, to adjust main lobe width and
Sidelobe level attenuation ratio；β is bigger, and with regard to smaller, main lobe width can then be increase accordingly secondary lobe；
The voltage sample data y (n) that step (1) obtains is multiplied with the discrete series of Kaiser window function, obtains adding window
Treated data x (n), shown in formula specific as follows:
X (n)=y (n) (n)；
(3) it is equally spaced frequency range by the frequency partition in subsynchronous and supersynchronous range, adds to what is obtained from step (2)
Window treated data x (n) carries out the local Discrete Fourier transformation of each frequency range, calculates the peak in all frequency ranges in addition to fundamental frequency
It is worth frequency point spectral line component, the calculation formula of local Discrete Fourier transformation are as follows:
Wherein,K is integer, herein k=1,2 ..., that is, calculate the corresponding n-th k root of subsynchronous frequency
Spectral line component；J is imaginary unit, and e is the truth of a matter of natural exponential function；It is the n-th 1 spectral line components, is plural form,
It takes absolute value respectively and the amplitude of the n-th 1 spectral lines can be obtained:
(4) it is X (n1) that the n-th 1 peak value frequency point spectral lines that step (3) is calculated, which correspond to amplitude, and nearby sampling obtains
Amplitude maximum and time greatly be respectively X (nmax1) and X (nmax2)；Enable α=nm-nmax1- 0.5, in maximum value X (nmax1) and X
(nmax2), between carry out interpolation, the expression formula of interpolation are as follows:
Wherein λ is scale parameter, and w (2 π f) is Kaiser window function；Due to X (nmax1) and X (nmax2) be it is known, can be with
Calculate λ value；
(5) λ-expression obtained to step (4) carries out fitting of a polynomial using the method for curve matching and approaches, can obtain
α=F (λ)
Then the frequency of the n-th 1 spectral lines is fn1=(+n1+0.5) Δ f；
Wherein, resolution ratio
(6) spectral line is weighted and averaged, amplitude correction formula can be obtained:
(7) phase only pupil filter formula are as follows:
(8) frequency, phase and amplitude of the nk with spectral line are successively sought.
The present invention is acquired voltage signal and pre-processes first, then uses the method frequency-division section of adding window spectral line interpolation
The correcting value of subsynchronous and supersynchronous phasor is calculated, finally calculates and obtains subsynchronous and supersynchronous phasor frequency, amplitude and first phase
Position.
The present invention can obtain compared with the existing technology it is following the utility model has the advantages that
Eliminate the leakage and fence effect of frequency spectrum to a certain extent using the discrete Fourier transform of Kaiser window interpolation
It answers, can be improved the measurement accuracy of frequency；Using the spectral line phase value feature unrelated with frequency departure, subsynchronous phasor phase is estimated
The method of position is simple and easy.
Specific embodiment
Technical solution of the present invention is further stated in detail below.
Subsynchronous phasor measurement embodiment of the method for the invention, comprising the following steps:
(1) obtain voltage signal data: according to IEC 61000-4-7 standard, acquisition obtains 50 of PMU installation point continuously
The voltage sample data y (n) of power frequency period；N=0 in formula, 1 ... ..., N-1, wherein n is integer, indicates current sampling point；N
For positive integer, the total length of institute's sampled data is indicated；
(2) time-domain windowed is handled:
Time-domain windowed processing is carried out using Kaiser window function to the voltage sample data y (n) that step (1) obtains,
The time-domain expression of Kaiser window function w (n) are as follows:
Wherein, I0(β) is Bessel function, n=0,1,2 ..., N/2.β is adjustable parameter, adjustable main lobe width and
Sidelobe level attenuation ratio；β is bigger, and with regard to smaller, main lobe width can then be increase accordingly secondary lobe；
The voltage sample data y (n) that step (1) obtains is multiplied with the discrete series of Kaiser window function, can be obtained
Data x (n) after windowing process takes β=12 here, shown in formula specific as follows:
X (n)=Y (n) W (n)
3, local Discrete Fourier transformation:
It is equally spaced frequency range by the frequency partition in subsynchronous and supersynchronous range, to the adding window obtained from step (2)
Treated, and data x (n) carries out the local Discrete Fourier transformation of each frequency range, calculates the peak value in all frequency ranges in addition to fundamental frequency
Frequency point spectral line component, the calculation formula of local Discrete Fourier transformation are as follows:
Wherein,K is integer, herein k=1,2 ..., that is, calculate the corresponding n-th k root of subsynchronous frequency
Spectral line component；J is imaginary unit, and e is the truth of a matter of natural exponential function；It is the n-th 1 spectral line components, is plural form,
It takes absolute value respectively and the amplitude of the n-th 1 spectral lines can be obtained:
(4) it is X (n1) that the n-th 1 peak value frequency point spectral lines that step (3) is calculated, which correspond to amplitude, and nearby sampling obtains
Amplitude maximum and time greatly be respectively X (nmax1) and X (nmax2), enable α=nm-nmax1- 0.5, in maximum value X (nmax1) and X
(nmax2), between carry out interpolation, the expression formula of interpolation are as follows:
Wherein λ is scale parameter, and w (2 π f) is Kaiser window function.Due to X (nmax1) and X (nmax2) be it is known, can be with
Calculate λ value；
(5) λ-expression obtained to step (4) carries out fitting of a polynomial using the method for curve matching and approaches, can obtain
α=F (λ)
Then the frequency of the n-th 1 spectral lines is fn1=(α+n1+0.5) Δ f.
Wherein, Δ f is resolution ratio.
(6) spectral line is weighted and averaged, amplitude correction formula can be obtained:
(7) phase only pupil filter formula are as follows:
(8) the n-th k root spectral line is successively asked to correspond to frequency, the phase and amplitude of subsynchronous or supersynchronous ingredient.
Claims (5)
1. a kind of subsynchronous/supersynchronous phasor measurement method of power system voltage harmonic wave, it is characterised in that comprise the steps of:
1) it is acquired line number Data preprocess of going forward side by side to mains voltage signal, time domain is carried out to obtained voltage sample data
The processing of Kaiser window；
2) the local Discrete Fourier transformation that the data after adding window are successively carried out with frequency-division section, calculates the breadth of spectral line in each frequency range
Value；
3) it counts and arranges the spectral line amplitude in addition to fundamental frequency in each frequency range, obtain the spectral line of amplitude and specific；
4) to obtained spectral line amplitude, double spectral line interpolation is carried out between two neighboring maximum amplitude, frequency correction system is calculated
Number, so that it is determined that the frequency of subsynchronous and supersynchronous component；
5) frequency correction coefficients are utilized, the correction formula for acquiring amplitude and phase is calculated, obtain subsynchronous and supersynchronous component
Amplitude and phase.
2. subsynchronous/supersynchronous phasor measurement method of power system voltage harmonic wave according to claim 1, feature exist
In: the step 1) specifically:
Acquisition obtains the voltage sample data y (n) of 50 continuous power frequency periods of PMU installation point；N=0 in formula, 1 ... ..., N-
1；Wherein, n is integer, indicates current sampling point；N is positive integer, indicates the total length of institute's sampled data；
Time-domain windowed processing, Kaiser window function w (n) are carried out using Kaiser window function to obtained voltage sample data y (n)
Time-domain expression are as follows:
Wherein, I0(β) is Bessel function, n=0,1,2 ..., N/2；β is adjustable parameter, to adjust main lobe width and secondary lobe
Level attenuation ratio；β is bigger, and with regard to smaller, main lobe width can then be increase accordingly secondary lobe；
Obtained voltage sample data y (n) is multiplied with the discrete series of Kaiser window function, after windowing process can be obtained
Data X (n), shown in formula specific as follows:
X (n)=y (n) W (n).
3. subsynchronous/supersynchronous phasor measurement method of power system voltage harmonic wave according to claim 1, feature exist
In: the step 2) specifically:
It is equally spaced frequency range by the frequency partition in subsynchronous and supersynchronous range, the windowing process to being obtained from step (2)
Data x (n) afterwards carries out the local Discrete Fourier transformation of each frequency range, calculates the peak value frequency point in all frequency ranges in addition to fundamental frequency
Spectral line component, the calculation formula of local Discrete Fourier transformation are as follows:
Wherein,K is integer, herein k=1,2 ..., that is, calculate the corresponding n-th k root spectral line of subsynchronous frequency
Component；J is imaginary unit, and e is the truth of a matter of natural exponential function；It is the n-th 1 spectral line components, is plural form, respectively
It takes absolute value and the amplitude of the n-th 1 spectral lines can be obtained:
4. subsynchronous/supersynchronous phasor measurement method of power system voltage harmonic wave according to claim 1, feature exist
In: the step 3) and step 4) specifically:
It is X (n1) that the n-th 1 peak value frequency point spectral lines being calculated, which correspond to amplitude, the obtained amplitude maximum of nearby sample with it is secondary
It is greatly respectively X (nmax1) and X (nmax2)；
Enable α=nm-nmax1- 0.5, in maximum value X (nmax1) and X (nmax2), between carry out interpolation, the expression formula of interpolation are as follows:
Wherein λ is scale parameter, and w (2 π f) is Kaiser window function；Due to X (nmax1) and X (nmax2) be it is known, can calculate
λ value out；
To obtained λ-expression, fitting of a polynomial is carried out using the method for curve matching and is approached, can be obtained
α=F (λ)
Then the frequency of the n-th 1 spectral lines is fn1=(α+n1+0.5) Δ f；
Wherein, Δ f is resolution ratio.
5. subsynchronous/supersynchronous phasor measurement method of power system voltage harmonic wave according to claim 1, feature exist
In: the step 5) specifically:
Spectral line is weighted and averaged, amplitude correction formula can be obtained:
Phase only pupil filter formula are as follows:
The n-th k root spectral line is successively asked to correspond to frequency, the phase and amplitude of subsynchronous or supersynchronous ingredient.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910054291.5A CN109884390A (en) | 2019-01-21 | 2019-01-21 | A kind of subsynchronous/supersynchronous phasor measurement method of power system voltage harmonic wave |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910054291.5A CN109884390A (en) | 2019-01-21 | 2019-01-21 | A kind of subsynchronous/supersynchronous phasor measurement method of power system voltage harmonic wave |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109884390A true CN109884390A (en) | 2019-06-14 |
Family
ID=66926389
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910054291.5A Withdrawn CN109884390A (en) | 2019-01-21 | 2019-01-21 | A kind of subsynchronous/supersynchronous phasor measurement method of power system voltage harmonic wave |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109884390A (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111294018A (en) * | 2020-03-17 | 2020-06-16 | 厦门傅里叶电子有限公司 | LRA drive pulse waveform design method based on Kaiser window |
| CN111650436A (en) * | 2020-05-28 | 2020-09-11 | 电子科技大学 | Subsynchronous oscillation identification method for high-proportion renewable energy power system |
| CN112748284A (en) * | 2020-12-30 | 2021-05-04 | 清华大学 | Method and device for measuring synchronous waveform and broad-spectrum phasor of power system |
| CN117892065A (en) * | 2024-03-14 | 2024-04-16 | 江西飞尚科技有限公司 | Vibrating wire sensor signal correction method, system, computer and storage medium |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN201352236Y (en) * | 2008-12-02 | 2009-11-25 | 湖南海兴电器有限责任公司 | Harmonic power meter based on Kaiser window double-spectrum-line interpolation FFT |
| CN103575984A (en) * | 2012-08-02 | 2014-02-12 | 西安元朔科技有限公司 | Harmonic analysis method based on Kaiser window double-spectral-line interpolation FFT |
| CN105182073A (en) * | 2015-09-06 | 2015-12-23 | 中国电力科学研究院 | Dynamic phasor measurement system of phasor measurement unit (PMU), and measurement method thereof |
| CN105242225A (en) * | 2015-09-06 | 2016-01-13 | 中国电力科学研究院 | Dynamic phasor measurement device calibration system and calibration method thereof |
| US20160329709A1 (en) * | 2015-05-07 | 2016-11-10 | Industry-Academic Cooperaiton Foundation, Yonsei University | Inverter-based stand-alone microgrid control system using pmu |
| CN107561359A (en) * | 2017-08-31 | 2018-01-09 | 国网新疆电力公司 | A kind of m-Acetyl chlorophosphonazo extracting method based on Blackman window |
| CN108957129A (en) * | 2018-07-27 | 2018-12-07 | 中国南方电网有限责任公司 | A kind of new energy grid connection system time/supersynchronous oscillation self-adaptive monitoring method |
-
2019
- 2019-01-21 CN CN201910054291.5A patent/CN109884390A/en not_active Withdrawn
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN201352236Y (en) * | 2008-12-02 | 2009-11-25 | 湖南海兴电器有限责任公司 | Harmonic power meter based on Kaiser window double-spectrum-line interpolation FFT |
| CN103575984A (en) * | 2012-08-02 | 2014-02-12 | 西安元朔科技有限公司 | Harmonic analysis method based on Kaiser window double-spectral-line interpolation FFT |
| US20160329709A1 (en) * | 2015-05-07 | 2016-11-10 | Industry-Academic Cooperaiton Foundation, Yonsei University | Inverter-based stand-alone microgrid control system using pmu |
| CN105182073A (en) * | 2015-09-06 | 2015-12-23 | 中国电力科学研究院 | Dynamic phasor measurement system of phasor measurement unit (PMU), and measurement method thereof |
| CN105242225A (en) * | 2015-09-06 | 2016-01-13 | 中国电力科学研究院 | Dynamic phasor measurement device calibration system and calibration method thereof |
| CN107561359A (en) * | 2017-08-31 | 2018-01-09 | 国网新疆电力公司 | A kind of m-Acetyl chlorophosphonazo extracting method based on Blackman window |
| CN108957129A (en) * | 2018-07-27 | 2018-12-07 | 中国南方电网有限责任公司 | A kind of new energy grid connection system time/supersynchronous oscillation self-adaptive monitoring method |
Non-Patent Citations (2)
| Title |
|---|
| 王建飞等: "一种改进型加Kaiser窗电力谐波检测方法", 《辽宁工程技术大学学报（自然科学版）》 * |
| 高云鹏等: "基于Kaiser窗双谱线插值FFT的谐波分析方法", 《仪器仪表学报》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111294018A (en) * | 2020-03-17 | 2020-06-16 | 厦门傅里叶电子有限公司 | LRA drive pulse waveform design method based on Kaiser window |
| CN111650436A (en) * | 2020-05-28 | 2020-09-11 | 电子科技大学 | Subsynchronous oscillation identification method for high-proportion renewable energy power system |
| CN111650436B (en) * | 2020-05-28 | 2021-06-04 | 电子科技大学 | Subsynchronous oscillation identification method for high-proportion renewable energy power system |
| CN112748284A (en) * | 2020-12-30 | 2021-05-04 | 清华大学 | Method and device for measuring synchronous waveform and broad-spectrum phasor of power system |
| CN117892065A (en) * | 2024-03-14 | 2024-04-16 | 江西飞尚科技有限公司 | Vibrating wire sensor signal correction method, system, computer and storage medium |
| CN117892065B (en) * | 2024-03-14 | 2024-05-31 | 江西飞尚科技有限公司 | Vibrating wire sensor signal correction method, system, computer and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109884390A (en) | A kind of subsynchronous/supersynchronous phasor measurement method of power system voltage harmonic wave | |
| Wen et al. | Hanning self-convolution window and its application to harmonic analysis | |
| CN105223418B (en) | The measurement method and measuring device of subsynchronous and supersynchronous harmonic phasor | |
| CN101587146B (en) | Parameter detection method for harmonic wave and indirect harmonic wave | |
| CN107247182B (en) | Inter-harmonic component reduction method based on measured phasor data | |
| CN110837003B (en) | Double-window full-phase DFT (discrete Fourier transform) synchronous phasor measurement method and system based on triangular window | |
| CN106841778B (en) | PMU-based subsynchronous and supersynchronous harmonic parameter processing method | |
| CN103353550A (en) | Method for measuring signal frequency and harmonic parameters of electric power system | |
| CN103575984A (en) | Harmonic analysis method based on Kaiser window double-spectral-line interpolation FFT | |
| CN110095650A (en) | The complicated harmonic detecting analysis method of four spectral line interpolation FFTs based on five Rife-Vincent (I) windows | |
| CN109507480B (en) | Inter-harmonic detection method and device for adjacent fundamental waves/harmonic waves | |
| CN103308804A (en) | Method for extracting time-frequency parameters of power quality disturbance signals on basis of fast K-S (Kaiser-S) transformation | |
| CN106771591B (en) | A kind of method for parameter estimation of Complex Power harmonic wave | |
| CN108918964A (en) | A kind of harmonic analysis method of sparsity enhancing | |
| CN108535613B (en) | Voltage flicker parameter detection method based on combined window function | |
| CN107271774A (en) | A kind of APF harmonic detecting methods based on spectrum leakage correcting algorithm | |
| CN112505407B (en) | Power grid broadband oscillation monitoring method, system and equipment and readable storage medium | |
| CN105137181A (en) | Double-spectrum-line interpolation harmonic analysis algorithm based on Nuttall-Kaiser composite window | |
| CN110940877A (en) | Voltage flicker detection method based on windowed FFT interpolation algorithm | |
| CN103091545A (en) | Sinusoidal signal phasor half-wave computing method irrelevant to frequency | |
| CN105486921A (en) | Kaiser third-order mutual convolution window triple-spectrum-line interpolation harmonic wave and inter-harmonic wave detection method | |
| CN103018555A (en) | High-precision electric power parameter software synchronous sampling method | |
| CN108896944A (en) | A kind of synchronous measuring apparatus laboratory investment instrument and its synchronous phasor measuring method | |
| CN110954746A (en) | Six-interpolation FFT algorithm based on four-term Nuttall cosine window | |
| CN108776263A (en) | Harmonic detecting method based on the peaceful self-convolution window of the high-order Chinese and improvement interpolation algorithm |
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 | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20190614 |
|
| WW01 | Invention patent application withdrawn after publication |