CN110531299B - Compressed sensing indirect measurement method for electric energy value of electric energy meter dynamic test signal - Google Patents

Compressed sensing indirect measurement method for electric energy value of electric energy meter dynamic test signal Download PDF

Info

Publication number
CN110531299B
CN110531299B CN201810500638.XA CN201810500638A CN110531299B CN 110531299 B CN110531299 B CN 110531299B CN 201810500638 A CN201810500638 A CN 201810500638A CN 110531299 B CN110531299 B CN 110531299B
Authority
CN
China
Prior art keywords
electric energy
dynamic
signal
energy value
steady
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
Application number
CN201810500638.XA
Other languages
Chinese (zh)
Other versions
CN110531299A (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.)
Beijing University of Chemical Technology
Original Assignee
Beijing University of Chemical Technology
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 Beijing University of Chemical Technology filed Critical Beijing University of Chemical Technology
Priority to CN201810500638.XA priority Critical patent/CN110531299B/en
Publication of CN110531299A publication Critical patent/CN110531299A/en
Application granted granted Critical
Publication of CN110531299B publication Critical patent/CN110531299B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R35/00—Testing or calibrating of apparatus covered by the other groups of this subclass
    • G01R35/04—Testing or calibrating of apparatus covered by the other groups of this subclass of instruments for measuring time integral of power or current

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Remote Monitoring And Control Of Power-Distribution Networks (AREA)

Abstract

The invention discloses a compressed sensing indirect measurement method for an electric energy value of a dynamic test signal of an electric energy meter, which comprises a modal decomposition unit of the dynamic test signal, a steady-state and dynamic signal measurement unit and a dynamic electric energy value calculation unit; the dynamic electric energy value calculation unit comprises an indirect measurement model of the dynamic electric energy value and reading or accumulating of the electric energy value of the standard meter, the indirect measurement model is provided with three input auxiliary measurement parameters, and tracing from the dynamic test signal electric energy value to the steady-state test signal electric energy value is realized according to the indirect measurement model; the compressed sensing indirect measurement method for the electric energy value of the dynamic test signal solves the problem of accurate measurement theory of the electric energy value of the current dynamic test signal and the problem of testing the dynamic error of the electric energy meter under the condition of a pseudo-random dynamic test signal, is simple and convenient to realize, and can be used for comprehensively testing the dynamic error characteristic of the electric energy meter.

Description

