CN105738771B - A kind of one-end fault ranging method based on the route containing TCSC of distribution character along fault traveling wave - Google Patents
A kind of one-end fault ranging method based on the route containing TCSC of distribution character along fault traveling wave Download PDFInfo
- Publication number
- CN105738771B CN105738771B CN201610200478.8A CN201610200478A CN105738771B CN 105738771 B CN105738771 B CN 105738771B CN 201610200478 A CN201610200478 A CN 201610200478A CN 105738771 B CN105738771 B CN 105738771B
- Authority
- CN
- China
- Prior art keywords
- traveling wave
- wave
- catastrophe point
- fault
- line
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/085—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution lines, e.g. overhead
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Locating Faults (AREA)
Abstract
The present invention relates to a kind of one-end fault ranging methods based on the route containing TCSC of distribution character along fault traveling wave, belong to Relay Protection Technology in Power System field.Firstly, obtaining measuring end fault current traveling wave data by the end M high speed acquisition device, and voltage traveling wave is constructed using adjacent sound circuit current traveling wave and wave impedance；Secondly, obtaining line mode voltage traveling wave and line mould current traveling wave by the phase-model transformation operation of the phase containing failure.Again, according to line mould electric current and line mode voltage, material calculation takes 0.1m along the line, calculates being distributed along the line for voltage and current traveling wave mutation using Bei Jielong transmission equation；Finally, in [t0,t0+ l/ (2v)] and [t0+l/(2v),t0+ l/v] when window in, take absolute value to integrate again and can obtain range function f to traveling wave mutationuI(x) and fuII(x), and according to the regularity of distribution is mutated along range function fault localization is realized.
Description
Technical field
The present invention relates to a kind of one-end fault ranging method based on the route containing TCSC of distribution character along fault traveling wave,
Belong to Relay Protection Technology in Power System field.
Background technique
The installing compensation equipment on extra high voltage network, can be by different compensation ways come the stable state of enhancement line
Transmission power ability and the transient stability margin for improving system, optimization Line Flow and reduction line loss etc..Due to route
In the uniformity of power transmission line all fronts impedance is just destroyed containing compensation device, therefore it is past using the single-ended impedance method ranging of power frequency quantity
It is past to obtain correct abort situation.Since compensation device response operating condition variation needs the regular hour, line flashover
After failure, it is feasible for carrying out fault localization using the fault traveling wave data in its several milliseconds.
Summary of the invention
The technical problem to be solved by the present invention is to propose a kind of route containing TCSC based on distribution character along fault traveling wave
One-end fault ranging method, to solve the above problems.
The technical scheme is that a kind of one-end fault based on the route containing TCSC of distribution character along fault traveling wave
Distance measuring method, when line failure, firstly, obtaining measuring end fault current traveling wave number by measuring end high speed acquisition device
According to, and voltage traveling wave is constructed using adjacent sound circuit current traveling wave and wave impedance；Secondly, the phase moding for passing through the phase containing failure
Operation is changed to obtain line mode voltage traveling wave and line mould current traveling wave；Again, it according to line mould electric current and line mode voltage, is walked along line computation
Length takes 0.1m, calculates being distributed along the line for voltage and current traveling wave mutation using Bei Jielong transmission equation；Finally, in [t0,t0+l/
(2v)] and [t0+l/(2v),t0+ l/v] when window in, take absolute value to integrate again and can obtain range function f to traveling wave mutationuI
(x) and fuII(x), and according to the regularity of distribution is mutated along range function fault localization is realized.
Specific steps are as follows:
The first step reads traveling wave data:
The measuring end fault current traveling wave data obtained by measuring end high speed acquisition device, and intercept the initial traveling wave of failure and arrive
L/v time window length after being reached up to preceding l/ (2v) time window length and the initial traveling wave of failure, the i.e. traveling wave of 1.5l/v time window length in total
Data；Wherein, l is the overall length of the transmission line of electricity containing TCSC；
Second step constructs voltage traveling wave using adjacent sound circuit current traveling wave and wave impedance, it may be assumed that
uM=ik×Zc (1)
U in formula (1)MFor measuring end voltage, ikFor longest sound circuit measuring end electric current, ZcFor surge impedance of a line；
Third step calculates direction traveling wave along route distribution:
The current traveling wave and voltage traveling wave obtained according to the first step and second step is calculated using Bei Jielong formula in [t0,t0+
L/v] it is distributed along time window length voltage traveling wave and current traveling wave, wherein t0At the time of surveying end for the initial traveling wave amount of reach of failure,
That is:
In formula, subscript s indicates modulus, s=1,2..., uM,sTo measure end line mode voltage, iM,sTo measure end line mould electric current,
X is the distance of the amount of leaving side, rsThe line mould resistance of unit length, Zc,sFor line mould wave impedance, vsLine mould wave velocity；
4th step calculates direct wave and backward-travelling wave:
Forward voltage traveling wave are as follows:
u+ x,s=(ux,s+Zc,six,s)/2 (4)
Backward voltage traveling wave are as follows:
u- x,s=(ux,s-Zc,six,s)/2 (5)
In formula, u+ x,sFor the direct wave at the x of distance measuring end, u- x,sFor the backward-travelling wave that distance measuring end is at x, ux,s
For the voltage traveling wave at the x of distance measuring end, ix,sFor the current traveling wave at the x of distance measuring end；
5th step, the mutation for extracting direct wave and backward-travelling wave:
Firstly, being obtained using formula (6) and (7) calculus of differences
With
For the difference result of direct wave,For the difference result of backward-travelling wave, Δ t is the sampling interval；
Secondly, calculating difference resultIn the energy S of a period of time2u(x, t), it may be assumed that
In formula (8),For the energy of direct wave whithin a period of time, in formula (9),For backward-travelling wave
Energy whithin a period of time；
6th step, building range function:
It is obtained according to formula (8) and formula (9)WithIn [t0,t0+ l/ (2v)] and [t0+l/(2v),t0+l/
V] in time window length, range function f is obtained according to formula (10)uI(x)、fuII(x), it may be assumed that
7th step determines fault distance:
By [t0,t0+ l/ (2v)] when window and [0, l/2] range in range function fuI(x) catastrophe point is denoted as catastrophe point solution
Collect fuI=[xI1,xI2... ...], by [t0+l/(2v),t0+ l/v] when window and [l/2, l] interior range function fuII(x) catastrophe point
It is denoted as catastrophe point disaggregation fuII=[xII1,xII2,……]；
It is matched according to formula (11), obtains fault distance:
By range function fu(x) number containing catastrophe point is defined as range function fuLength, compare fuIAnd fuIILength
Degree, and by the longer f of lengthuIt, then will be successively prominent with reference function containing catastrophe point in another range function as benchmark
Height is matched according to formula (11), and uses Euclidean distance metrics match error；The minimum corresponding catastrophe point of matching error is
For the catastrophe point of faults position, and if x* IPolarity be negative, then measuring end x is left in fault point* IIf x* IPolarity be positive,
Then measuring end l-x is left in fault point* I；
x* I+x* II=l x* I∈fuI,x* II∈fuII (11)
Meet wire length x more than a pair of of catastrophe point if it exists* I+x* II=l constraint then as described below realizes failure mutation
The identification of point；
First by [t0, t0+ l/ (2v)] when window in the catastrophe point solution concentration that is distributed along range function meet formula x* I+x* II
The corresponding distance of the catastrophe point of=l is denoted as x=[x1,x2,……xn]；
Secondly, successively assuming x=[x1,x2,……xn] in xk=xf(k=1,2 ... ... n), and in traveling wave observation window
[t0,t0+tk] in calculate range function fuk(x) the mutation distribution in [0, l/2] range；
If range function fuk(x) only exist a catastrophe point in range along, then the catastrophe point be faults position or
The catastrophe point of antithesis abort situation, and if the mutation point mutation polarity be negative, fault distance xf=x1If the mutation point mutation
Polarity is positive, then fault distance xf=l-x1；
If the catastrophe point is not reflection event without catastrophe point or an extra catastrophe point in range along range function
Hinder the catastrophe point of position or antithesis abort situation；
Wherein tkValue range be [x1/v,(x1+lMk1)/v], lMk1For the overall length of the most short sound circuit of measuring end.
The beneficial effects of the present invention are:
(1) it is not necessarily to consider the acting characteristic of compensation device, location algorithm is more reliable effectively；
(2) there is the high-pass filter effect in wire length dimension using Bei Jielong circuit model, so that distance measuring method is more
Have robustness and universality, it is easy to accomplish the functionization of single end distance measurement.
Detailed description of the invention
Fig. 1 is the line assumption diagram in the embodiment of the present invention 1, embodiment 2 and embodiment 3, total track length 200km；
Fig. 2 (a) is measuring end (end M) range function in window is grown when [t0, t0+l/ (2v)] in the embodiment of the present invention 1
Mutation distribution within the scope of wire length；
Fig. 2 (b) is measuring end (end M) ranging in window is grown when [t0+l/ (2v), t0+l/v] in the embodiment of the present invention 1
Mutation distribution of the function within the scope of wire length.
Fig. 3 (a) is measuring end (end M) range function edge in window is grown when [t0, t0+x1/v] in the embodiment of the present invention 1
Mutation distribution within the scope of wire length；
Fig. 3 (b) is measuring end (end M) range function edge in window is grown when [t0, t0+x2/v] in the embodiment of the present invention 1
Mutation distribution within the scope of wire length.
Fig. 4 (a) is measuring end (end M) range function in window is grown when [t0, t0+l/ (2v)] in the embodiment of the present invention 2
Mutation distribution within the scope of wire length；
Fig. 4 (b) is measuring end (end M) ranging in window is grown when [t0+l/ (2v), t0+l/v] in the embodiment of the present invention 2
Mutation distribution of the function within the scope of wire length.
Fig. 5 (a) is measuring end (end M) range function edge in window is grown when [t0, t0+x1/v] in the embodiment of the present invention 2
Mutation distribution within the scope of wire length；
Fig. 5 (b) is measuring end (end M) range function edge in window is grown when [t0, t0+x2/v] in the embodiment of the present invention 2
Mutation distribution within the scope of wire length.
Specific embodiment
With reference to the accompanying drawings and detailed description, the invention will be further described.
A kind of one-end fault ranging method based on the route containing TCSC of distribution character along fault traveling wave, when route occurs
When failure, firstly, obtaining measuring end fault current traveling wave data by measuring end high speed acquisition device, and adjacent sound circuit is utilized
Current traveling wave and wave impedance construct voltage traveling wave；Secondly, obtaining line mode voltage by the phase-model transformation operation of the phase containing failure
Traveling wave and line mould current traveling wave；Again, according to line mould electric current and line mode voltage, material calculation takes 0.1m along the line, using Bei Jielong
Transmission equation calculates being distributed along the line for voltage and current traveling wave mutation；Finally, in [t0,t0+ l/ (2v)] and [t0+l/(2v),t0+
L/v] when window in, take absolute value to integrate again and can obtain range function f to traveling wave mutationuI(x) and fuII(x), and according to ranging
It is mutated the regularity of distribution along function and realizes fault localization.
Specific steps are as follows:
The first step reads traveling wave data:
The measuring end fault current traveling wave data obtained by measuring end high speed acquisition device, and intercept the initial traveling wave of failure and arrive
L/v time window length after being reached up to preceding l/ (2v) time window length and the initial traveling wave of failure, the i.e. traveling wave of 1.5l/v time window length in total
Data；Wherein, l is the overall length of the transmission line of electricity containing TCSC；
Second step constructs voltage traveling wave using adjacent sound circuit current traveling wave and wave impedance, it may be assumed that
uM=ik×Zc (1)
U in formula (1)MFor measuring end voltage, ikFor longest sound circuit measuring end electric current, ZcFor surge impedance of a line；
Third step calculates direction traveling wave along route distribution:
The current traveling wave and voltage traveling wave obtained according to the first step and second step is calculated using Bei Jielong formula in [t0,t0+
L/v] it is distributed along time window length voltage traveling wave and current traveling wave, wherein t0At the time of surveying end for the initial traveling wave amount of reach of failure,
That is:
In formula, subscript s indicates modulus, s=1,2..., uM,sTo measure end line mode voltage, iM,sTo measure end line mould electric current,
X is the distance of the amount of leaving side, rsThe line mould resistance of unit length, Zc,sFor line mould wave impedance, vsLine mould wave velocity；
4th step calculates direct wave and backward-travelling wave:
Forward voltage traveling wave are as follows:
u+ x,s=(ux,s+Zc,six,s)/2 (4)
Backward voltage traveling wave are as follows:
u- x,s=(ux,s-Zc,six,s)/2 (5)
In formula, u+ x,sFor the direct wave at the x of distance measuring end, u- x,sFor the backward-travelling wave that distance measuring end is at x, ux,s
For the voltage traveling wave at the x of distance measuring end, ix,sFor the current traveling wave at the x of distance measuring end；
5th step, the mutation for extracting direct wave and backward-travelling wave:
Firstly, being obtained using formula (6) and (7) calculus of differences
With
For the difference result of direct wave,For the difference result of backward-travelling wave, Δ t is the sampling interval；
Secondly, calculating difference resultIn the energy S of a period of time2u(x, t), it may be assumed that
In formula (8),For the energy of direct wave whithin a period of time, in formula (9),For backward-travelling wave
Energy whithin a period of time；
6th step, building range function:
It is obtained according to formula (8) and formula (9)WithIn [t0,t0+ l/ (2v)] and [t0+l/(2v),t0+l/
V] in time window length, range function f is obtained according to formula (10)uI(x)、fuII(x), it may be assumed that
7th step determines fault distance:
By [t0,t0+ l/ (2v)] when window and [0, l/2] range in range function fuI(x) catastrophe point is denoted as catastrophe point solution
Collect fuI=[xI1,xI2... ...], by [t0+l/(2v),t0+ l/v] when window and [l/2, l] interior range function fuII(x) catastrophe point
It is denoted as catastrophe point disaggregation fuII=[xII1,xII2,……]；
It is matched according to formula (11), obtains fault distance:
By range function fu(x) number containing catastrophe point is defined as range function fuLength, compare fuIAnd fuIILength
Degree, and by the longer f of lengthuIt, then will be successively prominent with reference function containing catastrophe point in another range function as benchmark
Height is matched according to formula (11), and uses Euclidean distance metrics match error；The minimum corresponding catastrophe point of matching error is
For the catastrophe point of faults position, and if x* IPolarity be negative, then measuring end x is left in fault point* IIf x* IPolarity be positive,
Then measuring end l-x is left in fault point* I；
x* I+x* II=l x* I∈fuI,x* II∈fuII (11)
Meet wire length x more than a pair of of catastrophe point if it exists* I+x* II=l constraint then as described below realizes failure mutation
The identification of point；
First by [t0, t0+ l/ (2v)] when window in the catastrophe point solution concentration that is distributed along range function meet formula x* I+x* II
The corresponding distance of the catastrophe point of=l is denoted as x=[x1,x2,……xn]；
Secondly, successively assuming x=[x1,x2,……xn] in xk=xf(k=1,2 ... ... n), and in traveling wave observation window
[t0,t0+tk] in calculate range function fuk(x) the mutation distribution in [0, l/2] range；
If range function fuk(x) only exist a catastrophe point in range along, then the catastrophe point be faults position or
The catastrophe point of antithesis abort situation, and if the mutation point mutation polarity be negative, fault distance xf=x1If the mutation point mutation
Polarity is positive, then fault distance xf=l-x1；
If the catastrophe point is not reflection event without catastrophe point or an extra catastrophe point in range along range function
Hinder the catastrophe point of position or antithesis abort situation；
Wherein tkValue range be [x1/v,(x1+lMk1)/v], lMk1For the overall length of the most short sound circuit of measuring end.
The measuring end can also be expressed as the end M according to Fig. 1.
Embodiment 1: by taking transmission line of electricity shown in FIG. 1 as an example, ground connection occurs on the left of the TCSC element at distance M end 60km therefore
Barrier.
According to step 1 in specification, the traveling wave data of 1.5l/v time window length are got in measuring end (end M)；According to step
Rapid two, u is obtained using adjacent sound circuit current traveling wave and wave impedance to construct voltage traveling waveM=ik×Zc；According to step 3,
It calculates voltage traveling wave and current traveling wave and is distributed u along routex,s(x, t) and ix,s(x,t)；According to step 4 calculating direct wave and instead
To traveling wave u+ x,sAnd u- x,s；According to step 5, the mutation for extracting direct wave and anti-traveling wave is calculatedWithAnd
EnergyWithAccording to step 6, range function is constructed.Obtain range function fuI(x)、fuII(x).According to step
Rapid seven, obtain fault distance.fuI=[29.8 60.4] km, fuII=[109.2 138.9 168.5] km, it is known that, xI1+xII3=
29.8+168.5=198.3 ≈ l, xI2+xII2=60.4+138.9=199.3 ≈ l.It similarly there are and meet more than a pair of of catastrophe point
Formula x* I+x* II=l, therefore according to above-mentioned steps, obtain x=[29.8 60.4].It is now assumed that xf=x1=29.8km, and in traveling wave
Window [t when observation0,t0+x1/ v] in calculate range function fu(x), range function only one catastrophe point along wire length, it is false
It is set as standing, fault distance xf=x1=59.7km.It is next assumed that xf=x2=60.4km, the window [t in traveling wave observation0,t0+x2/
V]] in calculate range function fu(x) traveling wave in completely long range is mutated distribution.Range function has multiple along wire length
Catastrophe point, it is known that xf≠x2, further illustrate, fault distance xf=x1=59.7km.
Embodiment 2: it by taking transmission line of electricity shown in FIG. 1 as an example, is grounded at distance M end 140km on the right side of the TCSC element
Failure.
According to step 1 in specification, the traveling wave data of 1.5l/v time window length are got in measuring end M；According to step
Two, u is obtained using adjacent sound circuit current traveling wave and wave impedance to construct voltage traveling waveM=ik×Zc；According to step 3, meter
It calculates voltage traveling wave and current traveling wave and is distributed u along routex,s(x, t) and ix,s(x,t)；According to step 4 calculating direct wave and reversely
Traveling wave u+ x,sAnd u- x,s；According to step 5, the mutation for extracting direct wave and anti-traveling wave is calculatedWithAnd energy
AmountWithAccording to step 6, range function is constructed.Obtain range function fuI(x)、fuII(x).According to step
Seven, obtain fault distance.According to step 7, fault distance is obtained.fuI=[29.8 60.4] km, fuII=[109.2 138.9
168.5] km, it is known that, xI1+xII3=29.8+168.5=198.3 ≈ l, xI2+xII2=60.4+138.9=199.3 ≈ l.Together
Sample exists and meets formula x more than a pair of of catastrophe point* I+x* II=l, therefore according to above-mentioned steps, obtain x=[29.8 60.4].Existing vacation
If xf=x1=29.8km, and the window [t in traveling wave observation0,t0+x1/ v] in calculate range function fu(x) in completely long range
Traveling wave be mutated distribution, range function without catastrophe point, therefore assumes invalid, x along wire lengthf≠x1.It is next assumed that xf=
x2=60.4km, the window [t in traveling wave observation0,t0+x2/ v]] in calculate range function fu(x) traveling wave in completely long range is prominent
Variation cloth, range function only one catastrophe point along wire length, and the corresponding distance of the catastrophe point are x2, it is mutated polarity
It is positive, it assumes that set up, and fault point is left at the 139.6km of the end M.
In conjunction with attached drawing, the embodiment of the present invention is explained in detail above, but the present invention is not limited to above-mentioned
Embodiment within the knowledge of a person skilled in the art can also be before not departing from present inventive concept
Put that various changes can be made.
Claims (1)
1. a kind of one-end fault ranging method based on the route containing TCSC of distribution character along fault traveling wave, it is characterised in that:
When line failure, firstly, obtaining measuring end fault current traveling wave data by measuring end high speed acquisition device, and phase is utilized
Adjacent sound circuit current traveling wave and wave impedance construct voltage traveling wave；Secondly, being obtained by the phase-model transformation operation of the phase containing failure
Line taking mode voltage traveling wave and line mould current traveling wave；Again, according to line mould electric current and line mode voltage, material calculation takes 0.1m along the line,
Being distributed along the line for voltage and current traveling wave mutation is calculated using Bei Jielong transmission equation；Finally, in [t0,t0+ l/ (2v)] and [t0+
l/(2v),t0+ l/v] when window in, take absolute value to integrate again and can obtain range function f to traveling wave mutationuI(x) and fuII(x),
And fault localization is realized according to the regularity of distribution is mutated along range function；
Specific steps are as follows:
The first step reads traveling wave data:
The measuring end fault current traveling wave data obtained by measuring end high speed acquisition device, and intercept before the initial traveling wave of failure reaches
L/v time window length after l/ (2v) time window length and the initial traveling wave of failure reach, i.e. traveling wave data of 1.5l/v time window length in total；
Wherein, l is the overall length of the transmission line of electricity containing TCSC；
Second step constructs voltage traveling wave using adjacent sound circuit current traveling wave and wave impedance, it may be assumed that
uM=ik×Zc (1)
U in formula (1)MFor measuring end voltage, ikFor longest sound circuit measuring end electric current, ZcFor surge impedance of a line；
Third step calculates direction traveling wave along route distribution:
The current traveling wave and voltage traveling wave obtained according to the first step and second step is calculated using Bei Jielong formula in [t0,t0+l/v]
It is distributed along time window length voltage traveling wave and current traveling wave, wherein t0At the time of surveying end for the initial traveling wave amount of reach of failure, it may be assumed that
In formula, subscript s indicates modulus, s=1,2..., uM,sTo measure end line mode voltage, iM,sTo measure end line mould electric current, x is
Leave the distance of measuring end, rsThe line mould resistance of unit length, Zc,sFor line mould wave impedance, vsLine mould wave velocity；
4th step calculates direct wave and backward-travelling wave:
Forward voltage traveling wave are as follows:
u+ x,s=(ux,s+Zc,six,s)/2 (4)
Backward voltage traveling wave are as follows:
u- x,s=(ux,s-Zc,six,s)/2 (5)
In formula, u+ x,sFor the direct wave at the x of distance measuring end, u- x,sFor the backward-travelling wave that distance measuring end is at x, ux,sFor away from
From the voltage traveling wave at measuring end x, ix,sFor the current traveling wave at the x of distance measuring end；
5th step, the mutation for extracting direct wave and backward-travelling wave:
Firstly, being obtained using formula (6) and (7) calculus of differences
With
For the difference result of direct wave,For the difference result of backward-travelling wave, Δ t is the sampling interval；
Secondly, calculating difference resultIn the energy S of a period of time2u(x, t), it may be assumed that
In formula (8),For the energy of direct wave whithin a period of time, in formula (9),It is backward-travelling wave one
Energy in the section time；
6th step, building range function:
It is obtained according to formula (8) and formula (9)WithIn [t0,t0+ l/ (2v)] and [t0+l/(2v),t0+ l/v] when
In window length, range function f is obtained according to formula (10)uI(x)、fuII(x), it may be assumed that
7th step determines fault distance:
By [t0,t0+ l/ (2v)] when window and [0, l/2] range in range function fuI(x) catastrophe point is denoted as catastrophe point disaggregation fuI
=[xI1, xI2... ...], by [t0+l/(2v),t0+ l/v] when window and [l/2, l] interior range function fuII(x) catastrophe point is denoted as
Catastrophe point disaggregation fuII=[xII1,xII2,……]；
It is matched according to formula (11), obtains fault distance:
By range function fu(x) number containing catastrophe point is defined as range function fuLength, compare fuIAnd fuIILength, and
By the longer f of lengthuAs benchmark, then will successively be pressed with the catastrophe point of reference function in another range function containing catastrophe point
Illuminated (11) is matched, and uses Euclidean distance metrics match error；The minimum corresponding catastrophe point of matching error is to reflect
The catastrophe point of abort situation, and if x* IPolarity be negative, then measuring end x is left in fault point* IIf x* IPolarity be positive, then failure
Point leaves measuring end l-x* I；
x* I+x* II=l x* I∈fuI,x* II∈fuII (11)
Meet wire length x more than a pair of of catastrophe point if it exists* I+x* II=l constraint then as described below realizes failure catastrophe point
Identification；
First by [t0, t0+ l/ (2v)] when window in the catastrophe point solution concentration that is distributed along range function meet formula x* I+x* II=l's
The corresponding distance of catastrophe point is denoted as x=[x1,x2,……xn]；
Secondly, successively assuming x=[x1,x2,……xn] in xk=xf(k=1,2 ... ... n), and in traveling wave observation window [t0,t0+
tk] in calculate range function fuk(x) the mutation distribution in [0, l/2] range；
If range function fuk(x) a catastrophe point is only existed in range along, then the catastrophe point is faults position or antithesis
The catastrophe point of abort situation, and if the mutation point mutation polarity be negative, fault distance xf=x1If the mutation point mutation polarity
It is positive, then fault distance xf=l-x1；
If in range without catastrophe point or an extra catastrophe point, the catastrophe point not being faults position along range function
It sets or the catastrophe point of antithesis abort situation；
Wherein tkValue range be [x1/v,(x1+lMk1)/v], lMk1For the overall length of the most short sound circuit of measuring end.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610200478.8A CN105738771B (en) | 2016-04-01 | 2016-04-01 | A kind of one-end fault ranging method based on the route containing TCSC of distribution character along fault traveling wave |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610200478.8A CN105738771B (en) | 2016-04-01 | 2016-04-01 | A kind of one-end fault ranging method based on the route containing TCSC of distribution character along fault traveling wave |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105738771A CN105738771A (en) | 2016-07-06 |
| CN105738771B true CN105738771B (en) | 2019-07-05 |
Family
ID=56252716
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610200478.8A Active CN105738771B (en) | 2016-04-01 | 2016-04-01 | A kind of one-end fault ranging method based on the route containing TCSC of distribution character along fault traveling wave |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105738771B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106646140B (en) * | 2017-01-25 | 2019-07-12 | 国网四川省电力公司电力科学研究院 | HVDC transmission line internal fault external fault recognition methods based on measurement wave impedance |
| CN110456226B (en) * | 2019-08-22 | 2021-01-29 | 清华大学 | Fault location system of power transmission line |
| CN114609479A (en) * | 2022-05-09 | 2022-06-10 | 广东电网有限责任公司珠海供电局 | Fault positioning method and device, electronic equipment and storage medium |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1141768C (en) * | 2001-12-28 | 2004-03-10 | 清华大学 | Digital travelling wave protection method for power transmission line and its relay and protection system |
| CN105116295B (en) * | 2015-09-22 | 2018-05-25 | 昆明理工大学 | It is a kind of that overhead line fault distance-finding method is directly matched somebody with somebody based on the calibration of traveling wave mutation distance |
| 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 |
| CN105372559B (en) * | 2015-12-16 | 2019-04-12 | 昆明理工大学 | A kind of limited continuation method of chain type transmission line of electricity one-end fault ranging |
| CN105403813A (en) * | 2015-12-16 | 2016-03-16 | 昆明理工大学 | Aerial line single-end travelling wave distance measuring method based on direction traveling wave decomposition and distance calibration |
| CN105403812B (en) * | 2015-12-16 | 2018-11-27 | 昆明理工大学 | It is a kind of based on decomposing along fault traveling wave and the triangle loop grid one-end fault ranging method of distance calibration |
-
2016
- 2016-04-01 CN CN201610200478.8A patent/CN105738771B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105738771A (en) | 2016-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105372559B (en) | A kind of limited continuation method of chain type transmission line of electricity one-end fault ranging | |
| CN105116295B (en) | It is a kind of that overhead line fault distance-finding method is directly matched somebody with somebody based on the calibration of traveling wave mutation distance | |
| CN105738771B (en) | A kind of one-end fault ranging method based on the route containing TCSC of distribution character along fault traveling wave | |
| CN105182184B (en) | A kind of distribution cable mixing direct distribution lines fault distance-finding method decomposed based on direction traveling wave along the line | |
| CN104391224B (en) | A kind of distribution network failure data motor synchronizing method based on instantaneous amplitude change | |
| CN104865495B (en) | Feeder line fault travelling wave positioning method and system | |
| CN102175954B (en) | Circuit inter-phase fault single-end ranging method | |
| CN103762564B (en) | A kind of transient energy direction protection method of band SSSC transmission line of electricity | |
| CN105403812B (en) | It is a kind of based on decomposing along fault traveling wave and the triangle loop grid one-end fault ranging method of distance calibration | |
| CN102967779B (en) | Identifying method of distribution parameters of transmission line | |
| CN104459458B (en) | A kind of circuit of specific short window integrated value of utilization orientation traveling wave, which closes a floodgate, arrives fault recognition method | |
| CN105790239A (en) | High-voltage direct-current power transmission line non-unit protection method based on S transformation | |
| CN105891671A (en) | STATCOM-contained line single-end traveling wave distance measurement method based on fault traveling wave line distribution characteristics | |
| CN104833900B (en) | The faulty line selection method of low current singlephase earth fault | |
| CN105842584B (en) | A kind of T link multiterminal travelling wave ranging method based on distribution character along fault traveling wave | |
| CN107271842A (en) | A kind of parallel erected on same tower double-circuit line Fault Locating Method based on positive-sequence component in the same direction | |
| CN108896869A (en) | A kind of HVDC transmission line lightning stroke interference identification method based on S-transformation | |
| CN102508086B (en) | Parameter estimation method for electric power circuit in T-type connection | |
| CN105403813A (en) | Aerial line single-end travelling wave distance measuring method based on direction traveling wave decomposition and distance calibration | |
| CN105510770A (en) | Power transmission line one-terminal fault location method based on faulty traveling wave distribution characters along the line within two successive time windows | |
| CN106019077B (en) | A kind of current mode traveling wave ranging device Optimizing method | |
| CN103762563B (en) | A kind of overlapping sequences differential direction guard method of band SSSC transmission line of electricity | |
| CN110045226A (en) | A kind of compressed sensing based electrical power distribution network fault location method | |
| CN109164319A (en) | A kind of building user power utilization abnormality determination method | |
| CN104880648A (en) | Ranging method for overhead line power distribution network single-phase earth fault |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |