CN101074979A - Method for measuring power-transmission circuit double-ended distance by distributing parameter - Google Patents

Abstract

A method of using distribution parameter to realize distance-measurement of two ends on power transmission line includes obtaining relevant instant value of current and voltage at line local side by line protection device through sampling, using Fourier algorithm to derive out Fourier from of local side 3-phase current and 3-phase voltage, receiving phase form of line opposite side current and voltage obtained by protection synchrous sampling through communication network and utilizing equality of positive sequence voltage at fault two sides to obtain distance from protection erection to reference point.

Description

Utilize distribution parameter to realize the two ends of electric transmission line distance-finding method

Technical field

The present invention relates to field of power, relate more specifically to the method for electric power system fault range finding, is the method for utilizing the distribution parameter judgement transmission line malfunction position of the phasor current at circuit two ends and voltage and circuit itself.

Background technology

Transmission line of electricity is the basic equipment of electric system generating, conveying electricity etc., occupies important status in electric system.During transmission line malfunction, if can carry out localization of fault fast and accurately, not only help in time to repair faulty line, the assurance electric system is reliably powered, and all very important to the safety and stability and the economical operation of electric system.

At present, the method for measuring distance of transmission line fault can be divided into single end distance measurement and two kinds of methods of both-end distance measuring.Utilize single-ended measurement not to be subjected to the restriction of line channel condition apart from method, but be subjected to the influence of too many factor and cause dysmetria true, just be difficult to eliminate such as the influence of transition resistance, and circuit model is to adopt lumped parameter model that location algorithm is difficult to reach pinpoint accuracy from principle to the former.

At present, the both-end distance measuring algorithm can be divided into based on the location algorithm of lumped parameter model with based on the location algorithm of distributed parameter model.Algorithm based on lumped parameter model can be divided into the consideration distributed capacitance again and not consider two kinds of distributed capacitance.On principle, it is more accurate than adopting lumped parameter model to adopt distributed parameter model, particularly for the high pressure long transmission line.The problem of bringing is the more complicated of finding the solution of equation, and calculated amount is big, for some iterative algorithm, also exists and finds the solution the problem of hyperbolic function, and result of calculation is very likely dispersed, thereby causes the range finding failure.Adopt in the middle of the algorithm of lumped parameter, consider that distributed capacitance is obviously accurate than not considering distributed capacitance, for short-term, can ignore the influence of distributed capacitance, but, not consider that distributed capacitance can bring bigger error for the middle or long line road.

Summary of the invention

For the extra-high voltage long transmission line, because the influence of line distribution capacitance electric current can not be ignored, therefore, traditional is foundation with the lumped parameter, and the algorithm that utilizes single-ended amount to carry out fault localization at the circuit head end can not satisfy the requirement of existing electric system to localization of fault, therefore, need a kind of new, simple Two-terminal Fault Location element is realized both-end distance measuring algorithm more accurately, to satisfy the demand of electric system under existence conditions.

Technical scheme of the present invention is as follows:

At first utilize existing synchronized algorithm that two ends electric current and voltage are adjusted to synchronously, select corresponding positive sequence or negative phase-sequence location algorithm for use according to system failure situation again, simultaneously, can be as the case may be, both are comprehensive, utilize positive sequence, negative phase-sequence combined amount hybrid ranging.

It is as follows to get fault localization by forward-order current and positive sequence voltage:

$X=\frac{1}{2r_1}\ln \frac{(I_{n1}{Z}_{c1}+U_{n1})+{(I_{m1}{Z}_{c1}-U_{m1})e}^{r_{1}l}}{(I_{n1}{Z}_{c1}-U_{n1})+{(I_{m1}{Z}_{c1}+U_{m1})e}^{-r_{1}l}}---\left(1\right)$

To asymmetric capable fault, can also can calculate abort situation by the positive sequence amount by the negative phase-sequence amount;

It is as follows to calculate fault localization by negative-sequence current and voltage:

$X=\frac{1}{2r_1}\ln \frac{(I_{n2}{Z}_{c1}+U_{n2})+{(I_{m2}{Z}_{c1}-U_{m2})e}^{r_{1}l}}{(I_{n2}{Z}_{c1}-U_{n2})+{(I_{m2}{Z}_{c1}+U_{m2})e}^{-r_{1}l}}---\left(2\right)$

As follows by positive sequence and negative phase-sequence comprehensive utilization range finding formula:

