CN110967661A - Electrical data calibration method based on curve fitting - Google Patents

Electrical data calibration method based on curve fitting Download PDF

Info

Publication number
CN110967661A
CN110967661A CN201911327267.0A CN201911327267A CN110967661A CN 110967661 A CN110967661 A CN 110967661A CN 201911327267 A CN201911327267 A CN 201911327267A CN 110967661 A CN110967661 A CN 110967661A
Authority
CN
China
Prior art keywords
value
model
error
data
electrical
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
CN201911327267.0A
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.)
Ningxia Kaichen Electric Group Co ltd
Original Assignee
Ningxia Kaichen Electric Group Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Ningxia Kaichen Electric Group Co ltd filed Critical Ningxia Kaichen Electric Group Co ltd
Priority to CN201911327267.0A priority Critical patent/CN110967661A/en
Publication of CN110967661A publication Critical patent/CN110967661A/en
Pending legal-status Critical Current

Links

Images

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)
  • Measurement Of Current Or Voltage (AREA)

Abstract

The invention relates to an electrical data calibration method based on curve fitting. The method is characterized by comprising the following steps: (1) collecting a measured value and an actual value of the electrical measurement data, and calculating an error value; (2) performing visualization processing on the obtained measurement value and the error value, and defining an x axis as the measurement value and a y axis as the error value; (3) selecting a proper model according to the obtained image trend; (4) calculating various parameters of the model according to the error minimum principle; (5) calculating a residual error according to the original data brought by the model, analyzing whether the residual error meets the requirement, and returning to the step (3) to reselect a higher-order model if the residual error does not meet the requirement; (6) and if the residual error meets the requirement, calibrating the measured value according to the model, and taking the original measured value of the model as the input quantity of the model, wherein the output quantity is the calibrated value. The method can solve the precision problem caused by the nonlinear error caused by the difference of sensor hardware.

Description

