CN105116295B - 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 - Google Patents
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 Download PDFInfo
- Publication number
- CN105116295B CN105116295B CN201510606519.9A CN201510606519A CN105116295B CN 105116295 B CN105116295 B CN 105116295B CN 201510606519 A CN201510606519 A CN 201510606519A CN 105116295 B CN105116295 B CN 105116295B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- traveling wave
- mfrac
- 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
Landscapes
- Locating Faults (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
Overhead line fault distance-finding method is directly matched somebody with somebody based on the calibration of traveling wave mutation distance the present invention relates to a kind of, belongs to Relay Protection Technology in Power System field.When singlephase earth fault occurs for direct distribution lines, under sample rate 1MHz, fault feeder and the adjacent most long faulted phase current for perfecting feeder line initiating terminal are sampled, busbar terminal voltage traveling wave is obtained by adjacent most long feeder line initiating terminal current traveling wave and its wave impedance of perfecting, and most long perfected the current traveling wave that feeder line origin or beginning current traveling wave and fault feeder measuring end obtain using adjacent and calculated construction current traveling wave.It asks for being distributed along the voltage traveling wave of measuring end using obtained voltage traveling wave and current traveling wave and current traveling wave is distributed along the line, then Directional Decomposition is carried out to being distributed along traveling wave, obtain the direction traveling wave being distributed along the line, and extract the direct wave mutation being distributed along the line and backward-travelling wave mutation, finally by the two be multiplied then at traveling wave observe when window in integrated to construct range function realization fault localization.Theory analysis and simulation result show that this method works well.
Description
Technical field
The present invention for it is a kind of based on the calibration of traveling wave mutation distance directly with overhead line fault distance-finding method, belong to electric system after
Electric protection technical field.
Background technology
Under definite Power Network Neutral, power network topology and busbar outlet connection type, protection installation place is surveyed
The time domain traveling wave ripple of point reflects traveling wave attribute to moment, wave head (mutation) amplitude and polarity, containing fault location information, it with it is former
There is corresponding relation in barrier position, and the initial traveling wave of time-domain transient, trouble point back wave or the opposite end that can utilize single-ended observation are anti-
Ejected wave ripple arrival time difference carries out Single Terminal Traveling Wave Fault Location, and the initial traveling wave traveling wave of failure observed using both-end reaches the bilateral absolute moment
Difference, carry out both-end travelling wave ranging.The identification of trouble point back wave and examination are more complicated in Single Terminal Traveling Wave Fault Location method, to traveling wave width
Value, steepness are more demanding.Both-end travelling wave ranging does not need to recognize trouble point back wave, conveniently forms the machine of travelling wave ranging
Device is realized and automation, but both-end travelling wave ranging requirement circuit both ends cycle accurate is synchronous, requirement participates in what fault distance calculated
Line length engineering exhales title value l, and close to its " true value ", unavoidably there are objective range errors.In any case, current row
Ripple distance measuring method is based on fault traveling wave temporal signatures and in being observed, portray and wave head mark to traveling wave on time shaft mostly
Fixed and fault distance calculating.Therefore, it is badly in need of proposing a kind of new fault distance-finding method, it is effective from fault traveling wave wave head
Recognize the influence to fault location accuracy synchronous with ranging cycle accurate.
The content of the invention
The purpose of the present invention is overcome conventional Time-domain travelling wave ranging requirement fault traveling wave effectively identification and ranging cycle accurate
Synchronous limitation, proposition is a kind of directly to match somebody with somebody overhead line fault distance-finding method based on the calibration of traveling wave mutation distance, on solving
State problem.
The technical scheme is that: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,
When singlephase earth fault occurs for direct distribution lines, under sample rate 1MHz, to fault feeder and adjacent most long perfect feeder line initiating terminal
Faulted phase current sampled, obtain busbar end electricity by adjacent most long feeder line initiating terminal current traveling wave and its wave impedance of perfecting
Traveling wave is pressed, structure is calculated using the adjacent most long current traveling wave for perfecting feeder line origin or beginning current traveling wave and the acquisition of fault feeder measuring end
Make current traveling wave;Distribution and current traveling wave along the voltage traveling wave of measuring end are asked for using obtained voltage traveling wave and current traveling wave
It is distributed along the line, then carries out Directional Decomposition to being distributed along traveling wave, obtain the direction traveling wave being distributed along the line, and be distributed along extraction
Direct wave mutation and backward-travelling wave mutation, finally by the two be multiplied then at traveling wave observe when window in integrated to construct survey
Fault localization is realized away from function.
It is as follows:
(1) when singlephase earth fault occurs for direct distribution lines, under sample rate 1MHz, to the faulted phase current of fault feeder and
The adjacent most long feeder line initiating terminal faulted phase current that perfects is sampled, and is respectively obtained faulted phase current sampled value sequence, is denoted as:i
(k), i ' (k), wherein k expression sampled point, k=1,2 ...;
(2) busbar terminal voltage is obtained according to formula (1) and formula (2) respectively and constructs the discrete series u of current traveling waveM(k) and iM
(k):
uM(k)=i ' (k) × Zc (1)
iM(k)=i (k)-i ' (k) (2)
In formula, uM(k) terminal voltage for being busbar M, iM(k) the construction current traveling wave of fault feeder, ZcFor feeder line wave impedance;
(3) calculating of distribution along:Calculate voltage's distribiuting and the edge along the line of fault feeder respectively using formula (3) and formula (4)
Line current is distributed:
In formula, x is any point along the line to the distance of measuring end;vsFor the wave velocity of circuit;ZcFor feeder line wave impedance;rs
For circuit resistance per unit length;uM(k) terminal voltage for being busbar M;iM(k) it is the construction current traveling wave of fault feeder;uM,x(x,k)
It it is the k moment away from the voltage at measuring end x;iM,x(x, k) is the k moment away from the electric current at measuring end x;
(4) direct wave and backward-travelling wave being distributed along the line are calculated:Fault feeder is calculated according to formula (5) and formula (6) respectively
The forward voltage traveling wave being distributed along the line, the backward voltage traveling wave being distributed along the line, i.e.,:
u+ M,x=(uM,x+ZciM,x)/2 (5)
u- M,x=(uM,x-ZciM,x)/2 (6)
(5) calculating of the direct wave gradient and backward-travelling wave gradient of distribution along:Utilize the forward voltage being distributed along the line
The forward voltage gradient that the difference construction of the two neighboring sampled value of traveling wave is distributed along the line, i.e.,:
c+ M,dif_u(k)=u+ k,x(k)-u+ k,x(k-1) (7)
The backward voltage ladder being distributed along the line using the difference construction for the two neighboring sampled value of backward voltage traveling wave being distributed along the line
Degree, i.e.,:
c- M,dif_u(k)=u- k,x(k)-u- k,x(k-1) (8)
(6) the direct wave mutation being distributed along the line and backward-travelling wave mutation are calculated:It is extracted according to formula (9) along faulty line
The forward voltage traveling wave mutation of distribution, i.e.,:
It extracts the backward voltage traveling wave being distributed along faulty line according to formula (10) to be mutated, i.e.,:
In formula, R is taken as 3;
(7) construction of range function:Using formula (11) and formula (12), direct wave mutation that step (6) is obtained and anti-
Multiplication and window [k when traveling wave is observed are mutated to traveling wave0,k0+l/(2vs)] and [k0+l/(2vs),k0+l/vs] in accumulated
Point, obtain range function fuI(x) and fuII(x) traveling wave mutation along;
In formula, k0Represent the initial traveling wave arrival moment of failure that measuring end M is detected;L is fault feeder line length;
(8) construction of fault location criterion:[k is calculated according to step (7)0,k0+l/(2vs)] and [k0+l/(2vs),
k0+l/vs] two in succession when window in, range function fuI(x) and fuII(x) distribution catastrophe point, respective distances are remembered respectively along
For [xI1,xI2...] and [xII1,xII2... ...], if [xI1,xI2...] and in mutation distance x* I[xII1,xII2,……]
In mutation distance x* IIMeet the line length constraints shown in formula (13), and x* IMutation point-polarity be negative, x* IICatastrophe point
Polarity is negative, x* ICatastrophe point amplitude be less than x* IICatastrophe point amplitude, then failure be located at half line length within, trouble point distance measurements
The distance for surveying end is x* I;If [xI1,xI2...] and in mutation distance x* I[xII1,xII2...] and in mutation distance x* IIIt is full
Line length constraints shown in sufficient formula (13), and x* IMutation point-polarity be negative, x* IIMutation point-polarity be negative, x* IMutation
Point amplitude is more than x* IICatastrophe point amplitude, then failure be located at outside half line length, the distance at trouble point distance measuring end is x* II;
x* I+x* II=l (13).
The beneficial effects of the invention are as follows:
For the present invention for directly fault location is carried out with overhead transmission line, principle is simple, does not need calibration fault traveling wave ripple wave head,
And from the influence of the factors such as failure instantaneity, fault resistance variation, distance measurement result is accurately and reliably.
Description of the drawings
Fig. 1 is embodiment 1, the straight of embodiment 2 matches somebody with somebody overhead system structure chart;
Fig. 2 is [k within half line length under failure0,k0+L1/(2vs)] when window in range function fu(x) mutation distribution knot
Fruit;
Fig. 3 is [k within half line length under failure0+L1/(2vs),k0+L1/vs] when window in range function fu(x) mutation distribution
As a result;
Fig. 4 is [k outside half line length under failure0,k0+L1/(2vs)] when window in range function fu(x) mutation distribution knot
Fruit;
Fig. 5 is [k outside half line length under failure0+L1/(2vs),k0+L1/vs] when window in range function fu(x) mutation distribution
As a result.
Specific embodiment
Below in conjunction with the drawings and specific embodiments, the invention will be further described.
When singlephase earth fault occurs for direct distribution lines, under sample rate 1MHz, to fault feeder and adjacent most long sound feedback
The faulted phase current of line initiating terminal is sampled, and is obtained by adjacent most long feeder line initiating terminal current traveling wave and its wave impedance of perfecting
Busbar terminal voltage traveling wave, and utilize the adjacent most long electric current row for perfecting feeder line origin or beginning current traveling wave and fault feeder measuring end and obtaining
Ripple calculates construction current traveling wave.It asks for being distributed along the voltage traveling wave of measuring end using obtained voltage traveling wave and current traveling wave
It is distributed with along current traveling wave, then carries out Directional Decomposition to being distributed along traveling wave, obtained the direction traveling wave being distributed along the line, and carry
The direct wave that is distributed along the line is taken to be mutated and backward-travelling wave mutation, finally by the two be multiplied then at traveling wave observe when window in accumulate
Divide to construct range function realization fault localization.
Embodiment 1:
Circuit is fed out using 35kV as shown in Figure 1 directly with overhead system, fault feeder L more1=20km, and its whole story two
The connection type for the combination of " I-III " class busbar is held, measuring end is arranged at beginning M.Perfect feeder line L2=5km, L3=15km, and it is strong
Full line end is group iii busbar connection type.Assuming that feeder line L1A phases occur within half line length at distance M end 8km and are grounded event
Barrier, the initial phase angle of failure are 90 °, and transition resistance is set to 0.01 Ω, sample rate 1MHz, using voltage traveling wave and construction electric current row
Ripple, along the line material calculation take 0.1km, window [k when choosing two successive traveling waves observation respectively0,k0+L1/(2vs)] and [k0+L1/
(2vs),k0+L1/vs] expert's wave number evidence, it can be calculated range function fuI(x) and fuII(x) along completely long traveling wave mutation distribution
As shown in Figures 2 and 3.As shown in Figure 2, [k0,k0+l/(2vs)] when window in, fuI(x) catastrophe point A (x)=8km, and polarity is
It is negative;From the figure 3, it may be seen that [k0+l/(2vs),k0+l/vs] when window in, fuII(x) catastrophe point B (x)=12km, and polarity is negative.Cause
For A (x)+B (x)=8+12=20km, meet the line length constraints shown in formula (13), and the amplitude of A (x) is less than the width of B (x)
Value, so failure is located within half line length, the distance of trouble point distance measuring end M is 8km.
Embodiment 2:
Circuit is fed out using 35kV as shown in Figure 1 directly with overhead system, fault feeder L more1=20km, and its whole story two
The connection type for the combination of " I-III " class busbar is held, measuring end is arranged at beginning M.Perfect feeder line L2=5km, L3=15km, and it is strong
Full line end is group iii busbar connection type.Assuming that feeder line L1A phases occur outside half line length at distance M end 14km to be grounded
Failure, the initial phase angle of failure are 90 °, and transition resistance is set to 0.01 Ω, sample rate 1MHz, using voltage traveling wave and construction electric current
Traveling wave, along the line material calculation take 0.1km, window [k when choosing two successive traveling waves observation respectively0,k0+L1/(2vs)] and [k0+
L1/(2vs),k0+L1/vs] expert's wave number evidence, it can be calculated range function fuI(x) and fuII(x) along completely long traveling wave mutation point
Cloth is as shown in Figure 4 and Figure 5.As shown in Figure 4, [k0,k0+l/(2vs)] when window in, fuI(x) catastrophe point B (x)=6km, and polarity
It is negative;As shown in Figure 5, [k0+l/(2vs),k0+l/vs] when window in, fuII(x) catastrophe point A (x)=14km, and polarity is negative.
Because A (x)+B (x)=14+6=20km, meet the line length constraints shown in formula (13), and the amplitude of B (x) is more than A's (x)
Amplitude, so failure is located at outside half line length, the distance of trouble point distance measuring end M is 14km.
Claims (2)
1. a kind of directly match somebody with somebody overhead line fault distance-finding method based on the calibration of traveling wave mutation distance, it is characterised in that:Direct distribution lines are sent out
It is mutually electric with the adjacent most long failure for perfecting feeder line initiating terminal to fault feeder under sample rate 1MHz during raw singlephase earth fault
Stream is sampled, and busbar terminal voltage traveling wave, profit are obtained by adjacent most long feeder line initiating terminal current traveling wave and its wave impedance of perfecting
Construction electric current row is calculated with the adjacent most long current traveling wave for perfecting feeder line origin or beginning current traveling wave and the acquisition of fault feeder measuring end
Ripple;It asks for being distributed along the voltage traveling wave of measuring end using obtained voltage traveling wave and current traveling wave and current traveling wave divides along the line
Then cloth carries out Directional Decomposition to being distributed along traveling wave, obtain the direction traveling wave being distributed along the line, and extract the forward direction being distributed along the line
Traveling wave be mutated and backward-travelling wave mutation, finally by the two be multiplied then at traveling wave observe when window in integrated to construct range function
Realize fault localization.
2. according to claim 1 directly match somebody with somebody overhead line fault distance-finding method, feature based on the calibration of traveling wave mutation distance
It is to be as follows:
(1) when singlephase earth fault occurs for direct distribution lines, under sample rate 1MHz, faulted phase current to fault feeder and adjacent
The most long feeder line initiating terminal faulted phase current that perfects is sampled, and is respectively obtained faulted phase current sampled value sequence, is denoted as:i(k)、
I ' (k), wherein k expression sampled point, k=1,2 ...;
(2) busbar terminal voltage is obtained according to formula (1) and formula (2) respectively and constructs the discrete series u of current traveling waveM(k) and iM(k):
uM(k)=i ' (k) × Zc (1)
iM(k)=i (k)-i ' (k) (2)
In formula, uM(k) terminal voltage for being busbar M, iM(k) the construction current traveling wave of fault feeder, ZcFor feeder line wave impedance;
(3) calculating of distribution along:Calculated respectively using formula (3) and formula (4) fault feeder along voltage's distribiuting and along the line electricity
Flow distribution:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mi>s</mi>
</msub>
<mi>x</mi>
<mo>/</mo>
<mn>4</mn>
</mrow>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msub>
<mi>u</mi>
<mi>M</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mfrac>
<mi>x</mi>
<msub>
<mi>v</mi>
<mi>s</mi>
</msub>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>i</mi>
<mi>M</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mfrac>
<mi>x</mi>
<msub>
<mi>v</mi>
<mi>s</mi>
</msub>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>r</mi>
<mi>s</mi>
</msub>
<mi>x</mi>
</mrow>
<mn>4</mn>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mi>s</mi>
</msub>
<mi>x</mi>
<mo>/</mo>
<mn>4</mn>
</mrow>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msub>
<mi>u</mi>
<mi>M</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mfrac>
<mi>x</mi>
<msub>
<mi>v</mi>
<mi>s</mi>
</msub>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>i</mi>
<mi>M</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mfrac>
<mi>x</mi>
<msub>
<mi>v</mi>
<mi>s</mi>
</msub>
</mfrac>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mi>s</mi>
</msub>
<mi>x</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>r</mi>
<mi>s</mi>
</msub>
<mi>x</mi>
<mo>/</mo>
<mn>4</mn>
</mrow>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<msub>
<mi>u</mi>
<mi>M</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>r</mi>
<mi>s</mi>
</msub>
<mi>x</mi>
</mrow>
<mn>4</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mi>s</mi>
</msub>
<mi>x</mi>
<mo>/</mo>
<mn>4</mn>
</mrow>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mi>s</mi>
</msub>
<mi>x</mi>
<mo>/</mo>
<mn>4</mn>
</mrow>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>i</mi>
<mi>M</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>i</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mi>s</mi>
</msub>
<mi>x</mi>
<mo>/</mo>
<mn>4</mn>
</mrow>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msub>
<mi>u</mi>
<mi>M</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>x</mi>
<mo>/</mo>
<msub>
<mi>v</mi>
<mi>s</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>i</mi>
<mi>M</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mi>x</mi>
<mo>/</mo>
<msub>
<mi>v</mi>
<mi>s</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mi>s</mi>
</msub>
<mi>x</mi>
<mo>/</mo>
<mn>4</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mi>s</mi>
</msub>
<mi>x</mi>
<mo>/</mo>
<mn>4</mn>
</mrow>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msub>
<mi>u</mi>
<mi>M</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>x</mi>
<mo>/</mo>
<msub>
<mi>v</mi>
<mi>s</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>i</mi>
<mi>M</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>x</mi>
<mo>/</mo>
<msub>
<mi>v</mi>
<mi>s</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mi>s</mi>
</msub>
<mi>x</mi>
<mo>/</mo>
<mn>4</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<mfrac>
<mrow>
<msub>
<mi>r</mi>
<mi>s</mi>
</msub>
<mi>x</mi>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<msub>
<mi>u</mi>
<mi>M</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>i</mi>
<mi>M</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>r</mi>
<mi>s</mi>
</msub>
<mi>x</mi>
<mo>/</mo>
<mn>4</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, x is any point along the line to the distance of measuring end;vsFor the wave velocity of circuit;ZcFor feeder line wave impedance;rsFor line
Road resistance per unit length;uM(k) terminal voltage for being busbar M;iM(k) it is the construction current traveling wave of fault feeder;uM,x(x, k) is k
Moment is away from the voltage at measuring end x;iM,x(x, k) is the k moment away from the electric current at measuring end x;
(4) direct wave and backward-travelling wave being distributed along the line are calculated:It is calculated respectively along fault feeder according to formula (5) and formula (6)
The forward voltage traveling wave of distribution, the backward voltage traveling wave being distributed along the line, i.e.,:
u+ M,x=(uM,x+ZciM,x)/2 (5)
u- M,x=(uM,x-ZciM,x)/2 (6)
(5) calculating of the direct wave gradient and backward-travelling wave gradient of distribution along:Utilize the forward voltage traveling wave being distributed along the line
The forward voltage gradient that the difference construction of two neighboring sampled value is distributed along the line, i.e.,:
c+ M,dif_u(k)=u+ k,x(k)-u+ k,x(k-1) (7)
The backward voltage gradient being distributed along the line using the difference construction for the two neighboring sampled value of backward voltage traveling wave being distributed along the line,
I.e.:
c- M,dif_u(k)=u- k,x(k)-u- k,x(k-1) (8)
(6) the direct wave mutation being distributed along the line and backward-travelling wave mutation are calculated:It is extracted according to formula (9) and is distributed along faulty line
Forward voltage traveling wave mutation, i.e.,:
<mrow>
<msub>
<msup>
<mi>S</mi>
<mo>+</mo>
</msup>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mn>2</mn>
<mi>u</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mi>k</mi>
<mo>-</mo>
<mi>R</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<msup>
<mi>c</mi>
<mo>+</mo>
</msup>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mi>d</mi>
<mi>i</mi>
<mi>f</mi>
<mo>_</mo>
<mi>u</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>3</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
It extracts the backward voltage traveling wave being distributed along faulty line according to formula (10) to be mutated, i.e.,:
<mrow>
<msub>
<msup>
<mi>S</mi>
<mo>-</mo>
</msup>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mn>2</mn>
<mi>u</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mi>k</mi>
<mo>-</mo>
<mi>R</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<msup>
<mi>c</mi>
<mo>-</mo>
</msup>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mi>d</mi>
<mi>i</mi>
<mi>f</mi>
<mo>_</mo>
<mi>u</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>3</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, R is taken as 3;
(7) construction of range function:Using formula (11) and formula (12), the direct wave mutation that step (6) is obtained and reversed row
Ripple mutation is multiplied and window [k when traveling wave is observed0,k0+l/(2vs)] and [k0+l/(2vs),k0+l/vs] in integrated,
Obtain range function fuI(x) and fuII(x) traveling wave mutation along;
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>u</mi>
<mi>I</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>k</mi>
<mn>0</mn>
</msub>
<mrow>
<msub>
<mi>k</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mi>l</mi>
<mo>/</mo>
<mn>2</mn>
<mi>v</mi>
</mrow>
</msubsup>
<msubsup>
<mi>S</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mn>2</mn>
<mi>u</mi>
</mrow>
<mo>+</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<msubsup>
<mi>S</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mn>2</mn>
<mi>u</mi>
</mrow>
<mo>-</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>k</mi>
<mo>,</mo>
<mi>x</mi>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mn>0</mn>
<mo>,</mo>
<mi>l</mi>
<mo>/</mo>
<mn>2</mn>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>u</mi>
<mi>I</mi>
<mi>I</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mrow>
<msub>
<mi>k</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mi>l</mi>
<mo>/</mo>
<mn>2</mn>
<mi>v</mi>
</mrow>
<mrow>
<msub>
<mi>k</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mi>l</mi>
<mo>/</mo>
<mi>v</mi>
</mrow>
</msubsup>
<msubsup>
<mi>S</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mn>2</mn>
<mi>u</mi>
</mrow>
<mo>+</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<msubsup>
<mi>S</mi>
<mrow>
<mi>M</mi>
<mo>,</mo>
<mn>2</mn>
<mi>u</mi>
</mrow>
<mo>-</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>k</mi>
<mo>,</mo>
<mi>x</mi>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mi>l</mi>
<mo>/</mo>
<mn>2</mn>
<mo>,</mo>
<mi>l</mi>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, k0Represent the initial traveling wave arrival moment of failure that measuring end M is detected;L is fault feeder line length;
(8) construction of fault location criterion:[k is calculated according to step (7)0,k0+l/(2vs)] and [k0+l/(2vs),k0+l/
vs] two in succession when window in, range function fuI(x) and fuII(x) distribution catastrophe point, respective distances are denoted as respectively along
[xI1,xI2...] and [xII1,xII2... ...], if [xI1,xI2...] and in mutation distance x* I[xII1,xII2...] in
Mutation distance x* IIMeet the line length constraints shown in formula (13), and x* IMutation point-polarity be negative, x* IIMutation point pole
Property is negative, x* ICatastrophe point amplitude be less than x* IICatastrophe point amplitude, then failure be located at half line length within, trouble point distance measuring
The distance at end is x* I;If [xI1,xI2...] and in mutation distance x* I[xII1,xII2...] and in mutation distance x* IIMeet
Line length constraints shown in formula (13), and x* IMutation point-polarity be negative, x* IIMutation point-polarity be negative, x* ICatastrophe point
Amplitude is more than x* IICatastrophe point amplitude, then failure be located at outside half line length, the distance at trouble point distance measuring end is x* II;
x* I+x* II=l (13).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510606519.9A CN105116295B (en) | 2015-09-22 | 2015-09-22 | 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 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510606519.9A CN105116295B (en) | 2015-09-22 | 2015-09-22 | 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 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105116295A CN105116295A (en) | 2015-12-02 |
CN105116295B true CN105116295B (en) | 2018-05-25 |
Family
ID=54664330
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510606519.9A Active CN105116295B (en) | 2015-09-22 | 2015-09-22 | 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 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105116295B (en) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105548802B (en) * | 2015-12-04 | 2019-02-19 | 昆明理工大学 | A kind of asynchronous fault distance-finding method in three end of T link based on distribution character along fault traveling wave |
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 |
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 |
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 |
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 |
CN105738771B (en) * | 2016-04-01 | 2019-07-05 | 昆明理工大学 | A kind of one-end fault ranging method based on the route containing TCSC of distribution character along fault traveling wave |
CN105891672B (en) * | 2016-04-01 | 2020-02-07 | 昆明理工大学 | UPFC-containing line traveling wave distance measurement method based on fault traveling wave distribution characteristics along line |
CN105891671B (en) * | 2016-04-01 | 2020-04-03 | 昆明理工大学 | STATCOM-containing line single-end traveling wave distance measurement method based on fault traveling wave distribution characteristics along line |
CN105842584B (en) * | 2016-04-01 | 2019-07-16 | 昆明理工大学 | A kind of T link multiterminal travelling wave ranging method based on distribution character along fault traveling wave |
CN105866624B (en) * | 2016-04-01 | 2019-05-14 | 昆明理工大学 | A kind of both-end distance measuring method using distribution character along fault traveling wave |
CN109739215B (en) * | 2019-01-21 | 2021-05-04 | 济南大学 | Method for realizing fault initial angle control and automatic traversal in PSCAD |
EP3723224B1 (en) | 2019-04-08 | 2023-08-23 | Hitachi Energy Switzerland AG | Time domain distance protection based on polarity comparison for power transmission lines |
CN113406436B (en) * | 2021-06-17 | 2022-08-26 | 山东大学 | Traveling wave fault location method and system for alternating-current and direct-current transmission line based on 5G communication |
CN113466623B (en) * | 2021-06-29 | 2022-04-19 | 华中科技大学 | Traction network fault location method, system and medium based on transient sudden change current |
CN115015686A (en) * | 2022-01-26 | 2022-09-06 | 昆明理工大学 | LCC-VSC hybrid high-voltage direct-current transmission line fault distance measurement method and system |
CN115047284B (en) * | 2022-01-26 | 2024-07-26 | 昆明理工大学 | High-voltage direct-current transmission line fault location method and system |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6798211B1 (en) * | 1997-10-30 | 2004-09-28 | Remote Monitoring Systems, Inc. | Power line fault detector and analyzer |
CN101252274A (en) * | 2008-04-03 | 2008-08-27 | 昆明理工大学 | Same tower double back transmission line fault distance measuring time domain method using single end current flow |
CN103018627A (en) * | 2012-10-25 | 2013-04-03 | 福建省电力有限公司 | Adaptive fault type fault line detection method for non-effectively earthed system |
CN103809082A (en) * | 2014-02-17 | 2014-05-21 | 四川大学 | Distance measurement method for power distribution network single-phase earth fault on the basis of aerial mode traveling wave mutation |
CN104133158A (en) * | 2014-08-04 | 2014-11-05 | 昆明理工大学 | Distribution network fault line selection method based on zero-mode current multi-order difference transformation |
CN104155575A (en) * | 2014-08-04 | 2014-11-19 | 昆明理工大学 | Cable mixed circuit failure zone identification method adopting line model current for PCA cluster analysis |
CN104459458A (en) * | 2015-01-09 | 2015-03-25 | 昆明理工大学 | Method for recognizing circuit switching on and failures through specific short window integral value of directional traveling waves |
-
2015
- 2015-09-22 CN CN201510606519.9A patent/CN105116295B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6798211B1 (en) * | 1997-10-30 | 2004-09-28 | Remote Monitoring Systems, Inc. | Power line fault detector and analyzer |
CN101252274A (en) * | 2008-04-03 | 2008-08-27 | 昆明理工大学 | Same tower double back transmission line fault distance measuring time domain method using single end current flow |
CN103018627A (en) * | 2012-10-25 | 2013-04-03 | 福建省电力有限公司 | Adaptive fault type fault line detection method for non-effectively earthed system |
CN103809082A (en) * | 2014-02-17 | 2014-05-21 | 四川大学 | Distance measurement method for power distribution network single-phase earth fault on the basis of aerial mode traveling wave mutation |
CN104133158A (en) * | 2014-08-04 | 2014-11-05 | 昆明理工大学 | Distribution network fault line selection method based on zero-mode current multi-order difference transformation |
CN104155575A (en) * | 2014-08-04 | 2014-11-19 | 昆明理工大学 | Cable mixed circuit failure zone identification method adopting line model current for PCA cluster analysis |
CN104459458A (en) * | 2015-01-09 | 2015-03-25 | 昆明理工大学 | Method for recognizing circuit switching on and failures through specific short window integral value of directional traveling waves |
Non-Patent Citations (2)
Title |
---|
一种基于行波测距的输电线路接地故障距离保护方案;司大军等;《电工技术学报》;20030831;第18卷(第04期);第65-69页 * |
一种基于高频量衰减特性的特高压直流输电线路故障测距方法;陈仕龙等;《电力系统保护与控制》;20140516;第42卷(第10期);第77-83页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105116295A (en) | 2015-12-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
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 | |
CN105182184B (en) | A kind of distribution cable mixing direct distribution lines fault distance-finding method decomposed based on direction traveling wave along the line | |
CN101672883B (en) | Overhead and cable mixed power transmission line fault locating method and device | |
CN101718833B (en) | Method of single end distance measurement of power transmission line malfunction based on traveling wave inherent frequency extraction | |
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 | |
CN103869221B (en) | Based on the double-circuit line singlephase earth fault distance measurement method of SV network samples | |
CN103869220B (en) | Based on directly adopting straight jumping communication mode double-circuit line method for locating single-phase ground fault | |
CN106959401A (en) | A kind of distributed travelling wave ranging method | |
CN104155582A (en) | Distribution network line fault section positioning method based on full-waveform information | |
CN104569744A (en) | Comprehensive single-end fault positioning method applicable to power distribution network lines | |
CN101706533B (en) | Method for measuring zero sequence impedance parameters of transmission lines by utilizing current of line ends | |
CN102129014B (en) | Method for realizing line phase-to-phase fault single-end distance measurement by utilizing distribution parameter model | |
CN102967779B (en) | Identifying method of distribution parameters of transmission line | |
CN103235237A (en) | High-voltage direct-current grounding electrode line high-resistance fault positioning method | |
CN103399209A (en) | Method for measuring power frequency parameters of ultra-high voltage bipolar direct current (DC) transmission line | |
CN104898021A (en) | Fault line selecting method based on k-means cluster analysis for power distribution network | |
CN105137289B (en) | Suitable for the traveling wave fault positioning method of high voltage overhead lines-cable hybrid line | |
CN104090210B (en) | The non-same famous prime minister's cross-line earth fault method of single end distance measurement of double-circuit line | |
CN104865495A (en) | Travelling wave positioning method and system for electric power line fault | |
CN103760461A (en) | Bus protection method based on sudden-change direction of fault current containing power frequency bands | |
CN104655981A (en) | Comparison type directional protection method for travelling wave polarity by utilizing crossing sequence differential transformation | |
CN107329045A (en) | Distribution Network Failure least square location algorithm | |
CN109444665A (en) | The unrelated Fault Locating Method of power distribution network multiple-limb overhead transmission line parameter and system | |
CN105842584B (en) | A kind of T link multiterminal travelling wave ranging method based on distribution character along fault traveling wave | |
CN107247215B (en) | Distribution Network Failure population location algorithm based on Multipoint synchronous measurement data |
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 |