The content of the invention
It is an object of the invention to for above-mentioned problems of the prior art, it is proposed that a kind of multiple-loop line power transmission line
Zero-sequence mutual inductance impedance computation method on road, reduces the mistake caused by rough estimation is carried out to zero sequence mutual impedance parameter as far as possible
Difference, improves the stability of protective relaying device and the accuracy of Electrical power system analysis and computing.
To reach the purpose of foregoing invention, the present invention is achieved through the following technical solutions:
A kind of zero-sequence mutual inductance impedance computation method of the invention disclosed on multiple-loop line transmission line of electricity, including step are as follows:
Step 1, the basic data for being related to each loop line, base shaft tower, transmission pressure, ground wire is obtained;
Step 2, soil it will be equivalent to, with the equally distributed soil resistivity of sectional, make 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 identical soil resistivity, it is the homogeneous of circuit to make the soil resistivity
Soil resistivity ρj, use ρjConstruct ρt、ρfIterative initial value, if having v loop lines with tower section, wherein v is between 3 to 6, then
With tower section soil resistivity ρ on r loop linestWith 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, each loop line is write on ρ by rowt、ρfLinear algebra
Equation group, and iterative method and above-mentioned iterative initial value are used, the soil resistivity solved with tower section is ρtWith the soil of non-same tower section
Earth resistivity is ρf;
Step 5, each parameter value in basic data, by the soil resistivity ρ of same tower section and non-same tower sectiont、ρfRespectively
In the shaft tower for substituting into correspondence model, the impedance matrix of each model shaft tower unit length of each loop line is obtained, and utilize symmetrical component method
The zero-sequence mutual inductance impedance of each model shaft tower unit length in each loop line can be obtained from impedance matrix;
Step 6, zero-sequence mutual inductance impedance of each base shaft tower of every loop line in its span is calculated, by these zero-sequence mutual inductances
Impedance adds up, you can obtain the zero-sequence mutual inductance impedance of whole piece circuit.
Compared with prior art, patent of the present invention has the advantage that as follows:
Circuit zero sequence mutual impedance is calculated using this method, line impedance estimated without rule of thumb coefficient
Calculate, the result calculated considers the influence of tower structure and soil resistivity, with very high accuracy.Using this method meter
The zero-sequence mutual inductance impedance calculated, can greatly improve protection seting calculating, short circuit calculation, Load flow calculation, system operation mode
The levels of precision of the work such as selection.In addition, basic data used in calculating process can be obtained from power department, to base
The processing of plinth data can be carried out using computer, and actual operation process is relatively simple, and the realizability of this method is strong.
Embodiment
Clear, complete description is carried out to the technical scheme in the embodiment of the present invention with reference to the accompanying drawings and examples, is shown
So, described embodiment is only a part of embodiment of the invention, rather than whole embodiments.
For uniform the transposition symmetrical circuit of triphase parameter, i.e. balanced transmission line road during theoretical calculation, its unit length
Positive and negative zero sequence impedance can be calculated (unit is Ω/km) by following formula respectively:
In above formula:R represents D.C. resistance, and 0.05 is used to characterize ground resistance (ground resistance=ZL=π2f×10-4(Ω/
Km), when taking power frequency 50Hz, its value about 0.05), r ' expression wire equivalence radius, DgRepresent that ground returns to the equivalent depth in road, Dave
Represent the geometric mean distance between power transmission line (when considering single loop lineConsider the geometric mean distance between double loop
When,dmnRepresent the space length between correspondence conductor.
In practice, though to same back transmission line, due to along in erection process can across different topography and geomorphologies and
Building etc., soil resistivity along the line is not necessarily worth, but changes with the change of geological conditions, therefore along the line
Soil resistivity is change;In addition, the transmission tower for the construction of line is nor changeless, but with leap side
Formula and it is spanned thing and accordingly changes, therefore the position between conductor can change.Influenceed by geological conditions, soil resistivity
ρ leads to be changed in the range of 100~4000 Ω/m, thus lgDgExcursion 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 in interior fluctuation in a big way.
Because the co-ordination that the measurement of zero-sequence mutual inductance is related to is sufficiently complex, and the erection situation of circuit is limited to, because
Zero-sequence mutual inductance impedance in this engineering generally using empirical method to circuit is estimated, i.e., using a certain fixed coefficient and solid wire
Zero sequence impedance is multiplied to replace the zero sequence mutual impedance between same tower circuit.But from the point of view of accurate calculate, result of calculation and actual value it
Between error there is very big randomness, it is and inadvisable.
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, it includes:
Step 1, obtain and be related to each loop line, each base shaft tower, transmission pressure, the basic data of ground wire;
Step 2, soil it will be equivalent to, with the equally distributed soil resistivity of sectional, make 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 identical soil resistivity, it is the homogeneous of circuit to make the soil resistivity
Soil resistivity ρj, use ρjConstruct ρt、ρfIterative initial value, if having v loop lines with tower section, wherein v is between 3 to 6, then
With tower section soil resistivity ρ on r loop linestWith 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, each loop line is write on ρ by rowt、ρfLinear algebra
Equation group, and iterative method and above-mentioned iterative initial value are used, the soil resistivity solved with tower section is ρtWith the soil of non-same tower section
Earth resistivity is ρf;
Step 5, each parameter value in basic data, by the soil resistivity ρ of same tower section and non-same tower sectiont、ρfRespectively
In the shaft tower for substituting into correspondence model, the impedance matrix of each model shaft tower unit length of each loop line is obtained, and utilize symmetrical component method
The zero-sequence mutual inductance impedance of each model shaft tower unit length in each loop line can be obtained from 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, you can to obtain the zero-sequence mutual inductance value of whole piece circuit.
The basic data, as the specific embodiment of step 1, including:Each loop line is every used in erection process
The actual measurement zero sequence resistance of the phase sequence arrangement mode of wire, each loop line on the structural parameters and its span of one base shaft tower, each base shaft tower
Conductor parameter, the total length of circuit and the line length with tower section of the anti-, radical of transmission pressure parameter, ground wire and ground wire, it is described
Transmission pressure parameter includes division number, line footpath and D.C. resistance.
Tower structure and soil resistivity are to influence the principal element of 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 be surveyed one by one.Therefore in practice, influence zero sequence is mutual
The principal element that sense impedance is accurately calculated is soil resistivity.The step 2 and step 3 are at the equalization to soil resistivity
Reason, the soil resistivity by along carries out equalization processing, is the committed step for calculating all fronts zero-sequence mutual inductance impedance, i.e., by zero sequence
Mutual inductance is thought to be distributed along shaft tower piecewise uniform.
Specifically, after the influence of tower structure is considered, equalization can be expressed as:In the same tower section of transmission line of electricity,
It is equally distributed to think zero-sequence mutual inductance;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 zero-sequence mutual inductance with tower section and non-same tower section, the difference by the soil resistivity of same tower section and non-same tower section Lai
Embody.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 obtains circuit.
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 on ρ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,
Calculation times, calculating time and amount of calculation required for causing in solution procedure are all greatly increased, and practicality is low.According to zero sequence
The experience estimate of mutual inductive impedance, can calculate a corresponding experience soil resistivity ρ0, but ρ0And the real soil of logicalnot circuit
Resistivity, the ρ calculated in some cases0May obviously it not conform in this, as initial value outside 100~4000 Ω/m scopes
It is suitable.Therefore need further to seek ρt、ρfRelatively reasonable iterative initial value.
The concept of homogeneous soil resistivity is introduced below to illustrate how to obtain ρt、ρfIterative initial value.
The influence of tower structure is taken into account, it is believed that completely with identical soil resistivity, calculated with the resistivity
Zero sequence impedance out is equal to the measured value of circuit zero sequence impedance, and the soil resistivity is called the homogeneous soil resistivity of circuit
ρj。
Each parameter value in basic data, soil resistivity is uniformly taken as ρj, each model shaft tower of the loop line can be obtained
The impedance matrix of unit length, and it is long using symmetrical component method each model shaft tower unit in the loop line can be obtained from impedance matrix
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 lines;The non-same tower section numbering s of r loop lines shaft tower correspondence unit
The zero sequence impedance and span of length are expressed as Zrfs(ρjr) and Lrfs, Zrfs(ρjr) it is on the homogeneous electric resistance of soil of r loop lines
Rate ρjrExpression formula;R loop lines are expressed as with the zero sequence impedance and span of tower section numbering s shaft tower correspondence unit length
Zrts(ρjr) and Lrts, Zrfs(ρjr) it is on the homogeneous soil resistivity ρ of r loop linesjrExpression formula.
Constrained by the equation of formula (1), therefore can will be used as ρ0ρjrInitial value is iterated calculating, finally obtains homogeneous soil
Earth electricalresistivityρjExact value.To be better described, illustrate how to solve homogeneous 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 homogeneous soil resistivity is not limited thereto.
1) zero sequence impedance discreet value is calculated
If the experience estimate of the zero-sequence mutual inductance impedance on r loop lines road is ZrM0', and Tm is set as fed distance most long bar
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 try to achieve experience soil resistivity ρ0。
By ρ0It is used as ρrjInitial value substitute into formula (1), obtain the discreet value Z of circuit zero sequence impedance during the 0th corrected Calculationr0
(0), if the zero sequence impedance measured value of r loop lines is Zr0, note deviation ε0=Z0(0)-Z0.According to the actual conditions of circuit, the limits of error
σ is taken as, when | ε0|<During σ, it is believed that ρ0As need the homogeneous soil resistivity solved.|ε0|>During σ, then it is assumed that ρ0Value not
Appropriately, it need to be modified.
2) soil resistivity is corrected
The excursion 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 ρ0Necessarily interval [a,
B] in.The soil resistivity for remembering n-th corrected Calculation is ρn, the circuit A of n-th corrected Calculation zero sequence impedance value 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) zero sequence mutual impedance is proportional to, illustrates that now zero sequence mutual impedance is bigger than normal, should
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.Now, by ρ0B is assigned to, a keeps constant.
If ε0<σ, i.e. Z0(0)<Z0, due to Z0(0) zero sequence mutual impedance is proportional to, illustrates that now zero sequence mutual impedance is less than normal, should
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.Now, by ρ0A is assigned to, b keeps constant.
Make ρ1=(a+b)/2, repeat the above steps calculating Z0(1)、ε1If, | ε1|<σ then stops calculating, by ρ1As homogeneous
Soil resistivity is used to calculate homogeneous zero sequence mutual impedance, and otherwise repeatedly said process calculates ρ respectively2、ρ3、ρ4... until | εj|<
Stop calculating during σ, by ρ nowjIt is used as homogeneous soil resistivity.
3) ρ is obtainedt、ρfIterative initial value
If having v loop lines (6 >=v >=3) with tower section, the homogeneous soil of each loop line is calculated respectively using the above method
Electricalresistivityρj1…ρjv, the iterative initial value on r loop lines with tower section and non-same tower section soil resistivity can then be taken as respectively:(ρj1
+…+ρjv)/v、ρjr。
Referring to Fig. 4, the step 4 calculates non-same tower section or as described below with the specific method of tower section resistivity:
As one embodiment, if multiple-loop line line is using the whole erection mode with tower, with tower section and non-same tower section
Soil resistivity ρt、ρfIt is equal and all be homogeneous soil resistivity ρj, without solving again.By ρjSubstituting into formula can calculate
The zero-sequence mutual inductance impedance of each loop line.
As another embodiment, the transmission line of electricity set up for local multiple-loop line is logical with the inlet wire of tower section and outlet
Often constituted in units 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 generally made up of two groups of double-circuit line on same pole;With the inlet wire of the back transmission line of tower six
(or outlet) is generally made up of three groups of double-circuit line on same pole;Situation with three times erections of tower is less, and now inlet wire (or outlet) is usual
It is made up 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 also inevitable identical, the ρ of the soil resistivity along this two loop linet=ρt2, ρf1=ρf2.Set up jointly in multiple-loop line line
Part, geological conditions along the line is identical, therefore the equal ρ of soil resistivity on each loop line of this section of circuitt1=ρt2=...=
ρtv。
By taking the same loop line of tower four as an example, if I loop lines and II loop lines are one group of double-circuit line on same pole, whole process is joint used;III loop lines
It is two single loop lines with IV loop lines.The zero sequence impedance value surveyed according to each loop line, can arrange and write equation below 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, the zero sequence impedance of I~IV loop lines is represented respectively;ρf1~ρf4, respectively represent I~
The soil resistivity of the non-same tower section of IV loop lines;ρt1~ρt4, soil resistivity of I~IV loop lines with tower section is represented respectively;ρf1~
ρf4, the soil resistivity of the non-same tower section of I~IV loop lines is represented respectively;The non-same tower section numbering s of r loop lines shaft tower correspondence unit is long
The zero sequence impedance and span of degree are expressed as Zrfs(ρfr) and Lrfs, Zrfs(ρfr) it is on the non-same tower section electric resistance of soil of the loop line
Rate ρfrExpression formula;R loop lines are expressed as with the zero sequence impedance and span of tower section numbering s shaft tower correspondence unit length
Zrfs(ρt) and Lrfs, Zrfs(ρt) it is with tower section soil resistivity ρ on the loop linetExpression formula.
To be better described, below illustrate how to solve each loop line with tower section and non-same tower section exemplified by gaussian iteration method
Soil resistivity, calculation process is as shown in Figure 2.But each loop line of solution of the present invention is with the soil resistivity of tower section and non-same tower section
Method be not limited thereto.
Because I loop lines and II loop lines double circuits on same tower are set up, therefore there is ρf1=ρf2., can be by formula (2) by simple conversion
~(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, its value depend on tower section and non-same tower section bar structural parameters and
The corresponding span of each shaft tower.
Formula (6)~(9) are rewritten as to ρ=A ρ+B form, are:
Gaussian iteration form corresponding to formula (9) is:
ρ(k+1)=A ρ(k)+ B, k=0,1,2 ... (11)
ρ initial value is expressed as
Make ρt0=(ρj1+ρj2+ρj3+ρtj4)/4, to assigning 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 working as d<During δ
Stop iteration, now can obtain each loop line with tower section and the accurate solution ρ of non-same tower section soil resistivityt、ρf1(=ρf2)、ρf3、
ρf4。
The step 5 is each parameter value in basic data, asks for the zero sequence mutual impedance between circuit, specific to solve such as
Under, remember respectively:The zero-sequence mutual inductance impedance of the non-same tower section numbering s of r loop lines shaft tower correspondence unit length is Zrmfs(ρfr), Zrmfs
(ρfr) it is on the non-same tower section soil resistivity ρ of the loop linefrExpression formula;R loop lines are single with tower section numbering s shaft tower correspondence
The zero-sequence mutual inductance impedance of bit length is Zrmfs(ρt), Zrmfs(ρt) it is with tower section soil resistivity ρ on the loop linetExpression formula.
The then zero-sequence mutual inductance impedance related to r loop lines be:
Zrm0=Zrmf1(ρfr)Lrf1+…+Zrmfp(ρfr)Lrfp+Zrmt1(ρt)Lrt1…+ZImtq(ρt)Lrtq
The zero-sequence mutual inductance impedance value of each loop line and other loop lines is calculated, the numerical value of the basic data of specific circuit is substituted into
Above formula.
The invention provides a kind of zero-sequence mutual inductance impedance computation method on multiple-loop line transmission line of electricity, using present invention side
Method is calculated the zero sequence mutual impedance of transmission line of electricity, according to the actual conditions founding mathematical models of circuit, the zero sequence calculated
Mutual inductive impedance has considered the actual erection situation and topography variation situation of transmission line of electricity, has higher accuracy, without root
Line impedance is estimated according to empirical coefficient.The zero-sequence mutual inductance impedance that this method is calculated is applied to protect with soil resistivity
Protect in the research work such as adaptive setting, short circuit calculation, Load flow calculation, Computer Simulation, the essence that analysis is calculated can be greatly improved
True degree, is reduced due to the error that experience estimation is caused.And the calculating section of the inventive method can be complete in computer platform
It is simple into, actual mechanical process, with it is very strong can practicality.
Above-described embodiment only not limits technical scheme described in the invention to illustrate the present invention;Therefore, although
This specification is with reference to each above-mentioned embodiment to present invention has been detailed description, still, the ordinary skill of this area
Personnel should be appreciated that still can modify or equivalent substitution to the present invention;And all do not depart from the present invention spirit and
The technical scheme of scope and its improvement, it all should cover among scope of the presently claimed invention.