CN107390022A - Electric energy gauging method based on discrete spectrum correction - Google Patents

Electric energy gauging method based on discrete spectrum correction Download PDF

Info

Publication number
CN107390022A
CN107390022A CN201710729394.8A CN201710729394A CN107390022A CN 107390022 A CN107390022 A CN 107390022A CN 201710729394 A CN201710729394 A CN 201710729394A CN 107390022 A CN107390022 A CN 107390022A
Authority
CN
China
Prior art keywords
mrow
harmonic
msub
wave
spectrum
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
Application number
CN201710729394.8A
Other languages
Chinese (zh)
Other versions
CN107390022B (en
Inventor
周杰文
任智仁
汪龙峰
黄杰
杨辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wasion Group Co Ltd
Original Assignee
Wasion Group Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wasion Group Co Ltd filed Critical Wasion Group Co Ltd
Priority to CN201710729394.8A priority Critical patent/CN107390022B/en
Publication of CN107390022A publication Critical patent/CN107390022A/en
Application granted granted Critical
Publication of CN107390022B publication Critical patent/CN107390022B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R22/00Arrangements for measuring time integral of electric power or current, e.g. electricity meters
    • G01R22/06Arrangements for measuring time integral of electric power or current, e.g. electricity meters by electronic methods
    • G01R22/10Arrangements for measuring time integral of electric power or current, e.g. electricity meters by electronic methods using digital techniques

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Current Or Voltage (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

The invention discloses a kind of electric energy gauging method based on discrete spectrum correction, including the voltage and current signal of digital quantity is converted into using a length of N of window tetra- three rank windows progress time domain truncations of Nuttall;The spectrum barycenter and fundamental frequency of fundamental wave are calculated using spectrum barycenter formula;Harmonic spike position is estimated, calculates harmonic spectrum barycenter and the maximum spectrum peak position and offset of fundamental wave and each harmonic;Principal wave harmonic wave amplitude correction coefficient is determined, and calculates tested voltage, the fundamental wave of electric current, harmonic amplitude and phase;The electrical parameter of fundamental wave harmonic is calculated, completes electric energy metrical.Therefore the inventive method need not carry out peak value searching, offset need not be calculated according to the type of institute's windowed function, and mutually eliminate zero-crossing examination, it is directly calculated using fft analysis method, eliminate the window spectrum computing of complexity, calculate and analyze obtained frequency accuracy to greatly improve, and program realizes that simply real-time is good.

Description

Electric energy gauging method based on discrete spectrum correction
Technical field
Present invention relates particularly to a kind of electric energy gauging method based on discrete spectrum correction.
Background technology
With the development and the improvement of people's living standards of economic technology, electric energy has become the daily production of people and life Essential secondary energy sources in work, production and life to people bring endless facility.
Harmonious Waves in Power Systems has significant effect to the accuracy of electric energy metrical, using metering fundamental energy harmonic respectively The metering of electric energy is a kind of relatively rational mode, and with the development of technology, requirement of the market to principal wave harmonic wave measuring accuracy is more next It is higher.According to the current development in electric energy metrical field, the design of harmonic wave metering mainly has two kinds:A kind of is special metering Chip, another kind are the special AD+MCU of DSP+.Former scheme cost is relatively low, is mainly used in the design of common list, three-phase meter In;The latter's cost is of a relatively high, and main application is mainly used in regional grid, provincial power network, light with Source of Gateway Meter platform, Source of Gateway Meter The stoichiometric point that volt, wind-powered electricity generation, water power and electric railway even load frequently change, the field changeable in face of the big load of this harmonic content Close, the program is more flexible compared to the former function, it is easier to meets the needs of complicated.
Source of Gateway Meter does not require nothing more than higher electric energy metrical precision, and require lower fortune as a ripe product Row power consumption.And many algorithms have been proposed currently for the high-precision measuring of fundamental wave, the amplitude of harmonic wave, phase and frequency, this its In have quite a few algorithm because realizing complicated difficult to accomplish requirement of real-time in embedded systems, and realize it is more complicated, System operation power consumption can be improved accordingly, thus finds a kind of measuring accuracy height, realized simply, and time and spatial complexity It can meet that the algorithm of existing hardware platform is particularly important.
The content of the invention
It is an object of the invention to provide a kind of precision is of a relatively high, real-time is good and algorithm is simply based on discrete spectrum The electric energy gauging method of correction.
This electric energy gauging method based on discrete spectrum correction provided by the invention, comprises the following steps:
S1. it is digital quantity by sample conversion by voltage and current signal, using a length of N of window tetra- three ranks of Nuttall Window carries out time domain truncation to the digital quantity after conversion, obtains blocking the spectrum signal Y (k) of rear signal;The spectrum signal Y (k) Including voltage spectrum value YuAnd current spectrum value Y (k)i(k);
S2. the voltage spectrum value Y obtained according to step S1u(k) the spectrum barycenter of fundamental wave, is calculated using spectrum barycenter formula ksc1With fundamental frequency f0
S3. the fundamental frequency f obtained using step S20Estimate harmonic spike positionAnd calculated using spectrum barycenter formula Obtain harmonic spectrum barycenter ksch, and the maximum spectrum peak position k of fundamental wave and each harmonic is further calculatedhWith offset λh
S4. principal wave harmonic wave amplitude correction coefficient g (λ are determined according to the offset obtained in step S3h), and be calculated by Survey voltage, the fundamental wave of electric current, harmonic amplitude and phase;
S5. tested voltage, the fundamental wave of electric current, harmonic amplitude and the phase calculation fundamental wave harmonic obtained according to step S4 Electrical parameter, so as to complete the metering of electric energy.
The spectrum signal Y (k) for obtaining blocking rear signal described in step S1, after specially being blocked using fft analysis The spectrum signal Y (k) of signal.
Tetra- three rank windows of Nuttall of a length of N of window described in step S1, the specially main lobe of tetra- three rank windows of Nuttall Width isThere are 8 spectral lines in main lobe.
The a length of N of described window tetra- three rank windows of Nuttall meet following formula:
N=1,2 in formula ..., N-1;bmFor window function coefficient, meet
Fundamental wave spectrum barycenter k is calculated described in step S2sc1With fundamental frequency f0, specially calculated using following steps Obtain:
A. fundamental wave spectrum barycenter k is calculated according to spectrum barycenter formulasc1
In formulaAndFor downward bracket function, fmin is The minimum value of signal fundamental frequency, fmax be signal fundamental frequency maximum, Δ f=fs/ N, fs are the sampling of analog-to-digital conversion Frequency, N are that the analysis of fft algorithm is counted;
B. the fundamental wave obtained according to step A composes barycenter ksc1Measure fundamental frequency f0
f0=ksc1*Δf。
Harmonic spectrum barycenter k is calculated described in step S3sch, fundamental wave and each harmonic maximum spectrum peak position khWith it is inclined Shifting amount λh, specially it is calculated using following steps:
A. using following formula estimation harmonic spectrum peak
Wherein h value is 2,3 ..., M;M is the maximum overtone order of signal;
B. obtained according to step aHarmonic spectrum barycenter k is calculated using following formulasch
C. the harmonic spectrum barycenter k obtained according to step bsch, fundamental wave and each harmonic are calculated most using equation below Big spectrum peak position kh
kh=round (ksch)
Round () is the function that rounds up in formula;
D. the maximum spectrum peak position k obtained according to step chOffset λ is calculatedh
λh=kh-ksch
Determination principal wave harmonic wave amplitude correction coefficient g (λ described in step S4h), specially calculated using equation below:
Wherein p1~p5It is constant.
Tested voltage, the fundamental wave of electric current, harmonic amplitude and phase is calculated described in step S4, specially using as follows Updating formula is corrected to amplitude and phase:
A=| Y (kh)|g(λh)
A is revised amplitude in formula, and θ is revised phase, angle (Y (kh)) represent Y (kh) phase value.
The electrical parameter of calculating fundamental wave harmonic described in step S5, specially calculates fundamental wave harmonic using following formula Electrical parameter:
Define principal wave harmonic wave voltage magnitude Uh, voltage-phase θ uh, current amplitude Ih, current phase θ ih, h=1,2 ..., 63;
Phase difference:θh=θ uh-θih
Principal wave harmonic wave active power:Ph=UhIh cos(θuh-θih)
Principal wave harmonic wave reactive power:Qh=UhIh sin(θuh-θih)
Full wave voltage:
Full-wave electric current:
All-wave is active:
All-wave is idle:
Fundamental active electric flux:EP=P1T, t are the time;
Fundamental wave reactive power electric flux:EQ=Q1T, t are the time;
Harmonic wave active energy amount:T is the time;
Harmonic wave reactive energy amount:T is the time;
All-wave active energy amount:EP=Pt, t are the time;
All-wave reactive energy amount:EQ=Qt, t are the time.
This electric energy gauging method based on discrete spectrum correction provided by the invention, is accurately looked for using spectrum centroid method first Go out signal spectrum peak value, to eliminate caused spectral leakage and fence response under non-synchronous sampling;Further to voltage, electric current Spectrum value be corrected using amplitude correction formula and phase only pupil filter formula, finally try to achieve all-wave, fundamental wave harmonic electric parameter And electric flux;Therefore the inventive method need not carry out peak value searching, it is not necessary to calculated and offset according to the type of institute's windowed function Amount, and zero-crossing examination is mutually eliminated, it is directly calculated using fft analysis method, eliminates the window spectrum computing of complexity, meter Calculate and analyze obtained frequency accuracy to greatly improve, and program realizes that simply real-time is good.
Brief description of the drawings
Fig. 1 is the method flow diagram of the inventive method.
Fig. 2 is the spectrum centroid calculation schematic diagram of the inventive method.
Fig. 3 is the correction coefficient matched curve of the inventive method.
Fig. 4 is the harmonic frequency error precision of the inventive method.
Fig. 5 is the harmonic amplitude error precision of the inventive method.
Fig. 6 is the harmonic phase error precision of the inventive method.
Embodiment
It is the method flow diagram of the inventive method as shown in Figure 1:It is provided by the invention it is this based on discrete spectrum correction Electric energy gauging method, comprise the following steps:
S1. it is digital quantity by sample conversion by voltage and current signal, using a length of N of window tetra- three ranks of Nuttall Window carries out time domain truncation to the digital quantity after conversion, obtains blocking the spectrum signal Y (k) of rear signal using fft analysis;The frequency Spectrum signal Y (k) includes voltage spectrum value YuAnd current spectrum value Y (k)i(k);
In the specific implementation, the main lobe width of preferably tetra- three rank windows of Nuttall isThere are 8 spectral lines in main lobe;
The a length of N of described window tetra- three rank windows of Nuttall meet following formula:
N=1,2 in formula ..., N-1;bmFor window function coefficient, meetb0= 0.338946, b1=0.481973, b2=0.161054, b3=0.018027;
Described application window function is to do following processing to discrete voltage u (n), current signal i (n):U ' (n)=u (n) * w (n), i ' (n)=i (n) * w (n);U ' (n) and i ' (n) obtains discrete spectrum Y after fft analysisuAnd Y (k)i(k), for ease of point Analysis, it is used uniformly Y (k) expressions;
S2. the voltage spectrum value Y obtained according to step S1u(k), using spectrum barycenter formula to calculate, (Fig. 2 is the inventive method Spectrum centroid calculation schematic diagram) obtain the spectrum barycenter k of fundamental wavesc1With fundamental frequency f0;Specially calculated using following steps Arrive:
A. fundamental wave spectrum barycenter k is calculated according to spectrum barycenter formulasc1
In formulaAndFor downward bracket function, fmin is The minimum value of signal fundamental frequency, specially 45Hz;Fmax be signal fundamental frequency maximum, specially 65Hz;Δ f= fs/ N, fs are the sample frequency of analog-to-digital conversion, and N is that the analysis of fft algorithm is counted;
B. the fundamental wave obtained according to step A composes barycenter ksc1Measure fundamental frequency f0
f0=ksc1*Δf
S3. the fundamental frequency f obtained using step S20Estimate harmonic spike positionAnd calculated using spectrum barycenter formula Obtain harmonic spectrum barycenter ksch, and the maximum spectrum peak position k of fundamental wave and each harmonic is further calculatedhWith offset λh;Tool Body is to be calculated using following steps:
A. using following formula estimation harmonic spectrum peak
Wherein h value is 1,2,3 ..., M;M is the maximum overtone order of signal;
B. obtained according to step aHarmonic spectrum barycenter k is calculated using following formulasch
C. the harmonic spectrum barycenter k obtained according to step bsch, fundamental wave and each harmonic are calculated most using equation below Big spectrum peak position kh
kh=round (ksch)
Round () is the function that rounds up in formula;
D. offset λ is calculated in the maximum spectrum peak position kh obtained according to step ch
λh=kh-ksch
S4. principal wave harmonic wave amplitude correction coefficient g (λ are determined according to the offset obtained in step S3h), and be calculated by Survey voltage, the fundamental wave of electric current, harmonic amplitude and phase;
Calculating g (λ are carried out using equation belowh):
Wherein p1~p5It is constant;Preferable p1=0.2073, p2=-0.01546, p3=0.9941, p4=- 0.000251, p5=2.9503;
Amplitude and phase are corrected using following updating formula:
A=| Y (kh)|g(λh)
A is revised amplitude in formula, and θ is revised phase, angle (Y (kh)) represent Y (kh) phase value;
The derivation of above-mentioned amplitude correction coefficient described briefly below, amplitude correction formula and phase only pupil filter formula:
Consider single-frequency signals y (n) with sample frequency f firstsThe discrete-time signal that uniform sampling obtains is:
A is amplitude in formula, f0For frequency,For initial phase;
Y is obtained to y (n) plus tetra- three rank windows of Nuttallw(n)=y (n) w (n), yw(n) continuous fourier transform is
Discrete sampling is carried out to above formula, and ignores negative frequency point-f0Going out the secondary lobe at frequency peak influences, and obtains signal after adding window The expression formula of discrete Fourier transform is:
Make f0=kscΔ f, then kscFor the spectrum barycenter of signal;As shown in Fig. 2 k*For close to spectrum barycenter kscMaximum spectral peak Position, therefore can obtain:
Make λ=k*-ksc, λ is offset, λ ∈ [- 0.5,0.5];
Further derive that the discrete spectrum approximation to function of tetra- three rank windows of Nuttall is as follows:
Amplitude correction formula, which can to sum up be obtained, is:
In formulaFor window function offset amplitude;Curve map it is as shown in Figure 3;Therefrom It can be seen that the curve can use fitting of a polynomial, to reduce operand, first approached using fourth order polynomial fitting formulaCurve, fitting formula are as follows:
G (λ)=p1λ4+p2|λ|3+p3λ2+p4|λ|+p5
Wherein p1=0.2073, p2=-0.01546, p3=0.9941, p4=-0.000251, p5=2.9503;
Therefore amplitude formula can be obtained is:
A=| Y (k*)|g(λ)
Phase only pupil filter formula is:
S5. tested voltage, the fundamental wave of electric current, harmonic amplitude and the phase calculation fundamental wave harmonic obtained according to step S4 Electrical parameter, so as to complete the metering of electric energy;The electrical parameter of fundamental wave harmonic is specially calculated using following formula:
Define principal wave harmonic wave voltage magnitude Uh, voltage-phase θ uh, current amplitude Ih, current phase θ ih, h=1,2 ..., 63;
Phase difference:θh=θ uh-θih
Principal wave harmonic wave active power:Ph=UhIh cos(θuh-θih)
Principal wave harmonic wave reactive power:Qh=UhIh sin(θuh-θih)
Full wave voltage:
Full-wave electric current:
All-wave is active:
All-wave is idle:
Fundamental active electric flux:EP=P1T, t are the time;
Fundamental wave reactive power electric flux:EQ=Q1T, t are the time;
Harmonic wave active energy amount:T is the time;
Harmonic wave reactive energy amount:T is the time;
All-wave active energy amount:EP=Pt, t are the time;
All-wave reactive energy amount:EQ=Qt, t are the time.
The electric energy gauging method provided below by way of verification experimental verification invention:
The inventive method be applied to Source of Gateway Meter (ammeter specifications and models are DTSD341-MA2,3 × 57.7V, 3 × 1.5 (6) A, 20000imp/kWh, 50Hz), program is implemented in meter platform Blackfin BF533 DSP process chips.
The fundamental active pulse elementary error precision of this method is as shown in Figure 1:
The fundamental active pulse elementary error precision of table 1
Harmonic pulse Setup Experiments:
Fundamental voltage Un=57.7V;Fundamental current In is 0.5*Imax, as In=2*Ib=3A;Phasor power factor For 1.0;Fundamental frequency is 50Hz;Single harmonic component is superimposed on the basis of fundamental wave, overtone order highest is added to 48 times;Harmonic voltage Amplitude 0.05*Un;Harmonic current 0.4*In;Harmonic power factor is 0.5L;Harmonic pulse constant 1000000imp/kWh;
It is superimposed below for Source of Gateway Meter under single harmonic component, the error that harmonic wave goes out pulse is as shown in table 2:
The harmonic wave of table 2 goes out the error signal table of pulse
From upper table it can be found that Source of Gateway Meter measured data greatly differs from each other with theoretical value, this is due to that experiment condition is limited, is made Obtaining experimental data can not reach best, and nonetheless, upper table error information has met GB/T17215.302-2013 pair harmonic wave Electric energy metrical precision peak demand, and entered every strict test, fully meet standard requirement.Harmonic pulse error is main Influenceed by caused by hardware sample circuit and transformer than difference and difference, in table harmonic error data be ammeter through amplitude and Data are surveyed after phasing.After experiment condition permission, above-mentioned data precision will have the lifting of matter.
MATLAB simulation results as shown in Fig. 3~Fig. 6, respectively correction coefficient matched curve, harmonic frequency error precision, Harmonic amplitude error precision harmonic phase error precision;It can be seen that the present invention, which carries algorithm, can realize that harmonic wave is believed Number high-precision measuring.

Claims (9)

1. a kind of electric energy gauging method based on discrete spectrum correction, comprises the following steps:
S1. it is digital quantity by sample conversion by voltage and current signal, using a length of N of window tetra- three rank windows pair of Nuttall Digital quantity after conversion carries out time domain truncation, obtains blocking the spectrum signal Y (k) of rear signal;The spectrum signal Y (k) includes Voltage spectrum value YuAnd current spectrum value Y (k)i(k);
S2. the voltage spectrum value Y obtained according to step S1u(k) the spectrum barycenter k of fundamental wave, is calculated using spectrum barycenter formulasc1With Fundamental frequency f0
S3. the fundamental frequency f obtained using step S20Estimate harmonic spike positionAnd it is calculated using spectrum barycenter formula Harmonic spectrum barycenter ksch, and the maximum spectrum peak position k of fundamental wave and each harmonic is further calculatedhWith offset λh
S4. principal wave harmonic wave amplitude correction coefficient g (λ are determined according to the offset obtained in step S3h), and tested electricity is calculated Pressure, the fundamental wave of electric current, harmonic amplitude and phase;
S5. the electricity ginseng of tested voltage, the fundamental wave of electric current, harmonic amplitude and the phase calculation fundamental wave harmonic obtained according to step S4 Number, so as to complete the metering of electric energy.
2. the electric energy gauging method according to claim 1 based on discrete spectrum correction, it is characterised in that described in step S1 The spectrum signal Y (k) for obtaining blocking rear signal, specially obtain blocking the spectrum signal Y (k) of rear signal using fft analysis.
3. the electric energy gauging method according to claim 2 based on discrete spectrum correction, it is characterised in that described in step S1 The a length of N of window tetra- three rank windows of Nuttall, the main lobe width of specially tetra- three rank windows of Nuttall isHave in main lobe 8 spectral lines.
4. the electric energy gauging method according to claim 3 based on discrete spectrum correction, it is characterised in that described window length Meet following formula for N tetra- three rank windows of Nuttall:
<mrow> <mi>w</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>0</mn> </mrow> <mn>3</mn> </munderover> <msup> <mrow> <mo>(</mo> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>m</mi> </msup> <msub> <mi>b</mi> <mi>m</mi> </msub> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mn>2</mn> <mi>&amp;pi;</mi> <mi>n</mi> <mo>*</mo> <mi>m</mi> <mo>/</mo> <mi>N</mi> <mo>)</mo> </mrow> </mrow>
N=1,2 in formula ..., N-1;bmFor window function coefficient, meet
5. the electric energy gauging method according to claim 4 based on discrete spectrum correction, it is characterised in that described in step S2 Fundamental wave spectrum barycenter k is calculatedsc1With fundamental frequency f0, specially it is calculated using following steps:
A. fundamental wave spectrum barycenter k is calculated according to spectrum barycenter formulasc1
<mrow> <msub> <mi>k</mi> <mrow> <mi>s</mi> <mi>c</mi> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <msub> <mi>k</mi> <mi>min</mi> </msub> <mo>-</mo> <mi>B</mi> <mo>/</mo> <mn>2</mn> </mrow> <mrow> <msub> <mi>k</mi> <mi>max</mi> </msub> <mo>+</mo> <mi>B</mi> <mo>/</mo> <mn>2</mn> </mrow> </munderover> <msup> <mi>kY</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </mrow> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <msub> <mi>k</mi> <mi>min</mi> </msub> <mo>-</mo> <mi>B</mi> <mo>/</mo> <mn>2</mn> </mrow> <mrow> <msub> <mi>k</mi> <mi>max</mi> </msub> <mo>+</mo> <mi>B</mi> <mo>/</mo> <mn>2</mn> </mrow> </munderover> <msup> <mi>Y</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow>
In formulaAndFor downward bracket function, fmin is signal base The minimum value of wave frequency rate, fmax be signal fundamental frequency maximum, Δ f=fs/ N, fs be analog-to-digital conversion sample frequency, N Counted for the analysis of fft algorithm;
B. the fundamental wave obtained according to step A composes barycenter ksc1Measure fundamental frequency f0
f0=ksc1*Δf。
6. the electric energy gauging method according to claim 5 based on discrete spectrum correction, it is characterised in that described in step S3 Harmonic spectrum barycenter k is calculatedsch, fundamental wave and each harmonic maximum spectrum peak position khWith offset λh, specially using such as Lower step is calculated:
A. using following formula estimation harmonic spectrum peak
Wherein h value is 2,3 ..., M;M is the maximum overtone order of signal;
B. obtained according to step aHarmonic spectrum barycenter k is calculated using following formulasch
<mrow> <msub> <mi>k</mi> <mrow> <mi>s</mi> <mi>c</mi> <mi>h</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <msubsup> <mi>k</mi> <mi>h</mi> <mo>*</mo> </msubsup> <mo>-</mo> <mi>B</mi> <mo>/</mo> <mn>2</mn> </mrow> <mrow> <msubsup> <mi>k</mi> <mi>h</mi> <mo>*</mo> </msubsup> <mo>+</mo> <mi>B</mi> <mo>/</mo> <mn>2</mn> </mrow> </munderover> <msup> <mi>kY</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </mrow> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <msubsup> <mi>k</mi> <mi>h</mi> <mo>*</mo> </msubsup> <mo>-</mo> <mi>B</mi> <mo>/</mo> <mn>2</mn> </mrow> <mrow> <msubsup> <mi>k</mi> <mi>h</mi> <mo>*</mo> </msubsup> <mo>+</mo> <mi>B</mi> <mo>/</mo> <mn>2</mn> </mrow> </munderover> <msup> <mi>Y</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow>
C. the harmonic spectrum barycenter k obtained according to step bsch, fundamental wave is calculated using equation below and the maximum of each harmonic is composed Peak position kh
kh=round (ksch)
Round () is the function that rounds up in formula;
D. the maximum spectrum peak position k obtained according to step chOffset λ is calculatedh
λh=kh-ksch
7. the electric energy gauging method according to claim 6 based on discrete spectrum correction, it is characterised in that described in step S4 Determination principal wave harmonic wave amplitude correction coefficient g (λh), specially calculated using equation below:
<mrow> <mi>g</mi> <mrow> <mo>(</mo> <msub> <mi>&amp;lambda;</mi> <mi>h</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>p</mi> <mn>1</mn> </msub> <msubsup> <mi>&amp;lambda;</mi> <mi>h</mi> <mn>4</mn> </msubsup> <mo>+</mo> <msub> <mi>p</mi> <mn>2</mn> </msub> <msubsup> <mi>&amp;lambda;</mi> <mi>h</mi> <mn>3</mn> </msubsup> <mo>+</mo> <msub> <mi>p</mi> <mn>3</mn> </msub> <msubsup> <mi>&amp;lambda;</mi> <mi>h</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msub> <mi>p</mi> <mn>4</mn> </msub> <msub> <mi>&amp;lambda;</mi> <mi>h</mi> </msub> <mo>+</mo> <msub> <mi>p</mi> <mn>5</mn> </msub> </mrow>
Wherein p1~p5It is constant.
8. the electric energy gauging method according to claim 7 based on discrete spectrum correction, it is characterised in that described in step S4 Tested voltage, the fundamental wave of electric current, harmonic amplitude and phase is calculated, specially using following updating formula to amplitude and phase Position is corrected:
A=| Y (kh)|g(λh)
<mrow> <mi>&amp;theta;</mi> <mo>=</mo> <mi>a</mi> <mi>n</mi> <mi>g</mi> <mi>l</mi> <mi>e</mi> <mrow> <mo>(</mo> <mi>Y</mi> <mo>(</mo> <msub> <mi>k</mi> <mi>h</mi> </msub> <mo>)</mo> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&amp;lambda;</mi> <mi>h</mi> </msub> <mi>&amp;pi;</mi> <mo>+</mo> <mfrac> <mi>&amp;pi;</mi> <mn>2</mn> </mfrac> </mrow>
A is revised amplitude in formula, and θ is revised phase, angle (Y (kh)) represent Y (kh) phase value.
9. the electric energy gauging method according to claim 8 based on discrete spectrum correction, it is characterised in that described in step S5 Calculating fundamental wave harmonic electrical parameter, specially using following formula calculate fundamental wave harmonic electrical parameter:
Define principal wave harmonic wave voltage magnitude Uh, voltage-phase θ uh, current amplitude Ih, current phase θ ih, h=1,2 ..., 63;
Phase difference:θh=θ uh-θih
Principal wave harmonic wave active power:Ph=UhIhcos(θuh-θih)
Principal wave harmonic wave reactive power:Qh=UhIhsin(θuh-θih)
Full wave voltage:
Full-wave electric current:
All-wave is active:
All-wave is idle:
Fundamental active electric flux:EP=P1T, t are the time;
Fundamental wave reactive power electric flux:EQ=Q1T, t are the time;
Harmonic wave active energy amount:T is the time;
Harmonic wave reactive energy amount:T is the time;
All-wave active energy amount:EP=Pt, t are the time;
All-wave reactive energy amount:EQ=Qt, t are the time.
CN201710729394.8A 2017-08-23 2017-08-23 Electric energy metering method based on discrete spectrum correction Active CN107390022B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710729394.8A CN107390022B (en) 2017-08-23 2017-08-23 Electric energy metering method based on discrete spectrum correction

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710729394.8A CN107390022B (en) 2017-08-23 2017-08-23 Electric energy metering method based on discrete spectrum correction

Publications (2)

Publication Number Publication Date
CN107390022A true CN107390022A (en) 2017-11-24
CN107390022B CN107390022B (en) 2020-04-14

Family

ID=60354344

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710729394.8A Active CN107390022B (en) 2017-08-23 2017-08-23 Electric energy metering method based on discrete spectrum correction

Country Status (1)

Country Link
CN (1) CN107390022B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110795070A (en) * 2019-05-24 2020-02-14 杭州海兴电力科技股份有限公司 Virtual gateway table platform and construction method
CN112230054A (en) * 2019-12-20 2021-01-15 青岛鼎信通讯股份有限公司 Harmonic content and power calculation method
CN113358930A (en) * 2021-08-09 2021-09-07 南京派格测控科技有限公司 Harmonic wave test system, device and method based on signal deviation
CN114167198A (en) * 2021-10-18 2022-03-11 国网山东省电力公司平原县供电公司 Method and platform for measuring synchronous line loss data
CN114184839A (en) * 2021-12-10 2022-03-15 江苏集萃智能集成电路设计技术研究所有限公司 ADC frequency spectrum testing method for inhibiting signal source jitter

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010144857A3 (en) * 2004-12-21 2011-03-24 Montana Tech Of The University Of Montana Method of detecting system function by measuring frequency response
CN102445595A (en) * 2011-09-23 2012-05-09 重庆大学 Real-time measuring method for time-varying power of electrical power system
CN104391178A (en) * 2014-12-05 2015-03-04 国家电网公司 Time shift phase difference steady harmonic signal correction method based on Nuttall window
CN105445541A (en) * 2015-12-24 2016-03-30 合肥工业大学 Method for adaptively calculating power under arbitrary frequencies
CN106405229A (en) * 2016-08-30 2017-02-15 威胜集团有限公司 Fundamental wave and harmonic wave electric energy metering method

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010144857A3 (en) * 2004-12-21 2011-03-24 Montana Tech Of The University Of Montana Method of detecting system function by measuring frequency response
CN102445595A (en) * 2011-09-23 2012-05-09 重庆大学 Real-time measuring method for time-varying power of electrical power system
CN104391178A (en) * 2014-12-05 2015-03-04 国家电网公司 Time shift phase difference steady harmonic signal correction method based on Nuttall window
CN105445541A (en) * 2015-12-24 2016-03-30 合肥工业大学 Method for adaptively calculating power under arbitrary frequencies
CN106405229A (en) * 2016-08-30 2017-02-15 威胜集团有限公司 Fundamental wave and harmonic wave electric energy metering method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
曾金芳: "谱质心在电力谐波分析中的应用", 《中国电机工程学报》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110795070A (en) * 2019-05-24 2020-02-14 杭州海兴电力科技股份有限公司 Virtual gateway table platform and construction method
CN110795070B (en) * 2019-05-24 2023-07-25 杭州海兴电力科技股份有限公司 Virtual gateway table platform and construction method
CN112230054A (en) * 2019-12-20 2021-01-15 青岛鼎信通讯股份有限公司 Harmonic content and power calculation method
CN113358930A (en) * 2021-08-09 2021-09-07 南京派格测控科技有限公司 Harmonic wave test system, device and method based on signal deviation
CN114167198A (en) * 2021-10-18 2022-03-11 国网山东省电力公司平原县供电公司 Method and platform for measuring synchronous line loss data
CN114167198B (en) * 2021-10-18 2024-03-01 国网山东省电力公司平原县供电公司 Method and platform for measuring synchronous line loss data
CN114184839A (en) * 2021-12-10 2022-03-15 江苏集萃智能集成电路设计技术研究所有限公司 ADC frequency spectrum testing method for inhibiting signal source jitter
CN114184839B (en) * 2021-12-10 2024-03-01 江苏集萃智能集成电路设计技术研究所有限公司 ADC frequency spectrum testing method for inhibiting signal source jitter

Also Published As

Publication number Publication date
CN107390022B (en) 2020-04-14

Similar Documents

Publication Publication Date Title
CN107390022A (en) Electric energy gauging method based on discrete spectrum correction
CN102288807B (en) Method for measuring electric network voltage flicker
CN109633262A (en) Three phase harmonic electric energy gauging method, device based on composite window multiline FFT
CN103454497B (en) Based on the method for measuring phase difference improving windowed DFT
CN103344825B (en) Based on the electric energy metered system of AC sampling
CN106405229B (en) A kind of fundamental wave and harmonic wave electric energy gauging method
CN108318852B (en) Square wave influence test method for intelligent electric energy meter
CN104360235B (en) Method for positioning voltage sag source of complex power distribution network including DGs
CN104267258B (en) A kind of harmonic wave instantaneous power computational methods utilizing incomplete S-transformation
CN103399204A (en) Rife-Vincent (II) window interpolation FFT (Fast Fourier Transform)-based harmonic and inter-harmonic detection method
CN105486921A (en) Kaiser third-order mutual convolution window triple-spectrum-line interpolation harmonic wave and inter-harmonic wave detection method
CN103389412B (en) The electrical network electric parameter measurement method and device that a kind of small generator is constituted
CN102095929A (en) Method for rapidly measuring frequency of alternating-current signals
CN109061345A (en) Effective value measurement method and system suitable for electric system
CN104375014B (en) Single-frequency overlapping method for power capacitor noise test
CN107167658B (en) A kind of jamproof electric system fundamental frequency of high-precision and Method for Phase Difference Measurement
CN104407197A (en) Signal phasor measurement method based on trigonometric function iteration
CN106771552B (en) A kind of distortion power metering method
CN108051682B (en) Verification method of single-phase rectifier system impedance model
CN102749488A (en) Power grid harmonic wave real-time on-line monitor and method for detecting harmonic wave using same
CN104360198A (en) High-performance power quality analyzer
CN104865480B (en) A kind of method for obtaining load nonlinear function
CN103105532A (en) Harmonic electric energy measuring system of improved Fourier transform and control method thereof
CN203012015U (en) Harmonic energy metering system improving Fourier transformation
CN106771584A (en) A kind of harmonic detecting method for being applied to power distribution network

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant