CN110031788A - A kind of hollow coil current transformer error environment correlation analysis - Google Patents

A kind of hollow coil current transformer error environment correlation analysis Download PDF

Info

Publication number
CN110031788A
CN110031788A CN201910270802.7A CN201910270802A CN110031788A CN 110031788 A CN110031788 A CN 110031788A CN 201910270802 A CN201910270802 A CN 201910270802A CN 110031788 A CN110031788 A CN 110031788A
Authority
CN
China
Prior art keywords
matrix
current transformer
coil current
hollow coil
parameter
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.)
Pending
Application number
CN201910270802.7A
Other languages
Chinese (zh)
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.)
Huazhong University of Science and Technology
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
Electric Power Research Institute of State Grid Jiangxi Electric Power Co Ltd
Yangzhou Power Supply Co of Jiangsu Electric Power Co
Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
Huazhong University of Science and Technology
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
Electric Power Research Institute of State Grid Jiangxi Electric Power Co Ltd
Yangzhou Power Supply Co of Jiangsu Electric Power Co
Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huazhong University of Science and Technology, State Grid Corp of China SGCC, State Grid Jiangsu Electric Power Co Ltd, Electric Power Research Institute of State Grid Jiangxi Electric Power Co Ltd, Yangzhou Power Supply Co of Jiangsu Electric Power Co, Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd filed Critical Huazhong University of Science and Technology
Priority to CN201910270802.7A priority Critical patent/CN110031788A/en
Publication of CN110031788A publication Critical patent/CN110031788A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R35/00Testing or calibrating of apparatus covered by the other groups of this subclass
    • G01R35/02Testing 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Complex Calculations (AREA)

Abstract

The invention discloses a kind of hollow coil current transformer error environment correlation analysis, include the following steps: in assessment time window according to acquisition data building original matrix;The original matrix is extended based on Kalman filter, establishes higher-dimension random matrix;The higher-dimension random matrix is standardized, so that it is converted to row vector mean value 0, the non-Hermite Matrix that variance is 1;Influence amount relevance evaluation matrix is obtained according to the non-Hermite Matrix;Hollow coil current transformer error environment relevance evaluation index is obtained according to the influence amount relevance evaluation matrix, and the correlation between hollow coil current transformer error and environment parameter is assessed according to the influence amount relevance evaluation matrix and the relevance evaluation index.Advantage: can obtain the correlation degree of mutual inductor kinematic error Yu one or more environment parameter in real time, be conducive to control and assess the running error state stability of mutual inductor.

Description