$X=\frac{1}{2{\mathrm{\γ}}_1}\ln \frac{({(I_{n1}+I_{n2})Z_{c1}+(U_{n1}+U_{n2}))+((I_{m1}+I_{m2})Z_{c1}-(U_{m1}+U_{m2}))e}^{r_{1}l}}{((I_{n1}+I_{n2})Z_{c1}-\left({U_{n1}+U}_{n2}\right))+{((I_{m1}+I_{m2})Z_{c1}+{(U_{m1}+U}_{m2}\left)\right)e}^{-r_{1}l}}---\left(3\right)$

In the following formula: the positive sequence and the negative sequence component of subscript 1,2 expression institute calculated amount, n, m represent the circuit both sides;

Z _C1Be circuit positive sequence wave impedance,, therefore,, can replace with the positive sequence wave impedance using negative phase-sequence to calculate in the range finding because circuit positive sequence is identical with the negative phase-sequence parameter;

γ ₁Be the propagation constant of circuit, the positive-negative sequence propagation constant of circuit is identical;

I _n, I _m, U _n, U _mThe voltage of representing electric current, n side and the m side of n side and m side respectively;

L represents the protection domain Route Length, need be given in advance, and unit is km;

X is the distance that circuit n end is arrived in the trouble spot, and unit is km;

Wherein circuit being carried out fault localization can calculate range finding respectively by positive sequence amount (formula 1), negative phase-sequence amount (2), also can calculate range finding by the comprehensive amount (3) of positive sequence and negative phase-sequence.

Description of drawings

Fig. 1 has shown route protection and variable position synoptic diagram.

Specific embodiments

As shown in Figure 1, line protective devices installation site and calculating expression amount have been shown.

Wherein: X is the distance that the n end is arrived in the trouble spot, and unit is km, and the positive sequence voltage and the forward-order current at circuit n, m two ends are respectively U _N1, U _M1, I _N1, I _M1, the negative sequence voltage and the negative-sequence current at circuit n, m two ends are respectively U _N2, U _M2, I _N2, I _M2, l represents the protection domain Route Length, is unit with km.

At first utilize ripe sampling markers to adjust synchronized algorithm with the two ends data sync, such as utilize relatively know at present the sampling label with footwork realize the both sides protective device synchronously;

Obtain the Fu Shi form of this side three-phase current and three-phase voltage by fourier algorithm, and the phasor form that receives offside protection synchronized sampling and the electric current that calculates through filtering, voltage by the optical-fibre communications net;

Utilize strain sequence arithmetic to calculate the positive sequence voltage and the electric current U at circuit two ends _N1, U _M1, I _N1, I _M1, negative sequence voltage and electric current U _N2, U _M2, I _N2, I _M2, circuit just (is being born) preface wave impedance Z _C1Just (bearing) preface propagation constant γ ₁Can calculate by following formula according to prior given parameter in advance,

$Z_{c1}=\sqrt{(r_1+\mathrm{j\ω}{L}_1)/(g_1+\mathrm{j\ω}{C}_1)}$ ${\mathrm{\γ}}_1=\sqrt{(r_1+\mathrm{j\ω}{L}_1)*(g_1+\mathrm{j\ω}{C}_1)};$ In the formula: r ₁Be circuit resistance per unit length, L ₁Be circuit unit length induction reactance, g ₁For the electricity of lead unit length is over the ground led (generally smaller, as can to ignore), C ₁Be the circuit capacitance per unit length; Because of the positive sequence and the negative phase-sequence element of circuit are equal, therefore, in negative phase-sequence range finding computing formula, available positive sequence wave impedance Z _C1With positive sequence propagation constant γ ₁

It is as follows to utilize the positive sequence amount to calculate the abort situation formula:

$X=\frac{1}{2r_1}\ln \frac{(I_{n1}{Z}_{c1}+U_{n1})+{(I_{m1}{Z}_{c1}-U_{m1})e}^{r_{1}l}}{(I_{n1}{Z}_{c1}-U_{n1})+{(I_{m1}{Z}_{c1}+U_{m1})e}^{-r_{1}l}}---\left(1\right)$

It is as follows to utilize the negative phase-sequence amount to calculate the abort situation formula:

$X=\frac{1}{2r_1}\ln \frac{(I_{n2}{Z}_{c1}+U_{n2})+{(I_{m2}{Z}_{c1}-U_{m2})e}^{r_{1}l}}{(I_{n2}{Z}_{c1}-U_{n2})+{(I_{m2}{Z}_{c1}+U_{m2})e}^{-r_{1}l}}---\left(2\right)$

In the following formula, remove positive sequence voltage and electric current U _N1, U _M1, I _N1, I _M1, negative sequence voltage and electric current U _N2, U _M2, I _N2, I _M2, wave impedance Z _C1With propagation constant γ ₁Outward, l is a Route Length in the protection domain;