Electrical data calibration method based on curve fitting
Technical Field
The invention relates to an electrical data calibration method based on curve fitting.
Background
High-voltage complete equipment (high-voltage distribution cabinet) refers to indoor and outdoor alternating-current switch equipment which operates in a power system with voltage of 3kV or above and frequency of 50Hz or below. The power system protection device is mainly used for controlling and protecting power systems (including users of power plants, substations, power transmission and distribution lines, industrial and mining enterprises and the like), can be used for not only putting a part of power equipment or lines into or out of operation according to the operation requirements of a power grid, but also quickly cutting off the fault part from the power grid when the power equipment or lines have faults, thereby ensuring the normal operation of the fault-free part in the power grid and the safety of equipment and operation maintenance personnel. In the whole link, voltage and current are required to be collected from high-voltage electricity (3KV and above) for judgment, the voltage and the current on the high-voltage side are called primary data, the data is subjected to transformation ratio through a mutual inductor, the transformation ratio of the current (CT transformation ratio) is generally 2000 times, and the transformation ratio of the voltage (PT transformation ratio) is generally 100 times. The transformed data is referred to as secondary data.
The error value of the secondary data is amplified by the transformation ratio of the transformer, so that the precision of the secondary value is particularly important, and the comprehensive protection algorithm and the measurement algorithm are all based on the secondary value data, so the quality and the safety of the electrical equipment and the comprehensive protection equipment are influenced by the precision of the electrical data, and the quality of the power grid is also concerned.
Disclosure of Invention
The invention aims to provide an electrical data calibration method based on curve fitting, which can improve the accuracy of secondary values of electrical data as much as possible.
A method for calibrating electrical data based on curve fitting is characterized by comprising the following steps:
(1) collecting a measured value and an actual value of the electrical measurement data, and calculating an error value;
(2) performing visualization processing on the obtained measurement value and the error value, and defining an x axis as the measurement value and a y axis as the error value;
(3) selecting a proper model according to the obtained image trend;
(4) calculating various parameters of the model according to the error minimum principle;
(5) calculating a residual error according to the original data brought by the model, analyzing whether the residual error meets the requirement, and returning to the step (3) to reselect a higher-order model if the residual error does not meet the requirement;
(6) and if the residual error meets the requirement, calibrating the measured value according to the model, taking the original measured value of the model as the input quantity of the model, wherein the output quantity is the calibrated value, and ending the whole calibration process of the electrical data.
The electrical measurement data is specifically current, voltage, electrical power, electrical energy, phase, frequency, resistance, capacitance or inductance.
The method has the beneficial effects that: 1. the accuracy problem caused by non-linear errors due to differences in sensor hardware can be solved. 2. The accuracy of the secondary values of the electrical data can be improved as much as possible.
Drawings
FIG. 1 is a graph of the error of the A-phase voltage versus the measured voltage;
FIG. 2 is a fitting curve and an error curve of the A-phase voltage;
FIG. 3 is a graph of the fitting of the A-phase voltage and the residual after the fitting of the A-phase voltage.
Detailed Description
The invention relates to a solution for electric data requiring high precision, which is based on the minimum error theory, and different models can be selected to match the error of the whole precision, thereby establishing a set of voltage and current precision fitting compensation method system suitable for the electric industry. The device is suitable for microcomputer comprehensive protection and enterprise units needing high-precision current and voltage measurement.
The invention aims to provide a method for compensating the accuracy of the original value of the electrical data, which is mainly based on the two ideas that the sum of all errors of the created model and the electrical data value is minimum or the sum of the squares of the errors is minimum. The sample regression function has the smallest total error of the values of the explained variables from the known sample data point values, i.e. in the coordinate graph the perpendicular distance of the points of all the explained variables on the fitting function from the known sample data points is smallest with the same value of the explained variable.
The principle of the method is as follows:
for a given set of data points { (x)i,yi): if the curve model is fitted, i is 1, 2, 3, …, n
Figure BDA0002328696630000031
The error distance at the i-th position is
Figure BDA0002328696630000032
Then the fitted curve model value F (x) at the i-th positioni) With the actual given data point YiThe sum of the squares of the differences is
Figure BDA0002328696630000033
Further, find out
Figure BDA0002328696630000034
To obtain a fitted curve
Figure BDA0002328696630000035
The method firstly determines a model for fitting a curve, such as a linear model, a quadratic model, a cubic model, a fourth-order polynomial, a fifth-order polynomial and the like. After the model is determined, the parameters of the selected model are determined according to the concept of minimum error. After obtaining parameters of the model, the original data to be fitted is taken as a variable of the model to be input, and the obtained output is the data after fitting.
Regarding the principle of linear fitting, for a given data point (x)i,yi): 1, 2, 3, …, n is obtained as a linear equation of y ═ a + bx, with the total error being calculated
Figure BDA0002328696630000036
To minimize, the parameters a, b that make Q extremely good should be satisfied
Figure BDA0002328696630000037
This is true.
Figure BDA0002328696630000038
The parameters of the linear model are then found from the set of equations above.
With respect to the principle of polynomial fitting, sometimes it is not appropriate to fit a given data point with a straight line, and the use of polynomial fitting is considered. For a given data point (x)i,yi): 1, 2, 3, …, n is obtained as a polynomial of degree m (m < n)
Figure BDA0002328696630000042
Make the total error
Figure BDA0002328696630000043
Is minimal. Since Q can be considered as relating to aj(j is 0, 1, …, m), so the problem of fitting a polynomial can be reduced to the problem of the extremum of the multivariate function.
Figure BDA0002328696630000044
To obtain
Figure BDA0002328696630000045
Namely, the coefficient a is obtainedjThe system of linear equations is commonly referred to as a system of regular equations.
Example 1:
in a first step, measured and actual values of the electrical data are collected and an error value is calculated.
The first step of the method is to collect electrical data meeting the requirements, and the mode of collecting the electrical data can refer to the following scheme that high voltage and high current, namely voltage and current of a primary side are transmitted to a secondary transformer after being transformed by the transformer, 3kv and above are converted into 110V, 10kA is converted into a range of 5A, the alternating current data is sampled by a high-precision AD sampler, the effective value of the sampled data is calculated by a fourier transform or root mean square method, and a data interval required to be subjected to precision compensation is subjected to value taking, wherein the smaller the value step length of the interval is, the higher the compensation precision is. The error values obtained by subtracting the measured values from the actual values are shown in the following table.
Figure BDA0002328696630000041
Figure BDA0002328696630000051
And secondly, performing visualization processing on the data: the data is visualized, the horizontal axis (x) is measurement data, the vertical axis (y) is an error value, and x (after compensation) is x + y; and f (x), we need to find f (x) in the following, see fig. 1 for details.
Thirdly, selecting a proper function model and calculating model parameters:
an appropriate functional model is selected based on the error and the graphical trend of the measurement data, where we select a linear model, and we obtain the parameters of the model based on the principle of error minimization and based on the inventive principles described above. The linear model is y ═ 0.089 x-0.18. This is visualized as shown in the following figure, which is described in detail in figure 2.
Fourthly, calculating a residual error, and judging whether the requirement is met:
and calculating a residual error between the compensation curve and the error curve, wherein the residual error is the error fluctuation of the compensation curve, and the error fluctuation after compensation is basically within 0.2V, if the error fluctuation does not meet the requirement, performing high-order model compensation, and particularly referring to an attached figure 3.

