CN115508603A - Tower structure parameter inversion method adopting space electric field fixed point integral - Google Patents

Tower structure parameter inversion method adopting space electric field fixed point integral Download PDF

Info

Publication number
CN115508603A
CN115508603A CN202211227996.0A CN202211227996A CN115508603A CN 115508603 A CN115508603 A CN 115508603A CN 202211227996 A CN202211227996 A CN 202211227996A CN 115508603 A CN115508603 A CN 115508603A
Authority
CN
China
Prior art keywords
tower
electric field
insulator
integral
fixed point
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
CN202211227996.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.)
Guizhou Power Grid Co Ltd
Original Assignee
Guizhou Power Grid 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 Guizhou Power Grid Co Ltd filed Critical Guizhou Power Grid Co Ltd
Priority to CN202211227996.0A priority Critical patent/CN115508603A/en
Publication of CN115508603A publication Critical patent/CN115508603A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R19/00Arrangements for measuring currents or voltages or for indicating presence or sign thereof
    • G01R19/0084Arrangements for measuring currents or voltages or for indicating presence or sign thereof measuring voltage only
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/08Measuring electromagnetic field characteristics
    • G01R29/0807Measuring electromagnetic field characteristics characterised by the application
    • G01R29/0814Field measurements related to measuring influence on or from apparatus, components or humans, e.g. in ESD, EMI, EMC, EMP testing, measuring radiation leakage; detecting presence of micro- or radiowave emitters; dosimetry; testing shielding; measurements related to lightning
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/08Measuring electromagnetic field characteristics
    • G01R29/0864Measuring electromagnetic field characteristics characterised by constructional or functional features
    • G01R29/0892Details related to signal analysis or treatment; presenting results, e.g. displays; measuring specific signal features other than field strength, e.g. polarisation, field modes, phase, envelope, maximum value
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Electromagnetism (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Insulators (AREA)

Abstract

The invention relates to the technical field of transmission tower measurement data processing, in particular to a tower structure parameter inversion method adopting space electric field fixed point integral, which comprises the following steps: judging the structure type of the tower insulator, if the structure type of the tower insulator is judged to be the first type, determining an integral path of the tower insulator, selecting a first set number of sensors to be arranged on the tower insulator, and executing the next step; if the structure type of the tower insulator is judged to be the second type, determining an integral path of the tower insulator, selecting a second set number of sensors to be arranged on the tower insulator, and executing the next step; and acquiring electric field change data of the tower, performing inversion calculation based on the structure type of the tower to obtain the voltage value of the tower, and completing the inversion calculation of the voltage value of the transmission tower.

Description

Fixed point integration by adopting space electric field Pole tower structure parameter inversion method
Technical Field
The invention relates to the technical field of transmission tower measurement data processing, in particular to a tower structure parameter inversion method adopting space electric field fixed point integral.
Background
For the moment, voltage measurement is the most important and fundamental measurement content in grid sensing measurement. At the present stage, the monitoring of the electric energy state of the transmission and distribution side is mainly concentrated in a transformer substation, and in contrast, the difficulty of monitoring the voltage of a transmission tower is very high, and measurement and calculation are usually performed under the condition of contacting with the high voltage at the primary side of a transmission line, so that the installation difficulty of a sensor is high, and the installation safety of the sensor cannot be guaranteed; in addition, due to the fact that the tower structure is complex, the problem of solving the over-determined equation set is involved in the inverse problem calculation process of reversely calculating the field source parameters after electric field signals around the complex electromagnetic field are obtained, and therefore the problem of solving the inverse problem of the lead voltage is very difficult. Based on the above, in order to solve the above problems, a pole tower structure parameter inversion method adopting space electric field fixed point integration is designed.
Disclosure of Invention
The invention aims to provide a tower structure parameter inversion method adopting space electric field fixed point integral, which is used for solving the technical problem.
The embodiment of the invention is realized by the following technical scheme:
a tower structure parameter inversion method adopting space electric field fixed point integral comprises the following steps:
judging the structure type of the tower insulator, if the structure type of the tower insulator is judged to be the first type, determining an integral path of the tower insulator, selecting a first set number of sensors to be arranged on the tower insulator, and executing the next step; if the structure type of the tower insulator is judged to be the second type, determining an integral path of the tower insulator, selecting a second set number of sensors to be arranged on the tower insulator, and executing the next step;
and acquiring electric field change data of the tower, performing inversion calculation based on the structure type of the tower to obtain the voltage value of the tower, and completing the inversion calculation of the voltage value of the transmission tower.
Optionally, the tower insulator has two specific structure types.
Optionally, the structure types of the pole tower insulator include a glass insulator and a silica gel insulator, the glass insulator is the first structure type, and the silica gel insulator is the second structure type.
Optionally, the inversion calculation specifically uses a gaussian-legendre integral.
Optionally, the sensor is specifically an electric field sensor.
The technical scheme of the embodiment of the invention at least has the following advantages and beneficial effects:
according to the embodiment, the structure of the transmission tower is subjected to inversion calculation, and the voltage value of the transmission tower can be measured through inversion calculation of the structure of the tower by measuring the change of an electric field at one side of the structure of the transmission tower without contacting the high voltage at the primary side of the transmission line, so that the installation process and the difficulty of a tower sensor are greatly simplified, and the installation safety of the sensor is improved; moreover, inversion calculation is carried out through Gaussian-Legendre integration, calculation is faster compared with the inverse problem, and meanwhile the requirement of measurement accuracy of the voltage of the transmission tower can be met under most conditions.
Drawings
Fig. 1 is a schematic flow diagram of a pole tower structure parameter inversion method using space electric field fixed point integration according to the present invention;
FIG. 2 is a schematic view of an integration path of a glass insulator according to the present invention;
fig. 3 is a schematic diagram of an integral path of a silica gel insulator according to the present invention;
fig. 4 is a schematic diagram of the electric field distribution after the irregular part is cut off.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. The components of embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations.
As shown in fig. 1, the present invention provides one of the embodiments: a tower structure parameter inversion method adopting space electric field fixed point integral comprises the following steps:
judging the structure type of a pole tower insulator, selecting Gaussian Legendre integral under a 4-sensor as an integral scheme when the pole tower insulator is a glass insulator and the integral path of the pole tower glass insulator is shown in figure 2, and arranging an electric field sensor on the pole tower insulator to perform inversion by adopting the Gaussian Legendre integral; when the insulator of the tower is a silica gel insulator, see fig. 3, and the silica gel insulator is selected as an integral path in the tension tower, the Gaussian Legendre integral under the 5 sensor is selected as an integral scheme, and the electric field sensor is arranged on the insulator of the tower to perform inversion by adopting the Gaussian Legendre integral.
In this embodiment, the fundamental principle of the electric field inversion voltage is shown as follows:
Figure BDA0003880800720000041
wherein, U bc The potential difference between the point b and the point a is shown, and the minus sign before integration shows that the direction of the electric field intensity is opposite to the direction of the voltage rise. E represents the value of the electric field along the integration path. The superposition of the integral equation approximately equal to N discrete points is the description of integral physical meaning, and the larger N is, the smaller the error after the dispersion is. With this inversion scheme, only a certain number of sensors need to be placed on the electric field integration path.
Based on the numerical integration scheme of the above integration method inversion, the gaussian-legendre integration is adopted in the present embodiment.
In the present embodiment, the integration scheme is a gaussian-legendre integral, i.e., a gaussian integral variant, and is a special gaussian integral in the case where the weight function in the gaussian integral is set to ρ (x) =1 and the integration interval is changed to [ -1,1], such as the following formula, where the weight and position are determined using a legendre polynomial having an orthogonal property.
Figure BDA0003880800720000042
Wherein the legendre polynomial for solving the weights and the integral node locations in the gaussian-legendre integral is shown as follows:
Figure BDA0003880800720000043
when the integral constants under the N integral nodes are solved, the zero points corresponding to the N-order Legendre polynomial are corresponding integral nodes. Since the gaussian legendre integral range can only be [ -1,1], and the actual electric field integral range is [ a, b ], the original integral object dx needs to be converted, and the conversion formula is as follows:
Figure BDA0003880800720000051
and putting the Gaussian Legendre integral formula into an actual integral formula of the electric field solving voltage to obtain the Gaussian Legendre integral formula for the electric field inversion voltage as follows:
Figure BDA0003880800720000052
the actual gaussian-legendre polynomial integral node locations and weights can be found as shown in table 1. In practice, when the electric field inversion voltage is performed, the number of integration nodes is determined, and then the corresponding electric field sensor placement position is determined according to table 1 for inversion.
TABLE 1
Figure BDA0003880800720000061
As shown in fig. 4, this embodiment further provides an application example, in which a voltage-electric field distribution curve on the silicone rubber insulator is obtained through calculation, and irregular parts in the voltage-electric field distribution curve are removed, as shown in the following diagram. After the irregular part is deleted, the voltage difference between the two ends of the integration path is 666090V, the voltage is used as an inversion reference voltage value, and the integration length is 14385mm.
In this application example, based on comparison, chebyshev integral inversion was performed on the silicone rubber insulator integration path in the case of 3 sensors, 4 sensors, and 5 sensors, respectively, as shown in tables 2, 3, and 4.
TABLE 2
Figure BDA0003880800720000071
TABLE 3
Figure BDA0003880800720000072
TABLE 4
Figure BDA0003880800720000073
The accuracy is obviously improved under the condition that the number of sensors is increased for the Chebyshev integral on the silicon rubber insulator. Further, in the case of sensors 3, 4, and 5, the gaussian legendre integrals were performed as shown in tables 5, 6, and 7.
TABLE 5
Figure BDA0003880800720000081
TABLE 6
Figure BDA0003880800720000082
TABLE 7
Figure BDA0003880800720000083
It can be analyzed from the above table that, under the condition that the silicone rubber insulator is used as an integral path for inversion, the gaussian legendre integral precision is higher than the chebyshev integral precision. Therefore, when the tension tower selects a silicon rubber insulator as an integral path, the Gaussian Legendre integral under a sensor of 5 is selected as an optimal integral scheme, and the effectiveness of the invention is proved.
The present invention has been described in terms of the preferred embodiment, and it is not intended to be limited to the embodiment. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (5)

1. A tower structure parameter inversion method adopting space electric field fixed point integral is characterized by comprising the following steps:
judging the structure type of the pole tower insulator, if the structure type of the pole tower insulator is judged to be the first type, determining an integral path of the pole tower insulator, selecting a first set number of sensors to be arranged on the pole tower insulator, and executing the next step; if the structure type of the tower insulator is judged to be the second type, determining an integral path of the tower insulator, selecting a second set number of sensors to be arranged on the tower insulator, and executing the next step;
and acquiring electric field change data of the tower, performing inversion calculation based on the structure type of the tower to obtain the voltage value of the tower, and completing the inversion calculation of the voltage value of the transmission tower.
2. The pole tower structure parameter inversion method adopting space electric field fixed point integration according to claim 1, wherein two types of pole tower insulators are specifically adopted.
3. The pole tower structure parameter inversion method adopting space electric field fixed point integral as claimed in claim 2, wherein the structure types of the pole tower insulators comprise glass insulators and silica gel insulators, the glass insulators are of the first structure type, and the silica gel insulators are of the second structure type.
4. The tower structure parameter inversion method adopting space electric field fixed point integration according to claim 1, characterized in that the inversion calculation specifically adopts Gaussian-Legendre integration.
5. The pole tower structure parameter inversion method adopting space electric field fixed point integration according to claim 1, wherein the sensor is specifically an electric field sensor.
CN202211227996.0A 2022-10-09 2022-10-09 Tower structure parameter inversion method adopting space electric field fixed point integral Pending CN115508603A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211227996.0A CN115508603A (en) 2022-10-09 2022-10-09 Tower structure parameter inversion method adopting space electric field fixed point integral

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211227996.0A CN115508603A (en) 2022-10-09 2022-10-09 Tower structure parameter inversion method adopting space electric field fixed point integral

Publications (1)

Publication Number Publication Date
CN115508603A true CN115508603A (en) 2022-12-23

Family

ID=84509036

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211227996.0A Pending CN115508603A (en) 2022-10-09 2022-10-09 Tower structure parameter inversion method adopting space electric field fixed point integral

Country Status (1)

Country Link
CN (1) CN115508603A (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102103650A (en) * 2011-03-25 2011-06-22 江苏南大先腾信息产业有限公司 Three-dimensional simulation arrangement-along construction method for power transmission line
CN102694352A (en) * 2012-06-07 2012-09-26 甘肃省电力公司电力科学研究院 Method for improving insulator potential distribution in overhead power transmission line
CN104361168A (en) * 2014-11-12 2015-02-18 国家电网公司 Method for simulating electric field of 500kV alternating-current straight line tower
CN105425054A (en) * 2015-12-24 2016-03-23 国网重庆市电力公司电力科学研究院 Noncontact potential measurement method and device
CN107543953A (en) * 2017-09-18 2018-01-05 重庆大学 The phases line voltage detection method of transmission line of electricity based on modified Gauss Legendre integration
CN108387772A (en) * 2018-03-08 2018-08-10 清华大学 A kind of measurement method of transmission line of electricity overvoltage
CN112285495A (en) * 2020-09-21 2021-01-29 国网辽宁省电力有限公司营口供电公司 Electric field distribution based method for judging deterioration of insulator of power transmission line

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102103650A (en) * 2011-03-25 2011-06-22 江苏南大先腾信息产业有限公司 Three-dimensional simulation arrangement-along construction method for power transmission line
CN102694352A (en) * 2012-06-07 2012-09-26 甘肃省电力公司电力科学研究院 Method for improving insulator potential distribution in overhead power transmission line
CN104361168A (en) * 2014-11-12 2015-02-18 国家电网公司 Method for simulating electric field of 500kV alternating-current straight line tower
CN105425054A (en) * 2015-12-24 2016-03-23 国网重庆市电力公司电力科学研究院 Noncontact potential measurement method and device
CN107543953A (en) * 2017-09-18 2018-01-05 重庆大学 The phases line voltage detection method of transmission line of electricity based on modified Gauss Legendre integration
CN108387772A (en) * 2018-03-08 2018-08-10 清华大学 A kind of measurement method of transmission line of electricity overvoltage
CN112285495A (en) * 2020-09-21 2021-01-29 国网辽宁省电力有限公司营口供电公司 Electric field distribution based method for judging deterioration of insulator of power transmission line

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
周仿荣: "不同类型绝缘子在鸟粪影响下电场分布的仿真分析及防范建议", Retrieved from the Internet <URL:https://wwww.doc88.com/p-01899558913557.html> *

Similar Documents

Publication Publication Date Title
CN107589351B (en) Low and zero value insulator detection method for insulator detection robot
CN107561397B (en) Method and system for detecting electrical contact state of GIS contact system
CN203133171U (en) Testing device for testing wafer of electronic circuit
CN108108530B (en) Fatigue life calibration method suitable for structural connecting piece
CN110907881A (en) Current calibration method and system for power distribution network overhead line monitoring device
CN111521915A (en) High-voltage direct-current line corona onset field strength determination method and system
CN114740303B (en) Fault monitoring system of wireless passive high-voltage switch cabinet
CN116070068B (en) Stay cable damage identification method, device and terminal based on wavelet transformation of primary derivative of girder deflection difference
CN108548997B (en) Transformer substation space partial discharge positioning method and system
CN115508603A (en) Tower structure parameter inversion method adopting space electric field fixed point integral
CN112180177B (en) Power frequency electromagnetic field evaluation method and system fusing measured data
CN113218315A (en) Thickness measuring method, device and system
CN117074782A (en) Direct-current resistance measurement method and system for wires and cables
CN115640732B (en) Power distribution network arc fault positioning method based on magnetic field distribution
CN112784218B (en) Kurtosis estimation method for non-Gaussian random vibration signal
KR100316314B1 (en) Amorphous Region Determination Method and Apparatus
CN115524527A (en) Method for inverting voltage of internal through-flow conductor through cable accessory space electric field
CN116127352A (en) Switch load simulation test method based on data model
CN107729593B (en) Method for correcting measured voltage value of insulator detection robot based on support vector machine
CN108931808B (en) Improved ion beam vertical angle measuring method
CN116294947A (en) Perpendicular double-axis displacement measurement device, system and method based on magnetic induction
CN108763830A (en) The flicker noise model and its extracting method of semiconductor devices
CN111337732B (en) Voltage measurement method based on electric field inversion
CN113064023B (en) Fault identification method and system for same-tower multi-circuit power transmission line
CN111177641B (en) Wire temperature rise assessment method under high wind speed

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