A kind of hollow coil current transformer error environment correlation analysis
Technical field
Invention belongs to equipment for power transmission and distribution status assessment field, more particularly, to a kind of sky based on higher dimensional matrix theory Wire-core coil current transformer error environmental dependence analysis method.
Background technique
Voltage current transformer is that the important measuring device of current information, maneuverability are provided for substation's electric energy metered system The accuracy of metering device can be related to.Not only insulation system is complicated, volume is big for the conventional electromagnetic electric current being widely used at present, Cost is high, there is also magnetic saturation, dynamic range is small the disadvantages of, it is difficult to meet the technical need of electric system.Electronic current is mutual Sensor has many advantages, such as that dynamic range is big, Hz-KHz is wide, small in size, light weight, has complied with electric power digital, intelligence With the developing direction of networking.Hollow coil current transformer is one kind of electronic current mutual inductor, recently as intelligence The construction of substation is rapidly developed, and by the exploration and practice of many years, achieves a series of achievements.
It is numerous comprising link however, hollow coil current transformer structure is complex, in operational process by temperature, The influence of the environment parameters such as humidity, vibration, magnetic field and a load current.From the point of view of live operation problem, air core coil electricity The accuracy problems of current transformer still occupy biggish ratio.The fairness for affecting electric energy trade settlement, leads to tubular wire The popularization and application of loop current mutual inductor are hindered.It discloses the inner link of transformer error and various environment parameters and influences to advise Rule, specifies main influence amount, guidance is provided for the design and processes of hollow coil current transformer, to hollow coil current mutual inductance The control of device error robustness and assessment are of great significance.
The prior art includes the correlation analysis based on model, according to environment parameter to hollow coil current transformer The mechanism model of effect, the mechanism and rule that analysis environment parameter influences transformer error, this method height rely on model Accuracy, various hypothesis and premise also will affect analysis as a result, needing to establish different mechanism moulds for different environment parameters Type, versatility is poor, and is unable to get correlation quantitative assessing index.
The prior art further includes the analysis method based on data-driven, does not need to construct accurate mechanism model, passes through digging The error information and environment parameter data of pick, processing and analysis mutual inductor, it is related to environment parameter to obtain transformer error Property.However, fortune hollow coil current transformer error and environment parameter between relationship present it is multiple coupled and it is high at random Feature determines that the degree of correlation of transformer error and environment parameter is more difficult, and above-mentioned analysis method is simultaneously not suitable for.
Summary of the invention
The technical problem to be solved by the present invention is to overcoming the deficiencies of existing technologies, a kind of mechanism model that do not depend on is provided Hollow coil current transformer error environment correlation analysis based on higher dimensional matrix theory can be hollow coil current The control method of transformer error provides reference.
In order to solve the above technical problems, the present invention provides a kind of hollow coil current transformer error environment correlation analysis Method, which comprises the steps of:
S1: acquisition environment parameter data and hollow coil current transformer error information, the basis in assessment time window The environment parameter data and error information of acquisition construct original matrix;
S2: the original matrix is extended based on Kalman filter, establishes higher-dimension random matrix;
S3: being standardized the higher-dimension random matrix, so that it is converted to row vector mean value 0, variance is 1 Non- Hermite Matrix;
S4: influence amount relevance evaluation matrix is obtained according to the non-Hermite Matrix;
S5: hollow coil current transformer error environmental dependence is obtained according to the influence amount relevance evaluation matrix and is commented Estimate index, and according to the influence amount relevance evaluation matrix and the relevance evaluation index to hollow coil current transformer Correlation between error and environment parameter is assessed.
Further, in step sl, pass through the environment parameter data constructing environment parameter matrix of acquisitionWherein, element PijIt indicates that environment parameter can be surveyed in the measured value at j moment i-th, i is It can survey the serial number of environment parameter, i=1,2 ... ... M, M are the number of environment parameter, and j be the serial number measured, j=1,2 ... ... T, T is pendulous frequency;Error state matrix is constructed by the error information of acquisitionWherein, Element QijIndicating measured value of i-th of transformer error parameter at the j moment, i is the serial number of transformer error parameter, i=1, 2 ... ... N, N are the number of transformer error parameter, and j is the serial number of measurement, and the original matrix of j=1,2 ... ... T, building areWherein, k=M+N.
Further, in step s 2, the higher-dimension random matrix that obtains is after extensionK' is the number of the state parameter after extension, and the value range of N' is full Sufficient k'/T ∈ (0,1], T is pendulous frequency.
Further, in step s3, the non-Hermite Matrix isWherein Indicate sample xiBe averaged Value, σ (xij) indicate sample xiStandard deviation, xiFor higher-dimension random matrix D3Row vector, xi=(xi1,xi2,...,xiT),1≤i ≤ k', k' are the number of the state parameter after extension, and T is pendulous frequency, yijFor higher-dimension random matrix D3In variable xijBy The new variable obtained after the standardized way.
Further, the step S4 specifically:
Calculate the singular value equivalent matrice D of the non-Hermite Matrixu
According to the singular value equivalent matrice DuCalculating matrix product Z;
Error state evaluating matrix Z is obtained according to the matrix product Z2
Further, the matrix productL=1, institute State error state evaluating matrixWherein,ziFor square The row vector of battle array Z, zi=(zi1,zi2,...,ziT), 1≤i≤k',For matrix Z2Row to Amount, σ (zi) indicate ziStandard deviation, k' be extension after state parameter number, T is pendulous frequency.
Further, in step s 5, the relevance evaluation index includes dMSRAnd IMSR, wherein dMSRevref, dMSRIntegral to the time is IMSR:Wherein, t1And t2Indicate initial time and the finish time of assessment,λiFor the characteristic value of corresponding original matrix, λwiFor correspondence The characteristic value of R-matrix, n and n2The characteristic value number of respectively corresponding original matrix and R-matrix, E () indicate feature It is worth sample expectation, the R-matrixR-matrix is by air core coil error state matrix and Gauss white noise Sound matrix is constituted, wherein D1For environment parameter matrix;DNFor noise matrix, dimension is identical with environment parameter matrix, and element is The stochastic variable of standardized normal distribution is obeyed, the white Gaussian noise amplitude being superimposed in amplitude and matrix-expand is identical.
A kind of hollow coil current transformer error state monitoring system, which is characterized in that construct mould including original matrix Block, higher-dimension random matrix building module, standardization module, influence amount relevance evaluation matrix module and relevance evaluation Module;
The original matrix building module is for acquiring environment parameter data and hollow coil current transformer margin of error According to according to the environment parameter data and error information of acquisition building original matrix in assessment time window;
The higher-dimension random matrix building module is used to be extended the original matrix based on Kalman filter, builds Vertical higher-dimension random matrix;
The standardization module for being standardized to the higher-dimension random matrix, make its be converted to row to Measure mean value be 0, the non-Hermite Matrix that variance is 1;
The influence amount relevance evaluation matrix module is used to obtain influence amount correlation according to the non-Hermite Matrix Evaluating matrix;
The relevance evaluation module is used to obtain hollow coil current according to the influence amount relevance evaluation matrix mutual Sensor error environment relevance evaluation index, and according to the influence amount relevance evaluation matrix and the relevance evaluation index Correlation between hollow coil current transformer error and environment parameter is assessed.
Further, the original matrix building module is used for the environment parameter data constructing environment parameter square by acquisition Battle arrayWherein, element PijIndicate can to survey for i-th environment parameter in the measured value at j moment, I is the serial number that can survey environment parameter, i=1, and 2 ... ... M, M are the number of environment parameter, and j be the serial number measured, j=1, 2 ... ... T, T are pendulous frequency;Error state matrix is constructed by the error information of acquisitionWherein, element QijIndicate measured value of i-th of transformer error parameter at the j moment, i For the serial number of transformer error parameter, i=1,2 ... ... N, N are the number of transformer error parameter, and j is the serial number of measurement, j= The original matrix of 1,2 ... ... T, building isWherein, k=M+N.
Further, the higher-dimension random matrix of the higher-dimension random matrix building module foundation isK' is the number of the state parameter after extension, and the value range of N' is full Sufficient k'/T ∈ (0,1], T is pendulous frequency.
Further, the non-Hermite Matrix that the standardization resume module obtains is that the non-Hermite Matrix isWherein Indicate sample xiBe averaged Value, σ (xij) indicate sample xiStandard deviation, xiFor higher-dimension random matrix D3Row vector, xi=(xi1,xi2,...,xiT),1≤i ≤ k', k' are the number of the state parameter after extension, and T is pendulous frequency, yijFor higher-dimension random matrix D3In variable xijBy The new variable obtained after the standardized way.
Further, the influence amount relevance evaluation matrix module is used to calculate the singular value of the non-Hermite Matrix Equivalent matrice Du, according to the singular value equivalent matrice DuCalculating matrix product Z obtains error state according to the matrix product Z Evaluating matrix Z2
The matrix productThe error shape State evaluating matrixWherein,ziFor the row of matrix Z Vector, zi=(zi1,zi2,...,ziT), 1≤i≤k',For matrix Z2Row vector, σ (zi) indicate ziStandard deviation, k' be extension after state parameter number, T is pendulous frequency
Further, the relevance evaluation index includes dMSRAnd IMSR, wherein dMSRevref, dMSRTo the product of time It is divided into IMSR:Wherein, t1And t2Indicate initial time and the finish time of assessment,λiFor the characteristic value of corresponding original matrix, λwiFor correspondence The characteristic value of R-matrix, n and n2The characteristic value number of respectively corresponding original matrix and R-matrix, E () indicate feature It is worth sample expectation, the R-matrixR-matrix is by air core coil error state matrix and Gauss white noise Sound matrix is constituted, wherein D1For environment parameter matrix;DNFor noise matrix, dimension is identical with environment parameter matrix, and element is The stochastic variable of standardized normal distribution is obeyed, the white Gaussian noise amplitude being superimposed in amplitude and matrix-expand is identical
Advantageous effects of the invention:
The present invention is without establishing any physical model, without assumed condition and simplified condition, according only to hollow coil current The inner link of transformer error and environment parameter is quantified as relativity evaluation by data of transformer error and environment parameter data Index, according to the size and variation tendency of relativity evaluation index, can obtain in real time mutual inductor kinematic error and one or The correlation degree of multiple environment parameters is conducive to control and assess the running error state stability of mutual inductor.
Detailed description of the invention
Fig. 1 is estimation flow schematic diagram of the invention;
Fig. 2 is hollow coil current transformer error state on-line monitoring platform schematic diagram;
Fig. 3 is influence amount relevance evaluation matrix Dev1Distribution figure of characterized values;
Fig. 4 is influence amount relevance evaluation matrix Dev2Distribution figure of characterized values;
Fig. 5 is corresponding influence amount relevance evaluation matrix Dev1Relevance evaluation index trend chart;
Fig. 6 is corresponding influence amount relevance evaluation matrix Dev2Relevance evaluation index trend chart.
Wherein, 1 is hollow coil current transformer, and 2 be electromagnetic current transducer, and 3 be environmental monitoring unit, and 4 be light Fine remote transmission unit, 5 be signal acquisition unit, and 6 be data processing unit, and 7 be time synchronization unit, and 8 be interchanger, and 9 be service Device.
Specific embodiment
The invention will be further described below in conjunction with the accompanying drawings.Following embodiment is only used for clearly illustrating the present invention Technical solution, and not intended to limit the protection scope of the present invention.
The present invention only needs to establish higher-dimension according to the error information and environment parameter data of hollow coil current transformer random Matrix obtains relevance evaluation index by higher-dimension random matrix, and relevance evaluation index reflects higher-dimension random matrix element Statistical Distribution can be used to characterize the correlation between transformer error state and environment parameter, accordingly can be to mutual inductance The correlation of device error state is analyzed.
In the present invention, the hollow coil current transformer error environment correlation analysis based on higher dimensional matrix theory, The following steps are included:
Step 1: the on-line monitoring of hollow coil current transformer error and environment parameter, in assessment time window, building Environment parameter matrix, error state matrix and original matrix;
There is time varying characteristic due to monitoring obtained time series data on-line, using the time slip-window side of processing in real time The termination time of method, i.e. previous moment time window is the initial time of subsequent time time window, obtains current time and history The environment parameter and error information at moment, the data of each sampling instant are arranged according to time series, current time and history The data at moment may be embodied in original matrix.Time slip-window length value range is 100~∞ (s), it is preferable that sliding Time window length can choose as 1800s.
It includes non-electrical parameter that the error pattern of hollow coil current transformer, which includes than difference and angular difference, environment parameter type, And electric parameter, electric parameter can be divided into magnetic field parameter and a load current parameter (referred to as load parameter), non-electrical Parameter can be divided into temperature parameters, humidity parameter and Vibration Parameter.In the assessment time window of each interception, to M environment Parameter measurement T times, N number of data of transformer error measured T times.The value range of M is 1~5, it is preferable that M is selected as 5; The value range of N is 1~2, it is preferable that N is selected as 2;The value range of T is 300~∞, it is preferable that T is selected as 360.It is all Measurement data may be constructed original matrix D:
Wherein, k=M+N, xijIndicating the value of i-th of parameter jth time measurement, i is the serial number of parameter, i=1,2 ... ... k, J is the serial number of pendulous frequency, j=1,2 ... ... T.
Step 2: the Kalman filter based on single quantity of state is extended original matrix D, obtains higher-dimension random matrix D3
Since the error information of hollow coil current transformer and the type of environment parameter are less, even if combining the two Afterwards, the dimension of constructed influence amount relevance evaluation matrix is still less, is unable to satisfy the building condition of higher-dimension random matrix. In order to solve this problem, using the Matrix extension method of the Kalman filter equation based on single quantity of state.
Based on Kalman filter equation, the accurate measurements of measuring system are estimated are as follows: Wherein, xkFor current time system state space amount, xk+1For subsequent time system state space amount, ykFor system measurement;ξk For 0 mean value model noise;ηkNoise is measured for 0 mean value.
Change measurement noise ηkValue, the output valve of available multiple groups Kalman filter, ηkValue range 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 become that k ' is a from k, the dimension of technology transform, the value range of k ' need to guarantee k'/T ∈ (0,1], Preferably, k ' is selected as 20, constructs higher-dimension random matrix D accordingly3:
Step 3: to higher-dimension random matrix D3It is standardized, it is made to be converted to row vector mean value 0, variance 1 Non- Hermite Matrix;
To matrix D3Become non-Hermite Matrix D after carrying out following normalizing operationstd:Wherein Indicate sample xiAverage value, σ (xij) indicate sample xiStandard deviation, xi=(xi1,xi2,...,xiT), 1≤i≤N' is matrix D3Row to Amount.So that the matrix D after normalizing operationstd=(yij)k'×TMeetWherein, yi=(yi1, yi2,...,yiT),1≤i≤N’。
Step 4: being calculated by singular value equivalent matrice, matrix product calculates and evaluating matrix calculates link, foundation influence Measure relevance evaluation matrix;
Firstly, seeking the singular value equivalent matrice D of non-Hermite Matrixu:
Wherein,Representing matrix DstdConjugation means, U be Ha Er unitary matrice.
Then, calculating matrix product Z: Wherein, Du,iIndicate each independent singular value equivalent matrice, the value range of L is 1~∞, it is preferable that L is taken as 1.
Finally, being based on matrix product Z, influence amount relevance evaluation matrix Z is obtained2:
Wherein, zi=(zi1,zi2,...,ziT), 1≤i≤k' is the row vector of matrix Z,For matrix Z2's Row vector, σ (zi) indicate ziStandard deviation.
Step 5: establishing hollow coil current transformer error environment relevance evaluation index, evaluation index linear character Valued Statistics indicate;According to influence amount relevance evaluation matrix Z2Hollow coil current transformer is missed with relevance evaluation index Difference is analyzed with the correlation between environment parameter.
Linear character Valued Statistics are able to reflect the characteristic value distribution situation an of random matrix, for a random matrix For, single feature value can not reflect the statistical law of matrix element in assessment time window, and the mark of matrix is able to reflect matrix The statistical nature of element.R-matrix may be constructed by air core coil error state matrix and white Gaussian noise matrixWherein, DNFor noise matrix, dimension is identical with environment parameter matrix, and element is to obey standardized normal distribution Stochastic variable, the white Gaussian noise amplitude being superimposed in amplitude and matrix-expand is identical.The characteristic value of R-matrix may be constructed Characteristic value sampleError state evaluating matrix Z2Characteristic value may be constructed characteristic value sample V={ λ1, λ2,…λn, the central moment of calculating matrixWith Wherein, λwiFor the characteristic value of corresponding R-matrix, λiFor the characteristic value of corresponding original matrix, n and n2Respectively correspond to original square The characteristic value number of battle array and R-matrix.E () indicates the expectation of characteristic value sample.Define relevance evaluation index dMSR: dMSR= εevref, dMSRIntegral to the time is IMSR:Wherein t1And t2Indicate initial time and the end of assessment Moment.
If there are correlation between the error of hollow coil current transformer and environment parameter, by data of transformer error and The singular value equivalent matrice of the random matrix of environment parameter data building will meet monocycle theorem, and characteristic value can be evenly distributed on tool In the ring for having specific interior outer radius;Otherwise feature Distribution value can change, and be distributed no longer uniform.
To further understand the present invention, monocycle theorem in the present invention is briefly illustrated below:
In practical applications, if matrixFor non-Hermite Matrix, and the row vector of matrix A meet it is equal Value is 0, variance 1.For multiple non-Hermite Matrix Ai, define matrix productWherein, Au,iFor AiIt is unusual It is worth equivalent matrice.By matrix A2It is standardized as Astd, it is made to meet σ2(ai)=1/n, wherein aiFor matrix AstdRow vector, then AstdLimit Spectral structure convergence with probability one to probability density function are as follows: In formula (12), and c=m/n ∈ (0,1], m, n → ∞.AstdCharacteristic value in the distribution of complex plane be an annulus, inner ring Radius is (1-c)L/2, the radius of outer ring is 1.In the normal situation of equipment state, matrix AstdMeet following property: singular value The standardization product matrix that equivalent matrice is converted by Ha Er unitary matrice should meet monocycle theorem.
The present invention will be further explained below with reference to the attached drawings and specific examples.Embodiment is exemplary, it is intended to be 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 according to the following steps between hollow coil current transformer error and environment parameter Correlation is analyzed:
(1) hollow coil current transformer error state monitoring platform as shown in Figure 2 is built, platform includes: environment prison Survey unit 3, optical fiber remote transmission unit 4, signal acquisition unit 5, data processing unit 6, time synchronization unit 7.It is equipped in platform One 0.2 grade of hollow coil current transformer 1 and one 0.2 grade of electromagnetic current transducer 2.With electromagnetic type Current Mutual Inductance The output of device 2 is standard signal, the comparison result of available 1 error of hollow coil current transformer.Environmental monitoring unit 3 can be right The environment parameter of mutual inductor installation place is acquired, including the parameters such as temperature, humidity, vibration, magnetic field;Optical fiber remote transmission unit 4 is then By the data normalization of environmental monitoring unit 3, it is sent to a data processing unit 6;Data are passed through friendship by data processing unit 6 It changes planes and 8 is transferred to server 9, monitoring data are stored in server 9;Signal acquisition unit 5 can acquire digitlization electricity The output data of magnetic-type current transformer;Data processing unit 6 receives the output data and tubular wire of signal acquisition unit 5 simultaneously The sampling value message data of loop current mutual inductor 1.It is constructed according to the error information of environment parameter and hollow coil current transformer 1 Original random matrix D;Clock synchronization unit 7 constructs the synchronized clock system of whole system, is responsible for synchronization optical fiber remote transmission unit 4, data processing unit 6 and signal acquisition unit 5.
(2) it is based on Kalman filter, forms higher-dimension random matrix D after being extended to matrix D3.Utilize air core coil The ratio difference data of current transformer, angular difference data, non-electrical parametric data, electric parameter data constitute original matrix, using base Each original matrix is extended in the matrix-expand method of Kalman filter, the scale of original matrix and extended matrix such as table Shown in 1.Hollow coil current transformer error and the every 10min of environment parameter are calculated once, and time slip-window length is for 24 hours, altogether It calculates 144 times, this is the columns of original matrix, after extending using Kalman filter, constitutes one 20 × 144 extension Matrix.
Due to only calculating electronic mutual inductor phase, the line number of original matrix is 1, and the every 10s of electronic mutual inductor phase It calculates 1 time, time slip-window length is 1h, amounts to and calculates 360 times, this is the columns of original matrix, utilizes Kalman filter After extension, one 150 × 360 extended matrix is constituted.
1 higher dimensional matrix scale of table
(3) using formula (4) to matrix D3It is standardized and is converted to matrix Dstd
(4) influence amount relevance evaluation matrix Z is sought using formula (6)~formula (8)2
(5) error state matrix is constituted using the ratio difference data of hollow coil current transformer, utilizes temperature parameters data Error state matrix and environment parameter matrix are merged into influence amount relevance evaluation matrix D by composing environment parameter matrixev1, square Battle array scale is 40 × 144;It is related to humidity parametric data composition influence amount using the ratio difference data of hollow coil current transformer Property evaluating matrix Dev2, matrix size is similarly 40 × 144.It is respectively D that time slip-window, which is chosen for 1800s, Fig. 3 and Fig. 4,ev1With Dev2Feature Distribution value, comparison diagram 3 and Fig. 4, it can be seen that Dev1Singular value equivalent matrice feature Distribution value it is more dispersed, Partial feature value has exceeded the limitation of annulus;Dev2The feature Distribution value of singular value equivalent matrice more concentrate, be substantially distributed in In one annulus.
Relevance evaluation index is calculated according to formula (9)~formula (10), Fig. 5 and Fig. 6 are respectively according to Dev1And Dev2It obtains Relevance evaluation index trend chart, it can be seen that for evaluating matrix Dev1For, evaluation index dMSRMaximum value on It is raised near 0.35, IMSRReach 273.15;For evaluating matrix Dev2For, evaluation index dMSRIt remains near 0, IMSR It is 43.8, is much smaller than matrix Dev1IMSR, this demonstrate the correlation of the ratio difference of hollow coil current transformer and temperature is stronger, It is weaker with the correlation of humidity.
It should be understood by those skilled in the art that, embodiments herein can provide as method, system or computer program Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the application Apply the form of example.Moreover, it wherein includes the computer of computer usable program code that the application, which can be used in one or more, The computer program implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) produces The form of product.
The application is referring to method, the process of equipment (system) and computer program product according to the embodiment of the present application Figure and/or block diagram describe.It should be understood that can be realized by computer program instructions each in flowchart and/or the block diagram The combination of process and/or box in process and/or box and flowchart and/or the block diagram.It can provide these computer journeys Sequence instruct to general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices processor with A machine is generated, so that the instruction generation executed by computer or the processor of other programmable data processing devices is used for Realize the dress for the function of specifying in one or more flows of the flowchart and/or one or more blocks of the block diagram It sets.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates, Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one The step of function of being specified in a box or multiple boxes.The above is only a preferred embodiment of the present invention, it is noted that For those skilled in the art, without departing from the technical principles of the invention, if can also make Dry improvement and deformation, those modifications and variations should also be regarded as the protection scope of the present invention.

