CN100361123C - Whole line and multi-parameter integrated optimizing method for determining lightningproof performance of power transmission line - Google Patents
Whole line and multi-parameter integrated optimizing method for determining lightningproof performance of power transmission line Download PDFInfo
- Publication number
- CN100361123C CN100361123C CNB2005101166946A CN200510116694A CN100361123C CN 100361123 C CN100361123 C CN 100361123C CN B2005101166946 A CNB2005101166946 A CN B2005101166946A CN 200510116694 A CN200510116694 A CN 200510116694A CN 100361123 C CN100361123 C CN 100361123C
- Authority
- CN
- China
- Prior art keywords
- lightning
- parameter
- shaft tower
- tripping rate
- lightning protection
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000005540 biological transmission Effects 0.000 title claims abstract description 38
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 12
- 238000005457 optimization Methods 0.000 claims abstract description 9
- 230000005611 electricity Effects 0.000 claims description 26
- 238000013461 design Methods 0.000 abstract description 8
- 230000008901 benefit Effects 0.000 abstract description 3
- APTZNLHMIGJTEW-UHFFFAOYSA-N pyraflufen-ethyl Chemical compound C1=C(Cl)C(OCC(=O)OCC)=CC(C=2C(=C(OC(F)F)N(C)N=2)Cl)=C1F APTZNLHMIGJTEW-UHFFFAOYSA-N 0.000 description 27
- 238000004364 calculation method Methods 0.000 description 13
- 238000009413 insulation Methods 0.000 description 13
- 230000008878 coupling Effects 0.000 description 11
- 238000010168 coupling process Methods 0.000 description 11
- 238000005859 coupling reaction Methods 0.000 description 11
- 239000012212 insulator Substances 0.000 description 11
- 208000025274 Lightning injury Diseases 0.000 description 10
- 230000000694 effects Effects 0.000 description 8
- 238000009826 distribution Methods 0.000 description 7
- 238000005070 sampling Methods 0.000 description 7
- 238000012876 topography Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 239000004567 concrete Substances 0.000 description 4
- 239000004020 conductor Substances 0.000 description 3
- 241000209094 Oryza Species 0.000 description 2
- 235000007164 Oryza sativa Nutrition 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 208000028659 discharge Diseases 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 235000009566 rice Nutrition 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000002459 sustained effect Effects 0.000 description 2
- 230000009885 systemic effect Effects 0.000 description 2
- 241000282326 Felis catus Species 0.000 description 1
- 244000080575 Oxalis tetraphylla Species 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000011150 reinforced concrete Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005309 stochastic process Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Landscapes
- Suspension Of Electric Lines Or Cables (AREA)
Abstract
The present invention relates to a whole line, multiple parameter and integrated optimization method for judging lightning protection performance of a transmission line, which belongs to the technical field of protecting lightning for a transmission line in an electric power system. The method comprises: selecting a parameter of each rod of a transmission line which is required to be judged; using the Monte Carlo method for calculating the tower tripping rate with lightning strikes of each rod; then calculating the tripping rate with lightning strikes of the line; when the tripping rate is smaller than a set value which can satisfy the lightning protection performance requirements, the lightning protection performance is in accordance with the requirements, the judgment is end, otherwise towers which have high tripping rates are found and regarded as towers which have poor lightning protection performance; adjustable parameters are selected; a plurality of points are respectively selected for calculating the tripping rate with lightning strikes of each parameter, and a multiple parameter integrated optimization curve is matched; a suitable optimization point is selected, the tripping rate with lightning strikes of the line is recalculated till the performance requirements of lightning protection is satisfied. The present invention can enable lightning protection measures of a transmission line to have the advantages of pertinence, scientificalness and systematicality in the phases of design and operation.
Description
Technical field
The invention belongs to power system transmission line shocking preventing technology field, be particularly related to a kind of method that is used to judge the transmission line of electricity lightning protection properties, this method can be carried out the emulation of all fronts road, multi-parameter integrated optimizing to the lightning protection properties of transmission line of electricity, in order to the lightning protection work in scientific guidance transmission line of electricity design and running stage.
Background technology
According to power department statistics, thunderbolt is the major failure type of transmission line of electricity, and wherein the line is busy that road tripping operation ratio can be up to 60% because of lightning stroke trip.Therefore, no matter still move department for line design department, analyzing the quality of judging the transmission line of electricity lightning protection properties all is a very important job, and the leading indicator of reflection transmission line of electricity lightning protection properties is the tripping rate with lightning strike of circuit.In actual lightening protection engineering, generally be according to concrete line parameter circuit value such as shaft tower model, stake resistance, insulation configuration etc., select suitable tripping rate with lightning strike method for use, calculate the theoretical value of this lightning outage rate, be higher than setting value if find this value, then need adopt corresponding lightning protection measures, make tripping rate with lightning strike be controlled in the acceptable scope.
At present, transmission line of electricity being carried out tripping rate with lightning strike, to calculate general in the world be the method that the method introduced in " Guidefor Improving the Lighting Performance of the Transmission Line " of IEEE1997 and CIGRE1991 introduce in " Guide to Procedures for Estimating the LightningPerformance of Transmission Lines ".
Utilize above-mentioned two kinds of methods to judge that the basic ideas of transmission line of electricity lightning protection properties method are similarly, promptly represent the lightning protection properties of circuit with the tripping rate with lightning strike of typical shaft tower in the circuit, idiographic flow as shown in Figure 1, each step is described below:
1) at first to select representative typical shaft tower (go back neither one science and unified index system at present, and more will rely on on-site experience) for the measurement of typical shaft tower;
2) choose required shaft tower parameter (being provided by power department), parameter mainly comprises the topography and geomorphology of shaft tower major dimension, stake resistance, insulation configuration, position etc.;
3) with in a series of analytical expressions in aforementioned IEEE of parameter substitution or the CIGRE computing method, final calculation result is the tripping rate with lightning strike value of this shaft tower;
4) with the foundation of this tripping rate with lightning strike value as the quality of the lightning protection properties of judging the whole piece circuit;
5) if this trip-out rate value less than the setting value that satisfies the lightning protection properties requirement, then this line thunder protection performance is qualified, then skips to step 7), judge to finish, otherwise enters next step.
6) adjusting some influences the shaft tower parameter of lightning protection properties or installs lightning protection measures additional, jumps to step 2)
7) judge end.
From top step as can be seen, this method to the analysis of line thunder protection performance based on single shaft tower, though method is simplified, but in actual track, numerous condition difference such as topography and geomorphology at the model of each shaft tower, place are very big along the line, by this short-cut method analysis all fronts road lightning protection properties, not only computational accuracy can't be guaranteed, and the lightning protection properties of each section of road, each basic shaft tower can't be grasped in detail completely.
Simultaneously when analyzing various lightning protection measures effect, these methods can only one-parameter be calculated, relatively and adjust at present, can't analyze the lightning protection properties variation tendency that multi-parameter comprehensive changes the back circuit, are difficult to find optimum point when taking lightning protection measures especially.Finally cause transmission line of electricity in the design and running stage, lightning protection measures uses and lacks science and systemic.
Summary of the invention
The objective of the invention is in order to overcome the deficiency that existing transmission line of electricity lightning protection decision method can't carry out road emulation completely, multi-parameter integrated optimizing, a kind of all fronts road, multi-parameter integrated optimizing method of transmission line lightning stroke tripping operation emulation have been proposed, this method can be carried out the emulation of all fronts road, multi-parameter integrated optimizing to the lightning protection properties of transmission line of electricity, make transmission line of electricity in the design and running stage, the lightning protection measures of taking has more specific aim, science and systemic.
A kind of all fronts road, multi-parameter integrated optimizing method of judging the transmission line of electricity lightning protection properties that the present invention proposes is characterized in that this method specifically may further comprise the steps:
1) the selected transmission line of electricity (can be that the whole piece circuit also can be certain section circuit of relatively being concerned about) that needs judgement;
2) be that unit pursues choosing of bar parameter with the shaft tower to this circuit, the parameter of choosing comprises adjustable parameter and non-adjustable parameter (adjustable parameter can comprise: insulation configuration, stake resistance, lightning protection measures, lightning conducter shielding angle, non-adjustable parameter can comprise: lightening activity parameter, shaft tower parameter, line parameter circuit value, sea level elevation, topography and geomorphology etc.);
3) according to selected parameter, utilize Monte Carlo method (handling a kind of mathematical method of stochastic problem), calculate the tripping rate with lightning strike of shaft tower by bar;
4) based on each shaft tower tripping rate with lightning strike, with shaft tower both sides horizontal span sum 1/2nd as weighting coefficient, the tripping rate with lightning strike of computational scheme, computing formula as the formula (1);
R is the tripping rate with lightning strike (1/100km.yr) of circuit in the formula
N is the overhead line structures number
R
iIt is the tripping rate with lightning strike (1/100km.yr) of i base shaft tower
L
iBe half of i base shaft tower both sides horizontal span sum
L is total line length (km);
5) lightning outage rate that obtains according to step 4), if this trip-out rate less than the setting value that satisfies the lightning protection properties requirement lightning protection properties meet the requirements, jump to step 10) and judge to finish, otherwise enter step 6);
6) according to each shaft tower tripping rate with lightning strike along the line, finding out the higher shaft tower of trip-out rate is the shaft tower of lightning protection properties weakness;
7) according to the adjustable extent of field condition to the selected adjustable parameter of this weakness shaft tower;
8) according to the adjustable extent of adjustable parameter or the lightning protection measures of employing, in range of adjustment, each parameter is selected access point respectively, calculate the tripping rate with lightning strike of selected point according to the method for step 4), and the tripping rate with lightning strike of selected point is simulated the multi-parameter integrated optimizing curve;
9) according to multi-parameter integrated optimizing curve (at the scene under Cao Zuo feasibility and the economy prerequisite), choose the proper optimization point, utilize the tripping rate with lightning strike after Monte Carlo method recomputates this base shaft tower improvement, jump to step 4);
10) judge end.
The characteristics of the technical solution adopted for the present invention to solve the technical problems are: not only cause the whole tripping rate with lightning strike computational accuracy of circuit to improve by the accurate modeling of bar, and can also help the concrete lightning protection properties of grasping each shaft tower along the line, have more specific aim when making next step take lightning protection measures.When pointedly weak shaft tower being taked lightning protection measures, at first determine the adjustable parameters in the actual motion, calculate when each adjustable parameter changes within the specific limits the trend curve that the shaft tower lightning protection properties changes then.So just can help the user according to actual conditions, determine the optimum point when multi-parameter comprehensive is adjusted.
The invention has the beneficial effects as follows, by the transmission line of electricity lightning protection properties is carried out all fronts road, multi-parameter integrated optimizing, the transmission line of electricity design and running of scientific guidance power department has improved the power supply reliability of transmission line of electricity effectively, has produced tangible economic benefit and social benefit.
Description of drawings
Fig. 1 is existing transmission line of electricity lightning protection properties decision method process flow diagram.
Fig. 2 is all fronts road, the multi-parameter integrated optimizing method overview flow chart of judgement transmission line of electricity lightning protection properties of the present invention.
Fig. 3 is the embodiment particular flow sheet of the Monte Carlo method that adopts of the present invention.
Fig. 4 is the electric geometric model schematic diagram that the present invention introduces.
Fig. 5 is that the electric geometric model that the present invention introduces calculates the risk of shielding failure synoptic diagram.
Fig. 6 is that volt-second characteristic intersects method judgement insulation flashover synoptic diagram.
Fig. 7 is the trip-out rate curve of certain bar transmission line of electricity shaft tower along the line of the inventive method calculating.
Fig. 8 is the lightning protection properties multi-parameters optimization curve of certain basic shaft tower of the inventive method calculating.
Embodiment
The present invention is further described below in conjunction with drawings and Examples.
All fronts road, the multi-parameter integrated optimizing method of the judgement transmission line of electricity lightning protection properties that the present invention proposes as shown in Figure 2, specifically may further comprise the steps:
1) selected circuit or the part of path that needs to judge (can be imported the general status about this circuit in this step.Mainly comprise line name, line length, overhead line structures number etc.Line name can be used as and calculate the filename that respective lines parameter and result of calculation file are preserved in the back, and line length and overhead line structures number are the call parameters in the computation process);
2) circuit is pursued choosing of bar parameter, concrete parameter comprises: insulation configuration, stake resistance, lightning protection measures (install the coupling ground wire of certain altitude additional; Install leakage conductor additional; Install negative angle protection pin etc. additional), lightning conducter shielding angle, lightening activity parameter, shaft tower parameter (shaft tower model; Exhale the journey height), line parameter circuit value (electric pressure: exchange 110kV, 220kV, 500kV, 750kV, 1000kV, direct current ± 500kV, ± 800kV; Shaft tower both sides horizontal span; ), sea level elevation, topography and geomorphology etc.Wherein insulation configuration, stake resistance, lightning protection measures, lightning conducter shielding angle are adjustable parameter, and lightening activity parameter, shaft tower parameter, line parameter circuit value, sea level elevation, topography and geomorphology are non-adjustable parameter;
3), adopt Monte Carlo method to calculate the tripping rate with lightning strike of each basic shaft tower according to selected parameter; Because the transmission line lightning stroke tripping operation itself is to have stronger property of probability, handles it with probabilistic method and more can embody its essential characteristic.Monte Carlo method is as a kind of stochastic problem computing method of classics, after its selected target function, can regulate and change influence factor as required, utilize the method for mathematics to produce the sample of random variable sequence of various different distributions, and provide the probability statistics model of simulate given problem thus, give the progressive statistic estimated value of the numerical solution that goes wrong then.The outstanding feature of Monte Carlo method is that Consideration is comprehensive, and the controlling factor choice easily, and computation model is simple and clear.
In the present invention, with the function of transmission line lightning stroke trip-out rate N as some stochastic variables.
N＝f(I，U，B…) (2)
Amplitude of lightning current I wherein, thunderbolt moment lead power-frequency voltage value U, the position B during the thunderbolt circuit etc. is a stochastic variable, and has different probability distribution f (I), f (U) and f (B).
When utilizing Monte Carlo method to calculate tripping rate with lightning strike, the situation that needs analog line repeatedly to be struck by lightning, promptly adopt the method for random sampling from the probability distribution of the stochastic variable of each lightning stroke trip influence factor f (I), f (U) and f (B), to sample, corresponding random number I when obtaining the k time thunderbolt
k, U
kAnd B
k, cause insulation flashover when utilizing corresponding criterion to judge these this time thunderbolts then.
y
k＝f(I
k，U
k，B
k)，k＝1，2，3…N (3)
If cause flashover, then y
k=1, otherwise y
k=0 through N sampling, just can obtain the probability estimate value:
The embodiment flow process of the Monte Carlo method that the present invention adopts specifically may further comprise the steps as shown in Figure 3:
The 1st step: read in the calculating desired parameters, main parameter comprises: insulation configuration, stake resistance, lightning protection measures (install the coupling ground wire of certain altitude additional; Install leakage conductor additional; Install negative angle protection pin etc. additional), lightning conducter shielding angle, lightening activity parameter, shaft tower parameter (shaft tower model; Exhale the journey height), line parameter circuit value (electric pressure: exchange 110kV, 220kV, 500kV, 750kV, 1000kV, direct current ± 500kV, ± 800kV; Shaft tower both sides horizontal span; ), sea level elevation, topography and geomorphology etc.
The 2nd step: simulation thunderbolt number of times is set
At first utilize formula (5) to calculate this model shaft tower under the thunderbolt density of appointment according to the shaft tower parameter of first step input, thunderbolt number of times on the 100 kilometers circuits in 1 year inherence (because transmission line of lightning strike is a stochastic process, just can make numerical convergence satisfy accuracy requirement thereby therefore must calculate enough behind the thunderbolt number of times.Consider the speed of calculating and the accuracy requirement of engineering calculation, the thunderbolt number of times of each simulation calculation can be equivalent to the thunderbolt number of times that suffers of 100 kilometers circuits of 10000-30000).
N
s＝N
g(4h
s+b)/10 (5)
H wherein
sBe in lightning conducter average height (rice) that b is two lightning conducter spacings (rice).
h
tBe the shaft tower height, f
sIt is the lightning conducter sag.
The 3rd step random sampling lightning current waveform
According to amplitude of lightning current probability distribution formula (7), wave head, wave rear be value 2.6/50 μ s respectively.
i＝-88logr
1 (7)
R wherein
1Be [0,1] equally distributed random number.
Power-frequency voltage instantaneous value on the 4th step random sampling lead
It is generally acknowledged that most of thunders are negative polarity, then be struck by lightning moment stack power-frequency voltage for more greatly on the occasion of the time situation more serious.Calculate and produce [0, a 1] equally distributed random number r
2, because the power frequency phase angle is equally distributed, the random number that produces can be transformed to the random number of phase angle thus on [0,2 π]:
＝2πr
2 (8)
After obtaining the value, just can calculate the three-phase main-frequency instantaneous voltage:
U
a＝U
msin() (9)
U wherein
mBe the phase voltage peak value;
The 5th step is according to electric geometric model decision thunderbolt position
In the online pass thunder width, thunderbolt circuit position mainly is divided into two kinds: counterattack and shielding.Can be divided into thunderbolt shaft tower and thunderbolt mid span again for counterattack.For thunderbolt span central authorities, not it is generally acknowledged to cause insulator arc-over.For the calculating of shielding, counterattack probability, the present invention has introduced electric geometric model.
Electric geometric model is meant that the physical dimension with the flash-over characteristic of thunder and lightning and circuit connects and a kind of geometrical calculation model of setting up.Its ultimate principle is based upon on following notion and the hypothesis basis:
(1) the thunder and lightning guide arrives and to be hit that object is critical to be hit apart from being uncertain before, arrives hitting apart from interior promptly to its discharge of which object;
(2) hit apart from being the function of lightning current, size is relevant with amplitude of lightning current.According to theoretical research and test, hitting distance and amplitude of lightning current has following relation:
r
s＝kI
p (12)
R in the formula
sExpression is hit apart from (m), and I represents amplitude of lightning current (kA).
K, p are two constants.The formula that IEEE recommended in 1997 is:
r
c＝10I
0.65 (13)
H wherein
dFor lead on average to ground level.r
cWith r
gBe respectively that lead (comprising lightning conducter) hits distance and distance is hit on ground.
(3) do not consider the to be struck by lightning form effect of object and the influence of other factors.
Incident angle was obeyed certain distribution function when (4) guide was near ground, from the statistics angle, and the thunder and lightning distribution density maximum of its vertical drop, and the thunder and lightning density that level is attacked drops to zero.In general calculating, can think that the equal vertical ground of all thunders falls.
With ground elevation θ
g=0 is example.For a certain amplitude of lightning current I
i, can calculate corresponding hitting apart from r according to formula (13)
CiAnd r
Gi, as shown in Figure 4.
In Fig. 4, with r
CiBeing radius, is that circular arc B is made in the center of circle with lightning conducter s and lead c respectively
iC
iAnd C
iD
i, two circular arcs meet at C
iAt r overhead
GiThe place makes a horizontal line D
iE
i, meet at D with circular arc
iBy camber line B
iC
i, C
iD
iWith straight line D
iE
iForm curve road along the line direction move, form locating surface.When amplitude of lightning current is I
iPilot head portion be positioned at locating surface when above, its discharge development is not subjected to line influence.Have only when pilot head portion arrives the location curved surface, ability is influenced hits the thing development to quilt.If be introduced into B
iC
iCambered surface, then lightning conducter is hit, and lead is protected, so B
iC
iBe the protection arc.As be introduced into C
iD
iCambered surface is then hit lead, i.e. the shielding of lightning conducter was lost efficacy shielding, C have been taken place
iD
iFor exposing arc.If guide's head enters D
iE
iThe earth is then hit on the plane, so D
iE
iFor the earth is caught the thunder face.
Suppose that the thunder and lightning guide evenly vertically develops from the high-altitude earthward.As can be seen from Figure 5, amplitude is I
iHit in protection arc B
iC
iWith exposure arc C
iD
iProbability, respectively with this two arc projection F of surface level again
iG
iAnd G
iD
iProportional.So G
iD
iAnd F
iD
iRatio be exactly lightning current I
iRisk of shielding failure.As seen from the figure, risk of shielding failure is relevant with amplitude of lightning current.The lightning current that amplitude is low has bigger risk of shielding failure.Under the ideal conditions, after amplitude reached maximum shielding electric current, risk of shielding failure was 0.But, under the meteorological adverse condition, do not get rid of the possibility that big electric current shielding takes place in the mountain area.
With the once thunderbolt in the computation process is example.Calculate risk of shielding failure P according to electric geometric model
a, computing machine produces [0, a 1] equally distributed random number r subsequently
3If, r
3＜p
a, think that then thunderbolt this time is shielding; Otherwise think and be counterattack.Further again, if strike back, then think thunderbolt shaft tower (promptly hitting bar rate g) obedience probability distribution as shown in table 1 to produce [0, a 1] equally distributed random number r again
4If, r
4＜g then thinks the thunderbolt shaft tower takes place, otherwise is thunderbolt span central authorities.Behind the thunderbolt position of having determined thunderbolt this time, the mid span that is struck by lightning in this way just can be thought lightning stroke flashover not take place.
Shaft tower hits the bar rate under the various topographic conditions of table 1
| The lightning conducter radical | 1 | 2 |
| The |
1/4 | 1/6 |
| The |
1/3 | 1/4 |
The 6th step was calculated Lightning Over-voltage
The Lightning Over-voltage is here calculated the two kinds of situations that comprise, i.e. shielding superpotential and counterattack superpotential are specially shielding and still strike back by the decision of the 5th step.
When the shielding situation took place, the superpotential that bear at the insulator two ends was calculated suc as formula (15)
Wherein I is an amplitude of lightning current, and Zc is the lead wave impedance, gets 400 Ω in the engineering calculation usually.When counterattack took place, if when coupling ground wire not being installed, the superpotential that bear at the insulator two ends was calculated as (16)
(16) meaning of parameters in the formula is as shown in table 2.
Each meaning of parameters of table 2 (13) formula
| I is the random sampling (kV) of amplitude of lightning current | h aFor lead to ground level (m) |
| A is lightning current wave head steepness (kA/ μ S) | h tFor lightning conducter to ground level (m) |
| U is the random sampling (kV) of power-frequency voltage instantaneous value | h gBe lightning conducter average height (m) |
| K is the coupling coefficient of lead and lightning conducter | h cBe lead average height (m) |
| β is the shaft tower diverting coefficient | L tBe shaft tower inductance (μ H) |
| R iBe shaft tower impulse earthed resistance (Ω) |
If lightning protection measures is that coupling ground wire has been installed, the superpotential amplitude that bear at the insulator two ends is calculated according to (17) formula.
L
Gt a, L
Gt bBe the shaft tower equivalent inductance (μ H) of cat head to coupling ground wire hitch point section, coupling ground wire hitch point to the column foot section.
L
B1, L
B2Inductance value in parallel (μ H) for shaft tower both sides adjacent shelves lightning conducter and coupling ground wire.
R
ChBe shaft tower impulse earthed resistance (Ω).
C
1, C
2Be the coupling coefficient between lead and ground wire, coupling ground wire.
β is the shaft tower diverting coefficient
U is the random sampling (kV) of power-frequency voltage instantaneous value
A is lightning current wave head steepness (μ S)
The 7th step flashover is judged
For the counterattack flashover, available formula (18) approximate treatment voltage-time curve utilizes volt-second characteristic to intersect the method (see figure 6) then and judges.
V wherein
FOBe flashover voltage value (kV) that t is the time (μ s) before the flashover, W is the shortest clearance or insulator dry arcing distance (m).
The insulator arc-over family curve of curve 1 for calculating by (18) among Fig. 6, curve 2 and curve 3 are respectively two Lightning Over-voltage curves.Curve 2 and curve 1 intersect constantly at t1, can think that this moment, flashover took place in insulation.Curve 3 and curve 1 are non-intersect, and flashover does not take place in insulation.
Judge for the shielding flashover, if the shielding superpotential that (15) formula is calculated is U
U
CFO＝533l+132 (19)
Wherein l is insulator dry arcing distance or the shortest clearance (m)
If U＞U
CFO, i.e. shielding causes insulator arc-over, otherwise shielding does not cause insulator arc-over.
The 8th step cycle calculations
The 3rd step to the 7th step is the emulation of a lightning stroke process, and the thunderbolt emulation of every basic shaft tower should be carried out according to the number of times that second step calculated.Therefore need be to the 3rd step and the 7th continuous double counting of step, up to satisfying the calculation times requirement.
The 9th step was calculated the tripping operation ratio
The probability of sustained arc that multiply by formula (20) according to the lightning stroke flashover accident number of times of adding up is previously compared with total thunderbolt number of times, can obtain the ratio that the lightning stroke trip accident accounts for total thunderbolt number of times.
η＝(4.5E
0.75-14)×100
-2 (20)
E is the average working voltage gradient (kV/m) of insulator chain in the formula.
To effectively earthed system:
To neutral point insulation, arc suppression coil earthing system:
More than U in two formulas
nBe system's rated voltage (kV); l
1Be insulator length (m); l
2Be the wire spacing (m) of wooden cross arm circuit, to iron cross arm and reinforced concrete cross-arm circuit, l
2=0.
If E≤6kV/m, probability of sustained arc is very little, can be similar to think η=0.
The 10th step tripping rate with lightning strike
At present, generally adopt 1/100km.yr for tripping rate with lightning strike unit in the engineering calculation.Therefore the result of calculation in the 9th step need be carried out unit conversion.Suppose that it is P that the 9th lightning stroke trip accident that obtains of step accounts for total thunderbolt number of times ratio, and the thunderbolt number of times that per 100 kilometers of computational schemes suffered in a year to calculate gained by formula (4) be N
L, then the tripping rate with lightning strike of this circuit is N
LP/100km.yr.
4) weighted comprehensive computational scheme tripping rate with lightning strike.This step is calculated as the basis with previous step shaft tower tripping rate with lightning strike fully.Can the tripping rate with lightning strike of shaft tower be sorted, trip-out rate is higher to be some shaft tower of lightning protection properties weakness.The tripping rate with lightning strike of circuit calculates by formula (23)
R is the tripping rate with lightning strike (1/100km.yr) of circuit in the formula
N is the overhead line structures number
R
iIt is the tripping rate with lightning strike (1/100km.yr) of i base shaft tower
L
iBe half of i base shaft tower both sides horizontal span sum
L is total line length (km)
5) judge whether the line thunder protection performance needs to improve.If the trip-out rate value that formula (23) is calculated is less than satisfying the setting value (this setting value generally is made as electrical network competent authorities to this lightning outage rate examination value) that lightning protection properties requires, then this line thunder protection performance is qualified, then skip to step 10), judge and finish, otherwise need to improve the line thunder protection performance, change step 6) over to;
6) the relatively weaker shaft tower of lightning protection in the selected circuit.Weak shaft tower is generally chosen the higher shaft tower (consider the randomness of lightening activity, generally do not recommend to select weak shaft tower according to lightning strike accident situation in the actual motion) of trip-out rate that ordering obtains in the step 4);
7) adjustable extent of adjustable parameter in the selected shaft tower.Main adjustable parameter comprises stake resistance, insulation configuration, lightning protection measures (comprise the coupling ground wire that installs certain altitude additional, install leakage conductor additional, install negative angle protection pin etc. additional), also is adjustable parameter at design phase lightning conducter shielding angle.Whether these parameters are adjustable should judge according to the circuit actual conditions;
8) calculate the multi-parameter integrated optimizing curve.After adjustable parameter is determined, adjustable extent in routine, each parameter is selected access point respectively, calculate the tripping rate with lightning strike of selected point according to the method for step 4), tripping rate with lightning strike with these selected points simulates multi-parameter integrated optimizing curve (in order to guarantee computational accuracy and computing velocity simultaneously, 5 left-right dots of choosing of general every curve are used for calculating) then.
9) according to 8) the multi-parameters optimization curve of result of calculation, at the scene under Cao Zuo feasibility and the economy prerequisite, choose the proper optimization point, to the parameter substitution original parameter of this base shaft tower with the optimization point, utilize the tripping rate with lightning strike after Monte Carlo method recomputates this base shaft tower improvement, and jump to step 3);
10) judge end.
The effect of the application example of employing the inventive method is as follows:
Fig. 7 is the shaft tower tripping rate with lightning strike distribution schematic diagram along the line that adopts the inventive method to obtain to certain bar transmission line of electricity.Calculating simultaneously and calculating this lightning outage rate calculated value according to (20) formula is 1.41/100km.yr, and statistical value is 1.37/100km.yr, and the error of calculated value and actual value is 6%.As can be seen from the figure, this transmission line of electricity has 307 basic shaft towers.Wherein use the 101# of red mark, 102#, 105#, 106#, the 114# tripping rate with lightning strike is higher, belongs to the shaft tower that needs emphasis to strengthen lightning protection properties.
Fig. 8 is the highest 114# shaft tower of this a line tripping rate multi-parameter integrated optimizing curve map.As can be seen from the figure, reduce shielding angle and be to improve the effective measures of this base shaft tower lightning protection properties, should select the shaft tower of little shielding angle in this basic shaft tower position for use, simultaneously suitable reinforced insulation.Concrete optimum point should be based on this curve, and Electric Design and operations staff suitably choose from the feasibility of execute-in-place and the economy of transformation.
Claims (1)
1, a kind of all fronts road, multi-parameter integrated optimizing method of judging the transmission line of electricity lightning protection properties is characterized in that this method specifically may further comprise the steps:
1) the selected transmission line of electricity that needs judgement;
2) be that unit pursues choosing of bar parameter with the shaft tower to this circuit, the parameter of choosing comprises adjustable parameter and non-adjustable parameter;
3) according to selected parameter, utilize Monte Carlo method, calculate the tripping rate with lightning strike of shaft tower by bar;
4) based on each shaft tower tripping rate with lightning strike, with shaft tower both sides horizontal span sum 1/2nd as weighting coefficient, the tripping rate with lightning strike of computational scheme;
R is the tripping rate with lightning strike of circuit in the formula, and unit is 1/100km.yr, and N is the overhead line structures number, R
iThe tripping rate with lightning strike unit that is i base shaft tower is 1/100km.yr, L
iBe half of i base shaft tower both sides horizontal span sum, L is a total line length, and unit is km;
5) lightning outage rate that obtains according to step 4), if this trip-out rate less than the setting value that satisfies the lightning protection properties requirement lightning protection properties meet the requirements, jump to step 10) and judge to finish, otherwise enter step 6);
6) according to each shaft tower tripping rate with lightning strike along the line, finding out the higher shaft tower of trip-out rate is the shaft tower of lightning protection properties weakness;
7) according to the adjustable extent of field condition to the selected adjustable parameter of this weakness shaft tower;
8) according to the adjustable extent of adjustable parameter or the lightning protection measures of employing, in range of adjustment, choose the different numerical points of each parameter respectively, calculate the tripping rate with lightning strike of selected point according to the method for step 4), and the tripping rate with lightning strike of selected point is simulated the multi-parameter integrated optimizing curve;
9) under the feasibility and economy prerequisite of operating at the scene according to the multi-parameter integrated optimizing curve, choose the proper optimization point, utilize Monte Carlo method to recomputate tripping rate with lightning strike after this base shaft tower improves, jump to step 4);
10) judge end.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2005101166946A CN100361123C (en) | 2005-10-28 | 2005-10-28 | Whole line and multi-parameter integrated optimizing method for determining lightningproof performance of power transmission line |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2005101166946A CN100361123C (en) | 2005-10-28 | 2005-10-28 | Whole line and multi-parameter integrated optimizing method for determining lightningproof performance of power transmission line |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1770140A CN1770140A (en) | 2006-05-10 |
| CN100361123C true CN100361123C (en) | 2008-01-09 |
Family
ID=36751447
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2005101166946A Expired - Fee Related CN100361123C (en) | 2005-10-28 | 2005-10-28 | Whole line and multi-parameter integrated optimizing method for determining lightningproof performance of power transmission line |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN100361123C (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012145876A1 (en) * | 2011-04-29 | 2012-11-01 | 山西省电力公司 | Method for determining spatial location of conducting wire and aerial earth wire of power transmission line |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100498357C (en) * | 2007-09-14 | 2009-06-10 | 清华大学 | High voltage electricity transmission system thunder-proof property parameter test method |
| CN101900779B (en) * | 2010-07-28 | 2012-10-17 | 西南交通大学 | Method for measuring lightning shielding failure trip-out rate of high-voltage power transmission line by combining topography |
| CN102072992B (en) * | 2010-10-26 | 2013-01-23 | 国网电力科学研究院武汉南瑞有限责任公司 | Detailed terrain data-based transmission line failure-shielding and lightning-protection performance evaluation method |
| CN102435921B (en) * | 2011-09-26 | 2015-10-28 | 山西省电力公司忻州供电分公司 | The decision method of same tower double back transmission line insulation and resistance to lightning impulse performance |
| CN103675510B (en) * | 2012-09-20 | 2016-12-21 | 西南交通大学 | The method of overhead section of electric railway tripping rate with lightning strike of single line under evaluation and test AT mode |
| CN103236666B (en) * | 2013-04-02 | 2015-12-09 | 国家电网公司 | Transmission line located lightening arresting method |
| CN103488815A (en) * | 2013-08-26 | 2014-01-01 | 国家电网公司 | Transmission line lightning shielding failure risk evaluation method |
| CN103915810B (en) * | 2014-04-17 | 2017-02-15 | 国家电网公司 | Lightning protection method for 500 kV circuit JG type pole and tower |
| CN109586226A (en) * | 2018-12-04 | 2019-04-05 | 金华八达集团有限公司科技信息分公司 | The method and device of steel tower tower head overvoltage when a kind of inhibition transmission line lightning stroke |
| CN109521343B (en) * | 2018-12-29 | 2020-11-10 | 广东电网有限责任公司 | Method for evaluating lightning-induced tower protection range |
| CN110879333B (en) * | 2019-12-05 | 2022-02-01 | 国家电网有限公司 | High-voltage direct-current transmission line lightning stroke restart rate calculation and evaluation method |
| CN112487670A (en) * | 2020-12-25 | 2021-03-12 | 国网山东省电力公司电力科学研究院 | Generation method of lightning overvoltage simulation model of large-span high tower |
| CN113777441B (en) * | 2021-08-20 | 2024-04-09 | 云南电网有限责任公司楚雄供电局 | Lightning strike same-jump evaluation method and platform considering height of coupling ground wire |
| CN113572143B (en) * | 2021-09-26 | 2022-01-25 | 广东电网有限责任公司 | Lightning protection method and device for power transmission line under multiple lightning strikes |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH1080058A (en) * | 1996-09-02 | 1998-03-24 | Morinaga Denshi Kk | Surge absorber unit for electronic equipment |
| CN2336493Y (en) * | 1998-02-10 | 1999-09-01 | 江西省电力通信公司 | Lightning protection device for transmission lines |
-
2005
- 2005-10-28 CN CNB2005101166946A patent/CN100361123C/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH1080058A (en) * | 1996-09-02 | 1998-03-24 | Morinaga Denshi Kk | Surge absorber unit for electronic equipment |
| CN2336493Y (en) * | 1998-02-10 | 1999-09-01 | 江西省电力通信公司 | Lightning protection device for transmission lines |
Non-Patent Citations (3)
| Title |
|---|
| Line surge arrester selection to improve lightningperformance of transmission lines. J.A.Tarchini and W.Gimenez.IEEE. 2003 * |
| Transmission Line Lightning Performance Based Design. A.P.Sakis Meliopoulos, Ron Cooper.IEEE. 2000 * |
| 架空输电线路雷击跳闸率影响因素研究现状. 孙鹞鸿等.高电压技术，第30卷第12期. 2004 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012145876A1 (en) * | 2011-04-29 | 2012-11-01 | 山西省电力公司 | Method for determining spatial location of conducting wire and aerial earth wire of power transmission line |
| US20140052393A1 (en) * | 2011-04-29 | 2014-02-20 | Shanxi Electric Power Research Institute | Method for determining spatial location of conducting wire and aerial earth wire of power transmission line |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1770140A (en) | 2006-05-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN100361123C (en) | Whole line and multi-parameter integrated optimizing method for determining lightningproof performance of power transmission line | |
| CN102435921B (en) | The decision method of same tower double back transmission line insulation and resistance to lightning impulse performance | |
| CN102411105B (en) | Evaluation method of lightning disturbance risk of power transmission line section | |
| CN101315400B (en) | Electric power line lightning protection performance estimation method based on thunder and lightning parametric statistics | |
| CN103488815A (en) | Transmission line lightning shielding failure risk evaluation method | |
| CN103530453B (en) | Transmission line of lightning strike causes the analysis method of extra-high voltage direct-current system generation locking | |
| CN105242133B (en) | A kind of improvement distribution line lightning trip-out rate calculation method | |
| CN109753703A (en) | A kind of distribution line lightning protection grade appraisal procedure | |
| CN101261616A (en) | Electric network thunderstorm damage distribution determination method | |
| CN107992962A (en) | A kind of Lightning stroke Protection Measures for Over-Head Lines optimum choice method based on entropy assessment | |
| CN104836223A (en) | Power grid parameter error and bad data coordinated identification and estimation method | |
| CN110309527A (en) | A kind of overhead distribution damage to crops caused by thunder methods of risk assessment based on electric geometry method | |
| CN114236288B (en) | Fault positioning method based on power transmission line | |
| CN109917175A (en) | It is a kind of for high anti-back-out when overvoltage method for quick predicting | |
| CN105406471A (en) | Bad data identification and estimation method for power grid | |
| CN110954767A (en) | Accurate analysis method for power line loss | |
| CN111695775A (en) | Power transmission line lightning stroke risk assessment method based on optimal combined weighting of preference ratio-entropy method | |
| CN104392397A (en) | Entropy weight method based short-circuit current suppression scheme evaluation method | |
| CN106443253A (en) | Power transmission line parameter identification method based on PMU (phasor measurement unit) data | |
| CN105912509A (en) | Statistical method of cumulative probability distribution of multi-section lightning current amplitude | |
| CN110504707A (en) | Multi-infeed DC network voltage temporarily stabilizes horizontal appraisal procedure and device | |
| CN106872861A (en) | Super UHV Transmission Engineering circuit attacks the lightning current inversion method and system of wire directly | |
| CN110298053A (en) | A kind of calculation method of the transmission line lightning stroke trip-out rate based on the protection of inductive type lightning rod | |
| CN105929264B (en) | A kind of 750kV transmission line thunderbolts performance estimating method | |
| CN102968668A (en) | Urban power distribution network evaluation system and method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20080109 Termination date: 20171028 |