CN108896874A - A kind of fault positioning method for transmission line of end connection short-term road - Google Patents
A kind of fault positioning method for transmission line of end connection short-term road Download PDFInfo
- Publication number
- CN108896874A CN108896874A CN201810756249.3A CN201810756249A CN108896874A CN 108896874 A CN108896874 A CN 108896874A CN 201810756249 A CN201810756249 A CN 201810756249A CN 108896874 A CN108896874 A CN 108896874A
- Authority
- CN
- China
- Prior art keywords
- wave
- transmission line
- head
- line
- current traveling
- 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/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
-
- 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
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 discloses a kind of fault positioning method for transmission line of end connection short-term road.This method needs respectively to arrange a current traveling wave measurement point in route head end, midpoint and end, but not seek precise synchronization.Firstly, in conjunction with wave head polarity, determining that failure occurs in short-term road or transmission line of electricity using the absolute difference of route both ends modulus traveling-waves time difference；Secondly, when failure occurs to determine fault section using the ratio of route both ends modulus traveling-waves reaching time-difference in transmission line of electricity；Finally, calculating three measurement dotted line mould current traveling wave the first two wave head reaching time-differences, and compare with theoretical value, the source of second wave head is determined, using the relational implementation precision ranging of the first two wave head reaching time-difference and fault distance.The present invention is not necessarily to each measurement point precise synchronization, is not influenced by short-term road, abort situation, fault resstance, failure initial phase angle, precision with higher, reliability.
Description
Technical field
The present invention relates to transmission line malfunction fields of measurement, especially a kind of transmission line malfunction of end connection short-term road
Distance measuring method.
Background technique
With fast development economic in recent years, the scale and quantity of power load are being continuously increased, the use of certain areas
Electric load may concentrate in several lesser regions, this makes that several can be connected to the transmission line of electricity end of its long distance powedr transmission
Shorter route.End is connected to the quick and precisely fault location of the transmission line of electricity of short-term road for reducing economic loss and improving system
System reliability is of great significance.
Currently, power transmission line fault locating method can be divided into impedance method and traveling wave method.Impedance method utilizes fault distance and meter
The relationship for calculating impedance carries out fault location, although method is simple, the influence vulnerable to factors such as fault resstance, electric arcs.Traveling wave
Method generally by detection initial transient traveling wave Mintrop wave head reach transmission line of electricity both ends measurement point at the time of, using these moment with
Internal logical relationship between fault distance carries out fault location, and positioning accuracy is high, is not influenced by factors such as electric arcs, using more
Come wider.In recent years, with the development of signal processing technology, traveling wave method, which has been got back, further to be developed.
Currently, end be connected to the measuring distance of transmission line fault of short-term road there are the problem of mainly have：(1) it is difficult to judge failure
Occur on transmission line of electricity or short-term road；(2) it is difficult to judge the source of second traveling wave reflection wave head；(3) for proximal end event
Barrier is difficult to precision ranging.Therefore research simple, practical, high reliablity the end of principle is connected to the transmission line of electricity event of short-term road
Barrier distance measuring method not only has theoretical research value, but also has important practical significance for engineering practice.
Summary of the invention
Goal of the invention：In order to solve the problems of the measuring distance of transmission line fault of current end connection short-term road, this
Invention proposes a kind of fault positioning method for transmission line of end connection short-term road.
To achieve the above object, the technical solution adopted by the present invention is：
A kind of fault positioning method for transmission line of end connection short-term road, this method include the steps that successively executing (1)
To (7)：
(1) measurement point is respectively set in head end, midpoint and the end of transmission line of electricity (not including end short-term road), remembers
Head end, midpoint, end measurement point be respectively Ps、Ph、Pr；Assuming that failure occurs in F1Point, it is dynamic in the protective device of transmission line of electricity
After work, it is utilized respectively measurement point Ps、Ph、PrThe current transformer of place's installation extracts current traveling wave signal, then to the electric current of extraction
Travelling wave signal carries out phase-model transformation, obtains corresponding line mould current traveling wave signal and zero mould current traveling wave signal；
(2) arrival time for obtaining each measurement point line mould, zero mould current traveling wave Mintrop wave head, remember measurement point PsPlace's line mould,
The reaching time-difference of zero mould current traveling wave Mintrop wave head is Δ ts_F1, measurement point PrLocate the arrival of line mould, zero mould current traveling wave Mintrop wave head
Time difference is Δ tr_F1；
(3) according to the calculated result of step (2), Δ t is calculateds_F1With Δ tr_F1The absolute value delta t of differenceF1：
ΔtF1=| Δ ts_F1-Δtr_F1|
(4) when the arbitrary point F in short-term road occurs for failure2When, theoretically measurement point PsAnd PrLocate line mould, zero mould electric current row
The difference of wave arrival time is Δ tF2：
Wherein, v1And v0Respectively line mould velocity of wave and zero mould velocity of wave, Ls-rFor the length of transmission line of electricity；
(5) determine that position, including step occur for failure：
(5-1) judges Δ tF1With Δ tF2Whether formula is met：
|ΔtF1-ΔtF2| < δ
Wherein, δ is preset threshold value.
If satisfied, then determine that failure occurs to execute step (5-3) in transmission line of electricity, it is no to then follow the steps (5-2)；
(5-2) judges measurement point PsAnd PrWhether the line line wave Mintrop wave head polar relationship measured meets：
Ws=-Wr
Wherein, WsAnd WrRespectively indicate measurement point PsAnd PrThe initial wave head polarity factor of line line wave at place；If satisfied, then
Determine that failure occurs to be transferred to step (5-3) in transmission line of electricity；If not satisfied, then determining that failure occurs in short-term road, the failure
Distance measuring method terminates；
(5-3) calculates Δ t when judgement is out of order generation on the transmission lines_F1With Δ tr_F1Ratio, according to calculating
As a result judge fault section：
Further, after judgement is out of order generation on the transmission line, event is carried out according to the calculated result of step (5)
Hinder precision ranging, including step：
(a) remember Ln-sAnd Lm-rRespectively shortest line length in the connected short-term road of transmission line of electricity head end and end, sentences
It is disconnected whether to meet Ln-s≤Lm-r；If so, being transferred to step (b), otherwise, it is transferred to step (c)；
(b) when failure occurs in front half section route, failure precision ranging equation is as follows：
Wherein, xestPositional distance transmission line head end bus B occurs for failuresDistance；ε is a small positive number；ts、
thAnd trRespectively measurement point Ps、PhAnd PrThe time difference that line mould current traveling wave the first two wave head that place measures reaches；
When failure occurs in second half section route, failure precision ranging equation is as follows：
(c) when failure occurs in front half section route, failure precision ranging equation is as follows：
When failure occurs in second half section route, failure precision ranging equation is as follows：
Further, known in the step (2) by the wave head that Teager energy operator and wavelet transform combine
Other method obtains the arrival time of each measurement point line mould, zero mould current traveling wave Mintrop wave head, including step：
1. obtaining the line mould current traveling wave signal of measurement point；
2. using db6 small echo, 4 layers of wavelet decomposition are carried out to obtained line mould current traveling wave signal, extract d1Layer details system
Number；
3. using same morther wavelet to d1Layer detail coefficients carry out wavelet reconstruction, obtain reconstructing detail coefficients accordingly,
It is denoted as：D [j]=[d1, d2..., dk], wherein j is the length of original line mode voltage travelling wave signal；
4. calculating the Teager energy vectors Te of the reconstruction coefficients：
Te ([d [j]])=d [j]2-d[j-1]·d[j+1]
In vector T e, as Mintrop wave head arrival time T at the time of maximum value element corresponds to1, second local maximum
As second wave head arrival time T at the time of element corresponds to2。
Further, in the step (5), the value of threshold value δ is calculated by following formula and sampling precision：
Further, in the step (5), steps are as follows for the calculating of the current traveling wave Mintrop wave head polarity factor：
1. obtaining transmission line of electricity first and last end line mould current traveling wave signal；
2. with line mould current traveling wave Mintrop wave head arrival time T1For initial time, N is taken backwardsThe signal of a data window length
As original traveling wave Mintrop wave head signal；
3. the route head end and end line mould current traveling wave Mintrop wave head signal kept in mind are respectively Is(k) and Ir(k) (k=1,
2 ..., Ns)；Calculate reference signal Iref(k) it is：
Wherein, e (k) is and Is(k) or Ir(k) the identical unit vector of length；
4. calculating Iref(k) corresponding difference signal：
Wherein, Isd(k) or IrdIt (k) is respectively route head end and end line mould current traveling wave Mintrop wave head difference signal；
5. calculating the line mould current traveling wave Mintrop wave head polarity factor of corresponding measuring point：
Wherein, w is the polarity factor of corresponding measuring point；IdIt (k) is the difference signal of corresponding measuring point；Sign () is symbol letter
Number.
Further, the ε takes the 1 μ s of minimum sampling interval under sample frequency 1MHz
Beneficial effect：The present invention is not necessarily to each measurement point precise synchronization, is not influenced by transition resistance and fault initial angle, energy
Quickly, fault zone is recognized accurately, and carries out failure precise positioning, reliability with higher and practical engineering value.
Detailed description of the invention
Fig. 1 is flow chart of the present invention；
Fig. 2 is the topology diagram for the transmission line of electricity that end connects short-term road；
Fig. 3 is the transmission line travelling wave propagation path figure that end connects short-term road；
Fig. 4 is l1-2And l2-4Near terminal fault route first and last end measuring point difference signal figure occurs.
Specific embodiment
The present invention will be further explained with reference to the accompanying drawing.
The topological structure of the transmission line of electricity of one typical end connection short-term road is as shown in Fig. 2, the invention proposes one
The fault positioning method for transmission line of kind end connection short-term road, overall flow is as shown in Figure 1, include the following steps：
(1) measurement point is respectively set in head end, midpoint and the end of transmission line of electricity (not including end short-term road), remembers
Head end, midpoint, end measurement point be respectively Ps、Ph、Pr；Assuming that failure occurs in F1Point, it is dynamic in the protective device of transmission line of electricity
After work, it is utilized respectively measurement point Ps、Ph、PrThe current transformer of place's installation extracts current traveling wave signal, then to the electric current of extraction
Travelling wave signal carries out phase-model transformation, obtains corresponding line mould current traveling wave signal and zero mould current traveling wave signal；
(2) the wave head recognition methods combined using wavelet transform and Teager energy operator, obtains each measurement point
Line mould, zero mould current traveling wave Mintrop wave head arrival time calculate measurement point PsAnd PrLocate modulus current traveling wave reaching time-difference Δ ts_F1
With Δ tr_F1：
Δts_F1=ta_s-tz_s
Δtr_F1=ta_r-tz_r
Wherein, ta_s、tz_sAnd ta_r、tz_rIt respectively represents in measurement point PsAnd PrThe reality of the line mould of acquisition, zero mode voltage traveling wave
Border arrival time.
(3) according to the calculated result of step (2), Δ t is calculateds_F1With Δ tr_F1The absolute value delta t of differenceF1：
ΔtF1=| Δ ts_F1-Δtr_F1|
Wherein, Δ ts_F1With Δ tr_F1It is P respectivelysAnd PrLocate modulus current traveling wave reaching time-difference.
(4) when failure occurs in short-term road arbitrary point F2When, theoretically measurement point PsAnd PrLocate line mould, zero mould current traveling wave
The difference of arrival time is Δ tF2：
Wherein, v1And v0For line mould and zero mould velocity of wave, Ls-rIt is the length of transmission line of electricity.
(5) determine that failure occurs in transmission line of electricity rather than the criterion of both ends short-term road is as follows：
(a) compare Δ tF1With Δ tF2Between difference：
|ΔtF1-ΔtF2| < δ
Wherein, δ is preset positive number, and value is calculated by following formula and sampling precision, is taken as 44ps in the present embodiment.
Wherein, v1And v0For line mould and zero mould velocity of wave, value takes 2.98 × 10 respectively8M/s and 2.8 × 108m/s。Ls-rIt is
The length of transmission line of electricity takes 200km in emulation.
(b) measurement point P is calculatedsAnd PrThe line line wave Mintrop wave head polarity factor at place, judges the initial traveling wave wave of electric current at two
Whether the polarity of head meets following formula relationship：
Ws=-Wr
Wherein, WsAnd WrRepresent measuring point PsAnd PrLocate the line mould current traveling wave Mintrop wave head polarity factor.
If measurement point PsAnd PrThe modulus time difference Δ t at placeF1With Δ tF2The relationship for meeting criterion (a), judges failure
On the transmission line, if not satisfied, criterion (b) is recycled, according to measurement point PsAnd PrLocate the polarity of line mould current traveling wave Mintrop wave head
Abort situation is judged, if the two polarity is on the contrary, failure occurs on the transmission line, conversely, failure occurs on short-term road.
(6) when judging that failure occurs on the transmission line, Δ t is calculateds_F1With Δ tr_F1Ratio, utilize following criterion
Judge fault section：
(7) assume in transmission line of electricity both ends connection short-term road, bus Bn-s, bus Bm-rWith bus BsBetween short-term road be
Shortest two in all short-term roads, length is respectively Ln-sAnd Lm-r, and Ln-s≤Lm-r.P on front half section route1And P2It is two
A critical point, midpoint P1With measurement point PsThe distance between be equal to Ln-s；Point P2With measurement point PsThe distance between be equal to Lm-r。
Front half section route is divided into three track sections by the two critical points.
When failure occurs in point P2With measurement point PhBetween any point F on track section3(F3With measurement point PsBetween away from
From being denoted as xa) when, due to xa> Lm-r, measure PsAnd PrLocating second line mould current traveling wave wave head is through bus Bn-sWith Bm-rIt is anti-
Ejected wave.At this point, the measurement point P among routehSecond wavefront that place detects is through bus BsBack wave, it is available
PhThe first two line mould wavefront that place detects calculates abort situation；When failure occurs in point P1And P2Between on track section
Any point F2(F2With measurement point PsThe distance between be denoted as xb) when, due to Lm-r> xb> Ln-s, measurement point PsWhat place detected
Second wavefront is through bus Bn-sBack wave；Measurement point PrSecond wavefront that place detects is by opposite end
Bus BsBack wave, using measurement point PrThe first two line mould wavefront that place detects calculates abort situation；When failure is sent out
Life is in measuring point PsWith point P1Between any point F on track section1(F1With measurement point PsThe distance between be denoted as xc) when, due to xc
> Ln-s, measurement point PsSecond wavefront that place detects is through fault point F1Transmitted wave；Measurement point PrWhat place detected
Second wavefront is opposite end bus BsBack wave, measurement point PsAnd PrThe first two line mould wavefront that place detects is all
It can be used to calculate fault distance.
So failure occurs in front half section route, failure precision ranging equation is as follows：
Wherein, xestPositional distance transmission line head end bus B occurs for failuresDistance；ε is a small positive number, is taken
1 μ s of minimum sampling interval under sample frequency 1MHz；ts、thAnd trRespectively measurement point Ps、PhAnd PrThe line mould electric current row that place measures
The time difference that wave the first two wave head reaches；
When failure occurs in second half section route, failure precision ranging equation is as follows：
If Ln-s≥Lm-r, in front half section route, failure precision ranging equation is as follows for failure generation：
If Ln-s≥Lm-r, in second half section route, failure precision ranging equation is as follows for failure generation：
In above scheme, wavefront, which is identified by the wave head that Teager energy operator and wavelet transform combine, to be known
Other method realizes, Teager energy operator be it is nonlinear, it can effectively reflect signal amplitude, frequency and instantaneous energy
Significant changes.Therefore, the present invention combines wavelet transform and Teager energy operator, obtains electric current modulus traveling-waves letter
Number arrival time, this with using noise-containing single wavelet transformation compared with, have better effect.The detailed step of this method
It is as follows：
1. obtaining measurement dotted line mould current traveling wave signal；
2. using db6 small echo, 4 layers of wavelet decomposition are carried out to obtained signal, extract d1Layer detail coefficients；
3. using same morther wavelet to d1Layer detail coefficients carry out wavelet reconstruction, obtain reconstructing detail coefficients accordingly,
It is denoted as：D [j]=[d1, d2..., dk], wherein j is the length of original line mode voltage travelling wave signal；
4. calculating the Teager energy vectors Te of the reconstruction coefficients：
Te ([d [j]])=d [j]2-d[j-1]·d[j+1]
In vector T e, as Mintrop wave head arrival time T at the time of maximum value element corresponds to1, second local maximum
As second wave head arrival time T at the time of element corresponds to2。
The calculation method of the current traveling wave Mintrop wave head polarity factor is：
In order to determine the polarity of route both ends measurement point current traveling wave Mintrop wave head, propose that the wavefront polarity factor calculates
Method.The algorithm goes out the symbol of the polarity factor according to the difference signal Integration Solving of traveling wave Mintrop wave head signal to determine wave head polarity.
Steps are as follows for specific method：
1. obtaining route first and last end line mould current traveling wave signal；
2. with current traveling wave Mintrop wave head arrival time T1For initial time, N is taken backwardsThe signal conduct of a data window length
Original traveling wave Mintrop wave head signal.N in emulationsTake 50.
8. the route head end and end traveling wave Mintrop wave head signal kept in mind are respectively Is(k) and Ir(k) (k=1,2 ...,
Ns).Reference signal I is calculated according to the following formularef(k)：
Wherein, e (k) is and Is(k) or Ir(k) the identical unit vector of length.
4. calculating corresponding difference signal：
Wherein, Isd(k) or IrdIt (k) is respectively route head end and end traveling wave Mintrop wave head difference signal.
5. calculating the traveling wave Mintrop wave head polarity factor of corresponding measuring point：
Wherein, w is the polarity factor of corresponding measuring point；IdIt (k) is the difference signal of corresponding measuring point；Sign () is symbol letter
Number.Simulating, verifying
In order to examine effectiveness of the invention and reliability, the defeated of end connection short-term road is built on PSCAD/EMTDC
Electric line simulation model, the transmission line travelling wave propagation path that end connects short-term road are as shown in Figure 3.Transmission line of electricity overall length is
200km, system fundamental frequency are 60Hz, voltage class 500kV.Route left and right ends 2 short-term roads of each band, length are respectively：
L1-2=10km, L3-2=30km, L4-5=40km, L4-6=25km.Transmission line of electricity, which uses, meets actual frequency dependent character model.It leads
The model of line and lightning conducter is respectively：4 × LGJ-400/35 and GJ-80.Power supply and transformer parameter are as shown in the figure.In emulation,
Three current traveling wave measurement points are set at route head end bus 2, midpoint and end bus 4 and (are expressed as P2、PmidAnd P4),
Its sample frequency is 1MHz.Transmission line of electricity is divided into two sections by these three measurement points, first segment P2And PmidBetween part of path,
Second segment is PmidAnd P4Between part of path.Failure is with respect to range error erIt is defined as follows：
Wherein, xrFor physical fault distance.Ls-rFor transmission line of electricity total length, value 200km.
In order to verify the validity of proposition method, respectively between bus 1 and bus 2 short-term road apart from bus 2 be 2km,
Apart from bus 1 to simulate singlephase earth fault at 2km, 15km and 50km on transmission line of electricity, fault resstance is 10 Ω, failure first phase
Angle is 90 °.Table 1 is the faulty line definitive result after above four positions are broken down, l1-2And l2-4Near terminal fault line occurs
Road first and last end measuring point difference signal is as shown in Figure 4, wherein FLN indicates faulty line；Δt2With Δ t4Respectively at bus 2 and 4
The modulus current traveling wave reaching time-difference being calculated；l1-2And l2-4Respectively indicate the short-term road and bus 2 between bus 1 and 2
And the transmission line of electricity between 4；FLDR representing fault line location result；FD is fault point distance transmission route survey point P2Away from
From；CS1For determine failure occur in transmission line of electricity rather than the criterion symbol of both ends short-term road；W is polarity because of subvector, includes
Two elements, respectively measuring point P2And P4The line mould current traveling wave Mintrop wave head polarity factor.Table 2 is transmission line of electricity l2-4On break down
When fault localization result.Wherein, LS1Indicate first part of path.t2、tmidAnd t4Respectively indicate measurement point P2、PmidAnd P4Locate line
Mould current traveling wave the first two wave head reaching time-difference；CS2For the selecting criterion symbol of failure precision ranging equation；FDC is to calculate
Fault distance.Proposed method can accurately calculate fault distance for the failure of different location, not by route both ends
The influence of short-term road.
1 faulty line definitive result of table
2 failure precision ranging result of table
The above is only a preferred embodiment of the present invention, it should be pointed out that：For the ordinary skill people of the art
For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also answered
It is considered as protection scope of the present invention.
Claims (6)
1. a kind of fault positioning method for transmission line of end connection short-term road, which is characterized in that include the steps that successively executing
(1) to (7)：
(1) measurement point is respectively set in the head end of transmission line of electricity, midpoint and end, remembers the measurement point of head end, midpoint, end
Respectively Ps、Ph、Pr；Assuming that failure occurs in F1Point is utilized respectively measurement point P after the protective device movement of transmission line of electricitys、
Ph、PrThe current transformer of place's installation extracts current traveling wave signal, then carries out phase-model transformation to the current traveling wave signal of extraction,
Obtain corresponding line mould current traveling wave signal and zero mould current traveling wave signal；
(2) arrival time for obtaining each measurement point line mould, zero mould current traveling wave Mintrop wave head, remember measurement point PsLocate line mould, zero mould electricity
The reaching time-difference of popular wave Mintrop wave head is Δ ts_F1, measurement point PrLocate the reaching time-difference of line mould, zero mould current traveling wave Mintrop wave head
For Δ tr_F1；
(3) according to the calculated result of step (2), Δ t is calculateds_F1With Δ tr_F1The absolute value delta t of differenceF1：
ΔtF1=| Δ ts_F1-Δr_F1|
(4) when the arbitrary point F in short-term road occurs for failure2When, measurement point P in computational theorysAnd PrPlaceLine mould, zero mould current traveling wave
The difference of arrival time is Δ tF2：
Wherein, v1And v0Respectively line mould velocity of wave and zero mould velocity of wave, Ls-rFor the length of transmission line of electricity；
(5) determine that position, including step occur for failure：
(5-1) judges Δ tF1With Δ tF2Whether formula is met：
|AtF1-ΔtF2| < δ
Wherein, δ is preset threshold value；
If satisfied, then determine that failure occurs to execute step (5-3) in transmission line of electricity, it is no to then follow the steps (5-2)；
(5-2) judges measurement point PsAnd PrWhether the initial wave head polar relationship of line line wave measured meets：
Ws=-Wr
Wherein, WsAnd WrRespectively indicate measurement point PsAnd PrThe line mould current traveling wave Mintrop wave head polarity factor at place；If satisfied, then determining
Failure occurs to be transferred to step (5-3) in transmission line of electricity；If not satisfied, then determining that failure occurs in short-term road, the fault localization
Method terminates；
(5-3) calculates Δ t when judgement is out of order generation on the transmission lines_F1With Δ tr_F1Ratio, according to calculated result
Judge fault section：
2. a kind of fault positioning method for transmission line of end connection short-term road according to claim 1, which is characterized in that
After judgement is out of order generation on the transmission line, failure precision ranging, including step are carried out according to the calculated result of step (5)
Suddenly：
(a) remember Ln-sAnd Lm-rRespectively shortest line length in the connected short-term road of transmission line of electricity head end and end, judgement be
It is no to meet Ln-s≤Lm-r；If so, being transferred to step (b), otherwise, it is transferred to step (c)；
(b) when failure occurs in front half section route, failure precision ranging equation is as follows：
Wherein, xestPositional distance transmission line head end bus B occurs for failuresDistance；ε is a small positive number；ts、thWith
trRespectively measurement point Ps、PhAnd PrThe time difference that line mould current traveling wave the first two wave head that place measures reaches；
When failure occurs in second half section route, failure precision ranging equation is as follows：
(c) when failure occurs in front half section route, failure precision ranging equation is as follows：
When failure occurs in second half section route, failure precision ranging equation is as follows：
3. a kind of fault positioning method for transmission line of end connection short-term road according to claim 2, which is characterized in that
The wave head recognition methods combined in the step (2) by Teager energy operator and wavelet transform obtains each measurement point
Locate the arrival time of line mould, zero mould current traveling wave Mintrop wave head, including step：
1. obtaining the line mould current traveling wave signal of measurement point；
2. using db6 small echo, 4 layers of wavelet decomposition are carried out to obtained line mould current traveling wave signal, extract d1Layer detail coefficients；
3. using same morther wavelet to d1Layer detail coefficients carry out wavelet reconstruction, obtain reconstructing detail coefficients accordingly, be denoted as：d
[j]=[d1, d2..., dk], wherein j is the length of original line mode voltage travelling wave signal；
4. calculating the Teager energy vectors Te of the reconstruction coefficients：
Te ([d [j]])=d [j]2-d[j-1]·d[j+1]
In vector T e, as Mintrop wave head arrival time T at the time of maximum value element corresponds to1, second local maximum element pair
As second wave head arrival time T at the time of answering2。
4. a kind of fault positioning method for transmission line of end connection short-term road according to claim 3, which is characterized in that
In the step (5), the value of threshold value δ is calculated by following formula and sampling precision：
5. a kind of fault positioning method for transmission line of end connection short-term road according to claim 4, which is characterized in that
In the step (5), steps are as follows for the calculating of the current traveling wave Mintrop wave head polarity factor：
1. obtaining transmission line of electricity first and last end line mould current traveling wave signal；
2. with line mould current traveling wave Mintrop wave head arrival time T1For initial time, N is taken backwardsThe signal conduct of a data window length
Original traveling wave Mintrop wave head signal；
8. the route head end and end line mould current traveling wave Mintrop wave head signal kept in mind are respectively Is(k) and Ir(k) (k=1,
2 ..., Ns)；Calculate reference signal Iref(k) it is：
Wherein, e (k) is and Is(k) or Ir(k) the identical unit vector of length；
4. calculating Iref(k) corresponding difference signal：
Wherein, Isd(k) or IrdIt (k) is respectively route head end and end line mould current traveling wave Mintrop wave head difference signal；
5. calculating the line mould current traveling wave Mintrop wave head polarity factor of corresponding measuring point：
Wherein, w is the polarity factor of corresponding measuring point；IdIt (k) is the difference signal of corresponding measuring point；Sign () is sign function.
6. a kind of fault positioning method for transmission line of end connection short-term road according to claim 5, it is characterised in that：
The ε takes the 1 μ s of minimum sampling interval under sample frequency 1MHz.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810756249.3A CN108896874B (en) | 2018-07-10 | 2018-07-10 | A kind of fault positioning method for transmission line of end connection short-term road |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810756249.3A CN108896874B (en) | 2018-07-10 | 2018-07-10 | A kind of fault positioning method for transmission line of end connection short-term road |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108896874A true CN108896874A (en) | 2018-11-27 |
| CN108896874B CN108896874B (en) | 2019-09-24 |
Family
ID=64348930
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810756249.3A Active CN108896874B (en) | 2018-07-10 | 2018-07-10 | A kind of fault positioning method for transmission line of end connection short-term road |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108896874B (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109521330A (en) * | 2018-12-06 | 2019-03-26 | 中国矿业大学 | A kind of transmission line malfunction travelling wave ranging method based on the prediction of ARIMA wave head |
| CN109901016A (en) * | 2019-03-13 | 2019-06-18 | 清华四川能源互联网研究院 | A kind of both-end distance measuring method eliminating velocity of wave and influencing |
| CN110426592A (en) * | 2019-08-16 | 2019-11-08 | 南京国电南自电网自动化有限公司 | Aerial and cable hybrid line Earth design method based on the traveling wave time difference |
| CN110514967A (en) * | 2019-10-12 | 2019-11-29 | 国网江苏省电力有限公司徐州供电分公司 | A kind of Fault Locating Method based on NARNN model prediction wavefront arrival time |
| CN110703012A (en) * | 2019-09-26 | 2020-01-17 | 国电南瑞科技股份有限公司 | Distributed fault diagnosis method for power transmission line |
| CN110854828A (en) * | 2019-12-03 | 2020-02-28 | 西安交通大学 | Single-ended adaptive traveling wave ultra-high speed protection system and method for dealing with high-resistance fault |
| CN111239553A (en) * | 2020-02-17 | 2020-06-05 | 国网浙江省电力有限公司建设分公司 | Line wave velocity and line length calculation method suitable for hybrid power transmission line |
| CN112505472A (en) * | 2020-10-19 | 2021-03-16 | 国网辽宁省电力有限公司电力科学研究院 | Three-phase hybrid power supply network fault location method based on traveling wave |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6798211B1 (en) * | 1997-10-30 | 2004-09-28 | Remote Monitoring Systems, Inc. | Power line fault detector and analyzer |
| CN101776725A (en) * | 2010-01-13 | 2010-07-14 | 上海交通大学 | Fault positioning method for transmission line |
| CN101923139A (en) * | 2010-04-19 | 2010-12-22 | 昆明理工大学 | Intelligent method for single-ended traveling wave fault location of power transmission line |
| CN102096021A (en) * | 2010-12-08 | 2011-06-15 | 西南交通大学 | Traveling wave natural frequency-based power transmission network failure networking positioning and distance measurement method |
| CN102435908A (en) * | 2011-09-07 | 2012-05-02 | 兖州煤业股份有限公司 | On-line ranging method for single-phase earth fault based on zero mode and phase mode current traveling waves |
| CN103809079A (en) * | 2014-02-17 | 2014-05-21 | 华北电力大学 | Double-end high frequency impedance type fault ranging method suitable for direct current distribution network |
| CN104965151A (en) * | 2015-05-28 | 2015-10-07 | 国家电网公司 | Fault distance detecting method based on voltage fault component of fault point |
| CN106771861A (en) * | 2016-12-08 | 2017-05-31 | 中国矿业大学 | Complex electric network Fault Locating Method based on wide area traveling wave energy and time difference |
| CN106771843A (en) * | 2015-11-19 | 2017-05-31 | 中国石油化工股份有限公司 | A kind of fault travelling wave ranging method of single-core power cables |
| CN107505538A (en) * | 2017-09-22 | 2017-12-22 | 中国矿业大学 | Half-wavelength power transmission line asynchronous fault positioning method based on line mode voltage traveling wave amplitude attenuation characteristic |
| CN107632238A (en) * | 2017-08-25 | 2018-01-26 | 国网辽宁省电力有限公司 | A kind of multiterminal transmission line failure distance-finding method based on WAMS systems |
| CN108008251A (en) * | 2017-11-29 | 2018-05-08 | 国网江苏省电力公司无锡供电公司 | The fault distance-finding method of mixed power transmission line unknown parameters |
| CN108120899A (en) * | 2017-12-21 | 2018-06-05 | 国网宁夏电力公司中卫供电公司 | A kind of single-ended Section Location of one-phase earthing failure in electric distribution network |
-
2018
- 2018-07-10 CN CN201810756249.3A patent/CN108896874B/en active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6798211B1 (en) * | 1997-10-30 | 2004-09-28 | Remote Monitoring Systems, Inc. | Power line fault detector and analyzer |
| CN101776725A (en) * | 2010-01-13 | 2010-07-14 | 上海交通大学 | Fault positioning method for transmission line |
| CN101923139A (en) * | 2010-04-19 | 2010-12-22 | 昆明理工大学 | Intelligent method for single-ended traveling wave fault location of power transmission line |
| CN102096021A (en) * | 2010-12-08 | 2011-06-15 | 西南交通大学 | Traveling wave natural frequency-based power transmission network failure networking positioning and distance measurement method |
| CN102435908A (en) * | 2011-09-07 | 2012-05-02 | 兖州煤业股份有限公司 | On-line ranging method for single-phase earth fault based on zero mode and phase mode current traveling waves |
| CN103809079A (en) * | 2014-02-17 | 2014-05-21 | 华北电力大学 | Double-end high frequency impedance type fault ranging method suitable for direct current distribution network |
| CN104965151A (en) * | 2015-05-28 | 2015-10-07 | 国家电网公司 | Fault distance detecting method based on voltage fault component of fault point |
| CN106771843A (en) * | 2015-11-19 | 2017-05-31 | 中国石油化工股份有限公司 | A kind of fault travelling wave ranging method of single-core power cables |
| CN106771861A (en) * | 2016-12-08 | 2017-05-31 | 中国矿业大学 | Complex electric network Fault Locating Method based on wide area traveling wave energy and time difference |
| CN107632238A (en) * | 2017-08-25 | 2018-01-26 | 国网辽宁省电力有限公司 | A kind of multiterminal transmission line failure distance-finding method based on WAMS systems |
| CN107505538A (en) * | 2017-09-22 | 2017-12-22 | 中国矿业大学 | Half-wavelength power transmission line asynchronous fault positioning method based on line mode voltage traveling wave amplitude attenuation characteristic |
| CN108008251A (en) * | 2017-11-29 | 2018-05-08 | 国网江苏省电力公司无锡供电公司 | The fault distance-finding method of mixed power transmission line unknown parameters |
| CN108120899A (en) * | 2017-12-21 | 2018-06-05 | 国网宁夏电力公司中卫供电公司 | A kind of single-ended Section Location of one-phase earthing failure in electric distribution network |
Non-Patent Citations (2)
| Title |
|---|
| WANG CHONGLIN等: "Fault Location for Single-Phase-To-Earth Faults Based on Transient Traveling Wave Method and Artificial Pulse Signal Injection Method", 《2010 INTERNATIONAL CONFERENCE ON ELECTRICAL AND CONTROL ENGINEERING》 * |
| 梁睿 等: "基于广域行波初始波头时差关系矩阵的配电网故障选线及测距", 《高电压技术》 * |
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109521330A (en) * | 2018-12-06 | 2019-03-26 | 中国矿业大学 | A kind of transmission line malfunction travelling wave ranging method based on the prediction of ARIMA wave head |
| CN109521330B (en) * | 2018-12-06 | 2020-06-09 | 中国矿业大学 | Power transmission line fault traveling wave distance measurement method based on ARIMA wave head prediction |
| 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 |
| CN110426592A (en) * | 2019-08-16 | 2019-11-08 | 南京国电南自电网自动化有限公司 | Aerial and cable hybrid line Earth design method based on the traveling wave time difference |
| CN110703012A (en) * | 2019-09-26 | 2020-01-17 | 国电南瑞科技股份有限公司 | Distributed fault diagnosis method for power transmission line |
| CN110703012B (en) * | 2019-09-26 | 2021-10-15 | 国电南瑞科技股份有限公司 | Distributed fault diagnosis method for power transmission line |
| CN110514967A (en) * | 2019-10-12 | 2019-11-29 | 国网江苏省电力有限公司徐州供电分公司 | A kind of Fault Locating Method based on NARNN model prediction wavefront arrival time |
| CN110514967B (en) * | 2019-10-12 | 2021-05-18 | 国网江苏省电力有限公司徐州供电分公司 | Fault positioning method for predicting arrival time of traveling wave head based on NARNN model |
| CN110854828B (en) * | 2019-12-03 | 2021-01-19 | 西安交通大学 | Single-ended adaptive traveling wave ultra-high speed protection system and method for dealing with high-resistance fault |
| CN110854828A (en) * | 2019-12-03 | 2020-02-28 | 西安交通大学 | Single-ended adaptive traveling wave ultra-high speed protection system and method for dealing with high-resistance fault |
| CN111239553A (en) * | 2020-02-17 | 2020-06-05 | 国网浙江省电力有限公司建设分公司 | Line wave velocity and line length calculation method suitable for hybrid power transmission line |
| CN112505472A (en) * | 2020-10-19 | 2021-03-16 | 国网辽宁省电力有限公司电力科学研究院 | Three-phase hybrid power supply network fault location method based on traveling wave |
| CN112505472B (en) * | 2020-10-19 | 2023-01-24 | 国网辽宁省电力有限公司电力科学研究院 | Three-phase hybrid power supply network fault location method based on traveling wave |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108896874B (en) | 2019-09-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108896874B (en) | A kind of fault positioning method for transmission line of end connection short-term road | |
| CN110187235B (en) | Distributed power line fault positioning system and method based on traveling wave speed dynamic measurement | |
| Jiao et al. | A new method to improve fault location accuracy in transmission line based on fuzzy multi-sensor data fusion | |
| CN107505538B (en) | The asynchronous Fault Locating Method of half-wave power transmission circuit based on line mode voltage traveling wave amplitude attenuation characteristic | |
| CN101299538B (en) | Cable-aerial mixed line fault travelling wave ranging method | |
| CN105353268B (en) | One kind is used for the judgement of transmission line of electricity distribution traveling wave fault and localization method | |
| Liao et al. | Optimal estimate of transmission line fault location considering measurement errors | |
| CN110596539B (en) | Power distribution network fault positioning method based on transient waveform correlation | |
| CN107621591B (en) | A kind of transmission line of electricity iteration distance measuring method based on zero mould traveling wave speed variation characteristic | |
| CN106771861B (en) | Complex electric network Fault Locating Method based on wide area traveling wave energy and time difference | |
| CN104569744A (en) | Comprehensive single-end fault positioning method applicable to power distribution network lines | |
| JP2008157940A (en) | Multi-ended fault location specifying system | |
| CN103592575A (en) | Self-adaptation weighting data fusion fault distance measurement method based on multi-sensor system | |
| CN105652151B (en) | Both-end distance measuring method based on line parameter circuit value detection and the asynchronous verification of data | |
| Xinzhou et al. | Optimizing solution of fault location | |
| CN105137293A (en) | Positioning method of fault points in power distribution network mixed circuits | |
| CN108693446A (en) | A kind of Fault Locating Method of non-synchronous sampling power grid transient state travelling wave modulus time difference | |
| CN109521330A (en) | A kind of transmission line malfunction travelling wave ranging method based on the prediction of ARIMA wave head | |
| CN107632236A (en) | A kind of single outgoing-feeder line one-end fault ranging method based on the identification of opposite end bus back wave | |
| CN108627740A (en) | Half-wavelength power transmission line fault location method considering traveling wave speed change and arrival time compensation | |
| CN110514967A (en) | A kind of Fault Locating Method based on NARNN model prediction wavefront arrival time | |
| CN106841919A (en) | Based on the high-precision transmission line method of single end distance measurement that triangle is calculated | |
| CN108845233A (en) | Electric distribution network overhead wire single-phase earthing double-end monitor localization method | |
| CN103969553A (en) | Cable and overhead line mixed line double-end traveling wave fault location algorithm based on piecewise compensation principle | |
| CN103424667A (en) | Distance measurement method for post measurement simulation match of cable-hybrid line fault wave arriving time sequence |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |