CN104635094A - Method for improving PMU (power management unit) synchronous phasor measurement precision - Google Patents
Method for improving PMU (power management unit) synchronous phasor measurement precision Download PDFInfo
- Publication number
- CN104635094A CN104635094A CN201510093656.7A CN201510093656A CN104635094A CN 104635094 A CN104635094 A CN 104635094A CN 201510093656 A CN201510093656 A CN 201510093656A CN 104635094 A CN104635094 A CN 104635094A
- Authority
- CN
- China
- Prior art keywords
- phasor
- frequency
- theta
- dft
- pmu
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
The invention discloses a method for improving the PMU (power management unit) synchronous phasor measurement precision under the dynamic conditions such as out-of-band interference, system low-frequency oscillation, power oscillation, system out of step, short circuit or broken circuit fault, and belongs to the technical field of electric automation. The method comprises the following steps that the high-precision synchronous sampling is carried out; DFT (discrete Fourier transform) phasor calculation and pre-filtering are carried out; the high-precision frequency measurement is carried out, and the DFT calculation phasor under the frequency deviation condition is corrected; on the basis of the phasor, the frequency and the frequency change rate worked out in the step SS3, the filtering is carried out through a rear-mounted up-supply filter, and the data is sent to a back end data concentrator; test and verification are carried out. The method has the advantages that the inhabitation on the spectrum leakage and picket fence effect during the frequency offset by a DFT algorithm is realized, the spectrum aliasing phenomenon during the calculation of components with the oscillation frequency being higher than the Nyquist frequency is avoided, and the dynamic measuring performance of the PMU under the conditions of modulation, phase step and the like is improved.
Description
Technical field
The present invention relates to a kind of under the dynamic conditions such as the outer interference of band, low frequency oscillations and oscillation of power, system step-out, short circuit or disconnection fault, promote the method for PMU synchronous phasor measurement precision, belong to technical field of electric power automation.
Background technology
Synchronous phasor measurement unit (phasor measurement unit, PMU) as WAMS (wide area measurement system, WAMS) raw data source, its measuring accuracy for Power Network Status Estimation, the important in inhibitings such as electricity net safety stable monitoring.
The definition of synchronized phasor as shown in Figure 1, simulating signal
corresponding phasor form is
when the maximal value of v (t) appears at pulse per second (PPS), the angle of phasor is 0 degree, and when v (t) positive going zeror crossing point is synchronous with pulse per second (PPS), the angle of phasor is-90 to spend.
When phasor amplitude is constant, the phase place of phasor and the frequency of simulating signal should meet following relation:
Namely, when the frequency of phasor equals 50Hz, the angle of phasor is constant; When the frequency of phasor is greater than 50Hz, the angle of phasor increases gradually, and when the frequency of phasor is less than 50Hz, the angle of phasor reduces gradually.
It is typical modulated process that electric system occurs when vibrating, and modulation frequency range is at 0.1Hz-2.5Hz.In order to realize the accurate identification of vibrating, requiring that PMU can follow the tracks of amplitude and the phase angle of phasor fast and accurately, good amplitude versus frequency characte can be had at passband simultaneously.During the fault that is short-circuited, likely there is transition in amplitude, the phase angle of voltage signal, is step process.In order to reflect the state before and after electric network fault strictly according to the facts, requiring that PMU can tracking signal transition, making error fall into accuracy rating as early as possible, and controlling the overshoot of response.
The Measurement Algorithm existence two large class problems that current PMU is general:
(1) algorithm that PMU adopts usually is at present based on DFT (discrete Fourier transform, DFT) algorithm, DFT algorithm has spectral leakage and fence effect when frequency deviation, equalization effect can affect dynamic accuracy, therefore in modulated process, the measured value of PMU can exceed accuracy rating, and its concrete reason is analyzed as follows:
The DFT formula calculating fundamental phasors is as follows:
In formula,
for fundamental phasors, x (k) is dis-crete sample values, and N is every cycle sampling number.
Sampled signal in time window is multiplied with synchronous orthogonal coefficient by being in the nature that DFT calculates, and then averaging obtains phasor; Under steady state conditions, DFT can filtering harmonic wave effectively, obtains fundamental phasors accurately.
Suppose that signal generation amplitude is vibrated, its expression formula is:
In formula, X
mfor fundamental phasors amplitude, f is fundamental frequency,
for modulating part initial phase angle.
Original sample value is multiplied by orthogonal coefficient:
Real imaginary part when all-round DFT calculates is f by constant, frequency
acomponent and high fdrequency component composition.When adopting all-round DFT to carry out phasor calculation, frequency is that the high fdrequency component of 2f can by filtering, and frequency is 2f-f
aand 2f+f
ahigh fdrequency component then can exist residual, thus result in the generation of spectral leakage.Spectral leakage makes the phasor measurement result of PMU there is deviation, affect the various analytical applications based on PMU data, making systems axiol-ogy to there is not the higher-order of oscillation, wide-area control system mistake time serious, may be made to send out steering order, totally unfavorable to the safe and stable operation of wide area system.
(2) at synchronized phasor by substation in the process of master station transmission, transfer rate Fs limits signal bandwidth, according to sampling thheorem, exceedes the interference of nyquist frequency Fs/2, then can cause spectral aliasing if existed in signal.Specifically, if electric system fundamental frequency is expressed as f
0, phasor uploading rate is expressed as F
s, then nyquist frequency is F
s/ 2, main website can accurate measurement electric system letter frequency band be [f
0-(F
s/ 2), f
0+ (F
s/ 2) frequency], beyond this frequency band range is called as the outer frequency of band.Therefore, when carrying out the transmission of different rates phasor, suppress to disturb outward with filter out-band, it is as follows that the outer interference of band calculates to synchronized phasor the analysis of causes producing error:
Suppose to there is the out-of-band signaling relevant to real-time transfer rate in input signal, expression formula is as follows:
In formula, X
mfor fundamental phasors amplitude, f is fundamental frequency,
for phasor initial phase angle, X
dfor the which amplitude modulation degree of depth, f
dthe outer frequency of band,
it is modulating part initial phase angle.
Original sample value is multiplied by orthogonal coefficient:
When input signal exists out-of-band signaling, the real imaginary part that DFT calculates is except constant, and also there is frequency is | f
d-f| and f
dtwo components of+f.When adopting all-round DFT to carry out phasor calculation, frequency is | f
d-f| and f
dtwo components of+f cannot filtering completely, easily causes alarm by mistake, and high-frequency signal residual in spectral leakage, then increase the weight of the impact of spectral aliasing further.
Summary of the invention
Order of the present invention is to provide a kind of method promoting PMU synchronous phasor measurement precision, the spectral leakage of DFT algorithm when frequency deviation and fence effect are suppressed, to avoid calculating in phasor containing the aliasing occurred during the component of oscillation frequency higher than nyquist frequency, improve the kinetic measurement performance of PMU under the conditions such as modulation, step.
The present invention adopts following technical scheme: a kind of method promoting PMU synchronous phasor measurement precision, is characterized in that, comprise the steps:
SS1 high-precise synchronization sample: PMU first by input PT/CT analog signals through anti-aliasing analog filtering (analogue low pass filtering, passband is 3.1KHz) after carry out A/D conversion, when analog to digital conversion, sampled value is stamped accurate absolute time mark, and stored in data buffer;
SS2DFT phasor calculation and pre-filtering: the sampled value of the data buffer produced by step SS1 carries out DFT conversion, then low-pass filtering is carried out to the real imaginary part of DFT calculating phasor, prefilter adopts FIR lowpass digital filter; The ripple method designs such as described FIR lowpass digital filter employing, group delay is 50 milliseconds, and passband edges frequency is 5Hz, and passband gain is 0.0002dB, and stop band gain is-80dB;
SS3 high accuracy frequency measurement DFT under frequency of amendment drift condition calculates phasor: the result of calculation based on step SS2 adopts the phasor measurement multiprecision arithmetic based on equal interval sampling, calculated the iteration of phasor by continuous 3 DFT, in 5ms, frequency, frequency change rate DFT under frequency of amendment drift condition calculates phasor can be calculated accurately;
The phasor that SS4 calculates based on step SS3, frequency and frequency change rate are again after sending filter filtering on one rearmounted, on deliver to the data concentrator of rear end, described rearmounted on send wave filter to be FIR lowpass digital filter, the ripple method designs such as described FIR lowpass digital filter employing, and select different parameters according to different phasor uploading rate;
SS5 test and validation: clock source provides common benchmark to PMU and test source, theoretical value for generating waveform playback, and compares with the simulation main website off-line files of calling and carries out error analysis; When carrying out error analysis, clip the data of each a second before and after test and carry out point-by-point comparison and get maximum error, avoid analog quantity to apply and exit steps jumping over the impact of journey; Test event comprises ± 5Hz rated frequency deviation test, and 10%2-13 subharmonic, outside band, step, 5Hz amplitude phase angle is modulated simultaneously, 1Hz/s frequency slope.
Preferably, described step SS1 specifically comprises: signal when PMU device can be accessed by the B code pair of GPS/ big dipper clock, is decoded by fpga chip, produces 1PPS signal, and producing the 4K sampled signal synchronous with 1PPS signal, this signal enabling AD sampling A/D chip carries out analog to digital conversion; PT/CT analog signals accesses A/D conversion chip after anti-aliasing analog filtering (analogue low pass filtering passband is 3.1KHz), and analog-to-digital result is put into data buffer and carried out DFT phasor calculation.
Preferably, the low-pass filtering described in described step SS2 specifically comprises: the higher-order of oscillation component that filtering spectral leakage and the outer frequency of band produce, and eliminate interference, improve phasor calculation precision, the phasor calculation formula adding prefilter is as follows:
In formula (8), (9), f
0for mains frequency; L is for there being limit for length's unit impact response digital filter exponent number; W
kfor low-pass filter coefficients, N is sampling number, and x (i+k) is i+k moment phasor value.
Preferably, described step SS3 specifically comprises: the time domain sinusoidal signal that is sampled of case of external input is
A kth sampled point in r window can be expressed as:
Phasor computing algorithm adopts recurrence DFT formula to be expressed as:
Again because
Order
Namely obtaining by formula (11), (12), (13), (14) the original phasor that DFT calculates is:
In formula:
Formula (15) gives recurrence Fu Shi algorithm r and walks the original phasor that there is spectrum leakage error calculating and obtain
phasor true in theory
between relational expression,
for the phasor that there is spectrum leakage error that Fu Shi algorithm calculates, C
rthe frequency leakage item that initial phase angle corresponding to data window is relevant,
for the true phasor of theory; Can solving equation (15) be passed through thus, (16) obtain revised free from error true phasor
if every m subsynchronous phasor of sampling calculating 1
adopt 3 equally spaced phasors
set up Simultaneous Equations:
Solve 3 simultaneous equations, obtain:
In formula:
Make Φ
m=e
-jm θ, then
In the hope of θ, then θ can be substituted into formula (14) by formula (20), the side-play amount of ongoing frequency relative system rated frequency can be tried to achieve, and then can in the hope of system frequency and frequency change rate;
be the phasor of step SS2 result of calculation, through type (21) can calculate C
r, then the substantial synchronization phasor in frequency shift (FS) situation can be calculated by formula (24)
Preferably, described step SS4 Selecting parameter is: when transfer rate is 25Hz, group delay 250 milliseconds, passband edges 5Hz, passband gain 0.0002dB, stop band gain-40dB; When transfer rate is 50Hz, group delay 125 milliseconds, passband edges 5Hz, passband gain 0.0002dB, stop band gain-60dB; When transfer rate is 100Hz, without the need to postfilter.
The beneficial effect that the present invention reaches: the spectral leakage of DFT algorithm when frequency deviation and fence effect are suppressed, to avoid calculating in phasor containing the aliasing occurred during the component of oscillation frequency higher than nyquist frequency, improve the kinetic measurement performance of PMU under the conditions such as modulation, step
Accompanying drawing explanation
Fig. 1 is the transformational relation figure between waveform signal and synchronized phasor.
Fig. 2 is a kind of process flow diagram promoting the method for PMU synchronous phasor measurement precision of the present invention.
Fig. 3 is the amplitude-versus-frequency curve figure of equiripple filter of the present invention.
Fig. 4 is PMU test macro configuration diagram of the present invention.
Embodiment
Below in conjunction with accompanying drawing, the invention will be further described.Following examples only for technical scheme of the present invention is clearly described, and can not limit the scope of the invention with this.
As shown in Figure 1 be transformational relation figure between waveform signal and synchronized phasor, Fig. 2 is a kind of process flow diagram promoting the method for PMU synchronous phasor measurement precision of the present invention, the present invention proposes a kind of method promoting PMU synchronous phasor measurement precision, it is characterized in that, comprises the steps:
SS1 high-precise synchronization sample: PMU first by input PT/CT analog signals through anti-aliasing analog filtering (analogue low pass filtering, passband is 3.1KHz) after carry out A/D conversion, when analog to digital conversion, sampled value is stamped accurate absolute time mark, and stored in data buffer; Described step SS1 specifically comprises: signal when PMU device can be accessed by the B code pair of GPS/ big dipper clock, decoded by fpga chip, produce 1PPS signal, and produce the 4K sampled signal synchronous with 1PPS signal, this signal enabling AD sampling A/D chip carries out analog to digital conversion; PT/CT analog signals accesses A/D conversion chip after anti-aliasing analog filtering (analogue low pass filtering passband is 3.1KHz), and analog-to-digital result is put into data buffer and carried out DFT phasor calculation.
SS2DFT phasor calculation and pre-filtering: the sampled value of the data buffer produced by step SS1 carries out DFT conversion, then low-pass filtering is carried out to the real imaginary part of DFT calculating phasor, prefilter adopts FIR lowpass digital filter; The ripple method designs such as described FIR lowpass digital filter employing, group delay is 50 milliseconds, and passband edges frequency is 5Hz, and passband gain is 0.0002dB, and stop band gain is-80dB; Low-pass filtering described in described step SS2 specifically comprises: the higher-order of oscillation component that filtering spectral leakage and the outer frequency of band produce, and eliminate interference, improve phasor calculation precision, the phasor calculation formula adding prefilter is as follows:
In formula (8), (9), f
0for mains frequency; L is for there being limit for length's unit impact response digital filter exponent number; W
kfor low-pass filter coefficients, N is sampling number, and x (i+k) is i+k moment phasor value.
Usually, the exponent number of wave filter is higher, and stopband attenuation is larger, and passband error is less, and transitional zone is narrower, but step response can be slower, and time delay also can be longer.Relative to the 20ms time delay of complete cycle ripple DFT phasor calculation, the group delay of wave filter wants large many, and at this moment performance of filter still may can not meet requirement.Therefore, need the wave filter of devise optimum, with the filtering performance that low order realization is as far as possible best.
FIR filter design mainly contains window function metht, frequency sampling method and method such as ripple such as grade.The shortcoming of window function metht is: be not easy the wave filter designing cutoff frequency given in advance; The filter order gone out designed by when meeting same design objective is usually bigger than normal.The shortcoming of frequency sampling method is that the value of cutoff frequency is limited.And the approximate error of window function metht and frequency sampling method is not equally distributed on frequency band interval, larger near band edge error, less away from band edge error.Be a kind of optimal-design method Deng the best approximatioss of ripple, the filter freguency response designed in this way is minimum relative to ideal filter maximum error, and the amplitude versus frequency characte of equiripple filter as shown in Figure 3.
SS3 high accuracy frequency measurement DFT under frequency of amendment drift condition calculates phasor: the result of calculation based on step SS2 adopts the phasor measurement multiprecision arithmetic based on equal interval sampling, calculated the iteration of phasor by continuous 3 DFT, in 5ms, frequency, frequency change rate DFT under frequency of amendment drift condition calculates phasor can be calculated accurately; Described step SS3 specifically comprises: the time domain sinusoidal signal that is sampled of case of external input is
A kth sampled point in r window can be expressed as:
Phasor computing algorithm adopts recurrence DFT formula to be expressed as:
Again because
Order
Namely obtaining by formula (11), (12), (13), (14) the original phasor that DFT calculates is:
In formula:
Formula (15) gives recurrence Fu Shi algorithm r and walks the original phasor that there is spectrum leakage error calculating and obtain
phasor true in theory
between relational expression,
for the phasor that there is spectrum leakage error that Fu Shi algorithm calculates, C
rthe frequency leakage item that initial phase angle corresponding to data window is relevant,
for the true phasor of theory; Can solving equation (15) be passed through thus, (16) obtain revised free from error true phasor
if every m subsynchronous phasor of sampling calculating 1
adopt 3 equally spaced phasors
set up Simultaneous Equations:
Solve 3 simultaneous equations, obtain:
In formula:
Make Φ
m=e
-jm θ, then
In the hope of θ, then θ can be substituted into formula (14) by formula (20), the side-play amount of ongoing frequency relative system rated frequency can be tried to achieve, and then can in the hope of system frequency and frequency change rate;
be the phasor of step SS2 result of calculation, through type (21) can calculate C
r, then the substantial synchronization phasor in frequency shift (FS) situation can be calculated by formula (24)
The phasor that SS4 calculates based on step SS3, frequency and frequency change rate are again after sending filter filtering on one rearmounted, on deliver to the data concentrator of rear end, described rearmounted on send wave filter to be FIR lowpass digital filter, the ripple method designs such as described FIR lowpass digital filter employing, and select different parameters according to different phasor uploading rate; Low-frequency oscillation of electric power system frequency is 0.1Hz ~ 2.5Hz, according to previous analysis result, the high fdrequency component scope that spectral leakage causes is 97.5Hz ~ 102.5Hz, digital filter should the component of effective these frequency ranges of filtering, synchronous phasor measuring device detects test specification and defines the outer frequency range of band corresponding to different uploading rate, and ginseng is shown in Table 1.
Table 1 transfer rate and the outer frequency table of comparisons of band
Uploading rate | Be with outer frequency range |
25Hz | 10Hz~37.5Hz 62.5Hz~100Hz |
50Hz | 10Hz~25Hz 75~100Hz |
100Hz | 100Hz~150Hz |
According to emulation testing, the calculating phasor of prefilter to DFT is carried out filtering and is inhibit DFT spectrum leakage and the outer interference of part band, but can not inhibition zone disturb outward completely.Therefore also need after the phasor of step SS3, frequency and frequency change rate calculate again after sending wave filter on one rearmounted again on deliver to the data concentrator of rear end, postfilter selects different parameters according to different phasor uploading rate.During transfer rate 25Hz, group delay 250 milliseconds, passband edges 5Hz, passband gain 0.0002dB, stop band gain-40dB; During transfer rate 50Hz, group delay 125 milliseconds, passband edges 5Hz, passband gain 0.0002dB, stop band gain-60dB, during transfer rate 100Hz, without the need to arranging postfilter.
SS5 test and validation: in order to the actual effect of verification algorithm, has built PMU test macro configuration diagram of the present invention as shown in Figure 4 and has surveyed, and compared with not adopting the PMU of hereinbefore method.Clock source provides common benchmark to PMU and test source, and theoretical value for generating waveform playback, and compares with the simulation main website off-line files of calling and carries out error analysis.When carrying out error analysis, clip the data of each a second before and after test and carry out point-by-point comparison and get maximum error, avoid analog quantity to apply and exit steps jumping over the impact of journey.Test event comprises ± 5Hz rated frequency deviation test, and 10%2-13 subharmonic, outside band, step, 5Hz amplitude phase angle is modulated simultaneously, and 1Hz/s frequency slope is as shown in table 2.
Table 2 dynamic property accuracy test contrasts
The above is only the preferred embodiment of the present invention; it should be pointed out that for those skilled in the art, under the prerequisite not departing from the technology of the present invention principle; can also make some improvement and distortion, these improve and distortion also should be considered as protection scope of the present invention.
Claims (5)
1. promote a method for PMU synchronous phasor measurement precision, it is characterized in that, comprise the steps:
SS1 high-precise synchronization sample: PMU first by input PT/CT analog signals through anti-aliasing analog filtering (analogue low pass filtering, passband is 3.1KHz) after carry out A/D conversion, when analog to digital conversion, sampled value is stamped accurate absolute time mark, and stored in data buffer;
SS2 DFT phasor calculation and pre-filtering: the sampled value of the data buffer produced by step SS1 carries out DFT conversion, then low-pass filtering is carried out to the real imaginary part of DFT calculating phasor, prefilter adopts FIR lowpass digital filter; The ripple method designs such as described FIR lowpass digital filter employing, group delay is 50 milliseconds, and passband edges frequency is 5Hz, and passband gain is 0.0002dB, and stop band gain is-80dB;
SS3 high accuracy frequency measurement DFT under frequency of amendment drift condition calculates phasor: the result of calculation based on step SS2 adopts the phasor measurement multiprecision arithmetic based on equal interval sampling, calculated the iteration of phasor by continuous 3 DFT, in 5ms, frequency, frequency change rate DFT under frequency of amendment drift condition calculates phasor can be calculated accurately;
The phasor that SS4 calculates based on step SS3, frequency and frequency change rate are again after sending filter filtering on one rearmounted, on deliver to the data concentrator of rear end, described rearmounted on send wave filter to be FIR lowpass digital filter, the ripple method designs such as described FIR lowpass digital filter employing, and select different parameters according to different phasor uploading rate;
SS5 test and validation: clock source provides common benchmark to PMU and test source, theoretical value for generating waveform playback, and compares with the simulation main website off-line files of calling and carries out error analysis; When carrying out error analysis, clip the data of each a second before and after test and carry out point-by-point comparison and get maximum error, avoid analog quantity to apply and exit steps jumping over the impact of journey; Test event comprises ± 5Hz rated frequency deviation test, and 10%2-13 subharmonic, outside band, step, 5Hz amplitude phase angle is modulated simultaneously, 1Hz/s frequency slope.
2. a kind of method promoting PMU synchronous phasor measurement precision according to claim 1, it is characterized in that, described step SS1 specifically comprises: signal when PMU device can be accessed by the B code pair of GPS/ big dipper clock, decoded by fpga chip, produce 1PPS signal, and producing the 4K sampled signal synchronous with 1PPS signal, this signal enabling AD sampling A/D chip carries out analog to digital conversion; PT/CT analog signals accesses A/D conversion chip after anti-aliasing analog filtering (analogue low pass filtering passband is 3.1KHz), and analog-to-digital result is put into data buffer and carried out DFT phasor calculation.
3. a kind of method promoting PMU synchronous phasor measurement precision according to claim 1, it is characterized in that, low-pass filtering described in described step SS2 specifically comprises: the higher-order of oscillation component that filtering spectral leakage and the outer frequency of band produce, eliminate interference, improve phasor calculation precision, the phasor calculation formula adding prefilter is as follows:
4. a kind of method promoting PMU synchronous phasor measurement precision according to claim 1, it is characterized in that, described step SS3 specifically comprises: the time domain sinusoidal signal that is sampled of case of external input is
A kth sampled point in r window can be expressed as:
Phasor computing algorithm adopts recurrence DFT formula to be expressed as:
Again because
Order
Namely obtaining by formula (11), (12), (13), (14) the original phasor that DFT calculates is:
In formula:
Solve 3 simultaneous equations, obtain:
In formula:
Make Φ
m=e
-jm θ, then
In the hope of θ, then θ can be substituted into formula (14) by formula (20), the side-play amount of ongoing frequency relative system rated frequency can be tried to achieve, and then can in the hope of system frequency and frequency change rate;
be the phasor of step SS2 result of calculation, through type (21) can calculate C
r, then the substantial synchronization phasor in frequency shift (FS) situation can be calculated by formula (24)
5. a kind of method promoting PMU synchronous phasor measurement precision according to claim 1, is characterized in that, described step SS4 Selecting parameter is: when transfer rate is 25Hz, group delay 250 milliseconds, passband edges 5Hz, passband gain 0.0002dB, stop band gain-40dB; When transfer rate is 50Hz, group delay 125 milliseconds, passband edges 5Hz, passband gain 0.0002dB, stop band gain-60dB; When transfer rate is 100Hz, without the need to postfilter.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510093656.7A CN104635094A (en) | 2015-03-02 | 2015-03-02 | Method for improving PMU (power management unit) synchronous phasor measurement precision |
PCT/CN2015/090850 WO2016138764A1 (en) | 2015-03-02 | 2015-09-26 | Method for improving pmu synchronous phasor measurement precision |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510093656.7A CN104635094A (en) | 2015-03-02 | 2015-03-02 | Method for improving PMU (power management unit) synchronous phasor measurement precision |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104635094A true CN104635094A (en) | 2015-05-20 |
Family
ID=53214061
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510093656.7A Pending CN104635094A (en) | 2015-03-02 | 2015-03-02 | Method for improving PMU (power management unit) synchronous phasor measurement precision |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN104635094A (en) |
WO (1) | WO2016138764A1 (en) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105224811A (en) * | 2015-10-21 | 2016-01-06 | 中国科学院光电技术研究所 | PMU dynamic data processing method based on feedback iterative frequency tracking |
WO2016138764A1 (en) * | 2015-03-02 | 2016-09-09 | 国电南瑞科技股份有限公司 | Method for improving pmu synchronous phasor measurement precision |
CN107132500A (en) * | 2017-03-14 | 2017-09-05 | 国家电网公司 | A kind of synchronous phasor measurement unit on-line calibration method and apparatus |
CN108614155A (en) * | 2018-05-31 | 2018-10-02 | 许继集团有限公司 | A kind of synchronous phasor measuring method and system that Hamming window is added |
CN108896944A (en) * | 2018-03-16 | 2018-11-27 | 华北电力大学 | A kind of synchronous measuring apparatus laboratory investment instrument and its synchronous phasor measuring method |
CN109143137A (en) * | 2018-06-14 | 2019-01-04 | 国电南瑞科技股份有限公司 | A kind of generator PMU device measurement accuracy method for improving |
CN109644121A (en) * | 2016-12-23 | 2019-04-16 | 华为技术有限公司 | A kind of method and apparatus that clock is synchronous |
CN110412415A (en) * | 2019-06-14 | 2019-11-05 | 国网辽宁省电力有限公司 | It is a kind of based on DFT, it is multistage filtering and mutation criteria synchronized phasor calculation method |
CN110568309A (en) * | 2019-08-08 | 2019-12-13 | 中国农业大学 | Filter, synchronous phasor measurement system and method |
CN111505375A (en) * | 2020-04-17 | 2020-08-07 | 许继集团有限公司 | Frequency measurement method and device |
CN112162176A (en) * | 2020-09-30 | 2021-01-01 | 国网河南省电力公司洛阳供电公司 | Power distribution network interphase short circuit fault positioning method based on mu PMU measurement data |
CN112255457A (en) * | 2020-09-22 | 2021-01-22 | 天津电气科学研究院有限公司 | Phase angle difference measuring method suitable for automatic quasi-synchronization device |
CN112485524A (en) * | 2020-11-10 | 2021-03-12 | 广东电网有限责任公司广州供电局 | Phasor calculation method of high-precision calibrator for PMU (phasor measurement Unit) test |
CN112946337A (en) * | 2021-05-13 | 2021-06-11 | 深圳市鼎阳科技股份有限公司 | High-precision signal measuring method and digital oscilloscope |
CN113064021A (en) * | 2021-02-04 | 2021-07-02 | 国电南瑞科技股份有限公司 | Measurement and control device and method for realizing high-order harmonic suppression of power electronic grid |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1477401A (en) * | 2003-07-18 | 2004-02-25 | 清华大学 | High-accuracy synchronous phasor measuring method |
CN101561462A (en) * | 2009-04-27 | 2009-10-21 | 国网电力科学研究院 | Phasor amplitude measurement method of power system |
WO2011089421A1 (en) * | 2010-01-22 | 2011-07-28 | The University Of Liverpool | Apparatus and method for measuring a phasor of an electrical power system signal |
US20130193951A1 (en) * | 2012-01-27 | 2013-08-01 | General Electric Company | Convolution Integral for Synchronized Phasor |
CN103543335A (en) * | 2013-10-30 | 2014-01-29 | 国家电网公司 | Method for measuring synchronous phasor |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104635094A (en) * | 2015-03-02 | 2015-05-20 | 国电南瑞科技股份有限公司 | Method for improving PMU (power management unit) synchronous phasor measurement precision |
-
2015
- 2015-03-02 CN CN201510093656.7A patent/CN104635094A/en active Pending
- 2015-09-26 WO PCT/CN2015/090850 patent/WO2016138764A1/en active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1477401A (en) * | 2003-07-18 | 2004-02-25 | 清华大学 | High-accuracy synchronous phasor measuring method |
CN101561462A (en) * | 2009-04-27 | 2009-10-21 | 国网电力科学研究院 | Phasor amplitude measurement method of power system |
WO2011089421A1 (en) * | 2010-01-22 | 2011-07-28 | The University Of Liverpool | Apparatus and method for measuring a phasor of an electrical power system signal |
US20130193951A1 (en) * | 2012-01-27 | 2013-08-01 | General Electric Company | Convolution Integral for Synchronized Phasor |
CN103543335A (en) * | 2013-10-30 | 2014-01-29 | 国家电网公司 | Method for measuring synchronous phasor |
Non-Patent Citations (2)
Title |
---|
张超: "同步相量测量单元PMU的研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
杨卫华 等: "基于DSP的同步相量测量装置的研究", 《单片机与嵌入式系统与应用》 * |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016138764A1 (en) * | 2015-03-02 | 2016-09-09 | 国电南瑞科技股份有限公司 | Method for improving pmu synchronous phasor measurement precision |
CN105224811A (en) * | 2015-10-21 | 2016-01-06 | 中国科学院光电技术研究所 | PMU dynamic data processing method based on feedback iterative frequency tracking |
CN105224811B (en) * | 2015-10-21 | 2018-01-26 | 中国科学院光电技术研究所 | PMU dynamic data processing method based on feedback iterative frequency tracking |
CN109644121A (en) * | 2016-12-23 | 2019-04-16 | 华为技术有限公司 | A kind of method and apparatus that clock is synchronous |
CN107132500B (en) * | 2017-03-14 | 2019-10-15 | 国家电网公司 | A kind of synchronous phasor measurement unit on-line calibration method and apparatus |
CN107132500A (en) * | 2017-03-14 | 2017-09-05 | 国家电网公司 | A kind of synchronous phasor measurement unit on-line calibration method and apparatus |
CN108896944A (en) * | 2018-03-16 | 2018-11-27 | 华北电力大学 | A kind of synchronous measuring apparatus laboratory investment instrument and its synchronous phasor measuring method |
CN108896944B (en) * | 2018-03-16 | 2020-09-29 | 华北电力大学 | Laboratory calibrator of synchronous measuring device and synchronous phasor measuring method thereof |
CN108614155A (en) * | 2018-05-31 | 2018-10-02 | 许继集团有限公司 | A kind of synchronous phasor measuring method and system that Hamming window is added |
WO2019228540A1 (en) * | 2018-05-31 | 2019-12-05 | 许继集团有限公司 | Synchronized phasor measurement method and system having added hamming window |
CN109143137A (en) * | 2018-06-14 | 2019-01-04 | 国电南瑞科技股份有限公司 | A kind of generator PMU device measurement accuracy method for improving |
CN109143137B (en) * | 2018-06-14 | 2020-09-29 | 国电南瑞科技股份有限公司 | Method for improving measurement precision of power generator PMU device |
CN110412415B (en) * | 2019-06-14 | 2021-11-23 | 国网辽宁省电力有限公司 | Synchronous phasor calculation method based on DFT, multi-order filtering and mutation criterion |
CN110412415A (en) * | 2019-06-14 | 2019-11-05 | 国网辽宁省电力有限公司 | It is a kind of based on DFT, it is multistage filtering and mutation criteria synchronized phasor calculation method |
CN110568309A (en) * | 2019-08-08 | 2019-12-13 | 中国农业大学 | Filter, synchronous phasor measurement system and method |
CN110568309B (en) * | 2019-08-08 | 2020-11-06 | 中国农业大学 | Filter, synchronous phasor measurement system and method |
CN111505375A (en) * | 2020-04-17 | 2020-08-07 | 许继集团有限公司 | Frequency measurement method and device |
CN112255457A (en) * | 2020-09-22 | 2021-01-22 | 天津电气科学研究院有限公司 | Phase angle difference measuring method suitable for automatic quasi-synchronization device |
CN112255457B (en) * | 2020-09-22 | 2022-06-07 | 天津电气科学研究院有限公司 | Phase angle difference measuring method suitable for automatic quasi-synchronization device |
CN112162176A (en) * | 2020-09-30 | 2021-01-01 | 国网河南省电力公司洛阳供电公司 | Power distribution network interphase short circuit fault positioning method based on mu PMU measurement data |
CN112162176B (en) * | 2020-09-30 | 2022-04-22 | 国网河南省电力公司洛阳供电公司 | Power distribution network interphase short circuit fault positioning method based on mu PMU measurement data |
CN112485524A (en) * | 2020-11-10 | 2021-03-12 | 广东电网有限责任公司广州供电局 | Phasor calculation method of high-precision calibrator for PMU (phasor measurement Unit) test |
CN112485524B (en) * | 2020-11-10 | 2024-01-19 | 广东电网有限责任公司广州供电局 | High-precision calibrator phasor calculation method for PMU test |
CN113064021A (en) * | 2021-02-04 | 2021-07-02 | 国电南瑞科技股份有限公司 | Measurement and control device and method for realizing high-order harmonic suppression of power electronic grid |
CN113064021B (en) * | 2021-02-04 | 2023-08-25 | 国电南瑞科技股份有限公司 | Measurement and control device and method for realizing power electronic power grid higher harmonic suppression |
CN112946337A (en) * | 2021-05-13 | 2021-06-11 | 深圳市鼎阳科技股份有限公司 | High-precision signal measuring method and digital oscilloscope |
CN112946337B (en) * | 2021-05-13 | 2021-08-03 | 深圳市鼎阳科技股份有限公司 | High-precision signal measuring method and digital oscilloscope |
Also Published As
Publication number | Publication date |
---|---|
WO2016138764A1 (en) | 2016-09-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104635094A (en) | Method for improving PMU (power management unit) synchronous phasor measurement precision | |
CN103454497B (en) | Based on the method for measuring phase difference improving windowed DFT | |
CN103487652B (en) | Frequency self-adaptive real-time fractional harmonic wave detection method | |
CN103543335B (en) | A kind of synchronous phasor measuring method | |
CN106199183B (en) | A kind of PMU and method realizing sub-synchronous oscillation on-line identification alarm | |
CN104049144A (en) | Synchronous phasor measurement implementing method with filtered-out attenuation direct current components | |
CN103869162A (en) | Dynamic signal phasor measurement method based on time domain quasi-synchronization | |
Xie et al. | Improved synchrophasor measurement to capture sub/super‐synchronous dynamics in power systems with renewable generation | |
CN102081132A (en) | Two-end distance measuring method of transmission line fault under dynamic condition | |
CN109490630B (en) | Dynamic phasor measurement method based on matrix beam | |
CN106841778A (en) | The processing method of the subsynchronous and supersynchronous harmonic parameters realized based on PMU | |
CN105137175A (en) | Windowing short-time Fourier transform three-point interpolation dynamic frequency measurement method | |
CN103983847B (en) | Based on the adaptive frequency track measuring method of RLS in a kind of synchronous phasor measurement | |
CN104020352A (en) | Synchronous phasor measuring method suitable for M-class PMU unit | |
CN105487034A (en) | 0.05-level electronic transformer verification method and system | |
CN104242267A (en) | Distance protection method for wind power generation transmission line | |
CN102033177B (en) | Method and system for measuring power angle of electric power circuit of distribution network | |
CN102253282A (en) | Method for obtaining continuous frequency spectrum interpolation power harmonic parameter of Nuttall window function | |
CN103472333A (en) | Wind power integration electric energy quality overall performance detection method | |
CN108414838B (en) | Method for measuring line impedance of inverter parallel system | |
CN105261202A (en) | Transmission delay measurement method based on B spline interpolation wavelet singularity detection | |
Ai et al. | A single monitor method for voltage sag source location using Hilbert Huang transform | |
CN102508022A (en) | Method for detecting power grid frequency by using optimal multiplier Newton algorithm | |
Xie et al. | Measurement of sub-and supersynchronous phasors in power systems with high penetration of renewables | |
CN203204090U (en) | Substation zinc oxide arrester on-line monitoring device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150520 |