CN107505538B - The asynchronous Fault Locating Method of half-wave power transmission circuit based on line mode voltage traveling wave amplitude attenuation characteristic - Google Patents
The asynchronous Fault Locating Method of half-wave power transmission circuit based on line mode voltage traveling wave amplitude attenuation characteristic Download PDFInfo
- Publication number
- CN107505538B CN107505538B CN201710866659.9A CN201710866659A CN107505538B CN 107505538 B CN107505538 B CN 107505538B CN 201710866659 A CN201710866659 A CN 201710866659A CN 107505538 B CN107505538 B CN 107505538B
- Authority
- CN
- China
- Prior art keywords
- points
- wave
- mode voltage
- traveling wave
- amplitude
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 67
- 230000005540 biological transmission Effects 0.000 title claims abstract description 33
- 238000005259 measurement Methods 0.000 claims abstract description 18
- 230000007257 malfunction Effects 0.000 claims description 7
- 230000005611 electricity Effects 0.000 claims description 5
- 230000004807 localization Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 230000008901 benefit Effects 0.000 abstract description 3
- 230000000644 propagated effect Effects 0.000 abstract description 2
- 238000009331 sowing Methods 0.000 abstract 1
- 238000000354 decomposition reaction Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 4
- 238000000605 extraction Methods 0.000 description 4
- 101100499229 Mus musculus Dhrsx gene Proteins 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000002592 echocardiography Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000013598 vector Substances 0.000 description 3
- 239000000284 extract Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000011800 void material Substances 0.000 description 2
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 239000004411 aluminium Substances 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000004088 simulation Methods 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
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 asynchronous Fault Locating Methods of half-wave power transmission circuit based on line mode voltage traveling wave amplitude attenuation characteristic, the amplitude attenuation rule of sowing time high fdrequency component is uploaded in half-wave power transmission circuit based on line mode voltage traveling wave, it is proposed that the asynchronous fault location scheme of half-wave power transmission circuit for being combined Single Terminal Traveling Wave Fault Location with traveling wave attenuation characteristic, core concept are：After failure occurs, rough estimate abort situation selects corresponding fault distance-finding method to be positioned then according to line mode voltage traveling wave amplitude attenuation degree first.The attenuation characteristic range measurement principle of amplitude and single end distance measurement principle when being propagated based on line mode voltage traveling wave are combined and carry out half-wave power transmission line fault positioning by the present invention, have given full play to the advantage of two kinds of distance measuring methods.The present invention is not necessarily to each measurement point precise synchronization, has higher practical engineering value.
Description
Technical field
The present invention relates to a kind of half-wave power transmission circuits using line mode voltage traveling wave Mintrop wave head height frequency component amplitude characteristic
Fault Locating Method.
Background technology
A kind of ac transmission mode of the half-wave power transmission half-wavelength ac transmission as extra long distance is transmitted electricity with regular alternating current
It compares, has without installing the advantages such as the strong, good economy performance of reactive-load compensation equipment, conveying capacity.Since half-wave power transmission circuit is defeated
Electrical distance is long, and traveling wave dispersion and decaying are notable after failure occurs, therefore very big using traditional travelling wave ranging method error.
Currently, transmission line travelling wave ranging mainly has single-ended method and both-end method, single-ended method to need to detect without synchronizing
Multiple wave head arrival times.Both-end method needs to synchronize, but only that detection Mintrop wave head arrival time.Either which kind of ranging side
Method is used directly in if carrying out fault location on half-wave power transmission circuit, and error is larger (maximum up to 2% or so).Therefore, it grinds
Study carefully good reliability, positioning accuracy is high, and the half-wave power transmission line fault localization method of asynchronous measurement is very necessary.
Invention content
Goal of the invention：In order to solve the problems, such as that current traveling wave method half-wave power transmission line fault positioning exists, the present invention carries
Go out a kind of asynchronous Fault Locating Method of half-wave power transmission circuit based on line mode voltage traveling wave amplitude attenuation characteristic.
Technical solution：To achieve the above object, the technical solution adopted by the present invention is：
A kind of asynchronous Fault Locating Method of half-wave power transmission circuit based on line mode voltage traveling wave amplitude attenuation characteristic, packet
Include following steps：
(1) indicate that the measurement point at half-wave power transmission circuit head and end and midpoint, F indicate failure respectively with A, B and C point
Point；After half-wave power transmission line failure, head end A points, end B points and midpoint C point primary voltage travelling wave signals are extracted respectively,
Then phase-model transformation is carried out to it, is obtained original line mode voltage travelling wave signal, is denoted as respectivelyUnder wherein
Mark indicates that A points, B points and C points, subscript (1) indicate line mould；
(2) A points are obtained, each original line mode voltage travelling wave signal that B points and C points measureDifferent data
The amplitude of the signal component of multiple and different frequency components in the Mintrop wave head signal of window length；Meanwhile seeking arbitrary frequency using following formula
Rate fnOffline mode voltage traveling wave amplitude attenuation constant αn：
Wherein Zn=Rn+jωnL and Yn=Gn+jωnC is respectively frequency fnThe impedance of line mould and admittance of lower circuit；J is void
Number unit；RnAnd GnRespectively frequency fnThe line mould resistance and conductance of lower circuit；L and C be respectively circuit power frequency line mould inductance and
Capacitance；ωn=2 π fnFor corresponding angular frequency；Subscript n represents all electrical quantity all in frequency fnLower acquirement；
(3) the Mintrop wave head signal for being 82 in signal data window length using the line mode voltage travelling wave signal that A points and B points measure
Middle frequency is ωiThe signal amplitude of=22.727kHzWithIt brings following formula into and calculates its Amplitude Ration KAB：
In formula, x is represented apart from variable, and i represents the amount in frequencies omegaiLower acquirement.
(4) judge KABSize with 1, determines fault section：If KAB>1, then failure be happened on front half section circuit, go to
Step (5)；If KAB<1, then failure be happened on second half section circuit, go to step (10)；If KAB=1, failure is happened at circuit
Midpoint, fault distance 1500km；
(5) A point original line mode voltage travelling wave signals are utilizedBy Teager energy operator methods, finds and wherein correspond to
The Teager energy for reaching A point the first two wavefronts remembers that its value is E1AAnd E2A, and bring down two formulas into and calculate traveling wave Teager
Energy attenuation factor η1AAnd η2A：
And η2A=E2A
In formula, subscript A, 1 and 2 indicate second A points, first wavefront of A points, A points wavefront respectively；
(6) η that step (5) is calculated1AAnd η2AIt brings following criterion into, determines the pinpoint method of consequent malfunction：
If (a) η1A>0.4 ‰ and η2A>1100, then the traveling wave information of proximal end A points is utilized, event is determined using single end distance measurement
Hinder position, goes to step (7)；
If (b) η1A≤ 0.4 ‰ or η2A≤ 1100, then the traveling wave information of proximal end A points and midpoint C points is utilized, using both-end
Amplitude Ration method determines abort situation, goes to step (8)；
(7) Teager energy operator methods are utilized, A point original line mode voltage travelling wave signals are demarcatedThe first two wave head reaches
Moment is respectively T1 and T2, brings following Single Terminal Traveling Wave Fault Location formula into and calculates apart from head end A point failure distances x：
Wherein, v is traveling wave speed, and value is 2.95 × 108m/s；
(8) by the line mode voltage traveling wave Mintrop wave head signal intermediate frequency rate minute of A points in step (2) and C point different length data windows
Amount falls the corresponding amplitude of all signals in 20kHz~50kHz and chooses, and solve respective frequencies component signal amplitude it
Than；For each Amplitude Ration KAC, bring following both-end Amplitude Ration ranging formula into and calculate a failure apart from head end A points
Distance x：
Wherein, L represents circuit overall length as 3000km；For frequency fiUnder line mode voltage traveling wave amplitude attenuation constant；On
Table (1) indicates that line mould, subscript i are indicated in frequency fiUnder.
(9) arithmetic mean of instantaneous value of all fault distances in step (8) is calculated, which is finally apart from head end A points
Fault distance, so far, fault location terminate；
(10) B point original line mode voltage travelling wave signals are utilizedBy Teager energy operator methods, find wherein corresponding
In the Teager energy for reaching B point the first two wavefronts, remember that its value is E1BAnd E2B, and bring down two formulas into and calculate traveling wave
Teager energy attenuation factors η1BAnd η2B：
And η2B=E2B
In formula, subscript B, 1 and 2 indicate second B points, first wavefront of B points, B points wavefront respectively.
(11) η that step (10) is calculated1BAnd η2BIt brings following criterion into, determines the pinpoint side of consequent malfunction
Method：
If (a) η1B>0.4 ‰ and η2B>1100, then the traveling wave information of proximal end B points is utilized, event is determined using single end distance measurement
Hinder position, goes to step (12)；
If (b) η1B≤ 0.4 ‰ or η2B≤ 1100, then the traveling wave information of proximal end B points and midpoint C points is utilized, using both-end
Amplitude Ration method determines abort situation, goes to step (13)；
(12) Teager energy operator methods are utilized, B point original line mode voltage travelling wave signals are demarcatedThe first two wave head reaches
Moment is respectively T1 and T2, brings following Single Terminal Traveling Wave Fault Location formula into and calculates apart from end B point failure distances x：
Wherein, v is traveling wave speed, and value is 2.95 × 108m/s；
(13) by the line mode voltage traveling wave Mintrop wave head signal intermediate frequency rate minute of B points in step (2) and C point different length data windows
Amount falls the corresponding amplitude of all signals in 20kHz~50kHz and chooses, and solve respective frequencies component signal amplitude it
Than；For each Amplitude Ration KBC, bring following both-end Amplitude Ration ranging formula into and calculate a failure apart from end B points
Distance x：
Wherein, L represents circuit overall length as 3000km；For frequency fiUnder line mode voltage traveling wave amplitude attenuation constant；On
Table (1) indicates that line mould, subscript i are indicated in frequency fiUnder.
(14) arithmetic mean of instantaneous value of all fault distances in step (13) is calculated, which is final range line mould end
The fault distance of B points；So far, fault location terminates.
In step (2), method, line line wave height frequency component magnitude extraction method, change data are determined according to data window is calculated
Each original line mode voltage travelling wave signal that window frequency sweep method obtains A points, B points and C points measureDifferent data window
The amplitude of the signal component of multiple and different frequency components in the Mintrop wave head signal of length.
Advantageous effect：The present invention is by the attenuation characteristic range measurement principle of amplitude and single-ended survey when being propagated based on line mode voltage traveling wave
It is combined away from principle and carries out half-wave power transmission line fault positioning, given full play to the advantage of two kinds of distance measuring methods.The present invention without
Each measurement point precise synchronization is needed, there is higher practical engineering value.
Description of the drawings
Fig. 1 is half-wave power transmission line fault schematic diagram；
Fig. 2 is flow chart of the present invention.
Specific implementation mode
The present invention is further described below in conjunction with the accompanying drawings.
Half-wave power transmission circuit is as shown in Figure 1, the wherein measurement at A, B and C point difference transmission line of electricity head and end and midpoint
Point, F are fault point, U1bAnd U1fThe respectively reversed and initial traveling wave of positive line mode voltage.Dotted arrow indicates initial traveling wave in figure
The direction of propagation.A kind of asynchronous failure of half-wave power transmission circuit based on line mode voltage traveling wave amplitude attenuation characteristic of the present invention is fixed
Position method, includes the following steps：
(1) after half-wave power transmission line failure, head end A points, end B points and midpoint C point primary voltages are extracted respectively
Then travelling wave signal carries out it phase-model transformation, obtain original line mode voltage travelling wave signal, be denoted as respectivelyWherein subscript indicates that A points, B points and C points, subscript (1) indicate line mould；
(2) method, line line wave height frequency component magnitude extraction method, changing data window frequency sweep method are determined according to calculating data window
Obtain A points, each original line mode voltage travelling wave signal that B points and C points measureThe head of different data window length
The amplitude of the signal component of multiple and different frequency components in wave head signal；Meanwhile seeking optional frequency f using following formulanOffline mould
Voltage traveling wave amplitude attenuation constant αn：
Wherein, Zn=Rn+jωnL and Yn=Gn+jωnC is respectively frequency fnThe impedance of line mould and admittance of lower circuit；J is void
Number unit；RnAnd GnRespectively frequency fnThe line mould resistance and conductance of lower circuit；L and C be respectively circuit power frequency line mould inductance and
Capacitance；ωn=2 π fnFor corresponding angular frequency；Subscript n represents all electrical quantity all in frequency fnLower acquirement.
(3) the Mintrop wave head signal for being 82 in signal data window length using the line mode voltage travelling wave signal that A points and B points measure
Middle frequency is ωiThe signal amplitude of=22.727kHzWithIt brings following formula into and calculates its Amplitude Ration KAB：
Wherein, x is represented apart from variable, and i represents the amount in frequencies omegaiLower acquirement；
(4) judge KABSize with 1, determines fault section：If KAB>1, then failure be happened on front half section circuit, go to
Step (5)；If KAB<1, then failure be happened on second half section circuit, go to step (10)；If KAB=1, failure is happened at circuit
Midpoint, fault distance 1500km；
(5) A point original line mode voltage travelling wave signals are utilizedBy Teager energy operator methods, finds and wherein correspond to
The Teager energy for reaching A point the first two wavefronts remembers that its value is E1AAnd E2A, and bring down two formulas into and calculate traveling wave Teager
Energy attenuation factor η1AAnd η2A：
And η2A=E2A
Wherein, subscript A, 1 and 2 indicate second A points, first wavefront of A points, A points wavefront respectively.(6) will
The η that step (5) is calculated1AAnd η2AIt brings following criterion into, determines the pinpoint method of consequent malfunction：
If (a) η1A>0.4 ‰ and η2A>1100, then the traveling wave information of proximal end A points is utilized, event is determined using single end distance measurement
Hinder position, goes to step (7)；
If (b) η1A≤ 0.4 ‰ or η2A≤ 1100, then the traveling wave information of proximal end A points and midpoint C points is utilized, using both-end
Amplitude Ration method determines abort situation, goes to step (8)；
(7) Teager energy operator methods are utilized, A point original line mode voltage travelling wave signals are demarcatedThe first two wave head reaches
Moment is respectively T1 and T2, brings following Single Terminal Traveling Wave Fault Location formula into and calculates apart from head end A point failure distances x：
Wherein v is traveling wave speed, and value is 2.95 × 108m/s；
(8) by the line mode voltage traveling wave Mintrop wave head signal intermediate frequency rate minute of A points in step (2) and C point different length data windows
Amount falls the corresponding amplitude of all signals in 20kHz~50kHz and chooses, and solve respective frequencies component signal amplitude it
Than；For each Amplitude Ration KAC, bring following both-end Amplitude Ration ranging formula into and calculate a failure apart from head end A points
Distance x：
Wherein, L represents circuit overall length as 3000km；For frequency fiUnder line mode voltage traveling wave amplitude attenuation constant；On
Table (1) indicates that line mould, subscript i are indicated in frequency fiUnder；
(9) arithmetic mean of instantaneous value of all fault distances in step (8) is calculated, which is finally apart from head end A points
Fault distance, so far, fault location terminate；
(10) B point original line mode voltage travelling wave signals are utilizedBy Teager energy operator methods, find wherein corresponding
In the Teager energy for reaching B point the first two wavefronts, remember that its value is E1BAnd E2B, and bring down two formulas into and calculate traveling wave
Teager energy attenuation factors η1BAnd η2B：
And η2B=E2B
Wherein, subscript B, 1 and 2 indicate second B points, first wavefront of B points, B points wavefront respectively；
(11) η that step (10) is calculated1BAnd η2BIt brings following criterion into, determines the pinpoint side of consequent malfunction
Method：
If (a) η1B>0.4 ‰ and η2B>1100, then the traveling wave information of proximal end B points is utilized, event is determined using single end distance measurement
Hinder position, goes to step (12)；
If (b) η1B≤ 0.4 ‰ or η2B≤ 1100, then the traveling wave information of proximal end B points and midpoint C points is utilized, using both-end
Amplitude Ration method determines abort situation, goes to step (13)；
(12) Teager energy operator methods are utilized, B point original line mode voltage travelling wave signals are demarcatedThe first two wave head reaches
Moment is respectively T1 and T2, brings following Single Terminal Traveling Wave Fault Location formula into and calculates apart from end B point failure distances x：
Wherein, v is traveling wave speed, and value is 2.95 × 108m/s；
(13) by the line mode voltage traveling wave Mintrop wave head signal intermediate frequency rate minute of B points in step (2) and C point different length data windows
Amount falls the corresponding amplitude of all signals in 20kHz~50kHz and chooses, and solve respective frequencies component signal amplitude it
Than；For each Amplitude Ration KBC, bring following both-end Amplitude Ration ranging formula into and calculate a failure apart from end B points
Distance x：
Wherein, L represents circuit overall length as 3000km；For frequency fiUnder line mode voltage traveling wave amplitude attenuation constant；On
Table (1) indicates that line mould, subscript i are indicated in frequency fiUnder.
(14) arithmetic mean of instantaneous value of all fault distances in step (13) is calculated, which is final range line mould end
The fault distance of B points；So far, fault location terminates.
Used various methods are as follows in the present invention：
1. calculating data window determines method
After half-wave power transmission line failure, the voltage traveling wave that fault point generates can reach both ends measurement along circuit
Point.Since this method calculating process is all to be directed to line mode voltage Mintrop wave head signal as object, it is therefore desirable to which it is fixed for failure to determine
Position calculates line mode voltage traveling wave Mintrop wave head signal data window, this has a major impact consequent malfunction positioning step.Determine signal
The step of data window, is in two steps：First, substantially determining traveling wave Mintrop wave head process height；Two are to determine signal data window length.
Since db6 wavelet energies are concentrated the most, and there is good effect in terms of inflection point detection, small echo in the present invention
Transformation is all made of db6 small echos as morther wavelet.Half-wave power transmission circuit is longer, and under extreme fault condition, line mode voltage traveling wave is high
There is apparent decaying in frequency component, leads to the reconstruction coefficients mistake of d1 layers of detail coefficients of wavelet decomposition of failure distal wire mode voltage traveling wave
It is small and cause edge effect more apparent, interfere inflection point detection result.In order to solve this problem, following methods are taken to solve：
Higher sample frequency (1MHz or so) is kept, 4 layers of decomposition, selection are carried out to original line mode voltage travelling wave signal using db6 small echos
Single order edge-smoothing pattern reconstructs d3Layer detail coefficients, seek corresponding modulus maximum to determine catastrophe point.This method is substantially all
It is substantially to be recognized to the catastrophe point of traveling wave Mintrop wave head lower frequency components (about 62.5kHz~125kHz), due to compared with low frequency
The decaying of rate component is slower, therefore can ensure that catastrophe point determines not by line mode voltage traveling wave high fdrequency component amplitude to a certain extent
Decay excessive interference.
After determining line line wave sign mutation point, on the basis of the catastrophe point, the point preceding m and rear n point are chosen
Data as calculating entire line mode voltage traveling wave Mintrop wave head signal data window.It is counted with MATLAB by PSCAD emulation
It calculates, m is taken as 5 and is adapted to various situations.Due to using changing data window frequency sweep method, n that can take multiple values, specific value herein
See changing data window frequency sweep method.
2. line line wave height frequency component magnitude extraction method
The present invention determines that the element, core concept are using modulus maximum：(1) according to the frequency (S-transformation extracted
Decomposition scale), determine corresponding row vector U in complex matrix obtained by S-transformationST.(2) complex vector U is soughtSTIn each element
Modulus value, maximum modulus value is to represent the amplitude of original signal respective frequencies component.
According to above-mentioned analysis, before using S-transformation extraction original signal high fdrequency component amplitude, it is also necessary to which height is extracted in determination
Frequency (or the referred to as decomposition scale S of frequency componentlevel).The frequency should select to be expert in wave attenuation coefficient-frequency characteristic
More precipitous middle section.According to this paper PSCAD emulation and MATLAB result of calculations：Line mode voltage traveling wave Mintrop wave head is believed
Signal component comparison of computational results of the frequency between 20kHz~50kHz is accurate in number, specific S-transformation decomposition scale SlevelIt takes
Value is shown in changing data window frequency sweep method.
3. changing data window frequency sweep method
The Specific Principles proposed by the present invention for becoming data frequency sweep method are as follows：
(1) determination of data window length：N values take 70,71,72 successively in 1 trifle ..., 160；
(2) S-transformation decomposition scale：S-transformation decomposition scale S in 2 trifleslevelValue is 1,2,3 successively ..., 9；
(3) data processing method：First, after obtaining corresponding measurement point original line mode voltage travelling wave signal, 1 trifle is utilized
Method determine the line mode voltage traveling wave Mintrop wave head signal of different data window length on principle (1).Then, for each
The Mintrop wave head signal of specific data window length on principle (2), extracts the frequency point under scale 1~11 respectively using S-transformation
Amplitude and the corresponding frequency of each scale are measured, the amplitude of the frequency component between 20kHz~50kHz is grown in screening, and is sought
The real part of the line line Propagation Constants of circuit, brings formula under the Amplitude Ration of two measuring point corresponding frequencies components and the frequency
Calculate abort situation.Finally, seek calculating the arithmetic mean of instantaneous value of abort situation, which is final fault distance.
4.Teager energy operator methods
Teager energy operators method proposed by the present invention is as follows：
(1) corresponding measurement point original line mode voltage travelling wave signal is obtained；
(2) db6 small echos are used, 4 layers of wavelet decomposition are carried out to obtained signal, extract d3Layer detail coefficients；
(3) using same morther wavelet to d3Layer detail coefficients carry out wavelet reconstruction, obtain reconstructing detail coefficients accordingly,
It is denoted as：Dr2=[dr1,dr2,…,drm]=dr[m], wherein m are the length of original line mode voltage travelling wave signal；
(4) the Teager energy vectors that following formula calculates the reconstruction coefficients are brought into：
Ψe[dr[m]]=dr[m]2-dr[m-1]·dr[m+1]
In vectorial Ψ 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。
Simulating, verifying
In order to verify effectiveness of the invention and reliability, half-wave power transmission circuit model is built on PSCAD/EMTDC,
As shown in Figure 1.Using the circuit model for meeting practical line mould frequency dependent character, wherein half-wave power transmission overhead line structures use triangle tower
Model, conducting wire is using 8 division steel-cored aluminium strands.Circuit head end, midpoint and end are equipped with voltage traveling wave measuring device, exist respectively
Fault simulation emulation is carried out under the influence of different fault type, fault distance, fault resstance, failure initial phase angle factor.According to this
Literary method calculates relevant parameter and fault distance using MATLAB.Fault localization error e is defined by the formula：
In above formula, XcFor the fault distance being calculated, XrFor physical fault distance, Lt=3000km is total line length.
Fault location result of calculation is as shown in table 1 below.As space is limited, when table 1 only lists front half section line fault as a result, therefore failure
Distance is the distance of distance line head end A points.Actual emulation second half section line fault location Calculation result and front half section circuit
Result be almost symmetric.In following table, RFL is physical fault distance, and FT is fault type (Ag, ABg, ABCg, AB, ABC
Respectively represent A phases be grounded, AB two phase grounds, ABC three-phase grounds, AB line to line fault, ABC three phase short circuit faults), RfFor failure electricity
Resistance, θfFor failure initial phase angle.
Fault location result of calculation in the case of 1 different faults 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 (2)
1. a kind of asynchronous Fault Locating Method of half-wave power transmission circuit based on line mode voltage traveling wave amplitude attenuation characteristic, special
Sign is：Include the following steps：
(1) indicate that the measurement point at half-wave power transmission circuit head and end and midpoint, F indicate fault point respectively with A, B and C point；Half
After wavelength transmission line of electricity breaks down, head end A points, end B points and midpoint C point primary voltage travelling wave signals are extracted respectively, then
Phase-model transformation is carried out to it, is obtained original line mode voltage travelling wave signal, is denoted as respectivelyWherein subscript table
Show that A points, B points and C points, subscript (1) indicate line mould；
(2) A points are obtained, each original line mode voltage travelling wave signal that B points and C points measureDifferent data window is long
The amplitude of the signal component of multiple and different frequency components in the Mintrop wave head signal of degree；Meanwhile seeking optional frequency f using following formulan
Offline mode voltage traveling wave amplitude attenuation constant αn：
Wherein, Zn=Rn+jωnL and Yn=Gn+jωnC is respectively frequency fnThe impedance of line mould and admittance of lower circuit；J is imaginary number list
Position；RnAnd GnRespectively frequency fnThe line mould resistance and conductance of lower circuit；L and C is respectively the power frequency line mould inductance and electricity of circuit
Hold；ωn=2 π fnFor corresponding angular frequency；Subscript n represents all electrical quantity all in frequency fnLower acquirement；
(3) the line mode voltage travelling wave signal measured using A points and B points angle in the Mintrop wave head signal that signal data window length is 82
Frequency is ωiThe signal amplitude of=22.727kHzWithIt brings following formula into and calculates its Amplitude Ration KAB：
Wherein, x is represented apart from variable, and i is indicated in angular frequencyiUnder；
(4) judge KABSize with 1, determines fault section：If KAB>1, then failure be happened on front half section circuit, go to step
(5)；If KAB<1, then failure be happened on second half section circuit, go to step (10)；If KAB=1, failure is happened at circuit midpoint；
(5) A point original line mode voltage travelling wave signals are utilizedBy Teager energy operator methods, finds wherein to correspond to and reach
The Teager energy of A point the first two wavefronts remembers that its value is E1AAnd E2A, and bring down two formulas into and calculate traveling wave Teager energy
Decay factor η1AAnd η2A：
And η2A=E2A
Wherein, subscript A, 1 and 2 indicate second A points, first wavefront of A points, A points wavefront respectively；
(6) η that step (5) is calculated1AAnd η2AIt brings following criterion into, determines the pinpoint method of consequent malfunction：
If (a) η1A>0.4 ‰ and η2A>1100, then the traveling wave information of head end A points is utilized, fault bit is determined using single end distance measurement
It sets, goes to step (7)；
If (b) η1A≤ 0.4 ‰ or η2A≤ 1100, then the traveling wave information of head end A points and midpoint C points is utilized, using both-end amplitude
Abort situation is determined than method, goes to step (8)；
(7) Teager energy operator methods are utilized, A point original line mode voltage travelling wave signals are demarcatedThe first two wave head arrival time
Respectively T1 and T2 brings following Single Terminal Traveling Wave Fault Location formula into and calculates apart from head end A point failure distances x：
Wherein v is traveling wave speed, and value is 2.95 × 108m/s；
(8) the line mode voltage traveling wave Mintrop wave head signal intermediate frequency rate component of A points in step (2) and C point different length data windows is fallen
The ratio between choose in the corresponding amplitude of all signals of 20kHz~50kHz, and solve respective frequencies component signal amplitude；It is right
In each Amplitude Ration KAC, bring following both-end Amplitude Ration ranging formula into and calculate a fault distance apart from head end A points
x：
In above formula, L represents circuit overall length；For frequency fiUnder line mode voltage traveling wave amplitude attenuation constant；Subscript (1) indicates
Line mould, subscript i are indicated in frequency fiUnder；
(9) arithmetic mean of instantaneous value of all fault distances in step (8) is calculated, which is the final failure apart from head end A points
Distance, so far, fault location terminate；
(10) B point original line mode voltage travelling wave signals are utilizedBy Teager energy operator methods, finds and wherein correspond to
Up to the Teager energy of B point the first two wavefronts, remember that its value is E1BAnd E2B, and bring down two formulas into and calculate traveling wave Teager energy
Measure decay factor η1BAnd η2B：
And η2B=E2B
In formula, subscript B, 1 and 2 indicate second B points, first wavefront of B points, B points wavefront respectively；
(11) η that step (10) is calculated1BAnd η2BIt brings following criterion into, determines the pinpoint method of consequent malfunction：
If (a) η1B>0.4 ‰ and η2B>1100, then the traveling wave information of end B points is utilized, fault bit is determined using single end distance measurement
It sets, goes to step (12)；
If (b) η1B≤ 0.4 ‰ or η2B≤ 1100, then the traveling wave information of end B points and midpoint C points is utilized, using both-end amplitude
Abort situation is determined than method, goes to step (13)；
(12) Teager energy operator methods are utilized, B point original line mode voltage travelling wave signals are demarcatedThe first two wave head arrival time
Respectively T1 and T2 brings following Single Terminal Traveling Wave Fault Location formula into and calculates apart from end B point failure distances x1：
Wherein, v is traveling wave speed, and value is 2.95 × 108m/s；
(13) the line mode voltage traveling wave Mintrop wave head signal intermediate frequency rate component of B points in step (2) and C point different length data windows is fallen
The ratio between choose in the corresponding amplitude of all signals of 20kHz~50kHz, and solve respective frequencies component signal amplitude；It is right
In each Amplitude Ration KBC, bring following both-end Amplitude Ration ranging formula into and calculate a fault distance apart from end B points
x1：
Wherein, L represents circuit overall length；For frequency fiUnder line mode voltage traveling wave amplitude attenuation constant；Subscript (1) indicates line
Mould, subscript i are indicated in frequency fiUnder；
(14) arithmetic mean of instantaneous value of all fault distances in step (13) is calculated, which is final range line mould end B points
Fault distance；So far, fault location terminates.
2. the asynchronous event of the half-wave power transmission circuit according to claim 1 based on line mode voltage traveling wave amplitude attenuation characteristic
Hinder localization method, it is characterised in that：In step (2), determine that method, line line wave height frequency component amplitude carry according to data window is calculated
Take method, each original line mode voltage travelling wave signal that changing data window frequency sweep method obtains A points, B points and C points measureThe amplitude of the signal component of multiple and different frequency components in the Mintrop wave head signal of different data window length.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710866659.9A CN107505538B (en) | 2017-09-22 | 2017-09-22 | The asynchronous Fault Locating Method of half-wave power transmission circuit based on line mode voltage traveling wave amplitude attenuation characteristic |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710866659.9A CN107505538B (en) | 2017-09-22 | 2017-09-22 | The asynchronous Fault Locating Method of half-wave power transmission circuit based on line mode voltage traveling wave amplitude attenuation characteristic |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107505538A CN107505538A (en) | 2017-12-22 |
| CN107505538B true CN107505538B (en) | 2018-09-14 |
Family
ID=60697754
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710866659.9A Active CN107505538B (en) | 2017-09-22 | 2017-09-22 | The asynchronous Fault Locating Method of half-wave power transmission circuit based on line mode voltage traveling wave amplitude attenuation characteristic |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107505538B (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108627740B (en) * | 2018-06-22 | 2019-05-14 | 中国矿业大学 | Half-wavelength power transmission line fault location method considering traveling wave speed change and arrival time compensation |
| CN108896874B (en) * | 2018-07-10 | 2019-09-24 | 中国矿业大学 | A kind of fault positioning method for transmission line of end connection short-term road |
| CN109188210A (en) * | 2018-10-23 | 2019-01-11 | 国网江苏省电力有限公司徐州供电分公司 | A kind of urban electric power cable Two-terminal Fault Location method based on VMD-Hilbert transformation |
| CN109212384B (en) * | 2018-10-31 | 2019-07-19 | 海南电网有限责任公司电力科学研究院 | A kind of power transmission line fault locating method based on traveling wave amplitude attenuation theory |
| CN110286296A (en) * | 2019-05-20 | 2019-09-27 | 昆明理工大学 | A kind of half-wave power transmission route transient longitudinal protection method based on Simulation after test principle |
| CN110426592A (en) * | 2019-08-16 | 2019-11-08 | 南京国电南自电网自动化有限公司 | Aerial and cable hybrid line Earth design method based on the traveling wave time difference |
| CN114200246B (en) * | 2021-11-29 | 2022-06-21 | 昆明理工大学 | Fault location method and system for direct-current transmission line |
| CN115267419B (en) * | 2022-06-22 | 2023-07-04 | 天津大学 | Flexible direct-current line direction pilot protection method independent of line parameters and boundary elements |
| CN117192292B (en) * | 2023-11-07 | 2024-02-06 | 昆明理工大学 | Lightning grounding electrode line fault distance measurement method and system |
| CN117723895B (en) * | 2024-02-08 | 2024-05-03 | 南方电网数字电网研究院股份有限公司 | Distribution network fault section positioning method and device based on multi-terminal amplitude ratio matrix |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101509949B (en) * | 2009-03-20 | 2011-05-04 | 华南理工大学 | Double-end asynchronous and parameter self-adapting fault distance measuring time-domain method for direct current transmission line |
| CN102520315B (en) * | 2011-12-05 | 2013-10-16 | 西南交通大学 | Fault single end positioning method of power transmission line based on traveling wave multi-scale information |
| US9465067B2 (en) * | 2013-04-08 | 2016-10-11 | ECOLE POLYTECHNIQUE FéDéRALE DE LAUSANNE | Efficient method based on the electromagnetic time reversal to locate faults in power network |
| CN106443347B (en) * | 2016-10-09 | 2018-05-15 | 中国矿业大学 | A kind of electric network fault Section Location based on fault transient travelling wave attenuation components |
-
2017
- 2017-09-22 CN CN201710866659.9A patent/CN107505538B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN107505538A (en) | 2017-12-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107505538B (en) | The asynchronous Fault Locating Method of half-wave power transmission circuit based on line mode voltage traveling wave amplitude attenuation characteristic | |
| CN101718833B (en) | Method of single end distance measurement of power transmission line malfunction based on traveling wave inherent frequency extraction | |
| CN110596539B (en) | Power distribution network fault positioning method based on transient waveform correlation | |
| CN105353268B (en) | One kind is used for the judgement of transmission line of electricity distribution traveling wave fault and localization method | |
| CN102520315B (en) | Fault single end positioning method of power transmission line based on traveling wave multi-scale information | |
| CN108896874B (en) | A kind of fault positioning method for transmission line of end connection short-term road | |
| CN103323741B (en) | A kind of D molded line cable mixed line fault section compared based on false voltage initial row wave amplitude for strong fault sentences method for distinguishing | |
| CN104898021B (en) | A kind of distribution network fault line selection method based on k means cluster analyses | |
| CN104569744A (en) | Comprehensive single-end fault positioning method applicable to power distribution network lines | |
| CN109188210A (en) | A kind of urban electric power cable Two-terminal Fault Location method based on VMD-Hilbert transformation | |
| CN109001594A (en) | A kind of functional failure travelling wave positioning method | |
| CN110247420B (en) | Intelligent fault identification method for HVDC transmission line | |
| CN107632236A (en) | A kind of single outgoing-feeder line one-end fault ranging method based on the identification of opposite end bus back wave | |
| CN103018636A (en) | Cable single-ended travelling wave distance measurement method using fault signature frequency bands and TT conversion | |
| CN103941153B (en) | A kind of outgoing lines radiation network fault distance-finding method of k-NN algorithm based on waveform similarity | |
| CN105182187B (en) | Distribution feeder phase-to phase fault distance-finding method based on mapping function | |
| CN106019080A (en) | Line-side energy mutation based single-end travelling wave fault location method for double DC circuits on same tower | |
| CN106646121A (en) | Power distribution network fault traveling-wave range identification method | |
| CN108627740A (en) | Half-wavelength power transmission line fault location method considering traveling wave speed change and arrival time compensation | |
| CN103913676A (en) | Power transmission line single-ended fault location method based on variable traveling wave recognition time window | |
| CN115980514B (en) | Complex power distribution network fault positioning method based on multi-terminal traveling wave frequency matrix | |
| CN110514967A (en) | A kind of Fault Locating Method based on NARNN model prediction wavefront arrival time | |
| CN105445614B (en) | Double-end traveling wave fault positioning method and system based on wavelet analysis | |
| CN107179473B (en) | A kind of power transmission line fault locating method | |
| CN103941152B (en) | A kind of cable mixed line fault distance-finding method of k-NN algorithm based on waveform similarity |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 221008 Research Institute of China University of Mining and Technology,, Jiangsu Applicant after: China University of Mining & Technology Address before: 221008 Research Institute, China University of Mining and Technology, Xuzhou University, Jiangsu, China, Applicant before: China University of Mining & Technology |
|
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |