CN104076203A - Ultralow-frequency inter-harmonic detection method with negative frequency influence considered - Google Patents
Ultralow-frequency inter-harmonic detection method with negative frequency influence considered Download PDFInfo
- Publication number
- CN104076203A CN104076203A CN201410317995.4A CN201410317995A CN104076203A CN 104076203 A CN104076203 A CN 104076203A CN 201410317995 A CN201410317995 A CN 201410317995A CN 104076203 A CN104076203 A CN 104076203A
- Authority
- CN
- China
- Prior art keywords
- frequency
- harmonic
- formula
- ultralow
- harmonic wave
- 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.)
- Granted
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 34
- 230000003595 spectral effect Effects 0.000 claims abstract description 52
- 238000001228 spectrum Methods 0.000 claims abstract description 18
- 238000010183 spectrum analysis Methods 0.000 claims abstract description 15
- 238000000034 method Methods 0.000 claims abstract description 9
- 230000021615 conjugation Effects 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000012886 linear function Methods 0.000 claims description 3
- 238000006467 substitution reaction Methods 0.000 claims description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 8
- 238000012937 correction Methods 0.000 abstract description 4
- 238000012360 testing method Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 108010076504 Protein Sorting Signals Proteins 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
The invention provides an ultralow- frequency inter-harmonic detection method with a negative frequency influence considered. An algorithm is different from a traditional interpolation correction algorithm and considers the interference of negative frequency on ultralow- frequency inter-harmonic spectral lines, the three spectral lines are simultaneous to form a ternary equation set to solve ultralow-frequency inter-harmonic frequency, and the amplitude and phase of inter-harmonics are solved according to the frequency spectrum value and conjugate value of the maximum spectral line. Relevant feature parameters of the ultralow- frequency inter-harmonics can be extracted fast and accurately in a short time, and the requirements for real-time performance and validity of the algorithm in the actual detection process are guaranteed. Meanwhile, a parameter correction formula generated when the inter-harmonic frequency is smaller than one frequency resolution delta f is derived, the inter-harmonic frequency which can be solved actually is lowered below 03delta f, the situation that that a traditional spectral analysis algorithm based on DFT is limited by the frequency resolution is greatly broken through, the length of an analysis signal is shortened, and other algorithms cannot break through the situation and cannot shorten the length of the analysis signals. The algorithm cannot be influenced by direct-current components, is not sensitive to system noise and guarantees accuracy and robustness.
Description
Technical field
The present invention relates to a kind of harmonic detecting method, particularly a kind of ultralow frequency harmonic detection method of considering negative frequency impact.
Background technology
Detect for harmonic wave between electric system medium and low frequency (lower than 50Hz subharmonic), the scheme that can realize at present has harmonic wave group and subgroup algorithm between window function and interpolation algorithm based on DFT (or FFT), IEC, based on Prony and expansion algorithm etc. thereof.And in electric system between harmonic wave there is the feature such as undulatory property, randomness, simultaneously IEC61000-4-7 and GB GB/T 24337-2009 for pair between the harmonic wave stability and the real-time that detect consider, regulation is used based on adding rectangular window DFT frequency spectrum analysis method and is carried out the measurement of a harmonic wave, and to limit Measuring Time window be 10 fundamental frequency cycles (being 200ms), the spectral line on frequency spectrum is spaced apart 5Hz.
Under IEC framework, for harmonic wave between the ultralow frequency that may occur in electric system (refering in particular to the subharmonic lower than 10Hz herein), the current detection algorithm based on DFT spectrum analysis lost efficacy substantially.When related algorithm based on DFT spectrum analysis detects between ultralow frequency harmonic wave under IEC framework, the analysis frequency resolution that is spectral line is 5Hz, adjacent spectral line is subject to the impact of spectrum leakage and fence effect serious, simultaneously because frequency is too low, can be for the spectral line number of analyzing limited and can not ignore the impact of negative frequency, the existence of negative frequency has further increased the weight of the interference to harmonic wave between ultralow frequency, makes to analyze to lose efficacy.For obtaining abundant spectral line information, tradition must increase sampling analysis based on DFT spectral analysis algorithm counts, but also to have increased sample window long simultaneously, thereby has affected the real-time of algorithm.
Prony and expansion algorithm thereof (expansion Prony method) are the methods in estimating based on modern spectrum, in the time detecting low frequency and ultra-low frequency component, be not subject in theory the restriction of frequency resolution, but must carry out large-scale matrix operation, also need signal to carry out determining accurately rank, cannot direct-detection for the signal of electrical network broadband multifrequency, affect real-time and the stability of algorithm.Meanwhile, the related algorithm based on Prony is very sensitive to system noise, and in actual detection, error is very large.
Can find out that prior art also exists a lot of deficiencies to the detection of harmonic wave between ultralow frequency.And correlative study also shows, occur when low-frequency oscillation and all may occur harmonic wave between ultralow frequency in emerging electric power energy (as wind energy, sun power etc.) system, can therefore accurately detect harmonic wave between ultralow frequency has practical significance to researchs such as follow-up detection, improvement and system stability.
Summary of the invention
The object of this invention is to provide a kind of ultralow frequency harmonic detection method of considering negative frequency impact, it is characterized in that, comprise the following steps:
Step 1: under operation of power networks normal condition, set up the unifrequency discrete signal model of power network signal:
(1) wherein
for the amplitude of power network signal,
for the frequency of power network signal and
for the initial phase of power network signal,
,
for the total sampling number of discrete signal intercepting,
for power network signal sample frequency;
Step 2: right
the normalization discrete Fourier transform (DFT) (Discrete Fourier Transform, DFT) that adds rectangular window obtains
(2), wherein k=0,1 ... N-1;
Step 3: the frequency resolution of establishing spectral line is
,
, and will in formula (1) substitution formula (2), obtain
X（k）
（3），
If normalized frequency is
,
, formula (3) is further reduced to:
(4), in formula
;
Step 3: get adjacent three spectral lines
,
with
spectrum value
,
with
, here
for integer and to meet its corresponding spectral line be maximum spectral line, here
for integer and to meet its corresponding spectral line be maximum spectral line, obtain through calculating abbreviation
(5)
Note
,
, have according to formula (4):
(6), wherein
=1,0,1, three unknown numbers in formula are
,
with
, when
=1,0,1, can become a ternary linear function group by simultaneous, therefrom can obtain
,
with
, wherein solve frequency normalization to be asked
expression formula be
(7)
Step 4: obtain the highest spectral line by formula (4)
equation:
(8),
In formula
,
,
, formula (9) is got to conjugation and obtains:
(9) obtain
after, simultaneous formula (8) and formula (9) composition system of equations, this system of equations the inside only has
with
be unknown quantity, all the other are known quantity, therefore can solve:
(10)
Step 5: formula (7) and formula (10) are considers humorous wave frequency, amplitude and phase estimation formula between negative frequency and interpolation ultralow.
Detection in full frequency band time harmonic frequency comprises the following steps:
1) signal electrical network being collected, intercepts one section with rectangular window and makes FFT(DFT) spectrum analysis;
2) for harmonic component between more than 50Hz, the regulation detecting according to harmonic wave between IEC61000-4-7 standard pair, with the form of a harmonic wave group or subgroup, measure frequency range between two harmonic waves and between first-harmonic and second harmonic between harmonic parameters, and frequency using the mean value of two harmonic wave sequence numbers as corresponding harmonic wave group or subgroup;
3) must consider the impact of negative frequency for harmonic wave between 0 ~ 50Hz frequency range; Between the ultralow frequency of 0 ~ 10Hz frequency range, harmonic wave calculates harmonic parameters amplitude between ultralow frequency according to ultralow frequency simple harmonic wave detection method
, frequency
and phase place
, and 10 ~ 50Hz frequency range between harmonic wave calculate all spectral lines of residue except first-harmonic spectral line according to the method for IEC, remove the root-mean-square value of harmonic amplitude after the intersegmental harmonic wave leakage contribution of ultralow frequency
;
4) in step 3)
, be 0 ~ 50Hz subharmonic frequency range between harmonic parameters, according to the regulation of IEC standard, be designated as harmonic wave between 0.5 time;
5) comprehensive step 2) and step 4) obtain full frequency band between harmonic parameters.
In an embodiment of the present invention, occur what harmonic wave between ultralow frequency produced
spectral line drops on frequency just
upper
,
spectral line drops on direct current spectral line, carries out every straight processing filtering DC component before carrying out FFT spectrum analysis
,
, make formula (6) can not be write as the secondly form of system of equations of ternary, now can select
be 1, corresponding
, wherein " * " represent get conjugation, at this moment, the spectral line sequence number of utilization be (1,1,2), correspondence spectrum value be (
,
,
), formula
can be reduced to:
(11).
The present invention is different from traditional interpolation correction algorithm, this algorithm has been considered the interference of negative frequency to Harmonic lines between ultralow frequency, become equation with three unknowns group to solve harmonic frequency between ultralow frequency by three spectral line value simultaneous, and solve amplitude and the phase place of a harmonic wave according to the spectrum value of maximum spectral line and conjugate thereof; Realize the relevant feature parameters of extracting fast and more exactly harmonic wave between ultralow frequency in short time, real-time and the validity requirement of algorithm while having ensured actual detection.Between the present invention has simultaneously derived and worked as, harmonic frequency is less than 1 frequency resolution
time parameter correction formula, actual can solve between harmonic frequency be low to moderate
below, greatly broken through traditional restriction that is subject to frequency resolution based on DFT harmonic analysis, shortened the length of analytic signal, this is that other algorithm institute is irrealizable.And algorithm of the present invention is not subject to the impact of DC component and to insensitive its accuracy and the robustness of having ensured of system noise.
Brief description of the drawings
The interference phenomenon of Fig. 1 negative frequency to frequency spectrum.
Harmonic wave overhaul flow chart between the full frequency band of Fig. 2 based on IEC and ultralow frequency detection algorithm.
The amplitude error that between Fig. 3 IEC61000-4-7 subharmonic, harmonic wave group detects.
The amplitude detection result of 10 cycle code book motion algorithms of Fig. 4 IEC regulation.
Harmonic frequency variation diagram between the minimum that the each algorithm of Fig. 5 ± 5% range of tolerable variance can detect.
Fig. 6 Prony algorithm (left side) and herein algorithm (right side) accuracy of detection changing trend diagram under noise circumstance.
Embodiment
Voltage and the frequency of considering electrical network in practical power systems are always changing, and the detection of a harmonic wave is always closely related with frequency with the voltage of electrical network, according to the voltage of electrical network and frequency change time speed, signal model in actual electric network is classified: stable state while 1) length, the voltage magnitude of power network signal and constant stationary random signal of frequency expectation value long period; 2) stable in short-term, the voltage magnitude of power network signal and constant stationary random signal of frequency expectation value short period; 3) unsteady state, the non-stationary signal becoming when the voltage magnitude of power network signal and frequency are when fault.And the foundation of model of the present invention is aimed at and analyzes under system normal operating condition, the constant expectation value of voltage magnitude and frequency voltage magnitude and frequency in actual electric network within a period of time, under a period of time steady state (SS)
,
; Consider that the noise signal in actual electric network is the white Gaussian noise signal separate with each frequency component of electrical network simultaneously
.
Set up the power network signal mathematical model containing harmonic wave between ultralow frequency:
(12)
In formula
,
with
for amplitude, frequency and the initial phase of each component in power network signal,
for the direct-flow offset weight that may comprise in actual electric network signal,
for the white Gaussian noise signal separate with each frequency component of electrical network, M is the frequency component number comprising in power network signal; Time
, the intercepting time period that T is signal analysis.
Consider that the signal collecting in actual electric network is discrete series
, therefore change the mathematical model in formula (12) into discrete form:
(13)
In formula
for the sample frequency of power network signal,
, N is the total sampling number of discrete signal intercepting.
It is pointed out that when the sampled signal of actual electric network is carried out to spectrum analysis, generally can first carry out carrying out filtering DC component every straight processing
, object be improve computing machine integer FFT do spectrum precision, prevent Frequency and interfering frequency that windowing (except rectangular window) produces.While carrying out the mathematical analysis of model, first ignore system noise part
, after analyze again the situation that affects on testing result when the system noise; Notice, the linear superposition that the signal in formula (12) is each component, carries out DFT(FFT to it simultaneously) still meet this overlaying relation when spectrum analysis.Therefore, for without loss of generality, analysis list frequency-distributed signal model:
(14)
Can find out from formula (14), the actual signal gathering from electrical network is discrete real sequence
,
, not single-frequency complex exponential signal sequence, the normalization discrete Fourier transform (DFT) (Discrete Fourier Transform, DFT) that it is added to rectangular window obtains
(15)
In formula
=0,1,, N 1.
If the frequency resolution of spectral line is
,
, by formula (14) substitution formula (15):
(16)
In actual conditions, general
be worth larger, therefore
.If normalized frequency is
,
, formula (15) is further reduced to:
(17)
In formula
.
From formula (17), can find out, actual electric network sampled signal is carried out after DFT spectrum analysis, there will be two frequency components on frequency spectrum
with
, wherein Section 2 is the negative frequency components composition of frequency spectrum.In the algorithm model of tradition based on DFT spectrum analysis, ignore the impact of negative frequency on spectral line, directly take Section 1 in formula (17), do so certain reason, reason be between reality harmonic component when enough words far away of distance 0 frequency spectrum, negative frequency is very little on the impact of harmonic component between reality, but for harmonic wave between ultralow frequency, negative frequency is very large on its impact, this is also that general algorithm is detecting the reason losing efficacy for ultralow harmonic wave, when detecting harmonic wave between power network signal ultralow frequency, reality must consider the impact of negative frequency on it.Negative frequency affects schematic diagram as shown in Figure 1 to frequency spectrum.
Consider the discrete spectrum of negative frequency impact, get adjacent three spectral lines
,
with
spectrum value
,
with
, here
for integer and to meet its corresponding spectral line be maximum spectral line.So
(18)
In like manner can verify
(19)
Therefore, note
,
, have according to formula (15):
(20)
In formula, three unknown numbers are
,
with
, wherein
=1,0,1, can become a ternary linear function group by simultaneous, and there is unique solution in this equation, therefrom can obtain
,
with
.For the convenience solving, formula (20) is considered as to system of homogeneous linear equations, untrivialo solution vector be
,
, 1}
t, the unique conditional of its existence is the matrix of coefficients of system of equations
determinant be zero
.
Matrix of coefficients
for:
(21)
By
can solve frequency normalization to be asked
explicit Expression formula be
(22)
For the accuracy that ensures to calculate, need in guarantee formula (22)
spectral line is the highest spectral line,
with
be respectively the adjacent spectral line in its left and right, one of them must be time high spectral line so, so the frequency estimation of harmonic wave is between the ultralow frequency of consideration negative frequency impact
(23)
In formula
for frequency resolution.
The impact of taking into account system DC component: probably occur what harmonic wave between ultralow frequency produced in reality
spectral line drops on frequency just
upper
, like this
spectral line drops on direct current spectral line.Owing to originally may there being the interference of DC component in system, before carrying out FFT spectrum analysis, carry out every straight processing filtering DC component
,
, make formula (20) can not be write as the secondly form of system of equations of ternary, now can select
be 1, corresponding
, " * " represents to get conjugation.At this moment, the spectral line sequence number of utilization is (1,1,2), corresponding spectrum value be (
,
,
), formula (22) can abbreviation be:
(24)
In order to ensure signal to noise ratio (S/N ratio) and the computational solution precision of actual signal spectral line, carry out harmonic amplitude and phase estimation between ultralow frequency time and other interpolation algorithm similar, need to select the highest spectral line in part to calculate.By the highest spectral line of formula (17)
equation:
(25)
In formula
,
。
Wherein
.
Formula (25) is got to conjugation to be obtained:
(26)
Obtained by formula (22)
after, simultaneous formula (25) and formula (26) composition system of equations, this system of equations the inside only has
with
be unknown quantity, all the other are known quantity, therefore can solve:
(27)
Formula (22) and formula (27) are considers humorous wave frequency, amplitude and phase estimation formula between negative frequency and interpolation ultralow.It should be noted that above-mentioned model is based upon signal and intercepts by rectangular window (i.e. not windowing), this is also the testing requirement that meets harmonic wave between IEC61000-4-7 pair.
Harmonic wave between the ultralow frequency having detected for above-mentioned theory, and the regulation detecting in conjunction with harmonic wave between 24337-2009 couple of IEC61000-4-7 and GB GB/T, can form a whole set of be applicable to Whole frequency band between harmonic wave detection scheme, its FB(flow block) is as shown in Figure 2.Its detailed step is:
1) signal electrical network being collected, intercepts one section (as 10 cycles of IEC regulation, 200ms, respective frequencies resolution is 5Hz) with rectangular window and makes FFT(DFT) spectrum analysis.
2) for harmonic component between more than 50Hz, the regulation detecting according to harmonic wave between IEC61000-4-7 pair, with the form of a harmonic wave group or subgroup, measure frequency range between two harmonic waves and between first-harmonic and second harmonic between harmonic parameters, and frequency using the mean value of two harmonic wave sequence numbers as corresponding harmonic wave group or subgroup (as harmonic wave between between 2 subharmonic and 3 subharmonic be designated as 2.5 between harmonic wave group or subgroup).
3) for the necessary impact of considering negative frequency of harmonic wave between 0 ~ 50Hz frequency range (being subharmonic), the frequency of subharmonic is less, and the impact that is subject to negative frequency is more serious; Consider same frequency range between harmonic wave quantity can be not too many, there is the fluctuation of frequency and amplitude in electrical network simultaneously, the 5Hz frequency resolution specifying taking IEC is example, and between the ultralow frequency of 0 ~ 10Hz frequency range, harmonic wave calculates harmonic parameters (amplitude between ultralow frequency according to algorithm mentioned in this article
, frequency
and phase place
), and 10 ~ 50Hz frequency range between harmonic wave calculate all spectral lines of residue (containing first-harmonic spectral line) according to the method for IEC and remove the root-mean-square value after the intersegmental harmonic wave leakage contribution of ultralow frequency
.
4) the 3rd) in step
, be 0 ~ 50Hz subharmonic frequency range between harmonic parameters, according to the regulation of IEC standard, be designated as harmonic wave between 0.5 time.
5) comprehensive step 2) and step 4) obtain full frequency band between harmonic parameters.
The harmonic detection method proposing for the present invention carries out Digital Simulation, simulated environment is Matlab 2012b, the 0.2s sampled signal length (synchronized sampling) specifying with IEC61000-4-7 in emulation embodiment, corresponding 5Hz frequency resolution, if humorous wave frequency is 3.5Hz between ultralow frequency, amplitude coefficient with respect to fundamental frequency is 1%, and under noise circumstance, the testing result of various algorithms is as shown in table 1:
The accuracy of detection of each algorithm in table 1 noise circumstance
Note: while having noise, carry out 100 Monte Carlo simulations and test the mean value of asking for result; "---" represents that corresponding algorithm is without correct testing result.
As can be seen from Table 1, the window function and interpolation algorithm based on DFT lost efficacy in the time detecting between ultralow frequency harmonic wave, and reason is the restriction that is subject to frequency resolution; Prony algorithm has high accuracy of detection in the time of noiseless, but trickle noise can cause sharply dropping to of precision to detect, and its robustness is bad; The method of carrying herein can have higher accuracy of detection in larger noise circumstance, has ensured robustness and the real-time of algorithm in real system detects.
The comparison of the present invention and other algorithms:
(1) algorithm and the comparison of IEC canonical algorithm herein
Between IEC61000-4-7 regulation, harmonic wave group algorithm detects subharmonic, supposes that between low frequency, the relative amplitude of harmonic wave is
, frequency is
, phase place is 0.Sampling time is 0.2s, when
while variation from 0 ~ 40Hz with the step-length of 0.05Hz, corresponding testing result as shown in Figure 3
As seen from Figure 3, IEC algorithm is 5Hz in frequency resolution, range of allowable error is in ± 5% time, between the minimum that can detect, harmonic frequency is about 7.5Hz, consider the phase angle variations of harmonic wave and the stability of algorithm between ultralow frequency, between 0 ~ 10Hz frequency range, harmonic wave IEC algorithm cannot detect exactly simultaneously.
And herein algorithm shows as Fig. 4 the intersegmental harmonic wave testing result of ultralow frequency:
As seen from Figure 4, for the sampled signal of 0.2s length, doing frequency resolution corresponding to spectrum analysis is 5Hz, algorithm herein can IEC standard regulation ± detect in 5% error range low tone harmonic frequency be 1Hz(now between correspondence harmonic phase be 0), greatly shake off the restriction of frequency resolution, realized and in the short period, detected rapidly and accurately harmonic wave between ultralow frequency.
Under noiseless disturbed condition, in specify ± 5% range of tolerable variance of IEC, between the phase angle variations scope of harmonic wave be 0 ~ 2 π, obtain harmonic frequency between minimum that herein algorithm and IEC canonical algorithm can detect in the time of unlike signal length, the numerical value of enumerating and changing trend diagram thereof are respectively as shown in table 2 and Fig. 5.
Harmonic frequency between the minimum that in table 2 range of tolerable variance, algorithm can be identified
Can be found out by table 2 and Fig. 5, herein algorithm can detect harmonic wave between ultralow frequency in a short period of time, and without increasing signal length, thereby provides new way for harmonic wave fast detecting between ultra-low frequency.
(2) algorithm and Prony algorithm comparison herein
10 signal cycles that specify taking IEC standard are as example, between harmonic frequency change at 0 ~ 10Hz, its amplitude coefficient (with respect to fundamental frequency) is 1%, carries out 100 Monte Carlo simulation experiments as shown in table 3 with the contrast of Prony algorithm testing result under different noise circumstances:
Each algorithm accuracy of detection contrast situation when table 3 noise
Note: "---" represents that corresponding algorithm cannot detect result.
In different noise circumstances, the accuracy of detection variation tendency of Prony algorithm and herein algorithm is as shown in 6:
Can be found out by table 3 and Fig. 6, Prony algorithm, in noise-free environment, has quite high accuracy of detection and (reaches
the order of magnitude), and along with the increase (being that signal to noise ratio (S/N ratio) diminishes) of noise, its accuracy of detection sharply drops to and cannot detect (noise that is approximately higher than 50dB intensity).And in practical power systems because the voltage of electrical network and frequency are always in continuous variation, sampled signal is introduced random noise unavoidablely simultaneously, Prony algorithm is very sensitive to noise signal, is not suitable for the detection of harmonic wave between electric system ultralow frequency.It can also be seen that, algorithm noise immunity is better herein simultaneously, and in the time that noise strengthens, the variation of testing result precision not quite (is stabilized in
~
the order of magnitude), can meet the demand of the actual detection of electric system.
In sum, the short time fast detecting that can realize harmonic wave between ultralow frequency of the present invention, this algorithm is not subject to the impact of DC component and insensitive to system noise simultaneously, can ensure robustness and the real-time of this algorithm in reality detects.
The foregoing is only preferred embodiment of the present invention, all equalizations of doing according to the present patent application the scope of the claims change and modify, and all should belong to covering scope of the present invention.
Claims (3)
1. a ultralow frequency harmonic detection method of considering negative frequency impact, is characterized in that, comprises the following steps:
Step 1: under operation of power networks normal condition, set up the unifrequency discrete signal model of power network signal:
(1) wherein
for the amplitude of power network signal,
for the frequency of power network signal and
for the initial phase of power network signal,
,
for the total sampling number of discrete signal intercepting,
for power network signal sample frequency;
Step 2: right
the normalization discrete Fourier transform (DFT) (Discrete Fourier Transform, DFT) that adds rectangular window obtains
(2), wherein k=0,1 ... N-1;
Step 3: the frequency resolution of establishing spectral line is
,
, and will in formula (1) substitution formula (2), obtain
X（k）
（3），
If normalized frequency is
,
, formula (3) is further reduced to:
(4), in formula
;
Step 3: get adjacent three spectral lines
,
with
spectrum value
,
with
, here
for integer and to meet its corresponding spectral line be maximum spectral line, here
for integer and to meet its corresponding spectral line be maximum spectral line, obtain through calculating abbreviation
(5)
Note
,
, have according to formula (4):
(6), wherein
=-1,0,1, three unknown numbers in formula are
,
with
, when
=-1,0,1, can become a ternary linear function group by simultaneous, therefrom can obtain
,
with
, wherein solve frequency normalization to be asked
expression formula be
(7)；
Step 4: obtain the highest spectral line by formula (4)
equation:
(8),
In formula
,
,
, formula (9) is got to conjugation and obtains:
(9) obtain
after, simultaneous formula (8) and formula (9) composition system of equations, this system of equations the inside only has
with
be unknown quantity, all the other are known quantity, therefore can solve:
(10);
Step 5: formula (7) and formula (10) are considers humorous wave frequency, amplitude and phase estimation formula between negative frequency and interpolation ultralow.
2. ultralow frequency simple harmonic wave detection method according to claim 1, is characterized in that: the detection in full frequency band time harmonic frequency comprises the following steps:
1) signal electrical network being collected, intercepts one section with rectangular window and makes FFT(DFT) spectrum analysis;
2) for harmonic component between more than 50Hz, the regulation detecting according to harmonic wave between IEC61000-4-7 standard pair, with the form of a harmonic wave group or subgroup, measure frequency range between two harmonic waves and between first-harmonic and second harmonic between harmonic parameters, and frequency using the mean value of two harmonic wave sequence numbers as corresponding harmonic wave group or subgroup;
3) must consider the impact of negative frequency for harmonic wave between 0 ~ 50Hz frequency range; Between the ultralow frequency of 0 ~ 10Hz frequency range, harmonic wave calculates harmonic parameters between ultralow frequency according to described ultralow frequency simple harmonic wave detection method: amplitude
, frequency
and phase place
, and 10 ~ 50Hz frequency range between harmonic wave calculate all spectral lines of residue except first-harmonic spectral line according to the method for IEC, obtain removing the root-mean-square value of amplitude after the intersegmental harmonic wave leakage contribution of ultralow frequency
;
4)
, be 0 ~ 50Hz subharmonic frequency range between harmonic amplitude parameter, according to the regulation of IEC standard, be designated as harmonic wave between 0.5 time;
5) comprehensive step 2) and step 4) obtain full frequency band between harmonic parameters.
3. ultralow frequency harmonic detection method according to claim 1, is characterized in that: between ultralow frequency, harmonic wave produces
spectral line drops on frequency just
when upper, can select
for-1, corresponding
, wherein " * " represent get conjugation, at this moment, the spectral line sequence number of utilization be (1,1,2), correspondence spectrum value be (
,
,
), formula
can be reduced to:
(11).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410317995.4A CN104076203B (en) | 2014-07-07 | 2014-07-07 | A kind of intrasonic harmonic detection method considering that negative frequency affects |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410317995.4A CN104076203B (en) | 2014-07-07 | 2014-07-07 | A kind of intrasonic harmonic detection method considering that negative frequency affects |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104076203A true CN104076203A (en) | 2014-10-01 |
| CN104076203B CN104076203B (en) | 2016-11-02 |
Family
ID=51597599
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410317995.4A Active CN104076203B (en) | 2014-07-07 | 2014-07-07 | A kind of intrasonic harmonic detection method considering that negative frequency affects |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104076203B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104820131A (en) * | 2015-04-29 | 2015-08-05 | 江苏大学 | Method for accurately identifying ultralow frequency signal through dual calculation |
| CN112448397A (en) * | 2020-12-08 | 2021-03-05 | 国网福建省电力有限公司电力科学研究院 | Method and system for rapidly treating inter-harmonic waves of power supply system with variable fundamental frequency |
| CN113608023A (en) * | 2021-08-18 | 2021-11-05 | 国网湖南省电力有限公司 | Sampling frequency self-adaptive harmonic detection method and system |
| CN113702702A (en) * | 2021-08-23 | 2021-11-26 | 中国电力科学研究院有限公司 | Multi-inter-harmonic identification method for improved spectrum superposition algorithm |
| CN115166355A (en) * | 2022-06-17 | 2022-10-11 | 江苏电力信息技术有限公司 | Power system harmonic signal frequency estimation method based on unbiased DFT |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100161259A1 (en) * | 2007-08-29 | 2010-06-24 | Powertron Engineering Co., Ltd. | Aging status diagnostic apparatus for power conversion system, and method their of |
| CN203149014U (en) * | 2013-02-28 | 2013-08-21 | 上海同盛电力有限公司 | Inter-harmonic wave detection apparatus for electric power distribution system |
| CN103575987A (en) * | 2013-11-25 | 2014-02-12 | 国家电网公司 | Inter-harmonic wave detection instrument based on DSP and detection method thereof |
| CN103592514A (en) * | 2012-08-15 | 2014-02-19 | 王景芳 | Novel harmonic high-precision detection method |
| CN103743949A (en) * | 2014-01-06 | 2014-04-23 | 国家电网公司 | Detection method of harmonic and inter-harmonic based on single-channel FastICA (Fast Independent Component Analysis) |
-
2014
- 2014-07-07 CN CN201410317995.4A patent/CN104076203B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100161259A1 (en) * | 2007-08-29 | 2010-06-24 | Powertron Engineering Co., Ltd. | Aging status diagnostic apparatus for power conversion system, and method their of |
| CN103592514A (en) * | 2012-08-15 | 2014-02-19 | 王景芳 | Novel harmonic high-precision detection method |
| CN203149014U (en) * | 2013-02-28 | 2013-08-21 | 上海同盛电力有限公司 | Inter-harmonic wave detection apparatus for electric power distribution system |
| CN103575987A (en) * | 2013-11-25 | 2014-02-12 | 国家电网公司 | Inter-harmonic wave detection instrument based on DSP and detection method thereof |
| CN103743949A (en) * | 2014-01-06 | 2014-04-23 | 国家电网公司 | Detection method of harmonic and inter-harmonic based on single-channel FastICA (Fast Independent Component Analysis) |
Non-Patent Citations (1)
| Title |
|---|
| 梅永 等: "电力参数间谐波测量算法的研究", 《计量技术》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104820131A (en) * | 2015-04-29 | 2015-08-05 | 江苏大学 | Method for accurately identifying ultralow frequency signal through dual calculation |
| CN104820131B (en) * | 2015-04-29 | 2017-01-18 | 江苏大学 | Method for accurately identifying ultralow frequency signal through dual calculation |
| CN112448397A (en) * | 2020-12-08 | 2021-03-05 | 国网福建省电力有限公司电力科学研究院 | Method and system for rapidly treating inter-harmonic waves of power supply system with variable fundamental frequency |
| CN113608023A (en) * | 2021-08-18 | 2021-11-05 | 国网湖南省电力有限公司 | Sampling frequency self-adaptive harmonic detection method and system |
| CN113702702A (en) * | 2021-08-23 | 2021-11-26 | 中国电力科学研究院有限公司 | Multi-inter-harmonic identification method for improved spectrum superposition algorithm |
| CN113702702B (en) * | 2021-08-23 | 2024-01-23 | 中国电力科学研究院有限公司 | Multi-inter-harmonic identification method for improving frequency spectrum superposition algorithm |
| CN115166355A (en) * | 2022-06-17 | 2022-10-11 | 江苏电力信息技术有限公司 | Power system harmonic signal frequency estimation method based on unbiased DFT |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104076203B (en) | 2016-11-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103245832B (en) | Based on harmonic wave time-frequency characteristic method for parameter estimation and the analyser of quick S-transformation | |
| Wen et al. | Simple interpolated FFT algorithm based on minimize sidelobe windows for power-harmonic analysis | |
| CN107247182B (en) | Inter-harmonic component reduction method based on measured phasor data | |
| CN104076203A (en) | Ultralow-frequency inter-harmonic detection method with negative frequency influence considered | |
| CN105723229A (en) | Synchronous phasor measurement method applicable to P-type phasor measurement unit (PMU) | |
| CN103018555B (en) | High-precision electric power parameter software synchronous sampling method | |
| Jain et al. | An adaptive time-efficient technique for harmonic estimation of nonstationary signals | |
| Zhang et al. | Improved smoothing frequency shifting and filtering algorithm for harmonic analysis with systematic error compensation | |
| CN110389312B (en) | Calibrator phasor measurement method suitable for field PMU test | |
| CN104316768B (en) | A kind of negative sequence impedance parameter evaluation method of three-phase imbalance disturbance source locating | |
| Chen et al. | Nonstationary signal denoising using an envelope-tracking filter | |
| CN109507480A (en) | A kind of harmonic detection method and device of neighbouring fundamental wave/harmonic wave | |
| CN104502707A (en) | Synchronized phasor measurement method for electrical power system based on cubic spline interpolation | |
| Ma et al. | Harmonic and interharmonic analysis of mixed dense frequency signals | |
| Xu et al. | An improved Taylor weighted least squares method for estimating synchrophasor | |
| Zhang et al. | Power system dynamic frequency measurement based on novel interpolated STFT algorithm | |
| CN109444537B (en) | Self-adaptive synchronous phasor measurement method considering out-of-band interference | |
| CN109581045A (en) | A kind of m-Acetyl chlorophosphonazo power measurement method meeting IEC standard frame | |
| Bertocco et al. | Frequency tracking for efficient phasor measurement based on a CSTFM model | |
| CN114487589A (en) | Power grid broadband signal self-adaptive measurement method, device and system | |
| CN115859083A (en) | Decoupling method for multi-vibration-source coupled vibration signals in hydroelectric generating set | |
| CN112034232A (en) | Power supply system voltage sag detection method | |
| Jia et al. | A novel adaptive linear prediction-based parameter estimation method for monitoring sub-/inter-harmonics during SSI events | |
| Belega et al. | Estimation of fundamental power quantities in single-phase systems under dynamic conditions | |
| Men et al. | A new time domain filtering method for calculating the RMS value of vibration signals |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20170503 Address after: 100031 Xicheng District West Chang'an Avenue, No. 86, Beijing Co-patentee after: State Grid Fujian Electric Power Co., Ltd. Patentee after: State Grid Corporation of China Co-patentee after: Electric Power Research Institute, State Grid Fujian Electric Power Co., Ltd. Co-patentee after: Fujian EPRI Power Commissioning Co., Ltd. Address before: 100031 Xicheng District West Chang'an Avenue, No. 86, Beijing Co-patentee before: State Grid Fujian Electric Power Co., Ltd. Patentee before: State Grid Corporation of China Co-patentee before: Electric Power Research Institute, State Grid Fujian Electric Power Co., Ltd. |