Claims (2)

1. A method for calibrating electrical data based on curve fitting is characterized by comprising the following steps:
(1) collecting a measured value and an actual value of the electrical measurement data, and calculating an error value;
(2) performing visualization processing on the obtained measurement value and the error value, and defining an x axis as the measurement value and a y axis as the error value;
(3) selecting a proper model according to the obtained image trend;
(4) calculating various parameters of the model according to the error minimum principle;
(5) calculating a residual error according to the original data brought by the model, analyzing whether the residual error meets the requirement, and returning to the step (3) to reselect a higher-order model if the residual error does not meet the requirement;
(6) and if the residual error meets the requirement, calibrating the measured value according to the model, taking the original measured value of the model as the input quantity of the model, wherein the output quantity is the calibrated value, and ending the whole calibration process of the electrical data.
2. A method of electrical data calibration based on curve fitting according to claim 1, wherein: the electrical measurement data is specifically current, voltage, electrical power, electrical energy, phase, frequency, resistance, capacitance or inductance.
CN201911327267.0A 2019-12-20 2019-12-20 Electrical data calibration method based on curve fitting Pending CN110967661A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911327267.0A CN110967661A (en) 2019-12-20 2019-12-20 Electrical data calibration method based on curve fitting

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911327267.0A CN110967661A (en) 2019-12-20 2019-12-20 Electrical data calibration method based on curve fitting

Publications (1)

Publication Number Publication Date
CN110967661A true CN110967661A (en) 2020-04-07

Family

ID=70035715

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911327267.0A Pending CN110967661A (en) 2019-12-20 2019-12-20 Electrical data calibration method based on curve fitting

Country Status (1)

Country Link
CN (1) CN110967661A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112953538A (en) * 2021-02-23 2021-06-11 青岛鼎信通讯股份有限公司 ADC calibration method applied to direct current meter calibrating device
CN114114131A (en) * 2021-11-09 2022-03-01 航天亮丽电气有限责任公司 Method for calibrating metering precision of electric energy meter
CN114778995A (en) * 2022-06-22 2022-07-22 南京亿高微波系统工程有限公司 High-frequency electrotome precision automatic measurement method and device
CN115993567A (en) * 2023-03-22 2023-04-21 深圳市北汉科技有限公司 Calibration method and system for bidirectional feedback power supply data

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040001277A1 (en) * 2002-06-27 2004-01-01 International Business Machines Corporation Apparatus and method to calibrate servo sensors in a noisy environment
JP2013072643A (en) * 2011-09-26 2013-04-22 Fujikura Ltd Magnetic measurement data calibration device and azimuth angle measurement device
CN104360303A (en) * 2014-11-04 2015-02-18 江阴长仪集团有限公司 Method for verifying self-calibration error along with voltage fluctuation on basis of intelligent ammeters
CN107193782A (en) * 2017-04-18 2017-09-22 中国船舶重工集团公司第七〇九研究所 A kind of method of abnormal value removing and correction fitted based on multinomial
CN107273336A (en) * 2017-05-10 2017-10-20 北京航天控制仪器研究所 A kind of curve-fitting method based on least square method
CN110242589A (en) * 2019-06-25 2019-09-17 江苏大学 A kind of centrifugal pump performance fitting modification method
CN110568254A (en) * 2019-09-27 2019-12-13 宁夏凯晨电气集团有限公司 Method for accurately detecting attenuated direct-current component parameters in fault current

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040001277A1 (en) * 2002-06-27 2004-01-01 International Business Machines Corporation Apparatus and method to calibrate servo sensors in a noisy environment
JP2013072643A (en) * 2011-09-26 2013-04-22 Fujikura Ltd Magnetic measurement data calibration device and azimuth angle measurement device
CN104360303A (en) * 2014-11-04 2015-02-18 江阴长仪集团有限公司 Method for verifying self-calibration error along with voltage fluctuation on basis of intelligent ammeters
CN107193782A (en) * 2017-04-18 2017-09-22 中国船舶重工集团公司第七〇九研究所 A kind of method of abnormal value removing and correction fitted based on multinomial
CN107273336A (en) * 2017-05-10 2017-10-20 北京航天控制仪器研究所 A kind of curve-fitting method based on least square method
CN110242589A (en) * 2019-06-25 2019-09-17 江苏大学 A kind of centrifugal pump performance fitting modification method
CN110568254A (en) * 2019-09-27 2019-12-13 宁夏凯晨电气集团有限公司 Method for accurately detecting attenuated direct-current component parameters in fault current

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
B.W. SILVERMAN: "Some Aspects of the Spline Smoothing Approach to Non-parametric Regression Curve Fitting", 《JOURNAL OF THE ROYAL STATISTICAL SOCIETY: SERIES B (METHODOLOGICAL)》 *
BERNARD B.T.KEE等: "Aresidual based error estimator using radial basis functions", 《FINITE ELEMENTS IN ANALYSIS AND DESIGN》 *
岳鹏等: "用MatLab曲线拟合工具箱处理物理实验数据", 《大学物理实验》 *
陈岚峰: "基于MATLAB的最小二乘曲线拟合仿真研究", 《沈阳师范大学学报(自然科学版)》 *
高志钰等: "基于曲线拟合的AR模型在地铁施工监测中的应用", 《测绘工程》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112953538A (en) * 2021-02-23 2021-06-11 青岛鼎信通讯股份有限公司 ADC calibration method applied to direct current meter calibrating device
CN114114131A (en) * 2021-11-09 2022-03-01 航天亮丽电气有限责任公司 Method for calibrating metering precision of electric energy meter
CN114778995A (en) * 2022-06-22 2022-07-22 南京亿高微波系统工程有限公司 High-frequency electrotome precision automatic measurement method and device
CN114778995B (en) * 2022-06-22 2022-09-06 南京亿高微波系统工程有限公司 High-frequency electrotome precision automatic measurement method and device
CN115993567A (en) * 2023-03-22 2023-04-21 深圳市北汉科技有限公司 Calibration method and system for bidirectional feedback power supply data
CN115993567B (en) * 2023-03-22 2023-05-19 深圳市北汉科技有限公司 Calibration method and system for bidirectional feedback power supply data

