CN109711088A - A kind of distribution line lightning induced voltage calculation method of Electromagnetic field - Google Patents
A kind of distribution line lightning induced voltage calculation method of Electromagnetic field Download PDFInfo
- Publication number
- CN109711088A CN109711088A CN201910033658.5A CN201910033658A CN109711088A CN 109711088 A CN109711088 A CN 109711088A CN 201910033658 A CN201910033658 A CN 201910033658A CN 109711088 A CN109711088 A CN 109711088A
- Authority
- CN
- China
- Prior art keywords
- lightning
- current
- model
- line
- indicate
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention provides a kind of distribution line lightning induced voltage calculation method of Electromagnetic field, the external electromagnetic field that lightning stroke generates in space is calculated using lightning return stroke model and Cooray-Rubinstein formula, further by considering Agrawal field wire coupling model, distribution line lightning induced voltage is calculated using FDTD method, improve calculating speed, and the influence of non-ideal the earth, leader process and corona phenomenon is considered, can be preferably applied in Practical Project.
Description
Technical field
The present invention relates to distribution line lightning protection field, the distribution line inductive lightning of specifically a kind of Electromagnetic field is excessively electric
Press calculation method.
Background technique
The transmission line of electricity and distribution line of China's electric system mostly use overhead line to transmit electric energy, are easily struck by lightning
To cause power failure, damage power supply unit and household electrical appliance etc..For voltage class from 6kV to 35kV power distribution network, generally
There is no ground wire protection and line insulation level is lower, in addition network structure is also intricate, therefore not only direct lightning strike can cause thunder
Evil accident, inductive lightning can also cause biggish harm.According to the operation investigation to power supply department, in distribution line, lightning stroke
Trip-out rate accounts for 80% or more of its total trip-out rate, and line tripping caused by lightning induced voltage accounts for majority among these.
The first half in 20th century, when illustrating thunder and lightning induction voltage Producing reason, most commonly seen theory is constraint
Charge theory；Between nineteen sixty to 1980, people think that guide's development process of thunder and lightning is to lead to distribution line overvoltage again
The main reason for, the lightning Protection Design of this theoretical distribution line to Japan produces very big influence.But current research
Show both theoretical corresponding lightning induced voltage ingredient very littles.With the gradually development of science and technology, measuring instrument also becomes
Must be more and more accurate, this provides the experimental data of high-quality to probe into the mechanism of thunder and lightning induction voltage generation, to be born
Thunder and lightning induction voltage based on Fields of Lightning Return Stroke process is theoretical.The present age thunder and lightning induction voltage theory think, non-direct lightning strike
What is generated, which threatens the overvoltage main component of distribution line insulation, is generated during the counterattack of thunder and lightning.In addition, in recent years
The electric field of further observation to generate to(for) lightning current shows when observation point is closer apart from thunder and lightning circulation road, fights back every time
Leader process before process can also generate biggish electric field, and then generate apparent thunder and lightning induction voltage on the line, no
The leader process for the entire lightning stroke process generation before crossing leader process here and referring to before is not same concept, development
Speed is also only than fighting back slow 10 times of process or so.
After the producing cause of clear Fields of Lightning Return Stroke process, the calculation method for the external electromagnetic field that thunder and lightning generates also is ground therewith
Study carefully and increases；Further, the coupling model between external electromagnetic field and open-wire transmission line, i.e. field wire coupling model, and therewith
Corresponding method for solving is there has also been more rapidly development, and for current, FDTD (finite time-domain calculus of finite differences) is commonly to solve
Method.
Summary of the invention
The purpose of the present invention is to propose to a kind of distribution line lightning induced voltage calculation methods of Electromagnetic field, utilize thunder and lightning
Return-stroke model calculates the external electromagnetic field that lightning stroke generates in space, further by considering Agrawal field wire coupling model, uses
Distribution line lightning induced voltage is calculated in FDTD method, improves calculating speed, and consider non-ideal the earth, guide's mistake
The influence of journey and corona phenomenon can preferably be applied in Practical Project.
A kind of distribution line lightning induced voltage calculation method of Electromagnetic field, comprising the following steps:
(1) waveform of lightning current the bottom of the channel return stroke current used by determining；
(2) the suitable lightning current of the waveform selection of the lightning current the bottom of the channel return stroke current determined according to step (1) is fought back
Engineering model；
(3) lightning current determined according to step (2) fights back engineering model and obtains the relevant parameter of thunder and lightning circulation road, and then makes
The external electromagnetic field that lightning current generates is calculated with maxwell equation group；
(4) Agrawal method is used, field wire coupled problem is handled as electromagnetic scattering problems, tangential electric field on transmission line
Intensity generates spread voltage source, and then determines Agrawal model equation；
(5) Agrawal model equation is solved using FDTD method, obtains the thunder and lightning induction voltage on open-wire transmission line
Value.
Further, the waveform of lightning current the bottom of the channel return stroke current includes triangle wave pattern, two fingers in the step (1)
Exponential model and Heidler model, wherein Heidler model expression be
I in formulapeakIt is the peak value of lightning current, τ1It is related to the rate of climb of lightning current, and τ2With the rate of decay of lightning current
Steepness before correlation, n and lightning current wave is related, and the wavefront of n more large lightning current is steeper, and η is current peak correction factor, generally takes
Further, it includes transmission line model (TL), current source that the lightning current that the step (2) determines, which fights back engineering model,
Propagation model (TCS), transmission line improved model (MTL), Bruce-Golde model (BG), each model use a unified public affairs
Formula is expressed as follows:
Wherein z is height of the observed point apart from ground, u (t-z/vf) it is Heaviside function, as t > z/vfWhen u=
1, otherwise u=0, P (z) are corrected parameter related with height z, vfIt is the development speed of thunder and lightning circulation road, v is electric current in channel
Spread speed.
Further, the external electromagnetic field that lightning current generates is calculated using maxwell equation group in the step (3) to have
Body are as follows:
When considering ideal the earth, vertical electric field EZ, horizontal component of electric field EP, horizontal magnetic field HPApproximate formula be respectively
Wherein r is to be observed a little the distance between with thunder and lightning circulation road, and z is the height being observed a little, z0For thunder and lightning circulation
The height of Dao Shangzhen's that point calculated, R, which is positive, that point calculated and to be observed the distance between a little, and H is lightning current
The total length in channel；
The case where for non-ideal the earth, it is only necessary to influence of the non-ideal the earth to horizontal component of electric field is considered, using Cooray
The calculation method of horizontal component of electric field caused by the non-ideal the earth lightning current of calculating proposed with Rubinstein, such as following formula institute
Show:
Wherein, σgIndicate the conductivity of the earth, εrgIndicate the relative permitivity of the earth.
Further, the Agrawal model equation in the step (4) is expressed as follows:
Wherein h indicates height of the route apart from ground, ErIndicate the electric field strength of horizontal direction, L and C indicate the electricity of route
Sense and capacitor；
The total voltage u that lightning current incudes generation on the line is calculated using following formula
Wherein us(x, t) is the spread voltage that horizontal component of electric field generates,It is that vertical electric field produces
Raw induced voltage；
Above calculating is carried out in the case where ideal the earth, when considering the earth and line loss,
Expression formula of the Agrawal model equation under time domain becomes again:
Wherein G indicates the conductance of route, and ε (t) indicates the transient resistance of soil, they are counted using following formula
It calculates:
Wherein
τg=h2μ0σg
σ in formulaairIndicate the conductivity of air, F-1(jw) inverse Fourier transform is indicated.
Further, the step (5) carries out solving Agrawal model equation using FDTD method specifically:
First spatially for voltage and current progress discrete processes, the voltage and current node chosen on overhead line,
It is required to meet design conditions Δ x >=v Δ t, wherein v is electromagnetic field in the speed for calculating spatial, in most cases
Take v=c；
Voltage and current distribution along known n time Δt, it is desirable that the distribution along (n+1) time Δt voltage and current, it is first
The voltage value at route both ends is first calculated according to the following formula:
Then route intermediate node is calculated in the voltage value of (n+1) time Δt:
Finally calculate the current distribution on (n+1) time Δt route:
Wherein HkIndicate height of the route apart from ground, G at distance line head end (k-1) Δ x1Indicate route head end
Port conductance, GN+1Indicate the port conductance of line end,Indicate that distance line head end (k-1) Δ x is in n time Δt
Vertical electric field strength,It indicates to be in the horizontal component of electric field intensity of (n+1/2) time Δt in distance line head end (k-1/2) Δ x
Component in direction along the line,Indicate that distance line head end (k-1) Δ x is in the voltage value of n time Δt,It indicates in distance
Route head end (k-1/2) Δ x is in the current value of (n+1/2) time Δt, and Δ x is the section unit length of route, C and L difference
Indicate the capacitance matrix and inductance matrix of route.
Further, when considering guide's development process, the model of guide's development is as follows:
U (t) is Heaviside function in formula, and δ (t) is Dirac function, and integral in the time domain is equal to 1, and works as t ≠ 0
When δ (t)=0, v indicate guide's speed of development, ρ0It is charge density constant.
Further, when considering corona phenomenon, originally changeless line capacitance C becomes a dynamic capacity Cdyn,
CdynIt is calculated using formula below:
Cdyn(x, t)=Cwhen (u (x, t) < uth(x, t))
Wherein uthThe discharge inception voltage for indicating conducting wire, according to bloom electric field strength EthIt is calculated, in practical projects generally
Take Eth=29kV/cm, k1、k2It is the parameter for needing to set in this dynamic model, k1It indicates that corona phenomenon, capacitor once occurs
Increase k1Times, k2It then indicates also to will continue to increase with the lasting enhancing of corona phenomenon, capacitance.
Distribution line lightning induced voltage calculation method proposed by the present invention has following excellent compared with original algorithm
Point:
1. the present invention calculates non-ideal the earth lightning current institute using lightning return stroke model and Cooray-Rubinstein formula
The horizontal component of electric field of generation, accurate and calculation amount is small, and reference is facilitated to use existing model parameter.
2. lightning induced voltage calculation method proposed by the present invention uses Agrawal field wire coupling model and the side FDTD
Method can quickly calculate distribution line lightning induced voltage, and consider non-ideal the earth, leader process and corona phenomenon
It influences, is suitably applied Practical Project.
Detailed description of the invention
Fig. 1 is the Electromagnetic Calculation schematic diagram that lightning current generates；
Fig. 2 is the equivalent circuit diagram of Agrawal field wire coupling model；
Fig. 3 is FDTD voltage and current node selection schematic diagram；
Fig. 4 is FDTD schematic diagram of calculation flow；
Fig. 5 is to consider guide and the electromagnetic field schematic diagram that counterattack process generates simultaneously.
Specific embodiment
Below in conjunction with the attached drawing in the present invention, the technical solution in the present invention is clearly and completely described.
Distribution line lightning induced voltage calculation method proposed by the present invention, it is main comprising based on lightning current return-stroke model
The FDTD of spatial electromagnetic field computation and Agrawal field wire coupling model solves two large divisions, comprising the following steps:
(1) waveform of lightning current the bottom of the channel return stroke current used by determining.In actual measurement, in existing technology item
It under part, is only possible to measure the counterattack lightning current waveform of lightning channel bottom, so existing lightning current is fought back in engineering model, return
The spatial and temporal distributions calculating for hitting lightning current is all based on the counterattack lightning current of the bottom of the channel.Currently used lightning current channel floor
Portion's return stroke current is triangle wave pattern, biexponential model and Heidler model.Since Heidler model is more in line with people
For the understanding and expection of lightning current, and there is no biexponential model derivative is not zero at zero point the case where, be more suitable use
In simulation calculation, therefore more recommended use Heidler model, expression formula are
I in formulapeakIt is the peak value of lightning current, τ1It is related to the rate of climb of lightning current, and τ2With the rate of decay of lightning current
Steepness before correlation, n and lightning current wave is related, and the wavefront of n more large lightning current is steeper, and η is current peak correction factor, generally
It takes
(2) the suitable lightning current of the waveform selection of the lightning current the bottom of the channel return stroke current determined according to step (1) is fought back
Engineering model.It is more famous at present to have: transmission line model (TL), current source propagation model (TCS), transmission line improved model
(MTL), Bruce-Golde model (BG), Diendorfer-Uman model (DU) and Diendorfer-Uman improved model
(MDU), wherein MTL model can be subdivided into again as MTLL model and MTLE model, these engineering models are to a certain degree
On can simulate actual lightning current well.
A unified formula expression generally can be used in these engineering models (other than DU model and MDU model)
It is as follows:
Wherein z is height of the observed point apart from ground, u (t-z/vf) it is Heaviside function, as t > z/vfWhen u=
1, otherwise u=0, P (z) are corrected parameter related with height z, vfIt is the development speed of thunder and lightning circulation road, v is electric current in channel
Spread speed, the value of the P (z) and v of various models is as shown in the table.
The value of table 1 various model P (z) and v
(3) lightning current determined according to step (2) fights back engineering model and obtains the relevant parameter of thunder and lightning circulation road, and then makes
The external electromagnetic field that lightning current generates is calculated with maxwell equation group.Usually regard thunder and lightning circulation road as one stand upright on
The antenna of big ground, considers its mirror image, using maxwell equation group can in the hope of the magnetic distribution around lightning current,
Wherein the relevant parameter of thunder and lightning circulation road determines that lightning current return-stroke model is provided by step (2).When considering ideal the earth, hang down
Straight electric field EZ, horizontal component of electric field EP, horizontal magnetic field HPApproximate formula be respectively
As shown in Figure 1, wherein r is to be observed a little the distance between with thunder and lightning circulation road, z is the height being observed a little, z0
For the height of that point just calculated on thunder and lightning circulation road, R, which is positive, that point calculated and to be observed the distance between a little,
H is the total length of thunder and lightning circulation road.
And the case where for non-ideal the earth, mirror image lightning current will generate distortion, and in the lower electric current of different frequency
Distortion degree is different and related with the conductivity of the earth and electric current source position.Traditional Sommerfeld method is on frequency domain
Calculation method, since thunder and lightning frequency spectrum is wide, so calculation amount is very big, the degree that can permit far beyond engineering.Actual
In engineering calculation, it is believed that the wherein electric field E in direction perpendicular to the groundZIt is little by being influenced, therefore only need to consider non-ideal big
Influence of the ground to horizontal component of electric field.Caused by the non-ideal the earth lightning current of calculating proposed using Cooray and Rubinstein
The calculation method of horizontal component of electric field, is shown below.
Wherein, σgIndicate the conductivity of the earth, εrgIndicate the relative permitivity of the earth.
(4) Agrawal method is used, field wire coupled problem is handled as electromagnetic scattering problems, tangential electric field on transmission line
Intensity generates spread voltage source, and then determines Agrawal model equation.Agrawal model equation can be expressed as follows:
Wherein h indicates height of the route apart from ground, ErIndicate the electric field strength of horizontal direction, L and C indicate the electricity of route
Sense and capacitor.
The total voltage u that lightning current incudes generation on the line can be used following formula and be calculated
Wherein us(x, t) is the spread voltage that horizontal component of electric field generates,It is that vertical electric field produces
Raw induced voltage.The equivalent circuit diagram of Agrawal model is as shown in Figure 2.
Above calculating is carried out in the case where ideal the earth, when considering the earth and line loss,
Expression formula of the Agrawal model under time domain becomes again:
Wherein G indicates the conductance of route, and ε (t) indicates the transient resistance of soil, they can be used following formula and carry out
It calculates:
Wherein
τg=h2μ0σg
σ in formulaairIndicate the conductivity of air, F-1(jw) inverse Fourier transform is indicated.
(5) Agrawal model equation is solved using FDTD (finite time-domain calculus of finite differences) method, obtained on open-wire transmission line
Thunder and lightning induction voltage value.Above-mentioned equation needs to be solved using certain numerical computation method, and is directed to such calculating
Scale is generally solved using FDTD method.
Discrete processes spatially are carried out for voltage and current first, the voltage and current node chosen on overhead line is such as
Shown in Fig. 3, it is desirable that meet design conditions Δ x >=v Δ t, wherein v is electromagnetic field in the speed for calculating spatial, absolutely mostly
V=c is taken in number situation.
If Fig. 4 is the schematic diagram of FDTD calculation process, it is known that the voltage and current distribution along n time Δt, it is desirable that (n+1)
Distribution along time Δt voltage and current can calculate the voltage value at route both ends according to the following formula first.
Then route intermediate node is calculated in the voltage value of (n+1) time Δt.
Finally (current distribution for calculating current time is to calculate to the current distribution on calculating (n+1) time Δt route
The voltage's distribiuting of subsequent time).
Wherein HkIndicate height of the route apart from ground, G at distance line head end (k-1) Δ x1Indicate route head end
Port conductance, GN+1Indicate the port conductance of line end,Indicate that distance line head end (k-1) Δ x is in n time Δt
Vertical electric field strength,It indicates to be in the horizontal component of electric field intensity of (n+1/2) time Δt in distance line head end (k-1/2) Δ x
Component in direction along the line,Indicate that distance line head end (k-1) Δ x is in the voltage value of n time Δt,It indicates in distance
Route head end (k-1/2) Δ x is in the current value of (n+1/2) time Δt, and Δ x is the section unit length of route, C and L difference
Indicate the capacitance matrix and inductance matrix of route.
On the basis of above, guide's development process is further considered.As it was noted above, the elder generation before each counterattack process
Biggish electric field can also be generated by leading process, therefore can increase consideration leader process in simulation calculation.Basic ideas such as Fig. 5 institute
Show, has accumulated the charge that total charge dosage is Q in thundercloud first, the guide of formation is with constant speed vLDevelopment downwards reaches in ground
After face, occurrence and development speed fights back process faster immediately.The model of guide's development is as follows:
U (t) is Heaviside function in formula, and δ (t) is Dirac function, and integral in the time domain is equal to 1, and works as t ≠ 0
When δ (t)=0, v indicate guide's speed of development, ρ0It is charge density constant.
On the basis of above, corona phenomenon is considered.When amplitude of lightning current is very big or lightning current channel distance overhead line
It is close, biggish induced overvoltage will be generated on the line, it is possible to generate corona phenomenon, in this case, originally
Changeless line capacitance C has reformed into a dynamic capacity Cdyn, CdynFormula below can be used to be calculated.
Cdyn(x, t)=Cwhen (u (x, t) < uth(x, t))
Wherein uthThe discharge inception voltage for indicating conducting wire, can be according to bloom electric field strength EthIt is calculated, in practical projects
Generally take Eth=29kV/cm, k1、k2It is the parameter for needing to set in this dynamic model, k1It indicates that corona phenomenon once occurs,
Capacitor just will increase k1Times, k2It then indicates also to will continue to increase with the lasting enhancing of corona phenomenon, capacitance.Consider that corona is existing
As when use the changeless line capacitance C in above-mentioned dynamic capacity replacement step (4).
The present invention calculates non-ideal the earth lightning current using lightning return stroke model and Cooray-Rubinstein formula and is produced
Raw horizontal component of electric field, accurate and calculation amount is small, and reference is facilitated to use existing model parameter；The lightning induced voltage calculating side of proposition
Method uses Agrawal field wire coupling model and FDTD method, can quickly calculate distribution line lightning induced voltage, and consider
The influence of non-ideal the earth, leader process and corona phenomenon, is suitably applied Practical Project.
Claims (8)
1. a kind of distribution line lightning induced voltage calculation method of Electromagnetic field, it is characterised in that the following steps are included:
(1) waveform of lightning current the bottom of the channel return stroke current used by determining；
(2) the suitable lightning current of the waveform selection of the lightning current the bottom of the channel return stroke current determined according to step (1) fights back engineering
Model；
(3) lightning current determined according to step (2) fights back engineering model and obtains the relevant parameter of thunder and lightning circulation road, and then uses wheat
Ke Siwei equation group calculates the external electromagnetic field that lightning current generates；
(4) Agrawal method is used, field wire coupled problem is handled as electromagnetic scattering problems, tangential electric field strength on transmission line
Spread voltage source is generated, and then determines Agrawal model equation；
(5) Agrawal model equation is solved using FDTD method, obtains the thunder and lightning induction voltage value on open-wire transmission line.
2. the distribution line lightning induced voltage calculation method of Electromagnetic field as described in claim 1, it is characterised in that: described
The waveform of lightning current the bottom of the channel return stroke current includes triangle wave pattern, biexponential model and Heidler mould in step (1)
Type, wherein Heidler model expression be
I in formulapeakIt is the peak value of lightning current, τ1It is related to the rate of climb of lightning current, and τ2With the rate of decay phase of lightning current
It closes, n is related with the steepness before lightning current wave, and the wavefront of n more large lightning current is steeper, and η is current peak correction factor, generally takes
3. the distribution line lightning induced voltage calculation method of Electromagnetic field as described in claim 1, it is characterised in that: described
It includes transmission line model (TL), current source propagation model (TCS), transmission line that the lightning current that step (2) determines, which fights back engineering model,
Improved model (MTL), Bruce-Golde model (BG), each model are expressed as follows using a unified formula:
Wherein z is height of the observed point apart from ground, u (t-z/vf) it is Heaviside function, as t > z/vfWhen u=1, it is no
Then u=0, P (z) are corrected parameter related with height z, vfIt is the development speed of thunder and lightning circulation road, v is the biography of electric current in channel
Broadcast speed.
4. the distribution line lightning induced voltage calculation method of Electromagnetic field as described in claim 1, it is characterised in that: described
The external electromagnetic field that lightning current generates is calculated using maxwell equation group in step (3) specifically:
When considering ideal the earth, vertical electric field EZ, horizontal component of electric field EP, horizontal magnetic field HPApproximate formula be respectively
Wherein r is to be observed a little the distance between with thunder and lightning circulation road, and z is the height being observed a little, z0For on thunder and lightning circulation road just
The height of that point calculated, R, which is positive, that point calculated and to be observed the distance between a little, and H is thunder and lightning circulation road
Total length；
The case where for non-ideal the earth, it is only necessary to consider influence of the non-ideal the earth to horizontal component of electric field, using Cooray and
The calculation method of horizontal component of electric field caused by the non-ideal the earth lightning current of calculating that Rubinstein is proposed, is shown below:
Wherein, σgIndicate the conductivity of the earth, εrgIndicate the relative permitivity of the earth.
5. the distribution line lightning induced voltage calculation method of Electromagnetic field as described in claim 1, it is characterised in that: described
Agrawal model equation in step (4) is expressed as follows:
Wherein h indicates height of the route apart from ground, ErIndicate the electric field strength of horizontal direction, L and C indicate route inductance and
Capacitor；
The total voltage u that lightning current incudes generation on the line is calculated using following formula
Wherein us(x, t) is the spread voltage that horizontal component of electric field generates,It is that vertical electric field generates
Induced voltage；
Above calculating is carried out in the case where ideal the earth, when considering the earth and line loss, Agrawal mould
Expression formula of the type equation under time domain becomes again:
Wherein G indicates the conductance of route, and ε (t) indicates the transient resistance of soil, they are calculated using following formula:
Wherein
σ in formulaairIndicate the conductivity of air, F-1(jw) inverse Fourier transform is indicated.
6. the distribution line lightning induced voltage calculation method of Electromagnetic field as described in claim 1, it is characterised in that: described
Step (5) carries out solving Agrawal model equation using FDTD method specifically:
Discrete processes, the voltage and current node chosen on overhead line spatially are carried out for voltage and current first, it is desirable that
Meet design conditions Δ x >=v Δ t, wherein v is electromagnetic field in the speed for calculating spatial, takes v=in most cases
c；
Voltage and current distribution along known n time Δt, it is desirable that the distribution along (n+1) time Δt voltage and current, first root
The voltage value at route both ends is calculated according to following formula:
Then route intermediate node is calculated in the voltage value of (n+1) time Δt:
Finally calculate the current distribution on (n+1) time Δt route:
Wherein HkIndicate height of the route apart from ground, G at distance line head end (k-1) Δ x1Indicate the port of route head end
Conductance, GN+1Indicate the port conductance of line end,Indicate that distance line head end (k-1) Δ x is in the vertical of n time Δt
Electric field strength,It indicates to be in the horizontal component of electric field intensity of (n+1/2) time Δt on edge in distance line head end (k-1/2) Δ x
The component in line direction,Indicate that distance line head end (k-1) Δ x is in the voltage value of n time Δt,It indicates in distance line
Head end (k-1/2) Δ x is in the current value of (n+1/2) time Δt, and Δ x is the section unit length of route, and C and L are respectively indicated
The capacitance matrix and inductance matrix of route.
7. the distribution line lightning induced voltage calculation method of Electromagnetic field as described in claim 1, it is characterised in that: consider
When guide's development process, the model of guide's development is as follows:
U (t) is Heaviside function in formula, and δ (t) is Dirac function, and integral in the time domain is equal to 1, and the δ as t ≠ 0
(t)=0, v indicates guide's speed of development, ρ0It is charge density constant.
8. the distribution line lightning induced voltage calculation method of Electromagnetic field as described in claim 1, it is characterised in that: consider
When corona phenomenon, originally changeless line capacitance C becomes a dynamic capacity Cdyn, CdynIt is counted using formula below
It calculates:
cdyn(x, t)=C when (u (x, t) < uth(x, t))
Wherein uthThe discharge inception voltage for indicating conducting wire, according to bloom electric field strength EthIt is calculated, generally takes E in practical projectsth
=29kV/cm, k1、k2It is the parameter for needing to set in this dynamic model, k1Indicate that capacitor increases once corona phenomenon occurs
k1Times, k2It then indicates also to will continue to increase with the lasting enhancing of corona phenomenon, capacitance.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910033658.5A CN109711088B (en) | 2019-01-15 | 2019-01-15 | Distribution line inductive lightning overvoltage calculation method combining field and line |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910033658.5A CN109711088B (en) | 2019-01-15 | 2019-01-15 | Distribution line inductive lightning overvoltage calculation method combining field and line |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109711088A true CN109711088A (en) | 2019-05-03 |
| CN109711088B CN109711088B (en) | 2022-12-23 |
Family
ID=66261329
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910033658.5A Active CN109711088B (en) | 2019-01-15 | 2019-01-15 | Distribution line inductive lightning overvoltage calculation method combining field and line |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109711088B (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110516362A (en) * | 2019-08-28 | 2019-11-29 | 哈尔滨工程大学 | A kind of FDTD method for solving of multi-conductor transmission lines magnetic distribution |
| CN112036005A (en) * | 2020-07-21 | 2020-12-04 | 广东电网有限责任公司广州供电局 | Method and device for calculating lightning induction voltage |
| CN112131736A (en) * | 2020-09-21 | 2020-12-25 | 中国电子科技网络信息安全有限公司 | Lightning electromagnetic environment simulation modeling method |
| CN112364472A (en) * | 2019-07-26 | 2021-02-12 | 天津大学青岛海洋技术研究院 | Matlab-based transient electromagnetic field propagation simulation method |
| CN112749500A (en) * | 2021-01-11 | 2021-05-04 | 华能新能源股份有限公司 | Method for estimating lightning induction overvoltage of current collecting circuit of wind power plant in high-altitude mountain area |
| CN113255975A (en) * | 2021-05-12 | 2021-08-13 | 云南电网有限责任公司丽江供电局 | Spatial position optimization method for distribution line erection overhead ground wire |
| CN113252971A (en) * | 2021-05-12 | 2021-08-13 | 云南电网有限责任公司丽江供电局 | FDTD-based distribution line inductive lightning overvoltage calculation method |
| CN113255976A (en) * | 2021-05-12 | 2021-08-13 | 云南电网有限责任公司丽江供电局 | Grounding optimization design method for distribution line erection overhead ground wire |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103245826A (en) * | 2013-05-16 | 2013-08-14 | 广州供电局有限公司 | Lightning induction voltage determining method and system on distributing lines |
| JP2017181039A (en) * | 2016-03-28 | 2017-10-05 | 一般財団法人電力中央研究所 | Estimation method, estimation device, and estimation program for lightening strike current characteristics value |
-
2019
- 2019-01-15 CN CN201910033658.5A patent/CN109711088B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103245826A (en) * | 2013-05-16 | 2013-08-14 | 广州供电局有限公司 | Lightning induction voltage determining method and system on distributing lines |
| JP2017181039A (en) * | 2016-03-28 | 2017-10-05 | 一般財団法人電力中央研究所 | Estimation method, estimation device, and estimation program for lightening strike current characteristics value |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112364472A (en) * | 2019-07-26 | 2021-02-12 | 天津大学青岛海洋技术研究院 | Matlab-based transient electromagnetic field propagation simulation method |
| CN110516362A (en) * | 2019-08-28 | 2019-11-29 | 哈尔滨工程大学 | A kind of FDTD method for solving of multi-conductor transmission lines magnetic distribution |
| CN112036005A (en) * | 2020-07-21 | 2020-12-04 | 广东电网有限责任公司广州供电局 | Method and device for calculating lightning induction voltage |
| CN112131736A (en) * | 2020-09-21 | 2020-12-25 | 中国电子科技网络信息安全有限公司 | Lightning electromagnetic environment simulation modeling method |
| CN112749500A (en) * | 2021-01-11 | 2021-05-04 | 华能新能源股份有限公司 | Method for estimating lightning induction overvoltage of current collecting circuit of wind power plant in high-altitude mountain area |
| CN112749500B (en) * | 2021-01-11 | 2023-09-22 | 华能新能源股份有限公司 | Estimation method for lightning induction overvoltage of wind power plant collector line in high-altitude mountain area |
| CN113255975A (en) * | 2021-05-12 | 2021-08-13 | 云南电网有限责任公司丽江供电局 | Spatial position optimization method for distribution line erection overhead ground wire |
| CN113252971A (en) * | 2021-05-12 | 2021-08-13 | 云南电网有限责任公司丽江供电局 | FDTD-based distribution line inductive lightning overvoltage calculation method |
| CN113255976A (en) * | 2021-05-12 | 2021-08-13 | 云南电网有限责任公司丽江供电局 | Grounding optimization design method for distribution line erection overhead ground wire |
| CN113255976B (en) * | 2021-05-12 | 2023-01-24 | 云南电网有限责任公司丽江供电局 | Grounding optimization design method for distribution line erection overhead ground wire |
| CN113255975B (en) * | 2021-05-12 | 2023-10-13 | 云南电网有限责任公司丽江供电局 | Space position optimization method for overhead ground wire erected on distribution line |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109711088B (en) | 2022-12-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109711088A (en) | A kind of distribution line lightning induced voltage calculation method of Electromagnetic field | |
| Grcev | Time-and frequency-dependent lightning surge characteristics of grounding electrodes | |
| Grcev | Modeling of grounding electrodes under lightning currents | |
| Javor et al. | A channel-base current function for lightning return-stroke modeling | |
| CN103245826B (en) | Distribution line lightning induced voltage determines method and system | |
| CN105160049B (en) | A kind of 35kV and following distribution line direct lightning strike trip-out rate computational methods | |
| JP5282478B2 (en) | Distribution line direct lightning strike rate calculation device, distribution line direct lightning strike rate calculation method, distribution line lightning strike rate calculation device, distribution line lightning strike rate calculation method, and program | |
| CN104850738B (en) | Lightning induced voltage computational methods at overhead power transmission line line pole tower | |
| Cao et al. | Lightning surge analysis of transmission line towers with a hybrid FDTD-PEEC method | |
| Shariatinasab et al. | The effect of grounding system modeling on lightning-related studies of transmission lines | |
| Araneo et al. | Utilization of underbuilt shield wires to improve the lightning performance of overhead distribution lines hit by direct strokes | |
| Alipio et al. | Analysis of the lightning impulse and low-frequency performance of wind farm grounding systems | |
| Tossani et al. | Estimation of the influence of direct strokes on the lightning performance of overhead distribution lines | |
| Malicki et al. | Shielding failure analysis of a hybrid transmission line with AC and DC systems on the same tower | |
| Visacro et al. | Measurements for qualifying the lightning response of tower-footing electrodes of transmission lines | |
| Alipio et al. | Influence of a lossy ground on the lightning performance of overhead transmission lines | |
| Javor et al. | Analytical representation of measured lightning currents and its application to electromagnetic field estimation | |
| Piantini et al. | Lightning interaction with medium-voltage overhead power distribution systems | |
| Vasconcellos et al. | The influence of the frequency-dependent behavior of ground electrical parameters on the lightning performance of transmission lines | |
| Martinez et al. | Influence of the stroke angle on the flashover rate of an overhead transmission line | |
| Malicki et al. | Estimating the lightning performance of a multi-circuit transmission tower | |
| Ioannidis et al. | Estimating the shielding failure flashover rate of single-circuit overhead lines with horizontal phase configuration via stochastic lightning attachment simulations | |
| Alípio et al. | Lightning Overvoltages Including Frequency-Dependent Soil Parameters in the Transmission Line Model | |
| Lertwanitrot et al. | Characteristics and behaviour of transient current during lightning strikes on transmission tower | |
| Zhang et al. | Validation and revision of far-field-current relationship for the lightning strike to electrically short objects |
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 |