CN108846193A - A kind of calculation method for the HVDC transmission line formate field intensity and ion current density considering real-time weather conditions - Google Patents
A kind of calculation method for the HVDC transmission line formate field intensity and ion current density considering real-time weather conditions Download PDFInfo
- Publication number
- CN108846193A CN108846193A CN201810587694.1A CN201810587694A CN108846193A CN 108846193 A CN108846193 A CN 108846193A CN 201810587694 A CN201810587694 A CN 201810587694A CN 108846193 A CN108846193 A CN 108846193A
- Authority
- CN
- China
- Prior art keywords
- formula
- haze
- calculated
- density
- field intensity
- 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
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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
The present invention proposes that nominal electric field strength E is calculated in a kind of calculation method of HVDC transmission line formate field intensity and ion current density for considering real-time weather conditions, including optimizing application Analogue charge method；The formate field intensity and charge density of HVDC transmission line when not considering haze are calculated using flux collimation method, bloom field strength is calculated using the correction formula for considering climatic effect；The carrying capacity for calculating suspension droplet and haze particle, obtains total charge density；Calculate the conductive line surfaces charge density and electric field strength amplitude A value considered under the conditions of haze；Total formate field intensity and ion current density are calculated according to flux collimation method, ionic mobility is calculated using the correction formula for considering climatic effect.Algorithm synthesis of the invention considers influence of the weather conditions such as temperature, humidity, haze index to HVDC transmission line corona ion flow field, formate field intensity and ion current density of the more realistic calculating of energy under different weather conditions.
Description
Technical field
The invention belongs to electrical engineering transmission line of electricity electromagnetic environment calculating fields, particularly a kind of to consider real-time weather
The HVDC transmission line formate field intensity of condition and the calculation method of ion current density.
Background technique
D.C. high voltage transmission is widely applied in China's power transmission engineering, corona loss, radio interference, ground
The synthesis factors such as field and ion current density needs take in and control in construction, and this aspect, which has also been made, in the country largely grinds
Study carefully, and achieves great successes, since HVDC transmission line traverse distance is remote, the weather conditions for various regions of passing through are different,
Temperature, humidity differ larger in different geographical, and haze weather is multiple in recent years in addition, and new complicated weather conditions are to direct current transportation
Engineering brings new challenge, and under haze weather, there are a large amount of suspension droplet and suspended particulate substances in air, due to power transmission line
Charged ion is filled in space around road, the suspended particulate in haze can adsorb charged ion, new charged particle is formed,
It cannot ignore because quantity is big, this keeps the calculating of Ion Flow Field more complicated.Comprehensively consider the gas such as temperature, humidity, haze index
The calculation method in the corona ion flow field of the HVDC transmission line under the conditions of time has in actual project planning an design
Important actual reference.
In the calculating of Ion Flow Field, FInite Element or flux collimation method are generallyd use, FInite Element can preferably realize list
The calculating of pole and bipolar DC Ion Flow Field, but computational efficiency is not high, especially in the case where model meshes complexity, calculating speed
It is relatively slow；Flux collimation method premised on Deutsch assumes, solve equation it is simple, speed is fast, be not necessarily to mesh generation, calculated result with
Actual measured results are closer to, and meet certain accuracy requirement, the present invention calculated using flux collimation method formate field intensity and from
Subflow is able to satisfy engineering demand.
Summary of the invention
A kind of calculation method for the HVDC transmission line formate field intensity and ion current density considering real-time weather conditions,
Specific step is as follows：
Step 1：Nominal electric field strength E is calculated in optimizing application Analogue charge method；
Step 2：The formate field intensity and electricity of HVDC transmission line when not considering haze are calculated using flux collimation method
Lotus density, bloom field strength are calculated using the correction formula for considering climatic effect；
Step 3：The carrying capacity for calculating suspension droplet and haze particle, obtains total charge density；
Step 4：Calculate the conductive line surfaces charge density and electric field strength amplitude A value considered under the conditions of haze；
Step 5：Total formate field intensity and ion current density are calculated according to flux collimation method, ionic mobility is using consideration
The correction formula of climatic effect calculates.
Further, in step 2, (a) is for space any point P (x, y), first optimizing application Analogue charge method meter
The electric field line for passing through the point is calculated, to determine the relative position of P (x, y), if taking anode between positive wire and the earth
Each parameter of conducting wire is calculated, conversely, taking each parameter of cathode conductor；
(b) conductive line surfaces field strength is maintained at bloom field intensity value after corona, when conducting wire voltage-to-ground is U, conducting wire discharge inception voltage
For U0, then the A value A of conductive line surfacesiFor：
Ai=U0/U (5)
Since field strength is directly proportional with voltage, AiIt can also be denoted as：
Consider that the bloom field strength of temperature and humidity and haze index can be modified to：
E in formulaonFor the bloom field strength of conductive line surfaces；M is the coefficient of roughness for reacting conductive line surfaces situation, according to power industry
Standard DL/T《High voltage direct current aerial power transmission line fire protection technology》When calculating corona loss, it is leading under haze that m, which takes 0.47, m',
The line coefficient of roughness, is generally taken as 0.42；R is sub-conductor radius；δ is the relative density of air, takes 1 under the conditions of normal weather,
Consider to calculate when temperature influences by formula (7)：
Wherein t is actual temperature when conducting wire is run, t020 DEG C are taken as reference temperature, p is gas pressure intensity, usual situation
Under, p=p0For a standard atmospheric pressure,
The bloom field strength E of equivalent conducting wire can be calculated using above formulaonr,For the flat of each sub-conductor maximum surface field strength
Mean value；A is calculated separately out using formula (14) and (5)iAnd U0；
Substitution formula (26) calculates ρm；
(c) two ρ are selected in conductive line surfacesiInitial value,
ρi1=f1ρm (27)
ρi2=f2ρm (28)
If it is positive wire, f1Take 2, f2Take 3；If it is cathode conductor, f1Take 1.5, f2Take 3；
(d) it is calculated according to formula (24) and corresponds respectively to ρi1And ρi2Charge density value ρ1And ρ2, formula (25) then are substituted into, are asked
Obtain corresponding ρm1And ρm2,
Then further ρiInitial value be calculated as follows,
(e) by ρi3Substitution formula (24), acquires charge density ρ3, by ρ3Substitution formula (25) obtains corresponding to ρi3Average electricity
Lotus density pm3, ρi3It whether is ρiTrue value can be judged with following formula:
If above formula is set up, then it is assumed that ρi3It is ρiTrue value, otherwise enable ρi1=ρi2, ρi2=ρi3, and repeat step (c)
(d), until formula (30) are set up；
(f) ρ is acquirediTrue value after, calculate A according to formula (23), conjunction when A substitution formula (4) is not considered haze
At field strength Es,
Es=AE (4)
A in formulaiIt is the A value on polar conductor surface, ρiFor the charge density on polar conductor surface, ρmTo pass through space arbitrary point
Average charge density on power line, E are nominal field strength, EsFor the total electric field strength for considering haze.
Further, the formate field intensity E in step 3, when will not consider hazesIt substitutes into following formula (21), can be obtained
The saturation carrying capacity q of suspension droplet and haze particles,
Consider haze total electrical charge volume density be：
ρ=ρe+ρf+ρp=ρe+(Nf+Np)·qs (22)
ρfFor the charge density of suspension droplet, ρpFor the space charge density of suspension haze particle, ρeIt is close for ionic space-charge
Degree.
Further, in step 4, in the situation known to the total electrical charge volume density ρ for considering haze, conductive line surfaces can be obtained
Charge density ρiCalculation formula：
If required point is located at positive polarity conductor area, ρiIt is calculated by formula (31), conversely, ρiIt is calculated by formula (32)；Obtain conducting wire
Surface charge density ρiAfterwards, formula (23) are substituted into and obtains the A value under the conditions of haze.
Further, in step 5：The A value that step 4 is acquired substitutes into formula (4) and (3), can acquire formate field intensity with
Ion current density,
J=KqρeEs+ρfvf+ρpvp (3)
Consider that the ionic mobility formula under the weather conditions such as temperature, humidity, haze index can be modified to：
In formula, KqFor ionic mobility, J is ion current density.
Detailed description of the invention
Fig. 1 is that the present invention considers the HVDC transmission line formate field intensity of real-time weather conditions and the meter of ion current density
Calculate method flow diagram；
Fig. 2 is the ground total electric field distribution under the different haze situations of the embodiment of the present invention；
Fig. 3 is the ground ion flow Density Distribution under the different haze situations of the embodiment of the present invention；
Fig. 4 is that the measured value of route ground formate field intensity is somebody's turn to do with what is calculated according to the method for the present invention in existing technical literature
Route ground formate field intensity comparison diagram；
Fig. 5 be in existing technical literature the measured value of route ground ion flow density with calculate according to the method for the present invention
The comparison diagram of the route ground ion flow density；
Fig. 6 is the HVDC transmission line formate field intensity and ion current density of the real-time weather conditions of consideration of the invention
The calculating surface chart of computing system；
Fig. 7 is the HVDC transmission line formate field intensity and ion current density of the real-time weather conditions of consideration of the invention
The calculation flow chart of computing system.
Specific embodiment
With reference to the accompanying drawings and examples, specific embodiments of the present invention will be described in further detail.
A kind of calculation method for the HVDC transmission line formate field intensity and ion current density considering real-time weather conditions,
Real-time weather conditions herein include temperature, humidity and haze situation.Theoretical basis of the invention is as follows：
The governing equation of Ion Flow Field
When haze weather, haze particle can be considered to be made of suspension droplet and subtle suspended solid particles object, haze particle
In corona bloom, Ion Flow Field formation and stabilization process, gas discharge can be had an impact, be mainly manifested in droplet and haze
The charged effect of particle, charged haze particle increases space charge density, accordingly, it is considered to the high-voltage dc transmission of haze particle
The governing equation of electric line Ion Flow Field is similar under normal weather, can be described with following equation:
J=KqρeEs+ρfvf+ρpvp (3)
In formula, EsFor total electric field strength, V/m；ρeFor ionic space-charge density, C/m3；ρfFor the charge of suspension droplet
Density, C/m3；ρpFor the space charge density of suspension haze particle, C/m3；KqFor ionic mobility, m2/Vs；vfFor suspension droplet
Movement speed, m/s；vpFor the movement speed of suspension haze particle, m/s；J is ion current density, A/m2。
Describe have space charge synthesis field equation be it is nonlinear, for convenience of solve, often make the following assumptions：
Space charge only influences amplitude of field strength, assumes without influencing its direction and Deutsch：
Es=AE (4)
A is scalar function in formula；
Conductive line surfaces field strength is maintained at bloom field intensity value after corona, and when conducting wire voltage-to-ground is U, conducting wire discharge inception voltage is
U0, then the A value A of conductive line surfacesiFor：
Ai=U0/U (5)
Ionic mobility is unrelated with field strength, is constant.The thickness of corona layer is ignored.Do not consider that the diffusion of ion is made
With.
Based on assumed above, the calculating of formate field intensity and ion current density can be calculated with flux collimation method.
Analysis of Influential Factors of the weather conditions to corona ion flow field
1. impact analysis of the weather conditions to bloom field strength
Influence of the temperature to corona discharge intensity
Document (Peek F.W.Dielectric phenomena in high-voltage engineering [M]:
Journal of the Franklin Institute.McGraw-Hill Book Company, Incorporated, 1930)
The Peek formula that data obtain through a large number of experiments is widely used in the calculating of direct current and transmission line of alternation current bloom field strength, needle
To DC power transmission line, the form of Peek formula is：
Wherein, EcFor the bloom field strength of conducting wire, E0For with reference to field intensity value, K is empirical, 0.301, r is usually takencTo lead
Line radius, m is the coefficient of roughness of conductive line surfaces, and for DC power transmission line, for value between 0.3~0.8, δ is week
The relative density for enclosing air, its calculation formula is：
Wherein t is actual temperature when conducting wire is run, t020 DEG C are taken as reference temperature, p is gas pressure intensity, usual situation
Under, p=p0For a standard atmospheric pressure.By above formula it is found that as the temperature rises, atmospheric density gradually decreases, lead
Causing the bloom field strength on transmission line of electricity surface reduces, i.e., transmission line of electricity is easier bloom, and the strength of discharge of corona can be with temperature
It increases and enhances.
Influence of the humidity to corona discharge intensity
A large amount of experiment and research the result shows that, there is the positive-negative polarity corona characteristic of direct current and influence in the humidity of air,
Wherein, Davies on the basis of Peek formula, by air humidity, take into account by this factor, proposes one to Peek formula
Correction formula:
Wherein, H is the absolute humidity of air, g/m3, δ is air relative density.Under normal circumstances, meteorological condition is shown
For relative humidity, the transformational relation of relative humidity and absolute humidity is：
Wherein h is relative humidity, and H is hydrometeor density, i.e. absolute humidity contained by reality, H in unit volumes of air0For equality of temperature
Degree is lower to be saturated hydrometeor density, and rule of thumb formula is expressed as saturation hydrometeor density：
Wherein H0For saturated steam density, e is natural logrithm, and T is absolute temperature.The temperature applicable range of the formula is
0~240 degree Celsius, with actual measured value positive and negative deviation less than 5%.Therefore, under the conditions of known to the relative humidity and temperature, absolutely
Humidity H can be calculated.
Influence of the haze to conductive surface bloom field strength
Long-term haze can be such that the pollution level of conducting wire increases, and influence the roughness of conducting wire, and then influence the bloom electricity of conducting wire
Pressure.Air humidity under the conditions of haze is an important factor for influencing conducting wire bloom field strength, and existing document studies this, and
Propose calculation formula.Bloom field strength under haze weather can be expressed as：
E'onFor the bloom field strength under haze, EonFor the bloom field strength under normal weather, m' is that the conducting wire under haze is coarse
Coefficient, being generally taken as 0.42, m is the coefficient of roughness under normal weather, is generally taken as 0.47.
To sum up, by the analysis of the influence factor to temperature, humidity and haze exponent pair bloom field strength, consider temperature and
Humidity and the bloom field strength of haze index can be modified to：
E in formulaonFor the bloom field strength (kV/cm) of conductive line surfaces；M is the coefficient of roughness for reacting conductive line surfaces situation, according to
Power industry standard DL/T《High voltage direct current aerial power transmission line fire protection technology》When calculating corona loss, it is mist that m, which takes 0.47, m',
The conducting wire coefficient of roughness under haze, is generally taken as 0.42；R is sub-conductor radius (cm)；δ is the relative density of air, in normal day
1 is taken under the conditions of gas, considers to calculate when temperature influences by formula (7).
To simplify the calculation, DC power transmission line can calculate equivalent conducting wire using above formula mostly with the replacement of equivalent conducting wire
Bloom field strength Eonr, comparison is calculated and is surveyed according to previous, takes the average value of each sub-conductor maximum surface field strength hereinRoot
According to formula (5) it is found that the value A on polar conductor surfaceiFor：
2. influence of the weather conditions to ionic mobility
Influence of the temperature to ionic mobility
Temperature will affect the average heat movement speed of air and the free travel of ion, and then influence the mobility of ion,
Ionic mobility formula in same gas is：
Wherein α is constant, and between 0.6~1.0, e is basic amount of electrons, and m is the quality of ion, and r is gas molecule
Radius, k are Boltzmann constant, and p is air pressure, and T is Kai Erwen temperature.
Under a certain atmospheric pressure, it is assumed that in the case where temperature is t, ionic mobility Kt, then relative to standard day
Under the conditions of ionic mobility K have：
Wherein t0It is 20 degrees Celsius, it can be seen that, as the temperature rises, ionic mobility is in rising trend.
Influence of the humidity to ionic mobility
Professor Yang Jinji when the behavior in the electric field of theoretical and test result analysis charged particle by pointing out, it is believed that
When the mass conservation of charged particle, the mobility of ion according to gas molecule mass M-0.5Variation, utilizes the conclusion, it is assumed that no
It is K with the ionic mobility under relative humidityHr, the ionic mobility under the conditions of normal weather is K, then has：
Wherein, M is the relative mass of air molecule under normal circumstances, is taken as 29, MHrFor the sky under different relative humidity
Gas molecule relative mass, is calculated as follows：
maFor the atmospheric density under normal temperature and pressure, ma=1.2kg/m3, H is the vapor content under different humidity, by formula above
(9) it calculates.
Influence of the haze to ionic mobility
A large amount of haze particle and suspension droplet will affect the mobility of ion in air, according to the calculating of ionic mobility
Formula can obtain：
Wherein KhFor the ionic mobility under the conditions of haze, K be it is normal under the conditions of ionic mobility, MhFor haze condition
Air molecule relative mass, M be it is normal under the conditions of air molecule relative mass, maFor the atmospheric density under normal temperature and pressure, H
For the density of droplet content in air, air vapor content when relative humidity is 100% is H=23g/m3, mpFor haze in air
The density of particle, as haze index select the concentration of PM2.5 for 310ug/m when serious pollution3, by each parameter values generation
Ionic mobility when room temperature seriously pollutes can be obtained by entering above formula.
To sum up, consider that the ionic mobility formula under the weather conditions such as temperature, humidity, haze index can be modified to：
3. influence of the weather conditions to space charge density
Suspension droplet when the influence of space charge density is mainly manifested in haze weather after transmission line of electricity bloom, in air
Can occur with haze particle it is charged, in general, ion is attached to droplet and haze particle by diffusional charging and electric-field charge two ways
On, it was found that electric-field charge plays a major role, due to the partial size of droplet and haze particle when particle size is greater than 0.5um
It is all larger than this value, diffusional charging can be ignored, and only consider electric-field charge herein.It is generated after particle adion occurs charged
When electric field and external electric field balance, the saturation carrying capacity of charged saturation, particle is：
R in formula0For particle radius, εrFor the relative dielectric constant of particle.
Haze is the droplet to suspend and the general designation of haze particle, and the constituent of haze is complicated, only studies PM2.5 herein to ion
The influence of stream.The partial size of suspension droplet is generally greater than 1um, takes 3.2um herein, and the average diameter of suspended particulate substance is
1.3um is taken herein, and after not considering that formate field intensity when haze calculates, substituting into formula (21) can be obtained suspension droplet and haze
The saturation carrying capacity q of particles, it is assumed that the density of suspended particulate substance is 1g/cm3, in haze exponent mpUnder conditions of given, unit
Haze number of particles N in volumepIt can calculate, the droplet number N of unit volumefTake 50 drops/cm respectively according to the class of pollution3To 300
Drop/cm3, specifically see the table below.To obtain the total electrical charge volume density for considering haze.
ρ=ρe+ρf+ρp=ρe+(Nf+Np)·qs (22)
Division according to China to haze grade herein carries out corresponding droplet and haze particle unit volume quantity
It calculates, as shown in the table.The computation model of this paper assumes the haze of each grade all with certain mist, i.e., in following table accordingly
Droplet and the sum of haze particle object quantity be total Charge particle quantity.
Droplet and haze number of particles under the different haze classes of pollution of table 1
Based on above theoretical and it is assumed that it is found that different weather conditions, such as temperature, humidity, haze index, to power transmission line
Bloom field strength, the ionic mobility on road have an impact, and haze has an impact to space charge density, these influence factors directly affect
The formate field intensity and ion current density in corona ion flow field will consider these influence factors, specifically in calculating process
Steps are as follows for calculating：
(1) nominal electric field strength is calculated in optimizing application Analogue charge method.
(2) formate field intensity and charge density for not considering transmission line of electricity when haze is charged is calculated using flux collimation method,
In calculating process, bloom field strength is calculated using the correction formula for considering climatic effect；
(3) carrying capacity for calculating suspension droplet and haze particle, obtains total charge density；
(4) the conductive line surfaces charge density and A value considered under the conditions of haze is calculated；
(5) total formate field intensity and ion current density are calculated according to flux collimation method, ionic mobility is using consideration weather
The correction formula of influence calculates.
Fig. 1 show the real-time weather conditions of consideration of the invention HVDC transmission line formate field intensity and ion stream it is close
The calculation method flow chart of degree, specific step is as follows：
(1) nominal electric field strength is calculated in optimizing application Analogue charge method.
(2) formate field intensity and charge density of transmission line of electricity when not considering haze, bloom are calculated using flux collimation method
Field strength is calculated using the correction formula for considering climatic effect；
Flux collimation method asks the calculation formula of synthesis field to have：
A in formulaiIt is the A value on polar conductor surface, ρiFor the charge density on polar conductor surface.ρmTo pass through space arbitrary point
Average charge density on power line, E are nominal field strength, and the thinking of formate field intensity is calculated according to flux collimation method, is needed following several
Step：
(a) for space any point P (x, y), optimizing application Analogue charge method first calculates the electric field line for passing through the point,
To determine the relative position of P (x, y), if taking each parameter of positive wire to be counted between positive wire and the earth
It calculates, conversely, taking each parameter of cathode conductor.
(b) A is calculated separately out using formula (14) and (5)iAnd U0, substitute into formula (26) and calculate ρm。
(c) two ρ are selected in conductive line surfacesiInitial value, taken with following formula, so that the energy in iterative process below
Reach 10 as early as possible-6The convergence precision of grade,
ρi1=f1ρm (27)
ρi2=f2ρm (28)
10 can be reached in the iterative process below as early as possible to make to calculate-6The convergence precision of grade, rule of thumb, if it is just
Polar conductor, f1Take 2, f2Take 3；If it is cathode conductor, f1Take 1.5, f2Take 3.
(d) it is calculated according to formula (24) and corresponds respectively to ρi1And ρi2Charge density value ρ1And ρ2, formula (25) then are substituted into, are asked
Obtain corresponding ρm1And ρm2, then further ρiInitial value be calculated as follows.
(e) ibid, by ρi3Substitution formula (24), acquires charge density ρ3, by ρ3Substitution formula (25) obtains corresponding to ρi3It is flat
Equal charge density ρm3, ρi3It whether is ρiTrue value can be judged with following formula:
If above formula is set up, then it is assumed that ρi3It is ρiTrue value, otherwise enable ρi1=ρi2, ρi2=ρi3, and (c) (d) is repeated,
Until above formula is set up.
(f) ρ is acquirediTrue value after, calculate A according to formula (23), conjunction when A substitution formula (4) is not considered haze
At field strength Es。
(3) carrying capacity for calculating suspension droplet and haze particle, obtains total charge density；
Formate field intensity E when haze is not consideredsAfter calculating, substitutes into formula (21), suspension droplet and haze can be obtained
The saturation carrying capacity q of particles.By qsSubstitute into the total electrical charge volume density that formula (22) obtain considering haze.
(4) the conductive line surfaces charge density and A value considered under the conditions of haze is calculated；
Total electric field at this time consists of three parts, the nominal electric field that conducting wire charge generates, charged droplet and haze particle
The electric field of generation, the electric field that space corona ion generates.After finding out the total charge density ρ for considering haze, because of the electricity of conductive line surfaces
Lotus density piWith A and AiAlso changed because increasing haze carrying capacity, need to recalculate.According to formula (24), in ρ
In known situation, conductive line surfaces charge density ρ can be obtainediCalculation formula：
If required point is located at positive polarity conductor area, ρiIt is calculated by formula (31), conversely, ρiIt is calculated by formula (32).Obtain conducting wire
Surface charge density ρiAfterwards, formula (23) are substituted into and obtains the A value under the conditions of haze.
(5) total formate field intensity and ion current density are calculated, ionic mobility is public using the amendment for considering climatic effect
Formula calculates.
After A value acquires, the formate field intensity for considering each weather conditions can be acquired by substituting into formula (4), which is substituted into
Formula (3) can acquire the ion current density for considering each weather conditions.
Calculation method of the invention is illustrated with specific embodiment below.
Using certain ± 800kV HVDC transmission line as model, calculating in temperature is 20 degrees Celsius, humidity 80%, no
With the formate field intensity and ion current density in the case of haze, line parameter circuit value is as follows, and wire type is 6 × LGJ-630/45, and son is led
Line radius is 1.68cm, and interpolar distance is 22m, and conducting wire distance away the ground is 21m, and division spacing is 0.45m, positive sub-conductor bloom
Voltage is taken as 604kV, and cathode sub-conductor discharge inception voltage is taken as -582kV, the conjunction under the conditions of several different hazes being calculated
It is as shown in Figures 2 and 3 at field strength and ion current density.
Find out from the calculated result of attached drawing 2 and attached drawing 3, ground formate field intensity and ion current density all phases under haze weather
Normal weather increased, and rise with the aggravation of the pollution level of haze weather.In the 6th grade of heavily contaminated and companion
With under the weather of heavy fog, when haze concentration is 310ug/m3, droplet number is 300 drops/cm3, formate field intensity maximum value is 28kV/m,
Than improving 19% under normal weather, have been approached China ground formate field intensity control standard 30kV/m has with the exacerbation of haze
It can exceed that control standard, ion stream reach 37.5nA/m under the conditions of haze2, than being higher by 44% under normal weather, so,
Formate field intensity characteristic distributions under the conditions of being considered as in the design of transmission line of electricity and construction in haze.
In order to confirm the validity and accuracy of calculation method of the invention, below with measured result in the prior art into
Row compares.
Existing technical literature：Johnson G.Electrical fields and ion current s of a 400kV
HVDC test line [J] .IEEE Transactions on Power Apparatus and Systems, 1983,102
(8)：2559-2568, under the bipolar direct current transmission line an of ± 400k formate field intensity and ion current density incited somebody to action
Nearly measurement in 7 months, obtains the route in normal weather, the greasy weather, the rainy day, the formate field intensity and ion current density under snowy day
A large amount of measured datas.In order to which the validity to this paper algorithm is verified, now using method of the invention to the prior art text
Earth's surface total electric field and ion current density of this route in fair weather and greasy weather in offering carries out simulation calculation, fair weather
When, temperature is set as 20 degrees Celsius, relative humidity 50%, and haze index takes 50ug/m3, and when the greasy weather, temperature is 20 degrees Celsius, phase
It is 90% to humidity, haze index takes 100ug/m3, and the actual measured results of calculated result and existing technical literature are carried out pair
Than.Line related parameter is：Double oidiospore wire radius are 1.91cm, and division spacing is 45.7cm, and polar conductor distance away the ground is
15.2m, interpolar is away from for 12.2m.
Such as Fig. 4, Fig. 5, the calculated value of method of the invention is shown compared between the measured value of existing technical literature.
It can be seen from the figure that algorithm of the invention and measurement result have preferable consistency, calculation method of the present invention is being demonstrated just
True property.The calculating peak value of formate field intensity and ion current density is less than normal than measured value near negative polarity, it may be possible to due to the present invention
Calculation method consider a certain proportion of haze particle, it is higher when bloom field strength is compared to the greasy weather to cause.
Based on the above calculation method, consideration temperature is developed, the HVDC transmission line of humidity, haze index synthesizes field
Strong and ion current density computing system, which can obtain the temperature, humidity and haze index in somewhere in real time, in given voltage
Grade, division number divide spacing, and interpolar after the line parameter circuit values such as sub-conductor radius, can calculate certain away from, conducting wire distance away the ground
The formate field intensity and ion current density of point.It is as shown in Figure 6 to calculate interface.
Fig. 7 shows the calculation flow chart of the software.After starting software, inputted in the input frame of city calculative
The title in city, such as " Beijing " click " updating weather " button, and software judges whether to have input city name, if not defeated
Enter, is then set to default city.If having input city name, the city name of software inquiry input whether there is in Urban Data library,
If it does not exist, prompting frame information is popped up, correct city name is re-entered, otherwise the automatic temperature for obtaining setting city from network
Degree, humidity and AQI value, and be arranged into following parameter box.Next remaining each parameter is inputted：Voltage class, is divided division number
Spacing, interpolar are split away from, the coordinate (X, Y, Z) of conducting wire distance away the ground, sub-conductor radius and target point." calculating " button is clicked, it is soft
Part first judges whether parameter inputs completely, if not inputted, pops up prompt information, supplements complete parameter, then judge each ginseng
Whether several formats is correct, such as erroneous input letter or other characters, pops up prompt information, again corrected parameter.All nothings
After accidentally, software transfer algorithm calculates the formate field intensity and ion current density of target point, and shows.
Algorithm synthesis of the invention considers the weather conditions such as temperature, humidity, haze index to HVDC transmission line
The influence in corona ion flow field, formate field intensity and ion stream of the more realistic calculating of energy under different weather conditions are close
Degree, the algorithm is simple and easy to do, while precision is also able to satisfy engine request.The computing system of algorithm establishment can be in real time according to the present invention
Ground obtains the weather conditions in somewhere, can calculate the formate field intensity of the point quickly according to transmission line parameter and the setting for calculating point
And ion current density, it can consider to mention when the influence factors such as each weather conditions in planning and design for HVDC transmission line
For effective reference frame.
Claims (5)
1. a kind of calculation method for the HVDC transmission line formate field intensity and ion current density for considering real-time weather conditions, tool
Steps are as follows for body：
Step 1：Nominal electric field strength E is calculated in optimizing application Analogue charge method；
Step 2：It is close that the formate field intensity of HVDC transmission line and charge when not considering haze are calculated using flux collimation method
Degree, bloom field strength are calculated using the correction formula for considering climatic effect；
Step 3：The carrying capacity for calculating suspension droplet and haze particle, obtains total charge density；
Step 4：Calculate the conductive line surfaces charge density and electric field strength amplitude A value considered under the conditions of haze；
Step 5：Total formate field intensity and ion current density are calculated according to flux collimation method, ionic mobility is using consideration weather
The correction formula of influence calculates.
2. the HVDC transmission line formate field intensity according to claim 1 for considering real-time weather conditions and ion stream are close
The calculation method of degree, the step 2 include：
(a) for space any point P (x, y), optimizing application Analogue charge method first calculates the electric field line for passing through the point, with true
The relative position of P (x, y) is determined, if taking each parameter of positive wire to be calculated, instead between positive wire and the earth
It, takes each parameter of cathode conductor；
(b) conductive line surfaces field strength is maintained at bloom field intensity value after corona, and when conducting wire voltage-to-ground is U, conducting wire discharge inception voltage is U0,
The then A value A of conductive line surfacesiFor：
Ai=U0/U (5)
Since field strength is directly proportional with voltage, AiIt can also be denoted as：
Consider that the bloom field strength of temperature and humidity and haze index can be modified to：
E in formulaonFor the bloom field strength of conductive line surfaces；M is the coefficient of roughness for reacting conductive line surfaces situation, according to power industry standard
DL/T《High voltage direct current aerial power transmission line fire protection technology》When calculating corona loss, m takes 0.47, m' thick for the conducting wire under haze
Rough coefficient, is generally taken as 0.42；R is sub-conductor radius；δ is the relative density of air, and 1 is taken under the conditions of normal weather, is considered
It is calculated when temperature influences by formula (7)：
Wherein t is actual temperature when conducting wire is run, t020 DEG C are taken as reference temperature, p is gas pressure intensity, it is generally the case that p
=p0For a standard atmospheric pressure,
The bloom field strength E of equivalent conducting wire can be calculated using above formulaonr,For being averaged for each sub-conductor maximum surface field strength
Value；A is calculated separately out using formula (14) and (5)iAnd U0；
Substitution formula (26) calculates ρm；
(c) two ρ are selected in conductive line surfacesiInitial value,
ρi1=f1ρm (27)
ρi2=f2ρm (28)
If it is positive wire, f1Take 2, f2Take 3；If it is cathode conductor, f1Take 1.5, f2Take 3；
(d) it is calculated according to formula (24) and corresponds respectively to ρi1And ρi2Charge density value ρ1And ρ2, formula (25) then are substituted into, are acquired pair
The ρ answeredm1And ρm2,
Then further ρiInitial value be calculated as follows,
(e) by ρi3Substitution formula (24), acquires charge density ρ3, by ρ3Substitution formula (25) obtains corresponding to ρi3Mean charge it is close
Spend ρm3, ρi3It whether is ρiTrue value can be judged with following formula:
If above formula is set up, then it is assumed that ρi3It is ρiTrue value, otherwise enable ρi1=ρi2, ρi2=ρi3, and repeat step (c) and
(d), until formula (30) are set up；
(f) ρ is acquirediTrue value after, calculate A according to formula (23), formate field intensity when A substitution formula (4) is not considered haze
Es,
Es=AE (4)
A in formulaiIt is the A value on polar conductor surface, ρiFor the charge density on polar conductor surface, ρmTo pass through the electric power of space arbitrary point
Average charge density on line, E are nominal field strength, EsFor the total electric field strength for considering haze.
3. the HVDC transmission line formate field intensity according to claim 2 for considering real-time weather conditions and ion stream are close
The calculation method of degree, the step 3 include：Formate field intensity E when will not consider hazesIt substitutes into following formula (21), can obtain
To the saturation carrying capacity q of suspension droplet and haze particles,
qs=12 π ε0Esr0 2εr/(εr+2) (21)
Consider haze total electrical charge volume density be：
ρ=ρe+ρf+ρp=ρe+(Nf+Np)·qs (22)
ρfFor the charge density of suspension droplet, ρpFor the space charge density of suspension haze particle, ρeFor ionic space-charge density.
4. the HVDC transmission line formate field intensity according to claim 3 for considering real-time weather conditions and ion stream are close
The calculation method of degree, the step 4 include：In the situation known to the total electrical charge volume density ρ for considering haze, wire table can be obtained
Density of surface charge ρiCalculation formula：
If required point is located at positive polarity conductor area, ρiIt is calculated by formula (31), conversely, ρiIt is calculated by formula (32)；Obtain conductive line surfaces
Charge density ρiAfterwards, formula (23) are substituted into and obtains the A value under the conditions of haze.
5. the HVDC transmission line formate field intensity according to claim 4 for considering real-time weather conditions and ion stream are close
The calculation method of degree, the step 5 include：The A value that the step 4 is acquired substitutes into formula (4) and (3), can acquire synthesis
Field strength and ion current density,
J=KqρeEs+ρfvf+ρpvp (3)
Consider that the ionic mobility formula under the weather conditions such as temperature, humidity, haze index can be modified to：
In formula, KqFor ionic mobility, J is ion current density.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810587694.1A CN108846193B (en) | 2018-06-08 | 2018-06-08 | Method for calculating composite field intensity and ion current density of high-voltage direct-current transmission line considering real-time climate conditions |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810587694.1A CN108846193B (en) | 2018-06-08 | 2018-06-08 | Method for calculating composite field intensity and ion current density of high-voltage direct-current transmission line considering real-time climate conditions |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108846193A true CN108846193A (en) | 2018-11-20 |
| CN108846193B CN108846193B (en) | 2022-08-30 |
Family
ID=64210723
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810587694.1A Active CN108846193B (en) | 2018-06-08 | 2018-06-08 | Method for calculating composite field intensity and ion current density of high-voltage direct-current transmission line considering real-time climate conditions |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108846193B (en) |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109387755A (en) * | 2018-12-12 | 2019-02-26 | 国网山东省电力公司电力科学研究院 | A kind of high voltage direct current corona discharge characteristic experimental apparatus system |
| CN109948252A (en) * | 2019-03-20 | 2019-06-28 | 河海大学 | A kind of ultra-high/extra-high voltage AC corona loss numerical computation method based on stream theory |
| CN112016808A (en) * | 2020-07-31 | 2020-12-01 | 山西晋浙环保科技有限公司 | High-efficiency electrostatic dust collector |
| CN112182920A (en) * | 2020-08-24 | 2021-01-05 | 中国电力科学研究院有限公司 | Iterative method for determining field intensity value of direct current transmission line synthetic electric field |
| CN112649676A (en) * | 2020-12-18 | 2021-04-13 | 国网西藏电力有限公司电力科学研究院 | Calculation method of high-altitude area mixed field intensity considering wind speed influence |
| CN113325239A (en) * | 2021-05-07 | 2021-08-31 | 中国电力科学研究院有限公司 | Method and system for calculating corona starting field intensity of wire based on synthetic electric field |
| CN113376444A (en) * | 2021-05-07 | 2021-09-10 | 中国电力科学研究院有限公司 | Method and system for calculating corona onset field intensity of wire based on ion current density |
| CN113625067A (en) * | 2021-08-12 | 2021-11-09 | 华北电力大学 | Device and method for measuring charge characteristics of suspended liquid drops in ion flow field |
| CN113721084A (en) * | 2021-07-22 | 2021-11-30 | 中国电力科学研究院有限公司 | Method and system for determining direct-current corona-starting field intensity of lead under different temperature and humidity conditions |
| CN114113815A (en) * | 2021-11-19 | 2022-03-01 | 华北电力大学 | Device and method for measuring charge characteristics of particulate matters in ion flow field |
| CN114137468A (en) * | 2021-10-21 | 2022-03-04 | 国网浙江省电力有限公司电力科学研究院 | Method and system for correcting field intensity measured by direct current electric field intensity detector |
| CN112182920B (en) * | 2020-08-24 | 2024-05-10 | 中国电力科学研究院有限公司 | Iterative method for determining field intensity value of DC transmission line composite electric field |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105203021A (en) * | 2015-10-21 | 2015-12-30 | 华北电力大学 | Measuring device and method of roughness coefficient of high-voltage conducting wire |
| CN107704655A (en) * | 2017-09-05 | 2018-02-16 | 国网浙江省电力公司电力科学研究院 | A kind of direct current transportation wire filth appraisal procedure and system |
-
2018
- 2018-06-08 CN CN201810587694.1A patent/CN108846193B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105203021A (en) * | 2015-10-21 | 2015-12-30 | 华北电力大学 | Measuring device and method of roughness coefficient of high-voltage conducting wire |
| CN107704655A (en) * | 2017-09-05 | 2018-02-16 | 国网浙江省电力公司电力科学研究院 | A kind of direct current transportation wire filth appraisal procedure and system |
Non-Patent Citations (5)
| Title |
|---|
| FU JIN ET AL.: "The Research on Influence Law of Different Environment on the Ion Flow Field of HVDC Transmission Line and Insulator Contamination", 《2016 IEEE 11TH CONFERENCE ON INDUSTRIAL ELECTRONICS AND APPLICATIONS (ICIEA)》 * |
| 崔翔 等: "高压直流输电线路离子流场计算方法研究进展", 《中国电机工程学报》 * |
| 杨帆 等: "雾霾天气下的直流输电线路离子流场分布特性及其影响因素", 《电工技术学报》 * |
| 谭震宇 等: "雾霾对高压直流输电线路电晕离子流场的影响", 《高电压技术》 * |
| 赵永生 等: "雾对高压直流输电线路离子流场的影响", 《中国电机工程学报》 * |
Cited By (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109387755A (en) * | 2018-12-12 | 2019-02-26 | 国网山东省电力公司电力科学研究院 | A kind of high voltage direct current corona discharge characteristic experimental apparatus system |
| CN109948252B (en) * | 2019-03-20 | 2022-08-26 | 河海大学 | Ultra/extra-high voltage alternating current corona loss numerical value calculation method based on streamer theory |
| CN109948252A (en) * | 2019-03-20 | 2019-06-28 | 河海大学 | A kind of ultra-high/extra-high voltage AC corona loss numerical computation method based on stream theory |
| CN112016808A (en) * | 2020-07-31 | 2020-12-01 | 山西晋浙环保科技有限公司 | High-efficiency electrostatic dust collector |
| CN112182920A (en) * | 2020-08-24 | 2021-01-05 | 中国电力科学研究院有限公司 | Iterative method for determining field intensity value of direct current transmission line synthetic electric field |
| CN112182920B (en) * | 2020-08-24 | 2024-05-10 | 中国电力科学研究院有限公司 | Iterative method for determining field intensity value of DC transmission line composite electric field |
| CN112649676A (en) * | 2020-12-18 | 2021-04-13 | 国网西藏电力有限公司电力科学研究院 | Calculation method of high-altitude area mixed field intensity considering wind speed influence |
| CN112649676B (en) * | 2020-12-18 | 2024-04-09 | 国网西藏电力有限公司电力科学研究院 | Calculation method of mixed field intensity of high-altitude area considering wind speed influence |
| CN113376444A (en) * | 2021-05-07 | 2021-09-10 | 中国电力科学研究院有限公司 | Method and system for calculating corona onset field intensity of wire based on ion current density |
| CN113376444B (en) * | 2021-05-07 | 2023-12-05 | 中国电力科学研究院有限公司 | Method and system for calculating corona onset field intensity of wire based on ion current density |
| CN113325239B (en) * | 2021-05-07 | 2023-12-05 | 中国电力科学研究院有限公司 | Method and system for calculating conductor corona onset field intensity based on composite electric field |
| CN113325239A (en) * | 2021-05-07 | 2021-08-31 | 中国电力科学研究院有限公司 | Method and system for calculating corona starting field intensity of wire based on synthetic electric field |
| CN113721084A (en) * | 2021-07-22 | 2021-11-30 | 中国电力科学研究院有限公司 | Method and system for determining direct-current corona-starting field intensity of lead under different temperature and humidity conditions |
| CN113721084B (en) * | 2021-07-22 | 2024-02-13 | 中国电力科学研究院有限公司 | Method and system for determining direct current corona onset field intensity of lead under different temperature and humidity conditions |
| CN113625067A (en) * | 2021-08-12 | 2021-11-09 | 华北电力大学 | Device and method for measuring charge characteristics of suspended liquid drops in ion flow field |
| CN113625067B (en) * | 2021-08-12 | 2022-06-10 | 华北电力大学 | Device and method for measuring charge characteristics of suspended liquid drops in ion flow field |
| CN114137468A (en) * | 2021-10-21 | 2022-03-04 | 国网浙江省电力有限公司电力科学研究院 | Method and system for correcting field intensity measured by direct current electric field intensity detector |
| CN114113815A (en) * | 2021-11-19 | 2022-03-01 | 华北电力大学 | Device and method for measuring charge characteristics of particulate matters in ion flow field |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108846193B (en) | 2022-08-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108846193A (en) | A kind of calculation method for the HVDC transmission line formate field intensity and ion current density considering real-time weather conditions | |
| CN106680603B (en) | A kind of characteristic line method calculating DC power transmission line three-dimensional Ion Flow Field | |
| CN108984819A (en) | A kind of method in HVDC transmission line corona ion flow field under the conditions of calculating haze | |
| Liu et al. | An efficient model to simulate stable glow corona discharges and their transition into streamers | |
| Xiao et al. | Calculation of 3-D ion-flow field at the crossing of HVdc transmission lines by method of characteristics | |
| Sheridan | The plasma sheath around large discs and ion collection by planar Langmuir probes | |
| Zhen et al. | 3-D finite-element method for calculating the ionized electric field and the ion current of the human body model under the UHVDC lines | |
| Hedtke et al. | Audible noise of hybrid AC/DC overhead lines: Comparison of different prediction methods and conductor arrangements | |
| CN102175930B (en) | Wind speed influence treatment method in measurement of total electric field of extraordinarily high-voltage DC (direct current) transmission line | |
| Mokrov et al. | Development of a positive corona from a long grounded wire in a growing thunderstorm field | |
| Chuan et al. | A hybrid approach for corona discharge in needle electrode configuration: in a large-scale space | |
| He et al. | Experiments and analysis of corona inception voltage under combined AC-DC voltages at various air pressure and humidity in rod to plane electrodes | |
| CN112649676A (en) | Calculation method of high-altitude area mixed field intensity considering wind speed influence | |
| Yi et al. | Conductor surface conditions effects on the ion-flow field of long-term operating conductors of the HVDC transmission line | |
| CN115825584A (en) | Lightning positioning optimization method and system for correcting influence of terrain and geological parameters on lightning electromagnetic wave transmission and medium | |
| Jaiswal et al. | Finite element modelling of ionized field quantities around a monopolar HVDC transmission line | |
| Liu et al. | Corona onset gradient of the bundle conductor on AC power lines under sand and dust weather condition at 2,200 m altitude | |
| Latham | Space charge generated by wind tunnel fires | |
| Donglai et al. | Measurement of surface charges on the dielectric film based on field mills under the HVDC corona wire | |
| Wang et al. | Ion flow field distribution near the crossing of two circuit UHVDC transmission lines | |
| Si et al. | Ion flow field of AC and DC parallel lines under different weather conditions | |
| Ma et al. | The effect of humidity and temperature on the total electric field of unipolar positive and negative HVDC lines in the natural environment | |
| Michelarakis et al. | An electrostatic approach for the calculation of breakdown voltage for quasi-homogeneous sphere-plane gaps | |
| Zhang et al. | Calculation of ion flow fields of HVDC transmission lines by coupling flux tracing method with finite element method | |
| Li et al. | Research on calculation method of ion flow field under multi-circuit HVDC transmission lines |
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 |