CN108344923A - A kind of high-adaptability fault positioning method for transmission line and system - Google Patents
A kind of high-adaptability fault positioning method for transmission line and system Download PDFInfo
- Publication number
- CN108344923A CN108344923A CN201810051826.9A CN201810051826A CN108344923A CN 108344923 A CN108344923 A CN 108344923A CN 201810051826 A CN201810051826 A CN 201810051826A CN 108344923 A CN108344923 A CN 108344923A
- Authority
- CN
- China
- Prior art keywords
- back wave
- transmission line
- wave
- fault
- electricity
- 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
- 230000005540 biological transmission Effects 0.000 title claims abstract description 124
- 238000000034 method Methods 0.000 title claims abstract description 54
- 230000005611 electricity Effects 0.000 claims abstract description 79
- 230000002596 correlated effect Effects 0.000 claims abstract description 48
- 230000002441 reversible effect Effects 0.000 claims abstract description 45
- 230000009466 transformation Effects 0.000 claims abstract description 45
- 230000004807 localization Effects 0.000 claims abstract description 39
- 230000007257 malfunction Effects 0.000 claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 230000011514 reflex Effects 0.000 claims description 25
- 230000017105 transposition Effects 0.000 claims description 22
- 238000000605 extraction Methods 0.000 claims description 16
- 238000004458 analytical method Methods 0.000 claims description 15
- 230000001568 sexual effect Effects 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims description 13
- 230000001052 transient effect Effects 0.000 claims description 8
- 238000012360 testing method Methods 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 16
- 238000013519 translation Methods 0.000 description 11
- 238000010276 construction Methods 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 3
- 230000007704 transition Effects 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000001934 delay Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008447 perception Effects 0.000 description 2
- 230000008521 reorganization Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 238000012300 Sequence Analysis Methods 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000008188 pellet Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000028527 righting reflex Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
Abstract
The present invention provides a kind of high-adaptability fault positioning method for transmission line and systems, including：The three-phase current after transmission line malfunction is extracted, carrying out modal transformation to the three-phase current obtains line modulus；Multi-scale transform, the failure initial time at positioning transmission line of electricity both ends are carried out to line modulus；According to the failure initial time calculation delay at transmission line of electricity both ends；The positive cross-correlation sequence and reverse mutual correlated series at transmission line of electricity both ends are determined according to wavelet transform dimension and the time delay；Fault point back wave and/or busbar back wave are identified according to the positive cross-correlation sequence and reverse mutual correlated series at transmission line of electricity both ends；Fault localization calculating is carried out according to the fault point back wave and/or busbar back wave identified.Technical solution provided by the invention is improved the recognition success rate of fault point/busbar back wave by the time-frequency correlation of two end data of circuit, is modified on this basis to line length, time service error, improves the adaptability and precision of fault localization.
Description
Technical field
The invention belongs to measuring distance of transmission line fault fields, and in particular to a kind of high-adaptability measuring distance of transmission line fault side
Method and system.
Background technology
After transmission line of electricity breaks down, no matter whether circuit successfully overlaps, and track walker is required for search fault point, according to
Damaged condition caused by failure judges whether to need interruption maintenance to remove a hidden danger, therefore, failure point of power transmission line it is accurate fixed
Position (also commonly referred to as fault localization) is an important link for ensureing power grid security.According to the different existing transmission lines of electricity of principle
Fault localization can be divided into：Impedance method, traveling wave method and fault analytical method etc..Since the nineties in last century, based on the defeated of traveling wave principle
Electric system obtains extensive use to line fault distance measuring device at home, and practical operating experiences show：Relative to impedance method, event
Hinder analytic approach etc., traveling wave fault location preferably meets user's requirement in terms of precision and reliability.
According to the difference of use electrical quantity, traveling wave fault location is divided into as single-ended method, both-end method, impulse method etc., scene reality
The essentially both-end method of border application.Both-end method utilizes the initial travelling wave signal that failure generates, and is arrived by calculating the initial traveling wave of failure
Abort situation is calculated up to the circuit both ends time difference, fault distance l is calculated as follows1：
In above formula；t1、t2Respectively traveling wave reaches the absolute time at circuit both ends, and L is total track length, and v is that traveling wave propagates speed
Degree.Both-end traveling wave method has the characteristics that principle is simple and reliable, but there is also following its range accuracies of influence in Practical Project
Factor：(1) both ends time service device error：The time service error of 1us can lead to the range error of about 150m.For the device that put into operation,
Both ends time service error is one of the main reason of fault localization failure, and caused fault localization failure case is close thus
30%.In addition, in Practical Project, when circuit both ends use different time service producer systems, two end devices often exist several micro-
The fixation time service error of second.In DC transmission engineering, during man-made short-circuit tests, there is consolidating for about 3us in both sides time service device
Determine time service error, you can lead to about 400 meters of fixation range error.(2) line length error：The line error meeting of 1km long can
Lead to the fixation range error of about 0.5km.In Practical Project, field condition is more complicated, design, unit of operation's often basis
Line length is estimated in geographical location, and unit in charge of construction generally calculates line length according to consumptive material, and the two often has differences.This
Outside, arc sag difference caused by circular temperature loads can also influence line length.Such as the DC transmission engineering in nearly 600km long
In, the difference between the actual measurement line length and design length of unit in charge of construction even reaches 20km, this is to range accuracy sometimes up to 4km
Cause larger impact, it is however generally that, difference between the two is within 3%.(3) line construction error：Line construction error can shadow
Ring the velocity of wave setting of fault location device.It is more universal to having circuit reorganization and expansion situation in the power grid construction in some places,
Often occur circuit single loop line being changed to the change line construction such as double loop during reorganization and expansion.It is past in traditional travelling wave ranging method
Estimate that velocity of wave, line construction variation cause velocity of wave variation to generate additional error toward according to line construction.With 220kV circuits
For, there is about 3% velocity error in single loop line, double loop, velocity of wave influences and position of failure point, line length range accuracy
Correlation, when fault point is located at circuit both ends, velocity of wave error influences to maximize, and by taking 100km circuits as an example, limit influence about exists
1.5km left and right.(4) signal cable error in length in standing：Transmit similar, CT/PT secondary side signals on the line with transient state travelling wave
There is also time delays in signal cable transmission process in station.Initial failure range unit is installed in station in secondary small chamber, and two
Side signal time delay is close, is influenced on range accuracy smaller.With the development of secondary device site, part circuit side device peace
In secondary small chamber, and offside device site is installed, and secondary singal propagation delay time just impacts range accuracy, 300m
Left-right signal cable transmission time delay is about 2us, can lead to additional 300 meters of range errors.
The single-ended traveling wave method of fault localization, reason are completed using initial traveling wave and fault point/opposite end busbar back wave time difference
By it is upper not by line length error, time service device, stand in cable influenced.Therefore, some researchers propose to utilize single-ended row
Wave method coordinates both-end method to reduce the influence of above-mentioned factor, but also needs to solve the problems, such as that back wave identifies, in practical applications
Mainly face following difficulty：(1) theoretically, back wave recognition success rate can be improved by both-end method adding window, but in physical fault
In analysis, when window width setting it is relatively difficult, and signal resonance itself makes single adding window be difficult to solve back wave identification and ask
Topic.In addition, being influenced by factors such as time service errors by both-end method precision, non-cutting time window can not improve back wave and identify successfully
Rate.(2) being on the increase with power network line quantity, there is also the interference of reflected wave problem between the circuit adjacent with busbar, with
The polarity of back wave between the adjacent circuit of busbar is identical as fault point back wave, and which further increases the difficulties of perception reflex wave
Degree.(3) velocity of wave calculates：When measurement end of adjusting the distance localization of fault farther out, compared with both-end traveling wave method, failure in single-ended traveling wave method
Point/busbar reflective wave transmission is apart from longer so that the influence of velocity of wave error further expands.
In view of above-mentioned deficiency, single-ended traveling wave method is in physical fault analysis, by the single-ended method auxiliary both-end method correction of tradition
Drawback in length, time service error technology also relies primarily on engineering staff and has been manually done.
Invention content
The present invention is based on two end signal time-frequency correlations to disclose a kind of high-adaptability fault positioning method for transmission line.
A kind of high-adaptability fault positioning method for transmission line provided by the invention, including：
The three-phase current after transmission line malfunction is extracted, carrying out modal transformation to the three-phase current obtains line modulus；
Multi-scale transform is carried out to the line modulus, positions the failure initial time at the transmission line of electricity both ends；
According to the failure initial time calculation delay at the transmission line of electricity both ends；
The positive cross-correlation sequence at the transmission line of electricity both ends and reversed is determined according to wavelet transform dimension and the time delay
Cross-correlation sequence；
Fault point back wave is identified according to the positive cross-correlation sequence and reverse mutual correlated series at the transmission line of electricity both ends
And/or busbar back wave；
Fault localization calculating is carried out according to the fault point back wave and/or busbar back wave identified.
Three-phase current after the extraction transmission line malfunction carries out modal transformation to the three-phase current and obtains line mould
Amount, including：
Modal transformation is carried out as the following formula to the three-phase current, is transformed to line modulus：
In formula, Ia、Ib、IcFor three-phase current；IαFor the α modulus of three-phase current, IβFor the β modulus of three-phase current, I0It is three
Zero modulus of phase current.
It is described that multi-scale transform is carried out to line modulus, the failure initial time at the transmission line of electricity both ends is positioned, including：
Multi-scale transform is carried out to the line modulus with binary wavelet, according to wavelet transform dimension, passes through wavelet transformation
Line modulus is decomposed into approximation coefficient and detail coefficients；
Using the amplitude mean value of the detail coefficients of the wavelet transformation as standard, the frequency of access line both ends transient state amplitude maximum
Band is analyzed；
The wavelet transformation detail coefficients of the selected frequency band of extraction, and calculate modulus maximum sequence；
According to the modulus maximum sequence positioning failure initial time.
Time delay is calculated as follows according to the failure initial time at the transmission line of electricity both ends：
Δ t=tm-tn
In formula, Δ t is time delay, tmAnd tnFor the failure initial time at circuit both ends.
Determine the positive cross-correlation sequence at the transmission line of electricity both ends as the following formula according to wavelet transform dimension and the time delay：
W (t)=xm(t)×xn(t+Δt)
In formula, w (t) is the positive cross-correlation sequence at transmission line of electricity both ends, xm(t) and xn(t) it is respectively transmission line of electricity both ends
Wavelet conversion coefficient, Δ t be time delay.
The reverse mutual correlated series at the transmission line of electricity both ends are determined according to wavelet transform dimension and the time delay, including：
The analysis data that the either end in the transmission line of electricity both ends is calculated as follows carry out transposition：
y(t0+ i)=x (t0+ DL-i)
Wherein, DL=L/v × 1/fs
In formula, y (t) is the transmission line of electricity single ended data after transposition, and x (t) is the transmission line of electricity single ended data before transposition, t0
For wave head initial time, i=1,2...DL, DL is time window, and L is total track length, and v is velocity of wave, fsFor sample frequency；
The transmission line of electricity of transmission line of electricity single ended data and non-transposition after transposition is substituted into following formula to end data, is obtained described
The reverse mutual correlated series at transmission line of electricity both ends：
W (t) '=xm(t)×xn(t+Δt)
In formula, w (t) ' is the reverse mutual correlated series at transmission line of electricity both ends, xm(t) and xn(t) it is respectively transmission line of electricity two
The wavelet conversion coefficient at end, Δ t are time delay.
It is described to identify that fault point is anti-according to the positive cross-correlation sequence and reverse mutual correlated series at the transmission line of electricity both ends
Ejected wave and/or busbar back wave, including：
Include according to positive cross-correlation recognition sequence fault point back wave and busbar back wave, condition：
1) fault point back wave and busbar back wave are cathode sexual reflex；
2) the amplitude A of fault point back wave and busbar back wavefMeet 0.01 × A≤Af≤0.25×β1×β2× A,
In, A is initial wave head amplitude, β1、β2The respectively busbar reflectance factor at circuit both ends；With
3) fault point back wave and peer failure point reflection wave meet following formula：
0.97×L≤(t'1-t0)×v+(t'2-t0)×v≤1.03×L
In formula, t'1For fault point back wave moment, t'2For busbar back wave moment, t0For wave head initial time, L is line
Road overall length, v are velocity of wave.
Fault point back wave and/or busbar back wave are identified with reverse mutual correlated series, including：
It only can recognize that fault point back wave or busbar back wave, condition include according to the identification of reverse mutual correlated series：
1) fault point back wave or busbar back wave are cathode sexual reflex；
2) fault point back wave or busbar reflex amplitude AfThe following conditions should be met：0.01×A≤Af≤0.25×β1×
A or 0.01 × A≤Af≤0.25×β2×A。
Fault localization calculating is carried out according to the fault point back wave and busbar back wave identified, including：
If identifying fault point back wave and busbar back wave simultaneously, fault localization calculating is carried out as the following formula：
In formula, d is the distance value that test obtains, t'1For fault point back wave moment, t'2For busbar back wave moment, t0
For wave head initial time, L is total track length；
After completing fault localization calculating, the back wave of identification is corresponded respectively to by single ended data mould according to modulus maximum sequence
Maximum determines that back wave is fault point back wave or busbar back wave.
If only can recognize that fault point back wave or busbar back wave, fault localization calculating is carried out as the following formula：
D=(t'1-t0)×v/2
After completing fault localization calculating, the back wave of identification is corresponded to by single ended data modulus maxima according to modulus maximum sequence
Value, is determined as fault point or opposite end busbar back wave.
The present invention provides a kind of high-adaptability measuring distance of transmission line fault system, including：
Extraction module carries out modal transformation for extracting the three-phase current after transmission line malfunction to the three-phase current
Obtain line modulus；
Locating module, for carrying out multi-scale transform to the line modulus, when the failure at positioning transmission line of electricity both ends is initial
It carves；
Computing module, for the failure initial time calculation delay according to transmission line of electricity both ends；
Correlated series module, for determining the forward direction at transmission line of electricity both ends according to wavelet transform dimension and the time delay mutually
Close sequence and reverse mutual correlated series；
Identification module, for identifying failure according to the positive cross-correlation sequence and reverse mutual correlated series at transmission line of electricity both ends
Point reflection wave and/or busbar back wave；
Range finder module, for carrying out fault localization calculating according to the fault point back wave and/or busbar back wave that identify.
Compared with the latest prior art, technical solution provided by the invention has the advantages that：
Technical solution provided by the invention improves fault point/busbar by the time-frequency correlation of two end data of circuit and reflects
The recognition success rate of wave is on this basis modified line length, time service error, and then improves the adaptation of fault localization
Property and precision；
Technical solution provided by the invention carries out fault localization using fault point and/or busbar back wave, avoids circuit
Error in length, time service error, velocity of wave error, stand in influence of the cable to range accuracy, can be used as GPS and lose under magnitude abnormal conditions
Standby correction means；
Technical solution provided by the invention, using axial symmetry when the back wave of both ends and waveform similarity, and signal resonance,
With the orthogonal feature of the interference volumes such as busbar adjacent lines back wave, based on the positive, anti-of small echo correlation analysis the two poles of the earth signal
Clutter and interferer suppression is realized to cross correlation, solves the problems, such as most important back wave identification in Single Terminal Traveling Wave Fault Location；
Technical solution provided by the invention, relative to the method that existing single-ended/both-end method combines, the method for the invention is not
It is influenced by both-end method range error, window width offering question when not having to consider.
Description of the drawings
Fig. 1 is a kind of high-adaptability fault positioning method for transmission line flow chart provided by the invention；
Fig. 2 is the when axial symmetry schematic diagram of two end signal of transmission line of electricity；
Fig. 3 is δ and transition resistance variation tendency schematic diagram in the embodiment of the present invention；
Fig. 4 is circuit wave conversion schematic diagram small in ends in the embodiment of the present invention four；
(a) it is circuit both ends wavelet conversion coefficient schematic diagram, is (b) that circuit both ends wavelet modulus maxima sequence is illustrated
Figure；
Fig. 5 is the positive cross-correlation sequence diagram of circuit both ends data in the embodiment of the present invention four；
Fig. 6 is the reverse mutual correlated series schematic diagram of circuit both ends data in the embodiment of the present invention four；
Fig. 7 is circuit wave conversion schematic diagram small in ends in the embodiment of the present invention five；
(a) it is circuit both ends wavelet conversion coefficient schematic diagram, is (b) that circuit both ends wavelet modulus maxima sequence is illustrated
Figure；
Fig. 8 is the positive cross-correlation sequence diagram of circuit both ends data in the embodiment of the present invention five；
Fig. 9 is the reverse mutual correlated series schematic diagram of circuit both ends data in the embodiment of the present invention five.
Specific implementation mode
The present invention will be further described in detail below in conjunction with the accompanying drawings：
Embodiment one,
Fig. 1 is a kind of high-adaptability fault positioning method for transmission line flow chart provided by the invention, as shown in Figure 1, this
Inventing a kind of high-adaptability fault positioning method for transmission line provided includes：
The three-phase current after transmission line malfunction is extracted, carrying out modal transformation to the three-phase current obtains line modulus；
Multi-scale transform is carried out to the line modulus, positions the failure initial time at the transmission line of electricity both ends；
According to the failure initial time calculation delay at the transmission line of electricity both ends；
The positive cross-correlation sequence at the transmission line of electricity both ends and reversed is determined according to wavelet transform dimension and the time delay
Cross-correlation sequence；
Fault point back wave is identified according to the positive cross-correlation sequence and reverse mutual correlated series at the transmission line of electricity both ends
And/or busbar back wave；
Fault localization calculating is carried out according to the fault point back wave and/or busbar back wave identified.
Three-phase current after the extraction transmission line malfunction carries out modal transformation to the three-phase current and obtains line mould
Amount, including：
Modal transformation is carried out as the following formula to the three-phase current, is transformed to line modulus：
In formula, Ia、Ib、IcFor three-phase current；IαFor the α modulus of three-phase current, IβFor the β modulus of three-phase current, I0It is three
Zero modulus of phase current.
It is described that multi-scale transform is carried out to line modulus, the failure initial time at the transmission line of electricity both ends is positioned, including：
Multi-scale transform is carried out to the line modulus with binary wavelet, according to wavelet transform dimension, passes through wavelet transformation
Line modulus is decomposed into approximation coefficient and detail coefficients；
Using the amplitude mean value of the detail coefficients of the wavelet transformation as standard, the frequency of access line both ends transient state amplitude maximum
Band is analyzed；
The wavelet transformation detail coefficients of the selected frequency band of extraction, and calculate modulus maximum sequence；
According to the modulus maximum sequence positioning failure initial time.
Time delay is calculated as follows according to the failure initial time at the transmission line of electricity both ends：
Δ t=tm-tn
In formula, Δ t is time delay, tmAnd tnFor the failure initial time at circuit both ends.
Determine the positive cross-correlation sequence at the transmission line of electricity both ends as the following formula according to wavelet transform dimension and the time delay：
W (t)=xm(t)×xn(t+Δt)
In formula, w (t) is the positive cross-correlation sequence at transmission line of electricity both ends, xm(t) and xn(t) it is respectively transmission line of electricity both ends
Wavelet conversion coefficient, Δ t be time delay.
The reverse mutual correlated series at the transmission line of electricity both ends are determined according to wavelet transform dimension and the time delay, including：
The analysis data that the either end in the transmission line of electricity both ends is calculated as follows carry out transposition：
y(t0+ i)=x (t0+ DL-i)
Wherein, DL=L/v × 1/fs
In formula, y (t) is the transmission line of electricity single ended data after transposition, and x (t) is the transmission line of electricity single ended data before transposition, t0
For wave head initial time, i=1,2...DL, DL is time window, and L is total track length, and v is velocity of wave, fsFor sample frequency；
The transmission line of electricity of transmission line of electricity single ended data and non-transposition after transposition is substituted into following formula to end data, is obtained described
The reverse mutual correlated series at transmission line of electricity both ends：
W (t) '=xm(t)×xn(t+Δt)
In formula, w (t) ' is the reverse mutual correlated series at transmission line of electricity both ends, xm(t) and xn(t) it is respectively transmission line of electricity two
The wavelet conversion coefficient at end, Δ t are time delay.
It is described to identify that fault point is anti-according to the positive cross-correlation sequence and reverse mutual correlated series at the transmission line of electricity both ends
Ejected wave and/or busbar back wave, including：
Include according to positive cross-correlation recognition sequence fault point back wave and busbar back wave, condition：
1) fault point back wave and busbar back wave are cathode sexual reflex；
2) the amplitude A of fault point back wave and busbar back wavefMeet 0.01 × A≤Af≤0.25×β1×β2× A,
In, A is initial wave head amplitude, β1、β2The respectively busbar reflectance factor at circuit both ends；With
3) fault point back wave and peer failure point reflection wave meet following formula：
0.97×L≤(t'1-t0)×v+(t'2-t0)×v≤1.03×L
In formula, t'1For fault point back wave moment, t'2For busbar back wave moment, t0For wave head initial time, L is line
Road overall length, v are velocity of wave.
Fault point back wave and/or busbar back wave are identified with reverse mutual correlated series, including：
It only can recognize that fault point back wave or busbar back wave, condition include according to the identification of reverse mutual correlated series：
1) fault point back wave or busbar back wave are cathode sexual reflex；
2) fault point back wave or busbar reflex amplitude AfThe following conditions should be met：0.01×A≤Af≤0.25×β1×
A or 0.01 × A≤Af≤0.25×β2×A。
Fault localization calculating is carried out according to the fault point back wave and busbar back wave identified, including：
If identifying fault point back wave and busbar back wave simultaneously, fault localization calculating is carried out as the following formula：
In formula, d is the distance value that test obtains, t'1For fault point back wave moment, t'2For busbar back wave moment, t0
For wave head initial time, L is total track length；
After completing fault localization calculating, the back wave of identification is corresponded respectively to by single ended data mould according to modulus maximum sequence
Maximum determines that back wave is fault point back wave or busbar back wave.
If only can recognize that fault point back wave or busbar back wave, fault localization calculating is carried out as the following formula：
D=(t'1-t0)×v/2
After completing fault localization calculating, the back wave of identification is corresponded to by single ended data modulus maxima according to modulus maximum sequence
Value, is determined as fault point or opposite end busbar back wave.
Embodiment two,
It, can be with the present invention provides a kind of high-adaptability measuring distance of transmission line fault system based on identical inventive concept
Including：
Extraction module carries out modal transformation for extracting the three-phase current after transmission line malfunction to the three-phase current
Obtain line modulus；
Locating module, for carrying out multi-scale transform to the line modulus, when the failure at positioning transmission line of electricity both ends is initial
It carves；
Computing module, for the failure initial time calculation delay according to transmission line of electricity both ends；
Correlated series module, for determining the forward direction at transmission line of electricity both ends according to wavelet transform dimension and the time delay mutually
Close sequence and reverse mutual correlated series；
Identification module, for identifying failure according to the positive cross-correlation sequence and reverse mutual correlated series at transmission line of electricity both ends
Point reflection wave and/or busbar back wave；
Range finder module, for carrying out fault localization calculating according to the fault point back wave and/or busbar back wave that identify.
Three-phase current after the extraction transmission line malfunction carries out modal transformation to the three-phase current and obtains line mould
Amount, including：
Modal transformation is carried out to three-phase current, is transformed to line modulus；
Modal transformation is carried out as the following formula：
In formula, Ia、Ib、IcFor three-phase current；IαFor the α modulus of three-phase current, IβFor the β modulus of three-phase current, I0It is three
Zero modulus of phase current.
It is described to line modulus carry out multi-scale transform, the failure initial time at positioning line fault both ends, including：
Multi-scale transform is carried out to the line modulus with binary wavelet, according to wavelet transform dimension, passes through wavelet transformation
Line modulus is decomposed into approximation coefficient and detail coefficients；
Using the amplitude mean value of the detail coefficients of wavelet transformation as standard, the frequency band of access line both ends transient state amplitude maximum into
Row analysis；
The wavelet transformation detail coefficients of the selected frequency band of extraction, and calculate modulus maximum sequence；
According to modulus maximum sequence positioning failure initial time.
Time delay is calculated as follows according to the failure initial time at circuit both ends：
Δ t=tm-tn
In formula, tmAnd tnFor the failure initial time at circuit both ends.
The positive cross-correlation sequence at circuit both ends is obtained according to wavelet transform dimension and the time delay, including：
The positive cross-correlation sequence at circuit both ends is calculated as follows：
W (t)=xm(t)×xn(t+Δt)
In formula, w (t) is the positive cross-correlation sequence at circuit both ends, xm(t) and xn(t) be respectively circuit both ends small echo become
Coefficient is changed, Δ t is time delay.
The reverse mutual correlated series at circuit both ends are obtained according to wavelet transform dimension and the time delay, including：
The analysis data that the either end in circuit both ends is calculated as follows carry out transposition：
y(t0+ i)=x (t0+ DL-i)
Wherein, DL=L/v × 1/fs
In formula, y (t) is the single ended data after transposition, and x (t) is the single ended data before transposition, t0For wave head initial time, i
=1,2...DL, DL are time window, and L is total track length, and v is velocity of wave, fsFor sample frequency；
The circuit of single ended data and non-transposition after transposition is substituted into following formula to end data, obtains the reverse mutual at circuit both ends
Correlated series：
W (t) '=xm(t)×xn(t+Δt)
In formula, w (t) ' is the reverse mutual correlated series at circuit both ends, xm(t) and xn(t) be respectively circuit both ends small echo
Transformation coefficient, Δ t are time delay.
It is described according to the positive cross-correlation sequence and reverse mutual correlated series at circuit both ends identify fault point back wave and/or
Busbar back wave, including：
According to positive cross-correlation recognition sequence fault point back wave and busbar back wave, to meet following conditions：
Fault point back wave and busbar back wave are cathode sexual reflex；
The amplitude A of fault point back wave and busbar back wavefMeet 0.01 × A≤Af≤0.25×β1×β2× A, wherein
A is initial wave head amplitude, β1、β2The respectively busbar reflectance factor at circuit both ends；
Fault point back wave and peer failure point reflection wave meet following formula：
0.97×L≤(t'1-t0)×v+(t'2-t0)×v≤1.03×L
In formula, t'1For fault point back wave moment, t'2For busbar back wave moment, t0For wave head initial time, L is line
Road overall length, v are velocity of wave；
If going out to meet the fault point back wave and busbar back wave of above-mentioned condition according to positive cross-correlation sequence None- identified,
Then further identified with reverse mutual correlated series.
Fault point back wave and/or busbar back wave are further identified with reverse mutual correlated series, including：
Fault point back wave or busbar back wave are only can recognize that according to the identification of reverse mutual correlated series, under identification will meet
State condition：
Fault point back wave or busbar back wave are cathode sexual reflex；
Fault point back wave or busbar reflex amplitude AfThe following conditions should be met：0.01×A≤Af≤0.25×β1×A
Or 0.01 × A≤Af≤0.25×β2×A。
Fault localization calculating is carried out according to the fault point back wave and busbar back wave identified, including：
If identifying fault point back wave and busbar back wave simultaneously, fault localization calculating is carried out as the following formula：
In formula, d is the distance value that test obtains, t'1For fault point back wave moment, t'2For busbar back wave moment, t0
For wave head initial time, L is total track length；
After completing fault localization calculating, the back wave of identification is corresponded respectively to by single ended data mould according to modulus maximum sequence
Maximum determines that back wave is fault point back wave or busbar back wave.
If only can recognize that fault point back wave or busbar back wave, fault localization calculating is carried out as the following formula：
D=(t'1-t0)×v/2
After completing fault localization calculating, the back wave of identification is corresponded to by single ended data modulus maxima according to modulus maximum sequence
Value, is determined as fault point or opposite end busbar back wave.
Embodiment three,
The present invention is based on two end signal time-frequency correlations to disclose a kind of high-adaptability fault positioning method for transmission line.By
In fault point/busbar back wave and line length, time service device, stand in signal cable it is unrelated, and pass through and select transmission range
Closer back wave, which completes ranging, also the influence for reducing velocity of wave error to range accuracy, therefore, high-adaptability transmission line of electricity
Fault distance-finding method emphasis is the reliable recognition for solving fault point/busbar back wave.
Conventional single-ended traveling wave method main interference factors in back wave identification have following：
1) signal resonance, if signal resonance caused by busbar end equipment, then circuit both ends harmonic wave exists apparent
Difference, and it is orthogonal；If caused by fault point equivalent source and line parameter circuit value resonance, two end signal harmonic wave phases
Seemingly, polarity is identical, and being single-ended traveling wave method needs key problems-solving；
2) with busbar adjacent lines back wave, with busbar adjacent lines back wave, circuit both ends related to this side bus
Such interference volume in collected transient state travelling wave signal is orthogonal；
3) power electronic equipment, stand in the noise jamming that brings of the actuating of relay, it is equally, such dry in two end signal of circuit
The amount of disturbing is also orthogonal；
Relative to above-mentioned noise jamming, the transient state travelling wave that transmission line malfunction generates generates when reaching each impedance discontinuity point
Catadioptric ejected wave axial symmetry and Waveform Correlation when then having, as shown in Fig. 2.Circuit both ends (referred to as M, N-terminal), the ends M
Fault point back wave it is similar to axial symmetry when the opposite end busbar back wave of N-terminal and waveform, and the opposite end busbar back wave at the ends M with
The fault point back wave situation of N-terminal is also similar.
The present invention improves the recognition success rate of fault point/busbar back wave by the time-frequency correlation of two end data of circuit,
Line length, time service error are modified on this basis, and then improve the adaptability and precision of fault localization.
Fault positioning method for transmission line provided by the invention specifically includes following steps：
Step 1 modal transformation
When transmission line malfunction, non-faulting grade and failure level three-phase current after failure are extracted, three-phase current is converted
To be analyzed again after line modulus, to eliminate the influence that mould mixes.Using symmetrical component modal transformation, and choose α modulus conducts
Main Analysis Aerial mode component, transformation matrix are as follows：
Step 2 wavelet transformation and the positioning of failure initial time
The present invention carries out multi-scale transform using binary wavelet to the line modulus that step 1 obtains, with wavelet transformation details
The amplitude mean value of coefficient is standard, and the obvious frequency band of access line both ends transient characteristic (correspondence analysis scale) is analyzed.
Wavelet transformation includes：
Signal decomposition is that approximation coefficient (low frequency part) and detail coefficients are (high by a under given analysis scale, wavelet transformation
Frequency part), as follows.
Y (i)=l (i)+d (i) (2)
In formula 2, d (i) (i=0,1 ... it is n) detail coefficients, detail coefficients are suitable with preferable temporal resolution
It is calculated for failure initial time.The present invention identifies failure initial time using wavelet modulus maxima method, and steps are as follows：
(1) extraction wavelet transformation detail coefficients d (i) calculates modulus maximum sequence；
(2) according to modulus maximum positioning failure initial time；
Can obtain circuit both ends initial time by Wavelet Modulus Maxima Algorithm is respectively：tm、tn, adopted since circuit both ends are asynchronous
Sample, therefore, generally there are time delays between the initial time that both ends measure：Δ t=tm-tn。
Step 3 forward direction translation compensation and cross-correlation sequence calculate
When two end data of transmission line of electricity is considered as two independent signals, at given wavelet transform dimension a and time delay Δ t,
Small echo cross-correlation sequence calculation formula is as follows：
W (t)=xm(t)×xn(t+Δt) (3)
In above formula, xm(t) and respectively xn(t) it is the ends M, the wavelet conversion coefficient of N-terminal, if on the basis of M end datas, it need to be to N
End signal Wavelet transformation coefficient translation Δ t realizes delay compensation, and the positive cross-correlation sequence of two end signals is can be obtained based on formula 3
Row.
In the positive cross-correlation sequence of two end signals, since fault point back wave is positive polarity, circuit opposite end busbar is anti-
Ejected wave is negative polarity, and therefore, fault point back wave and busbar back wave are cathode sexual reflex in cross-correlation sequence
The reversed translation compensation of step 4 and cross-correlation sequence calculate
Similar with step 3, when two end data of transmission line of electricity is considered as two independent signals, same type reflects in two end signals
Reversed on axis when wave, waveform is similar, therefore, can carry out transposition to M or N-terminal single ended data：
y(t0+ i)=x (t0+ DL-i) (4)
Wherein, y (t) is single ended data after transposition, t0For wave head initial time, DL=L/v × 1/fs, i=1,2...DL,
DL is time window, and window size is by total track length L, velocity of wave v and sample frequency fsIt determines.By y (t) and circuit to end data x (t)
Substitute into formula 3, you can obtain the reverse mutual correlated series of two end datas.
Step 5 back wave identifies
Interference source in conventional single-ended traveling wave method includes：With busbar adjacent lines back wave, bus equipment and signal resonance
Interference, in cross-correlation sequence, following characteristics are presented in above-mentioned interference volume：
1) with busbar adjacent lines back wave, the interference signal caused by bus equipment resonance due to orthogonal,
It is suppressed in cross-correlation sequence；
2) in positive cross-correlation sequence, the then opposite enhancing of interference signal caused by two end signal resonance, amplitude generally compared with
Height, and due to waveform is similar, when axial symmetry, be positive sexual reflex；In reverse mutual correlated series, resonance signal is asymmetric,
Therefore, there is theoretically no interfered caused by signal resonance.
In positive cross-correlation sequence, fault point, busbar back wave are identified according to following three points：
1) fault point back wave and busbar back wave are cathode sexual reflex.
2) assume that initial wave head amplitude is A, fault point/busbar reflex amplitude should be approximately equal to A × α0×β0×β1×β2,
It is fault point refraction coefficient α to enable δ0With busbar refraction coefficient β0Product, β1、β2Respectively M, N-terminal busbar reflectance factor calculate public
Formula is as follows：
β1≈(n1- 1)/(n1+1)β2≈(n2- 1)/(n2+1) (6)
In formula, n1、n2Respectively M, N-terminal are the same as busbar adjacent lines quantity, Zc1For transmission line of electricity characteristic impedance, RgFor short circuit
Transition resistance.In physical fault analysis, β1、β2It is definite value, δ is presented as shown in Fig. 3 with transition resistance relationship, and δ is generally not
More than 0.25.In addition, to avoid white noise acoustic jamming, in physical fault analysis, a lower threshold 0.01 is set, when amplitude is small
In 0.01 × A, then it is assumed that white noise acoustic jamming.Amid all these factors, therefore, fault point/busbar reflex amplitude AfShould meet with
Lower condition：0.01×A≤Af≤0.25×β1×β2×A。
3) fault point back wave and opposite end busbar back wave meet line length-velocity of wave constraints, it is assumed that fault point, mother
The line reflection wave moment is t'1、t'2, then theoretically t'2、t'3Meet condition：(t'1-t0)×v+(t'2-t0) × v=L, but consider
There is the error no more than 3% in practice to analysis scale, velocity of wave, line error, then screening conditions are as follows：
0.97×L≤(t'1-t0)×v+(t'2-t0)×v≤1.03×L (7)
Therefore, if two back waves can be filtered out in righting reflex wave train and meet formula 7, step can be directly entered
6 calculate.If the back wave combination for meeting above-mentioned condition can not be found, further identified using reverse mutual correlated series.
Compared with positive mutually facies-suite, in reverse mutual correlated series, the fault point/busbar back wave at circuit both ends is anode
Property, and the interference of signal resonance levels off to and zero is conducive to back wave identification, and but it is more sensitive to line length, velocity of wave, when circuit is long
There are under error condition, reflex amplitude can be then substantially reduced for degree, velocity of wave.Since reversed cross-correlation only identifies fault point or busbar
Back wave therefore only need to be according to following two condition perception reflex waves：
1) fault point back wave and busbar back wave are cathode sexual reflex.
2) fault point/busbar reflex amplitude AfThe following conditions should be met：0.01×A≤Af≤0.25×β1× A or 0.01
×A≤Af≤0.25×β2×A。
On the basis of completing back wave identification again, 6 completion single end distance measurements are entered step.
Step 6 one-end fault ranging
(1) if identifying fault point and busbar back wave simultaneously by step 5, it is contemplated that fault point, busbar back wave are full
Sufficient condition：(t'1-t0)×v+(t'2-t0) × v=L therefore, can be by routine to avoid the influence of line length, velocity of wave error
One-end fault ranging formula is modified to：
Complete fault localization calculate after, then by back wave correspond to single ended data modulus maximum, be determined as fault point or right
Hold busbar back wave.
(2) if being only capable of identification fault point back wave or busbar back wave by step 5, therefore, still according to the single-ended row of tradition
Wave method is completed fault point and is calculated, and is shown below.
D=(t'1-t0)×v/2 (9)
Equally, back wave single ended data modulus maximum be need to be corresponded to, fault point or opposite end busbar back wave are determined as.
Example IV,
It is analyzed by taking the red 1 line physical fault in the ports Liaoning electric power grid 220kV as an example
The ports Liaoning electric power grid 220kV 1 line failure of pellet, total track length about 55.48km.In the secondary failure, the failure of scene operation
Range unit does not provide fail result at the first time because GPS loses star, and rear practical line walking result is：It alters an agreement apart from port city fault point
16.919km alters an agreement 38.5km away from Dandong, which is typical singlephase earth fault.
Step 1：To the three-phase current of fault moment into line modular transformation；
Step 2：Wavelet transformation and failure initial time calculate：Using selected spline function small echo and change of scale to line
Mold component is obtained after decomposing shown in wavelet conversion coefficient such as Fig. 4 (a).The present invention identifies the event at circuit both ends using modulus maximum
Hinder initial time, and extract failure initial time time delay Δ t, as shown in Fig. 4 (b), after failure initial time, becomes in small echo
There are a fairly large number of noise jamming, these clutters to be impacted to single end distance measurement in mold changing maximum.
Step 3：Positive translation compensation and cross-correlation sequence calculate：After two end data forward direction translation compensations, small echo is utilized
The both ends forward direction cross-correlation sequence of extraction is converted as shown in figure 5, being effectively suppressed relative to Fig. 4 (b) noise jamming, in Fig. 5
Encircled portion is actual fault point/circuit opposite end busbar back wave, and is cathode sexual reflex.
Step 4：Reversed translation compensation and cross-correlation sequence calculate：After the data back translation compensation of both ends, small echo is utilized
The both ends reverse mutual correlated series of extraction are converted as shown in fig. 6, encircled portion is that actual fault point, circuit opposite end busbar are anti-in figure
Ejected wave, compared with positive small echo cross-correlation sequence, fault point, circuit opposite end busbar back wave are positive sexual reflex.
Step 5：Back wave identifies：Positive cross-correlation sequence is analyzed first, according to polarity described previously, amplitude item
Part screens back wave, may recognize that two back waves.And the two back waves are not with initial traveling wave time difference：
(t'1-t0)=72, (t'2-t0)=163 are about 295m/us in conjunction with 220kV double loop velocities of wave, meet line length-velocity of wave about
Beam condition (t'1-t0)×v+(t'2-t0) × v=L can be directly entered step 6 and complete to calculate.
Step 6：Single-ended traveling wave fault location：By (t'1-t0), (t'2-t0) substitute into formula 8, d=Can be calculated fault point is：16.998km, 79m is differed with actual fault point, and entire algorithm calculated
The absolute time mark that GPS need not be used to provide in journey, can complete fault localization, and with single-ended method under both-end markers abnormal conditions
Based on, it is not influenced by Divergent line length of cable, is influenced also relative reduction by line length, structure.
Embodiment five,
It is analyzed by taking Liaoning electric power grid 220kV member center line physical faults as an example
The total track length about 79.91km of Liaoning electric power grid 220kV member center lines, the secondary failure result in about extremely because of time service device
1.5km measurement errors, practical line walking result is after failure：Fault point distance element dragon is altered an agreement 37.375km, is altered an agreement away from middle stockaded village
42.53km, the secondary failure is close to metallic short circuit failure.
Step 1：To fault moment acquisition three-phase current into line modular transformation, phase component is become using symmetry transformation matrix
It is changed to Aerial mode component.
Step 2：Wavelet transformation and failure initial time calculate：Using selected spline function small echo and change of scale to line
Mold component is obtained after decomposing shown in wavelet conversion coefficient such as Fig. 7 (a).The present invention identifies the event at circuit both ends using modulus maximum
Hinder initial time, and extract failure initial time time delay Δ t, as shown in Fig. 7 (b).
Step 3：Positive translation compensation and cross-correlation sequence calculate：After two end data forward direction translation compensations, small echo is utilized
The both ends forward direction cross-correlation sequence of extraction is converted as shown in figure 8, being effectively suppressed relative to Fig. 7 noise jamming, but due to this
For secondary failure close to metallicity failure, circuit opposite end busbar back wave is weaker, therefore, without obvious fault in positive cross-correlation sequence
Point/busbar back wave.
Step 4：Reversed translation compensation and cross-correlation sequence calculate：After the data back translation compensation of both ends, small echo is utilized
The both ends reverse mutual correlated series of extraction are converted as shown in figure 9, compared with positive cross-correlation sequence, in reverse mutual correlated series therefore
Hinder the opposite enhancing of point reflection wave, as shown in red circle, and is positive sexual reflex.
Step 5：Back wave identifies：The step is anti-to in positive cross-correlation sequence analysis, not finding qualified two
Ejected wave combines；Therefore, according to reverse mutual correlated series, amplitude maximum back wave such as Fig. 9 institutes are filtered out by polarity and amplitude principle
Show.
Step 6：Single-ended traveling wave fault location：Will likely back wave substitute into formula 9 calculate, (t'1-t0)=158, in conjunction with
220kV single loop line velocities of wave 293m/us obtains fault point and is：37.035m differs 340m with actual fault point.
It should be understood by those skilled in the art that, embodiments herein can be provided as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the application
Apply the form of example.Moreover, the application can be used in one or more wherein include computer usable program code computer
The computer program production implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.)
The form of product.
The application is with reference to method, the flow of equipment (system) and computer program product according to the embodiment of the present application
Figure and/or block diagram describe.It should be understood that can be realized by computer program instructions every first-class in flowchart and/or the block diagram
The combination of flow and/or box in journey and/or box and flowchart and/or the block diagram.These computer programs can be provided
Instruct the processor of all-purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine so that the instruction executed by computer or the processor of other programmable data processing devices is generated for real
The device for the function of being specified in present one flow of flow chart or one box of multiple flows and/or block diagram or multiple boxes.
These computer program instructions, which may also be stored in, can guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works so that instruction generation stored in the computer readable memory includes referring to
Enable the manufacture of device, the command device realize in one flow of flow chart or multiple flows and/or one box of block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device so that count
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, in computer or
The instruction executed on other programmable devices is provided for realizing in one flow of flow chart or multiple flows and/or block diagram one
The step of function of being specified in a box or multiple boxes.
Finally it should be noted that：Above example is only used to illustrate the technical scheme of the present invention rather than to its protection domain
Limitation, although the application is described in detail with reference to above-described embodiment, those of ordinary skill in the art should
Understand：Those skilled in the art read the specific implementation mode of application can still be carried out after the application various changes, modification or
Person's equivalent replacement, but these changes, modification or equivalent replacement, are applying within pending claims.
Claims (11)
1. a kind of high-adaptability fault positioning method for transmission line, which is characterized in that including：
The three-phase current after transmission line malfunction is extracted, carrying out modal transformation to the three-phase current obtains line modulus；
Multi-scale transform is carried out to the line modulus, positions the failure initial time at the transmission line of electricity both ends；
According to the failure initial time calculation delay at the transmission line of electricity both ends；
According to wavelet transform dimension and the time delay determine the transmission line of electricity both ends positive cross-correlation sequence and it is reversed mutually
Close sequence；
According to the positive cross-correlation sequence and reverse mutual correlated series at the transmission line of electricity both ends identify fault point back wave and/or
Busbar back wave；
Fault localization calculating is carried out according to the fault point back wave and/or busbar back wave identified.
2. high-adaptability fault positioning method for transmission line as described in claim 1, which is characterized in that the extraction power transmission line
Three-phase current after the failure of road carries out modal transformation to the three-phase current and obtains line modulus, including：
Modal transformation is carried out as the following formula to the three-phase current, is transformed to line modulus：
In formula, Ia、Ib、IcFor three-phase current；IαFor the α modulus of three-phase current, IβFor the β modulus of three-phase current, I0For three-phase electricity
Zero modulus of stream.
3. high-adaptability fault positioning method for transmission line as described in claim 1, which is characterized in that it is described to line modulus into
Row multi-scale transform positions the failure initial time at the transmission line of electricity both ends, including：
Multi-scale transform is carried out to the line modulus with binary wavelet, according to wavelet transform dimension, by wavelet transformation by line
Modulus is decomposed into approximation coefficient and detail coefficients；
Using the amplitude mean value of the detail coefficients of the wavelet transformation as standard, the frequency band of access line both ends transient state amplitude maximum into
Row analysis；
The wavelet transformation detail coefficients of the selected frequency band of extraction, and calculate modulus maximum sequence；
According to the modulus maximum sequence positioning failure initial time.
4. high-adaptability fault positioning method for transmission line as described in claim 1, which is characterized in that according to the power transmission line
Time delay is calculated as follows in the failure initial time at road both ends：
Δ t=tm-tn
In formula, Δ t is time delay, tmAnd tnFor the failure initial time at circuit both ends.
5. high-adaptability fault positioning method for transmission line as described in claim 1, which is characterized in that according to wavelet transformation ruler
Degree and the time delay determine the positive cross-correlation sequence at the transmission line of electricity both ends as the following formula：
W (t)=xm(t)×xn(t+Δt)
In formula, w (t) is the positive cross-correlation sequence at transmission line of electricity both ends, xm(t) and xn(t) it is respectively the small of transmission line of electricity both ends
Wave conversion coefficient, Δ t are time delay.
6. high-adaptability fault positioning method for transmission line as described in claim 1, which is characterized in that according to wavelet transformation ruler
Degree and the time delay determine the reverse mutual correlated series at the transmission line of electricity both ends, including：
The analysis data that the either end in the transmission line of electricity both ends is calculated as follows carry out transposition：
y(t0+ i)=x (t0+ DL-i)
Wherein, DL=L/v × 1/fs
In formula, y (t) is the transmission line of electricity single ended data after transposition, and x (t) is the transmission line of electricity single ended data before transposition, t0For wave
Head initial time, i=1,2...DL, DL is time window, and L is total track length, and v is velocity of wave, fsFor sample frequency；
The transmission line of electricity of transmission line of electricity single ended data and non-transposition after transposition is substituted into following formula to end data, obtains the transmission of electricity
The reverse mutual correlated series at circuit both ends：
W (t) '=xm(t)×xn(t+Δt)
In formula, w (t) ' is the reverse mutual correlated series at transmission line of electricity both ends, xm(t) and xn(t) it is respectively transmission line of electricity both ends
Wavelet conversion coefficient, Δ t are time delay.
7. high-adaptability fault positioning method for transmission line as described in claim 1, which is characterized in that described according to described defeated
Positive cross-correlation sequence and reverse mutual correlated series the identification fault point back wave and/or busbar back wave at electric line both ends, packet
It includes：
Include according to positive cross-correlation recognition sequence fault point back wave and busbar back wave, condition：
1) fault point back wave and busbar back wave are cathode sexual reflex；
2) the amplitude A of fault point back wave and busbar back wavefMeet 0.01 × A≤Af≤0.25×β1×β2× A, wherein A is
Initial wave head amplitude, β1、β2The respectively busbar reflectance factor at circuit both ends；With
3) fault point back wave and peer failure point reflection wave meet following formula：
0.97×L≤(t'1-t0)×v+(t'2-t0)×v≤1.03×L
In formula, t'1For fault point back wave moment, t'2For busbar back wave moment, t0For wave head initial time, L is that circuit is complete
Long, v is velocity of wave.
8. high-adaptability fault positioning method for transmission line as described in claim 1, which is characterized in that with reversed cross-correlation sequence
Row identification fault point back wave and/or busbar back wave, including：
It only can recognize that fault point back wave or busbar back wave, condition include according to the identification of reverse mutual correlated series：
1) fault point back wave or busbar back wave are cathode sexual reflex；
2) fault point back wave or busbar reflex amplitude AfThe following conditions should be met：0.01×A≤Af≤0.25×β1× A or
0.01×A≤Af≤0.25×β2×A。
9. high-adaptability fault positioning method for transmission line as claimed in claim 3, which is characterized in that according to the event identified
Hinder point reflection wave and busbar back wave carries out fault localization calculating, including：
If identifying fault point back wave and busbar back wave simultaneously, fault localization calculating is carried out as the following formula：
In formula, d is the distance value that test obtains, t'1For fault point back wave moment, t'2For busbar back wave moment, t0For wave
Head initial time, L is total track length；
After completing fault localization calculating, the back wave of identification is corresponded respectively to by single ended data modulus maxima according to modulus maximum sequence
Value determines that back wave is fault point back wave or busbar back wave.
10. high-adaptability fault positioning method for transmission line as claimed in claim 3, which is characterized in that if only can recognize that
Fault point back wave or busbar back wave then carry out fault localization calculating as the following formula：
D=(t1'-t0)×v/2
After completing fault localization calculating, the back wave of identification is corresponded to by single ended data modulus maximum according to modulus maximum sequence,
It is determined as fault point or opposite end busbar back wave.
11. a kind of high-adaptability measuring distance of transmission line fault system, which is characterized in that including：
Extraction module carries out modal transformation to the three-phase current and obtains for extracting the three-phase current after transmission line malfunction
Line modulus；
Locating module, for carrying out multi-scale transform, the failure initial time at positioning transmission line of electricity both ends to the line modulus；
Computing module, for the failure initial time calculation delay according to transmission line of electricity both ends；
Correlated series module, the positive cross-correlation sequence for determining transmission line of electricity both ends according to wavelet transform dimension and the time delay
Row and reverse mutual correlated series；
Identification module, for identifying that fault point is anti-according to the positive cross-correlation sequence and reverse mutual correlated series at transmission line of electricity both ends
Ejected wave and/or busbar back wave；
Range finder module, for carrying out fault localization calculating according to the fault point back wave and/or busbar back wave that identify.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810051826.9A CN108344923B (en) | 2018-01-19 | 2018-01-19 | High-adaptability power transmission line fault location method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810051826.9A CN108344923B (en) | 2018-01-19 | 2018-01-19 | High-adaptability power transmission line fault location method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108344923A true CN108344923A (en) | 2018-07-31 |
| CN108344923B CN108344923B (en) | 2024-01-16 |
Family
ID=62961130
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810051826.9A Active CN108344923B (en) | 2018-01-19 | 2018-01-19 | High-adaptability power transmission line fault location method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108344923B (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109239532A (en) * | 2018-11-06 | 2019-01-18 | 安徽康能电气有限公司 | A kind of line fault positioner based on the sampling of FPGA traveling wave |
| CN110672986A (en) * | 2019-08-27 | 2020-01-10 | 西安电子科技大学 | Cable fault positioning system for reducing positioning blind area and improving resolution |
| CN111123033A (en) * | 2019-12-23 | 2020-05-08 | 贵州电网有限责任公司 | Distribution line potential fault identification method |
| CN111123032A (en) * | 2019-12-23 | 2020-05-08 | 贵州电网有限责任公司 | Distribution lines latent fault identification system |
| CN111948493A (en) * | 2020-08-21 | 2020-11-17 | 兰州理工大学 | MMC-HVDC direct current transmission line fault positioning method |
| CN113625104A (en) * | 2021-07-25 | 2021-11-09 | 三峡大学 | Traveling wave fault positioning-oriented line length correction method |
| CN113848428A (en) * | 2021-09-29 | 2021-12-28 | 华南理工大学 | Power transmission line double-end fault distance measurement method, system, device and medium |
| CN113884821A (en) * | 2021-08-23 | 2022-01-04 | 华能国际电力江苏能源开发有限公司南通电厂 | Line fault distance determination method and equipment based on traveling wave method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102122815A (en) * | 2011-02-28 | 2011-07-13 | 清华大学 | Ultra high-speed traveling wave direction pilot protection method, device and system for high voltage transmission line |
| CN102157981A (en) * | 2011-02-28 | 2011-08-17 | 清华大学 | High-speed data acquisition and digital signal processing device |
| CN104569744A (en) * | 2014-11-26 | 2015-04-29 | 国家电网公司 | Comprehensive single-end fault positioning method applicable to power distribution network lines |
-
2018
- 2018-01-19 CN CN201810051826.9A patent/CN108344923B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102122815A (en) * | 2011-02-28 | 2011-07-13 | 清华大学 | Ultra high-speed traveling wave direction pilot protection method, device and system for high voltage transmission line |
| CN102157981A (en) * | 2011-02-28 | 2011-08-17 | 清华大学 | High-speed data acquisition and digital signal processing device |
| CN104569744A (en) * | 2014-11-26 | 2015-04-29 | 国家电网公司 | Comprehensive single-end fault positioning method applicable to power distribution network lines |
Non-Patent Citations (3)
| Title |
|---|
| FERNANDO MARINHO DE MAGALHÄES 等: "Using current traveling waves to implement directional elements in parallel lines", 《2017 WORKSHOP ON COMMUNICATION NETWORKS AND POWER SYSTEMS》, pages 1 - 5 * |
| 张帆: "基于单端暂态行波的接地故障测距与保护研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》, no. 1 * |
| 钱峰 等: "一种新型双端时域故障定位算法", 《高电压技术》, vol. 31, no. 2, pages 81 - 83 * |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109239532A (en) * | 2018-11-06 | 2019-01-18 | 安徽康能电气有限公司 | A kind of line fault positioner based on the sampling of FPGA traveling wave |
| CN110672986A (en) * | 2019-08-27 | 2020-01-10 | 西安电子科技大学 | Cable fault positioning system for reducing positioning blind area and improving resolution |
| CN110672986B (en) * | 2019-08-27 | 2021-11-05 | 西安电子科技大学 | Cable fault positioning system for reducing positioning blind area and improving resolution |
| CN111123033A (en) * | 2019-12-23 | 2020-05-08 | 贵州电网有限责任公司 | Distribution line potential fault identification method |
| CN111123032A (en) * | 2019-12-23 | 2020-05-08 | 贵州电网有限责任公司 | Distribution lines latent fault identification system |
| CN111123033B (en) * | 2019-12-23 | 2022-07-29 | 贵州电网有限责任公司 | Distribution line potential fault identification method |
| CN111948493A (en) * | 2020-08-21 | 2020-11-17 | 兰州理工大学 | MMC-HVDC direct current transmission line fault positioning method |
| CN113625104A (en) * | 2021-07-25 | 2021-11-09 | 三峡大学 | Traveling wave fault positioning-oriented line length correction method |
| CN113625104B (en) * | 2021-07-25 | 2024-03-12 | 三峡大学 | Line length correction method for traveling wave fault location |
| CN113884821A (en) * | 2021-08-23 | 2022-01-04 | 华能国际电力江苏能源开发有限公司南通电厂 | Line fault distance determination method and equipment based on traveling wave method |
| CN113848428A (en) * | 2021-09-29 | 2021-12-28 | 华南理工大学 | Power transmission line double-end fault distance measurement method, system, device and medium |
| CN113848428B (en) * | 2021-09-29 | 2022-06-14 | 华南理工大学 | Power transmission line double-end fault distance measurement method, system, device and medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108344923B (en) | 2024-01-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108344923A (en) | A kind of high-adaptability fault positioning method for transmission line and system | |
| Suonan et al. | A novel fault-location method for HVDC transmission lines | |
| CN105353268B (en) | One kind is used for the judgement of transmission line of electricity distribution traveling wave fault and localization method | |
| CN106771861B (en) | Complex electric network Fault Locating Method based on wide area traveling wave energy and time difference | |
| CN102435912B (en) | Method for positioning fault disturbance point in power grid | |
| CN108896874B (en) | A kind of fault positioning method for transmission line of end connection short-term road | |
| CN107621591B (en) | A kind of transmission line of electricity iteration distance measuring method based on zero mould traveling wave speed variation characteristic | |
| CN105676075B (en) | A kind of power distribution network traveling wave fault positioning method based on multiterminal data | |
| CN104931855B (en) | The device and method for being identified and being extracted based on powerline fault traveling wave head | |
| CN108693446B (en) | A kind of Fault Locating Method of non-synchronous sampling power grid transient state travelling wave modulus time difference | |
| CN103592575A (en) | Self-adaptation weighting data fusion fault distance measurement method based on multi-sensor system | |
| CN105738759B (en) | A kind of direct current transmission line fault localization method based on transient state recorder data | |
| CN108020754A (en) | Single-ended traveling wave fault location method based on Waveform Reconstructing | |
| CN106054169A (en) | Multi-station radar signal fusion detection method based on tracking information | |
| CN106019080B (en) | A kind of common-tower double-return DC line Single Terminal Traveling Wave Fault Location method based on energy jump along the line | |
| CN110927510A (en) | Frequency domain method for power transmission line double-end traveling wave fault location | |
| CN109375060A (en) | A kind of distribution network failure wave-form similarity calculation method | |
| Bo et al. | Positional protection of transmission systems using global positioning system | |
| CN103941159A (en) | Mixed line fault location method | |
| CN113092946A (en) | Method and device for positioning ground fault of multi-branch overhead-cable mixed line of power distribution network | |
| CN106597215B (en) | A kind of fault distance-finding method and system based on transient state travelling wave in short time-window | |
| CN109709442A (en) | Power distribution network travelling wave ranging method and system based on bus Injection Signal | |
| CN109387733A (en) | A kind of distribution circuit single-phase earth fault localization method and system | |
| Li et al. | Study on wide‐area traveling wave fault line selection and fault location algorithm | |
| CN107179473A (en) | A kind of power transmission line fault locating method |
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 |