CN106548021B - A kind of electronic mutual inductor error state appraisal procedure based on higher-dimension random matrix - Google Patents
A kind of electronic mutual inductor error state appraisal procedure based on higher-dimension random matrix Download PDFInfo
- Publication number
- CN106548021B CN106548021B CN201610953170.0A CN201610953170A CN106548021B CN 106548021 B CN106548021 B CN 106548021B CN 201610953170 A CN201610953170 A CN 201610953170A CN 106548021 B CN106548021 B CN 106548021B
- Authority
- CN
- China
- Prior art keywords
- matrix
- mutual inductor
- error state
- electronic mutual
- appraisal procedure
- 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
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/02—Testing or calibrating of apparatus covered by the other groups of this subclass of auxiliary devices, e.g. of instrument transformers according to prescribed transformation ratio, phase angle, or wattage rating
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
- Testing Electric Properties And Detecting Electric Faults (AREA)
Abstract
The invention discloses a kind of electronic mutual inductor error state appraisal procedure based on higher-dimension random matrix, it is intended to independent of standard utensil or physical model, exported only according to electronic mutual inductor, whether have anomalous variation, so as to carry out error state assessment if the statistical distribution of the measurement error of electronic mutual inductor is judged based on higher-dimension random matrix.Detailed process includes：By the output data of electronic mutual inductor, based on time slip-window, random matrix is built；Based on Kalman filter, random matrix is extended；Ask for non-Hermite Matrix；Evaluating matrix is asked for；The evaluation index for covering sample central moment and sample moment of the orign is calculated, and electronic mutual inductor error state is assessed.The present invention is realized on the premise of independent of standard and physical model, and the error state of electronic mutual inductor is effectively assessed, and assessment result is accurate and visual, and appraisal procedure has universality, validity and easy implementation.
Description
Technical field
The invention belongs to equipment for power transmission and distribution state estimation and fault diagnosis field, and higher-dimension is based on more particularly, to one kind
The electronic mutual inductor error state appraisal procedure of random matrix.
Background technology
Electric power mutual-inductor is protecting electrical power system, control, metering accurately and reliably electric current, voltage signal, its error are provided
And reliability and the safe and reliable operation of power system it is closely bound up.With the development of China's power industry, power system capacity
Constantly increase, voltage class also constantly raise, and traditional electromagnetic transformer exposes a series of shortcomings：(1) volume and weight
It is larger；(2) insulation system is complicated；(3) phenomena such as intrinsic magnetic saturation, ferromagnetic resonance；(4) frequency band is narrow, transient characterisitics are poor；(5)
100V or 5A analog output does not meet digitized growth requirement of signal of intelligent substation.
Electronic mutual inductor becomes a kind of preferable due to outstanding advantages of sensing principle is simple, measurement quality is superior
The substitute of electromagnetic transformer, is widely applied in intelligent substation, possesses and is come into operation in batches in engineering
Ability.However, electronic mutual inductor still has following prominent question in the process of running：(1) part electronic type is mutual
There is drift to sensor in error in the process of running, have impact on protection, control, metering and the use of measurement function；(2) electronic type
The error state of transformer is in the process of running unknown for user, interior only during maintenance, utilizes standard utensil pair
Electronic mutual inductor is verified offline, and it is online in the case of error state assess and be difficult to.It is mutual in order to improve electronic type
The operation management of sensor is horizontal, electronic mutual inductor error state is responded in time, in the constraint independent of standard utensil
Under the conditions of, the appraisal procedure of research electronic mutual inductor error state becomes a kind of active demand.
Prior art includes the electronic mutual inductor error state appraisal procedure based on signal transacting, and its key is to electronics
The output abnormality signal characteristic of formula transformer is classified, extracted and identified.When the measurement error generation of electronic mutual inductor is different
Often during change, its output signal has some forms of expression, as deviation increases suddenly, signal suddenly disappears, output signal distorts
Deng, therefore measurement error anomalous variation can be identified by analyzing the output of electronic mutual inductor.But this method is present
Following deficiency：(1) comparison " standard " of the appraisal procedure is actually the experience of people, that is, thinks the output of electronic mutual inductor
Signal should not be mutated under normal circumstances, and this condition is not to set up always in actual applications；(2) appraisal procedure
The influence that power network once disturbs can not be excluded, i.e. the premise of accurate judgement is that power network primary side does not break down；(3) the assessment side
Method is more effective for electronic mutual inductor error mutability failure, but the judgement for error gradual failure is then more weary
Power.
Prior art also includes the electronic mutual inductor error state appraisal procedure based on mathematics physics model, and its key is
Equipment is modeled, then according to model solution sensor output value and is used as theoretical true value, by theoretical true value and reality
Output valve, which is compared, realizes that error state is assessed.But this method has the following disadvantages：(1) assessment of error state is high
Degree relies on the accuracy for establishing model, and accurate mathematics physics model is generally difficult to build；(2) the hypothesis bar in modeling process
Part and simplified condition have impact on the accuracy and reliability of assessment result, and the accuracy requirement of electronic mutual inductor error evaluation
Gap be present.
Prior art also includes the online Compare System that Guangdong DianKeYuan and Jiangsu DianKeYuan are established, to electronic mutual inductor
Error state assessed.This method is by setting up a set of identical essence in the same interval where transformer substation electronic transducer
The conventional electromagnetic transformer of grade is spent, and electronic transducer calibration instrument is installed in control room, so as to realize that electronic type is mutual
The long-term monitoring of sensor error state and assessment.But this method has the following disadvantages：(1) need to add in transformer station to set
Standby, this does not possess universality in engineer applied；(2) premise of this method application needs the error of conventional electromagnetic transformer to expire
Sufficient accuracy requirement, but its error state is also unknown in the process of running.
Higher-dimension Random Matrices Theory is a kind of effective big data analysis means, in physics, radio communication, metering finance
Economics, field of biology are applied, and on the other hand, the output parameter of electronic mutual inductor has the spy of steady random distribution
Sign, possesses the set up the condition of higher-dimension random matrix.Although the Asymptotic Behavior For Some in Random Matrices Theory requires the random square of higher-dimension
The dimension convergence of battle array is infinite, but in matrix of the dimension from tens to hundreds of, can also be observed that accurate asymptotic convergence result.
This is the premise that higher-dimension Random Matrices Theory is applied to engineering problem.
The content of the invention
The defects of for prior art, it is an object of the invention to provide it is a kind of independent of standard utensil based on higher-dimension with
The electronic mutual inductor error state appraisal procedure of machine matrix, it is ensured that electronic mutual inductor error state appraisal procedure can
Realisation, validity and accuracy.
The invention provides a kind of electronic mutual inductor error state appraisal procedure based on higher-dimension random matrix, including under
State step：
S1：The output valve of electronic mutual inductor is acquired based on time slip-window, and built according to the data of collection
Original matrix；
S2：Kalman filter based on single quantity of state is extended to the original matrix, obtains higher-dimension random matrix；
S3：The higher-dimension random matrix is standardized, make its be converted to row vector average be 0, variance be 1
Non- Hermite Matrix；
S4：Error state evaluating matrix is obtained according to the non-Hermite Matrix；
S5：The evaluation index of electric mutual inductor error state is obtained according to the error state evaluating matrix, and according to institute
State error state evaluating matrix and the evaluation index is assessed the error state of electric mutual inductor.
The present invention exports only according to electronic mutual inductor independent of standard utensil, judges this using higher-dimension random matrix
Whether a little parameters meet known rules, so as to which whether the statistical distribution of the measurement error to electronic mutual inductor has anomalous variation to do
Go out to judge, carry out error state assessment accordingly.
Further, in step sl, original matrix D is built1Based on time slip-window, in the state estimation of interception
Between in window, T times to N number of measured state parameter measurement of equipment, using measured state parameter number N as space dimensionality, measurement
Number T is as time dimension, and forms are with step-length TWIt is mobile.The original matrix of structure isWherein, xijValue of i-th of parameter in j measurement is represented, i is can
Survey state parameter sequence number, i=1,2 ... ... N, N be state parameter number, j be measurement sequence number, j=1,2 ... ... T, T
For pendulous frequency.
Further, in step s 2, based on Kalman filter, the square built according to electronic mutual inductor output valve
Battle array D1, using single state Kalman filter equation, Kalman filter output is obtained, using the output of wave filter as matrix
OK, with the dimension of extended matrix, matrix D is obtained2.The higher-dimension random matrix obtained after extensionN' be extension after state parameter number, N' value model
Enclose satisfaction
Further, in step s3, non-Hermite MatrixWhereinxi=(xi1,xi2,…,
xiT), 1≤i≤N' is matrix D2Row vector,Represent sample xiAverage value, σ (xij) represent sample xiStandard deviation.
Further, step S4 is specially：Calculate the singular value equivalent matrice D of the non-Hermite Matrixu；According to institute
State singular value equivalent matrice calculating matrix product Z；Error state evaluating matrix Z is obtained according to the matrix product Z2。
Further, the error state evaluating matrixWherein,
zi=(zi1,zi2,…,ziT), 1≤i≤N ',Represent error state evaluating matrix Z2Element,σ
(zi) represent ziStandard deviation.
Further, in step s 5, the evaluation index includes central moment αkWith moment of the orign βk；The central momentThe moment of the orignWherein, n is characterized value number, and E (V) is characterized value sample
It is expected.
Further, using 2 rank central moment α2With 1 rank moment of the orign β1Assessment as electronic mutual inductor error state
Criterion.
Further, after step s 5, also comprise the following steps：Judge to assess whether the time terminates, if so, then obtaining
Obtain electronic mutual inductor error state assessment result；If it is not, then t=t+TWIt is back to step S1；Wherein TWMoved for time slip-window
Dynamic step-length, t are the assessment time.
The present invention, can be real only according to the data of electronic mutual inductor itself output without using standard utensil at the scene
" self-evaluating " of existing error state, has been greatly reduced the cost assessed, and the O＆M for being advantageous to improve electronic mutual inductor is horizontal.
The present invention need not establish any physical model, without assumed condition and simplified condition, when the parameter of power system changes
The assessment result of the present invention will not be had any impact.No matter undergoing mutation property of error state failure or occur gradually changeable therefore
Barrier, the present invention can make accurate and effective assessment according to the regularity of distribution of characteristic value.
Brief description of the drawings
Fig. 1 is the estimation flow schematic diagram of the present invention.
Fig. 2 is electronic mutual inductor error state Evaluation Platform schematic diagram.
Fig. 3 is the distribution figure of characterized values under electronic mutual inductor normal condition.Wherein, discrete point is that error state assesses square
The feature Distribution value of battle array.
Fig. 4 is the distribution figure of characterized values under electronic mutual inductor malfunction feature.Wherein, discrete point is commented for error state
Estimate the feature Distribution value of matrix.
Wherein, 1 is AC 220V power supplys, and 2 be isolating transformer, and 3 be auto-transformer, and 4 be voltage transformer, and 5 be data
Collector, 6 be pulse per second (PPS) equipment, and 7 be computer.
Embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, it is right below in conjunction with drawings and Examples
The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and
It is not used in the restriction present invention.
The present invention exports only according to electronic mutual inductor, judges whether these parameters meet using higher-dimension random matrix
Know rule, so as to which whether the statistical distribution of the measurement error to electronic mutual inductor has anomalous variation to judge, and error
Statistical distribution and error state are directly related, carry out error state assessment accordingly.
In the present invention, the electronic mutual inductor error state appraisal procedure based on higher-dimension random matrix, comprise the following steps：
Step 1：The on-line monitoring of electronic mutual inductor output valve, based on time slip-window, built using output valve original
Matrix；
Electronic mutual inductor output is the Serial No. of time coherence, and measured state parameter contains the width of fundamental wave harmonic
The information such as value, phase, frequency, wherein, harmonic wave highest number is 13 times.In order to realize the real-time assessment of state, one can be used
The method of kind time slip-window, i.e. previous moment time window terminate the initial time that the time is subsequent time time window, constantly
The output valve at current time is obtained from source data, builds original matrix, time slip-window length span is 1000~∞
(s), it is preferable that time slip-window length can be selected as 3600s.
In the state estimation time window of a certain interception, to N number of measured state parameter measurement of electronic mutual inductor output
T times, using measured state parameter number N as space dimensionality, pendulous frequency T is as time dimension, and forms are with step-length TWIt is mobile.
State parameter N span is 1~39, it is preferable that N can be selected as 1；Pendulous frequency T span is 100~∞,
Preferably, T can be selected as 360；In the time window of interception, all measurement data may be constructed original matrix D1：
Wherein, xijRepresent i-th of parameter in j
The value of secondary measurement, i be measured state parameter sequence number, i=1,2 ... ... N, j be measurement sequence number, j=1,2 ... ... T.
Step 2：Kalman filter based on single quantity of state is to original matrix D1It is extended, obtains higher-dimension random matrix
D2；
When electronic mutual inductor error state based on higher-dimension Random Matrices Theory is assessed, due to characterizing electronic mutual inductor
The measured state parameter of error state is less, it is impossible to directly using measurement data structure higher-dimension random matrix.In order to overcome this
Problem, employ the Matrix extension method of the Kalman filter equation based on single quantity of state.
Estimate the accurate measurements of measuring system：Wherein, xkFor current time system
State space amount, xk+1For subsequent time system state space amount, ykFor system measurement；ξkIt is 0 mean value model noise；ηkIt is 0
Average measurement noise.
Change measurement noise ηkValue, the output valve of multigroup Kalman filter, η can be obtainedkSpan can be
0.1Vrms~10Vrms, wherein, VrmsFor the virtual value of Kalman filter input signal.Using the output of Kalman filter as
Row matrix, state parameter are changed into N' from N number of, and the dimension of technology transform, N' span needs to ensureIt is preferred that
Ground, N' can be selected as 150, so as to build higher-dimension random matrix D2：
Step 3：Higher-dimension random matrix is standardized, make its be converted to row vector average be 0, variance be 1
Non- Hermite Matrix；
To matrix D2It is changed into non-Hermite Matrix D after carrying out following normalizing operationstd：
Whereinxi
=(xi1,xi2,…,xiT), 1≤i≤N' is matrix D2Row vector,Represent sample xiAverage value, σ (xij) represent sample
xiStandard deviation.
So that the matrix D after normalizing operationstd=(yij)N'×TMeetWherein, yi=
(yi1,yi2,…,yiT),1≤i≤N’。
Step 4：Calculated by the calculating of singular value equivalent matrice, matrix product and evaluating matrix calculates link, obtain error
State estimation matrix；
First, the singular value equivalent matrice D of non-Hermite Matrix is asked foru：
Wherein,Representing matrix DstdConjugation means, U is Haar unitary matrice.
Then, calculating matrix product Z：
Wherein, Du,iIndependent singular value equivalent matrice is represented, L span is 1~∞, it is preferable that L is taken as 1.
Finally, based on matrix product Z, error state evaluating matrix Z is asked for2：
Wherein, zi=(zi1,
zi2,…,ziT), 1≤i≤N ',Represent error state evaluating matrix Z2Element,σ(zi) represent zi
Standard deviation.
Step 5：The evaluation index of electric mutual inductor error state is obtained, evaluation index includes central moment αkWith moment of the orign βk；
And according to error state evaluating matrix Z2The error state of electric mutual inductor is assessed with evaluation index；
The state of assessment equipment is judged according to the higher-dimension random matrix of structure, it is necessary to which whether test matrix meets some property
Matter, checkout procedure need to rely on some specific statistics.For monocyclic theorem, square statistic characteristic feature value can be used
The regularity of distribution, error state evaluating matrix Z2Characteristic value may be constructed characteristic value sample V={ V1,V2,…Vn, the k ranks of the sample
Central moment αkFor：Wherein, n is characterized value number, and E (V) is characterized the value sample phase
Hope, αkAs evaluation index J1, equally can be with the k rank moment of the origns β of defined feature value samplekAs evaluation index J2：
Preferably, 2 rank central moment α can be used2With 1 rank moment of the orign β1Assessment as electronic mutual inductor error state
Criterion.
Specific appraisal procedure：When electronic mutual inductor error state is normal, i.e., its error is in standard GB/T 20840.7-
2007《The part of transformer the 7th：Electronic type voltage transformer》With standard GB/T 20840.7-2007《The part of transformer the 8th：Electricity
Type current transformer》During defined electronic mutual inductor error limit requirement, measured state parameter keeps substantially constant constant,
Only make stable random fluctuation near average is horizontal, conversely, when electronic mutual inductor error state exception, measured state ginseng
There are unusual fluctuations in the measured value of amount, is no longer steady random fluctuation.According to monocyclic theorem, when electronic mutual inductor is in normal
During state, the singular value equivalent matrice for the random matrix that its state parameter measured value is formed will meet monocyclic theorem, characteristic value meeting
It is evenly distributed in the ring with specific interior outer radius；If abnormal, the square containing abnormal data occurs for electronic mutual inductor state
The feature Distribution value of battle array can change, and distribution is no longer uniform.Using electronic mutual inductor output valve, extract contain survey shape
State property, higher-dimension random matrix is built, obtains error state evaluating matrix, it is final to identify real-time electronic according to evaluation index
Formula transformer error state.After the completion of real-time electronic formula transformer error state estimation step, judge to assess whether the time ties
Beam, if so, then obtaining electronic mutual inductor error state assessment result；Continued executing with if it is not, being then back to step 1, during assessment
Between span be 1~7 day, it is preferable that assess the time can select as 1 day.
To further understand the present invention, relative theory in the present invention is briefly illustrated below：
(1) electronic mutual inductor output parameter Statistical Distribution Characteristics
Electronic mutual inductor gathers power network primary voltage and current information, by sampled data can be calculated frequency,
Phase, amplitude, power factor, harmonic wave, tri-phase unbalance factor etc., and can also be obtained between these parameters by mathematical operation
New parameter, i.e., almost all of operation of power networks parameter can directly or indirectly be observed by electronic mutual inductor.These electricity
The rule of gas parameter is accurate Analysis expression or can to obey known distribution etc., for example, the changing rule of node phase can parse
ForWherein, ω0For operation of power networks frequency,For the first phase of node phase.According to the defeated of electronic mutual inductor
Go out value, it can be determined that whether the parameter meets known rules, so as to which the statistical distribution of the measurement error to electronic mutual inductor is
No anomalous variation judges, and the statistical distribution and error performance of error are directly related, carries out error state assessment accordingly, this is
It is the general principle independent of standard utensil and physical model assessment electronics formula transformer error state.
(2) monocyclic theorem
Assuming that random matrix A=(aij)n×nIt is n rank Haar unitary matrice that A=PTQ, wherein P and Q, which can be decomposed into, and T is pair
Angle battle array, diagonal entry are A singular values.When certain conditions are met, the limit Spectral structure of matrix A is by the general of its singular value
Rate, which is estimated, to be uniquely determined, and characteristic value converges to annulus on a complex plane, and the interior outer radius of annulus is respectively：Wherein, v is the probability measure of the singular value of matrix A, and this is referred to as monocyclic theorem.
In actual applications, matrix is consideredFor non-Hermite Matrix, and the row vector of matrix A meets
Average is 0, variance 1.For multiple non-Hermite Matrix Ai, define matrix product A2ForWherein, Au,iFor Ai
Singular value equivalent matrice.By matrix A2It is standardized as Astd, it is met σ2(ai)=1/N, wherein, aiFor matrix AstdRow to
Measure, then AstdLimit Spectral structure convergence in (with)probability be to probability density function： Spectral distribution function.In formula (12),p,
M→∞。AstdCharacteristic value in the distribution of complex plane be an annulus, inner ring radius is (1-c)L/2, outer shroud radius is 1.Setting
In the case of standby state is normal, matrix AstdMeet following property：Singular value equivalent matrice converts what is obtained by Haar unitary matrice
Standardization product matrix should meet monocyclic theorem.
The invention will be further described with specific embodiment below in conjunction with the accompanying drawings.Embodiment is exemplary, it is intended to is used
It is of the invention in explaining, and be not considered as limiting the invention.
As shown in figure 1, the present invention is assessed the error state of electronic mutual inductor according to following steps：
(1) electronic mutual inductor error state Evaluation Platform as shown in Figure 2 is built, platform includes：AC 220V power supplys 1,
Isolating transformer 2, auto-transformer 3, voltage transformer 4, data acquisition unit 5, pulse per second (PPS) equipment 6 and computer 7.AC 220V
Power supply, isolating transformer, auto-transformer are used to provide signals to voltage transformer, because the error of electronic mutual inductor is not easy
Change, the simulation electronic formula transformer by the way of voltage transformer+data acquisition unit, pulse per second (PPS) equipment carries to data acquisition unit
Sampled for pps pulse per second signal with triggering.Voltage transformer data are passed through by being sent to computer, computer after data acquisition unit
Ethernet interface captures electronic mutual inductor packet, and original random matrix D is built according to electronic current mutual inductor output1。
(2) Kalman filter is based on, to matrix D1Higher-dimension random matrix D is formed after being extended2.Original matrix and expansion
The scale for opening up matrix is as shown in table 1.Due to only calculating electronic mutual inductor phase, the line number of original matrix is 1, and electronic type is mutual
Sensor phase calculates 1 time per 10s, and time slip-window length is 1h, calculates 360 times altogether, and this is the columns of original matrix, utilizes
After Kalman filter extension, the extended matrix of one 150 × 360 is constituted.
The higher dimensional matrix scale of table 1
(3) using formula (4) to matrix D2It is standardized and is converted to matrix Dstd。
(4) evaluating matrix Z is asked for using formula (6)~formula (8)2。
(5) method by being blocked in a computer to electronic mutual inductor data bits, simulation electronic formula mutual inductance
The normal condition and abnormal conditions of device error state, Fig. 3 are the characteristic value regularity of distribution for not blocking evaluating matrix in the case of digit,
Fig. 4 is the characteristic value regularity of distribution for blocking evaluating matrix in the case of digit is 10.From Fig. 3 and Fig. 4 contrasts as can be seen that when electricity
When minor transformer is in normal condition, the feature Distribution value of evaluating matrix is concentrated in an annulus；And work as electronic type mutual inductance
When device is in abnormality, characteristic value distribution starts discretization, and partial feature value fall into annulus.
Evaluation index is calculated according to formula (9)~formula (10), result of calculation is as shown in table 2.As can be seen from Table 2, when
When electronic mutual inductor is in abnormality, evaluation index J1And J2Value significantly become big, can be effective therefore according to evaluation index
Judge the error state of electronic mutual inductor.
The evaluation index of table 2 contrasts
| Evaluation index | Normal condition | Abnormality |
| J1 | 7.9128 | 10.7242 |
| J2 | 0.6841 | 1.1660 |
In order to as it will be easily appreciated by one skilled in the art that the foregoing is merely illustrative of the preferred embodiments of the present invention, not
To limit the present invention, all any modification, equivalent and improvement made within the spirit and principles of the invention etc. all should
Within protection scope of the present invention.
Claims (9)
1. a kind of electronic mutual inductor error state appraisal procedure based on higher-dimension random matrix, it is characterised in that including following
Step：
S1：The output valve of electronic mutual inductor is acquired based on time slip-window, and it is original according to the data of collection structure
Matrix；
S2：Kalman filter based on single quantity of state is extended to the original matrix, obtains higher-dimension random matrix；
S3：The higher-dimension random matrix is standardized, it is converted to the non-strategic point that row vector average is 0, variance is 1
The special matrix of rice；
S4：Error state evaluating matrix is obtained according to the non-Hermite Matrix；
S5：The evaluation index of electric mutual inductor error state is obtained according to the error state evaluating matrix, and according to the mistake
Poor state estimation matrix and the evaluation index are assessed the error state of electric mutual inductor.
2. electronic mutual inductor error state appraisal procedure as claimed in claim 1, it is characterised in that in step sl, structure
The original matrix built isWherein, xijRepresent i-th of parameter in j
The value of secondary measurement, i are the sequence number of measured state parameter, and i=1,2 ... ... N, N are the number of state parameter, and j is the sequence of measurement
Number, j=1,2 ... ... T, T are pendulous frequency.
3. electronic mutual inductor error state appraisal procedure as claimed in claim 2, it is characterised in that in step s 2, warp
The higher-dimension random matrix obtained is crossed after extensionN'
For the number of the state parameter after extension, N' span meets
4. electronic mutual inductor error state appraisal procedure as claimed in claim 3, it is characterised in that in step s3, institute
State non-Hermite MatrixWhereinxi=
(xi1,xi2,...,xiT), 1≤i≤N' is matrix D2Row vector,Represent sample xiAverage value, σ (xij) represent sample xi
Standard deviation.
5. the electronic mutual inductor error state appraisal procedure as described in claim any one of 1-4, it is characterised in that step S4
Specially：
Calculate the singular value equivalent matrice D of the non-Hermite Matrixu；
According to the singular value equivalent matrice calculating matrix product Z；
Error state evaluating matrix Z is obtained according to the matrix product Z2。
6. electronic mutual inductor error state appraisal procedure as claimed in claim 5, it is characterised in that the error state is commented
Estimate matrixWherein, zi=(zi1,zi2,...,ziT), 1≤i≤N ', zi
Representing matrix product Z row vector,Represent error state evaluating matrix Z2Row vector,σ(zi) table
Show ziStandard deviation, N' is the number of the state parameter after extension, and N' span meetsT is pendulous frequency.
7. electronic mutual inductor error state appraisal procedure as claimed in claim 1, it is characterised in that in step s 5, institute
Stating evaluation index includes central moment αkWith moment of the orign βk；The central momentThe moment of the orignWherein, away from the exponent number with moment of the orign, V={ V centered on k1,V2,…VnIt is error state evaluating matrix Z2's
Characteristic value sample, n are characterized value number,It is characterized the expectation of value sample.
8. electronic mutual inductor error state appraisal procedure as claimed in claim 7, it is characterised in that using 2 rank central moment α2
With 1 rank moment of the orign β1Evaluation criteria as electronic mutual inductor error state.
9. electronic mutual inductor error state appraisal procedure as claimed in claim 1, it is characterised in that after step s 5,
Also comprise the following steps：
Judge to assess whether the time terminates, if so, then obtaining electronic mutual inductor error state assessment result；If it is not, then t=t+
TWIt is back to step S1；Wherein TWFor the step-length of time slip-window movement, t is the assessment time.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610953170.0A CN106548021B (en) | 2016-11-03 | 2016-11-03 | A kind of electronic mutual inductor error state appraisal procedure based on higher-dimension random matrix |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610953170.0A CN106548021B (en) | 2016-11-03 | 2016-11-03 | A kind of electronic mutual inductor error state appraisal procedure based on higher-dimension random matrix |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106548021A CN106548021A (en) | 2017-03-29 |
| CN106548021B true CN106548021B (en) | 2018-01-05 |
Family
ID=58393574
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610953170.0A Active CN106548021B (en) | 2016-11-03 | 2016-11-03 | A kind of electronic mutual inductor error state appraisal procedure based on higher-dimension random matrix |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106548021B (en) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107132454A (en) * | 2017-05-04 | 2017-09-05 | 国网上海市电力公司 | The abnormal quick determination method of power network based on random matrix spectral radius method |
| CN109033031B (en) * | 2018-07-24 | 2022-08-23 | 南通大学 | Bearing state detection method based on high-dimensional random matrix |
| CN109359822B (en) * | 2018-09-21 | 2022-02-15 | 国网江苏省电力有限公司电力科学研究院 | Electronic voltage transformer measurement state evaluation method and system |
| CN110031788A (en) * | 2019-04-04 | 2019-07-19 | 国网江苏省电力有限公司电力科学研究院 | A kind of hollow coil current transformer error environment correlation analysis |
| CN110095744B (en) * | 2019-04-04 | 2021-05-11 | 国网江苏省电力有限公司电力科学研究院 | Error prediction method for electronic transformer |
| CN111474510B (en) * | 2020-04-25 | 2021-02-12 | 华中科技大学 | Error evaluation method and system for voltage transformer with non-stable output |
| CN112485747A (en) * | 2020-11-28 | 2021-03-12 | 国网江西省电力有限公司供电服务管理中心 | Error state evaluation method for capacitor voltage transformer |
| CN113050018B (en) * | 2021-03-04 | 2023-07-07 | 国网湖南省电力有限公司 | Voltage transformer state evaluation method and system based on data driving evaluation result change trend |
| CN115542064B (en) * | 2022-12-01 | 2023-04-07 | 威胜集团有限公司 | Real-time monitoring method and device for error state of mutual inductor |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN201145737Y (en) * | 2008-01-04 | 2008-11-05 | 宁波三维电测设备有限公司 | Apparatus for measuring error of current transformer based on parameter measurement |
| CN101726716A (en) * | 2009-12-17 | 2010-06-09 | 江苏省电力公司镇江供电公司 | Method for measuring error of voltage transformer by nomography |
| CN104933012A (en) * | 2015-05-08 | 2015-09-23 | 西南交通大学 | Method for online identifying measurement deviation fault for traction substation instrument transformer |
| CN105699804A (en) * | 2016-01-22 | 2016-06-22 | 吉林大学 | Big data fault detection and positioning method for power distribution network |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2001029572A1 (en) * | 1999-10-21 | 2001-04-26 | Siemens Power Transmission & Distribution, Inc. | External transformer correction in an electricity meter |
-
2016
- 2016-11-03 CN CN201610953170.0A patent/CN106548021B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN201145737Y (en) * | 2008-01-04 | 2008-11-05 | 宁波三维电测设备有限公司 | Apparatus for measuring error of current transformer based on parameter measurement |
| CN101726716A (en) * | 2009-12-17 | 2010-06-09 | 江苏省电力公司镇江供电公司 | Method for measuring error of voltage transformer by nomography |
| CN104933012A (en) * | 2015-05-08 | 2015-09-23 | 西南交通大学 | Method for online identifying measurement deviation fault for traction substation instrument transformer |
| CN105699804A (en) * | 2016-01-22 | 2016-06-22 | 吉林大学 | Big data fault detection and positioning method for power distribution network |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106548021A (en) | 2017-03-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106548021B (en) | A kind of electronic mutual inductor error state appraisal procedure based on higher-dimension random matrix | |
| CN107016236B (en) | Power grid false data injection attack detection method based on nonlinear measurement equation | |
| CN106353623B (en) | Low-frequency oscillation of electric power system mode on-line identification method based on random response signal | |
| Kong et al. | Online smart meter measurement error estimation based on EKF and LMRLS method | |
| CN112485747A (en) | Error state evaluation method for capacitor voltage transformer | |
| Wang et al. | Transmission lines positive sequence parameters estimation and instrument transformers calibration based on PMU measurement error model | |
| CN111474510A (en) | Error evaluation method and system for voltage transformer with non-stable output | |
| CN110208647B (en) | Power distribution network single-phase earth fault section positioning method based on random matrix | |
| CN105098763A (en) | Transient voltage stable on-line assessment method based on WAMS and EMS | |
| CN109284933B (en) | Electronic transformer state evaluation system and method based on mathematical statistics | |
| CN105406461B (en) | A kind of load self-adapting dynamic monitoring method of power distribution network power-off event | |
| Liu et al. | Remote malfunctional smart meter detection in edge computing environment | |
| CN110348114B (en) | Non-precise fault identification method for power grid completeness state information reconstruction | |
| CN105044649B (en) | A kind of current transformer error horizontal dynamic detection method and system | |
| Luan et al. | Energy theft detection via integrated distribution state estimation based on AMI and SCADA measurements | |
| CN111551887A (en) | Multidimensional identification voltage transformer metering performance online monitoring platform | |
| CN106841857A (en) | A kind of equipment for monitoring power quality reliability estimation method | |
| Wang et al. | Augmented state estimation of line parameters in active power distribution systems with phasor measurement units | |
| CN106407955B (en) | A kind of low-frequency oscillation of electric power system mode on-line identification method | |
| Li et al. | Dynamic topology awareness in active distribution networks under DG uncertainties using GMM-PSEs and KL divergence | |
| CN105068035A (en) | Voltage transformer error level dynamic detection method and system | |
| CN109902735A (en) | A kind of two norm Independent Component Analysis for realizing that threephase potential transformer kinematic error is assessed in substation | |
| CN107832959B (en) | Voltage stability evaluation method considering load characteristics and power supply constraints | |
| CN109613324A (en) | A kind of detection method and device of Harmonics amplification | |
| Li et al. | Real time evaluation algorithm for measurement performance of substation voltage transformer based on artificial neural 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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20211123 Address after: Room 1803-1805, building 2-07, guanggu.core center, 303 Guanggu Avenue, Donghu New Technology Development Zone, Wuhan City, Hubei Province, 430000 Patentee after: WUHAN GLORY ROAD INTELLIGENT TECHNOLOGY Co.,Ltd. Address before: 430074 Hubei Province, Wuhan city Hongshan District Luoyu Road No. 1037 Patentee before: HUAZHONG University OF SCIENCE AND TECHNOLOGY |
|
| CP01 | Change in the name or title of a patent holder | ||
| CP01 | Change in the name or title of a patent holder |
Address after: Room 1803-1805, building 2-07, guanggu.core center, 303 Guanggu Avenue, Donghu New Technology Development Zone, Wuhan City, Hubei Province, 430000 Patentee after: Wuhan Gelanruo Intelligent Technology Co.,Ltd. Address before: Room 1803-1805, building 2-07, guanggu.core center, 303 Guanggu Avenue, Donghu New Technology Development Zone, Wuhan City, Hubei Province, 430000 Patentee before: WUHAN GLORY ROAD INTELLIGENT TECHNOLOGY Co.,Ltd. |