Compressed sensing indirect measurement method for electric energy value of electric energy meter dynamic test signal
1. Field of the invention
The invention relates to a compressed sensing indirect measurement method for a dynamic test signal reference electric energy value, in particular to an indirect method for measuring a standard reference electric energy value of an electric energy meter dynamic error test signal, which can be used for testing the dynamic error characteristic of an intelligent electric energy meter.
2. Background of the invention
With the continuous development of smart grid construction, from the perspective of power supply and load, the smart grid presents two new characteristics: firstly, due to popularization and use of new power sources such as wind energy, solar energy, tidal energy and the like, the output power of the novel power source of the power grid presents unsteady state characteristics; secondly, heavy industrial users are increasingly loaded with high-power dynamics, such as: electric load power of steel-making electric arc furnace, electrolytic aluminum, high-speed rail electric locomotive, etc. shows frequent random dynamic change characteristics. The new characteristics of the smart grid make the metering of the electric energy generated and consumed and the dynamic error testing of the smart electric energy meter face new challenges. At present, the production, test and detection of intelligent electric energy meters, electronic electric energy meters and induction electric energy meters are all carried out under a steady-state condition, which leads to the fact that the dynamic error characteristics of the electric energy meters cannot be tested and evaluated. The method for accurately measuring and tracing the reference electric energy value of the dynamic test signal of the electric energy meter is complex, and a perfect reference electric energy value measurement model and an accurate measurement and tracing method of the dynamic test signal are not formed at present, so that the progress of the dynamic error characteristic test technology of the electric energy meter is indirectly influenced.
3. Summary of the invention
In order to avoid the defects of the existing method, the invention can simply and conveniently realize the tracing of the dynamic test signal reference electric energy value, measure and determine the dynamic electric energy value of the intelligent electric energy meter, and play an important role in solving the problem of testing the dynamic error characteristic of the electric energy meter under the high-power dynamic load.
The invention aims to solve the problems of measuring the reference electric energy value of a dynamic test signal of an intelligent electric energy meter and tracing the standard reference electric energy value, and provides a compressed sensing indirect measurement model and method of the standard reference electric energy value of the dynamic test signal, so that the accuracy of a dynamic error test method of the electric energy meter is improved. The indirect measurement model is a steady-state electric energy which is easy to measure and a positive characteristic value cumulant M+And negative eigenvalue cumulant M-And as an input auxiliary measurement parameter, determining an indirect measurement model of the measured dynamic reference electric energy value, thereby providing a tracing method of the auxiliary steady-state test signal electric energy value from the dynamic reference electric energy value to the indirect measurement model.
The invention relates to a compressed sensing indirect measurement method for an electric energy value of a dynamic test signal of an electric energy meter, which adopts a dynamic test voltage and current signal (1), a dynamic test signal modal decomposition preprocessing unit (2), a steady and dynamic signal measurement unit (3) and a dynamic electric energy value calculation unit (4) to complete the measurement and tracing of the reference electric energy value of the dynamic test signal. The dynamic test voltage and current signal (1) is a three-phase/single-phase steady-state test voltage signal uk(t) and the steady-state current signal ik(t) three-phase/single-phase dynamic test voltage signal
Figure GDA0002773420280000012
Dynamic test current signal modulated with two/three amplitudes
Figure GDA0002773420280000013
Dynamic test voltage and current signal (1) inputThe dynamic test signal modal decomposition and processing unit (2) obtains two intrinsic modal signals of a steady-state power signal and a dynamic current; the two intrinsic mode signals are respectively measured by a steady-state and dynamic signal measuring unit (3), and the steady-state power signal is measured by a standard electric energy meter to obtain the electric energy value ES0For a dynamic current eigenmode signal, a positive eigenvalue cumulant M of the signal is obtained by adopting a compressed sensing measurement matrix to measure+And negative eigenvalue cumulative quantity M-(ii) a Steady state electric energy value ES0Positive eigenvalue cumulant M+And negative eigenvalue cumulant M-As three auxiliary measurement parameters, the three parameters are input into a dynamic electric energy value indirect measurement model in a dynamic electric energy value calculation unit (4), and a reference electric energy value E of a dynamic test signal is calculated by adopting a dynamic measurement algorithmX0The indirect measurement model realizes the reference electric energy value E of the dynamic test signalX0To steady state test signal electric energy value ES0Tracing to the source by using a reference electric energy value EX0And the electric energy value E output by the electric energy meter to be testedZAnd calculating the dynamic error of the electric energy meter.
The steady state test voltage signal uk(t) and a steady-state test current signal ik(t) which are respectively expressed as:
uk(t)＝Usin(Ω1t+φk)
Figure GDA0002773420280000011
in the formula, omega1＝2πf1Is the angular frequency, f1At power frequency, T1＝1/f1Power frequency period of steady state signal, phikFor the three-phase steady-state voltage signal phase,
Figure GDA0002773420280000021
for three-phase steady-state current signal phases, k ═ a, b, c represent phase a, phase b, and phase c signals, respectively; from signal uk(t) and signal ik(t) multiplying to obtain a steady state power signal
Figure GDA0002773420280000022
Figure GDA0002773420280000023
The dynamic test voltage signal
Figure GDA0002773420280000024
And two-amplitude/three-amplitude dynamic test current signal
Figure GDA0002773420280000025
The expressions are respectively:
Figure GDA0002773420280000026
Figure GDA0002773420280000027
Figure GDA0002773420280000028
Figure GDA0002773420280000029
from signals
Figure GDA00027734202800000210
Sum signal
Figure GDA00027734202800000211
Multiplying to obtain dynamic test power signal
Figure GDA00027734202800000212
Figure GDA00027734202800000213
Figure GDA00027734202800000214
Figure GDA00027734202800000215
Figure GDA00027734202800000216
In the formula, AkL＝[αk0,αk1,…,αk(L-1)]，αk(L+n)＝αknFor dynamically testing the current
Figure GDA00027734202800000217
Characteristic vector of (B)L＝[β0,β1,…,βL-1]1×L，βk(L+n)＝βknFor dynamically testing voltage
Figure GDA00027734202800000218
L ∈ [2, ∞) is a bounded positive integer, T ═ LT1For dynamic test signal periods, G (t-nT) is a window function, Gn(t)＝[g(t-T)，g(t-2T)，…g(t-nT)]Is a window function vector.
The modal decomposition and processing unit of the dynamic test signal is used for decomposing and processing the dynamic test voltage
Figure GDA00027734202800000219
And dynamic test current
Figure GDA00027734202800000220
Multiplying to obtain dynamic test power signal
Figure GDA00027734202800000221
Simultaneously dynamically testing current signals
Figure GDA00027734202800000222
And
Figure GDA00027734202800000223
decomposing to obtain dynamic current characteristic signal
Figure GDA00027734202800000224
And steady state
Figure GDA00027734202800000225
Using steady-state test voltage uk(t) and the steady-state test current ik(t) multiplying to obtain steady state power
Figure GDA00027734202800000226
By using
Figure GDA00027734202800000227
Instead of the former
Figure GDA00027734202800000228
Sending into a steady-state and dynamic signal measuring unit, namely, keeping steady-state power
Figure GDA00027734202800000229
The stability of the dynamic test signal is convenient for realizing the tracing of the electric energy value of the dynamic test signal.
The steady-state and dynamic signal measuring unit is composed of a standard electric energy meter and a compressed sensing measuring matrix: at a measurement time L' T1Measuring steady-state power by internally adopting standard electric energy meter
Figure GDA00027734202800000230
Accumulated steady state test signal electric energy value ES0(ii) a At the same measuring time L' T1Internal, to dynamic current signal
Figure GDA00027734202800000231
Or dynamic current characteristic signal
Figure GDA00027734202800000232
Measuring matrix phi using compressed sensing matched filterBMeasurement calculation
Figure GDA00027734202800000233
Positive eigenvalue cumulative quantity M of+And negative eigenvalue cumulative quantity M-(ii) a Compressed sensing matched filter measurement matrix phiBComprises the following steps:
Figure GDA00027734202800000234
in the formula (I), the compound is shown in the specification,
Figure GDA0002773420280000031
positive eigenvalue cumulative quantity M+And negative eigenvalue cumulative quantity M-The calculation algorithm is as follows:
Figure GDA0002773420280000032
in the formula (I), the compound is shown in the specification,
Figure GDA0002773420280000033
respectively taking a binary sequence or a ternary sequence.
In the dynamic electric energy value calculation unit, the input auxiliary measurement parameter of the dynamic electric energy value indirect measurement model is a steady-state test signal electric energy value ES0Positive eigenvalue cumulant M+And negative eigenvalue cumulative quantity M-The output is a reference electric energy value E of the measured dynamic test signalX0(ii) a By reading the electric energy value of the standard electric energy meter or accumulating the standard electric energy meter at the measuring time L' T1The number of the internally output electric energy pulses is calculated to obtain a steady-state electric energy value ES0(ii) a Wherein, the indirect measurement model of the dynamic electric energy value is as follows:
Figure GDA0002773420280000034
in the formula, EX0For dynamic test signalsL (|) is a likelihood function,
Figure GDA0002773420280000035
ES0the steady state test signal electric energy value;
the reference electric energy value E of the dynamic test signal is calculatedX0The dynamic measurement algorithm is as follows:
Figure GDA0002773420280000036
this equation gives the reference electrical energy value E from the dynamic test signalX0To steady state test signal electric energy value ES0The source tracing relationship.
The dynamic load current signal characteristic vector AkL′＝[αk0,αk1,…,αk(L′-1)]，αkn＝αk(L+n)Can be a binary sequence of "0" and "+ 1", or a ternary sequence of "+ 1", "0" and "-1", and can be a deterministic sequence or a pseudo-random sequence.
The invention has the beneficial effects that:
the method can accurately measure and obtain the reference electric energy value of the dynamic test signal by utilizing the indirect measurement model of the dynamic electric energy value under the condition of the dynamic test signal of the electric energy meter according to the test application requirement of the dynamic error characteristic of the electric energy meter, effectively solves the measurement and source tracing problems of the electric energy value of the determined and pseudo-random dynamic test signal, and obviously improves the accuracy of the electric energy meter dynamic error test method.
4. Drawings and description of the drawings
FIG. 1 is a compressed sensing indirect measurement method for the electric energy value of the dynamic test signal of the electric energy meter according to the invention.
5. Detailed description of the preferred embodiments
The first embodiment is as follows: the present embodiment will be described with reference to fig. 1. The compressed sensing indirect measurement method for the electric energy value of the dynamic test signal of the electric energy meter in the embodiment comprises the following steps: the dynamic test signal modal decomposition unit (2), the steady-state and dynamic signal measurement unit (3) and the dynamic electric energy value calculation unit (4). Wherein the input signals are: three-phase/single-phase steady-state voltage and current signals (1) generated using a steady-state power source, and two-amplitude/three-amplitude dynamic test current signals (1) generated using a signal generation unit or device to modulate the steady-state voltage and current signals (1).
According to the compressed sensing indirect measurement method for the electric energy value of the dynamic test signal of the electric energy meter, according to the measurement requirement of the dynamic error characteristic of the electric energy meter; aiming at three-phase/single-phase dynamic test current signals and three-phase/single-phase steady-state test voltage and current signals, a standard electric energy meter is adopted to measure steady-state power
Figure GDA0002773420280000037
Cumulative steady state test signal power magnitude ES0While measuring matrix phi by using compressed sensing matched filterBMeasurement calculation
Figure GDA0002773420280000038
Positive eigenvalue cumulative quantity M of+And negative eigenvalue cumulative quantity M-And measuring to obtain the reference electric energy value of the measured dynamic test signal by adopting an indirect measurement model of the dynamic electric energy value and a dynamic measurement algorithm of the electric energy value of the dynamic test signal.
The second embodiment is as follows: in this embodiment, the first embodiment will be described further, and the three-phase dynamic test voltage signal in this embodiment will be described
Figure GDA0002773420280000041
Three-phase dynamic test current signal
Figure GDA0002773420280000042
Respectively as follows:
Figure GDA0002773420280000043
wherein the content of the first and second substances,
Figure GDA0002773420280000044
Figure GDA0002773420280000045
wherein the content of the first and second substances,
Figure GDA0002773420280000046
dynamic test power signal
Figure GDA0002773420280000047
Is composed of signals
Figure GDA0002773420280000048
Sum signal
Figure GDA0002773420280000049
Multiplication results in:
Figure GDA00027734202800000410
Figure GDA00027734202800000411
Figure GDA00027734202800000412
Figure GDA00027734202800000413
in the formula uk(t)＝Usin(Ω1t+φk)，
Figure GDA00027734202800000414
Ω1＝2πf1Is the angular frequency, f1At power frequency, T1＝1/f1Power frequency period of steady state signal, phikFor the three-phase steady-state voltage signal phase,
Figure GDA00027734202800000415
the three-phase steady-state current signal phase is defined as k, a, b and c respectively represent a phase, a phase and a phase; a. thekL＝[αk0,αk1,…,αk(L-1)]，αk(L+n)＝αknFor dynamically testing the current
Figure GDA00027734202800000416
Characteristic vector of (B)L＝[β0,β1,…,βL-1]1×L，βk(L+n)＝βknFor dynamically testing voltage
Figure GDA00027734202800000417
L ∈ [2, ∞) is a bounded positive integer, T ═ LT1For dynamic test signal periods, G (t-nT) is a window function, Gn(t)＝[g(t-T)，g(t-2T)，g(t-nT)]Is a window function vector;
Figure GDA00027734202800000418
for steady state power signals:
Figure GDA00027734202800000419
the third concrete implementation mode: the present embodiment is further described with respect to the first and second embodiments, and the dynamic test signal modal decomposition and processing unit of the present embodiment is configured to perform a dynamic test on a dynamic test power signal
Figure GDA00027734202800000420
And dynamic test current signal
Figure GDA00027734202800000421
Processing and decomposing to obtain dynamic current characteristic signal
Figure GDA00027734202800000422
And steady state
Figure GDA00027734202800000423
Dynamic current signature
Figure GDA00027734202800000424
Comprises the following steps:
Figure GDA00027734202800000425
Figure GDA00027734202800000426
and
Figure GDA00027734202800000427
having the same feature vector AkL＝[αk0,αk1,…,αk(L-1)]，αk(L+n)＝αkn。
Using steady-state test voltage uk(t) and the steady-state test current ik(t) multiplying to obtain steady state power
Figure GDA00027734202800000428
Instead of the former
Figure GDA00027734202800000429
Ensuring dynamic test power signals
Figure GDA00027734202800000430
To steady state power signal
Figure GDA00027734202800000431
The source tracing is carried out on the original data,
Figure GDA00027734202800000432
and
Figure GDA00027734202800000433
sending into a steady state and dynamic signal measuring unit for measurement
Figure GDA00027734202800000434
Cumulative steady state test signal power magnitude ES0And
Figure GDA00027734202800000435
positive eigenvalue cumulative quantity M of+And negative eigenvalue cumulant M-。
The fourth concrete implementation mode: the present embodiment is further described with respect to the first, second and third embodiments, and the steady-state and dynamic signal measuring unit of the present embodiment measures the time L' T1Measuring steady-state power by using standard electric energy meter in internal (L' > L)
Figure GDA0002773420280000051
Cumulative steady state test signal power magnitude ES0(ii) a At the same measuring time L' T1Internal, using a compressed sensing matched filter to measure the matrix phiBMeasurement calculation
Figure GDA0002773420280000052
Or
Figure GDA0002773420280000053
Positive eigenvalue cumulative quantity M of+And negative eigenvalue cumulative quantity M-；
The compressed sensing matched filter measurement matrix phiBComprises the following steps:
Figure GDA0002773420280000054
in the formula (I), the compound is shown in the specification,
Figure GDA0002773420280000055
positive eigenvalue cumulative quantity M+And negative eigenvalue cumulative quantity M-The measurement calculation algorithm is as follows:
Figure GDA0002773420280000056
in the formula (I), the compound is shown in the specification,
Figure GDA0002773420280000057
respectively taking a binary sequence or a ternary sequence.
The fifth concrete implementation mode: the present embodiment is further described with respect to the first, second, third, or fourth embodiment, and the dynamic electric energy value calculation unit according to the present embodiment includes two parts, namely, a dynamic electric energy value indirect measurement model and a read/accumulated electric energy value. The dynamic electric energy value indirect measurement model adopts a steady-state test signal electric energy value ES0Positive eigenvalue cumulant M+And negative eigenvalue cumulative quantity M-Reference electric energy value E as input auxiliary measuring parameter through dynamic testing signalX0The reference electric energy value E of the measured dynamic test signal is obtained by calculationX0(ii) a By reading the electric energy value of the standard electric energy meter or accumulating the standard electric energy meter at the measuring time L' T1The number of the internally output electric energy pulses is calculated to obtain the electric energy value E of the steady-state test signalS0；
Wherein, the indirect measurement model of the dynamic electric energy value is as follows:
Figure GDA0002773420280000058
in the formula, EX0For the reference power value of the dynamic test signal, L (-) is a likelihood function,
Figure GDA0002773420280000059
ES0the steady state test signal electric energy value;
reference electric energy value E of the dynamic test signalX0The dynamic measurement algorithm is as follows:
Figure GDA00027734202800000510
in the formula, EX0To ES0The calculation relationship gives the principle of tracing from the electric energy value of the dynamic test signal to the electric energy value of the steady test signal。
It will be evident to those skilled in the art that the invention is not limited to the details of the foregoing illustrative embodiments, and that the present invention may be embodied in other specific forms without departing from the spirit or essential attributes thereof. The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. Any reference sign in a claim should not be construed as limiting the claim concerned.
Furthermore, it should be understood that although the present description refers to embodiments, not every embodiment may contain only a single embodiment, and such description is for clarity only, and those skilled in the art should integrate the description, and the embodiments may be combined as appropriate to form other embodiments understood by those skilled in the art.

Claims (4)

1. A compressed sensing indirect measurement method for an electric energy value of an electric energy meter dynamic test signal is characterized in that an input dynamic voltage and current test signal (1) is a three-phase/single-phase steady-state test voltage signal uk(t) and the steady-state current signal ik(t), three-phase/single-phase dynamic test voltage signal
Figure FDA0002773420270000011
Dynamic test current signal modulated with two/three amplitudes
Figure FDA0002773420270000012
Respectively called as a steady-state periodic eigenmode signal and a dynamic eigenmode signal modulated by two amplitudes/three amplitudes; the input dynamic test signal obtains two modes of a steady-state power signal and a dynamic current intrinsic mode signal in the dynamic test signal mode decomposition and processing unit (2); the two modal signals are measured by the steady state and dynamic signals respectivelyThe unit (3) measures the steady-state power signal by adopting a standard electric energy meter to obtain the electric energy value ES0For a dynamic current eigenmode signal, a positive eigenvalue cumulant M of the signal is obtained by adopting a compressed sensing measurement matrix to measure+And negative eigenvalue cumulative quantity M-(ii) a Steady state test signal electric energy value ES0Positive eigenvalue cumulant M+And negative eigenvalue cumulant M-As three auxiliary measurement parameters, the three parameters are input into a dynamic electric energy value indirect measurement model in a dynamic electric energy value calculation unit (4), and a reference electric energy value E of a dynamic test signal is calculated by adopting a dynamic measurement algorithmX0The indirect measurement model realizes the reference electric energy value E of the dynamic test signalX0To steady state test signal electric energy value ES0Tracing to the source, using reference electric energy value E of the dynamic test signalX0And the electric energy value E output by the electric energy meter to be testedZCalculating the dynamic error of the electric energy meter;
wherein the steady-state periodic eigenmode signal is a steady-state test voltage signal uk(t) and a steady-state test current signal ik(t) which are respectively expressed as:
uk(t)＝Usin(Ω1t+φk)
Figure FDA0002773420270000013
in the formula, omega1＝2πf1Is the angular frequency, f1At power frequency, T1＝1/f1Power frequency period of steady state signal, phikFor the three-phase steady-state voltage signal phase,
Figure FDA0002773420270000014
the three-phase steady-state current signal phase is defined as k, a, b and c respectively represent a phase, a phase and a phase; from signal uk(t) and signal ik(t) multiplying to obtain a steady state power signal
Figure FDA0002773420270000015
Figure FDA0002773420270000016
The dynamic eigenmode signal input with the two-amplitude/three-amplitude modulation is dynamic test voltage
Figure FDA0002773420270000017
And dynamic test current signal
Figure FDA0002773420270000018
The expressions are respectively:
Figure FDA0002773420270000019
Figure FDA00027734202700000110
Figure FDA00027734202700000111
Figure FDA00027734202700000112
from signals
Figure FDA00027734202700000113
Sum signal
Figure FDA00027734202700000114
Multiplying to obtain dynamic test power signal
Figure FDA00027734202700000115
Figure FDA00027734202700000116
Figure FDA00027734202700000117
Figure FDA00027734202700000118
Figure FDA00027734202700000119
In the formula, AkL＝[αk0,αk1,…,αk(L-1)]，αk(L+n)＝αknFor dynamically testing the current
Figure FDA00027734202700000120
Characteristic vector of (B)L＝[β0,β1,…,βL-1]1×L，βk(L+n)＝βknFor dynamically testing voltage
Figure FDA00027734202700000121
L ∈ [2, ∞) is a bounded positive integer, T ═ LT1For dynamic test signal periods, G (t-nT) is a window function, Gn(t)＝[g(t-T)，g(t-2T)，…g(t-nT)]Is a window function vector;
the input auxiliary measurement parameter of the dynamic electric energy value indirect measurement model in the dynamic electric energy value calculation unit (4) is a steady-state test signal electric energy value ES0Positive eigenvalue cumulant M+And negative eigenvalue cumulative quantity M-The output is a reference electric energy value E of the measured dynamic test signalX0(ii) a By reading the electric energy value of the standard electric energy meter or accumulating the standard electric energy meter at the measuring time L' T1The number of the internally output electric energy pulses is calculated to obtain the electric energy value E of the steady-state test signalS0(ii) a Calculating reference electric energy value E of dynamic test signal by adopting dynamic measurement algorithmX0；
Wherein, the indirect measurement model of the dynamic electric energy value is as follows:
Figure FDA0002773420270000021
in the formula, EX0Reference electric energy value for dynamic test signal, ES0For steady state test signal power magnitude, L (-) is a likelihood function,
Figure FDA0002773420270000022
the reference electric energy value E of the dynamic test signal is calculatedX0The dynamic measurement algorithm is as follows:
Figure FDA0002773420270000023
this formula gives EX0To ES0The source tracing relationship.
2. The method of claim 1, wherein the modal decomposition and processing unit of the dynamic test signal transforms the dynamic test voltage into a dynamic test voltage
Figure FDA0002773420270000024
And dynamic test current
Figure FDA0002773420270000025
Multiplying to obtain dynamic test power signal
Figure FDA0002773420270000026
Simultaneously dynamically testing current signals
Figure FDA0002773420270000027
And
Figure FDA0002773420270000028
decomposing to obtain dynamic current characteristic signal
Figure FDA0002773420270000029
And steady state
Figure FDA00027734202700000210
Using steady-state test voltage uk(t) and the steady-state test current ik(t) multiplying to obtain steady state power
Figure FDA00027734202700000211
By using
Figure FDA00027734202700000212
Instead of the former
Figure FDA00027734202700000213
Sending into a steady state and dynamic signal measurement unit to reserve steady state power
Figure FDA00027734202700000214
The stability and the characteristic of being convenient for tracing the electric energy value are realized; dynamic current signature
Figure FDA00027734202700000215
Comprises the following steps:
Figure FDA00027734202700000216
3. the method of claim 1, wherein the steady-state and dynamic signal measurements are performed by a compressed sensing indirect measurement method of the magnitude of the dynamic test signal power of the electric energy meterThe unit cell is a unit cell, which is composed of a plurality of unit cells,
Figure FDA00027734202700000217
and
Figure FDA00027734202700000218
sending into a steady state and dynamic signal measuring unit for measuring time L' T1Measuring steady-state power by standard electric energy meter
Figure FDA00027734202700000219
Cumulative steady state test signal power magnitude ES0(ii) a At the same measuring time L' T1Internal, to dynamic current signal
Figure FDA00027734202700000220
Measuring matrix phi using compressed sensing matched filterBMeasurement calculation
Figure FDA00027734202700000221
Or
Figure FDA00027734202700000222
Positive eigenvalue cumulative quantity M of+And negative eigenvalue cumulative quantity M-；
Wherein the compressed sensing matched filter measures the matrix phiBComprises the following steps:
Figure FDA00027734202700000223
in the formula (I), the compound is shown in the specification,
Figure FDA00027734202700000224
positive eigenvalue cumulative quantity M+And negative eigenvalue cumulative quantity M-The calculation is as follows:
Figure FDA00027734202700000225
in the formula (I), the compound is shown in the specification,
Figure FDA0002773420270000031
respectively taking a binary sequence or a ternary sequence.
4. The method according to any one of claims 1 to 3, wherein the method is applied to the dynamic load current signal characteristic vector AkL′＝[αk0,αk1,…,αk(L′-1)]，αkn＝αk(L+n)Can be a binary sequence of "0" and "+ 1", or a ternary sequence of "+ 1", "0" and "-1", and can be a deterministic sequence or a pseudo-random sequence.
CN201810500638.XA 2018-05-23 2018-05-23 Compressed sensing indirect measurement method for electric energy value of electric energy meter dynamic test signal Active CN110531299B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810500638.XA CN110531299B (en) 2018-05-23 2018-05-23 Compressed sensing indirect measurement method for electric energy value of electric energy meter dynamic test signal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810500638.XA CN110531299B (en) 2018-05-23 2018-05-23 Compressed sensing indirect measurement method for electric energy value of electric energy meter dynamic test signal

Publications (2)

Publication Number Publication Date
CN110531299A CN110531299A (en) 2019-12-03
CN110531299B true CN110531299B (en) 2021-01-01

Family

ID=68656761

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810500638.XA Active CN110531299B (en) 2018-05-23 2018-05-23 Compressed sensing indirect measurement method for electric energy value of electric energy meter dynamic test signal

Country Status (1)

Country Link
CN (1) CN110531299B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103091659A (en) * 2013-01-22 2013-05-08 黑龙江省电力科学研究院 Electric energy meter dynamic measurement cycle power sequence proportional allocation calibration method
CN103630871A (en) * 2013-12-09 2014-03-12 国家电网公司 Testing device for dynamic characteristics of electronic energy meter
CN103809148A (en) * 2012-11-12 2014-05-21 深圳供电局有限公司 Method and system for measuring dynamic characteristics of electric energy meter
CN104345294A (en) * 2014-11-04 2015-02-11 国家电网公司 Electric energy meter dynamic performance test device
CN104569900A (en) * 2013-10-14 2015-04-29 北京化工大学 Dynamic error testing method and device for intelligent ammeter
CN106501755A (en) * 2016-10-18 2017-03-15 国网山东省电力公司烟台供电公司 A kind of intelligent electric energy meter dynamic error measuring method based on dynamic load model

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8803546B2 (en) * 2011-11-27 2014-08-12 Mediatek Singapore Pte. Ltd. Method for performing power detection, and associated apparatus

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103809148A (en) * 2012-11-12 2014-05-21 深圳供电局有限公司 Method and system for measuring dynamic characteristics of electric energy meter
CN103091659A (en) * 2013-01-22 2013-05-08 黑龙江省电力科学研究院 Electric energy meter dynamic measurement cycle power sequence proportional allocation calibration method
CN104569900A (en) * 2013-10-14 2015-04-29 北京化工大学 Dynamic error testing method and device for intelligent ammeter
CN103630871A (en) * 2013-12-09 2014-03-12 国家电网公司 Testing device for dynamic characteristics of electronic energy meter
CN104345294A (en) * 2014-11-04 2015-02-11 国家电网公司 Electric energy meter dynamic performance test device
CN106501755A (en) * 2016-10-18 2017-03-15 国网山东省电力公司烟台供电公司 A kind of intelligent electric energy meter dynamic error measuring method based on dynamic load model

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
dynamical waveforms and the dynamical source for electricity meter dynamical experiment;Hongtao Huang et.al;《2016 conference on precision electromagnetic measurements》;20160811;正文第1-2页 *
压缩感知伪随机动态功率信号的电能测量方法;王学伟 等;《电子与信息学报》;20170331;第39卷(第3期);第640-646页 *
电能表动态误差特性实验研究;王学伟 等;《电测与仪表》;20131231;第50卷(第12期);第1-4页 *