Similar Documents

Publication Publication Date Title
CN110967661A (en) Electrical data calibration method based on curve fitting
Moutis et al. Digital twin of distribution power transformer for real-time monitoring of medium voltage from low voltage measurements
CN109031000B (en) A kind of method and system based on non-faulting disturbance In situ Measurement grid short circuit capacity
CN103869183A (en) Online monitoring system and method of zinc oxide arrester
CN109490813A (en) A kind of current transformer characteristic appraisal procedure and system
CN103809010B (en) For determining the voltage of every phase and the device and method of power of middle-voltage network
KR20110015964A (en) Distribution automation system and its voltage control method for reactive power compensation
Špelko et al. CIGRE/CIRED JWG C4. 42: Overview of common methods for assessment of harmonic contribution from customer installation
CN106054023B (en) The method of two sides system impedance is estimated in a kind of transmission line of electricity single end distance measurement
US10901026B2 (en) Passive harmonic filter power quality monitor and communications device
CN112285634B (en) Method for identifying abnormal measurement error of high-voltage transformer in double-bus structure
CN107831378B (en) Device and method for detecting compensation effect of arc suppression coil
Leinakse et al. Impact of distributed generation on estimation of exponential load models
CN105678640B (en) Method for predicting direct current distribution of alternating current power grid in consideration of transformer saturation influence
CN105403808B (en) A kind of localization method and device of DC line earth fault
Blajszczak Resistive voltage divider for higher harmonics measurement in 400 kV network
CN204359923U (en) A kind of novel high voltage direct current generator calibration equipment
KR102027844B1 (en) Apparatus for correcting measurement error of measuring instrument using equipment information and method thereof
CN107171327A (en) A kind of Power Network Status Estimation method and apparatus
Yamashita et al. A development of dynamic load model parameter derivation method
Bracale et al. DC electric arc furnace modelling for power quality indices assessment
CN115372728A (en) Current transformer working state monitoring method and device and storage medium
Franković et al. Application of different least square methods for transmission line parameter estimation
CN109342817B (en) A kind of non-fully commutation three-phase transmission line method for parameter estimation measured based on PMU
Chiulan et al. A practical example of power transformer unit winding condition assessment by means of short-circuit impedance measurement

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: 20200407