CN106874565A - A kind of computational methods of rainy day transmission line of electricity lower section three-dimensional electric field - Google Patents
A kind of computational methods of rainy day transmission line of electricity lower section three-dimensional electric field Download PDFInfo
- Publication number
- CN106874565A CN106874565A CN201710034706.3A CN201710034706A CN106874565A CN 106874565 A CN106874565 A CN 106874565A CN 201710034706 A CN201710034706 A CN 201710034706A CN 106874565 A CN106874565 A CN 106874565A
- Authority
- CN
- China
- Prior art keywords
- dimensional
- charge
- electric
- points
- 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.)
- Granted
Links
- 230000005684 electric field Effects 0.000 title claims abstract description 87
- 230000005540 biological transmission Effects 0.000 title claims abstract description 39
- 230000005611 electricity Effects 0.000 title claims abstract description 38
- 238000000205 computational method Methods 0.000 title claims abstract description 20
- 238000009826 distribution Methods 0.000 claims abstract description 38
- 238000004088 simulation Methods 0.000 claims abstract description 26
- 238000000034 method Methods 0.000 claims abstract description 25
- 238000012795 verification Methods 0.000 claims description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 4
- 210000004027 cell Anatomy 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 7
- 229910000831 Steel Inorganic materials 0.000 description 5
- 239000004020 conductor Substances 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 5
- 239000010959 steel Substances 0.000 description 5
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 2
- 240000002853 Nelumbo nucifera Species 0.000 description 2
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/30—Polynomial surface description
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Geometry (AREA)
- Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Evolutionary Computation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Computer Hardware Design (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- General Engineering & Computer Science (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention relates to a kind of computational methods of rainy day transmission line of electricity lower section three-dimensional electric field, the method comprises the following steps：(1) three-dimensional system of coordinate is set up, the distribution situation according to transmission line of electricity and shaft tower sets up the threedimensional model of transmission line of electricity and shaft tower in three-dimensional system of coordinate using three-dimensional line unit simulation charge method；(2) raindrop model is set up using three-dimensional point charge simulation charge method in three-dimensional system of coordinate according to rain fall；(3) according to the coordinate of field strength region to be calculated, the electric-field intensity distribution of field strength region to be calculated is calculated according to field intensity computing formula and principle of stacking.Compared with prior art, result of calculation of the present invention is more true and reliable, the Electric Field Distribution during rainy day of all kinds of transmission line of electricity lower sections can be calculated, with applicability higher.
Description
Technical field
The present invention relates to a kind of computational methods of transmission line of electricity lower section three-dimensional electric field, more particularly, to a kind of rainy day power transmission line
The computational methods of road lower section three-dimensional electric field.
Background technology
With China's super-pressure and the high speed development of UHV transmission, the electromagnetic environment problem thus brought is more next
More by enjoying society to be paid close attention to its people.For actual transmission line of electricity, it is necessary to consider the influence of conducting wire sag, shaft tower, its
The calculating of electric field is complicated three dimensional field problem.Meanwhile, circuit actual motion under complicated weather environment, wherein with the rainy day and
Based on fine day, the electric field of rainy day is still not clear with the difference of the electric field of fine day, and at present to the rainy day when electric field calculating it is very few,
So being necessary for the three-dimensional electric field analysis under complicated weather condition.
The content of the invention
The purpose of the present invention is exactly to provide a kind of rainy day transmission line of electricity for the defect for overcoming above-mentioned prior art to exist
The computational methods of lower section three-dimensional electric field.
The purpose of the present invention can be achieved through the following technical solutions：
A kind of computational methods of rainy day transmission line of electricity lower section three-dimensional electric field, the method comprises the following steps：
(1) three-dimensional system of coordinate is set up, the distribution situation according to transmission line of electricity and shaft tower is in three-dimensional system of coordinate using three-dimensional
Line unit simulation charge method sets up the threedimensional model of transmission line of electricity and shaft tower；
(2) raindrop model is set up using three-dimensional point charge simulation charge method in three-dimensional system of coordinate according to rain fall；
(3) according to the coordinate of field strength region to be calculated, calculate to be calculated according to field intensity computing formula and principle of stacking
The electric-field intensity distribution of field strength region.
Three-dimensional line unit simulation charge method sets up transmission line of electricity and the threedimensional model of shaft tower is specially in step (1)：Will be defeated
The wire or bar of electric line are taken and be equivalent to respectively some sections of end to end compositions of three-dimensional line elementary charge, and electric charge is located at three-dimensional line list
First center, match point and checkpoint are arranged on three-dimensional line cell surface, according to match point and checkpoint carry out charge matching and
Verification determines the linear charge density of three-dimensional line unit.
Three-dimensional point charge simulation charge method is set up raindrop model and is specially in step (2)：Determine two according to rain fall
The distance between adjacent raindrop are d meters, determine the position of raindrop on zone boundary, x-axis, y-axis in three-dimensional system of coordinate, z-axis side
Upwards with d meters for the adjacent raindrop of interval setting, and rule completes the setting of raindrop position distribution in setting regions like this, in region
Each raindrop is correspondingly arranged a three-dimensional point electric charge, and electric charge is located at raindrop center, and match point and checkpoint are arranged on raindrop surface,
Charge matching is carried out according to match point and checkpoint and verification determines the quantity of electric charge of three-dimensional point electric charge.
The linear charge density of described three-dimensional line unit is tried to achieve especially by following manner：
(a1) quantity of electric charge for assuming 2 points of three-dimensional line elementary charge head and the tail is τ1And τ2, and obtain head and the tail 2 points corresponding
Current potential V with point1And V2And 2 points of current potential V of corresponding checkpoint of head and the tail3And V4；
(a2) it is P to calculate 2 points of coefficients of potential to corresponding match point of head and the tail1And P2, by τ1=V1/P1, τ2=V2/P2
Try to achieve the quantity of electric charge τ of 2 points of head and the tail1And τ2；
(a3) it is P to calculate 2 points of coefficients of potential to corresponding checkpoint of head and the tail respectively3And P4, it is calculated by following formula
The current potential V of checkpoint31And V41：V31=P3*τ1, V41=P4*τ2；
(a4) V is asked for respectively31With V3Difference and V41With V4Difference, when difference be less than setting value when perform step
(a5), otherwise return to step (a1)；
(a5) electric charge in three-dimensional line elementary charge is linear distribution, by 2 points of quantity of electric charge τ of head and the tail1And τ2Calculate
To the linear charge density of three-dimensional line unit.
The quantity of electric charge of three-dimensional point electric charge is tried to achieve by following manner：
(b1) quantity of electric charge for assuming three-dimensional point electric charge is q, and obtains the current potential V of the match point on three-dimensional point electric charge5；
(b2) coefficient of potential for calculating three-dimensional dot-charge pair match point is P5, by q=V/P5Try to achieve the electricity of three-dimensional point electric charge
Lotus amount；
(b3) coefficient of potential for calculating three-dimensional dot-charge pair checkpoint is P6, by V51=P5* q is calculated checkpoint
Current potential V51；
(b4) V is asked for51With V5Difference, the quantity of electric charge of three-dimensional point electric charge is determined when difference is less than setting value for q, otherwise
Return to step (b1).
Step (3) calculates the electric-field intensity of each point in field strength region to be calculated so as to obtain electric-field intensity respectively
Distribution, if the point for asking for electric-field intensity is P points, it is as follows that P points electric-field intensity calculates formula：
EPx=EPx lines+EPx points
EPy=EPy lines+EPy points
EPz=EPz lines+EPz points,
EPx、EPyAnd EPzCorrespond to component of the P point electric field strengths in x directions, y directions and z directions, EPx lines、EPy linesAnd EPz linesIt is right
Should be all three-dimensional line elementary charges in space P points electric-field intensity x directions, the component in y directions and z directions and, EPx points、
EPy pointsAnd EPz pointsCorrespond to all three-dimensional point electric charges in space P points electric-field intensity x directions, y directions and z directions component
With.
The distance between adjacent raindrop d is tried to achieve by equation below：
Wherein, D is the diameter of known raindrop, vbIt is terminal velocity, r is the quantum of rainfall in the t time periods, ρWaterIt is water
Density, ρIt is emptyIt is the density of air, g is acceleration of gravity, and T is a total time for rainfall, and R is the quantum of rainfall of T time section, t
It is the certain time period in T time section.
Certain three-dimensional line elementary charge P points electric-field intensity x directions, y directions and z directions component E 'Pxxian、
E′PyxianWith E 'PzxianIt is calculated by following manner：
If the line unit starting point of three-dimensional line elementary charge is P1(x1,y1,z1), terminal is P2(x2,y2,z2), if line unit is long
It is L to spend, and the linear charge density of three-dimensional line unit is τ (u)=au+b, and a, b are constant, and u=Lt (0≤t≤1), P point coordinates is P
(x, y, z), the then current potential that three-dimensional line elementary charge is produced in P pointsFor：
Wherein, ε0It is permittivity of vacuum,
A'=a'L=(- τ2+τ1)L
B'=b'=- τ1
τ (0)=τ1, τ (L)=τ2, b=τ1, a=(τ2-τ1)/L。
Certain three-dimensional point electric charge P points electric-field intensity x directions, y directions and z directions component E 'Pxdian、E′Pydian
With E 'PzdianIt is calculated by following manner：
If three-dimensional point electric charge coordinate is Q (x1,y1,z1), P point coordinates is P (x, y, z), then three-dimensional point electric charge is produced in P points
Current potentialFor：
Wherein, d is that three-dimensional point electric charge coordinate is the distance of Q to P points, and q is the quantity of electric charge of three-dimensional point electric charge, and ε is situated between for air
Electric constant；
And then,
Compared with prior art, the invention has the advantages that：
(1) present invention combines three-dimensional line unit simulation charge method and three-dimensional point charge simulation charge method, while using field intensity
Principle of stacking, it is proposed that a kind of computational methods for calculating transmission line of electricity lower section near-earth three-dimensional electric field during the rainy day, for calculating
Transmission line of electricity lower section near-earth electric field during the rainy day, so as to consider conducting wire sag, shaft tower, rainy day factor comprehensively so that result of calculation
It is more true and reliable, the Electric Field Distribution during rainy day of all kinds of transmission line of electricity lower sections can be calculated, with applicability higher；
(2) can to simulate heavy rain, moderate rain and light rain etc. by the change of the distance between adjacent raindrop d different for the present invention
Weather conditions so that result of calculation is more fitted actual condition, as a result relatively reliable.
Brief description of the drawings
Fig. 1 is the FB(flow block) of the computational methods of rainy day transmission line of electricity lower section of the invention three-dimensional electric field；
Fig. 2 is three-dimensional line unit simulation charge model；
Fig. 3 is three-dimensional point charge model；
Fig. 4 is the threedimensional model of shaft tower；
Fig. 5 is sag wire schematic diagram；
Fig. 6 be sunny weather under, shaft tower centre-to-centre spacing ground 1.5m at three-dimensional electric field distribution map；
Fig. 7 be light rain weather under, shaft tower centre-to-centre spacing ground 1.5m at three-dimensional electric field distribution map；
Fig. 8 is that fine day compares figure with light rain day Electric Field Distribution at shaft tower centre-to-centre spacing ground 1.5m；
Fig. 9 is that the rainy day Electric Field Distribution of different rainfall grades compares figure at shaft tower centre-to-centre spacing ground 1.5m；
Figure 10 is that below the phase of side at ground 1.5m, the Electric Field Distribution under different meteorology bars compares figure.
Specific embodiment
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.
Embodiment
As shown in figure 1, a kind of computational methods of rainy day transmission line of electricity lower section three-dimensional electric field, the method comprises the following steps：
(1) three-dimensional system of coordinate is set up, the distribution situation according to transmission line of electricity and shaft tower is in three-dimensional system of coordinate using three-dimensional
Line unit simulation charge method sets up the threedimensional model of transmission line of electricity and shaft tower；
(2) raindrop model is set up using three-dimensional point charge simulation charge method in three-dimensional system of coordinate according to rain fall；
(3) according to the coordinate of field strength region to be calculated, calculate to be calculated according to field intensity computing formula and principle of stacking
The electric-field intensity distribution of field strength region.
Three-dimensional line unit simulation charge method sets up transmission line of electricity and the threedimensional model of shaft tower is specially in step (1)：Will be defeated
The wire or bar of electric line are taken and be equivalent to respectively some sections of end to end compositions of three-dimensional line elementary charge, and electric charge is located at three-dimensional line list
First center, match point and checkpoint are arranged on three-dimensional line cell surface, according to match point and checkpoint carry out charge matching and
Verification determines the linear charge density of three-dimensional line unit.
Three-dimensional point charge simulation charge method is set up raindrop model and is specially in step (2)：Determine two according to rain fall
The distance between adjacent raindrop are d meters, determine the position of raindrop on zone boundary, x-axis, y-axis in three-dimensional system of coordinate, z-axis side
Upwards with d meters for the adjacent raindrop of interval setting, and rule completes the setting of raindrop position distribution in setting regions like this, in region
Each raindrop is correspondingly arranged a three-dimensional point electric charge, and electric charge is located at raindrop center, and match point and checkpoint are arranged on raindrop surface,
Charge matching is carried out according to match point and checkpoint and verification determines the quantity of electric charge of three-dimensional point electric charge.
If transmission line length is LlRice, numerical modeling setting is carried out using three-dimensional line unit simulation electric charge, and line unit is long
It is l meters to spend, and circuit is by Ll/ l sections of end to end the composition of line unit, electric charge is located at wire center, and match point and checkpoint are set
In conductive line surfaces.Shaft tower is made up of some steel columns, is taken wherein one section and is modeled explanation, due to the cylindrical use three-dimensional line of steel column
Unit simulation electric charge carries out numerical modeling setting, if this section of steel column length is LTM, line element length is lTM, it is by end to end
LT/lTSection line unit is constituted, and electric charge is located at steel column center, and match point and checkpoint are arranged on steel column surface.Raindrop spatial distribution
The distance between a length of x m in region, a width of y m, a height of z m, two adjacent raindrop are d m, first determine a rain on zone boundary
The position of drop, then adjacent raindrop are set with d m on x-axis, y-axis, z-axis direction, and rule completes raindrop position in region like this
Put distribution to set, each raindrop is correspondingly arranged a three-dimensional point electric charge in region, electric charge is located at raindrop center, match point and verification
Point is arranged on raindrop surface.By three-dimensional line unit simulation electric charge and the correlation computations formula of three-dimensional point electric charge, institute is obtained first
There is the coefficient of potential matrix [P] of match point, because the current potential matrix [V] of each match point is, it is known that can be by [V]=[P] [Q] premultiplication
[P]-1Try to achieve the electricity matrix [Q] of each electric charge.In order to determine amount of charge whether set rationally, quantities of charge calculate whether just
Really, it is necessary to be verified, the coefficient of potential matrix [P '] of checkpoint is tried to achieve, by [V ']=[P '] [Q], [V '] is checkpoint
Current potential matrix is calculated, each check point is calculated and is calculated current potential, each check point is calculated into the corresponding actual potential of current potential enters
Row error calculation, it is rational to illustrate that amount of charge is set if requirement is met, and it is correct that quantities of charge is calculated.Finally, may be used
Space any point P is tried to achieve in x-axis, y-axis, z-axis by the field intensity computing formula of three-dimensional line unit simulation electric charge and three-dimensional point electric charge
On field intensity component and close field intensity, and combine field intensity that field intensity superposition theorem forms all electric charges in space in P points according to arrow
Amount superposition obtains total intensity.
Specifically：The linear charge density of three-dimensional line unit is tried to achieve especially by following manner：
(a1) quantity of electric charge for assuming 2 points of three-dimensional line elementary charge head and the tail is τ1And τ2, and obtain head and the tail 2 points corresponding
Current potential V with point1And V2And 2 points of current potential V of corresponding checkpoint of head and the tail3And V4；
(a2) it is P to calculate 2 points of coefficients of potential to corresponding match point of head and the tail1And P2, by τ1=V1/P1, τ2=V2/P2
Try to achieve the quantity of electric charge τ of 2 points of head and the tail1And τ2；
(a3) it is P to calculate 2 points of coefficients of potential to corresponding checkpoint of head and the tail respectively3And P4, it is calculated by following formula
The current potential V of checkpoint31And V41：V31=P3*τ1, V41=P4*τ2；
(a4) V is asked for respectively31With V3Difference and V41With V4Difference, when difference be less than setting value when perform step
(a5), otherwise return to step (a1)；
(a5) electric charge in three-dimensional line elementary charge is linear distribution, by 2 points of quantity of electric charge τ of head and the tail1And τ2Calculate
To the linear charge density of three-dimensional line unit.
The quantity of electric charge of three-dimensional point electric charge is tried to achieve by following manner：
(b1) quantity of electric charge for assuming three-dimensional point electric charge is q, and obtains the current potential V of the match point on three-dimensional point electric charge5；
(b2) coefficient of potential for calculating three-dimensional dot-charge pair match point is P5, by q=V/P5Try to achieve the electricity of three-dimensional point electric charge
Lotus amount；
(b3) coefficient of potential for calculating three-dimensional dot-charge pair checkpoint is P6, by V51=P5* q is calculated checkpoint
Current potential V51；
(b4) V is asked for51With V5Difference, the quantity of electric charge of three-dimensional point electric charge is determined when difference is less than setting value for q, otherwise
Return to step (b1).
Step (3) calculates the electric-field intensity of each point in field strength region to be calculated so as to obtain electric-field intensity respectively
Distribution, if the point for asking for electric-field intensity is P points, it is as follows that P points electric-field intensity calculates formula：
EPx=EPx lines+EPx points
EPy=EPy lines+EPy points
EPz=EPz lines+EPz points,
EPx、EPyAnd EPzCorrespond to component of the P point electric field strengths in x directions, y directions and z directions, EPx lines、EPy linesAnd EPz linesIt is right
Should be all three-dimensional line elementary charges in space P points electric-field intensity x directions, the component in y directions and z directions and, EPx points、
EPy pointsAnd EPz pointsCorrespond to all three-dimensional point electric charges in space P points electric-field intensity x directions, y directions and z directions component
With.
The distance between adjacent raindrop d is tried to achieve by equation below：
Wherein, D is the diameter of known raindrop, vbIt is terminal velocity, r is the quantum of rainfall in the t time periods, ρWaterIt is water
Density, ρIt is emptyIt is the density of air, g is acceleration of gravity, and T is a total time for rainfall, and R is the quantum of rainfall of T time section, t
It is the certain time period in T time section.
Certain three-dimensional line elementary charge P points electric-field intensity x directions, y directions and z directions component E 'Pxxian、
E′PyxianWith E 'PzxianIt is calculated by following manner：
If the line unit starting point of three-dimensional line elementary charge is P1(x1,y1,z1), terminal is P2(x2,y2,z2), if line unit is long
It is L to spend, and the linear charge density of three-dimensional line unit is τ (u)=au+b, and a, b are constant, and u=Lt (0≤t≤1), P point coordinates is P
(x, y, z), the then current potential that three-dimensional line elementary charge is produced in P pointsFor：
Wherein, ε0It is permittivity of vacuum,
A '=a ' L=(- τ2+τ1)L
B '=b '-τ1
τ (0)=τ1, τ (L)=τ2, b=τ1, a=(τ2-τ1)/L。
Certain three-dimensional point electric charge P points electric-field intensity x directions, y directions and z directions component E 'Pxdian、E′Pydian
With E 'PzdianIt is calculated by following manner：
If three-dimensional point electric charge coordinate is Q (x1,y1,z1), P point coordinates is P (x, y, z), then three-dimensional point electric charge is produced in P points
Current potentialFor：
Wherein, d is that three-dimensional point electric charge coordinate is the distance of Q to P points, and q is the quantity of electric charge of three-dimensional point electric charge, and ε is situated between for air
Electric constant；
And then,
The three-dimensional electric field of computing electric power line lower section near-earth, traditional two-dimentional computational methods cannot meet requirement.Therefore, base
In the theoretical foundation of Analogue charge method, three-dimensional point charge simulation charge method is derived, and has combined three-dimensional line unit simulation charge method,
The threedimensional model of shaft tower, sag wire and raindrop is established, is distributed to simulate rainy day environment using raindrop are orderly in space, counted
The near-earth electric field of the rainy day let it pass under different rainfall intensities.Shaft tower and sag wire need to use three-dimensional line unit simulation when modeling
Electric charge, three-dimensional point electric charge is used for the foundation of raindrop model, with reference to the principle of stacking of field intensity, has obtained three-dimensional point electric charge and three-dimensional line
The field intensity formula of unit hybrid analog-digital simulation charge method.
Three-dimensional line unit simulation charge model as shown in Figure 2, it is assumed that P1It is line unit starting point, along P1P2Direction foundation office
Portion coordinate u, and it is L to set line element length, then unit inside any point Q (x3,y3,z3) coordinate can be obtained by following formula：
Wherein, L is line element length, u ∈ [0, L], linear charge density (i.e. quantities of charge), τ linear distributions in unit：
τ (u)=au+b (2)
Wherein, a, b are undetermined, then the current potential that any point P (x, y, z) is produced in fieldFor：
D is line unit internal field source point to P point distances, ε0It is permittivity of vacuum.Substitution of variable is made to the limit of integration, u=is made
Lt (0≤t≤1), then formula (3) can be expressed as：
If τ (0)=τ1, τ (L)=τ2, b=τ1, a=(τ2-τ1)/L, integration obtains the coefficient of potential and is：
P is the coefficient of potential.P1Electric charge and P2The coefficient of potential of electric charge is：
Due toGenerally known quantity, convolutionTherefore τ can be tried to achieve1And τ2Value, and then try to achieve a,
The value of b.
In order to meet boundary condition that ground potential is zero, it is necessary to add image charge, P1Corresponding image charge and P2It is right
The coefficient of potential of the image charge answered is：
L'=x2-x1
M'=y2-y1
N'=(- z2)-(-z1)=z1-z2
E'=l'2+m'2+n'2
F'=-2 [l (x-x1)+m(y-y1)+n(z-(-z1))]
G'=(x-x1)2+(y-y1)2+(z-(-z1))2 (9)
From electric field superposition principle, the actual electric-field intensity of space any point P points be non-mirror image and mirror image two parts it
With its expression formula is as follows：
In the present invention, needs are modeled according to the rainy day, derived three-dimensional point charge simulation charge method, it is assumed that Q is one in space
Point electric charge, P is any point in space, and two point coordinates are as shown in Figure 3.
Point charge Q can be calculated to the distance of spatial point P by following formula,
The current potential of P points is：
Wherein, q is quantities of charge.
Understand Q point electric charges the coefficient of potential be：
Due toGenerally known quantity, convolutionCan be in the hope of q.
The electric Field Calculation formula of P points is as follows：
In order that result of calculation is of universal significance, the present embodiment 500kV model of power transmission system common to China is carried out
Modeling, Tower Model uses wine glass-shaped shaft tower, wire type to use 4 × LGJ-400/35, and sub-conductor spacing is 0.45m, alternate
Away from being 12m, phase conductor suspension height is 31m.If circuit phase voltage presses the change of three-phase symmetrical sinusoidal rule, calculate using effective
Value, rated line voltage 500kV considers that calculating voltage takes 1.05 times of rated voltage, i.e. U with actual motionA=303.1kV, UB
=1.5 × 102-j2.6×102KV, UC=-1.5 × 102+j2.6×102kV.The threedimensional model of transmission tower is as shown in figure 4, arc
The wire schematic diagram that hangs down is as shown in Figure 5.
Meanwhile, in order to avoid the result of calculation that actual rainy day raindrop size, the isoparametric real-time change of raindrop spacing are brought
Particularity, is distributed to simulate the rain of different rainfall intensities in space in order using the three-dimensional point charge model of different radii herein
Its environment.The rainy day of different rainfall intensities, size, quantity, the spacing of raindrop all can produce different degrees of shadow to space electric field
Ring.
In formula, T is a total time for rainfall；R is the quantum of rainfall of T time section；T be T time section in sometime
Section；R is the quantum of rainfall in the t time periods；ρWaterIt is the density of water；ρIt is emptyIt is the density of air；G is acceleration of gravity；D is raindrop
Diameter；D is raindrop spacing distance；vbIt is terminal velocity.
Raindrop are arranged on using three-dimensional point charge model, charge concentration in the centre of sphere, current potential match point and current potential checkpoint
Spherical surface.With reference to data and formula (15)-(17) of table 1, the relevant parameter of raindrop under different rainfall intensity grades is calculated
As shown in table 2.
Raindrop average-size under the different rainfall grades of table 1
The relevant parameter of raindrop under the different rainfall grades of table 2
Fig. 6 be fair weather under, shaft tower center, [x (- 1.5,1.5), y (- 1.5, the 1.5)] square region at ground 1.5m
The Electric Field Distribution in domain.Fig. 7 be light rain weather under, shaft tower center, away from ground 1.5m at [x (- 1.5,1.5), y (- 1.5,1.5)] pros
The Electric Field Distribution in shape region.It can be found that there are many small spikes in electric field from Fig. 7, the presence of raindrop is illustrated so that electric field is sent out
Distortion is given birth to.
Fig. 8 be side phase lower section away from ground 1.5m at, the comparing of fine day and light rain weather Electric Field Distribution, from figure it can be found that
The electric field strength profile of light rain weather illustrates that the presence of light rain result in this phenomenon, the presence of light rain less than the curve of fine day
There is shielding action to the electric field of its inside.According to electrostatic screen principle, because the presence of light rain area of space, and rainwater are made
Be conductor, the inside for regarding conductor as that the inside of its area of space can be approximate, but because between raindrop gap presence, should
Interior volume is not completely enclosed conductor, therefore it has certain shielding action to electric field, but can not be complete
Full-shield electric field.
Fig. 9 is the comparing of the rainy day Electric Field Distribution of different rainfall grades at shaft tower centre-to-centre spacing ground 1.5m, and Figure 10 is under the phase of side
At ground 1.5m, the Electric Field Distribution under DIFFERENT METEOROLOGICAL CONDITIONS compares, therefrom it can be found that the corresponding Electric Field Distribution curve of light rain for side
Electric Field Distribution corresponding higher than moderate rain, Electric Field Distribution corresponding higher than heavy rain illustrates that shielding action of the heavy rain to electric field is most strong, its
Secondary is moderate rain, is then finally light rain.Because heavy rain, moderate rain, light rain are different in the density of spatial distribution, distribution of the heavy rain in space
Density more than moderate rain and light rain distribution density, the raindrop space that heavy rain is formed closer to an entirety, so to the screen of electric field
The effect of covering is more obvious.
Table 3 without ground wire, have shaft tower, have sag model field intensity comparative result under DIFFERENT METEOROLOGICAL CONDITIONS
By fine day and the numerical computations of rainy day, and protect what the committee (ICNIRP) formulated with reference to international Non-ionizing radiation
Directive/guide and standard on electromagnetic field exposure limits, corresponding data are listed in table 3.As shown in Table 3, the maximum field of fine day and rainy day
Respectively less than China's occupation exposes electric-field intensity to the open air and the public exposes electric-field intensity requirement to the open air by force.
Claims (9)
1. computational methods of a kind of rainy day transmission line of electricity lower section three-dimensional electric field, it is characterised in that the method comprises the following steps：
(1) three-dimensional system of coordinate is set up, the distribution situation according to transmission line of electricity and shaft tower uses three-dimensional line list in three-dimensional system of coordinate
First Analogue charge method sets up the threedimensional model of transmission line of electricity and shaft tower；
(2) raindrop model is set up using three-dimensional point charge simulation charge method in three-dimensional system of coordinate according to rain fall；
(3) according to the coordinate of field strength region to be calculated, electric field to be calculated is calculated according to field intensity computing formula and principle of stacking
The electric-field intensity distribution of intensity area.
2. a kind of computational methods of rainy day transmission line of electricity lower section according to claim 1 three-dimensional electric field, it is characterised in that step
Suddenly three-dimensional line unit simulation charge method sets up transmission line of electricity and the threedimensional model of shaft tower is specially in (1)：By leading for transmission line of electricity
Line or bar are taken and be equivalent to respectively some sections of end to end compositions of three-dimensional line elementary charge, and electric charge is located at three-dimensional line unit center,
Three-dimensional line cell surface is arranged on point and checkpoint, charge matching is carried out according to match point and checkpoint and verification is determined three
Tie up the linear charge density of line unit.
3. a kind of computational methods of rainy day transmission line of electricity lower section according to claim 2 three-dimensional electric field, it is characterised in that step
Suddenly three-dimensional point charge simulation charge method is set up raindrop model and is specially in (2)：According to rain fall determine two adjacent raindrop it
Between distance be d meters, determine the position of raindrop on zone boundary, be with d meters in the x-axis, y-axis, z-axis direction in three-dimensional system of coordinate
The adjacent raindrop of interval setting, and rule completes the setting of raindrop position distribution, each raindrop pair in region in setting regions like this
One three-dimensional point electric charge should be set, and electric charge is located at raindrop center, and match point and checkpoint are arranged on raindrop surface, according to match point
Charge matching is carried out with checkpoint and verification determines the quantity of electric charge of three-dimensional point electric charge.
4. a kind of computational methods of rainy day transmission line of electricity lower section according to claim 2 three-dimensional electric field, it is characterised in that institute
The linear charge density of the three-dimensional line unit stated is tried to achieve especially by following manner：
(a1) quantity of electric charge for assuming 2 points of three-dimensional line elementary charge head and the tail is τ1And τ2, and obtain 2 points of corresponding match points of head and the tail
Current potential V1And V2And 2 points of current potential V of corresponding checkpoint of head and the tail3And V4；
(a2) it is P to calculate 2 points of coefficients of potential to corresponding match point of head and the tail1And P2, by τ1=V1/P1, τ2=V2/P2Try to achieve
2 points of quantity of electric charge τ from beginning to end1And τ2；
(a3) it is P to calculate 2 points of coefficients of potential to corresponding checkpoint of head and the tail respectively3And P4, verification is calculated by following formula
The current potential V of point31And V41：V31=P3*τ1, V41=P4*τ2；
(a4) V is asked for respectively31With V3Difference and V41With V4Difference, when difference be less than setting value when perform step (a5),
Otherwise return to step (a1)；
(a5) electric charge in three-dimensional line elementary charge is linear distribution, by 2 points of quantity of electric charge τ of head and the tail1And τ2It is calculated three
Tie up the linear charge density of line unit.
5. computational methods of a kind of rainy day transmission line of electricity lower section according to claim 3 three-dimensional electric field, it is characterised in that three
The quantity of electric charge of dimension point electric charge is tried to achieve by following manner：
(b1) quantity of electric charge for assuming three-dimensional point electric charge is q, and obtains the current potential V of the match point on three-dimensional point electric charge5；
(b2) coefficient of potential for calculating three-dimensional dot-charge pair match point is P5, by q=V/P5Try to achieve the electric charge of three-dimensional point electric charge
Amount；
(b3) coefficient of potential for calculating three-dimensional dot-charge pair checkpoint is P6, by V51=P5* q is calculated the current potential of checkpoint
V51；
(b4) V is asked for51With V5Difference, determine that the quantity of electric charge of three-dimensional point electric charge, for q, is otherwise returned when difference is less than setting value
Step (b1).
6. a kind of computational methods of rainy day transmission line of electricity lower section according to claim 3 three-dimensional electric field, it is characterised in that step
Suddenly (3) calculate the electric-field intensity of each point in field strength region to be calculated so as to obtain electric-field intensity distribution respectively, if asking
The point of power taking field intensity is P points, and it is as follows that P points electric-field intensity calculates formula：
EPx=EPx lines+EPx points
EPy=EPy lines+EPy points
EPz=EPz lines+EPz points,
EPx、EPyAnd EPzCorrespond to component of the P point electric field strengths in x directions, y directions and z directions, EPx lines、EPy linesAnd EPz linesCorrespond to
In space all three-dimensional line elementary charges P points electric-field intensity x directions, the component in y directions and z directions and, EPx points、EPy points
And EPz pointsCorrespond to all three-dimensional point electric charges in space P points electric-field intensity x directions, the component in y directions and z directions and.
7. a kind of computational methods of rainy day transmission line of electricity lower section according to claim 3 three-dimensional electric field, it is characterised in that phase
The distance between adjacent raindrop d is tried to achieve by equation below：
Wherein, D is the diameter of known raindrop, vbIt is terminal velocity, r is the quantum of rainfall in the t time periods, ρWaterIt is the density of water,
ρIt is emptyIt is the density of air, g is acceleration of gravity, and T is a total time for rainfall, and R is the quantum of rainfall of T time section, when t is T
Between certain time period in section.
8. computational methods of a kind of rainy day transmission line of electricity lower section according to claim 6 three-dimensional electric field, it is characterised in that certain
Individual three-dimensional line elementary charge P points electric-field intensity x directions, y directions and z directions component E 'Pxxian、E′PyxianWith
E′PzxianIt is calculated by following manner：
If the line unit starting point of three-dimensional line elementary charge is P1(x1,y1,z1), terminal is P2(x2,y2,z2), if line element length is
L, the linear charge density of three-dimensional line unit is τ (u)=au+b, and a, b are constant, u=Lt (0≤t≤1), P point coordinates be P (x, y,
Z), then the current potential that three-dimensional line elementary charge is produced in P pointsFor：
Wherein, ε0It is permittivity of vacuum,
A'=a'L=(- τ2+τ1)L
B'=b'=- τ1
τ (0)=τ1, τ (L)=τ2, b=τ1, a=(τ2-τ1)/L。
9. computational methods of a kind of rainy day transmission line of electricity lower section according to claim 6 three-dimensional electric field, it is characterised in that certain
Individual three-dimensional point electric charge P points electric-field intensity x directions, y directions and z directions component E 'Pxdian、E′PydianWith E 'PzdianIt is logical
Following manner is crossed to be calculated：
If three-dimensional point electric charge coordinate is Q (x1,y1,z1), P point coordinates is P (x, y, z), the then electricity that three-dimensional point electric charge is produced in P points
PositionFor：
Wherein, d is that three-dimensional point electric charge coordinate is the distance of Q to P points, and q is the quantity of electric charge of three-dimensional point electric charge, and ε is that air dielectric is normal
Number；
And then,
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710034706.3A CN106874565B (en) | 2017-01-17 | 2017-01-17 | Method for calculating three-dimensional electric field below power transmission line in rainy days |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710034706.3A CN106874565B (en) | 2017-01-17 | 2017-01-17 | Method for calculating three-dimensional electric field below power transmission line in rainy days |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106874565A true CN106874565A (en) | 2017-06-20 |
| CN106874565B CN106874565B (en) | 2020-05-26 |
Family
ID=59158393
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710034706.3A Active CN106874565B (en) | 2017-01-17 | 2017-01-17 | Method for calculating three-dimensional electric field below power transmission line in rainy days |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106874565B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110488103A (en) * | 2019-08-28 | 2019-11-22 | 清华大学 | Spatial electronic field measurement system |
| CN110738013A (en) * | 2019-09-26 | 2020-01-31 | 国网江西省电力有限公司电力科学研究院 | Method and system for calculating surface field intensity distribution of split conductor |
| CN111667393A (en) * | 2019-08-26 | 2020-09-15 | 福建数***信息科技有限公司 | Method and terminal for simulating raining in virtual scene |
| CN113239557A (en) * | 2021-05-17 | 2021-08-10 | 江南大学 | ADSS optical cable hanging point position determination method based on analog charge method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101349720A (en) * | 2008-08-28 | 2009-01-21 | 重庆电力科学试验研究院 | High pressure overhead power line power frequency electric field shielding analysis method |
| CN104036058A (en) * | 2013-03-05 | 2014-09-10 | 辽宁省电力有限公司营口供电公司 | Electric transmission line galloping distribution diagram drawing method |
| CN105427190A (en) * | 2015-12-09 | 2016-03-23 | 国网重庆市电力公司电力科学研究院 | Three-dimensional power-frequency electric-field calculating method of ground below UHVAC power transmission line in complex terrain |
| CN106096110A (en) * | 2016-06-03 | 2016-11-09 | 国网电力科学研究院武汉南瑞有限责任公司 | Flashover of power transmission circuit caused by windage yaw method for early warning based on numerical weather forecast |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106168987A (en) * | 2016-04-18 | 2016-11-30 | 国网重庆市电力公司电力科学研究院 | Complex scene power frequency electric field three dimensional analysis method and system under a kind of power transmission line |
-
2017
- 2017-01-17 CN CN201710034706.3A patent/CN106874565B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101349720A (en) * | 2008-08-28 | 2009-01-21 | 重庆电力科学试验研究院 | High pressure overhead power line power frequency electric field shielding analysis method |
| CN104036058A (en) * | 2013-03-05 | 2014-09-10 | 辽宁省电力有限公司营口供电公司 | Electric transmission line galloping distribution diagram drawing method |
| CN105427190A (en) * | 2015-12-09 | 2016-03-23 | 国网重庆市电力公司电力科学研究院 | Three-dimensional power-frequency electric-field calculating method of ground below UHVAC power transmission line in complex terrain |
| CN106096110A (en) * | 2016-06-03 | 2016-11-09 | 国网电力科学研究院武汉南瑞有限责任公司 | Flashover of power transmission circuit caused by windage yaw method for early warning based on numerical weather forecast |
Non-Patent Citations (5)
| Title |
|---|
| HUICHUN XIE 等: "Statistical Analysis of Radio Interference of 1000 kV UHV AC Double-circuit Transmission Lines in Foul Weather", 《CSEE JOURNAL OF POWER AND ENERGY SYSTEMS》 * |
| 刘伟: "降雨对棒―板短空气间隙直流放电影响的仿真及试验研究", 《中国优秀硕士学位论文全文数据库》 * |
| 刘大为 等: "雨滴运动分析", 《科学技术与工程》 * |
| 彭迎 等: "模拟电荷法计算特高压架空线路3维工频电场", 《高电压技术》 * |
| 马爱清 等: "考虑环境因素时雨滴对直流输电导线电晕振动的影响", 《中国电机工程学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111667393A (en) * | 2019-08-26 | 2020-09-15 | 福建数***信息科技有限公司 | Method and terminal for simulating raining in virtual scene |
| CN111667393B (en) * | 2019-08-26 | 2023-07-07 | 福建数***信息科技有限公司 | Method and terminal for simulating raining in virtual scene |
| CN110488103A (en) * | 2019-08-28 | 2019-11-22 | 清华大学 | Spatial electronic field measurement system |
| CN110738013A (en) * | 2019-09-26 | 2020-01-31 | 国网江西省电力有限公司电力科学研究院 | Method and system for calculating surface field intensity distribution of split conductor |
| CN110738013B (en) * | 2019-09-26 | 2023-04-07 | 国网江西省电力有限公司电力科学研究院 | Method and system for calculating surface field intensity distribution of split conductor |
| CN113239557A (en) * | 2021-05-17 | 2021-08-10 | 江南大学 | ADSS optical cable hanging point position determination method based on analog charge method |
| CN113239557B (en) * | 2021-05-17 | 2024-03-22 | 江南大学 | ADSS optical cable hanging point position determining method based on analog charge method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106874565B (en) | 2020-05-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106874565A (en) | A kind of computational methods of rainy day transmission line of electricity lower section three-dimensional electric field | |
| AD | Ionospheric electrodynamics using magnetic apex coordinates | |
| CN102129484B (en) | Method and device for generating digitalized flat cross-section diagram of transmission line | |
| CN105404938B (en) | A kind of line walking method for optimizing route shortening the patrolling transmission line time | |
| CN102368286B (en) | Computation method of surface synthesis electric field of small grounding object under DC (direct current) power transmission line | |
| CN103995986B (en) | The method that the power frequency electric field of rolling ground overhead transmission line is determined based on Analogue charge method | |
| CN108205889A (en) | Freeway traffic flow Forecasting Methodology based on convolutional neural networks | |
| Liu et al. | Experiments on natural snow distribution around simplified building models based on open air snow-wind combined experimental facility | |
| CN103683086B (en) | Power transmission circuit caused by windage transient state analyzing method under the effect of motion thunderstorm wind | |
| CN102628913A (en) | Electric transmission line three-dimensional lightning calculating method based on electric transmission line structure and terrain | |
| CN106168987A (en) | Complex scene power frequency electric field three dimensional analysis method and system under a kind of power transmission line | |
| CN108205088A (en) | A kind of parallel radio interference of high voltage ac/dc circuit calculates and optimization method | |
| CN108961094A (en) | Wind leaning fault method for early warning based on transmission line of electricity minimum air void online measuring | |
| Zhang et al. | Power line simulation for safety distance detection using point clouds | |
| CN107045731A (en) | The method being modeled based on cloud data to substation bus bar | |
| Struzewska et al. | Impact of urban parameterization on high resolution air quality forecast with the GEM–AQ model | |
| Liu et al. | Road PV production estimation at city scale: A predictive model towards feasible assessing regional energy generation from solar roads | |
| CN103914737A (en) | Method for calculating whole circuit in-place weather information of electric transmission and transformation circuit | |
| CN103177300A (en) | Tower positioning method for transmission line | |
| CN111061897A (en) | Multi-scale automatic intercepting and classifying method suitable for high-resolution remote sensing image of power grid | |
| CN107391794A (en) | A kind of typhoon continuous stereo Wind-field Retrieval method | |
| CN105021903A (en) | A method for measuring and calculating electric field intensity nearby high voltage transmission lines on a hillside | |
| CN107704655A (en) | A kind of direct current transportation wire filth appraisal procedure and system | |
| CN106652031A (en) | 4D real dynamic display method for electric engineering design | |
| CN110866693B (en) | Overhead transmission line icing risk assessment method based on GIS model |
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 |