CN101915872B - Nonlinear-load electric energy measuring method - Google Patents
Nonlinear-load electric energy measuring method Download PDFInfo
- Publication number
- CN101915872B CN101915872B CN2010102732345A CN201010273234A CN101915872B CN 101915872 B CN101915872 B CN 101915872B CN 2010102732345 A CN2010102732345 A CN 2010102732345A CN 201010273234 A CN201010273234 A CN 201010273234A CN 101915872 B CN101915872 B CN 101915872B
- Authority
- CN
- China
- Prior art keywords
- formula
- current
- value
- electric weight
- nonlinear
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Measurement Of Current Or Voltage (AREA)
Abstract
The invention discloses a nonlinear-load measuring method having a high integral algebraic accuracy grade. Aiming at the traditional dot product measuring algorithm, the invention points out the algorithm accuracy deficiency under the condition of calculating nonlinear loads based on the analysis of algorithm characteristics, thereby providing a scheme to improve the accuracy and the adaptability of an electric energy meter on the modern nonlinear-load electric energy measurement by using a complex Newton-Cotes integral algorithm having a high-order measuring algebraic accuracy grade. Based on the theoretical analysis of the algorithm and the real-time operation characteristics of the current measuring system, the invention designs and provides a method for the real-time operation of the algorithm. Simulation and experiment results prove that the invention can effectively improve the accuracy of an electric meter in the electric energy measuring field, thereby improving the adaptability of the electric meter on wider nonlinear-load measurement.
Description
Technical field
The invention belongs to electric system electric energy metrical field.
Background technology
Electric energy metrical is the basic basis that carries out economic balance between electricity power enterprise, power transmission and transformation enterprise and the electricity consumption client, and its accuracy directly influences three's interests.In the process of modern chemical industry industry develop rapidly, electric system has inserted a large amount of nonlinear-loads, and consequent a large amount of harmonic currents, harmonic voltage have not only had a strong impact on power transmission efficiency, and the accuracy of electric energy metrical produced seriously influences.
For effectively solving the influence of mains by harmonics to the electric energy metrical accuracy, each experts and scholars has launched the research to the nonlinear load electric energy metrical one after another.At present, industry generally solves the problems such as accuracy, rationality of nonlinear-load metering through research part throttle characteristics, electrical network modeling and based on the methods such as frequency analysis of FFT.
At present, industry generally adopts dot-product operation to calculate to be had, capacity of idle power, has, reactive power and parameters such as voltage, current effective value.Wherein concrete active energy computing formula is as follows:
is instantaneous power in the formula;
is electric flux in the signal period;
is maximum sub-range length;
is the signal period, and
is the signal sampling time interval.It is thus clear that dot-product operation directly prolongs the basic definition formula of having used integral operation.But in calculating process, how many formulas need count to participate in calculating could be satisfied the accuracy of calculating, and formula does not provide guidance.So, when being input as standard pure sinusoid signal, the accuracy of some integration method is different under the unlike signal sampling rate.
Be calculated as example with active energy; If input voltage, current signal are:
,
; In the formula
;
is input signal cycle, can get:
(2)
Substitution formula (1), then the electric energy theoretical value is in the signal period:
(4)
Substitution formula (1) then with dot product mode signal calculated electric energy in the cycle, can get formula (5).
(5)
Contrast formula (5) and Theoretical Calculation be formula (3) result of calculation as a result, can be as drawing a conclusion:
(1) for the pure sinusoid periodic signal, be parameters such as the electric weight of accurate signal calculated, power, effective value, the signals sampling rate must guarantee that weekly the phase samples more than 3 at least;
(2) signal of the highest for containing
subharmonic; Want the parameters such as electric weight, power, effective value of accurate signal calculated, then the signals sampling rate must guarantee more than per primitive period signal sampling
point;
Notice that conclusion (2) is the inference of conclusion (1), it is a non-adequate condition of necessity.For example, for a sine signal source that contains 21 subharmonic, if will accurately calculate their electric weight, then signal sampling rate must guarantee fundamental signal ripple sampling more than 63 weekly.
Above-mentionedly know that the basic definition formula of having used the definite integral computing is just directly prolonged in the dot product of industry-wide adoption and computing.When sampling number was confirmed, the class of accuracy of the right integrating meter formula was then confirmed thereupon.Numerical analysis in modern times theory is pointed out: if quadrature formula is accurately set up the polynomial expression quadrature of all number of times
; But inaccurate for
order polynomial, claim that then this integral algorithm has
inferior algebraically accuracy.Obviously, the some integration method is a Riemann's ladder integral algorithm that has only 0 rank algebraically accuracy, and it is only accurately set up the integration of direct current measurement (0 order polynomial).Can know with reference to formula (4) and Fig. 1; Standard sine input signal instantaneous power is divided into constant component and positive string section; Because the odd symmetry property of sinusoidal signal; Make in the sinusoidal signal cycle more than 3 that evenly distribution each point sum is zero, the electric flux that therefore integral algorithm of 0 rank algebraically accuracy just can accurate Calculation two sinusoidal signals more than 3.But when input signal no longer be thereby that the standard sine signal is when having destroyed this symmetry; The error of calculation of algorithm will produce thereupon; This asymmetric non-standard sinusoidal input signal is very general in the nonlinear-load metering, causes existing ammeter error in dipping bigger.How could nonlinear-load accurately be measured, especially under asymmetric non-standard sinusoidal signal input condition, accurately to measure nonlinear-load, become present problem demanding prompt solution.
Summary of the invention
The object of the present invention is to provide and a kind ofly can improve the especially nonlinear-load electric energy measuring method of real time measure accuracy of the ammeter accuracy of measuring.
This nonlinear-load electric energy measuring method provided by the invention is at one-period
sampling
point; Gather current instantaneous voltage and momentary current; Calculate current instantaneous power [
]; Based on the principle of high-order newton-Ke Tesi Numerical Integral Formulas current sampling instant electric weight is carried out real-time integration and obtain integral and calculating value (
); And the real-time integration of electric weight is carried out the cycle accumulation obtain cycle electric weight accumulating value (
); Multiply by a coefficient and multiply by the cycle electric flux (
) that a factor obtains nonlinear-load with this accumulating value, the said factor is the ratio of periodic quantity
and sampling number
.
After obtaining above-mentioned electric weight accumulating value (
); When definite current input signal zero passage; Calculate the electric weight (
) of nonlinear-load, said integral and calculating value (
) is fit to following two kinds of formula:
In the above-mentioned formula
Following formula
is crossed the electric weight accumulating value in null cycle for signal; Preserve in be with the instantaneous electric weight of current calculating gained
and after electric weight accumulating value
adds up before this, deposit back
.
carries out real-time integration to current sampling instant electric weight to obtain real-time integral and calculating value;
is periodic quantity;
is sampling number,
be the periodic signal electric energy value after the execution cycle computing.
The inventive method utilizes the Newton-Cotes integral algorithm to have the characteristics of high-order accuracy of measuring grade; Current sampling instant electric weight is carried out real-time integration; And the real-time integration of electric weight carried out the cycle accumulation; This method has the accuracy of high-order numerical integration algebraically; Therefore it has higher algorithm accuracy than industry conventional point integration method, can be adapted to nonlinear-load metering widely, has effectively improved electric energy meter to the accuracy of Modern Nonlinear load electric energy metering and the adaptability that all kinds of nonlinear-loads are measured.The inventor shows that through emulation and experiment show this method provided by the invention is highly effective to improving ammeter accuracy aspect.
Description of drawings
Fig. 1 is the synoptic diagram of prior art metering.
Fig. 2 is the process flow diagram of a kind of embodiment of the inventive method.
Fig. 3 is the process flow diagram of the another kind of embodiment of the inventive method.
The figure as a result that Fig. 4 experimentizes with the inventive method.
Embodiment
Embodiment one: this is an embodiment with 3 rank accuracy metering grades, referring to Fig. 2.
One-period (
) some points of sampling (
) at input signal; Signal time in cycle is divided into
individual interval
; According to current instantaneous voltage and momentary current; [
] maybe can be abbreviated as
to calculate current instantaneous power; Utilize the Simpson integral algorithm that the active power of by stages
such as per 2 is carried out electric weight and calculate, be shown below:
In the formula; Because the complexification integral formula is to adopt two adjacent sectors integrations and the single interval integral formula that substitutes the initial point integration method; For explaining conveniently, adopted the index
different to represent that each computation interval originates in the even number index with former integral formula index
.Simultaneously, please note the difference of following formula and hot just physiology opinion quadrature formula, because algorithm is the continuous integrated value of calculating two intervals, so aforementioned calculation formula coefficient is 2 times of the theoretical calculation formula coefficient.Then, following formula in the summation that adds up of the interval upper integral value of each equal portions, can be got:
(2)
T: the signal period value
N: sampling number, i.e. signals sampling rate in the signal period
Above-mentioned formula (7) is a theoretical calculation formula, can not satisfy the real-time operation requirement of modern real-time digital signal processing system.In order in practical application, to utilize the Xin Fusheng algorithm to calculate in real time, need make suitable modification to following formula.Observe calculating formula (7), and launch and can get:
Formula (8) is real-time computing formula, promptly in practical application, when whenever obtaining a sampled point, then according to the current sampling point index data is accumulated: when index is even number, data are directly accumulated; When index is odd number, then data are taken advantage of 2 accumulations.If
is current sampling instant integral and calculating value; Above-mentioned formulate is become algorithm, and the formula of employing is:
In the formula
In the formula;
is current sampling index or the address offset amount of data in internal memory;
carries out mould 2 complementation computings to index
, and symbol
is the fixed-point data shift left operation.
Above-mentioned real-time integration type is carried out cycle accumulation
; The electric weight that obtains at last is multiplied by coefficient more respectively and multiply by the factor
; Run up to then in total electric weight, computing method are applicable to following formula:
(10)
The cycle electric weight that
obtains when being the input signal zero passage at last in the formula.
The electric weight that formula (5) is calculated is accumulated in total electric weight; Make ammeter accurately measure the electric weight that nonlinear load uses; And in this way; The signal of no matter importing is that the sinusoidal waveform or the asymmetrical waveform of symmetry can both accurately measure, and has effectively improved the accuracy of ammeter to the nonlinear-load metering.
Embodiment two: this is an embodiment with 5 rank accuracy metering grades, referring to Fig. 3.
The one-period of input signal (
) some points of sampling (
), signal time in the cycle is divided into
Individual interval
, according to current instantaneous voltage and momentary current, calculate current instantaneous power [
], utilize quadravalence newton-Ke Tesi integral formula to whenever
4Deng the by stages
Active power carry out electric weight and calculate, be shown below:
In like manner; Note the relation between following formula index
and the aforementioned index
, and the formula coefficient is 4 times of these difference of theoretical calculation formula coefficient.Then, following formula in the summation that adds up of the interval upper integral value of each equal portions, can be got:
: current period electric energy metrical value
In like manner, above-mentioned formula (12) need suitably to be revised the aforementioned calculation formula and be convenient to the engineering real-time operation in order in practical application, to utilize.Observe calculating formula (12), expansion can get:
Formula (13) is the real-time computing formula of 4 rank newton-Ke Tesi; Promptly in practical application; When whenever obtaining a sampled point; Then data are accumulated:, be added to the electric weight totalizer after current sampled data product enlarged 7 times when sampled voltage, when current data indicator memory mould 4 residual values are 0 according to the current sampling point index; When datarams desired value mould 4 residual values are 1, be added to the electric weight totalizer after then current sampled data product value being enlarged 16 times; When datarams desired value mould 4 residual values are 2, be added to the electric weight totalizer after then current sampled data product value being enlarged 6 times; When datarams desired value mould 4 residual values are 3, be added to the electric weight totalizer after then current sampled data product value being enlarged 16 times.So go round and begin again; Back and forth calculate; If
is current sampling instant integral and calculating value; Above-mentioned formulate is become algorithm, and the formula of employing is:
calculates the signal period electric weight of gained in the formula when being the input signal zero passage.This embodiment can obtain 5 grades of integration algebraically class of accuracies, can measure more accurately nonlinear-load.
In order to confirm that the inventive method possesses the higher accuracy of measuring than prior art, the inventor has done following experiment:
If the input voltage and input current signal is respectively
V;
A; Here for realizing the asymmetry of input signal, current signal adopts nonsinusoidal signal.Make
, its waveform can be referring to Fig. 4.Calculating should signal period electric flux theoretical value, can get:
The trapezoidal quadrature of dot product, complexification and the above-mentioned two kinds of embodiments of the present invention of using prior art respectively carry out the cycle active energy to
and calculate under the different mining sampling point.The error of calculation is listed in the table under each algorithm, and is as shown in table 1.
Table 1
Can know in the table, when sampling number increases, the also corresponding increase of the accuracy in computation of each integral formula.But under identical sampling number, high-order Newton-Cotes integral algorithm exceeds 2 one magnitude than the dot product accuracy.It is thus clear that the inventive method can further improve the class of accuracy of integral and calculating for asymmetric load under identical sampling number situation.
Below be the inventive method to be applied to ammeter (the ammeter specifications and models are DTSD341-MA1,3 * 57.7V/100V, 3 * 1.5 (6) A; 20000imp/kWh, 50Hz, 0.1s level (the mark Q/OKRW 013-2010 of enterprise)); Enumerate the meritorious three-phase and four-line Measuring and testing accuracy data of forward, as shown in table 2.
Test result: uncertainty/accuracy is 0.01% (k=2).
Table 2
Can find out from table 2, use ammeter down in the inventive method, its accuracy of measuring is lower than 20% of accuracy index that the highest domestic measurement standard GB/T17215.322-2008 stipulates, from minute having embodied the superiority of the inventive method.
Claims (1)
1. the electric energy gauging method of a nonlinear-load is characterized in that at one-period
Interior sampling
Individual, gather current instantaneous voltage and momentary current, calculate current instantaneous power
P(
k), based on the principle of high-order newton-Ke Tesi Numerical Integral Formulas current sampling instant electric weight is carried out real-time integration and obtain instantaneous integral and calculating value
I k , and the real-time integration of electric weight is carried out the cycle accumulation obtain the electric weight accumulating value
I, multiply by a coefficient and multiply by the cycle electric flux that a factor obtains nonlinear-load with this accumulating value
E, the said factor is a periodic quantity
TWith sampling number
NRatio, obtain the electric weight accumulating value
IAfter, when definite current input signal zero passage, calculate the electric weight of nonlinear-load
E, said integral and calculating value
I k Be fit to following two kinds of formula:
Formula two: based on quadravalence Ke Tesi formula
I k For:
The electric weight of said nonlinear-load
EBe fit to following two kinds of formula:
In the above-mentioned formula
I=I+I
k
Following formula
is crossed the electric weight accumulating value in null cycle for signal; Preserve in be with the instantaneous electric weight of current calculating gained
and after electric weight accumulating value
adds up before this, deposit back
;
In the formula;
is current sampling index or the address offset amount of data in internal memory;
carries out mould 2 complementation computings to index
, and symbol
is the fixed-point data shift left operation;
carries out real-time integration to current sampling instant electric weight to obtain real-time integral and calculating value;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102732345A CN101915872B (en) | 2010-09-07 | 2010-09-07 | Nonlinear-load electric energy measuring method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102732345A CN101915872B (en) | 2010-09-07 | 2010-09-07 | Nonlinear-load electric energy measuring method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101915872A CN101915872A (en) | 2010-12-15 |
CN101915872B true CN101915872B (en) | 2012-08-29 |
Family
ID=43323433
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010102732345A Active CN101915872B (en) | 2010-09-07 | 2010-09-07 | Nonlinear-load electric energy measuring method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101915872B (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103344824B (en) * | 2013-07-08 | 2016-05-04 | 国家电网公司 | The electric flux integration method measuring based on markers |
CN104316894B (en) * | 2014-10-16 | 2017-10-17 | 国家电网公司 | A kind of realistically displayed method of electric energy meter actual motion environment |
CN106033098A (en) * | 2015-03-19 | 2016-10-19 | 上海机电工程研究所 | AC power measuring method and device based on data collection |
CN104989381B (en) * | 2015-06-15 | 2017-12-05 | 西安华瑞网电设备有限公司 | A kind of flexible measurement method and device of beam pumping unit stroke cycle |
CN110687348B (en) * | 2019-10-29 | 2021-04-27 | 四川长虹电器股份有限公司 | Electric energy metering device and electric energy metering method |
CN113257638B (en) * | 2021-06-23 | 2021-11-02 | 江苏能电科技有限公司 | Universal type plastic shell circuit breaker |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1404196A (en) * | 2002-10-16 | 2003-03-19 | 武汉大学 | Harmonic, reactive and negative sequence comprehensive control method and device for complicated power system |
CN101055299A (en) * | 2007-01-31 | 2007-10-17 | 湖南大学 | Impact load electric energy metering method |
-
2010
- 2010-09-07 CN CN2010102732345A patent/CN101915872B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1404196A (en) * | 2002-10-16 | 2003-03-19 | 武汉大学 | Harmonic, reactive and negative sequence comprehensive control method and device for complicated power system |
CN101055299A (en) * | 2007-01-31 | 2007-10-17 | 湖南大学 | Impact load electric energy metering method |
Non-Patent Citations (1)
Title |
---|
粟时平 等.冲击负荷电能计量理论与算法研究.《电测与仪表》.2004,(第10期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN101915872A (en) | 2010-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101915872B (en) | Nonlinear-load electric energy measuring method | |
CN104502700B (en) | Sinusoidal parameter measurement method and system of power signal | |
CN102539892B (en) | Method for obtaining true effective values of AC (alternating current) signals | |
CN101833036B (en) | Method for measuring instantaneous phase of alternating current | |
EP0777125B1 (en) | Vector electricity meters and associated vector electricity metering methods | |
CN104535836B (en) | Electric signal fundamental frequency measuring method and system | |
CN102163844B (en) | Method for detecting state of power system based on phasor measurement unit (PMU) | |
CN102323494A (en) | A Multi-Harmonic Source Identification Method | |
CN107271768A (en) | A kind of least square fitting Dynamic Frequency Measurement method | |
CN102095934B (en) | Measuring method for phase difference of alternating current signals | |
CN101915874A (en) | Harmonic wave detection method based on Fourier transformation | |
CN102435844A (en) | Frequency-independent sinusoidal signal phasor calculation method | |
CN106026086B (en) | The method for dynamic estimation of operating status in a kind of power grid | |
CN105606900A (en) | Single-phase harmonic impedance measuring method based on square wave signals | |
CN102445595A (en) | A real-time measurement method of time-varying power in power system | |
CN102253281A (en) | Method for measuring power grid frequency by resisting harmonic interference and automatically changing phase based on software | |
CN108333426A (en) | Power system frequency measurement method based on fourier algorithm | |
CN106918737A (en) | A kind of method and reactive power detection device for detecting power network dynamic reactive power | |
CN104502775A (en) | Quantitative analysis method for influences on comprehensive energy consumption by electric energy quality | |
CN104569581B (en) | Multi-level set and single-cycle estimation method of power grid frequency measuring | |
CN102280877B (en) | Method for identifying parameter of poor branch of power system through a plurality of measured sections | |
CN101949974A (en) | Field tester of alternating current electric energy meter | |
CN114879125A (en) | Three-phase electric energy meter metering method, three-phase electric energy meter and readable storage medium | |
CN104158173A (en) | Method for calculating voltage-drop line loss rate based on load measurement and calculation | |
CN105067935B (en) | Load tap changer method for parameter estimation based on idle amount of unbalance |
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 | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20101215 Assignee: Shenzhen Wismei Technology Co.,Ltd. Assignor: Wasion Group Co., Ltd. Contract record no.: 2013440020316 Denomination of invention: Nonlinear-load electric energy measuring method Granted publication date: 20120829 License type: Exclusive License Record date: 20131112 |
|
LICC | Enforcement, change and cancellation of record of contracts on the licence for exploitation of a patent or utility model |