CN105699855B - Based on the single-ended traveling wave fault location calculation method not influenced by traveling wave speed and distance measuring method - Google Patents
Based on the single-ended traveling wave fault location calculation method not influenced by traveling wave speed and distance measuring method Download PDFInfo
- Publication number
- CN105699855B CN105699855B CN201610209799.4A CN201610209799A CN105699855B CN 105699855 B CN105699855 B CN 105699855B CN 201610209799 A CN201610209799 A CN 201610209799A CN 105699855 B CN105699855 B CN 105699855B
- Authority
- CN
- China
- Prior art keywords
- traveling wave
- wave
- fault
- section
- bus
- 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.)
- Expired - Fee Related
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
It is a kind of based on the single-ended traveling wave fault localization method not influenced by traveling wave speed degree, it is a kind of based on the single-ended traveling wave fault location calculation method not influenced by traveling wave speed degree；If it is t that the moment, which occurs, for failure0, initial current traveling wave is t at the time of reaching the end M1M, thereafter first or second class traveling wave reach the end M at the time of be t2MAnd first, second or third class traveling wave reach the end M at the time of be t3M、t4M..., if L is the overall length of route, x is distance of the fault point to bus M, when 0<x<When L/3, work as L/2<x<When 2L/3 and work as 2L/3<x<When L, three kinds of states are calculated, only relevant with the time, are no longer had relationship with traveling wave spread speed, therefore no longer limited by implementation condition, are improved the accuracy of ranging.
Description
Technical field
The invention patent relates to a kind of single-ended traveling wave fault location calculation methods, especially a kind of based on not by traveling wave speed
The single-ended traveling wave fault location calculation method and distance measuring method of influence.
Background technique
After power circuit breaks down, accurately determines the fault point of power circuit, be conducive to rapidly removing faults, it is extensive
Multiple 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
In fault localization calculation method, calculate fault distance using the traveling wave spread speed of estimation causes to survey due to the deviation of wave velocity
Away from accuracy it is low.
Summary of the invention
Object of the invention is a kind of based on the single-ended traveling wave fault location calculation method not influenced by traveling wave speed degree；
Object of the invention is a kind of based on the single-ended traveling wave fault location method not influenced by traveling wave speed degree.
In order to overcome the technical drawbacks described above, the purpose of the invention patent be to provide influenced based on not by traveling wave speed it is single-ended
Traveling wave fault location calculation method and distance measuring method, therefore no longer limited by implementation condition, improve the accuracy of ranging.
In order to achieve the above objectives, the invention patent adopts the technical scheme that：
It is a kind of based on the single-ended traveling wave fault location calculation method not influenced by traveling wave speed；Its step：Event occurs for route
When barrier, fault traveling wave is issued to M bus, N bus simultaneously from fault point.If it is t that the moment, which occurs, for failure0, initial current traveling wave arrives
It is t at the time of up to the end M1M, thereafter first or second class traveling wave reach the end M at the time of be t2MAnd first, second or
Three classes traveling wave is t at the time of reaching the end M3M、t4M..., if L is the overall length of route, x is distance of the fault point to bus M,
When 0<x<When L/3：
Generally, if tpThe transmitted wave at (p >=4) moment, the 1st back wave of opposite end bus reaches bus monitoring side,
t0=(3t1-t2)/2
Work as L/3<x<When L/2：
t0=(3t1-t2)/2
Work as L/2<x<When 2L/3：
t0=(3t1-t3)/2
Work as 2L/3<x<When L
t0=(3t1-tp)/2
Due to devising a point calculation method for different sections, pass relevant with the traveling wave surge arrival time of monitoring side
, no longer there is relationship in system with traveling wave spread speed, therefore no longer limited by implementation condition, improves the accuracy of ranging.
The present invention is devised based on the single-ended traveling wave fault location calculation method not influenced by traveling wave speed, step：
3.1 phase-model transformation
There are couplings to need first to eliminate the influence coupled between three-phase to traveling wave for three-phase high-voltage, medium-voltage line
Component carries out phase-model transformation, the phase component to intercouple is decoupled as mutually independent 0, α, β component, according to triumphant logical sequence Bell
(Karrenbauer) it converts, three-phase current decoupling is:
Wherein:Iα、Iβ、I0Respectively phase current Ia、Ib、Icα mold component, β mold component and zero mould under triumphant human relations Bell transformation
Component, wherein zero _exit is propagated between three-phase conductor and the earth, and α mold component is propagated between A phase and B phase line, β mould point
Amount 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 with frequency raising decaying, does not generally select zero _exit as the object of wavelet transformation；And
Line mould transient state travelling wave signal can guarantee enough sensitivity as measuring signal, can be used as measuring signal to carry out small wavelength-division
Analysis.Using the α mold component of fault current as measuring signal.
3.2 location algorithm
Fault current travelling wave signal is acquired at the end M.MN route different location failure, traveling wave process such as Fig. 2-Fig. 5 institute
Show.
When line failure, fault traveling wave is issued to M bus, N bus simultaneously from fault point.If the moment occurs for failure
For t0, initial current traveling wave is t at the time of reaching the end M1M, thereafter first or second class traveling wave reach the end M at the time of be t2M, with
And first, second or third class traveling wave reach the end M at the time of be t3M、t4M……
If L is the overall length of route, x is distance of the fault point to bus M.
According to the 2nd affiliated different type of traveling wave for reaching the end M, the overall length L of route can be divided into (0, L/2), (L/2, L)
Two sections；According to reach the end M the 3rd affiliated different type of traveling wave, (0, L/2) section be sub-divided into (0, L/3), (L/3,
L/2) two sub-segments, (L/2, L) section are sub-divided into (L/2,2L/3), (2L/3, L) two sub-segments.In conclusion route
Overall length L, can 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) for reaching the end M using wavefront calculate traveling wave speed and fault distance at a moment.
Each section location algorithm derives as follows.
3.2.1 0<x<L/3 section
0<x<L/3 section, fault transient travelling wave is as shown in Fig. 2, simultaneous following equation：
In formula, traveling wave propagates average speed v, moment t occurs for failure0, x be unknown number, t1、t2、t3With WAVELET TRANSFORM MODULUS pole
Big value correspond to the moment and acquires, and observes as it can be seen that the above simultaneous equations are mutually proportional, v, x can not be acquired, by grid chart 2 as it can be seen that introducing
Second class traveling wave arrival time t4, following equation group can be obtained：
Simultaneous above equation, is acquired：
t0=(3t1-t2)/2
Generally, if tpThe transmitted wave at (p >=4) moment, the 1st back wave of opposite end bus reaches bus monitoring side, then
In the presence of：
3.2.2 L/3<x<L/2 section
In L/3<x<L/2 section, fault transient travelling wave as shown in figure 3,
In the τ time interval of t≤3, equation is listed：
Simultaneous above equation, solves：
t0=(3t1-t2)/2
3.2.3 L/2<x<2L/3 section
In L/2<x<2L/3 section, fault transient travelling wave as shown in figure 4,
In the τ time interval of t≤3, equation is listed：
Simultaneous above equation, solves：
t0=(3t1-t3)/2
3.2.4 2L/3<x<L section
In 2L/3<x<L section, fault transient travelling wave as shown in figure 5,
In the τ time interval of t≤3, equation is listed：
Simultaneous above equation, no solution.
Utilize the 1st class traveling wave arrival time tp, obtain following equation：
Simultaneous above equation, solves：
t0=(3t1-tp)/2
According to above 4 sectional analyses as a result, in fault distance function x expression formula, independent variable is free of row wave speed
V, only with total track length and t1、t2、t3、t4(tp) etc. it is related, be used directly for the fault location ranging of route.
The present invention devises：It is a kind of based on the single-ended traveling wave fault location method not influenced by traveling wave speed,
Its step：
A, according to measurement end local terminal be three classes bus and opposite end is that a kind of bus carries out the starting of traveling wave recording, detect it is initial therefore
Hinder current traveling wave and carries out triumphant human relations Bell transformation；
B, to IαOne Dimension Continuous Wavelet Transform is carried out, t is extracted1、t2、t3、t4(tp)；
C, corresponding wavelet modulus maxima moment t1、t2、t3、t4(tp), wavelet coefficient polarity, size are read, determines row
Wave wave head polarity；
D, combinations of polarities is carried out according to following principle：Using the polarity of initial wavefront wavelet coefficient as reference polarity, if
It measures first arrival wave head wavelet coefficient to be negative, then second arrival wave head wavelet coefficient is then "+"；If measuring first to arrive
It is positive up to wave head wavelet coefficient, then second arrival wave head wavelet coefficient is then respectively "-"；
E, when if measuring first arrival wave head wavelet coefficient to be positive with second arrival wave head wavelet coefficient is "-",
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 measuring that first arrival wave head wavelet coefficient is negative, second arrival wave head wavelet coefficient is "+", 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 scheme, the sectional based on triumphant logical sequence Bell transformation carries out being calculated as important technology feature, is being based on
The single-ended traveling wave fault positioning not influenced by traveling wave speed degree and distance measuring method, have novelty, creativeness and practicability, at this
Term in technical solution is all that can be explained and be understood with patent document in the art.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with
It obtains other drawings based on these drawings.
Fig. 1 is the flow chart based on the single-ended traveling wave fault location method not influenced by traveling wave speed degree of invention；
Fig. 2 is invention 0<x<L/3 section, fault transient travelling wave figure；
Fig. 3 is invention in L/3<x<L/2 section, fault transient travelling wave figure；
Fig. 4 is invention in L/2<x<2L/3 section, fault transient travelling wave figure；
Fig. 5 is invention in 2L/3<x<L section, fault transient travelling wave figure；
Fig. 6 is the faulted phase current of 3 routes (containing including faulty line) on the same bus in the side Kang Jin in inventive embodiments
Transients fault compoment waveform diagram；
Fig. 7 is the wavelet transformation curve graph of invention figure 6；
Fig. 8 is that B phase ground fault occurs for Kang Sui A-wire in invention example, wherein the side Kang Jin faulty line current temporary state failure
Component waveform figure；
Fig. 9 is the wavelet transformation curve graph of invention figure 8.
Specific embodiment
Below with reference to embodiment, the present invention is further described, following embodiment is intended to illustrate invention rather than to this
Invention further limits.
It is a kind of based on the single-ended traveling wave fault location calculation method not influenced by traveling wave speed degree, step：
Single-ended traveling wave fault location algorithm
3.1 phase-model transformation
There are couplings to need first to eliminate the influence coupled between three-phase to traveling wave for three-phase high-voltage, medium-voltage line
Component carries out phase-model transformation, and the phase component to intercouple is decoupled as mutually independent 0, α, β component.According to triumphant logical sequence Bell
(Karrenbauer) it converts, three-phase current decoupling is:
Wherein:Iα、Iβ、I0Respectively phase current Ia、Ib、Icα mold component, β mold component and zero mould under triumphant human relations Bell transformation
Component.Wherein zero _exit is propagated between three-phase conductor and the earth, and α mold component is propagated between A phase and B phase line, β mould point
Amount 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 with frequency raising decaying, does not generally select zero _exit as the object of wavelet transformation；And
Line mould transient state travelling wave signal can guarantee enough sensitivity as measuring signal, can be used as measuring signal to carry out small wavelength-division
Analysis.Using the α mold component of fault current as measuring signal.
3.2 location algorithm
Fault current travelling wave signal is acquired at the end M.MN route different location failure, traveling wave process such as Fig. 2-Fig. 5 institute
Show.
By Fig. 2-Fig. 5 as it can be seen that when line failure, fault traveling wave is issued to M bus, N bus simultaneously from fault point.
If it is t that the moment, which occurs, for failure0, initial current traveling wave is t at the time of reaching the end M1M, thereafter first or second class traveling wave reach M
It is t at the time of end2MAnd first, second or third class traveling wave reach the end M at the time of be t3M、t4M……
If L is the overall length of route, x is distance of the fault point to bus M.
According to the 2nd affiliated different type of traveling wave for reaching the end M, the overall length L of route can be divided into (0, L/2), (L/2, L)
Two sections；According to reach the end M the 3rd affiliated different type of traveling wave, (0, L/2) section be sub-divided into (0, L/3), (L/3,
L/2) two sub-segments, (L/2, L) section are sub-divided into (L/2,2L/3), (2L/3, L) two sub-segments.In conclusion route
Overall length L, can 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.
Traveling wave speed and fault distance are calculated at a moment using three (four) for reaching the end M.
Each section location algorithm derives as follows：
3.2.1 0<x<L/3 section
0<x<L/3 section, fault transient travelling wave are as shown in Figure 2.Simultaneous following equation：
In formula, traveling wave propagates average speed v, moment t occurs for failure0, x be unknown number, t1、t2、t3With WAVELET TRANSFORM MODULUS pole
Big value corresponds to the moment and acquires.Observation can not acquire v, x as it can be seen that the above simultaneous equations are mutually proportional.By grid chart 2 as it can be seen that introducing
Second class traveling wave arrival time t4, following equation group can be obtained：
Simultaneous above equation, is acquired：
t0=(3t1-t2)/2 (11)
Generally, if tpThe transmitted wave at (p >=4) moment, the 1st back wave of opposite end bus reaches bus monitoring side, then
In the presence of：
3.2.2 L/3<x<L/2 section
In L/3<x<L/2 section, fault transient travelling wave are as shown in Figure 3.
In the τ time interval of t≤3, equation is listed：
Simultaneous above equation, solves：
t0=(3t1-t2)/2 (15)
3.2.3 L/2<x<2L/3 section
In L/2<x<2L/3 section, fault transient travelling wave are as shown in Figure 4.
In the τ time interval of t≤3, equation is listed：
Simultaneous above equation, solves：
t0=(3t1-t3)/2 (17)
3.2.4 2L/3<x<L section
In 2L/3<x<L section, fault transient travelling wave are as shown in Figure 5.
In the τ time interval of t≤3, equation is listed：
Simultaneous above equation, no solution.
Utilize the 1st class traveling wave arrival time tp, obtain following equation：
Simultaneous above equation, solves：
t0=(3t1-tp)/2 (19)
According to above 4 sectional analyses as a result, fault distance function x, independent variable is free of row wave speed v, only with line
Road overall length and t1、t2、t3、t4(tp) etc. it is related, be used directly for the fault location ranging of route.
3.3 fault localization processes
By taking " three is a kind of " bus (for measurement end local terminal as three classes bus, opposite end is a kind of bus) structure as an example, event is listed
Hinder ranging process, is illustrated in fig. 4 shown below.Wherein, combinations of polarities link is with reference to pole with the polarity of initial wavefront wavelet coefficient
First arrival wave head wavelet coefficient is negative in property namely figure, then second arrival wave head wavelet coefficient is then "+"；If measuring
First arrival wave head wavelet coefficient is positive, then second arrival wave head wavelet coefficient is then respectively "-".
Failure generates put in the line before it is different with second traveling wave classification after midpoint, reaching bus measurement end,
Constitute one of the criterion of following figure fault location.It, can be by wave head polarity and ranging during being subdivided into four fault zones
Formula combines, and verifies each other.
Fault localization process is as shown in Figure 1.
It is a kind of based on the single-ended traveling wave fault location method not influenced by traveling wave speed degree, first implementation of the invention
Example, step：
A, according to measurement end local terminal be three classes bus and opposite end is that a kind of bus carries out the starting of traveling wave recording, detect it is initial therefore
Hinder current traveling wave and carries out triumphant human relations Bell transformation；
B, to IαOne Dimension Continuous Wavelet Transform is carried out, t is extracted1、t2、t3、t4(tp)；
C, corresponding wavelet modulus maxima moment t1、t2、t3、t4(tp), wavelet coefficient polarity, size are read, determines row
Wave wave head polarity；
D, combinations of polarities is carried out according to principle once：Using the polarity of initial wavefront wavelet coefficient as reference polarity, if
It measures first arrival wave head wavelet coefficient to be negative, then second arrival wave head wavelet coefficient is then "+"；If measuring first to arrive
It is positive up to wave head wavelet coefficient, then second arrival wave head wavelet coefficient is then respectively "-"；
E, when being respectively if measuring first arrival wave head wavelet coefficient and being positive with second arrival wave head wavelet coefficient
When "-", 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, it is negative when measuring first arrival wave head wavelet coefficient, then when second arrival wave head wavelet coefficient is then "+",
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：Fault localization example one：5 days 14 April in 2002:33:07, Heilungkiang Suihua Electric Power Bureau is managed
B phase ground fault occurs for the 220kV Kang Sui A-wire for having jurisdiction over overall length 64.3km.Wherein 3 routes (contain failure on the same bus in the side Kang Jin
Including route) faulted phase current transients fault compoment waveform it is as shown in Figure 6.Faulty line both ends bus be all connected with it is a plurality of its
Its route is the 3rd class bus.
Triumphant human relations Bell transformation is carried out to fault current traveling wave, α mold component electric current is as shown in Fig. 7 the first row waveform.To α mould point
Electric current application One Dimension Continuous Wavelet Transform is measured, vanishing moment order is used to be decomposed for 10 db10 small echo, scale j=1~
100, then the figure as shown in the second row of Fig. 7 of the wavelet coefficient under each scale；Wavelet coefficient such as Fig. 7 third when scale j=50
Waveform shown in row；The maximum value of the wavelet coefficient of scale 1~100 waveform as shown in Fig. 7 fourth line, modulus maximum position are obvious.
Amplify along sampled point direction, it is as follows to read the corresponding sampling number of each wave head:
t1=62, t2=258, t3=327, t4=435
Failure is located at (L/3, L/2) section, brings into formula (16), distance measurement result 27.34km, with physical fault distance
27.4km error is 0.06km.
One of the present embodiment：16 days 4 April in 2002:29:39, above-mentioned one Kang Sui A-wire of fault diagnosis example occurs B phase and is grounded event
Barrier, wherein the side Kang Jin faulty line current temporary state fault component waveform is as shown in Figure 8.
Triumphant human relations Bell transformation is carried out to fault current traveling wave, α mold component electric current is as shown in Fig. 9 the first row waveform.To α mould point
Electric current application One Dimension Continuous Wavelet Transform is measured, vanishing moment order is used to be decomposed for 10 db10 small echo, scale j=1~
100, then the figure as shown in 9 second and third row respectively of the wavelet coefficient when wavelet coefficient under each scale, j=50.Scale j=1
The maximum value of~100 wavelet coefficients waveform as shown in 9 fourth lines reads the corresponding sampled point of each wave head along sampled point direction
Number is as follows:
t1=70, t2=133, t3=445, t4=454
Failure is located at (0, L/3) section, brings into formula (14), distance measurement result 9.062km, with physical fault distance
8.955km error is 0.107km, and compared with the position of failure point 9.2km that D type traveling wave range measurement principle obtains, distance measurement result is more
Accurately.
The invention patent has the characteristics that lower：
1, only relevant with wavelet coefficient polarity, wave head arrival time due to devising a point calculation method for different sections
Relationship no longer has relationship with traveling wave spread speed, therefore is no longer limited by implementation condition, improves the accuracy of ranging.
2, due to devising single-ended traveling wave fault location calculation method, no longer traveling wave spread speed is required, is extended
The use scope of Single Terminal Traveling Wave Fault Location device.
3, the restriction for carrying out numberical range to high-tension line length due to devising, makes numberical range the invention patent
Technical characteristic in technical solution is not the technical characteristic for being calculated by formula or being obtained by limited times test, experiments have shown that
The technical characteristic of the numberical range achieves good technical effect.
4, due to devising the technical characteristic of the invention patent, in the work of the independent and mutual set of technical characteristic
With, by experiments have shown that, the performance indexes of the invention patent is that existing performance indexes is at least 1.7 times, is led to
Crossing assessment has good market value.
Above-described embodiment is provided by the present invention based on the single-ended traveling wave fault location not influenced by traveling wave speed degree
A kind of way of realization of calculation method and distance measuring method, other deformations of provided scheme, increase or subtract according to the present invention
Few composition therein or step, or the present invention is used for other technical fields close with the present invention, belong to the present invention
Protection scope.
Claims (3)
1. a kind of based on the single-ended traveling wave fault location calculation method not influenced by traveling wave speed, it is characterised in that：Its step is：
When line failure, fault traveling wave is issued to M bus, N bus simultaneously from fault point, if it is t that the moment, which occurs, for failure0, just
Beginning current traveling wave is t at the time of reaching the end M1M, thereafter first or second class traveling wave reach the end M at the time of be t2M, Yi Ji
One, second or third class traveling wave reach the end M at the time of be t3M、t4M..., if L is the overall length of route, x is fault point to mother
The distance of line M,
When 0<x<When L/3：
Generally, if tpThe transmitted wave at (p >=4) moment, the 1st back wave of opposite end bus reaches bus monitoring side,
t0=(3t1-t2)/2
Work as L/3<x<When L/2：
t0=(3t1-t2)/2
Work as L/2<x<When 2L/3：
t0=(3t1-t3)/2
Work as 2L/3<x<When L
t0=(3t1-tp)/2
2. as described in claim 1 based on the single-ended traveling wave fault location calculation method not influenced by traveling wave speed, feature
It is：
Its step is：
3.1 phase-model transformation
There are couplings to need first to eliminate the influence coupled between three-phase to traveling-wave component for three-phase high-voltage, medium-voltage line
Phase-model transformation is carried out, the phase component to intercouple is decoupled as mutually independent 0, α, β component, according to triumphant logical sequence Bell
(Karrenbauer) it converts, three-phase current decoupling is:
Wherein:Iα、Iβ、I0Respectively phase current Ia、Ib、Icα mold component, β mold component and zero _exit under triumphant human relations Bell transformation,
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 in A
Mutually to be propagated between 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 with frequency raising decaying, does not generally select zero _exit as the object of wavelet transformation；And line mould
Transient state travelling wave signal can guarantee enough sensitivity as measuring signal, can be used as measuring signal to carry out wavelet analysis,
Using the α mold component of fault current as measuring signal；
3.2 location algorithm
Fault current travelling wave signal is acquired at the end M, MN route different location failure, traveling wave process is as Figure 2-Figure 5,
When line failure, fault traveling wave is issued to M bus, N bus simultaneously from fault point, if it is t that the moment, which occurs, for failure0,
Initial current traveling wave is t at the time of reaching the end M1M, thereafter first or second class traveling wave reach the end M at the time of be t2M, Yi Ji
One, second or third class traveling wave reach the end M at the time of be t3M、t4M……
If L is the overall length of route, x is distance of the fault point to bus M,
According to the 2nd affiliated different type of traveling wave for reaching the end M, the overall length L of route can be divided into (0, L/2), (L/2, L) two
Section；According to the 3rd affiliated different type of traveling wave for reaching the end M, (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 conclusion route is complete
Long L, can be divided into 4 sub-segments by L/3, L/2,2L/3 point, and respectively 0<x<L/3, L/3<x<L/2, L/2<x<2L/3,2L/3
<x<L,
Failure is calculated a moment using three (four) for reaching the end M, and moment and fault distance occurs,
Each section location algorithm derives as follows：
3.2.1 0<x<L/3 section
0<x<L/3 section, fault transient travelling wave is as shown in Fig. 2, simultaneous following equation：
In formula, traveling wave propagates average speed v, moment t occurs for failure0, x be unknown number, t1、t2、t3With wavelet modulus maxima
The corresponding moment is acquired, and observation can not acquire v, x, by grid chart 2 as it can be seen that introducing second as it can be seen that the above simultaneous equations are mutually proportional
Class traveling wave arrival time t4, following equation group can be obtained：
Simultaneous above equation, is acquired：
t0=(3t1-t2)/2
Generally, if tpThe transmitted wave at (p >=4) moment, the 1st back wave of opposite end bus reaches bus monitoring side, then exists：
3.2.2 L/3<x<L/2 section
In L/3<x<L/2 section, fault transient travelling wave as shown in figure 3,
In the τ time interval of t≤3, equation is listed：
Simultaneous above equation, solves：
t0=(3t1-t2)/2
3.2.3 L/2<x<2L/3 section
In L/2<x<2L/3 section, fault transient travelling wave as shown in figure 4,
In the τ time interval of t≤3, equation is listed：
Simultaneous above equation, solves：
t0=(3t1-t3)/2
3.2.4 2L/3<x<L section
In 2L/3<x<L section, fault transient travelling wave as shown in figure 5,
In the τ time interval of t≤3, equation is listed：
Simultaneous above equation, no solution；
Utilize the 1st class traveling wave arrival time tp, obtain following equation：
Simultaneous above equation, solves：
t0=(3t1-tp)/2
According to above 4 sectional analyses as a result, in fault distance function x expression formula, independent variable is free of row wave speed v, only
Only with total track length and t1、t2、t3、t4(tp) etc. it is related, be used directly for the fault location ranging of route.
3. it is a kind of based on the single-ended traveling wave fault location method not influenced by traveling wave speed,
Its step：
A, according to measurement end local terminal be three classes bus and opposite end is that a kind of bus carries out the starting of traveling wave recording, detection primary fault electricity
Popular wave and carry out triumphant human relations Bell transformation；
B, to IαOne Dimension Continuous Wavelet Transform is carried out, t is extracted1、t2、t3、t4(tp)；
C, corresponding wavelet modulus maxima moment t1、t2、t3、t4(tp), wavelet coefficient polarity, size are read, determines traveling wave wave
Head polarity；
D, combinations of polarities is carried out according to following principle：Using the polarity of initial wavefront wavelet coefficient as reference polarity, if measuring
First arrival wave head wavelet coefficient is "-", then second arrival wave head wavelet coefficient is then "+"；If measuring first arrival
Wave head wavelet coefficient is positive, then second arrival wave head wavelet coefficient is then respectively "-"；
If e, measure first arrival wave head wavelet coefficient, second arrival wave head wavelet coefficient polarity it is opposite when, 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 measure first arrival wave head wavelet coefficient, second arrival wave head wavelet coefficient polarity it is identical when, 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 CN105699855A (en) | 2016-06-22 |
| CN105699855B true 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) |
Families Citing this family (14)
| 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 |
| 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 |
| CN111433616B (en) * | 2017-11-17 | 2022-11-08 | 日立能源瑞士股份公司 | 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 |
| 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 |
| 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 |
| CN109901016B (en) * | 2019-03-13 | 2020-11-20 | 清华四川能源互联网研究院 | Double-end ranging method for eliminating wave velocity influence |
| CN110927510B (en) * | 2019-10-17 | 2021-10-01 | 东北大学 | Frequency domain method for power transmission line double-end traveling wave fault location |
| CN111766470B (en) * | 2020-06-24 | 2021-06-25 | 湖南大学 | Fault positioning method and system for high-voltage direct-current transmission line and direct-current transmission line |
| CN114113879B (en) * | 2021-11-04 | 2024-04-12 | 云南电网有限责任公司 | Low-voltage line fault positioning method and system based on wavelet transformation |
| CN114184884B (en) * | 2021-11-23 | 2024-05-24 | 昆明理工大学 | Automatic construction method of power grid fault traveling wave ranging equation |
| CN117233529B (en) * | 2023-09-20 | 2024-05-14 | 青岛鼎信通讯科技有限公司 | Single-ended distance measurement method suitable for inter-phase short-circuit fault of power distribution network |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
Family Cites Families (1)
| 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 |
-
2016
- 2016-04-06 CN CN201610209799.4A patent/CN105699855B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
|---|
| 一种不受波速影响的单端行波故障测距新算法;杨关春，彭丽;《山东理工大学学报》;20130131;第27卷(第1期);第59-62页 * |
| 输电线路故障暂态信号分析和行波测距研究;马丹丹;《中国优秀硕士论文全文数据库 工程科技Ⅱ辑》;20081215;论文正文第29-31页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105699855A (en) | 2016-06-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105699855B (en) | Based on the single-ended traveling wave fault location calculation method not influenced by traveling wave speed and distance measuring method | |
| Saleh et al. | Ultra-high-speed traveling-wave-based protection scheme for medium-voltage DC microgrids | |
| Suonan et al. | A novel fault-location method for HVDC transmission lines | |
| Ngu et al. | A combined impedance and traveling wave based fault location method for multi-terminal transmission lines | |
| Shu et al. | Single-ended fault location for hybrid feeders based on characteristic distribution of traveling wave along a line | |
| Saleh et al. | Fault detection and location in medium‐voltage DC microgrids using travelling‐wave reflections | |
| CN102074937B (en) | Ranging method for adaptive transmission line fault of model | |
| CN104898021B (en) | A kind of distribution network fault line selection method based on k means cluster analyses | |
| Azizi et al. | Fault location on multi-terminal DC systems using synchronized current measurements | |
| CN106019079B (en) | A kind of common-tower double-return DC line novel double end fault distance-finding method | |
| CN103777115B (en) | Electric transmission line single-terminal positioning method based on fault transient state and steady-state signal wave velocity difference | |
| CN110346688A (en) | A kind of high pressure complexity cable fault quick segmentation differentiates and localization method | |
| CN103235237A (en) | High-voltage direct-current grounding electrode line high-resistance fault positioning method | |
| CN109643890A (en) | Method and system for being positioned to the failure in mixed type transmission line of electricity | |
| CN106771868A (en) | A kind of flexible direct current distribution system fault locating method based on traveling wave | |
| CN102545177A (en) | Bergeron-model-based simulation-after-test method for fault phase selection of alternating current transmission line | |
| CN102565629B (en) | A kind of transmission line of alternation current Fault Phase Selection test simulation method based on lumped parameter Π model | |
| CN103869221A (en) | SV network sampling-based double-circuit line single-phase earth fault distance measurement method | |
| Jensen et al. | Online fault location on AC cables in underground transmission systems using sheath currents | |
| CN103954894A (en) | Partial discharge locating method for three-phase crossed and interconnected cables | |
| Xinzhou et al. | Optimizing solution of fault location | |
| CN104655981A (en) | Comparison type directional protection method for travelling wave polarity by utilizing crossing sequence differential transformation | |
| US20190094290A1 (en) | Traveling wave identification using distortions for electric power system protection | |
| CN106291262A (en) | The detection and location method of Distribution Network Frame ceases to be busy singlephase earth fault | |
| Yao et al. | A novel method to locate a fault of transmission lines by shielding failure |
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 |