Claims (13)

1. a kind of hollow coil current transformer error environment correlation analysis, which comprises the steps of:
S1: acquisition environment parameter data and hollow coil current transformer error information, according to acquisition in assessment time window Environment parameter data and error information construct original matrix;
S2: the original matrix is extended based on Kalman filter, establishes higher-dimension random matrix;
S3: being standardized the higher-dimension random matrix, so that it is converted to row vector mean value 0, the non-strategic point that variance is 1 The special matrix of rice;
S4: influence amount relevance evaluation matrix is obtained according to the non-Hermite Matrix;
S5: hollow coil current transformer error environment relevance evaluation is obtained according to the influence amount relevance evaluation matrix and is referred to Mark, and according to the influence amount relevance evaluation matrix and the relevance evaluation index to hollow coil current transformer error Correlation between environment parameter is assessed.
2. hollow coil current transformer error environment correlation analysis according to claim 1, which is characterized in that In step sl, pass through the environment parameter data constructing environment parameter matrix of acquisitionIts In, element PijIt indicating to survey environment parameter i-th in the measured value at j moment, i is the serial number that can survey environment parameter, i=1, 2 ... ... M, M are the number of environment parameter, and j is the serial number of measurement, and j=1,2 ... ... T, T are pendulous frequency;Pass through the mistake of acquisition Difference data constructs error state matrixWherein, element QijIndicate that i-th of mutual inductor misses Measured value of the poor parameter at the j moment, i are the serial number of transformer error parameter, and i=1,2 ... ... N, N are transformer error parameter Number, j is the serial number of measurement, and the original matrix of j=1,2 ... ... T, building is Wherein, k=M+N.
3. hollow coil current transformer error environment correlation analysis according to claim 1, which is characterized in that In step s 2, the higher-dimension random matrix that obtains is after extension K' be extension after state parameter number, the value range of N' meet k'/T ∈ (0,1], T is pendulous frequency.
4. hollow coil current transformer error environment correlation analysis according to claim 1, which is characterized in that In step s3, the non-Hermite Matrix isWherein Indicate sample xiAverage value, σ (xij) indicate sample xiStandard deviation, xiFor the random square of higher-dimension Battle array D3Row vector, xi=(xi1,xi2,...,xiT), 1≤i≤k', k' are the number of the state parameter after extension, and T is measurement time Number, yijFor the variable x in higher-dimension random matrixijThe new variable obtained after the standardized way.
5. hollow coil current transformer error environment correlation analysis according to claim 1, which is characterized in that The step S4 specifically:
Calculate the singular value equivalent matrice D of the non-Hermite Matrixu
According to the singular value equivalent matrice DuCalculating matrix product Z;
Error state evaluating matrix Z is obtained according to the matrix product Z2
6. hollow coil current transformer error environment correlation analysis according to claim 5, which is characterized in that The matrix productThe error state evaluating matrixWherein,ziFor the row vector of matrix Z, zi= (zi1,zi2,...,ziT), 1≤i≤k',For matrix Z2Row vector, σ (zi) indicate zi's Standard deviation, k' are the number of the state parameter after extension, and T is pendulous frequency.
7. hollow coil current transformer error environment correlation analysis according to claim 1, which is characterized in that In step s 5, the relevance evaluation index includes dMSRAnd IMSR, wherein dMSRevref, dMSRIntegral to the time is IMSR:Wherein, t1And t2Indicate initial time and the finish time of assessment,λiFor the characteristic value of corresponding original matrix, λwiFor correspondence The characteristic value of R-matrix, n and n2The characteristic value number of respectively corresponding original matrix and R-matrix, E () indicate feature It is worth sample expectation, the R-matrixR-matrix is by air core coil error state matrix and Gauss white noise Sound matrix is constituted, wherein D1For environment parameter matrix;DNFor noise matrix, dimension is identical with environment parameter matrix, and element is The stochastic variable of standardized normal distribution is obeyed, the white Gaussian noise amplitude being superimposed in amplitude and matrix-expand is identical.
8. a kind of hollow coil current transformer error state monitors system, which is characterized in that including original matrix building module, Higher-dimension random matrix constructs module, standardization module, influence amount relevance evaluation matrix module and relevance evaluation mould Block;
The original matrix building module is used to acquire environment parameter data and hollow coil current transformer error information, It assesses in time window and constructs original matrix according to the environment parameter data and error information of acquisition;
The higher-dimension random matrix building module is used to be extended the original matrix based on Kalman filter, establishes high Tie up random matrix;
The standardization module makes it be converted to row vector equal for being standardized to the higher-dimension random matrix Value is the non-Hermite Matrix that 0, variance is 1;
The influence amount relevance evaluation matrix module is used to obtain influence amount relevance evaluation according to the non-Hermite Matrix Matrix;
The relevance evaluation module is used to obtain hollow coil current transformer according to the influence amount relevance evaluation matrix Error environment relevance evaluation index, and according to the influence amount relevance evaluation matrix and the relevance evaluation index to sky Correlation between wire-core coil current transformer error and environment parameter is assessed.
9. hollow coil current transformer error state according to claim 8 monitors system, which is characterized in that the original Beginning matrix constructs the environment parameter data constructing environment parameter matrix that module is used to pass through acquisitionWherein, element PijIt indicates that environment parameter can be surveyed in the measured value at j moment i-th, i is It can survey the serial number of environment parameter, i=1,2 ... ... M, M are the number of environment parameter, and j be the serial number measured, j=1,2 ... ... T, T is pendulous frequency;Error state matrix is constructed by the error information of acquisitionWherein, Element QijIndicating measured value of i-th of transformer error parameter at the j moment, i is the serial number of transformer error parameter, i=1, 2 ... ... N, N are the number of transformer error parameter, and j is the serial number of measurement, and the original matrix of j=1,2 ... ... T, building areWherein, k=M+N.
10. hollow coil current transformer error state according to claim 8 monitors system, which is characterized in that described Higher-dimension random matrix constructs the higher-dimension random matrix that module is establishedk' For extension after state parameter number, the value range of N' meet k'/T ∈ (0,1], T is pendulous frequency.
11. hollow coil current transformer error state according to claim 8 monitors system, which is characterized in that described The non-Hermite Matrix that standardization resume module obtains is that the non-Hermite Matrix isWherein Indicate sample xiBe averaged Value, σ (xij) indicate sample xiStandard deviation, xiFor higher-dimension random matrix D3Row vector, xi=(xi1,xi2,…,xiT),1≤i ≤ k', k' are the number of the state parameter after extension, and T is pendulous frequency, and yij is the variable x in higher-dimension random matrixijBy this The new variable obtained after standardized way.
12. hollow coil current transformer error state according to claim 8 monitors system, which is characterized in that described Influence amount relevance evaluation matrix module is used to calculate the singular value equivalent matrice D of the non-Hermite Matrixu, according to the surprise Different value equivalent matrice DuCalculating matrix product Z obtains error state evaluating matrix Z according to the matrix product Z2
The matrix productThe error state assessment MatrixWherein,ziFor the row vector of matrix Z, zi =(zi1,zi2,…,ziT), 1≤i≤k',For matrix Z2Row vector, σ (zi) indicate zi Standard deviation, k' be extension after state parameter number, T is pendulous frequency
13. hollow coil current transformer error state according to claim 8 monitors system, which is characterized in that described Relevance evaluation index includes dMSRAnd IMSR, wherein dMSRevref, dMSRIntegral to the time is Wherein, t1And t2Indicate initial time and the finish time of assessment,λi For the characteristic value of corresponding original matrix, λwiFor the characteristic value of corresponding R-matrix, n and n2Respectively corresponding original matrix and The characteristic value number of R-matrix, E () indicate the expectation of characteristic value sample, the R-matrixR-matrix It is made of air core coil error state matrix and white Gaussian noise matrix, wherein D1For environment parameter matrix;DNFor noise matrix, Its dimension is identical with environment parameter matrix, and element is the stochastic variable for obeying standardized normal distribution, folds in amplitude and matrix-expand The white Gaussian noise amplitude added is identical.
CN201910270802.7A 2019-04-04 2019-04-04 A kind of hollow coil current transformer error environment correlation analysis Pending CN110031788A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910270802.7A CN110031788A (en) 2019-04-04 2019-04-04 A kind of hollow coil current transformer error environment correlation analysis

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910270802.7A CN110031788A (en) 2019-04-04 2019-04-04 A kind of hollow coil current transformer error environment correlation analysis

Publications (1)

Publication Number Publication Date
CN110031788A true CN110031788A (en) 2019-07-19

Family

ID=67237458

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910270802.7A Pending CN110031788A (en) 2019-04-04 2019-04-04 A kind of hollow coil current transformer error environment correlation analysis

Country Status (1)

Country Link
CN (1) CN110031788A (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111474510A (en) * 2020-04-25 2020-07-31 华中科技大学 Error evaluation method and system for voltage transformer with non-stable output
CN112014785A (en) * 2020-08-06 2020-12-01 三峡大学 Error compensation method for air-core coil current transformer based on elastic network
CN112485747A (en) * 2020-11-28 2021-03-12 国网江西省电力有限公司供电服务管理中心 Error state evaluation method for capacitor voltage transformer
CN112505386A (en) * 2020-08-25 2021-03-16 中国电力科学研究院有限公司 Method and system for detecting current value of direct current charging pile
CN112713876A (en) * 2020-12-16 2021-04-27 南京大学 Kalman filter-based data processing method applicable to ACCT
CN113792256A (en) * 2021-07-27 2021-12-14 国网河北省电力有限公司营销服务中心 Low-voltage current transformer metering performance evaluation method and device and electronic equipment
CN115480203A (en) * 2022-09-28 2022-12-16 武汉格蓝若智能技术有限公司 Online quantitative evaluation method and system for error state of current transformer
CN116295620A (en) * 2023-02-17 2023-06-23 南通科瑞环境科技有限公司 Environment monitoring, collecting and detecting method
CN116754830A (en) * 2023-06-13 2023-09-15 武汉阿迪克电子股份有限公司 Metering method and device of electric energy meter and storage medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104062617A (en) * 2014-07-08 2014-09-24 国家电网公司 Merging-unit metrological-characteristic on-site detecting system and method thereof
CN106548021A (en) * 2016-11-03 2017-03-29 华中科技大学 A kind of electronic mutual inductor error state appraisal procedure based on higher-dimension random matrix
CN107133684A (en) * 2016-02-26 2017-09-05 中国电力科学研究院 A kind of random matrix construction method towards GA for reactive power optimization
CN108828405A (en) * 2018-06-06 2018-11-16 西南交通大学 A kind of electric transmission line fault detection method based on random matrix

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104062617A (en) * 2014-07-08 2014-09-24 国家电网公司 Merging-unit metrological-characteristic on-site detecting system and method thereof
CN107133684A (en) * 2016-02-26 2017-09-05 中国电力科学研究院 A kind of random matrix construction method towards GA for reactive power optimization
CN106548021A (en) * 2016-11-03 2017-03-29 华中科技大学 A kind of electronic mutual inductor error state appraisal procedure based on higher-dimension random matrix
CN108828405A (en) * 2018-06-06 2018-11-16 西南交通大学 A kind of electric transmission line fault detection method based on random matrix

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
胡琛: "空心线圈电流互感器可靠性及长期稳定性关键问题研究", 《万方 HTTP://D.WANFANGDATA.COM.CN/CONFERENCE/7064537》 *
胡琛等: "基于随机矩阵理论的电子式互感器误差状态相关性分析方法", 《电力自动化设备》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111474510A (en) * 2020-04-25 2020-07-31 华中科技大学 Error evaluation method and system for voltage transformer with non-stable output
CN112014785A (en) * 2020-08-06 2020-12-01 三峡大学 Error compensation method for air-core coil current transformer based on elastic network
CN112505386A (en) * 2020-08-25 2021-03-16 中国电力科学研究院有限公司 Method and system for detecting current value of direct current charging pile
CN112505386B (en) * 2020-08-25 2022-09-02 中国电力科学研究院有限公司 Method and system for detecting current value of direct current charging pile
CN112485747A (en) * 2020-11-28 2021-03-12 国网江西省电力有限公司供电服务管理中心 Error state evaluation method for capacitor voltage transformer
CN112713876A (en) * 2020-12-16 2021-04-27 南京大学 Kalman filter-based data processing method applicable to ACCT
CN112713876B (en) * 2020-12-16 2023-12-05 南京大学 Kalman filter-based data processing method suitable for ACCT
CN113792256A (en) * 2021-07-27 2021-12-14 国网河北省电力有限公司营销服务中心 Low-voltage current transformer metering performance evaluation method and device and electronic equipment
CN115480203A (en) * 2022-09-28 2022-12-16 武汉格蓝若智能技术有限公司 Online quantitative evaluation method and system for error state of current transformer
CN116295620A (en) * 2023-02-17 2023-06-23 南通科瑞环境科技有限公司 Environment monitoring, collecting and detecting method
CN116754830A (en) * 2023-06-13 2023-09-15 武汉阿迪克电子股份有限公司 Metering method and device of electric energy meter and storage medium
CN116754830B (en) * 2023-06-13 2024-01-09 武汉阿迪克电子股份有限公司 Metering method and device of electric energy meter and storage medium

Similar Documents

Publication Publication Date Title
CN110031788A (en) A kind of hollow coil current transformer error environment correlation analysis
CN103198437B (en) A kind of electrical network metric data and electric network model error correction method and device
CN110095744A (en) A kind of electronic mutual inductor error prediction method
CN106338706B (en) A kind of methods, devices and systems of electric energy metering device global error detection
CN106443285B (en) Multiple-harmonic-source harmonic responsibility quantitative analysis method based on total least square method
CN106295911B (en) A kind of grid branch parameter evaluation method based on chromatographic assays
CN105158725B (en) A kind of electric energy meter metering accuracy appraisal procedure based on multidimensional influence amount
CN106548021B (en) A kind of electronic mutual inductor error state appraisal procedure based on higher-dimension random matrix
CN106338708A (en) Electric energy metering error analysis method by combining deep learning and recursive neural network
CN101807795A (en) Method for forming electric energy metering simulation system and device thereof
CN112485747A (en) Error state evaluation method for capacitor voltage transformer
CN107741577A (en) A kind of Source of Gateway Meter degree of accuracy on-line monitoring and analysis method and system
CN109307854A (en) A kind of electric energy metering device comprehensive estimation method and system
CN103324858A (en) Three-phase load flow state estimation method of power distribution network
Sheng et al. Online measurement-based estimation of dynamic system state matrix in ambient conditions
CN113743651A (en) Electronic voltage transformer error state prediction method based on machine learning and SVM algorithm
CN113837432A (en) Power system frequency prediction method driven by physics-data combination
CN109581271A (en) A kind of typical low pressure platform area electricity consumption data rapid simulation method
CN105068035A (en) Voltage transformer error level dynamic detection method and system
CN111999691B (en) Error calibration method and error calibration device for metering sensor device
CN106022970B (en) A kind of active distribution network measure configuration method counted and distributed generation resource influences
CN105426656B (en) A kind of error calculation method and system
CN109359822A (en) Electronic type voltage transformer measuring state appraisal procedure and system
CN107132500B (en) A kind of synchronous phasor measurement unit on-line calibration method and apparatus
CN109982365A (en) Antenna feeder problem check method and device based on emulation and MRO data

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20190719