Specific embodiment
The technical scheme in the embodiment of the invention is clearly and completely described with reference to the accompanying drawings and examples, shows
So, described embodiments are only a part of the embodiments of the present invention, rather than whole embodiments.
Route symmetrical for the triphase parameter that uniformly replaces, i.e. balanced transmission line road when theoretical calculation, unit length
Positive and negative zero sequence impedance can be calculated (unit is Ω/km) by following formula respectively:
In above formula: R indicates D.C. resistance, and 0.05 for characterizing ground resistance (ground resistance=ZL=π2f×10-4(Ω/
Km), when taking power frequency 50Hz, value about 0.05), r ' expression conducting wire equivalence radius, DgIndicate that ground returns to road equivalence depth, Dave
Indicate the geometric mean distance between power transmission line (when considering single loop lineConsider the geometric mean distance between double loop
When,dmnIndicate the space length between corresponding conductor.
In practice, even if to same back transmission line, due in erection process along the line can across different topography and geomorphology and
Building etc., soil resistivity along the line are not necessarily worth, but change with the variation of geological conditions, therefore along the line
Soil resistivity is variation;In addition, the transmission tower for the construction of line is also not changeless, but with leap side
It formula and is spanned object and accordingly changes, therefore the position between conductor can change.It is influenced by geological conditions, soil resistivity
ρ leads to be changed within the scope of 100~4000 Ω/m, thus lgDgVariation range be 2.9701~3.7711.For different structure
Shaft tower, DaveValue there is also difference, especially compact shaft tower and conventional type shaft tower, DaveDifference it is fairly obvious.Consider
After the influence of these two aspects factor, zero-sequence mutual inductance impedance may be fluctuated in a big way.
Since the co-ordination that the measurement of zero-sequence mutual inductance is related to is sufficiently complex, and it is limited to the erection situation of route, because
Empirical method is generallyd use in this engineering to estimate the zero-sequence mutual inductance impedance of route, that is, uses a certain fixed coefficient and solid wire
Zero sequence impedance is multiplied to replace the zero sequence mutual impedance between same tower route.But from the point of view of accurately calculating, calculated result and true value it
Between error there is very big randomness, and it is undesirable.
Based on this, the present invention proposes a kind of method for calculating zero-sequence mutual inductance impedance on multiple-loop line transmission line of electricity, referring to Fig. 1
And Fig. 2 comprising have:
Step 1, it obtains and is related to the basic data of each loop line, each base shaft tower, transmission pressure, ground wire;
Step 2, it soil will be equivalent to enable with the equally distributed soil resistivity of sectional with tower section electric resistance of soil along the line
Rate is ρt, non-same tower section soil resistivity be ρf;
Step 3, all fronts soil is equivalent to soil resistivity having the same, enabling the soil resistivity is the uniform of route
Soil resistivity ρj, use ρjConstruct ρt、ρfIterative initial value, if sharing v loop line with tower section, wherein v be in 3 to 6 between, then
With tower section soil resistivity ρ on r loop linetWith non-same tower section soil resistivity ρfIterative initial value can then be taken as respectively: (ρj1
+…+ρjv)/v and ρjr;
Step 4, the transmission line of electricity set up for local multiple-loop line writes each loop line about ρ by columnt、ρfLinear algebra
Equation group, and iterative method and above-mentioned iterative initial value are used, solution is ρ with the soil resistivity of tower sectiontWith the soil of non-same tower section
Earth resistivity is ρf;
Step 5, according to each parameter value in basic data, by the soil resistivity ρ of same tower section and non-same tower sectiont、ρfRespectively
It substitutes into the shaft tower of corresponding model, finds out the impedance matrix of each model shaft tower unit length of each loop line, and utilize symmetrical component method
It can be from the zero-sequence mutual inductance impedance for finding out each model shaft tower unit length in each loop line in impedance matrix;
Step 6, zero-sequence mutual inductance value of each base shaft tower of every loop line in its span is calculated, by these zero-sequence mutual inductance values
It is cumulative, it can be obtained the zero-sequence mutual inductance value of whole route.
The basic data, the specific embodiment as step 1, comprising: each loop line is used every in erection process
The actual measurement zero sequence resistance of the phase sequence arrangement mode of conducting wire, each loop line on the structural parameters and its span of one base shaft tower, each base shaft tower
Anti-, the radical of transmission pressure parameter, ground wire and ground wire conductor parameter, the overall length of route and the line length with tower section, it is described
Transmission pressure parameter includes division number, line footpath and D.C. resistance.
Tower structure and soil resistivity are the principal elements for influencing zero-sequence mutual inductance impedance, along the line the structure ginseng of each base shaft tower
Number can be obtained from circuit department, but soil resistivity along the line can not survey one by one.Therefore in practice, it is mutual to influence zero sequence
The principal element that sense impedance accurately calculates is soil resistivity.The step 2 and step 3 are at the equalization to soil resistivity
Soil resistivity along the line is carried out equalization processing, is the committed step for calculating the impedance of all fronts zero-sequence mutual inductance, i.e., by zero sequence by reason
Mutual inductance is thought to be distributed along shaft tower piecewise uniform.
Specifically, equalization can be stated after considering the influence of tower structure are as follows: in the same tower section of transmission line of electricity,
Think that zero-sequence mutual inductance is equally distributed;In the non-same tower section of transmission line of electricity, it is believed that zero-sequence mutual inductance is also equally distributed;Think
Had differences between tower section and the zero-sequence mutual inductance of non-same tower section, the difference by the soil resistivity of same tower section and non-same tower section Lai
It embodies.When calculating by soil resistivity segment processing, the soil resistivity with tower section is taken as ρt, the electric resistance of soil of non-same tower section
Rate is taken as ρf.Solving ρt、ρfAfterwards, it can arrange and write the zero-sequence mutual inductance impedance that Algebraic Equation set finds out route.
The zero sequence impedance measured value of each loop line is provided by basic data, can thus be arranged and be write each loop line about ρt、ρfLinear generation
Number equation, is solved with iterative method.And the convergence rate of iterative method and the selection of initial value are closely related, inappropriate initial value,
So that required calculation times, calculating time and calculation amount all greatly increase in solution procedure, practicability is low for meeting.According to zero sequence
The experience estimated value of mutual inductive impedance can calculate a corresponding experience soil resistivity ρ0, but ρ0And the true soil of logicalnot circuit
Resistivity, the ρ calculated in some cases0May not obviously it conform to outside 100~4000 Ω/m range in this, as initial value
It is suitable.Therefore it needs further to seek ρt、ρfRelatively reasonable iterative initial value.
Introduce the concept of uniform soil resistivity below to illustrate how to obtain ρt、ρfIterative initial value.
The influence of tower structure is taken into account, it is believed that completely soil resistivity having the same is calculated with the resistivity
Zero sequence impedance out is equal to the measured value of route zero sequence impedance, which is referred to as the uniform soil resistivity of route
ρj。
According to each parameter value in basic data, soil resistivity is uniformly taken as ρj, each model shaft tower of the loop line can be found out
The impedance matrix of unit length, and can be long from each model shaft tower unit in the loop line is found out in impedance matrix using symmetrical component method
The zero sequence impedance of degree.Then ρjMeet the constraint of following formula:
Zr0=Zrf1(ρjr)Lrf1+…+Zrfp(ρjr)Lrfp+Zrt1(ρjr)Lrt1+…+Zrtq(ρjr)Lrtq (1)
In above formula, Zr0For the zero sequence impedance measured value of r loop line;The shaft tower of the non-same tower section number s of r loop line corresponds to unit
The zero sequence impedance and span of length are expressed as Zrfs(ρjr) and Lrfs, Zrfs(ρjr) it is about the uniform electric resistance of soil of r loop line
Rate ρjrExpression formula;R loop line corresponds to the zero sequence impedance of unit length with the shaft tower of tower section number s and span is expressed as
Zrts(ρjr) and Lrts, Zrfs(ρjr) it is about the uniform soil resistivity ρ of r loop linejrExpression formula.
It is constrained by the equation of formula (1), therefore can will be used as ρ0ρjrInitial value is iterated calculating, finally finds out uniform soil
Earth electricalresistivityρjExact value.To be better described, illustrate how to solve uniform soil resistivity below by taking dichotomy as an example
ρj, calculation process is as shown in Figure 3.But the method that the present invention solves uniform soil resistivity is not limited thereto.
1) zero sequence impedance discreet value is calculated
If the experience estimated value of the zero-sequence mutual inductance impedance on r loop line road is ZrM0', and Tm is set as the longest bar of conveying distance
Tower number calculates the corresponding geometric mean distance D' of Tmave.According to ZrM0’、Dave' value, the calculating of the zero-sequence mutual inductance impedance of substitution
Formula can acquire experience soil resistivity ρ0。
By ρ0As ρrjInitial value substitute into formula (1), obtain the discreet value Z of route zero sequence impedance when the 0th corrected Calculationr0
(0), if the zero sequence impedance measured value of r loop line is Zr0, remember deviation ε0=Z0(0)-Z0.According to the actual conditions of route, the limits of error
It is taken as σ, when | ε0When | < σ, it is believed that ρ0The uniform soil resistivity for as needing to solve.|ε0When | > σ, then it is assumed that ρ0Value not
Appropriately, it need to be modified.
2) soil resistivity is corrected
The variation range of usual soil resistivity, might as well be by the lower limit 100 of soil resistivity between 100~4000 Ω/m
Ω/m is assigned to a as initial value, and 4000 Ω of the upper limit/m is assigned to b as initial value, then soil resistivity ρ0Inevitable section [a,
B] in.The soil resistivity for remembering n-th corrected Calculation is ρn, the zero sequence impedance value of the route A of n-th corrected Calculation is Z0(n),
The deviation of n-th corrected Calculation is εn.Then makeover process is as follows:
If ε0> σ, i.e. Z0(0)>Z0, due to Z0(0) it is proportional to zero sequence mutual impedance, illustrates that zero sequence mutual impedance is bigger than normal at this time, answers
Suitably reduce this value.And zero sequence mutual impedance is proportional to soil resistivity, therefore can be by the ρ that suitably reduces0Value carrys out adjusting circuit
Zero sequence impedance discreet value.At this point, by ρ0It is assigned to b, a is remained unchanged.
If ε0< σ, i.e. Z0(0)<Z0, due to Z0(0) it is proportional to zero sequence mutual impedance, illustrates that zero sequence mutual impedance is less than normal at this time, answers
Suitably increase this value.And zero sequence mutual impedance is proportional to soil resistivity, therefore can be by the ρ that suitably increases0Value carrys out adjusting circuit
Zero sequence impedance discreet value.At this point, by ρ0It is assigned to a, b is remained unchanged.
Enable ρ1=(a+b)/2, repeats the above steps and calculates Z0(1)、ε1If | ε1| < σ then stops calculating, by ρ1As uniform
Otherwise soil resistivity repeats the above process for calculating uniform zero sequence mutual impedance and calculates separately ρ2、ρ3、ρ4... until | εj|<
Stop calculating when σ, by ρ at this timejAs uniform soil resistivity.
3) ρ is obtainedt、ρfIterative initial value
If sharing v loop line (6 >=v >=3) with tower section, the uniform soil of each loop line is calculated separately out using the above method
Electricalresistivityρj1…ρjv, can be then taken as respectively with the iterative initial value of tower section and non-same tower section soil resistivity on r loop line: (ρj1
+…+ρjv)/v、ρjr。
Referring to Fig. 4, the step 4 calculates non-same tower section or that the specific method is as follows is described with tower section resistivity:
As one embodiment, if multiple-loop line line uses the whole erection mode with tower, with tower section and non-same tower section
Soil resistivity ρt、ρfIt is equal and all be uniform soil resistivity ρj, no longer need to solve.By ρjSubstituting into formula can calculate
The zero-sequence mutual inductance impedance of each loop line.
It is logical with the inlet wire of tower section and outlet for the transmission line of electricity that local multiple-loop line is set up as another embodiment
It is often formed as unit of double-circuit line on same pole, the inlet wire and outgoing line part of a small number of multiple-loop line lines can include single loop line.Specifically:
The inlet wire (or outlet) of transmission line with four-circuit on single tower is usually made of two groups of double-circuit line on same pole;With the inlet wire of six back transmission line of tower
(or outlet) is usually made of three groups of double-circuit line on same pole;The situation for returning erection with tower three is less, and inlet wire (or outlet) is usual at this time
It is made of one group of double-circuit line on same pole and a single loop line.
As another embodiment, for the double loop of parallel lines on same tower, due to whole parallel lines on same tower, Origin And Destination phase
Together, therefore the soil resistivity along this two loop line also must be the same, ρt=ρt2, ρf1=ρf2.It is set up jointly in multiple-loop line line
Part, geological conditions along the line is identical, therefore the soil resistivity on each loop line of this section of route is equal ρt1=ρt2=...=
ρtv。
By taking quadri-circuit lines on the same tower as an example, if I loop line and II loop line are one group of double-circuit line on same pole, whole process is joint used;III loop line
It is two single loop lines with IV loop line.According to the zero sequence impedance value that each loop line is surveyed, it can arrange and write following equation group:
ZI0=ZIf1(ρf1)LIf1+…+ZIfp(ρf1)LIfp+ZIt1(ρt)LIt1…+ZItq(ρt)LItq (2)
ZII0=ZIIf1(ρf2)LIIf1+…+ZIIfp(ρf2)LIIfp+ZIIt1(ρt)LIIt1…+ZIItq(ρt)LIitq (3)
ZIII0=ZIIIf1(ρf3)Lf1+…+ZIIIfp(ρf3)LIIIfp+ZIIIt1(ρt)LIIIt1…+ZIIItq(ρt)LIIItq (4)
ZIV0=ZIVf1(ρf4)LIvf1+…+ZIVfp(ρf4)LIvfp+ZIVt1(ρt)LIvt1…+ZIVtq(ρt)LIvtq (5)
In formula (2)~(5): ZI0~ZIV0, respectively indicate the zero sequence impedance of I~IV loop line;ρf1~ρf4, respectively indicate I~
The soil resistivity of the non-same tower section of IV loop line;ρt1~ρt4, I~IV loop line is respectively indicated with the soil resistivity of tower section;ρf1~
ρf4, respectively indicate the soil resistivity of the non-same tower section of I~IV loop line;It is long that the shaft tower of the non-same tower section number s of r loop line corresponds to unit
The zero sequence impedance and span of degree are expressed as Zrfs(ρfr) and Lrfs, Zrfs(ρfr) it is about the non-same tower section electric resistance of soil of the loop line
Rate ρfrExpression formula;R loop line corresponds to the zero sequence impedance of unit length with the shaft tower of tower section number s and span is expressed as
Zrfs(ρt) and Lrfs, Zrfs(ρt) it is about the loop line with tower section soil resistivity ρtExpression formula.
To be better described, below illustrate how to solve each loop line with tower section and non-same tower section for gaussian iteration method
Soil resistivity, calculation process are as shown in Figure 2.But each loop line of solution of the invention is the same as the soil resistivity of tower section and non-same tower section
Method be not limited thereto.
Since I loop line and II loop line double circuits on same tower are set up, there is ρf1=ρf2.It, can be by formula (2) by simply converting
~(5) are rewritten as following form:
ρt=a1ρt+a2ρf1+b1 (6)
ρf1=a3ρt+a4ρf1+b2 (7)
ρf3=a5ρf3+a6ρt+b3 (8)
ρf4=a7ρt+a8ρf4+b4 (9)
In formula (6)~(9), a1~a8 is constant, value depend on structural parameters with tower section and non-same tower section bar and
The corresponding span of each shaft tower.
Formula (6)~(9) are rewritten as to the form of ρ=A ρ+B, i.e., are as follows:
Gaussian iteration form corresponding to formula (9) are as follows:
ρ(k+1)=A ρ(k)+ B, k=0,1,2 ... (11)
The initial value of ρ is expressed as
Enable ρt0=(ρj1+ρj2+ρj3+ρtj4)/4, to tax initial value ρ(0)=[ρt0 ρj1 ρj2 ρj3]T,
Note d=max | ρ(k+1)-ρ(k)|, expansion is
The limits of error is taken as δ, constantly carries out Gaussian processes iteration using formula (11), calculates ρ(1)、ρ(2)..., until as d < δ
Stop iteration, each loop line can be obtained at this time with the accurate solution ρ of tower section and non-same tower section soil resistivityt、ρf1(=ρf2)、ρf3、
ρf4。
The step 5 is to seek the zero sequence mutual impedance between route according to each parameter value in basic data, specific to solve such as
Under, remember respectively: the zero-sequence mutual inductance impedance that the shaft tower of the non-same tower section number s of r loop line corresponds to unit length is Zrmfs(ρfr), Zrmfs
(ρfr) it is about the non-same tower section soil resistivity ρ of the loop linefrExpression formula;R loop line is corresponding single with the shaft tower of tower section number s
The zero-sequence mutual inductance impedance of bit length is Zrmfs(ρt), Zrmfs(ρt) it is about the loop line with tower section soil resistivity ρtExpression formula.
Then zero-sequence mutual inductance impedance relevant to r loop line are as follows:
Zrm0=Zrmf1(ρfr)Lrf1+…+Zrmfp(ρfr)Lrfp+Zrmt1(ρt)Lrt1…+ZImtq(ρt)Lrtq
The zero-sequence mutual inductance impedance value for calculating each loop line Yu other loop lines substitutes into the numerical value of the basic data of specific route
Above formula.
The present invention provides a kind of zero-sequence mutual inductance impedance computation methods on multiple-loop line transmission line of electricity, using present invention side
Method calculates the zero sequence mutual impedance of transmission line of electricity, according to the actual conditions founding mathematical models of route, calculated zero sequence
Mutual inductive impedance has comprehensively considered the practical erection situation and topography variation situation of transmission line of electricity, there is higher accuracy, is not necessarily to root
Line impedance is estimated according to empirical coefficient.The calculated zero-sequence mutual inductance impedance of this method and soil resistivity are applied to protect
It protects in the research work such as adaptive setting, short circuit calculation, Load flow calculation, Computer Simulation, the essence of analytical calculation can be greatly improved
True degree reduces the error as caused by experience estimation.And the calculating section of the method for the present invention can be complete in computer platform
At, actual mechanical process is simple, have it is very strong can practicality.
Above-described embodiment is only to illustrate the present invention and not limits the technical scheme described by the invention;Therefore, although
The present invention has been described in detail referring to above-mentioned each embodiment for this specification, still, the ordinary skill of this field
Personnel should be appreciated that and still can modify or equivalently replace the present invention;And all do not depart from spirit of the invention and
The technical solution and its improvement of range, are intended to be within the scope of the claims of the invention.