CN105699855A - Single-ended traveling fault location calculation method and location method insusceptible to traveling wave speed - Google Patents
Single-ended traveling fault location calculation method and location method insusceptible to traveling wave speed Download PDFInfo
- Publication number
- CN105699855A CN105699855A CN201610209799.4A CN201610209799A CN105699855A CN 105699855 A CN105699855 A CN 105699855A CN 201610209799 A CN201610209799 A CN 201610209799A CN 105699855 A CN105699855 A CN 105699855A
- Authority
- CN
- China
- Prior art keywords
- fault
- traveling wave
- wave
- moment
- arrive
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/085—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution lines, e.g. overhead
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Locating Faults (AREA)
Abstract
The invention provides a single-ended traveling fault location method insusceptible to traveling wave speed and a single-ended traveling fault location calculation method insusceptible to traveling wave speed. According to the invention, the fault occurrence time is set to be t0; the time when initial current traveling waves reach an M end is set to be T1M, the time when first or second type travelling waves reach the M end is set to be t2M, and the time when first, second or third type travelling waves reach the M end is set to be t2M, t4M...; L is set to be the total length of a line; x is set to be the distance between a fault point and a bus M; when 0<=x<=L/3, L/2<=x<=2L/3 and 2L/3<=x<=L, calculation is performed in said three conditions and only incidence relation with time exists while no incidence relation with traveling wave speed exists. Therefore, limits caused by implementation conditions are removed and location accuracy is improved.
Description
Technical field
Patent of the present invention relates to a kind of single-ended traveling wave fault location computational methods, especially a kind of based on the single-ended traveling wave fault location computational methods not affected by traveling wave speed and distance-finding method。
Background technology
After power circuit breaks down, judging the trouble point of power circuit exactly, be conducive to rapidly removing faults, recovering line powering, therefore single-ended traveling wave fault location is a kind of important high-voltage line protection method。In existing single-ended traveling wave fault location computational methods, utilize the traveling wave spread speed estimated to calculate fault distance, due to the deviation of wave velocity, cause that the degree of accuracy of range finding is low。
Summary of the invention
The object of the present invention is a kind of based on the single-ended traveling wave fault location computational methods not affected by traveling wave speed degree;
The object of the present invention is a kind of based on the single-ended traveling wave fault location method not affected by traveling wave speed degree。
In order to overcome above-mentioned technical disadvantages, the purpose of patent of the present invention is to provide based on the single-ended traveling wave fault location computational methods not affected by traveling wave speed and distance-finding method, is therefore no longer limited by implementation condition, improves the degree of accuracy of range finding。
For reaching above-mentioned purpose, patent of the present invention adopts the technical scheme that:
A kind of based on the single-ended traveling wave fault location computational methods not affected by traveling wave speed;Its step: during line failure, is sent fault traveling wave to M bus, N bus by trouble point simultaneously。If fault occurs the moment to be t0, the moment of initial current traveling wave arrival M end is t1M, thereafter first or second class traveling wave arrive M end moment be t2M, and first, second or the 3rd class traveling wave arrive M end moment be t3M、t4M..., if L is the total length of circuit, x is the distance that bus M is arrived in trouble point,
When 0 < x < during L/3:
Usually, if tpIn (p >=4) moment, the transmitted wave of the 1st echo of opposite end bus arrives bus monitoring side,
t0=(3t1-t2)/2
When L/3 < x < during L/2:
t0=(3t1-t2)/2
When L/2 < x < during 2L/3:
t0=(3t1-t3)/2
As 2L/3, < x is < during L
t0=(3t1-tp)/2
Owing to devising point computational methods for different sections, the time only arriving monitoring side with traveling wave surge is relevant, no longer has relation with traveling wave spread speed, is therefore no longer limited by implementation condition, improves the degree of accuracy of range finding。
The present invention devises based on the single-ended traveling wave fault location computational methods not affected by traveling wave speed, its step:
3.1 phase-model transformations
There is coupling in three-phase high-voltage, medium-voltage line, in order to eliminate the impact of coupling between three-phase, need first traveling-wave component to be carried out phase-model transformation, it is separate 0, α, β component by the phase component decoupling intercoupled, converting according to triumphant logical sequence Bel (Karrenbauer), three-phase current decoupling is:
Wherein: Iα、Iβ、I0Respectively phase current Ia、Ib、IcTriumphant human relations Bel convert under α mold component, β mold component and zero _exit, wherein zero _exit is propagated between three-phase conductor and the earth, and α mold component is propagated between A phase and B phase line, and β mold component is propagated between A phase and C phase line。α mold component and β mold component are only propagated between the conductors, so also known as Aerial mode component。
Zero _exit is very serious along with frequency raises decay, does not generally select zero _exit as the object of wavelet transformation;And line mould transient state travelling wave signal can ensure that enough sensitivity as measurement signal, it is possible to carry out wavelet analysis as measuring signal。Adopt the α mold component of fault current as measuring signal。
3.2 location algorithms
Fault current travelling wave signal is gathered at M end。MN circuit diverse location fault, traveling wave process is as Figure 2-Figure 5。
When line failure, trouble point send fault traveling wave to M bus, N bus simultaneously。If fault occurs the moment to be t0, the moment of initial current traveling wave arrival M end is t1M, thereafter first or second class traveling wave arrive M end moment be t2M, and first, second or the 3rd class traveling wave arrive M end moment be t3M、t4M……
If L is the total length of circuit, x is the distance that bus M is arrived in trouble point。
According to arrive M end the 2nd traveling wave belonging to dissimilar, the total length L of circuit can be divided into (0, L/2), (L/2, L) two sections;According to arrive M end the 3rd traveling wave belonging to dissimilar, (0, L/2) section is sub-divided into (0, L/3), (L/3, L/2) two sub-segments, (L/2, L) section is sub-divided into (L/2,2L/3), (2L/3, L) two sub-segments。In sum, the total length L of circuit, it is possible to be divided into 4 sub-segments by L/3, L/2,2L/3 point, respectively 0 < x < L/3, L/3 < x < L/2, L/2 < x < 2L/3,2L/3 < x < L。
Three (four) the individual moment that wavefront arrives M end are utilized to calculate traveling wave speed and fault distance。
Each section location algorithm is derived as follows。
3.2.10 < x < L/3 section
At 0 < x < L/3 section, fault transient travelling wave is as in figure 2 it is shown, simultaneous below equation:
In formula, traveling wave propagates average speed v, moment t occurs fault0, x be unknown number, t1、t2、t3Trying to achieve with the wavelet modulus maxima correspondence moment, observe visible, above simultaneous equations are mutually ratio, it is impossible to try to achieve v, x, from grid chart 2, introduce second class traveling wave due in t4, below equation group can be obtained:
Simultaneous above equation, tries to achieve:
t0=(3t1-t2)/2
Usually, if tpIn (p >=4) moment, the transmitted wave of the 1st echo of opposite end bus arrives bus monitoring side, then exist:
3.2.2L/3 < x < L/2 section
At L/3 < x < L/2 section, fault transient travelling wave as it is shown on figure 3,
At t≤3 τ time interval, list equation:
Simultaneous above equation, solves:
t0=(3t1-t2)/2
3.2.3L/2 < x < 2L/3 section
At L/2 < x < 2L/3 section, fault transient travelling wave as shown in Figure 4,
At t≤3 τ time interval, list equation:
Simultaneous above equation, solves:
t0=(3t1-t3)/2
3.2.42L/3 < x < L section
At 2L/3 < x < L section, fault transient travelling wave as it is shown in figure 5,
At t≤3 τ time interval, list equation:
Simultaneous above equation, without solving。
Utilize the 1st class traveling wave due in tp, obtain below equation:
Simultaneous above equation, solves:
t0=(3t1-tp)/2
Result according to above 4 sectional analyses, in fault distance function x expression formula, independent variable is without row wave speed v, only with total track length and t1、t2、t3、t4(tp) etc. relevant, it is possible to be directly used in circuit fault location range finding。
The present invention devises: a kind of based on the single-ended traveling wave fault location method not affected by traveling wave speed,
Its step:
A, it is three class buses and opposite end is that a class bus carries out traveling wave record ripple and starts according to measuring end local terminal, detection primary fault current traveling wave and carry out triumphant human relations Bel conversion;
B, to IαCarry out One Dimension Continuous Wavelet Transform, extract t1、t2、t3、t4(tp);
C, corresponding wavelet modulus maxima moment t1、t2、t3、t4(tp), read wavelet coefficient polarity, size, it is determined that wavefront polarity;
D, carrying out combinations of polarities according to following principle: with the polarity of initial wavefront wavelet coefficient for reference polarity, if recording first, to arrive wave head wavelet coefficient be negative, then second arrival wave head wavelet coefficient be then "+";If recording first to arrive wave head wavelet coefficient for just, then second arrives wave head wavelet coefficient then respectively "-";
If e, when record first arrive wave head wavelet coefficient for just with the second arrive wave head wavelet coefficient be "-" time, then x < L/2,
According to:
Formula one:
Determine whether to meet 0 < x < L/3,
Formula two:
Determine whether to meet L/3 < x < L/2;
F, when record first arrive wave head wavelet coefficient for negative, second arrive wave head wavelet coefficient for "+" time, then x > L/2,
According to:
Formula one:
Determine whether to meet L/2 < x < 2L/3,
Formula two:
Determine whether to meet 2L/3 < x < L。
In the technical program, sectional based on triumphant logical sequence Bel conversion is calculated as important technical characteristic, based on the single-ended traveling wave fault location not affected by traveling wave speed degree and distance-finding method, having novelty, creativeness and practicality, the term in the technical program is all can make an explanation with patent documentation in the art and understand。
Accompanying drawing explanation
In order to be illustrated more clearly that the embodiment of the present invention or technical scheme of the prior art, the accompanying drawing used required in embodiment or description of the prior art will be briefly described below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the premise not paying creative work, it is also possible to obtain other accompanying drawing according to these accompanying drawings。
Fig. 1 is the flow chart based on the single-ended traveling wave fault location method not affected by traveling wave speed degree of invention;
Fig. 2 is for inventing at 0 < x < L/3 section, fault transient travelling wave figure;
Fig. 3 is for inventing at L/3 < x < L/2 section, fault transient travelling wave figure;
Fig. 4 is for inventing at L/2 < x < 2L/3 section, fault transient travelling wave figure;
Fig. 5 is for inventing at 2L/3 < x < L section, fault transient travelling wave figure;
Fig. 6 is the faulted phase current transients fault compoment oscillogram of 3 circuits (containing faulty line interior) on the same bus in Kang Jin side in inventive embodiments;
Fig. 7 is the wavelet transformation curve chart of invention Fig. 6;
Fig. 8 is that in invention example, B phase earth fault, wherein Kang Jin side faulty line current temporary state fault component oscillogram occur Kang Sui A-wire;
Fig. 9 is the wavelet transformation curve chart of invention Fig. 8。
Detailed description of the invention
Below in conjunction with embodiment, the present invention being further described, following example are intended to illustrate the present invention rather than limitation of the invention further。
A kind of based on the single-ended traveling wave fault location computational methods not affected by traveling wave speed degree, its step:
Single-ended traveling wave fault location algorithm
3.1 phase-model transformations
There is coupling in three-phase high-voltage, medium-voltage line, in order to eliminate between three-phase the impact of coupling, it is necessary to first traveling-wave component is carried out phase-model transformation, is separate 0, α, β component by the phase component decoupling intercoupled。Converting according to triumphant logical sequence Bel (Karrenbauer), three-phase current decoupling is:
Wherein: Iα、Iβ、I0Respectively phase current Ia、Ib、IcTriumphant human relations Bel convert under α mold component, β mold component and zero _exit。Wherein zero _exit is propagated between three-phase conductor and the earth, and α mold component is propagated between A phase and B phase line, and β mold component is propagated between A phase and C phase line。α mold component and β mold component are only propagated between the conductors, so also known as Aerial mode component。
Zero _exit is very serious along with frequency raises decay, does not generally select zero _exit as the object of wavelet transformation;And line mould transient state travelling wave signal can ensure that enough sensitivity as measurement signal, it is possible to carry out wavelet analysis as measuring signal。Adopt the α mold component of fault current as measuring signal。
3.2 location algorithms
Fault current travelling wave signal is gathered at M end。MN circuit diverse location fault, traveling wave process is as Figure 2-Figure 5。
From Fig. 2-Fig. 5, when line failure, trouble point send fault traveling wave to M bus, N bus simultaneously。If fault occurs the moment to be t0, the moment of initial current traveling wave arrival M end is t1M, thereafter first or second class traveling wave arrive M end moment be t2M, and first, second or the 3rd class traveling wave arrive M end moment be t3M、t4M……
If L is the total length of circuit, x is the distance that bus M is arrived in trouble point。
According to arrive M end the 2nd traveling wave belonging to dissimilar, the total length L of circuit can be divided into (0, L/2), (L/2, L) two sections;According to arrive M end the 3rd traveling wave belonging to dissimilar, (0, L/2) section is sub-divided into (0, L/3), (L/3, L/2) two sub-segments, (L/2, L) section is sub-divided into (L/2,2L/3), (2L/3, L) two sub-segments。In sum, the total length L of circuit, it is possible to be divided into 4 sub-segments by L/3, L/2,2L/3 point, respectively 0 < x < L/3, L/3 < x < L/2, L/2 < x < 2L/3,2L/3 < x < L。
Three (four) the individual moment arriving M end are utilized to calculate traveling wave speed and fault distance。
Each section location algorithm is derived as follows:
3.2.10 < x < L/3 section
At 0 < x < L/3 section, fault transient travelling wave is as shown in Figure 2。Simultaneous below equation:
In formula, traveling wave propagates average speed v, moment t occurs fault0, x be unknown number, t1、t2、t3Try to achieve with the wavelet modulus maxima correspondence moment。Observing visible, above simultaneous equations are mutually ratio, it is impossible to try to achieve v, x。From grid chart 2, introduce second class traveling wave due in t4, below equation group can be obtained:
Simultaneous above equation, tries to achieve:
t0=(3t1-t2)/2(11)
Usually, if tpIn (p >=4) moment, the transmitted wave of the 1st echo of opposite end bus arrives bus monitoring side, then exist:
3.2.2L/3 < x < L/2 section
At L/3 < x < L/2 section, fault transient travelling wave is as shown in Figure 3。
At t≤3 τ time interval, list equation:
Simultaneous above equation, solves:
t0=(3t1-t2)/2(15)
3.2.3L/2 < x < 2L/3 section
At L/2 < x < 2L/3 section, fault transient travelling wave is as shown in Figure 4。
At t≤3 τ time interval, list equation:
Simultaneous above equation, solves:
t0=(3t1-t3)/2(17)
3.2.42L/3 < x < L section
At 2L/3 < x < L section, fault transient travelling wave is as shown in Figure 5。
At t≤3 τ time interval, list equation:
Simultaneous above equation, without solving。
Utilize the 1st class traveling wave due in tp, obtain below equation:
Simultaneous above equation, solves:
t0=(3t1-tp)/2(19)
Result according to above 4 sectional analyses, fault distance function x, independent variable is without row wave speed v, only with total track length and t1、t2、t3、t4(tp) etc. relevant, it is possible to be directly used in circuit fault location range finding。
3.3 fault localization flow processs
It is example with " 31 class " bus (for three class buses, opposite end is a class bus to measurement end local terminal) structure, lists fault localization flow process, be illustrated in fig. 4 shown below。Wherein, combinations of polarities link, with the polarity of initial wavefront wavelet coefficient for reference polarity, namely in figure first arrive wave head wavelet coefficient be negative, then second arrival wave head wavelet coefficient be then "+";If recording first to arrive wave head wavelet coefficient for just, then second arrives wave head wavelet coefficient then respectively "-"。
Fault produce in the line before with midpoint after, arriving bus, to measure second traveling wave classification of end different, one of criterion of composition figure below fault location。In the process being subdivided into four fault zones, it is possible to wave head polarity is combined with range finding formula, verifies each other。
Fault localization flow process is as shown in Figure 1。
A kind of based on the single-ended traveling wave fault location method not affected by traveling wave speed degree, first embodiment of the present invention, its step:
A, it is three class buses and opposite end is that a class bus carries out traveling wave record ripple and starts according to measuring end local terminal, detection primary fault current traveling wave and carry out triumphant human relations Bel conversion;
B, to IαCarry out One Dimension Continuous Wavelet Transform, extract t1、t2、t3、t4(tp);
C, corresponding wavelet modulus maxima moment t1、t2、t3、t4(tp), read wavelet coefficient polarity, size, it is determined that wavefront polarity;
D, carrying out combinations of polarities according to principle once: with the polarity of initial wavefront wavelet coefficient for reference polarity, if recording first, to arrive wave head wavelet coefficient be negative, then second arrival wave head wavelet coefficient be then "+";If recording first to arrive wave head wavelet coefficient for just, then second arrives wave head wavelet coefficient then respectively "-";
If e, when record first arrive wave head wavelet coefficient for just with the second arrive wave head wavelet coefficient, respectively during "-", then x < L/2, according to:
Formula one:
Determine whether to meet 0 < x < L/3,
Formula two:
Determine whether to meet L/3 < x < L/2;
F, when record first arrive wave head wavelet coefficient be negative, then second arrival wave head wavelet coefficient be then "+" time, then x > L/2, according to:
Formula one:
Determine whether to meet L/2 < x < 2L/3,
Formula two:
Determine whether to meet 2L/3 < x < L。
One of the present embodiment: B phase earth fault occurs the 220kV Kang Sui A-wire of fault localization example on April 5th, 1: 2002 14:33:07, Heilungkiang Suihua Electric Power Bureau administered total length 64.3km。Wherein on the same bus in Kang Jin side, the faulted phase current transients fault compoment waveform of 3 circuits (containing faulty line interior) is as shown in Figure 6。Faulty line two ends bus is all connected to other circuit a plurality of, is the 3rd class bus。
Fault current traveling wave carries out triumphant human relations Bel conversion, and α mold component electric current is such as shown in Fig. 7 the first row waveform。To α mold component electric current application One Dimension Continuous Wavelet Transform, the db10 small echo adopting vanishing moment exponent number to be 10 decomposes, yardstick j=1~100, the then such as figure shown in Fig. 7 the second row of the wavelet coefficient under each yardstick;Wavelet coefficient during yardstick j=50 is waveform shown in Fig. 7 the third line such as;The maximum of the wavelet coefficient of yardstick 1~100 is waveform shown in Fig. 7 fourth line such as, and modulus maximum position is obvious。Amplifying along sampled point direction, the sampling number reading each wave head corresponding is as follows:
t1=62, t2=258, t3=327, t4=435
Fault is positioned at (L/3, L/2) section, brings formula (16) into, and range measurement is 27.34km, is 0.06km with physical fault distance 27.4km error。
One of the present embodiment: on April 16th, 2002 4:29:39, above-mentioned fault diagnosis example one Kang Sui A-wire occur B phase earth fault, wherein Kang Jin side faulty line current temporary state fault component waveform is as shown in Figure 8。
Fault current traveling wave carries out triumphant human relations Bel conversion, and α mold component electric current is such as shown in Fig. 9 the first row waveform。To α mold component electric current application One Dimension Continuous Wavelet Transform, the db10 small echo adopting vanishing moment exponent number to be 10 decomposes, yardstick j=1~100, then the wavelet coefficient under each yardstick, j=50 time wavelet coefficient figure as shown in 9 second and third row respectively。The maximum of yardstick j=1~100 wavelet coefficient is waveform shown in 9 fourth lines such as, and the sampling number reading each wave head corresponding along sampled point direction is as follows:
t1=70, t2=133, t3=445, t4=454
Fault is positioned at (0, L/3) section, brings formula (14) into, range measurement is 9.062km, being 0.107km with physical fault distance 8.955km error, compared with the position of failure point 9.2km that D type traveling wave range measurement principle obtains, range measurement is more accurate。
Patent of the present invention has lower feature:
1, owing to devising point computational methods for different sections, only relevant for the time of advent with wavelet coefficient polarity, wave head, no longer there is relation with traveling wave spread speed, therefore no longer limited by implementation condition, improve the degree of accuracy of range finding。
2, owing to devising single-ended traveling wave fault location computational methods, no longer traveling wave spread speed is required, extend the use scope of Single Terminal Traveling Wave Fault Location device。
3, owing to devising the restriction that high-tension line length is carried out numerical range, technical characteristic in the technical scheme making numerical range be patent of the present invention, calculate or tested by limited number of time the technical characteristic drawn not by formula, test shows that the technical characteristic of this numerical range achieves good technique effect。
4, owing to devising the technical characteristic of patent of the present invention, effect in the set individually and each other of technical characteristic, by test show, the property indices of patent of the present invention be existing property indices be at least 1.7 times, by assess there is good market value。
Above-described embodiment is a kind of way of realization based on the single-ended traveling wave fault location computational methods not affected by traveling wave speed degree and distance-finding method provided by the present invention; other deformation according to scheme provided by the present invention; increase or reduce composition therein or step; or the present invention is used for other the technical field close with the present invention, belongs to protection scope of the present invention。
Claims (3)
1. one kind based on the single-ended traveling wave fault location computational methods not affected by traveling wave speed, it is characterised in that: the steps include:
During line failure, trouble point send fault traveling wave to M bus, N bus, if fault occurs the moment to be t simultaneously0, the moment of initial current traveling wave arrival M end is t1M, thereafter first or second class traveling wave arrive M end moment be t2M, and first, second or the 3rd class traveling wave arrive M end moment be t3M、t4M..., if L is the total length of circuit, x is the distance that bus M is arrived in trouble point,
When 0 < x < during L/3:
Usually, if tpIn (p >=4) moment, the transmitted wave of the 1st echo of opposite end bus arrives bus monitoring side,
t0=(3t1-t2)/2
When L/3 < x < during L/2:
t0=(3t1-t2)/2
When L/2 < x < during 2L/3:
t0=(3t1-t3)/2
As 2L/3, < x is < during L
t0=(3t1-tp)/2
2. as claimed in claim 1 based on the single-ended traveling wave fault location computational methods not affected by traveling wave speed, it is characterised in that:
The steps include:
3.1 phase-model transformations
There is coupling in three-phase high-voltage, medium-voltage line, in order to eliminate the impact of coupling between three-phase, need first traveling-wave component to be carried out phase-model transformation, it is separate 0, α, β component by the phase component decoupling intercoupled, converting according to triumphant logical sequence Bel (Karrenbauer), three-phase current decoupling is:
Wherein: Iα、Iβ、I0Respectively phase current Ia、Ib、IcTriumphant human relations Bel convert under α mold component, β mold component and zero _exit, wherein zero _exit is propagated between three-phase conductor and the earth, α mold component is propagated between A phase and B phase line, β mold component is propagated between A phase and C phase line, α mold component and β mold component are only propagated between the conductors, so also known as Aerial mode component
Zero _exit is very serious along with frequency raises decay, does not generally select zero _exit as the object of wavelet transformation;And line mould transient state travelling wave signal can ensure that enough sensitivity as measurement signal, it is possible to carry out wavelet analysis as measuring signal, adopt the α mold component of fault current as measuring signal;
3.2 location algorithms
M end gather fault current travelling wave signal, MN circuit diverse location fault, traveling wave process as Figure 2-Figure 5,
When line failure, trouble point send fault traveling wave to M bus, N bus, if fault occurs the moment to be t simultaneously0, the moment of initial current traveling wave arrival M end is t1M, thereafter first or second class traveling wave arrive M end moment be t2M, and first, second or the 3rd class traveling wave arrive M end moment be t3M、t4M……
If L is the total length of circuit, x is the distance that bus M is arrived in trouble point,
According to arrive M end the 2nd traveling wave belonging to dissimilar, the total length L of circuit can be divided into (0, L/2), (L/2, L) two sections;According to arrive M end the 3rd traveling wave belonging to dissimilar, (0, L/2) section is sub-divided into (0, L/3), (L/3, L/2) two sub-segments, (L/2, L) section is sub-divided into (L/2,2L/3), (2L/3, L) two sub-segments, in sum, the total length L of circuit, it is possible to be divided into 4 sub-segments by L/3, L/2,2L/3 point, respectively 0 < x < L/3, L/3 < x < L/2, L/2 < x < 2L/3,2L/3 < x < L
Utilize three (four) the individual moment arriving M end to calculate fault and moment and fault distance occur,
Each section location algorithm is derived as follows:
3.2.10 < x < L/3 section
At 0 < x < L/3 section, fault transient travelling wave is as in figure 2 it is shown, simultaneous below equation:
In formula, traveling wave propagates average speed v, moment t occurs fault0, x be unknown number, t1、t2、t3Trying to achieve with the wavelet modulus maxima correspondence moment, observe visible, above simultaneous equations are mutually ratio, it is impossible to try to achieve v, x, from grid chart 2, introduce second class traveling wave due in t4, below equation group can be obtained:
Simultaneous above equation, tries to achieve:
t0=(3t1-t2)/2
Usually, if tpIn (p >=4) moment, the transmitted wave of the 1st echo of opposite end bus arrives bus monitoring side, then exist:
3.2.2L/3 < x < L/2 section
At L/3 < x < L/2 section, fault transient travelling wave as it is shown on figure 3,
At t≤3 τ time interval, list equation:
Simultaneous above equation, solves:
t0=(3t1-t2)/2
3.2.3L/2 < x < 2L/3 section
At L/2 < x < 2L/3 section, fault transient travelling wave as shown in Figure 4,
At t≤3 τ time interval, list equation:
Simultaneous above equation, solves:
t0=(3t1-t3)/2
3.2.42L/3 < x < L section
At 2L/3 < x < L section, fault transient travelling wave as it is shown in figure 5,
At t≤3 τ time interval, list equation:
Simultaneous above equation, without solving。
Utilize the 1st class traveling wave due in tp, obtain below equation:
Simultaneous above equation, solves:
t0=(3t1-tp)/2
Result according to above 4 sectional analyses, in fault distance function x expression formula, independent variable is without row wave speed v, only with total track length and t1、t2、t3、t4(tp) etc. relevant, it is possible to be directly used in circuit fault location range finding。
3. based on the single-ended traveling wave fault location method not affected by traveling wave speed,
Its step:
A, it is three class buses and opposite end is that a class bus carries out traveling wave record ripple and starts according to measuring end local terminal, detection primary fault current traveling wave and carry out triumphant human relations Bel conversion;
B, to IαCarry out One Dimension Continuous Wavelet Transform, extract t1、t2、t3、t4(tp);
C, corresponding wavelet modulus maxima moment t1、t2、t3、t4(tp), read wavelet coefficient polarity, size, it is determined that wavefront polarity;
D, carry out combinations of polarities according to following principle: with the polarity of initial wavefront wavelet coefficient for reference polarity, if recording first to arrive wave head wavelet coefficient for " ", then second arrive wave head wavelet coefficient be then "+";If recording first to arrive wave head wavelet coefficient for just, then second arrives wave head wavelet coefficient then respectively " ";
If e records first arrival wave head wavelet coefficient, second arrival wave head wavelet coefficient opposite polarity, then x < L/2, according to:
Formula one:
Determine whether to meet 0 < x < L/3,
Formula two:
Determine whether to meet L/3 < x < L/2;
F, when record first arrive wave head wavelet coefficient, second arrive wave head wavelet coefficient polarity identical time, then x > L/2, according to:
Formula one:
Determine whether to meet L/2 < x < 2L/3,
Formula two:
Determine whether to meet 2L/3 < x < L。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610209799.4A CN105699855B (en) | 2016-04-06 | 2016-04-06 | Based on the single-ended traveling wave fault location calculation method not influenced by traveling wave speed and distance measuring method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610209799.4A CN105699855B (en) | 2016-04-06 | 2016-04-06 | Based on the single-ended traveling wave fault location calculation method not influenced by traveling wave speed and distance measuring method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105699855A true CN105699855A (en) | 2016-06-22 |
CN105699855B CN105699855B (en) | 2018-11-20 |
Family
ID=56218360
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610209799.4A Expired - Fee Related CN105699855B (en) | 2016-04-06 | 2016-04-06 | Based on the single-ended traveling wave fault location calculation method not influenced by traveling wave speed and distance measuring method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105699855B (en) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106093706A (en) * | 2016-07-05 | 2016-11-09 | 四川中光防雷科技股份有限公司 | A kind of transmission line malfunction interval determines method and system |
CN107037316A (en) * | 2016-11-23 | 2017-08-11 | 国网湖北省电力公司检修公司 | A kind of single ended waveform automatic identifying method for being adapted to transmission line travelling wave ranging |
CN108051705A (en) * | 2017-12-26 | 2018-05-18 | 华北电力大学 | The method for locating single-phase ground fault and system of power distribution network |
CN108120898A (en) * | 2017-12-18 | 2018-06-05 | 广东电网有限责任公司佛山供电局 | Low sampling rate traveling wave fault positioning method based on VMD and SDEO |
CN108169626A (en) * | 2018-01-16 | 2018-06-15 | 国网冀北电力有限公司秦皇岛供电公司 | A kind of T-shaped power transmission line fault locating method and system |
CN109188210A (en) * | 2018-10-23 | 2019-01-11 | 国网江苏省电力有限公司徐州供电分公司 | A kind of urban electric power cable Two-terminal Fault Location method based on VMD-Hilbert transformation |
CN109284830A (en) * | 2018-11-27 | 2019-01-29 | 广东电网有限责任公司 | A kind of combination single-ended traveling wave fault location algorithm |
CN109901016A (en) * | 2019-03-13 | 2019-06-18 | 清华四川能源互联网研究院 | A kind of both-end distance measuring method eliminating velocity of wave and influencing |
CN110927510A (en) * | 2019-10-17 | 2020-03-27 | 东北大学 | Frequency domain method for power transmission line double-end traveling wave fault location |
CN111433616A (en) * | 2017-11-17 | 2020-07-17 | Abb电网瑞士股份公司 | Parametric traveling wave based fault location for power transmission lines |
CN111766470A (en) * | 2020-06-24 | 2020-10-13 | 湖南大学 | Fault positioning method and system for high-voltage direct-current transmission line and direct-current transmission line |
CN114113879A (en) * | 2021-11-04 | 2022-03-01 | 云南电网有限责任公司 | Low-voltage line fault positioning method and system based on wavelet transformation |
CN114184884A (en) * | 2021-11-23 | 2022-03-15 | 昆明理工大学 | Automatic construction method for power grid fault traveling wave distance measurement equation |
CN117233529A (en) * | 2023-09-20 | 2023-12-15 | 青岛鼎信通讯科技有限公司 | Single-ended distance measurement method suitable for inter-phase short-circuit fault of power distribution network |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101050550B1 (en) * | 2009-12-29 | 2011-07-19 | 주식회사 효성 | Fault point detection method of transmission line using secondary difference of traveling wave signal |
CN102253310A (en) * | 2011-05-18 | 2011-11-23 | 昆明理工大学 | Method for identifying property of second wave head during fault location of alternating-current power transmission line |
CN102288869A (en) * | 2011-05-10 | 2011-12-21 | 山东大学 | Single-end traveling wave fault ranging method for power transmission line |
CN103278747A (en) * | 2013-06-03 | 2013-09-04 | 东南大学 | High-tension transmission line single-ended traveling wave fault distance detection method combined with time-frequency characteristics |
CN104793102A (en) * | 2015-04-08 | 2015-07-22 | 三峡大学 | Single-ended traveling-wave fault distance measurement method |
-
2016
- 2016-04-06 CN CN201610209799.4A patent/CN105699855B/en not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101050550B1 (en) * | 2009-12-29 | 2011-07-19 | 주식회사 효성 | Fault point detection method of transmission line using secondary difference of traveling wave signal |
CN102288869A (en) * | 2011-05-10 | 2011-12-21 | 山东大学 | Single-end traveling wave fault ranging method for power transmission line |
CN102253310A (en) * | 2011-05-18 | 2011-11-23 | 昆明理工大学 | Method for identifying property of second wave head during fault location of alternating-current power transmission line |
CN103278747A (en) * | 2013-06-03 | 2013-09-04 | 东南大学 | High-tension transmission line single-ended traveling wave fault distance detection method combined with time-frequency characteristics |
CN104793102A (en) * | 2015-04-08 | 2015-07-22 | 三峡大学 | Single-ended traveling-wave fault distance measurement method |
Non-Patent Citations (2)
Title |
---|
杨关春,彭丽: "一种不受波速影响的单端行波故障测距新算法", 《山东理工大学学报》 * |
马丹丹: "输电线路故障暂态信号分析和行波测距研究", 《中国优秀硕士论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106093706B (en) * | 2016-07-05 | 2018-11-09 | 四川中光防雷科技股份有限公司 | A kind of transmission line malfunction section determines method and system |
CN106093706A (en) * | 2016-07-05 | 2016-11-09 | 四川中光防雷科技股份有限公司 | A kind of transmission line malfunction interval determines method and system |
CN107037316A (en) * | 2016-11-23 | 2017-08-11 | 国网湖北省电力公司检修公司 | A kind of single ended waveform automatic identifying method for being adapted to transmission line travelling wave ranging |
US12044718B2 (en) | 2017-11-17 | 2024-07-23 | Hitachi Energy Ltd | Parameter free traveling wave based fault location for power transmission lines |
CN111433616A (en) * | 2017-11-17 | 2020-07-17 | Abb电网瑞士股份公司 | Parametric traveling wave based fault location for power transmission lines |
CN108120898B (en) * | 2017-12-18 | 2019-07-05 | 广东电网有限责任公司佛山供电局 | Low sampling rate traveling wave fault positioning method based on VMD and SDEO |
CN108120898A (en) * | 2017-12-18 | 2018-06-05 | 广东电网有限责任公司佛山供电局 | Low sampling rate traveling wave fault positioning method based on VMD and SDEO |
CN108051705A (en) * | 2017-12-26 | 2018-05-18 | 华北电力大学 | The method for locating single-phase ground fault and system of power distribution network |
CN108169626B (en) * | 2018-01-16 | 2021-05-14 | 国网冀北电力有限公司秦皇岛供电公司 | T-shaped power transmission line fault positioning method and system |
CN108169626A (en) * | 2018-01-16 | 2018-06-15 | 国网冀北电力有限公司秦皇岛供电公司 | A kind of T-shaped power transmission line fault locating method and system |
CN109188210A (en) * | 2018-10-23 | 2019-01-11 | 国网江苏省电力有限公司徐州供电分公司 | A kind of urban electric power cable Two-terminal Fault Location method based on VMD-Hilbert transformation |
CN109284830A (en) * | 2018-11-27 | 2019-01-29 | 广东电网有限责任公司 | A kind of combination single-ended traveling wave fault location algorithm |
CN109901016A (en) * | 2019-03-13 | 2019-06-18 | 清华四川能源互联网研究院 | A kind of both-end distance measuring method eliminating velocity of wave and influencing |
CN109901016B (en) * | 2019-03-13 | 2020-11-20 | 清华四川能源互联网研究院 | Double-end ranging method for eliminating wave velocity influence |
CN110927510A (en) * | 2019-10-17 | 2020-03-27 | 东北大学 | Frequency domain method for power transmission line double-end traveling wave fault location |
CN110927510B (en) * | 2019-10-17 | 2021-10-01 | 东北大学 | Frequency domain method for power transmission line double-end traveling wave fault location |
CN111766470A (en) * | 2020-06-24 | 2020-10-13 | 湖南大学 | Fault positioning method and system for high-voltage direct-current transmission line and direct-current transmission line |
CN114113879A (en) * | 2021-11-04 | 2022-03-01 | 云南电网有限责任公司 | Low-voltage line fault positioning method and system based on wavelet transformation |
CN114113879B (en) * | 2021-11-04 | 2024-04-12 | 云南电网有限责任公司 | Low-voltage line fault positioning method and system based on wavelet transformation |
CN114184884A (en) * | 2021-11-23 | 2022-03-15 | 昆明理工大学 | Automatic construction method for power grid fault traveling wave distance measurement equation |
CN114184884B (en) * | 2021-11-23 | 2024-05-24 | 昆明理工大学 | Automatic construction method of power grid fault traveling wave ranging equation |
CN117233529A (en) * | 2023-09-20 | 2023-12-15 | 青岛鼎信通讯科技有限公司 | Single-ended distance measurement method suitable for inter-phase short-circuit fault of power distribution network |
CN117233529B (en) * | 2023-09-20 | 2024-05-14 | 青岛鼎信通讯科技有限公司 | Single-ended distance measurement method suitable for inter-phase short-circuit fault of power distribution network |
Also Published As
Publication number | Publication date |
---|---|
CN105699855B (en) | 2018-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105699855A (en) | Single-ended traveling fault location calculation method and location method insusceptible to traveling wave speed | |
CN101299538B (en) | Cable-aerial mixed line fault travelling wave ranging method | |
CN102495336B (en) | Distributed single-phase earth fault ranging system and ranging method thereof | |
WO2019080238A1 (en) | Time-domain calculation method for voltage traveling wave differential protection of flexible direct-current transmission line | |
CN105223466B (en) | It is a kind of using modulus maximum than extra high voltage direct current transmission line method of single end distance measurement | |
CN102967779B (en) | Identifying method of distribution parameters of transmission line | |
CN105425107A (en) | Method and system for active power distribution network fault diagnosis and location | |
CN107390081A (en) | It is a kind of to be used for the device and method being monitored powered to deformation of transformer winding | |
CN103913676B (en) | Based on the transmitting line one-end fault localization method of window during variable row ripple identification | |
CN207181589U (en) | It is a kind of to be used for the device being monitored powered to deformation of transformer winding | |
EP3710842B1 (en) | Traveling wave based fault location using unsynchronized measurements for transmission lines | |
Dashtdar et al. | Fault location in the transmission network using a discrete wavelet transform | |
CN103954894A (en) | Partial discharge locating method for three-phase crossed and interconnected cables | |
CN111766470B (en) | Fault positioning method and system for high-voltage direct-current transmission line and direct-current transmission line | |
CN109387744A (en) | Distribution line fault point positioning method and device based on singular value decomposition | |
CN102545177A (en) | Bergeron-model-based simulation-after-test method for fault phase selection of alternating current transmission line | |
CN103439636A (en) | Cable open-circuit point position detecting method | |
US11038342B2 (en) | Traveling wave identification using distortions for electric power system protection | |
CN113253052A (en) | High-voltage direct-current transmission line fault distance measurement method based on improved SMMG | |
CN115902530A (en) | Earth electrode line fault distance measurement method and system | |
CN104655981A (en) | Comparison type directional protection method for travelling wave polarity by utilizing crossing sequence differential transformation | |
CN106841918A (en) | A kind of cable low resistance faults localization method combined using single both-end | |
CN112886547A (en) | Traveling wave protection method and system for high-voltage transmission line based on traveling wave waveform distance | |
CN101825676B (en) | Fault ranging method by utilizing red and green colors to express traveling wave of DC power transmission circuit malfunction voltage | |
CN117421976A (en) | Wind power plant current collection line fault positioning monitoring method and system based on ensemble learning |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181120 Termination date: 20190406 |