CN104931849B - A kind of supply line's arc grounding fault distance-finding method - Google Patents
A kind of supply line's arc grounding fault distance-finding method Download PDFInfo
- Publication number
- CN104931849B CN104931849B CN201510243171.1A CN201510243171A CN104931849B CN 104931849 B CN104931849 B CN 104931849B CN 201510243171 A CN201510243171 A CN 201510243171A CN 104931849 B CN104931849 B CN 104931849B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- msup
- mtr
- mtd
- 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.)
- Active
Links
Landscapes
- Locating Faults (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
The invention belongs to field of power, more particularly to a kind of supply line's arc grounding fault distance-finding method.Most of supply line's earth fault belongs to arc grounding, after faulty line power failure, with DC current generator to line charging, raise line-to-earth voltage, when reaching failure breakdown voltage, trouble point is punctured again, and trouble point is changed into rapidly low resistive state by high-impedance state, line distribution capacitance is discharged by breakdown point, while the voltage on circuit drastically declines.Regard circuit as a dynamical system, establish the dynamic model of circuit, with dynamic system parameter discrimination method, according to the electric current on circuit and voltage sample series, identification circuit dynamic model parameters, be out of order distance according to model parameter calculation.
Description
Technical field
The invention belongs to field of power, more particularly to a kind of supply line's arc grounding fault distance-finding method.
Background technology
Supply line's failure more than 80% belongs to earth fault.Earth fault is roughly divided into resistance eutral grounding and arc grounding two
Class.In mesohigh supply line, arc grounding ratio is far above the ratio of resistance eutral grounding.
Existing earth fault distance measurement method is on-line checking mostly, i.e. fault localization equipment is hung on power network always, inspection
After measuring earth fault, the distance measurement function of starting device, completes fault localization immediately.During online fault localization, faulty line is also
Without departing from network system.When electric network composition is more complicated, for example when branch and more user, distance measurement result is often insincere.It is special
It is not, in the faulty line power failure maintenance stage, the ground of the online completely useless force of fault localization equipment.
Existing fault localization principle can be divided into two kinds of impedance method and traveling wave method.Signal injection method, phasor approach, impedance method and
Modulus Analysis etc. all belongs to impedance method.This method calculates line impedance according to the electric current under malfunction and voltage characteristic, so as to
Calculate fault distance.Equipment is simple, easily realizes.But this method there is also some problems in actual applications：1st, alternating current
Press under breakdown condition, the impedance of trouble point is non-linear serious, influences measurement accuracy；2nd, user load produces very big to measurement process
Interference；3rd, faulty line and regular link link together, complicated, cause distance measurement result deviation very big；4th, ground connection event
It is arc light arcing fault to hinder majority of case, and trouble point can not be represented with simple linear resistance.These factors greatly shadow
Ring the ranging effect of impedance method.Various methods based on traveling wave principle, such as single-ended traveling wave, both-end traveling wave, wavelet analysis, small echo
Combined with neutral net, simulated annealing etc., fault distance is calculated by measuring fault traveling wave transmission time, the method
Succeed application on the transmission line of electricity of single loop long range.When the technology is applied to multipoint line such as distribution network systems event
When hindering ranging, effect is unsatisfactory.Traveling wave is influenceed very big, signal multiple reflections fusion in transmitting procedure by lines branch,
Distortion is serious, and the wave head for detecting traveling wave is relatively difficult.In addition, travelling wave ranging equipment is complicated, it is difficult to spreads to mesolow power supply system
System.
The content of the invention
Offline ranging and Dynamic Signal ranging are realized the technical problem to be solved in the present invention is to provide one kind, improve ranging essence
Degree and reliability, the safety of safeguard work personnel, solve to be difficult to maintain fail-steady problem during arc grounding, calculate journey
Supply line's arc grounding fault distance-finding method that sequence is simple, measuring speed is fast.
In order to solve the above technical problems, the present invention adopts the following technical scheme that：
Most of supply line's earth fault belongs to arc grounding.When the voltage on circuit is relatively low, it is impossible to breakdown fault point
When, high-impedance state is presented in trouble point.If after faulty line power failure, with DC current generator to line charging, make circuit pair
Ground voltage is raised, and when reaching failure breakdown voltage, trouble point is punctured again, and trouble point is changed into rapidly low-resistance by high-impedance state
State.Line distribution capacitance can fiercely be discharged by breakdown point, form very big discharge pulse current, while the voltage on circuit
Drastically decline.Failure breakdown process can generally repeat, and the amplitude of variation of electric current and voltage is all bigger in discharge process, these
All it is the factor beneficial to measurement.On the other hand, the failure breakdown process duration is shorter, generally only continues hundreds of microseconds to tens
Millisecond, so trouble point distance can not be calculated to analyze with steady-state model and steady-state analysis method.Circuit can be regarded as to one
Dynamical system, the dynamic model of circuit is established, with dynamic system parameter discrimination method, such as, least square method, Kalman's dynamic
The methods of filtering, according to circuit DC current input electric current and voltage sample series, identification circuit dynamic model parameters,
Then, it is out of order distance according to model parameter calculation.Thus propose that a kind of supply line's arc grounding fault distance-finding method is as follows：
A kind of supply line's arc grounding fault distance-finding method, includes high voltage direct current flow-generator, high-voltage direct current
The output end of generator is connected with current sampling unit and voltage sampling unit, comprises the following steps：
(1) broken down in supply line, after electric power system separation, with high voltage direct current flow-generator to faulty line
DC current is injected, line-to-earth voltage is gradually risen, trouble point is punctured again, current sampling unit and voltage sampling unit
Respectively to the electric current and voltage sample in breakdown process, electric current and voltage sample series, i (0), i (1) ... ... i (N), u are obtained
(0),u(1),……u(N)；
(2) according to the characteristics of faulty line, the dynamic model of faulty line is built, the parameter and fault distance of the model have
Proportional relation, electric current and voltage sample series according to acquired by (1) step, with Dynamical System Identification, estimates dynamic analog
The parameter of type；
To improve range accuracy, according to different line features, the dynamic analog of one or more faulty lines is constructed in advance
Type, it is determined that corresponding Dynamical System Identification, during for ranging selection use；
(3) according to the model parameter of estimation, and the relation of fault distance and model parameter, fault distance is calculated.
The method for building the dynamic model of faulty line is a lot, using the circuit model being made up of resistance, inductance and electric capacity,
Such as single order fault equivalence circuit, second order fault equivalence circuit or three rank fault equivalence circuits, or faults line current
The transfer function model of voltage relationship.
Dynamical System Identification is also a lot, and conventional has least square method, auxiliary variable method, Kalman filtering method.
Preferably, when the dynamic model of the faulty line is three rank fault equivalence circuit, using Dynamic System Identification side
Method and fault distance are calculated as follows：
Discrete dynamic system model such as following formula is built according to three rank fault equivalence circuits,
U (n)=- a1u(n-1)-a2u(n-2)-a3u(n-3)
+b0i(n)+b1i(n-1)+b2i(n-2)+b3i(n-3)
Parameter θ=[a of the dynamic system model is recognized with auxiliary variable method1a2a3b0b1b2b3]T,
The parameter of gained is substituted into following system of linear equations：
WhereinT is the sampling period, and solving equations obtain parameter ЬS0、ЬS1、ЬS2、аS1、аS2And аS3；
By the parameter Ь of gainedS0、ЬS1、ЬS2、аS1、аS2And аS3Substitute into below equation group：
In above equation group, L is Current injection points to the line equivalent inductance between failure breakdown point, with fault distance into
Direct ratio；R0It is that Current injection points are directly proportional to the line equivalent resistance between failure breakdown point, and fault distance；R1It is that failure is hit
Ground connection transition resistance after wearing, C0It is that line capacitance converts the equivalent capacity for arriving Current injection points, C1It is that event is arrived in line capacitance conversion
Hinder the equivalent capacity of point.
Solution above equation group obtains：
The distance of L divided by circuit unit length inductance, as Current injection points to trouble point；Or by above-mentioned equation
Group solves R0, R0Divided by the distance of the resistance of circuit unit length, as Current injection points to trouble point.
The present invention compared with prior art, has the advantage that and beneficial effect：
Compared to the prior art, the present invention is main proposes 2 essence improvement：It is improved to offline by online ranging
Ranging, Dynamic Signal ranging is improved to by steady-state signal ranging.
For the safety of safeguard work personnel, after line failure, when repairing failure, always first to have a power failure what is overhauled again,
Off-line type fault localization more meets the requirement of real work.After faulty line departs from power network, interference substantially reduces, and is advantageous to improve
Range accuracy.Off-line type ranging can be repeated several times, and particularly measurement point may move, by different place rangings, progressively forcing
Nearly trouble point, reliability are more preferable.
The present invention uses Dynamic Signal ranging, when overcoming arc grounding, it is difficult to maintains fail-steady problem, holds
Easily realize, calculation procedure is simple, and measuring speed is fast, and precision is high.
Brief description of the drawings
Fig. 1 is the operation principle schematic diagram that offline ranging is carried out using the inventive method.
In Fig. 1,1 is high voltage direct current flow-generator, and 2 be current sampling unit, and 3 be voltage sampling unit, and 4 be control meter
Display unit is calculated, 5 represent the faulty line after having a power failure, and 6 represent earth fault.
Fig. 2 is the single order fault equivalence circuit diagram of faulty line in the present invention.
Fig. 3 is the second order fault equivalence circuit diagram of faulty line in the present invention.
Fig. 4 is three rank fault equivalence circuit diagrams of faulty line in the present invention.
Embodiment
The present invention is elaborated below in conjunction with the drawings and specific embodiments, but not formed to the claims in the present invention
The limitation of protection domain.
As shown in figure 1, broken down in supply line, after electric power system separation, with high voltage direct current flow-generator 1 to
Faulty line 5 injects DC current, is stepped up line-to-earth voltage, trouble point 6 is punctured again, the He of current sampling unit 2
Voltage sampling unit 3 to the electric current and voltage sample in breakdown process, obtains electric current and voltage sample series, i (0), i respectively
(1),……i(N),u(0),u(1),……u(N)；Control calculating display unit 4 is serial according to electric current and voltage sample, according to
The faulty line dynamic model being previously set, using corresponding dynamic system parameter discrimination method, computation model parameter, finally by
Model parameter calculation fault distance, and result is shown.
Embodiment 1：
It is the single order fault equivalence circuit diagram of faulty line as shown in Figure 2.In figure, L is that Current injection points puncture to failure
Line equivalent inductance between point, and fault distance are directly proportional；R0Be Current injection points to the circuit between failure breakdown point etc.
Effect resistance, and fault distance are directly proportional；R1It is the ground connection transition resistance after failure breakdown.
Laplace transform is carried out to single order fault equivalence circuit, obtains its transmission function, after discretization, obtains discrete arteries and veins
Transmission function is rushed, with Dynamical System Identification, the parameter of pulsed transfer function is tried to achieve, continuous system is tried to achieve by solving equation
Load transfer function coefficient, and then equivalent inductance L, L in equivalent circuit divided by the inductance of faulty line unit length are calculated, produce
Fault distance.
Embodiment 2：
It is the second order fault equivalence circuit diagram of faulty line as shown in Figure 3.In figure, L is that Current injection points puncture to failure
Line equivalent inductance between point, and fault distance are directly proportional；R0Be Current injection points to the circuit between failure breakdown point etc.
Effect resistance, and fault distance are directly proportional；R1Be failure breakdown after ground connection transition resistance, C1It is that faulty line is converted to trouble point
Equivalent capacity.
Laplace transform is carried out to second order fault equivalence circuit, obtains its transmission function, after discretization, obtains discrete arteries and veins
Transmission function is rushed, with Dynamical System Identification, the parameter of pulsed transfer function is tried to achieve, continuous system is tried to achieve by solving equation
Load transfer function coefficient, and then calculate the equivalent inductance L in equivalent circuit or equivalent resistance R0, L divided by faulty line unit are grown
The inductance of degree, fault distance is produced, or, R0Divided by the resistance of faulty line unit length, also obtain fault distance.
Embodiment 3：
It is three rank fault equivalence circuit diagrams of faulty line as shown in Figure 4.In figure, L is that Current injection points puncture to failure
Line equivalent inductance between point, and fault distance are directly proportional；R0Be Current injection points to the circuit between failure breakdown point etc.
Effect resistance, and fault distance are directly proportional；R1Be failure breakdown after ground connection transition resistance, C0It is that line capacitance conversion is noted to electric current
The equivalent capacity of access point, C1It is equivalent capacity of the line capacitance conversion to trouble point.The line capacitance includes each branch line
The line capacitance between circuit and trouble point and line end between road, trouble point and Current injection points.
Laplace transform is carried out to three rank fault equivalence circuits, obtains its transmission function, after discretization, obtains discrete arteries and veins
Transmission function is rushed, with Dynamical System Identification, the parameter of pulsed transfer function is tried to achieve, continuous system is tried to achieve by solving equation
Load transfer function coefficient, and then calculate the equivalent inductance L in equivalent circuit or equivalent resistance R0, L divided by faulty line unit are grown
The inductance of degree, produces fault distance；Or R0Divided by the resistance of faulty line unit length, also obtain fault distance.
Fault equivalence circuit and common line equivalent circuit shown in Fig. 4 make a big difference, and add ground connection transition
Resistance R1, particularly C0And C1It is separate variable.Substantial amounts of emulation and Physical Experiment show, are no more than in circuit overall length
During 40 km, in direct current to 2000 frequency ranges, the frequency of fault equivalence circuit and physical fault circuit shown in Fig. 4
Characteristic error is no more than 2%, can meet the actual requirement of engineering.
Discretization method described in embodiment 1-3 is constant using first difference method, Bilinear transformation method, impulse response
Method or step response not political reform.
Dynamical System Identification described in embodiment 1-3 is using least square method, auxiliary variable method, Kalman filtering
Method or maximum-likelihood method.
The advantages of single order fault equivalence circuit is simple in construction, and unknown number is few, and calculating speed is fast, and shortcoming is that error is larger.
Three rank fault equivalence circuit structures are complicated, and parameter is more, and its advantage is precision height.
Below by taking three rank fault equivalence circuits as an example, illustrate the calculating process of the present invention.
Laplace changes are carried out to three rank fault equivalence circuits, obtain the functional relation of input voltage and input current such as
Shown in following formula (1)：
In formula (1)：
Discretization is carried out to formula (1) and obtains following formula (3)：
B in formula (3)3、b2、b1、b0、a3、a2、a1Meet below equation group (4)：
In equation group (4)T is the sampling period.
Shown in the difference equation such as following formula (5) for writing out input voltage u and input current i according to formula (3)：
According to formula (5) and current sample series i (0), i (1) ... i (N), voltage sample series u (0), u (1) ... u
(N), the relational expression of input and output is designated as shown in matrix form such as following formula (6)：
U=A θ (6)
θ=[a in formula (6)1 a2 a3 b0 b1 b2 b3]TFor parameter to be estimated,
U=[u (4) u (5) u (6) u (N+3)]TIt is serial for output,
To recognize equation calculation matrix.
According to formula (6), a kind of θ parameter Estimation is obtained using least square method.
When least-squares parameter estimation algorithm is applied to Dynamic System Identification, algorithm is simple, but result be usually have it is inclined.
The unbiased parameter estimation of formula (3), auxiliary variable method parameter Estimation formula such as following formula can be obtained using auxiliary variable method
(7) shown in：
WhereinFor systematic parameter θ unbiased estimator, B is auxiliary matrix of variables.Auxiliary variable battle array B is by output estimation valueFormed with input quantity i：
Auxiliary variable method parameter Estimation is an iterative process, and step is as follows：
1st, the initial value of auxiliary variable method estimation is determined with least square methodB, iterations K initial value are 0；
2nd, byTry to achieveOrderAuxiliary variable is tried to achieve according to formula (8)
Battle array B (K+1)；
3rd, new parameter Estimation is tried to achieve by formula (7)IfThenEstimate for required parameter
Evaluation, stop iteration.Otherwise makeB (K)=B (K+1) goes to step 2；
With least square method or auxiliary variable method, after obtaining the estimation parameter of formula (3), the estimation parameter is substituted into equation
Group (4), with the Solving Linear algorithm of standard, the parameter Estimation of formula (1) is obtained, by the parameter Estimation substitution side of formula (1)
Journey group (2) can solve line inductance L
The distance of inductance value L divided by circuit unit length inductance, as trouble point.
Similarly, if the parameter Estimation of equation group (1) is substituted into equation group (2), line resistance R is solved0, resistance value R0Divided by
The distance of the resistance of circuit unit length, as trouble point.
Claims (4)
1. a kind of supply line's arc grounding fault distance-finding method, include high voltage direct current flow-generator, high-voltage direct current hair
The output end of raw device connects current sampling unit and voltage sampling unit respectively, it is characterised in that comprises the following steps：
(1) break down in supply line, after electric power system separation, injected with high voltage direct current flow-generator to faulty line
DC current, line-to-earth voltage is raised, trouble point is punctured again, current sampling unit and voltage sampling unit are respectively to hitting
Electric current and voltage sample during wearing, respectively obtain current sample series i (0), i (1) ... ... i (N), and voltage sample system
U (0), u (1) are arranged ... u (N)；
(2) according to the characteristics of faulty line, the dynamic model of faulty line is built, the parameter and fault distance of the dynamic model have
Proportional relation, electric current and voltage sample series according to acquired by (1) step, the dynamic analog is estimated with Dynamical System Identification
The parameter of type；
(3) according to the model parameter of estimation, and the relation of fault distance and model parameter, fault distance is calculated.
2. supply line's arc grounding fault distance-finding method according to claim 1, it is characterised in that：The faulty line
Dynamic model, be faulty line equivalent circuit either faults line current voltage relationship transfer function model.
3. supply line's arc grounding fault distance-finding method according to claim 1 or 2, it is characterised in that：The dynamic
Model is single order fault equivalence circuit, second order fault equivalence circuit or the three rank fault equivalence circuits of faulty line.
4. supply line's arc grounding fault distance-finding method according to claim 3, it is characterised in that：The faulty line
Dynamic model when being three rank fault equivalence circuits, the Dynamical System Identification and fault distance are calculated as follows：
Discrete dynamic system model such as following formula is built according to three rank fault equivalence circuits,
U (n)=- a1u(n-1)-a2u(n-2)-a3u(n-3)
+b0i(n)+b1i(n-1)+b2i(n-2)+b3i(n-3)
Parameter θ=[a of the dynamic system model is recognized with auxiliary variable method1 a2 a3 b0 b1 b2 b3]T,
The parameter of gained is substituted into following system of linear equations：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>b</mi>
<mrow>
<mi>s</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>cb</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msub>
<mi>b</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<msub>
<mi>ca</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<msup>
<mi>c</mi>
<mn>3</mn>
</msup>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>3</mn>
<msub>
<mi>b</mi>
<mrow>
<mi>s</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>cb</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msub>
<mi>b</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<msub>
<mi>ca</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<msup>
<mi>c</mi>
<mn>3</mn>
</msup>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>3</mn>
<msub>
<mi>b</mi>
<mrow>
<mi>s</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>cb</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msub>
<mi>b</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<msub>
<mi>ca</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<msup>
<mi>c</mi>
<mn>3</mn>
</msup>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>b</mi>
<mrow>
<mi>s</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>cb</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msub>
<mi>b</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>3</mn>
</msub>
<msub>
<mi>ca</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>3</mn>
</msub>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>3</mn>
</msub>
<msup>
<mi>c</mi>
<mn>3</mn>
</msup>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>b</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mi>c</mi>
<mo>-</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msup>
<mi>c</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<mn>3</mn>
<msup>
<mi>c</mi>
<mn>3</mn>
</msup>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mn>3</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mi>c</mi>
<mo>+</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<msup>
<mi>c</mi>
<mn>3</mn>
</msup>
<mo>-</mo>
<mn>3</mn>
<msup>
<mi>c</mi>
<mn>3</mn>
</msup>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mn>3</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<mi>c</mi>
<mo>+</mo>
<mi>c</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<msup>
<mi>c</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<msup>
<mi>c</mi>
<mn>3</mn>
</msup>
<mo>+</mo>
<msup>
<mi>c</mi>
<mn>3</mn>
</msup>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
WhereinT is the sampling period, and solving equations obtain parameter bs0、bs1、bs2、аs1、аs2And аs3；
By the parameter b of gaineds0、bs1、bs2、аs1、аs2And аs3Substitute into below equation group：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mi>L</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mi>L</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mrow>
<mi>s</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<msub>
<mi>C</mi>
<mn>0</mn>
</msub>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mn>0</mn>
</msub>
<mi>L</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>C</mi>
<mn>0</mn>
</msub>
<msub>
<mi>R</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mn>0</mn>
</msub>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<msub>
<mi>C</mi>
<mn>0</mn>
</msub>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mi>L</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
In above equation group, L is Current injection points to the line equivalent inductance between failure breakdown point, with fault distance into just
Than；R0It is that Current injection points are directly proportional to the line equivalent resistance between failure breakdown point, and fault distance；R1It is failure breakdown
Ground connection transition resistance afterwards, C0It is that line capacitance converts the equivalent capacity for arriving Current injection points, C1It is that failure is arrived in line capacitance conversion
The equivalent capacity of point；
Solution above equation group obtains：
<mrow>
<mi>L</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>b</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>b</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>b</mi>
<mrow>
<mi>s</mi>
<mn>0</mn>
</mrow>
</msub>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>s</mi>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mfrac>
</mrow>
The distance of L divided by circuit unit length inductance, as Current injection points to trouble point；Or by above-mentioned solution of equations
Go out R0, R0Divided by the distance of the resistance of circuit unit length, as Current injection points to trouble point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510243171.1A CN104931849B (en) | 2015-05-13 | 2015-05-13 | A kind of supply line's arc grounding fault distance-finding method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510243171.1A CN104931849B (en) | 2015-05-13 | 2015-05-13 | A kind of supply line's arc grounding fault distance-finding method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104931849A CN104931849A (en) | 2015-09-23 |
| CN104931849B true CN104931849B (en) | 2018-01-26 |
Family
ID=54119107
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510243171.1A Active CN104931849B (en) | 2015-05-13 | 2015-05-13 | A kind of supply line's arc grounding fault distance-finding method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104931849B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106646143A (en) * | 2016-12-20 | 2017-05-10 | 国网河北省电力公司电力科学研究院 | Method and system for puncturing high-resistance fault point of cable based on charge superposition principle |
| CN107064734A (en) * | 2017-03-17 | 2017-08-18 | 北京交通大学 | A kind of flexible direct current Fault Location for Distribution Network method of utilization fault transient process |
| CN108957225B (en) * | 2018-06-08 | 2020-08-18 | 西安理工大学 | Direct-current distribution line single-end fault location method considering cable distribution capacitance |
| CN113552443B (en) * | 2021-06-08 | 2022-12-02 | 广西大学 | Hybrid connection high-voltage line fault identification method based on alternating current-direct current pulse and random forest |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102096019A (en) * | 2009-12-15 | 2011-06-15 | 黄洪全 | Method and device for locating single-phase grounding fault of low-current grounding system |
| CN102455399A (en) * | 2010-10-18 | 2012-05-16 | 黄洪全 | Tracking and locating device for earth fault of overhead line |
| CN103529356A (en) * | 2012-07-03 | 2014-01-22 | Abb技术有限公司 | Method and apparatus for determining the distance to phase-to-earth fault |
| CN103743996A (en) * | 2013-10-15 | 2014-04-23 | 昆明理工大学 | Pi-type equivalent circuit based direct current earth electrode line fault location method |
| CN103743995A (en) * | 2013-10-15 | 2014-04-23 | 昆明理工大学 | High-voltage direct current earth electrode line fault location method based on T-type equivalent model |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100821702B1 (en) * | 2006-09-21 | 2008-04-14 | 전명수 | Trolley line impedance detector and Method for locating fault distance using the same |
-
2015
- 2015-05-13 CN CN201510243171.1A patent/CN104931849B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102096019A (en) * | 2009-12-15 | 2011-06-15 | 黄洪全 | Method and device for locating single-phase grounding fault of low-current grounding system |
| CN102455399A (en) * | 2010-10-18 | 2012-05-16 | 黄洪全 | Tracking and locating device for earth fault of overhead line |
| CN103529356A (en) * | 2012-07-03 | 2014-01-22 | Abb技术有限公司 | Method and apparatus for determining the distance to phase-to-earth fault |
| CN103743996A (en) * | 2013-10-15 | 2014-04-23 | 昆明理工大学 | Pi-type equivalent circuit based direct current earth electrode line fault location method |
| CN103743995A (en) * | 2013-10-15 | 2014-04-23 | 昆明理工大学 | High-voltage direct current earth electrode line fault location method based on T-type equivalent model |
Non-Patent Citations (2)
| Title |
|---|
| 利用暂态信号的小电流接地故障距离识别技术;马士聪;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20090515(第05期);C042-70 * |
| 基于注入法的小电流接地系统故障定位新方法;张慧芬 等;《电力系统自动化》;20040210;第28卷(第3期);第64-66页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104931849A (en) | 2015-09-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104931849B (en) | A kind of supply line's arc grounding fault distance-finding method | |
| CN104898021B (en) | A kind of distribution network fault line selection method based on k means cluster analyses | |
| CN105223466B (en) | It is a kind of using modulus maximum than extra high voltage direct current transmission line method of single end distance measurement | |
| CN101587159B (en) | Power distribution network outgoing feeder fault route selecting method by S transform amplitude detection | |
| CN103278756B (en) | A kind of method assessing transformer oil paper insulation ageing state | |
| CN103941149B (en) | Electrified railway traction networks localization of fault method | |
| CN102005755B (en) | Extra-high voltage direct current transmission line boundary element forming method based on support vector machine | |
| CN105116295A (en) | Direct distribution overhead line fault range finding method based on traveling wave abrupt change distance calibration | |
| CN104880648B (en) | The distance-finding method of overhead line one-phase earthing failure in electric distribution network | |
| CN102967779B (en) | Identifying method of distribution parameters of transmission line | |
| CN108107319A (en) | A kind of multiterminal flexible direct current electric network fault localization method and system | |
| CN103777115B (en) | Electric transmission line single-terminal positioning method based on fault transient state and steady-state signal wave velocity difference | |
| CN103760460B (en) | A kind of shape high voltage DC earthing pole line time-domain fault distance-finding method based on Bei Jielong model | |
| CN104865495B (en) | Feeder line fault travelling wave positioning method and system | |
| CN101252275A (en) | Coupling double line down distance measuring method using six-phase system new phase mode transformation | |
| CN102508086B (en) | Parameter estimation method for electric power circuit in T-type connection | |
| CN106019079A (en) | Novel double end fault location method for double DC circuits on same tower | |
| CN109358268A (en) | A kind of accurate fault distance-finding method of double-circuit line on same pole independent of line parameter circuit value | |
| CN109100605A (en) | Utilize the single end positioning method of the high-tension cable singlephase earth fault of failure boundary condition | |
| CN105372559A (en) | Chain type power transmission line single terminal fault location limited continuation method based on characteristic of fault traveling wave distribution along line | |
| CN106019080A (en) | Line-side energy mutation based single-end travelling wave fault location method for double DC circuits on same tower | |
| CN107037325A (en) | Method, device and system for fault line selection of arc suppression coil grounding system | |
| CN103762594A (en) | Feeder line loop closing impact current calculation method based on clock synchronization data matrix | |
| CN112485601B (en) | Fault analysis method and system based on double-end line electrical quantity information | |
| CN105891672A (en) | UPFC-contained line single-end traveling wave distance measurement method based on fault traveling wave line distribution characteristics |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |