CN101509949B - Double-end asynchronous and parameter self-adapting fault distance measuring time-domain method for direct current transmission line - Google Patents

Double-end asynchronous and parameter self-adapting fault distance measuring time-domain method for direct current transmission line Download PDF

Info

Publication number
CN101509949B
CN101509949B CN2009100380777A CN200910038077A CN101509949B CN 101509949 B CN101509949 B CN 101509949B CN 2009100380777 A CN2009100380777 A CN 2009100380777A CN 200910038077 A CN200910038077 A CN 200910038077A CN 101509949 B CN101509949 B CN 101509949B
Authority
CN
China
Prior art keywords
line
wave
fault
domain
time
Prior art date
Application number
CN2009100380777A
Other languages
Chinese (zh)
Other versions
CN101509949A (en
Inventor
王钢
梁远升
李海锋
黄敏
朱革兰
黎小林
傅闯
Original Assignee
华南理工大学
南方电网技术研究中心
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 华南理工大学, 南方电网技术研究中心 filed Critical 华南理工大学
Priority to CN2009100380777A priority Critical patent/CN101509949B/en
Publication of CN101509949A publication Critical patent/CN101509949A/en
Application granted granted Critical
Publication of CN101509949B publication Critical patent/CN101509949B/en

Links

Abstract

The invention discloses a time-domain failure distance measuring method with two asynchronous ends of a direct-current power line and self-adaptive parameters, which comprises the following steps of: setting up an asynchronous time-domain failure observation equation of the direct-current line, which includes quantities of fault distance, distributed resistances of the line, transmission wave velocity, wave impedance, asynchronous time difference of data at the two ends and the like to be observed; taking the time of initial failure traveling waves which arrive at the two ends of the line respectively as reference time of data at each end; introducing a traveling wave characteristic equation; and eliminating the quantities to be observed of asynchronous time difference in the time-domain failure observing equation, thus obtaining a novel time-domain failure distance measuring method of the direct-current line so as to realize self-adaption of the line parameter self-adapting without synchronizing with a clock. The time-domain failure distance measuring method has higher transition resistance durability and can effectively improve the accuracy and reliability of failure distance measurement at the two ends of the direct-current line.

Description

The fault distance measuring time-domain method of direct current transmission line double-end asynchronous and parameter self-adapting
Technical field
The present invention relates to the method for electric system direct current transmission line fault range finding, particularly a kind ofly need not synchronous clock and parameter adaptive direct current transmission line double-end fault distance measuring time-domain method in conjunction with the row wave property.
Background technology
The fed distance of DC transmission system is long, the probability height of line failure, and the accurate range finding of line fault is significant for the reliability of minimizing interruption maintenance time, raising DC transmission system.At present, what the direct current transmission line fault range finding was adopted is the both-end traveling wave method, promptly realizes fault localization by the accurate moment at detection failure initial line mode voltage row ripple arrival circuit two ends.But under the situation of high resistance earthing fault, the more weak situation of transient state travelling wave may appear, especially on the over distance DC power transmission line, be subjected to the influence of row wave dispersion and decay even more serious, the both-end traveling wave method often occurs than mistake because of row ripple arrival circuit two ends accurately obtain than difficulty constantly.
The accuracy that row wave-wave speed detects and the clock synchronization of two end datas also have considerable influence to the both-end travelling wave ranging.For definite method of row wave-wave speed, can revise row wave-wave speed by the historical failure data, but these class methods can not satisfy the requirement that real time environment changes, the slight change of velocity of wave all can bring than mistake.The both-end traveling wave method two ends data synchronization that places one's entire reliance upon, though the precision of synchronous clock (GPS) can reach the requirement of both-end travelling wave ranging at present,, if both-end distance measuring can be broken away from the dependence to synchronous clock, then not only can reduce cost, more can improve the reliability of Two-terminal Fault Location.
Summary of the invention
The objective of the invention is to solve the range error that circuit parameter uncertainty and both-end data asynchronism are brought in the direct current transmission line fault ranging process, the accuracy and the reliability of raising DC line fault range finding under metallic earthing and high resistance earthing fault situation, provide a kind of in conjunction with the fault traveling wave characteristic, dependence can be broken away from, and the direct current transmission line fault distance-finding method of double-end asynchronous and parameter self-adapting can be realized synchronous clock.
For reaching above-mentioned purpose, the present invention adopts following technical scheme: a kind of fault distance measuring time-domain method of direct current transmission line double-end asynchronous and parameter self-adapting, be characterized in, and comprise the steps:
The first step is with sample frequency f sGather utmost point I, utmost point II DC line two ends M, the voltage of N, current instantaneous value u MI, i NI, u MII, i MII, u NI, u NII, i NII, i NII, convert line mode voltage, the line mould electric current u at M, N two ends to m, i m, u n, i n:
u m = 2 ( u mI - u mII ) i m = 2 ( i mI - i mII ) u n = 2 ( u nI - u nII ) i n = 2 ( i nI - i nII )
Second step, window scope [t during specified data Star, t Stop], the line mode voltage travelling wave signal u when decomposing data in the window with multi-scale wavelet mAnd u nCarry out layering, and calculate the frequency content E of each layer K, in frequency domain scope 4~30kHz, choose E KPeaked frequency domain (experimental results show that the frequency range of 4~30kHz can avoid the influence of direct current and low frequency component on the one hand as the centre frequency of wave filter, can reduce the influence of line parameter circuit value frequency dependent characteristic on the other hand), to described line mode voltage, line mould current traveling wave signal u m, i m, u n, i nFilter, obtain band signal best in the fault analysis signal;
Wherein, t StartAnd t StopThe initial sum of window stops constantly when being respectively data,
t stop=t start+3l/v min+t rdn
In the formula: l is the distance of DC line two ends M, N, v MinFor the line mould of DC line is propagated the minimum value of velocity of wave, v Min=0.97~0.99c, c are the light velocity; t RdnFor calculating redundancy time, t Rdn=p Δ t Cal, p is a fault time domain observation equation number, the p>actual number for the treatment of observed quantity, Δ t CalBe interval computing time, Δ t Cal>to the voltage of utmost point I or utmost point II DC line two ends M, N, current instantaneous value sampling time interval 10 times;
E K = 1 t stop - t start ∫ t = t start t stop S K 2 ( t ) dt
In the formula: s K(t) be described line mode voltage travelling wave signal u mAnd u nThe HFS of K layer signal after wavelet decomposition; The frequency band of described K layer correspondence is: f s/ 2 K+1~f s/ 2 K, centre frequency is 3f s/ 2 K+1, bandwidth is f s/ 2 K+1
In the 3rd step, set up fault time domain observation equation:
f Location ( d , r , v , Z C , t Δ , t ) = 1 2 ( 1 + rd 4 Z C ) 2 [ u m ( t - t Δ + d v ) - i m ( t - t Δ + d v ) ( Z C + rd 4 ) ]
+ 1 2 ( 1 - rd 4 Z C ) 2 [ u m ( t - t Δ - d v ) - i m ( t - t Δ - d v ) ( Z C - rd 4 ) ]
- ( rd 4 Z C ) 2 u m ( t - t Δ ) - rd 4 ( 1 + rd 4 Z C ) ( 1 - rd 4 Z C ) i m ( t - t Δ ) - - - ( 1 )
- 1 2 ( 1 + r ( l - d ) 4 Z C ) 2 [ u n ( t + l - d v ) - i n ( t + l - d v ) ( Z C + r ( l - d ) 4 ) ]
- 1 2 ( 1 - r ( l - d ) 4 Z C ) 2 [ u n ( t - l - d v ) - i n ( t - l - d v ) ( Z C - r ( l - d ) 4 ) ]
+ ( r ( l - d ) 4 Z C ) 2 u n ( t ) + r ( l - d ) 4 ( 1 + r ( l - d ) 4 Z C ) ( 1 - r ( l - d ) 4 Z C ) i n ( t ) = 0
Wherein, l is that distance, the d of DC line two ends M, N are distance, the t of trouble spot apart from DC line M end ΔFor DC line M, asynchronous mistiming of N two ends, r are that the line mould distributed resistance of DC line, the line mould that v is DC line are propagated velocity of wave, Z CLine modular character impedance for DC line;
Wherein, t Δ=t M0-t N0
In the formula, t M0For the reference of the M end data clock of DC line constantly, t N0Be the reference moment of DC line N end data clock;
Set up row wave property equation:
d=(v twΔt f+l)/2
Δt f=T fm-T fm=t fm-t fn+t Δ
In the formula, v TwBe the capable wave-wave speed of described line mode voltage, Δ t fBe the capable ripple arrival DC line M of described line mode voltage, the mistiming at N two ends, T FmAnd T FnBe respectively fault initial line mode voltage row ripple and arrive the moment that DC line M holds and N holds; t FmAnd t FnBe respectively fault initial line mode voltage row ripple and arrive the moment at DC line M, N two ends with respect to each self-reference moment t M0, t N0Time;
The 4th step is to described line mode voltage travelling wave signal u mAnd u nAsk for modulus maximum constantly, it is arrived moment of DC line M end and N end as fault initial line mode voltage row ripple, and be made as the reference moment t of M, N two ends data clock the moment that fault initial line mode voltage row ripple arrives DC line M, N two ends respectively M0, t N0, t then Fm=t Fn=0, t Δ=Δ t f, promptly M, asynchronous mistiming of N two ends equal the mistiming that fault initial line mode voltage row ripple arrives the circuit two ends, according to row wave property equation, then:
t Δ=Δt f=(2d-l)/k vv (2)
Wherein, k vRatio for travelling wave analysis line mould velocity of wave and time-domain analysis line mould velocity of wave;
In the 5th step,, get respectively and calculate t=t constantly formula (2) substitution formula (1) 1, t 2, t 3..., t P, foundation comprises the described actual asynchronous fault Optimization of Time Domain of the both-end equation for the treatment of observed quantity:
F ( d , r , v , Z C ) = min [ Σ k = 1 P | f Location ( d , r , v , Z C , 2 d - l k v v , t k ) | ] - - - ( 3 )
Formula (3) is carried out global optimizing find the solution, can realize the fault time domain range finding of direct current transmission line double-end asynchronous and parameter self-adapting;
Wherein, calculate t constantly 1, t 2..., t k..., t pFor:
t 1 = t start + 2 l / v min t p = t stop - l / v min t 2 = t 1 + ( t p - t 1 ) / p · · · t k = t 1 + ( k - 1 ) ( t p - t 1 ) / p
In the formula, 1<k<p, k are natural number.
More specifically, the collection utmost point I in the described first step, utmost point II DC line two ends M, the voltage of N, current instantaneous value are meant with the fault traveling wave wave recording device that is installed in DC line M, N two ends and gather DC line two ends M, the voltage of N, current instantaneous value.
Sample frequency f in the described first step sBe 1MHz or 500kHz, this is a preferred version.
Wave filter in described second step is meant the Butterworth digital band-pass filter.
Described in described second step and the 5th step actual treat observed quantity be meant the trouble spot apart from the line modular character impedance of the line mould propagation velocity of wave of the line mould distributed resistance of the distance of DC line M end, DC line, DC line and DC line totally four treat observed quantity.
Described the 4th the step in to described line mode voltage travelling wave signal u mAnd u nAsk for modulus maximum constantly, be meant and adopt wavelet analysis detection signal singular point theory described line mode voltage travelling wave signal u mAnd u nAsk for modulus maximum constantly.
In described the 5th step formula (3) being carried out global optimizing finds the solution and is meant that utilizing particle swarm optimization algorithm that formula (3) is carried out global optimizing finds the solution.
The ratio k of travelling wave analysis line mould velocity of wave and time-domain analysis line mould velocity of wave in described the 4th step vDefinite method be: (1) adopts LINE CONSTANTS program, according to the actual track parameter: the solid conductor internal diameter of transmission line of electricity, solid conductor external diameter, direct current resistance, horizontal range, suspension height, sag minimum altitude, split conductor number, heading spacing and lead angle, the line mould that calculates DC line is propagated the relation curve of velocity of wave and analysis frequency; (2), obtain described the 4th step employing wavelet analysis detection signal singular point theory to described line mode voltage travelling wave signal u according to described relation curve mAnd u nAsk for the line mould of analysis frequency range pairing DC line in curve that modulus maximum constantly adopted and propagate velocity of wave And the E that chooses in described second step KThe line mould of the DC line of peaked frequency domain correspondence is propagated velocity of wave v '; (3) can get as calculated Because And the value of v ' has certain scope, therefore the k that calculates vCertain scope is also arranged; This value is surveyed used frequency range and the required frequency range of time-domain analysis, then k as long as determined capable wave-wave head for a certain transmission line of electricity vThe span of value not for some reason the barrier situation change and change, be a scope of determining, k vOptional value in its span all can reach and reduce row wave-wave speed v TwThe purpose of the influence that brings with the difference of propagating velocity of wave v.
The present invention with respect to the major advantage and the beneficial effect of prior art is: the fault distance measuring time-domain method of direct current transmission line double-end asynchronous and parameter self-adapting of the present invention, can be in conjunction with the row wave property, realization need not the asynchronous and adaptive direct current transmission line fault range finding of line parameter circuit value of both-end of synchronous clock, have higher transition resistance tolerance, can improve the accuracy and the reliability of DC line Two-terminal Fault Location effectively.
Description of drawings
Fig. 1 is the structural representation that comprises the DC transmission system of utmost point I, utmost point II DC line; Wherein, utmost point I circuit is the transmission line of electricity that breaks down, and utmost point II circuit is the circuit of normal transmission of electricity;
Fig. 2 be utilize among embodiment one and the embodiment two LINE CONSTANTS program right ± 500kV and ± 800kV DC line model calculates the line mould velocity of wave of gained and the relation curve of frequency.
Embodiment
The present invention is described in further detail below in conjunction with embodiment and accompanying drawing, but embodiments of the present invention are not limited thereto.
Embodiment one
This example employing ± 500kV DC transmission system model, its transmission line of electricity total length is 936km, and wherein utmost point I circuit is the transmission line of electricity that breaks down, and utmost point II circuit is the transmission line of electricity of normal transmission of electricity, and the transmission line of electricity two ends are made as M end and N end respectively.The transmission line of electricity two ends adopt the fault traveling wave wave recording device to gather the voltage and current of M end and N end, and the data sampling frequency is 1MHz.
The first step converts the voltage and current instantaneous value that utmost point I, utmost point II DC line M hold and N holds that collects to M, the line mode voltage at N two ends, electric current:
u m = 2 ( u mI - u mII ) i m = 2 ( i mI - i mII ) u n = 2 ( u nI - u nII ) i n = 2 ( i nI - i nII )
Second step, window scope [t during specified data Star, t Stop], begin to calculate t from fault moment StartGet 0, v MinDetermine according to 97% of the light velocity, get 2.91 * 10 5Km/s, p gets 100, computing time interval of delta t CalGet 0.03ms, greater than 10 times of the sampling interval of data, so calculate redundancy time t Rdn=p Δ t Cal=3ms, t Stop=t Start+ 3l/v Min+ t Rdn=12.65ms; Utilize multi-scale wavelet to decompose, during to data in the window signal carry out layering, sample frequency is 1MHz, 4~30KHz frequency range correspondence be the frequency domain scope of the 7th layer to the 5th layer high-frequency signal, this frequency range can be avoided the influence of direct current and low frequency component on the one hand, can reduce the influence of line parameter circuit value frequency dependent characteristic on the other hand, utilize formula
E K = 1 t stop - t start ∫ t = t start t stop S K 2 ( t ) dt
The frequency content E of evaluation between from the 7th layer to the 5th layer K(if sample frequency is 500kHz, the frequency range of 4~30KHz correspondence is the frequency domain scope of the 6th layer to the 4th layer high-frequency signal, the frequency content between needing this moment to estimate from the 6th layer to the 4th layer.) choose E KPeaked frequency domain, i.e. f s/ 2 K+1~f s/ 2 K, as the logical frequency of band of Butterworth digital band-pass filter, to described line mode voltage, line mould current traveling wave signal u m, i m, u n, i nFilter, obtain band signal best in the fault analysis signal;
In the 3rd step, set up fault time domain observation equation:
f Location ( d , r , v , Z C , t Δ , t ) = 1 2 ( 1 + rd 4 Z C ) 2 [ u m ( t - t Δ + d v ) - i m ( t - t Δ + d v ) ( Z C + rd 4 ) ]
+ 1 2 ( 1 - rd 4 Z C ) 2 [ u m ( t - t Δ - d v ) - i m ( t - t Δ - d v ) ( Z C - rd 4 ) ]
- ( rd 4 Z C ) 2 u m ( t - t Δ ) - rd 4 ( 1 + rd 4 Z C ) ( 1 - rd 4 Z C ) i m ( t - t Δ ) - - - ( 1 )
- 1 2 ( 1 + r ( l - d ) 4 Z C ) 2 [ u n ( t + l - d v ) - i n ( t + l - d v ) ( Z C + r ( l - d ) 4 ) ]
- 1 2 ( 1 - r ( l - d ) 4 Z C ) 2 [ u n ( t - l - d v ) - i n ( t - l - d v ) ( Z C - r ( l - d ) 4 ) ]
+ ( r ( l - d ) 4 Z C ) 2 u n ( t ) + r ( l - d ) 4 ( 1 + r ( l - d ) 4 Z C ) ( 1 - r ( l - d ) 4 Z C ) i n ( t ) = 0
Wherein, l is that distance, the d of DC line two ends M, N are distance, the t of trouble spot apart from DC line M end ΔFor DC line M, asynchronous mistiming of N two ends, r are that the line mould distributed resistance of DC line, the line mould that v is DC line are propagated velocity of wave, Z CLine modular character impedance for DC line;
Wherein, t Δ=t M0-t N0
In the formula, t M0For the reference of the M end data clock of DC line constantly, t N0Be the reference moment of DC line N end data clock;
Set up row wave property equation:
d=(v twΔt f+l)/2
Δt f=T fm-T fn=t fm-t fn+t Δ
In the formula, v TwBe the capable wave-wave speed of described line mode voltage, Δ t fBe the capable ripple arrival DC line M of described line mode voltage, the mistiming at N two ends, T FmAnd T FnBe respectively fault initial line mode voltage row ripple and arrive the moment that DC line M holds and N holds; t FmAnd t FnBe respectively fault initial line mode voltage row ripple and arrive the moment at DC line M, N two ends with respect to each self-reference moment t M0, t N0Time;
In the 4th step, adopt wavelet analysis detection signal singular point theory to above-mentioned line mode voltage travelling wave signal u mAnd u nAnalyze, ask for described line mode voltage travelling wave signal u mAnd u nModulus maximum constantly and is made as the moment that fault initial line mode voltage row ripple arrives DC line M, N two ends respectively reference the moment t of M, N two ends data clock M0, t N0, t then Fm=t Fn=0, t Δ=Δ t f, promptly M, asynchronous mistiming of N two ends equal the mistiming that fault initial line mode voltage row ripple arrives the circuit two ends, according to row wave property equation, then:
t Δ=Δt f=(2d-l)/k vv (2)
Wherein, k vRatio for travelling wave analysis line mould velocity of wave and time-domain analysis line mould velocity of wave, according to the actual track parameter: the solid conductor internal diameter of transmission line of electricity, the solid conductor external diameter, direct current resistance, horizontal range, suspension height, the sag minimum altitude, the split conductor number, heading spacing and lead angle, the relation curve that utilizes LINECONSTANTS program computational scheme line mould to propagate velocity of wave and analysis frequency is determined, as shown in Figure 2 ± the 500kV curve is described relation curve, in time-domain analysis, the 5th layer to the 7th layer HFS of the corresponding wavelet decomposition of the frequency range of bandpass filter, sample frequency for 1MHz, the centre frequency scope is 5k to 20kHz, and corresponding velocity of wave scope is 2.963~2.966km/s; Above-mentioned employing wavelet analysis detection signal singular point theory is to described line mode voltage travelling wave signal u mAnd u nAsk for the analysis frequency range that modulus maximum is adopted constantly, centre frequency is about 250~500kHz, and corresponding row wave-wave speed is about 2.977km/s, so for this circuit, k vSpan be: 1.004~1.0047, get 1.004 herein.
In the 5th step,, get a plurality of calculating t=t constantly with formula (2) substitution formula (1) 1, t 2, t 3..., t P, foundation comprises fault distance, distributed resistance, propagation velocity of wave and characteristic impedance totally 4 the asynchronous fault Optimization of Time Domain of both-end equations for the treatment of observed quantity:
F ( d , r , v , Z C ) = min [ Σ k = 1 P | f Location ( d , r , v , Z C , 2 d - l k v v , t k ) | ] - - - ( 3 )
Formula (3) is carried out global optimizing find the solution, can realize that line parameter circuit value is adaptive and need not the synchronous DC power transmission line time domain fault localization of GPS;
Wherein, calculate t constantly 1, t 2..., t k..., t p:
t 1 = t start + 2 l / v min = 6.4 ms t 100 = t stop - l / v min = 9.4 ms t 2 = t 1 + ( t p - t 1 ) / p = 6.43 ms t 3 = t 2 + ( t p - t 1 ) / p = 6.46 ms · · ·
To optimizing equation (3), adopt PSO algorithm (particle swarm optimization algorithm) to find the solution.
At the transmission line of electricity diverse location metallic earthing fault and high resistance earthing fault are set respectively, the range finding of various faults the results are shown in table 1, as seen, adopts method disclosed by the invention to carry out fault localization, and error is little, accuracy rate is high.
Table 1
10 9.923 2.9726 -3.0820 -0.077 10.118 2.9749 -3.0783 0.118
100 100.332 2.9700 -2.4759 0.332 100.640 2.9705 -2.434 0.640
200 199.520 2.9724 -1.8065 -0.480 200.404 2.9724 -1.8001 0.404
300 300.189 2.9730 -1.1289 0.189 300.251 2.9722 -1.1288 0.251
400 399.886 2.9728 -0.4583 -0.114 400.267 2.9714 -0.4559 0.267
536 535.951 2.9711 0.4574 -0.049 535.884 2.9717 0.4569 -0.116
636 635.770 2.9725 1.1288 -0.230 635.808 2.9728 1.1289 -0.192
736 736.460 2.9726 1.8062 0.460 735.654 2.9717 1.8014 -0.346
836 835.803 2.9724 2.4748 -0.197 835.769 2.9718 2.4751 -0.231
926 926.255 2.9723 3.0836 0.255 925.385 2.9713 3.0787 -0.615
936 935.880 2.9759 3.1445 -0.120 935.223 2.9709 3.1453 -0.777
Embodiment two
This example employing ± 800kV DC transmission system model, its transmission line of electricity total length is 1438km, and wherein utmost point I circuit is the transmission line of electricity that breaks down, and utmost point II circuit is the transmission line of electricity of normal transmission of electricity, and the circuit two ends are made as M end and N end respectively.The circuit two ends adopt the fault traveling wave wave recording device to gather the voltage and current of M end and N end, and the data sampling frequency is 1MHz.
The first step converts the voltage and current instantaneous value that utmost point I, utmost point II DC line M hold and N holds that collects to M, the line mode voltage at N two ends, electric current:
u m = 2 ( u mI - u mII ) i m = 2 ( i mI - i mII ) u n = 2 ( u nI - u nII ) i n = 2 ( i nI - i nII )
Second step, window scope [t during specified data Star, t Stop], begin to calculate t from fault moment StartGet 0, v MinDetermine according to 97% of the light velocity, get 2.91 * 10 5Km/s, p gets 100, Δ t CalGet 0.03ms, greater than 10 times of the sampling interval of data, so calculate redundancy time t Rdn=p Δ t Cal=3ms, t Stop=t Start+ 3l/v Min+ t Rdn=17.82ms; Utilize multi-scale wavelet to decompose, during to data in the window signal carry out layering, sample frequency is 1MHz, the frequency range of 4~30KHz correspondence is the frequency domain scope of the 7th layer to the 5th layer high-frequency signal, this frequency range can be avoided the influence of direct current and low frequency component on the one hand, can reduce the influence of line parameter circuit value frequency dependent characteristic on the other hand, utilize formula
E K = 1 t stop - t start ∫ t = t start t stop S K 2 ( t ) dt
The frequency content E of evaluation between from the 7th layer to the 5th layer K(if sample frequency is 500kHz, the frequency range of 4~30KHz correspondence is the frequency domain scope of the 6th layer to the 4th layer high-frequency signal, the frequency content between needing this moment to estimate from the 6th layer to the 4th layer.) choose E KPeaked frequency domain, i.e. f s/ 2 K+1~f s/ 2 K, as the logical frequency of band of Butterworth digital band-pass filter, to described line mode voltage, line mould current traveling wave signal u m, i m, u n, i nFilter, obtain band signal best in the fault analysis signal;
In the 3rd step, set up fault time domain observation equation:
f Location ( d , r , v , Z C , t Δ , t ) = 1 2 ( 1 + rd 4 Z C ) 2 [ u m ( t - t Δ + d v ) - i m ( t - t Δ + d v ) ( Z C + rd 4 ) ]
+ 1 2 ( 1 - rd 4 Z C ) 2 [ u m ( t - t Δ - d v ) - i m ( t - t Δ - d v ) ( Z C - rd 4 ) ]
- ( rd 4 Z C ) 2 u m ( t - t Δ ) - rd 4 ( 1 + rd 4 Z C ) ( 1 - rd 4 Z C ) i m ( t - t Δ ) - - - ( 1 )
- 1 2 ( 1 + r ( l - d ) 4 Z C ) 2 [ u n ( t + l - d v ) - i n ( t + l - d v ) ( Z C + r ( l - d ) 4 ) ]
- 1 2 ( 1 - r ( l - d ) 4 Z C ) 2 [ u n ( t - l - d v ) - i n ( t - l - d v ) ( Z C - r ( l - d ) 4 ) ]
+ ( r ( l - d ) 4 Z C ) 2 u n ( t ) + r ( l - d ) 4 ( 1 + r ( l - d ) 4 Z C ) ( 1 - r ( l - d ) 4 Z C ) i n ( t ) = 0
Wherein, l is that distance, the d of DC line two ends M, N are distance, the t of trouble spot apart from DC line M end ΔFor DC line M, asynchronous mistiming of N two ends, r are that the line mould distributed resistance of DC line, the line mould that v is DC line are propagated velocity of wave, Z CLine modular character impedance for DC line;
Wherein, t Δ=t M0-t N0
In the formula, t M0For the reference of the M end data clock of DC line constantly, t N0Be the reference moment of DC line N end data clock;
Set up row wave property equation:
d=(v twΔt f+l)/2
Δt f=T fm-T fn=t fm-t fn+t Δ
In the formula, v TwBe the capable wave-wave speed of described line mode voltage, Δ t fBe the capable ripple arrival DC line M of described line mode voltage, the mistiming at N two ends, T FmAnd T FnBe respectively fault initial line mode voltage row ripple and arrive the moment that DC line M holds and N holds; t FmAnd t FnBe respectively fault initial line mode voltage row ripple and arrive the moment at DC line M, N two ends with respect to each self-reference moment t M0, t N0Time;
In the 4th step, adopt wavelet analysis detection signal singular point theory to above-mentioned line mode voltage travelling wave signal u mAnd u nAnalyze, ask for described line mode voltage travelling wave signal u mAnd u nModulus maximum constantly and is made as the moment that fault initial line mode voltage row ripple arrives DC line M, N two ends respectively reference the moment t of M, N two ends data clock M0, t N0, then be t Fm=t Fn=0, t Δ=Δ t f, M, asynchronous mistiming of N two ends equal the mistiming that fault initial line mode voltage row ripple arrives the circuit two ends, according to row wave property equation, then:
t Δ=Δt f=(2d-l)/k vv (2)
Wherein, k vRatio for travelling wave analysis line mould velocity of wave and time-domain analysis line mould velocity of wave, according to the actual track parameter: the solid conductor internal diameter of transmission line of electricity, the solid conductor external diameter, direct current resistance, horizontal range, suspension height, the sag minimum altitude, the split conductor number, heading spacing and lead angle, the relation curve that utilizes LINECONSTANTS program computational scheme line mould to propagate velocity of wave and frequency is determined, as shown in Figure 2 ± the 800kV curve is described relation curve, in time-domain analysis, the 5th layer to the 7th layer HFS of the corresponding wavelet decomposition of the frequency range of bandpass filter, sample frequency for 1MHz, the centre frequency scope is 5k to 20kHz, and corresponding velocity of wave scope is 2.958~2.961km/s; Above-mentioned employing wavelet analysis detection signal singular point theory is to described line mode voltage travelling wave signal u mAnd u nAsk for the analysis frequency range that modulus maximum is adopted constantly, centre frequency is about 250~500kHz, and corresponding row wave-wave speed is about 2.966km/s, so for this circuit, k vSpan be: 1.0017~1.0027, get 1.002 herein.
In the 5th step,, get a plurality of calculating t=t constantly with formula (2) substitution formula (1) 1, t 2, t 3..., t P, foundation comprises fault distance, distributed resistance, propagation velocity of wave and characteristic impedance totally 4 the asynchronous fault Optimization of Time Domain of both-end equations for the treatment of observed quantity:
F ( d , r , v , Z C ) = min [ Σ k = 1 P | f Location ( d , r , v , Z C , 2 d - l k v v , t k ) | ] - - - ( 3 )
Formula (3) is carried out global optimizing find the solution, can realize that line parameter circuit value is adaptive and need not the synchronous DC power transmission line time domain fault localization of GPS;
Wherein, calculate t constantly 1, t 2..., t k..., t p:
t 1 = t start + 2 l / v min = 9 . 9 ms t 100 = t stop - l / v min = 12.9 ms t 2 = t 1 + ( t p - t 1 ) / p = 9 . 93 ms t 3 = t 2 + ( t p - t 1 ) / p = 9 . 96 ms · · ·
To optimizing equation (3), adopt PSO algorithm (particle swarm optimization algorithm) to find the solution.
At the circuit diverse location metallic earthing fault and high resistance earthing fault are set respectively, the range finding of various faults the results are shown in table 2, as seen, adopts method disclosed by the invention to carry out fault localization, and error is little, accuracy rate is high.
Table 2
The foregoing description is a preferred implementation of the present invention; but embodiments of the present invention are not restricted to the described embodiments; other any do not deviate from change, the modification done under spirit of the present invention and the principle, substitutes, combination, simplify; all should be the substitute mode of equivalence, be included within protection scope of the present invention.

Claims (3)

1. the fault distance measuring time-domain method of a direct current transmission line double-end asynchronous and parameter self-adapting is characterized in that, comprises the steps:
The first step is with sample frequency f sGather utmost point I, utmost point II DC line two ends M, the voltage of N, current instantaneous value u MI, i MI, u MII, i MII, u NI, i NI, u NII, i NII, convert line mode voltage, the line mould electric current u at M, N two ends to m, i m, u n, i n:
u m = 2 ( u mI - u mII ) i m = 2 ( i mI - i mII ) u n = 2 ( u nI - u nII ) i n = 2 ( i nI - i nII )
Second step, window scope [t during specified data Star, t Stop], the line mode voltage travelling wave signal u when decomposing data in the window with multi-scale wavelet mAnd u nCarry out layering, and calculate the frequency content E of each layer K, in frequency domain scope 4~30kHz, choose E KPeaked frequency domain is as the centre frequency of wave filter, to described line mode voltage, line mould current traveling wave signal u m, i m, u n, i nFilter, obtain band signal best in the fault analysis signal; Described wave filter is meant the Butterworth digital band-pass filter;
Wherein, t StartAnd t StopThe initial sum of window stops constantly when being respectively data,
t stop=t start+3l/v min+t rdn
In the formula: l is the distance of DC line two ends M, N, v MinFor the line mould of DC line is propagated the minimum value of velocity of wave, v Min=0.97~0.99c, c are the light velocity; t RdnFor calculating redundancy time, t Rdn=p Δ t Cal, p is a fault time domain observation equation number, the p>actual number for the treatment of observed quantity, Δ t CalBe interval computing time, Δ t Cal>at interval 10 times of the voltage of utmost point I or utmost point II DC line two ends M, N, current instantaneous value acquisition time; Described actual treat observed quantity be meant propagate velocity of wave and DC line apart from the line mould of the line mould distributed resistance of the distance of DC line M end, DC line, DC line in the trouble spot the impedance of line modular character totally four treat observed quantity;
E K = 1 t stop - t start ∫ t = t start t stop S K 2 ( t ) dt
In the formula: s K(t) be described line mode voltage travelling wave signal u mAnd u nThe HFS of K layer signal after wavelet decomposition; The frequency band of described K layer correspondence is: f s/ 2 K+1~f s/ 2 K, centre frequency is 3f s/ 2 K+1, bandwidth is f s/ 2 K+1
In the 3rd step, set up fault time domain observation equation:
f Location ( d , r , v , Z C , t Δ , t ) = 1 2 ( 1 + rd 4 Z C ) 2 [ u m ( t - t Δ + d v ) - i m ( t - t Δ + d v ) ( Z C + rd 4 ) ]
+ 1 2 ( 1 - rd 4 Z C ) 2 [ u m ( t - t Δ - d v ) - i m ( t - t Δ - d v ) ( Z C - rd 4 ) ]
- ( rd 4 Z C ) 2 u m ( t - t Δ ) - rd 4 ( 1 + rd 4 Z C ) ( 1 - rd 4 Z C ) i m ( t - t Δ ) - - - ( 1 )
- 1 2 ( 1 + r ( l - d ) 4 Z C ) 2 [ u n ( t + l - d v ) - i n ( t + l - d v ) ( Z C + r ( l - d ) 4 ) ]
- 1 2 ( 1 - r ( l - d ) 4 Z C ) 2 [ u n ( t - l - d v ) - i n ( t - l - d v ) ( Z C - r ( l - d ) 4 ) ]
+ ( r ( l - d ) 4 Z C ) 2 u n ( t ) + r ( l - d ) 4 ( 1 + r ( l - d ) 4 Z C ) ( 1 - r ( l - d ) 4 Z C ) i n ( t ) = 0
Wherein, l is that distance, the d of DC line two ends M, N are distance, the t of trouble spot apart from DC line M end ΔFor DC line M, asynchronous mistiming of N two ends, r are that the line mould distributed resistance of DC line, the line mould that v is DC line are propagated velocity of wave, Z CLine modular character impedance for DC line;
Wherein, t Δ=t M0-t N0
In the formula, t M0For the reference of the M end data clock of DC line constantly, t N0Be the reference moment of DC line N end data clock;
Set up row wave property equation:
d=(v twΔt f+l)/2
Δt f=T fm-T fn=t fm-t fn+t Δ
In the formula, v TwBe the capable wave-wave speed of described line mode voltage, Δ t fBe the capable ripple arrival DC line M of described line mode voltage, the mistiming at N two ends, T FmAnd T FnBe respectively fault initial line mode voltage row ripple and arrive the moment that DC line M holds and N holds; t FmAnd t FnBe respectively fault initial line mode voltage row ripple and arrive the moment at DC line M, N two ends with respect to each self-reference moment t M0, t N0Time;
In the 4th step, adopt wavelet analysis detection signal singular point theory to described line mode voltage travelling wave signal u mAnd u nAsk for modulus maximum constantly, it is arrived moment of DC line M end and N end as fault initial line mode voltage row ripple, and be made as the reference moment t of M, N two ends data clock the moment that fault initial line mode voltage row ripple arrives DC line M, N two ends respectively M0, t N0, t then Fm=t Fn=0, t Δ=Δ t f, promptly M, asynchronous mistiming of N two ends equal the mistiming that fault initial line mode voltage row ripple arrives the circuit two ends, according to row wave property equation, then:
t Δ=Δt f=(2d-l)/k vv (2)
Wherein, k vRatio for travelling wave analysis line mould velocity of wave and time-domain analysis line mould velocity of wave;
The ratio k of described travelling wave analysis line mould velocity of wave and time-domain analysis line mould velocity of wave vDefinite method be: (1) adopts LINE CONSTANTS program, according to the actual track parameter: the solid conductor internal diameter of transmission line of electricity, solid conductor external diameter, direct current resistance, horizontal range, suspension height, sag minimum altitude, split conductor number, heading spacing and lead angle, the line mould that calculates DC line is propagated the relation curve of velocity of wave and analysis frequency; (2), obtain described the 4th step employing wavelet analysis detection signal singular point theory to described line mode voltage travelling wave signal u according to described relation curve mAnd u nAsk for the line mould of analysis frequency range pairing DC line in curve that modulus maximum constantly adopted and propagate velocity of wave And the E that chooses in described second step KThe line mould of the DC line of peaked frequency domain correspondence is propagated velocity of wave v '; (3) can get as calculated
In the 5th step,, get respectively and calculate t=t constantly formula (2) substitution formula (1) 1, t 2, t 3..., t P, foundation comprises the described actual asynchronous fault Optimization of Time Domain of the both-end equation for the treatment of observed quantity:
F ( d , r , v , Z C ) = min [ Σ k = 1 P | f Location ( d , r , v , Z C , 2 d - l k v v , t k ) | ] - - - ( 3 )
Utilize particle swarm optimization algorithm that formula (3) is carried out global optimizing and find the solution, can realize the fault time domain range finding of direct current transmission line double-end asynchronous and parameter self-adapting; Described actual treat observed quantity be meant propagate velocity of wave and DC line apart from the line mould of the line mould distributed resistance of the distance of DC line M end, DC line, DC line in the trouble spot the impedance of line modular character totally four treat observed quantity;
Wherein, calculate t constantly 1, t 2..., t k..., t pFor:
t 1 = t start + 2 l / v min t p = t stop - l / v min t 2 = t 1 + ( t p - t 1 ) / p · · · t k = t 1 + ( k - 1 ) ( t p - t 1 ) / p
In the formula, 1<k<p, k are natural number.
2. according to the fault distance measuring time-domain method of the described direct current transmission line double-end asynchronous and parameter self-adapting of claim 1, it is characterized in that the collection utmost point I in the described first step, utmost point II DC line two ends M, the voltage of N, current instantaneous value are meant with the fault traveling wave wave recording device that is installed in DC line M, N two ends gathers DC line two ends M, the voltage of N, current instantaneous value.
3. according to the fault distance measuring time-domain method of the described direct current transmission line double-end asynchronous and parameter self-adapting of claim 1, it is characterized in that the sample frequency f in the described first step sBe 1MHz or 500kHz.
CN2009100380777A 2009-03-20 2009-03-20 Double-end asynchronous and parameter self-adapting fault distance measuring time-domain method for direct current transmission line CN101509949B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100380777A CN101509949B (en) 2009-03-20 2009-03-20 Double-end asynchronous and parameter self-adapting fault distance measuring time-domain method for direct current transmission line

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100380777A CN101509949B (en) 2009-03-20 2009-03-20 Double-end asynchronous and parameter self-adapting fault distance measuring time-domain method for direct current transmission line

Publications (2)

Publication Number Publication Date
CN101509949A CN101509949A (en) 2009-08-19
CN101509949B true CN101509949B (en) 2011-05-04

Family

ID=41002409

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100380777A CN101509949B (en) 2009-03-20 2009-03-20 Double-end asynchronous and parameter self-adapting fault distance measuring time-domain method for direct current transmission line

Country Status (1)

Country Link
CN (1) CN101509949B (en)

Families Citing this family (39)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102010001198A1 (en) 2010-01-26 2011-07-28 Robert Bosch GmbH, 70469 Process for the identification of consumers or producers in a pneumatic, hydraulic or electrical network
CN101867178B (en) * 2010-03-30 2012-07-04 昆明理工大学 Fault location method using three primary colours to represent travel waves of single-phase earth fault current of transmission line
CN101825678B (en) * 2010-04-23 2012-02-01 华南理工大学 Parameter adaptive failure distance measurement method for transmission line with four-circuit on single tower
CN101825676B (en) * 2010-04-27 2012-01-11 昆明理工大学 Fault ranging method by utilizing red and green colors to express traveling wave of DC power transmission circuit malfunction voltage
CN101943738B (en) * 2010-08-04 2012-05-02 清华大学 Single-phase grounding distance measuring method and distance measuring device
CN102087332B (en) * 2010-12-20 2012-10-31 中国南方电网有限责任公司超高压输电公司检修试验中心 Direct current (DC) travelling wave fault location method based on wave velocity optimization
CN102495336B (en) * 2011-12-29 2013-09-11 上海交通大学 Distributed single-phase earth fault ranging system and ranging method thereof
RU2508556C1 (en) * 2012-10-24 2014-02-27 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Иркутский государственный технический университет" (ФГБОУ ВПО "ИрГТУ") Method for determination of short-circuit place on overhead transmission line under nonsynchronised measurements on its both ends
CN103116114B (en) * 2013-01-23 2015-08-19 中国南方电网有限责任公司超高压输电公司检修试验中心 Fault distance-finding method under a kind of DC de-icing device ground wire deicing pattern and system
CN103217622B (en) * 2013-01-28 2015-07-29 长沙理工大学 Based on the distribution network fault line selection method of multiterminal voltage traveling wave
CN103412240B (en) * 2013-07-24 2016-06-22 昆明理工大学 A kind of same tower double back transmission line single-ended traveling wave fault location method being independent of wave head identification
CN103809079A (en) * 2014-02-17 2014-05-21 华北电力大学 Double-end high frequency impedance type fault ranging method suitable for direct current distribution network
CN103760468B (en) * 2014-02-18 2016-02-24 国家电网公司 Based on the double-circuit line non-synchronous sampling Two-terminal Fault Location method of Jacobi matrix
CN103837800A (en) * 2014-03-18 2014-06-04 国家电网公司 Method for locating direct-current transmission line unipolar grounding fault
CN103941151B (en) * 2014-04-03 2016-09-21 昆明理工大学 A kind of utilize voltage, the magnitude of current coordinate be independent of both-end synchronize zero, line mould time difference radiation network Fault Locating Method
CN103969554A (en) * 2014-05-30 2014-08-06 智友光电技术发展有限公司 Online failure positioning device for high voltage cable line and positioning method thereof
CN104062549B (en) * 2014-06-20 2016-09-21 国家电网公司 The distribution joint line distance-finding method that distributed constant method of identification is combined with traveling wave method
CN104360222B (en) * 2014-12-09 2017-05-17 国网宁夏电力公司 Rapid calculation method of high-voltage/extra-high voltage direct current (DC) transmission line fault voltage traveling wave
CN105738760B (en) * 2014-12-12 2019-01-11 国家电网公司 A kind of high resistive fault distance measuring method of frequency domain method in conjunction with traveling wave method
CN104597376B (en) * 2015-01-07 2018-08-03 西安理工大学 A kind of high voltage direct current transmission line fault location method considering actual measurement velocity of wave
CN104749488A (en) * 2015-03-31 2015-07-01 华南理工大学 Direct-current circuit time domain fault distance measuring method based on continuous data window
CN105067949B (en) * 2015-07-17 2016-06-08 广东电网有限责任公司电力调度控制中心 Method for synchronizing time and system based on circuit both end voltage phase contrast
CN105301438B (en) * 2015-08-06 2018-09-21 徐振宇 Fault distance computational methods suitable for transmission line of electricity high resistance failure
CN105182184B (en) * 2015-09-29 2018-10-02 昆明理工大学 A kind of distribution cable mixing direct distribution lines fault distance-finding method decomposed based on direction traveling wave along the line
CN105510770A (en) * 2015-12-04 2016-04-20 昆明理工大学 Power transmission line one-terminal fault location method based on faulty traveling wave distribution characters along the line within two successive time windows
CN105652151B (en) * 2015-12-30 2018-08-10 武汉大学 Both-end distance measuring method based on line parameter circuit value detection and the asynchronous verification of data
CN105552863B (en) * 2016-01-28 2019-04-19 中国电力科学研究院 A kind of adjoint impedance protecting method for half-wave power transmission route
CN106093708B (en) * 2016-08-02 2019-01-04 昆明理工大学 A kind of modified common-tower double-return direct current transmission line double-end fault distance-finding method of velocity of wave
CN107797025B (en) * 2016-09-06 2020-05-12 西门子公司 Fault positioning method and device for power system
CN106443348A (en) * 2016-10-11 2017-02-22 中国南方电网有限责任公司超高压输电公司检修试验中心 Power transmission line fault travelling wave fault location automatic correction method
CN106291261A (en) * 2016-10-26 2017-01-04 四川大学 The localization method of series-parallel connection one-phase earthing failure in electric distribution network
CN106908691B (en) * 2017-02-28 2019-08-13 国网四川省电力公司电力科学研究院 Direct current transmission line fault distance measuring method based on dual data source
CN106959401B (en) * 2017-03-24 2019-11-15 国网山东省电力公司 A kind of distribution travelling wave ranging method
CN107422224B (en) * 2017-04-26 2018-11-02 中国矿业大学 A kind of asynchronous fault zone localization method of power grid
CN107505538B (en) * 2017-09-22 2018-09-14 中国矿业大学 The asynchronous Fault Locating Method of half-wave power transmission circuit based on line mode voltage traveling wave amplitude attenuation characteristic
CN107607836B (en) * 2017-09-22 2019-02-19 中国矿业大学 A kind of electric network fault area positioning method based on multidimensional information
CN107817402B (en) * 2017-10-27 2021-04-06 国网四川省电力公司电力科学研究院 Direct-current transmission line fault direction identification method based on measured wave impedance
CN107861024B (en) * 2017-10-29 2020-02-21 天津大学 Voltage traveling wave differential protection time domain calculation method for flexible direct current transmission line
CN109270406B (en) * 2018-11-02 2021-01-22 广东电网有限责任公司 Method for calculating transition resistance based on double-end asynchronous electrical quantity

Also Published As

Publication number Publication date
CN101509949A (en) 2009-08-19

Similar Documents

Publication Publication Date Title
Auzanneau Wire troubleshooting and diagnosis: Review and perspectives
Bahmanyar et al. A comparison framework for distribution system outage and fault location methods
US9002672B2 (en) Method and system for time synchronization of phase of signals from respective measurement devices
CN103675605B (en) A kind of power distribution network earth fault line selection method based on the correlation analysis of fault-signal transient state
Wu et al. Ultra-high-speed directional protection of transmission lines using mathematical morphology
CN101776725B (en) Fault positioning method for transmission line
CN103576002B (en) A kind of computing method of capacitive insulator arrangement dielectric loss angle
CN102495336B (en) Distributed single-phase earth fault ranging system and ranging method thereof
CN102822688B (en) Fault wave arrival determination
CN103344875B (en) Classification line selection method for single-phase earth fault of resonance earthing system
CN101291055B (en) Method for precisely marking arriving time of initial wave of fault generated traveling waves for electricity transmission line
He et al. Natural frequency-based line fault location in HVDC lines
CN104204824B (en) fault location in distribution system
CN102096021B (en) Traveling wave natural frequency-based power transmission network failure networking positioning and distance measurement method
Korkali et al. Traveling-wave-based fault-location technique for transmission grids via wide-area synchronized voltage measurements
CN101299538B (en) Cable-aerial mixed line fault travelling wave ranging method
EP3047416A1 (en) Power line parameter adjustment and fault location using traveling waves
Marx et al. Traveling wave fault location in protective relays: Design, testing, and results
CN102435912B (en) Method for positioning fault disturbance point in power grid
CN102565628B (en) Method for identifying properties of lightning faults of overhead line on basis of amplitude interval distribution of lightning current
CN102830328B (en) Distributed fault location method for T-circuit
CN105474022B (en) By calculating traveling wave arrival time fault location is carried out using traveling wave
CN102944814A (en) Power distribution network single-phase earth fault locating method based on transient state
CN102288869B (en) Single-end traveling wave fault ranging method for power transmission line
US20150081236A1 (en) Traveling wave validation using estimated fault location

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
ASS Succession or assignment of patent right

Owner name: RESEARCH INSTITUTE OF SOUTHERN POWER GRID CO., LTD

Free format text: FORMER OWNER: SOUTHERN POWER GRID TECHNOLOGY RESEARCH CENTER

Effective date: 20120308

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20120308

Address after: 510640 Tianhe District, Guangdong, No. five road, No. 381,

Co-patentee after: Research Institute of Southern Power Grid Co., Ltd.

Patentee after: South China University of Technology

Address before: 510640 Tianhe District, Guangdong, No. five road, No. 381,

Co-patentee before: Southern Power Grid Technology Research Center

Patentee before: South China University of Technology