Also Published As

Publication number Publication date
CN110531299A (en) 2019-12-03

Similar Documents

Publication Publication Date Title
CN102193031B (en) AC impedance measuring device
CN103513211B (en) Ac impedance measurement measuring device
CN109752584B (en) Method for measuring effective value of periodic signal
CN106569164A (en) Method and system for synchronization testing of electric quantity of double-core electric energy meter
CN108318852B (en) Square wave influence test method for intelligent electric energy meter
Perez et al. Sensitivity-based interval PDE observer for battery SOC estimation
WO2017113754A1 (en) Method for detecting dynamic reactive power of grid and reactive power detector
CN111142032B (en) Method, device and equipment for determining battery electric quantity and storage medium
CN102288817B (en) A kind of charging pile consumption detection system and detection method based on virtual load calibrating
JP2014052186A (en) Method for estimating capacity maintenance rate of secondary battery
MT Ibraheem Al-Naib et al. A Cost-Effective Method for Power Factor‎ Metering Systems‎
CN114200381A (en) System and method for detecting reliability of intelligent electric meter
CN109829652A (en) A kind of long time scale dynamic harmonic divisions of responsibility method
CN110531299B (en) Compressed sensing indirect measurement method for electric energy value of electric energy meter dynamic test signal
KR101386344B1 (en) The real-time monitoring device of fourier transform electrochemical impedance and method.
CN103293424A (en) Data center power usage efficiency (PUE) determination method
KR20110063184A (en) Apparatus for measuring generation power and load power for power generation system
CN109856550A (en) The evaluation method and system of the SOC of flow battery
Zheng et al. Estimation of harmonic impedance and harmonic contribution with harmonic complex power in the absence of harmonic phase angle
CN202393834U (en) Precise ohmic internal resistance measurer for battery
CN105158555A (en) Method for measuring and verifying saved electric power and electric quantity in boosting operation of power distribution line
CN105911346A (en) Bi-directional intelligent electric meter respectively used for metering fundamental wave and harmonic wave and metering method thereof
KR20150015200A (en) Device for measuring frequency of power based on zero-crossing and method for measuring frequency of power based on zero-crossing
Singh et al. Regression based approach for measurement of current in single-phase smart energy meter
CN204065387U (en) A kind of synchronous demodulator and comprise the power standard source of this synchronous demodulator

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