Consideration for improving the precision and the sensitivity of fault localization, can be considered positive sequence when the asymmetry fault, the negative phase-sequence comprehensive utilization, and formula is as follows:

$X=\frac{1}{2{\mathrm{\γ}}_1}\ln \frac{({(I_{n1}+I_{n2})Z_{c1}+(U_{n1}+U_{n2}))+((I_{m1}+I_{m2})Z_{c1}-(U_{m1}+U_{m2}))e}^{r_{1}l}}{((I_{n1}+I_{n2})Z_{c1}-\left({U_{n1}+U}_{n2}\right))+{((I_{m1}+I_{m2})Z_{c1}+{(U_{m1}+U}_{m2}\left)\right)e}^{-r_{1}l}}---\left(3\right)$

For symmetric fault,, therefore, can only calculate with the positive sequence amount owing to there is not the negative phase-sequence amount.

Compare with present used location algorithm, this algorithm is not subjected to line mutual-inductance, not affected by trouble point transition resistance etc., And be not subjected to the impact of system power supply impedance and load current, and for positive sequence and negative phase-sequence amount, after fault, need only fault not Excision is long-term existence.

Compare with present used both-end distance measuring algorithm and since formula in used be standard natural Exponents and logarithm, therefore, From existing Calculating Foundation, existing standard operation. Theory and practice shows that this algorithm has very big proposing to range accuracy High.

Claims (2)

1. One kind in electric system by two ends electric current and voltage and the method for utilizing Transmission Line Distributed Parameter that fault is found range, this method comprises the steps:

   Utilize synchronized algorithm that two ends electric current and voltage are adjusted to synchronously;

   Line protective devices obtain the electric current and voltage instantaneous value to the voltage current waveform sampling of mutual inductor;

   Obtain the phasor form of three-phase current and voltage by fourier algorithm;

   Calculate both sides m by conversion by both sides three-phase current and voltmeter, the negative-sequence current I of n _M2, I _N2, negative sequence voltage U _M2, U _N2With forward-order current I _M1, I _N2With positive sequence voltage U _M1, U _N1

   For symmetrical fault, can calculate abort situation by the positive sequence amount,

   It is as follows to calculate fault localization by forward-order current and positive sequence voltage:

   $X=\frac{1}{{2r}_1}\ln \frac{(I_{n1}{Z}_{c1}+U_{n1})+(I_{m1}{Z}_{c1}-U_{m1})e^{r_{1}l}}{(I_{n1}{Z}_{c1}-U_{n1})+(I_{m1}{Z}_{c1}+U_{m1})e^{-r_{1}l}}$

   For the asymmetry fault, can calculate abort situation by the positive sequence amount, or can calculate abort situation by the negative phase-sequence amount; It is as follows to calculate fault localization by negative-sequence current and negative sequence voltage:

   $X=\frac{1}{{2r}_1}\ln \frac{(I_{n2}{Z}_{c1}+U_{n2})+(I_{m2}{Z}_{\mathrm{cl}}-U_{m2})e^{r_{1}l}}{(I_{n2}{Z}_{c1}-U_{n2})+(I_{m2}{Z}_{\mathrm{cl}}+U_{m2})e^{-r_{1}l}}$

   Wherein: * being the distance of trouble spot to the n end, unit is km, the positive sequence voltage and the forward-order current at circuit n, m two ends are respectively U _N1, U _M1, I _N1, I _M1, the negative sequence voltage and the negative-sequence current at circuit n, m two ends are respectively U _N2, U _M2, I _N2, I _M2, Z _C1Be circuit positive sequence wave impedance, γ 1 is a circuit positive sequence propagation constant, and l represents the protection domain Route Length, is unit with km.
2. 2. method as claimed in claim 1 further comprises, to the asymmetry fault, comprehensive utilization positive sequence and negative phase-sequence are calculated abort situation, wherein:

   As follows by positive sequence and negative phase-sequence comprehensive utilization range finding formula:

   $X=\frac{1}{{2\mathrm{\γ}}_1}\ln \frac{((I_{n1}+I_{n2})Z_{c1}+(U_{n1}+U_{n2}))+((I_{m1}+I_{m2})Z_{\mathrm{cl}}-(U_{m1}+U_{m2}))e^{r_{1}l}}{((I_{n1}+I_{n2})Z_{c1}-(U_{n1}+U_{n2}))+((I_{m1}+I_{m2})Z_{c1}+(U_{m1}+U_{m2}))e^{-r_{1}l}}.$

Images