CN110738009B - Setting method of analog charges in conducting wire in electric field calculation of power transmission line - Google Patents
Setting method of analog charges in conducting wire in electric field calculation of power transmission line Download PDFInfo
- Publication number
- CN110738009B CN110738009B CN201910975391.1A CN201910975391A CN110738009B CN 110738009 B CN110738009 B CN 110738009B CN 201910975391 A CN201910975391 A CN 201910975391A CN 110738009 B CN110738009 B CN 110738009B
- Authority
- CN
- China
- Prior art keywords
- charge
- analog
- electric field
- potential
- simulated
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Computing Systems (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Molecular Biology (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
The invention discloses a setting method of analog charges in a wire in electric field calculation of a power transmission line, which comprises the following steps: a. uniformly arranging a plurality of analog line charges on a certain concentric circle of each equivalent conducting wire, and arranging corresponding matching points and check points on the surface of the conducting wire; b. calculating an analog charge matrix; c. calculating the potential error of the check point; d. and obtaining the radius of the simulated charge distribution circle which minimizes the potential error of the check point by using a particle swarm optimization simulated charge method (PSOCSM), calculating the coordinates and the charge quantity of the corresponding simulated charge group, and finally calculating the electric field around the transmission line by using the simulated charge group. The invention can effectively realize the optimal setting of the analog charge in the lead, and compared with the conventional setting method (a central single method) of the analog charge in the lead, the electric potential error of the check point is greatly reduced, and the electric field calculation precision is improved.
Description
Technical Field
The invention belongs to the field of electric field calculation, and particularly relates to a method for setting simulated charges in a wire in electric field calculation of a power transmission line.
Background
With the rapid development of the economy in China, various high-voltage transmission projects are successively constructed, so that the electromagnetic environment problem around the transmission line is also brought into close attention. The power frequency electric field generated by the high-voltage transmission line is a main factor in the electromagnetic environment problem. Compared with field measurement, the electric field numerical calculation can save a great deal of manpower, material resources and time, and the electric field distribution of the whole field can be obtained through calculation. Therefore, the method has great practical significance for predicting the influence of the power frequency electric field around the high-voltage transmission line on the human and the environment by calculating the power frequency electric field. Therefore, the numerical calculation and analysis of the power frequency electric field around the power transmission line gradually becomes a research hot spot. In this respect, most of researches are performed by considering a complex environment around a power transmission line and then performing electric field calculation and analysis by using CSM, but few methods for improving electric field calculation accuracy are mentioned, and high calculation accuracy is a precondition for accurately predicting electric field influence.
The accuracy of the potential verification and electric field calculation depends on the setting of the analog charge in the wire. The setting method of analog charges in the lead is generally two, one is a central single method, and the method is to set 1 infinite length analog line charge in the center of the lead; the other is an empirical formula method, which is to set a plurality of analog line charges in a conducting wire and determine the positions of the analog line charges according to the empirical formula. However, the method for setting 1 line charge has limited calculation accuracy and cannot well meet the potential boundary condition of the lead; the empirical formula method is not suitable for the situation that the simulated line charges in the lead are 2-5, and if too many line charges are arranged, the calculated amount is obviously increased, meanwhile, the position is determined by human experience, and the CSM is not guaranteed to be optimal. It is contemplated that optimized CSM may be used to achieve optimal placement of analog charge in the conductive line.
Disclosure of Invention
Based on the analysis, a plurality of analog line charges are uniformly arranged on concentric circles of each equivalent conducting wire, only the distribution circle radius of the analog line charges is used as an optimization variable, and PSOCSM is utilized to search the distribution circle radius of the analog charges which minimizes potential verification errors, so that the optimal arrangement of the analog charges of the conducting wires is realized. And then calculating an electric field by the corresponding simulated charge group, and comparing the check point potential error of the distributed circle radius optimization method and the center single method with the electric field calculation precision so as to embody the high precision of electric field calculation after the distributed circle radius optimization method is adopted.
The invention provides a setting method of analog charges in a wire in electric field calculation of a power transmission line, which comprises the following steps:
a. uniformly arranging a plurality of analog line charges on a certain concentric circle of each equivalent conducting wire, and arranging corresponding matching points and check points on the surface of the conducting wire;
b. calculating an analog charge matrix;
c. calculating the potential error of the check point;
d. and obtaining the radius of the simulated charge distribution circle which minimizes the potential error of the check point by using a particle swarm optimization simulated charge method, calculating the coordinates and the charge quantity of the corresponding simulated charge group, and finally calculating the electric field around the transmission line by using the simulated charge group.
Further, in step b, the analog charge τ matrix is:
in the middle ofFor equivalent conductor surface potential, p i,k Potential coefficients at the i-th matching point for the kth analog charge:
ε 0 is the dielectric constant of air, x i ,y i X is the abscissa, ordinate, x of the ith matching point k ,y k Is the abscissa of the kth analog line charge. Equation (2) already contains the effect of the x-axis mirror charge and takes any point on the x-axis as a reference potential.
Further, in step d, the adaptive value function of PSOCSM is:
the optimization objective is minimum { fitness }. Wherein n is the check point number,modulo the n-th checkpoint potential calculation value,/->The standard value of the potential of the nth check point is modulo. The average relative error of the modulus of the calculated value of the potential of the i check points and the modulus of the standard value is calculated in the formula (3).
In step d, the radius of the simulated charge distribution circle minimizing the potential error of the check point is obtained by using PSOCSM, the coordinates of the corresponding simulated charge group are obtained by mathematical relationship, and the simulated charge amount τ matrix is obtained by equation (1).
In step d, after obtaining the coordinates of the analog charge group and the electric quantity τ matrix that minimize the fitness, the magnitudes of the x-direction component and the y-direction component of the electric field intensity of a certain measuring point are calculated according to the following formula:
τ in i 、x i 、y i The electric quantity and the abscissa of the electric charge of the ith analog line are respectively, and x and y are the abscissas and ordinates of a certain measuring point. Equations (4) and (5) already contain the effect on the x-axis mirror charge and take any point on the x-axis as the reference potential. It can be seen that E x 、E y Plural number, the total electric field intensity of a certain measuring point is
The invention has the beneficial effects that: according to the method for setting the simulated charges in the conducting wire in the electric field calculation of the power transmission line, the optimal setting of the simulated charges in the conducting wire can be realized, and the precision of the electric field calculation around the power transmission line is improved.
Drawings
The invention is further described below with reference to the accompanying drawings and examples:
FIG. 1 is a flow chart of the optimization algorithm particle swarm optimization simulated charge method used in the present invention;
FIG. 2 is a schematic diagram showing the arrangement of analog charge, matching points and check points in each equivalent conducting wire of the present invention;
FIG. 3 is a schematic diagram showing the comparison of the relative error of the potential of the calibration points on the surface of the equivalent conductive wire with the conventional analog charge setting method;
FIG. 4 is a graph of the relative error of the electric field value at 1.5m above ground versus the COMSOL simulation for the methods herein and the conventional method.
Detailed Description
As shown in fig. 1-4, the method for setting the analog charges in the conducting wire in the electric field calculation of the power transmission line provided by the invention comprises the following steps:
a. uniformly arranging a plurality of analog line charges on a certain concentric circle of each equivalent conducting wire, and arranging corresponding matching points and check points on the surface of the conducting wire;
b. calculating an analog charge matrix;
c. calculating the potential error of the check point;
d. and obtaining the radius of the simulated charge distribution circle which minimizes the potential error of the check point by using a particle swarm optimization simulated charge method, calculating the coordinates and the charge quantity of the corresponding simulated charge group, and finally calculating the electric field around the transmission line by using the simulated charge group. The invention can effectively realize the optimal setting of the analog charge in the lead, and compared with the conventional setting method (a central single method) of the analog charge in the lead, the electric potential error of the check point is greatly reduced, and the electric field calculation precision is improved.
In this embodiment, in step a, regarding the setting of analog charges in the equivalent wire, it is conventional practice to set 1 infinite analog line charge in the geometric center of the equivalent wire, and use a total of 3 infinite analog line charges equivalent to a three-phase power line, referred to herein as a center single-phase method. In order to better meet the surface potential boundary condition of the equivalent wire without greatly increasing the calculated amount, 4 infinite line charges are uniformly arranged on the concentric circle of each equivalent wire, as shown in fig. 2.
In this embodiment, in step b, the analog charge τ matrix is:
in the middle ofFor equivalent conductor surface potential, p i,k Potential coefficients at the i-th matching point for the kth analog charge:
ε 0 is the dielectric constant of air, x i ,y i X is the abscissa, ordinate, x of the ith matching point k ,y k Is the abscissa of the kth analog line charge. Equation (2) already contains the effect of the x-axis mirror charge and takes any point on the x-axis as a reference potential.
In this embodiment, in step c, after the coordinate of the analog line charge corresponding to the radius of the analog charge distribution circle and the electric quantity τ matrix are obtained, the matrix inverse operation of the formula (1) is usedAnd (3) calculating the potential of the check point, and performing potential check by calculating the elements of the P matrix as shown in the formula (2). The fitness function of PSOCSM is:
the optimization objective is minimum { fitness }. Wherein n is the check point number,modulo the n-th checkpoint potential calculation value,/->The standard value of the potential of the nth check point is modulo. The average phase of the modulus of the calculated value of the electric potential of the 12 check points and the modulus of the standard value is calculated in the formula (3)Error is corrected.
Taking the power frequency electric field of an actual 500kV three-phase transmission line as an example, the calibration point potential relative error curve of the conventional method and the conventional method is shown in fig. 3.
In this embodiment, in step d, the radius of the simulated charge distribution circle that minimizes fitness is obtained by using a particle swarm optimization simulated charge method, the coordinates of the corresponding simulated charge set are obtained by a mathematical relationship, and the simulated charge amount τ matrix is obtained by equation (1).
In this embodiment, in step d, after obtaining the coordinates of the analog charge group and the electric quantity τ matrix that minimize the fitness, the magnitudes of the x-direction component and the y-direction component of the electric field intensity of a certain measurement point are calculated according to the following formula:
τ in i 、x i 、y i The electric quantity and the abscissa of the electric charge of the ith analog line are respectively, and x and y are the abscissas and ordinates of a certain measuring point. Equations (4) and (5) already contain the effect on the x-axis mirror charge and take any point on the x-axis as the reference potential. It can be seen that E x 、E y Plural number, the total electric field intensity of a certain measuring point is
The relative error plot of the calculated electric field value at 1.5m above ground versus the COMSOL simulation value for the conventional and the methods herein is shown in fig. 4.
Finally, it is noted that the above embodiments are only for illustrating the technical solution of the present invention and not for limiting the same, and although the present invention has been described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications and equivalents may be made thereto without departing from the spirit and scope of the technical solution of the present invention, which is intended to be covered by the scope of the claims of the present invention.
Claims (5)
1. The setting method of the analog charge in the conducting wire in the electric field calculation of the power transmission line is characterized by comprising the following steps:
a. uniformly arranging a plurality of analog line charges on a certain concentric circle of each equivalent conducting wire, and arranging corresponding matching points and check points on the surface of the conducting wire;
b. calculating an analog charge matrix;
c. calculating the potential error of the check point;
d. and obtaining the radius of the simulated charge distribution circle which minimizes the potential error of the check point by using a particle swarm optimization simulated charge method, calculating the coordinates and the charge quantity of the corresponding simulated charge group, and finally calculating the electric field around the transmission line by using the simulated charge group.
2. The method for setting analog charges in a conductor in electric field calculation of a power transmission line according to claim 1, wherein in step b, the analog charge τ matrix is:
in the middle ofFor equivalent conductor surface potential, p i,k Potential coefficients at the i-th matching point for the kth analog charge:
ε 0 is the dielectric constant of air, x i ,y i X is the abscissa, ordinate, x of the ith matching point k ,y k The abscissa of the charge for the kth analog line; equation (2) already contains the effect of the x-axis mirror charge generation, andany point on the x-axis is taken as a reference potential.
3. The method for setting analog charges in a conductor in electric field calculation of a power transmission line according to claim 1, wherein in step d, an adaptive value function of PSOCSM is:
the optimization objective is minimum { fitness }, where n is the checkpoint sequence number,modulo the n-th checkpoint potential calculation value,/->The average relative error of the modulus of the calculated value of the potential of the i check points and the modulus of the standard value is calculated as the modulus of the n-th check point potential standard value in the formula (3).
4. The method for setting simulated charges in a conductor in electric field calculation of a power transmission line according to claim 3, wherein in step d, a radius of a simulated charge distribution circle minimizing a check point potential error is obtained by using PSOCSM, coordinates of a corresponding simulated charge group are obtained by a mathematical relationship, and a simulated charge amount τ matrix is obtained by formula (1).
5. The method for setting analog charges in a conducting wire in electric field calculation of a power transmission line according to claim 3 or 4, wherein in step d, after obtaining coordinates of an analog charge group minimizing a fitness and an electric quantity τ matrix, the magnitudes of x-direction components and y-direction components of the electric field intensity of a certain measuring point are calculated by the following formula:
τ in i 、x i 、y i The electric quantity and the abscissa of the electric charge of the ith analog line are respectively, and x and y are the abscissas and ordinates of a certain measuring point;
equations (4) and (5) already contain the effect on the x-axis mirror line charge, and take any point on the x-axis as a reference potential;
know E x 、E y Plural number, the total electric field intensity of a certain measuring point is
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910975391.1A CN110738009B (en) | 2019-10-14 | 2019-10-14 | Setting method of analog charges in conducting wire in electric field calculation of power transmission line |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910975391.1A CN110738009B (en) | 2019-10-14 | 2019-10-14 | Setting method of analog charges in conducting wire in electric field calculation of power transmission line |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110738009A CN110738009A (en) | 2020-01-31 |
CN110738009B true CN110738009B (en) | 2023-08-04 |
Family
ID=69268920
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910975391.1A Active CN110738009B (en) | 2019-10-14 | 2019-10-14 | Setting method of analog charges in conducting wire in electric field calculation of power transmission line |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110738009B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111859713B (en) * | 2020-08-03 | 2023-11-24 | 国网重庆市电力公司电力科学研究院 | Indoor substation power frequency electric field optimizing device |
CN112765897B (en) * | 2021-01-30 | 2022-07-05 | 上海工程技术大学 | Two-dimensional magnetic induction intensity control method based on overhead transmission line |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107315903A (en) * | 2017-05-23 | 2017-11-03 | 浙江大学 | A kind of intelligent analysis of electric field system |
CN109344431A (en) * | 2018-08-24 | 2019-02-15 | 国网安徽省电力有限公司建设分公司 | The method for accurately calculating conductive line surfaces electric field strength based on Analogue charge method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5568419B2 (en) * | 2010-09-06 | 2014-08-06 | 株式会社リコー | Surface charge distribution measuring method and surface charge distribution measuring apparatus |
-
2019
- 2019-10-14 CN CN201910975391.1A patent/CN110738009B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107315903A (en) * | 2017-05-23 | 2017-11-03 | 浙江大学 | A kind of intelligent analysis of electric field system |
CN109344431A (en) * | 2018-08-24 | 2019-02-15 | 国网安徽省电力有限公司建设分公司 | The method for accurately calculating conductive line surfaces electric field strength based on Analogue charge method |
Non-Patent Citations (3)
Title |
---|
张庆辉等."基于模拟电荷法的特高压直流输电线路表面场强计算方法".《清远职业技术学院学报》.2015,第第8卷卷(第第1期期),第56-59页. * |
彭孟杰."绝缘子工频电场逆向检测及优化方法研究".《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》.2012,(第第5期期),摘要、正文第2-4章. * |
谭传武 ; 陈卫兵 ; 邹豪杰 ; 李长云 ; .一种新颖自偏置带隙基准电压源的设计.湖南工业大学学报.2010,(第03期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN110738009A (en) | 2020-01-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110738009B (en) | Setting method of analog charges in conducting wire in electric field calculation of power transmission line | |
CN110135019A (en) | A kind of loss of power cable and core temperature prediction technique | |
CN108387824B (en) | Oiled paper insulation state evaluation method for improving return voltage measurement mode | |
CN112287297B (en) | Electric energy meter quality consistency evaluation method and medium based on random sampling | |
CN109344431B (en) | Method for accurately calculating electric field intensity on surface of wire based on analog charge method | |
Van Lanen et al. | Simulation of interstrand coupling loss in cable-in-conduit conductors with JackPot-AC | |
CN113609749B (en) | Current calculation method applicable to multiple scenes and based on magnetic field signals | |
CN111337797A (en) | XLPE power cable main insulation non-uniform thermal aging state evaluation method based on modified Debye model | |
CN112285435B (en) | Equivalent simulation method of high-power magnetic field radiation source | |
Su et al. | Parameter extraction of photovoltaic single-diode model using integrated current–voltage error criterion | |
CN108761279B (en) | Polar coordinate grid searching method and system for positioning partial discharge source | |
CN112685958B (en) | SiC MOSFET blocking voltage determination method based on neural network | |
CN112001059B (en) | Method and device for establishing broadband model of flexible direct current converter valve submodule | |
CN110646707B (en) | Grounding grid fault diagnosis method based on micro-processing method and particle swarm optimization | |
CN117056875A (en) | Cable transmission performance analysis method and device based on test data | |
CN109799391B (en) | Conductor alternating current resistance measuring method, system and computer storage medium | |
CN114818330A (en) | Fitting method of time-thickness relation curve and paint film thickness calculation method | |
CN106169409A (en) | Use multiple inductance and capacitance apparatus to use multiple condition of plasma to determine the system and method for matching network model | |
CN115483884A (en) | Power performance evaluation method for batch photovoltaic modules and related equipment | |
CN115034111A (en) | GIS electrostatic field calculation method based on U-net convolution neural network | |
CN112765842A (en) | Optimization design method for combined insulator voltage-sharing structure | |
CN116451544B (en) | Intelligent optimization design method for high-power high-frequency transformer | |
CN112765897B (en) | Two-dimensional magnetic induction intensity control method based on overhead transmission line | |
Talaat et al. | Multidimensional Optimization of CSM With Different Algorithms for Investigating the Electric Field in the Vicinity of a Point Electrode | |
CN110398635A (en) | The computation model of ground resistance |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |