CN107219440A - The localization method of single-ended radiation type distribution network singlephase earth fault - Google Patents
The localization method of single-ended radiation type distribution network singlephase earth fault Download PDFInfo
- Publication number
- CN107219440A CN107219440A CN201710232520.9A CN201710232520A CN107219440A CN 107219440 A CN107219440 A CN 107219440A CN 201710232520 A CN201710232520 A CN 201710232520A CN 107219440 A CN107219440 A CN 107219440A
- Authority
- CN
- China
- Prior art keywords
- zero
- mrow
- phase difference
- feeder line
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/086—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Locating Faults (AREA)
Abstract
The invention discloses the localization method of single-ended radiation type distribution network singlephase earth fault, the distortion degree of zero line wave height frequency component is characterized using specific frequency Amplitude Ration in the initial zero mode voltage traveling wave Mintrop wave head signal of measurement point and phase difference, its core concept is, select location in distribution analogue system applies singlephase earth fault, extract the mode voltage traveling wave Mintrop wave head height frequency component amplitude of feeder line first and last end zero and phase sets up standard Amplitude Ration vector sum phase difference vector, compared by the Amplitude Ration and phase difference that obtain physical fault calculating with standard Amplitude Ration and phase difference and positioned come route selection and part of path, and branch's positioning is realized according to the energy of zero mould current traveling wave Mintrop wave head height frequency component.Emulation shows that this method is not influenceed by the factor such as electric arc and synchro measure error, with higher economy and stronger practical value.
Description
Technical field
The invention belongs to electric network fault field of locating technology, more particularly to single-ended radiation type distribution network singlephase earth fault
Localization method.
Background technology
With the propulsion of intelligent power distribution network construction, FLT searches failure for rapid, during fast recovery of power supply
Between, reduce each side's economic loss important in inhibiting.Power distribution network it is topological more complicated, its line connection is also a variety of many
Sample, various disequilibrium generally existings, these factors bring many challenges to fault location.
According to the scope of fault location, single-ended radial pattern Distribution Network Failure localization method can be broadly divided into two major classes：One class is
The fault zone substantially estimation technique, including route selection, branch's (section) and part of path localization method；Another kind of is that abort situation is accurate
Calculating method, can be divided into impedance method and traveling wave method again.Fault Locating Method based on traveling wave typically utilizes the transient state after failure generation
Traveling wave information, recognizes traveling wave Mintrop wave head using certain mathematical method or wave head reaches the moment thereafter, and opened up according to distribution
Information is flutterred, builds specific fault location criterion to solve.
The subject matter that traditional traveling wave method is applied to distribution network failure positioning presence has：(1) the complicated topology of distribution line
Structure causes traveling wave catadioptric rule complex；(2) it is accurate in distribution to determine traveling wave Mintrop wave head arrival time and traveling wave ripple
Speed is more difficult；(3) traveling wave method requires higher to measurement synchronism；(4) amplitude attenuation and delayed phase of traveling wave high fdrequency component
Rule is complicated, the more difficult accurate calculating of its value.
The content of the invention
In order to solve the technical problem that above-mentioned background technology is proposed, the present invention is intended to provide single-ended radiation type distribution network is single-phase
The localization method of earth fault, overcomes traditional traveling wave method to realize that complicated, synchronism requires high defect.
In order to realize above-mentioned technical purpose, the technical scheme is that：
The localization method of single-ended radiation type distribution network singlephase earth fault, comprises the following steps：
(1) before failure generation, the simulation model of single-ended radiation type distribution network system is pre-established, in system substation bus
Place's simulation singlephase earth fault, gathers each feeder line head end and the mode voltage traveling wave Mintrop wave head signal of end zero, and extract these letters
Frequency is the amplitude and phase of ω high-frequency signal in number；
(2) for any feeder line, the zero mode voltage traveling wave Mintrop wave head signal intermediate frequency rate for calculating the feeder line head end and end is
The Amplitude Ration and phase difference absolute value of ω high-frequency signal, obtain n dimension Amplitude Ration vector sum n dimension phase difference vectors, n is distribution
Feeder line number in net system；
(3) for mixing feeder line, it is made up of if certain mixes feeder line N sections of different wave impedance circuits, then the ripple of the mixing feeder line
The number of impedance discontinuity point is N+1, simulates singlephase earth fault at this N+1 point respectively, gathers the N+1 mixing feeder line
Head end and the mode voltage traveling wave Mintrop wave head signal of end zero, and extract width of this N+1 times signal intermediate frequency rate for ω high-frequency signal
Value and phase, calculate high-frequency signal of the zero mode voltage traveling wave Mintrop wave head signal intermediate frequency rate of N+1 signal head end and end for ω
Amplitude Ration and phase difference absolute value, obtain N+1 dimension Amplitude Ration vector sum N+1 dimension phase difference vectors；
(4) after actual power distribution network breaks down, each feeder line head end and the mould of end zero electricity are extracted according to the method for step (1)
It is the amplitude and phase of ω high-frequency signal to press traveling wave Mintrop wave head signal intermediate frequency rate, and obtains reality according to the method for step (2)
N dimension Amplitude Ration vector sum n dimension phase difference vectors after failure；
(5) ratio of the Amplitude Ration for the same feeder line that the Amplitude Ration that calculation procedure (2) is obtained is obtained with step (4), and
The ratio of the phase difference absolute value for the same feeder line that the phase difference absolute value that step (2) is obtained is obtained with step (4), obtains n dimensions
Amplitude Ration ratio vector sum n dimension phase difference ratio vectors；
(6) greatest member in n dimension Amplitude Ration ratios vector, n dimension phase difference ratio vectors is found out respectively, if two vectors
Greatest member it is corresponding be same feeder line, then the feeder line is faulty line, is if two vectorial greatest members are corresponding
Different feeder lines, the then corresponding feeder line of greatest member that n is tieed up to Amplitude Ration ratio vector is used as faulty line；
(7) if faulty line is mixing feeder line, the n that step (4) is obtained ties up Amplitude Ration vector sum n and tieed up in phase difference vector
The corresponding element of the faulty line obtains all elements in N+1 dimension Amplitude Ration vector sum N+1 dimension phase difference vectors with step (3) and entered
Row compares, and is tieed up from N+1 and the element for meeting following formula is found out during Amplitude Ration vector sum N+1 ties up phase difference vector：
SARm≥a≥SARm+1
SPDw≤p≤SPDw+1
In above formula, a is the corresponding element of faulty line in the n dimension Amplitude Ration vectors that step (4) is obtained, and p obtains for step (4)
To n dimension phase difference vector in the corresponding element of faulty line, SARmM-th of element in Amplitude Ration vector, SPD are tieed up for N+1w
W-th of element in phase difference vector, m, w ∈ [1, N], if SAR are tieed up for N+1mWith SPDwCorresponding is same part of path, then should
Part of path is the fault wire section of faulty line, if SARmWith SPDwCorresponding is different part of paths, then SARmCorresponding part of path
For the fault wire section of faulty line；
(8) if faulty line is branch feeder, extract zero mode voltage traveling wave Mintrop wave head at each branch point of the branch feeder and believe
Frequency is the complex vector of ω high-frequency signal in number；
(9) energy of each mode voltage traveling wave Mintrop wave head signal intermediate frequency rate of branch point zero for the complex vector of ω high-frequency signal is calculated
Amount, then the maximum branch of energy is the fault branch of the faulty line.
Further, when gathering zero mode voltage traveling wave Mintrop wave head signal, it is thus necessary to determine that calculate data window, its determination method
It is as follows：
(a) catastrophe point of zero mould travelling wave signal is looked for, if zero line ripple signal vector of measurement point collection is U0=[u1,
u2,...,uk-1,uk,uk+1,...,uL], wherein L is signal total length, should be if the line ripple catastrophe point of failure zero is k simultaneously
Meet following two formula：
In above formula, uk-2,uk-1,ukIt is not zero, ε is the threshold value of setting；
(b) 100 sampled points are chosen before catastrophe point k and X sampled point is as calculating data window thereafter, wherein X is by following formula
It is determined that：
In above formula, LminFor the length of most short feeder line in all feeder lines in power distribution network, v is the mould of high frequency zero electricity that frequency is ω
Press the velocity of wave of traveling wave.
Further, threshold epsilon=2.
Further, S-transformation is used to extract width of the zero mode voltage traveling wave Mintrop wave head signal intermediate frequency rate for ω high-frequency signal
Value, phase and complex vector.
Further, frequencies omega=45kHz of zero mode voltage traveling wave Mintrop wave head signal high frequency signal.
Further, in step (9), each mode voltage traveling wave Mintrop wave head signal intermediate frequency rate of branch point zero is believed for ω high frequency
Number complex vector energy calculation formula it is as follows：
In above formula, S is the complex vector for the high-frequency signal that certain mode voltage traveling wave Mintrop wave head signal intermediate frequency rate of branch point zero is ω,
E is S energy, and Amp (S) represents to seek S modulus value.
The beneficial effect brought using above-mentioned technical proposal：
The core of the present invention is to pre-establish bus by analog simulation and mix when feeder line respectively connects point failure respectively to present
The Amplitude Ration and phase difference array of the mode voltage traveling wave Mintrop wave head height frequency component of line first and end zero, recycle S-transformation quick, accurate
The phase of frequency component needed for extracting.The present invention implements simple, is easy to implement, due to need to only extract the amplitude and phase of a certain frequency
Position, therefore without zero mode voltage traveling wave speed and transient state travelling wave Mintrop wave head due in, and without precise synchronization measurement, with compared with
High practical engineering value.
Brief description of the drawings
Fig. 1 is typical single-ended radial pattern joint line power distribution network schematic diagram；
Fig. 2 is the flow chart of the inventive method；
Fig. 3 builds single-ended radial pattern joint line distribution network model schematic diagram by simulating, verifying；
Fig. 4 is arc current simulation waveform；
Fig. 5 is the end measurement point primary signal of feeder line 1 and considers 10 μ s, 30 μ s, the signal graph of 50 μ s synchro measure errors.
Embodiment
Below with reference to accompanying drawing, technical scheme is described in detail.
Typical single-ended radial pattern joint line power distribution network is illustrated in fig. 1 shown below, and the present invention proposes a kind of for single-ended spoke
The localization method of emitting one-phase earthing failure in electric distribution network, step is as shown in Figure 2.
Step 1：Before failure generation, the simulation model of single-ended radiation type distribution network system is pre-established, in system substation
Singlephase earth fault is simulated at bus, each feeder line head end and the mode voltage traveling wave Mintrop wave head signal of end zero are gathered, using S-transformation
Extract the amplitude and phase for the high-frequency signal that these signal intermediate frequency rates are ω.Preset ω=45kHz.
Step 2：For any feeder line, zero mode voltage traveling wave Mintrop wave head signal intermediate frequency rate of the feeder line head end and end is calculated
Amplitude Ration and phase difference absolute value for ω high-frequency signal, obtain n dimension Amplitude Ration vector AR and n dimension phase difference vectors PD, n
For the feeder line number in distribution network system.
Step 3：For mixing feeder line, if certain mixing feeder line is made up of N section difference wave impedance circuits, then the mixing feeder line
The number of wave impedance discontinuity point is N+1, simulates singlephase earth fault, collection N+1 mixing feedback at this N+1 point respectively
Line head end and the mode voltage traveling wave Mintrop wave head signal of end zero, and it is ω high-frequency signal to extract this N+1 times signal intermediate frequency rate
Amplitude and phase, calculate zero mode voltage traveling wave Mintrop wave head signal intermediate frequency rate of N+1 signal head end and end and believe for ω high frequency
Number Amplitude Ration and phase difference absolute value, obtain N+1 dimensions Amplitude Ration vector SAR and N+1 and tie up phase difference vector SPD.
Step 4：After actual power distribution network breaks down, each feeder line head end and the mould of end zero are extracted according to the method for step 1
Voltage traveling wave Mintrop wave head signal intermediate frequency rate is the amplitude and phase of ω high-frequency signal, and obtains reality according to the method for step 2
N dimension Amplitude Ration vector RAR and n dimension phase difference vectors RPD after failure.
Step 5：The ratio of the Amplitude Ration for the same feeder line that the Amplitude Ration that calculation procedure 2 is obtained is obtained with step 4, Yi Jibu
The ratio of the phase difference absolute value for the same feeder line that rapid 2 obtained phase difference absolute values are obtained with step 4, obtains n dimension Amplitude Rations
Ratio vector CAR and n dimension phase difference ratio vector CPD.
Step 6：The greatest member in CAR, CPD is found out respectively, if it is same feedback that two vectorial greatest members are corresponding
Line, then the feeder line is faulty line, if it is different feeder lines that two vectorial greatest members are corresponding, by the maximum in CAR
The corresponding feeder line of element is used as faulty line.
Step 7：If faulty line is mixing feeder line, the corresponding member of the faulty line in RAR and RPD that step 4 is obtained
All elements are compared in SAR and SPD that element is obtained with step 3, and the element for meeting following formula is found out from SAR and SPD：
SARm≥a≥SARm+1
SPDw≤p≤SPDw+1
In above formula, a is the corresponding element of faulty line in RAR, and p is the corresponding element of faulty line, SAR in RPDmFor
M-th of element in SAR, SPDwFor w-th of element in SPD, m, w ∈ [1, N], if SARmWith SPDwCorresponding is same line
Section, then the part of path is the fault wire section of faulty line, if SARmWith SPDwCorresponding is different part of paths, then SARmIt is right
The part of path answered is the fault wire section of faulty line.
Step 8：If faulty line is branch feeder, zero mode voltage at each branch point of the branch feeder is extracted by S-transformation
Traveling wave Mintrop wave head signal intermediate frequency rate is the complex vector of ω high-frequency signal.
Step 9：Calculate the complex vector for the high-frequency signal that each mode voltage traveling wave Mintrop wave head signal intermediate frequency rate of branch point zero is ω
Energy, then the maximum branch of energy is the fault branch of the faulty line：
In above formula, S is the complex vector for the high-frequency signal that certain mode voltage traveling wave Mintrop wave head signal intermediate frequency rate of branch point zero is ω,
E is S energy, and Amp (S) represents to seek S modulus value.
Carry out fault location, it is necessary first to accurate to extract zero line ripple Mintrop wave head signal.Distribution medium wave impedance discontinuity
Catadioptric occurs for the initial zero line ripple that the presence of point can to propagate to herein, if initial zero line ripple Mintrop wave head signal
If the length selection of calculating data window is improper, then the reflecting component of certain ingredients can be included in its signal.Present invention determine that meter
The method for calculating data window is as follows：
(1) catastrophe point of zero mould travelling wave signal is looked for, if zero line ripple signal vector of measurement point collection is U0=[u1,
u2,...,uk-1,uk,uk+1,...,uL], wherein L is signal total length, should be if the line ripple catastrophe point of failure zero is k simultaneously
Meet following two formula：
In above formula, uk-2,uk-1,ukIt is not zero, ε is the threshold value of setting, through emulation experiment, when ε values are 2, Ke Yibao
Card accurately obtains the line ripple sign mutation point of failure zero in all cases.
(2) 100 sampled points are chosen before catastrophe point k and X sampled point is as calculating data window thereafter, wherein X is by following formula
It is determined that：
In above formula, LminFor the length of most short feeder line in all feeder lines in power distribution network, v is the mould of high frequency zero electricity that frequency is ω
The velocity of wave of traveling wave is pressed, the value need not be calculated accurately, can be taken as the light velocity.According to this principle, X values are adapted to each for 10
The situation of kind.
Simulating, verifying
In order to verify the correctness of the present invention, the single-ended radial pattern joint line distribution of a 10Kv have been built on PSCAD
Net, as shown in Figure 3.For convenience's sake, all overhead lines use same configuration, and all cables are also matched somebody with somebody using identical
Put, overhead line and cable in practice is characterized using accurate frequency dependent character phase model.The measurement point arrangement of each feeder line
As it was previously stated, sample frequency is 1MHz.It is resulting by after simulated failure at bus and the specified point of Hybrid connections feeder line 1
Route selection and part of path localization criteria Amplitude Ration and phase difference array, it is as shown in table 1 below.Phase angle difference list in all forms herein
Position is radian.In order to verify influence of the various factors to context of methods, respectively in different neutral grounding modes, failure electricity
Emulated in the presence of resistance, failure initial phase angle, electric arc and synchro measure error.
The failure line selection of table 1 and part of path localization criteria Amplitude Ration and phase difference array
1. the influence of different grounding modes, fault resstance, fault distance
It is different with the overhead line of feeder line 3 in the cable of feeder line 1 respectively in order to verify influence of these factors to the present invention
At position, under two kinds of different earthing modes of isolated neutral and neutral by arc extinction coil grounding, fault resstance be 10 Ω,
When 200 Ω, 500 Ω, the emulation of A phases earth fault is carried out under 42 kinds of different situations altogether, shown in such as table 2 (a) and 2 (b).
Result of calculation in 1 cable diverse location of feeder line occurs for table 2 (a) physical fault
Result of calculation in 3 overhead line diverse location of feeder line occurs for table 2 (b) physical fault
Not Wei failure occur result of calculation when in 3 overhead line diverse location of the cable of feeder line 1 and feeder line, wherein failure away from
From the distance for referring to feeder line head end substation bus with a distance from trouble point.Result be can be seen that from table, and this method is not influenceed by factors above.
2. the influence in the presence of electric arc
In order to verify influence of the electric arc produced during single-phase earthing to context of methods, using best suiting actual field situation
Cybernetics Arc Modelling.Apply A phases in isolated neutral, at the diverse location of different feeder line and be grounded arc fault, survey
The arc current obtained is as shown in figure 4, it is seen that it is about 50ms that the moment, which occurs, for arc fault, instant of failure produces electric current
Spike signal, after after arc stability, positive and negative asymmetry is presented on the whole for arc current waveform and each cycle has two
Secondary " zero stops " phenomenon.Simulation result as shown in table 3, the fault distance in table refer at trouble point to feeder line head end substation bus away from
From.As can be seen from the table, even if occurring arc fault, this method remains able to reliably working.
Result of calculation during arc fault occurs at the different faults position of table 3
3. the influence of synchro measure error
In order to verify influence of the synchro measure error to this paper algorithms, zero mode voltage that each feeder line head end measurement point is measured
Travelling wave signal keep it is constant, and the zero mode voltage traveling wave Mintrop wave head signal that end measurement point is measured on the original basis respectively to
Postpone 10 afterwards, 30,50 points, and in the case of property point is earth-free in systems, respectively the cable of feeder line 1, the overhead line of feeder line 2,
The midpoint of the overhead line of feeder line 3 simulates single-phase electric arc type earth fault.When Fig. 5 is that arcing ground occurs for the cable midpoint of feeder line 1,
Signal graph obtained by zero mode voltage travelling wave signal backward delay 10 that the end measurement point of feeder line 1 is obtained, 30,50 points.
As shown in table 4, wherein SME is synchro measure error to final calculation result.The determination of signal data window is only relevant with catastrophe point, because
This actually synchro measure error on context of methods without influence, it is also seen that this point from table.
Table 4 considers result of calculation during synchro measure error
The technological thought of embodiment only to illustrate the invention, it is impossible to which protection scope of the present invention is limited with this, it is every according to
Technological thought proposed by the present invention, any change done on the basis of technical scheme, each falls within the scope of the present invention.
Claims (6)
1. the localization method of single-ended radiation type distribution network singlephase earth fault, it is characterised in that comprise the following steps：
(1) before failure generation, the simulation model of single-ended radiation type distribution network system is pre-established, the mould at the bus of system substation
Quasi-single- phase earth fault, gathers each feeder line head end and the mode voltage traveling wave Mintrop wave head signal of end zero, and extract in these signals
Frequency is the amplitude and phase of ω high-frequency signal；
(2) for any feeder line, it is ω's to calculate zero mode voltage traveling wave Mintrop wave head signal intermediate frequency rate of the feeder line head end and end
The Amplitude Ration and phase difference absolute value of high-frequency signal, obtain n dimension Amplitude Ration vector sum n dimension phase difference vectors, n is power distribution network system
Feeder line number in system；
(3) for mixing feeder line, it is made up of if certain mixes feeder line N sections of different wave impedance circuits, then the wave impedance of the mixing feeder line
The number of discontinuity point is N+1, simulates singlephase earth fault at this N+1 point respectively, gathers the N+1 mixing feeder line head end
With the mode voltage traveling wave Mintrop wave head signal of end zero, and extract this N+1 times signal intermediate frequency rate for the amplitude of ω high-frequency signal and
Phase, calculates width of the zero mode voltage traveling wave Mintrop wave head signal intermediate frequency rate of N+1 signal head end and end for ω high-frequency signal
Value ratio and phase difference absolute value, obtain N+1 dimension Amplitude Ration vector sum N+1 dimension phase difference vectors；
(4) after actual power distribution network breaks down, each feeder line head end and the mode voltage row of end zero are extracted according to the method for step (1)
Ripple Mintrop wave head signal intermediate frequency rate is the amplitude and phase of ω high-frequency signal, and obtains physical fault according to the method for step (2)
N dimension Amplitude Ration vector sum n dimension phase difference vectors afterwards；
(5) ratio of the Amplitude Ration for the same feeder line that the Amplitude Ration that calculation procedure (2) is obtained is obtained with step (4), and step
(2) ratio of the phase difference absolute value for the same feeder line that the phase difference absolute value obtained is obtained with step (4), obtains n dimension amplitudes
Than ratio vector sum n dimension phase difference ratio vectors；
(6) find out the greatest member in n dimension Amplitude Ration ratios vector, n dimension phase difference ratio vectors respectively, if two it is vectorial most
Corresponding big element is same feeder line, then the feeder line is faulty line, if it is different that two vectorial greatest members are corresponding
Feeder line, the then corresponding feeder line of greatest member that n is tieed up to Amplitude Ration ratio vector is used as faulty line；
(7) if faulty line is mixing feeder line, the n that step (4) is obtained ties up Amplitude Ration vector sum n and ties up the event in phase difference vector
The corresponding element of barrier circuit obtains all elements in N+1 dimension Amplitude Ration vector sum N+1 dimension phase difference vectors with step (3) to be compared
Compared with, from N+1 tie up Amplitude Ration vector sum N+1 dimension phase difference vector in find out and meet the element of following formula：
SARm≥a≥SARm+1
SPDw≤p≤SPDw+1
In above formula, a is the corresponding element of faulty line in the n dimension Amplitude Ration vectors that step (4) is obtained, and p is what step (4) was obtained
The corresponding element of faulty line, SAR in n dimension phase difference vectorsmM-th of element in Amplitude Ration vector, SPD are tieed up for N+1wFor N+1
Tie up w-th of element in phase difference vector, m, w ∈ [1, N], if SARmWith SPDwCorresponding is same part of path, then the part of path
For the fault wire section of faulty line, if SARmWith SPDwCorresponding is different part of paths, then SARmCorresponding part of path is failure
The fault wire section of circuit；
(8) if faulty line is branch feeder, extract at each branch point of the branch feeder in zero mode voltage traveling wave Mintrop wave head signal
Frequency is the complex vector of ω high-frequency signal；
(9) energy of each mode voltage traveling wave Mintrop wave head signal intermediate frequency rate of branch point zero for the complex vector of ω high-frequency signal is calculated,
Then the maximum branch of energy is the fault branch of the faulty line.
2. the localization method of single-ended radiation type distribution network singlephase earth fault according to claim 1, it is characterised in that：Adopting
When collecting zero mode voltage traveling wave Mintrop wave head signal, it is thus necessary to determine that calculate data window, its determination method is as follows：
(a) catastrophe point of zero mould travelling wave signal is looked for, if zero line ripple signal vector of measurement point collection is U0=[u1,
u2,...,uk-1,uk,uk+1,...,uL], wherein L is signal total length, should be if the line ripple catastrophe point of failure zero is k simultaneously
Meet following two formula：
<mrow>
<mfrac>
<msub>
<mi>u</mi>
<mi>k</mi>
</msub>
<msub>
<mi>u</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>&GreaterEqual;</mo>
<mi>&epsiv;</mi>
</mrow>
<mrow>
<mfrac>
<msub>
<mi>u</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>u</mi>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>2</mn>
</mrow>
</msub>
</mfrac>
<mo><</mo>
<mi>&epsiv;</mi>
</mrow>
In above formula, uk-2,uk-1,ukIt is not zero, ε is the threshold value of setting；
(b) 100 sampled points are chosen before catastrophe point k and X sampled point is as calculating data window thereafter, wherein X is determined by following formula：
<mrow>
<mfrac>
<msub>
<mi>L</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mi>V</mi>
</mfrac>
<mo>&le;</mo>
<mi>X</mi>
<mo>&le;</mo>
<mn>2</mn>
<mfrac>
<msub>
<mi>L</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mi>V</mi>
</mfrac>
</mrow>
In above formula, LminFor the length of most short feeder line in all feeder lines in power distribution network, v is the mode voltage row of high frequency zero that frequency is ω
The velocity of wave of ripple.
3. the localization method of single-ended radiation type distribution network singlephase earth fault according to claim 2, it is characterised in that：It is described
Threshold epsilon=2.
4. the localization method of single-ended radiation type distribution network singlephase earth fault according to claim 1, it is characterised in that：Using
S-transformation extracts amplitude, phase and complex vector of the zero mode voltage traveling wave Mintrop wave head signal intermediate frequency rate for ω high-frequency signal.
5. the localization method of single-ended radiation type distribution network singlephase earth fault according to claim 1, it is characterised in that：Zero mould
Frequencies omega=45kHz of voltage traveling wave Mintrop wave head signal high frequency signal.
6. the localization method of single-ended radiation type distribution network singlephase earth fault according to claim 1, it is characterised in that：In step
Suddenly in (9), each mode voltage traveling wave Mintrop wave head signal intermediate frequency rate of branch point zero is the meter of the energy of the complex vector of ω high-frequency signal
Calculate formula as follows：
<mrow>
<mi>E</mi>
<mo>=</mo>
<msqrt>
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>A</mi>
<mi>m</mi>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>S</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>2</mn>
</msub>
</mrow>
</msqrt>
</mrow>
In above formula, S is the complex vector for the high-frequency signal that certain mode voltage traveling wave Mintrop wave head signal intermediate frequency rate of branch point zero is ω, and E is S
Energy, Amp (S) represents to seek S modulus value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710232520.9A CN107219440B (en) | 2017-04-11 | 2017-04-11 | The localization method of single-ended radiation type distribution network singlephase earth fault |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710232520.9A CN107219440B (en) | 2017-04-11 | 2017-04-11 | The localization method of single-ended radiation type distribution network singlephase earth fault |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107219440A true CN107219440A (en) | 2017-09-29 |
| CN107219440B CN107219440B (en) | 2018-11-27 |
Family
ID=59928216
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710232520.9A Active CN107219440B (en) | 2017-04-11 | 2017-04-11 | The localization method of single-ended radiation type distribution network singlephase earth fault |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107219440B (en) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108008246A (en) * | 2017-11-22 | 2018-05-08 | 湖南湘能智能电器股份有限公司 | Circuit fault distance measurement and device based on time-domain adaptive window |
| CN108445354A (en) * | 2018-05-23 | 2018-08-24 | 南方电网科学研究院有限责任公司 | Fault location method, device, equipment and medium based on traveling wave differential current |
| CN108693446A (en) * | 2018-05-25 | 2018-10-23 | 中国矿业大学 | A kind of Fault Locating Method of non-synchronous sampling power grid transient state travelling wave modulus time difference |
| CN108919054A (en) * | 2018-07-05 | 2018-11-30 | 国网上海市电力公司 | The intelligent power distribution network short circuit fault localization method estimated based on load and branch line |
| CN109188189A (en) * | 2018-07-11 | 2019-01-11 | 天津大学 | Based on arc characteristic ultra-high/extra-high voltage transmission line of electricity permanent fault recognition methods |
| CN109815593A (en) * | 2019-01-25 | 2019-05-28 | 周耘立 | A kind of analysis method of the distribution successive failure coupling mechanism of transmission |
| CN110133445A (en) * | 2019-06-10 | 2019-08-16 | 集美大学 | A kind of submarine cable fault distance-finding method, terminal device and storage medium |
| CN111007427A (en) * | 2019-11-23 | 2020-04-14 | 清华大学 | Distribution line single-phase earth fault line selection method and computer readable storage medium |
| CN111426908A (en) * | 2020-03-31 | 2020-07-17 | 国电南瑞科技股份有限公司 | Single-phase earth fault protection method, device and system for small current earthing system |
| CN112881864A (en) * | 2021-01-25 | 2021-06-01 | 云南电网有限责任公司昆明供电局 | Fault traveling wave acquisition device and method for medium-voltage power distribution network |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103926511A (en) * | 2014-03-17 | 2014-07-16 | 四川大学 | Distance measurement method for power distribution network ground fault based on zero-mode travelling wave difference |
| US20140300370A1 (en) * | 2013-04-08 | 2014-10-09 | Ecole Polytechnique Federale De Lausanne (Epfl) | Efficient Method Based on the Electromagnetic Time Reversal to Locate Faults in Power Network |
| CN105004969A (en) * | 2015-07-15 | 2015-10-28 | 四川大学 | Method for positioning single-phase grounding fault of overhead power distribution network |
| CN105203923A (en) * | 2015-09-18 | 2015-12-30 | 国家电网公司 | Power distribution network cable failure monitoring method based on traveling wave amplitude measure |
| CN106291261A (en) * | 2016-10-26 | 2017-01-04 | 四川大学 | The localization method of series-parallel connection one-phase earthing failure in electric distribution network |
-
2017
- 2017-04-11 CN CN201710232520.9A patent/CN107219440B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140300370A1 (en) * | 2013-04-08 | 2014-10-09 | Ecole Polytechnique Federale De Lausanne (Epfl) | Efficient Method Based on the Electromagnetic Time Reversal to Locate Faults in Power Network |
| CN103926511A (en) * | 2014-03-17 | 2014-07-16 | 四川大学 | Distance measurement method for power distribution network ground fault based on zero-mode travelling wave difference |
| CN105004969A (en) * | 2015-07-15 | 2015-10-28 | 四川大学 | Method for positioning single-phase grounding fault of overhead power distribution network |
| CN105203923A (en) * | 2015-09-18 | 2015-12-30 | 国家电网公司 | Power distribution network cable failure monitoring method based on traveling wave amplitude measure |
| CN106291261A (en) * | 2016-10-26 | 2017-01-04 | 四川大学 | The localization method of series-parallel connection one-phase earthing failure in electric distribution network |
Non-Patent Citations (2)
| Title |
|---|
| 梁睿 等: "单双端行波特征综合考虑的辐射状电网组合测距技术", 《高电压技术》 * |
| 梁睿 等: "复杂输电网中基于全网行波信息的测距装置最优配置", 《电工技术学报》 * |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108008246B (en) * | 2017-11-22 | 2020-04-21 | 湖南湘能智能电器股份有限公司 | Line fault distance measurement method and device based on time domain adaptive window |
| CN108008246A (en) * | 2017-11-22 | 2018-05-08 | 湖南湘能智能电器股份有限公司 | Circuit fault distance measurement and device based on time-domain adaptive window |
| CN108445354A (en) * | 2018-05-23 | 2018-08-24 | 南方电网科学研究院有限责任公司 | Fault location method, device, equipment and medium based on traveling wave differential current |
| CN108445354B (en) * | 2018-05-23 | 2020-05-08 | 南方电网科学研究院有限责任公司 | Fault location method, device, equipment and medium based on traveling wave differential current |
| CN108693446A (en) * | 2018-05-25 | 2018-10-23 | 中国矿业大学 | A kind of Fault Locating Method of non-synchronous sampling power grid transient state travelling wave modulus time difference |
| CN108919054A (en) * | 2018-07-05 | 2018-11-30 | 国网上海市电力公司 | The intelligent power distribution network short circuit fault localization method estimated based on load and branch line |
| CN109188189A (en) * | 2018-07-11 | 2019-01-11 | 天津大学 | Based on arc characteristic ultra-high/extra-high voltage transmission line of electricity permanent fault recognition methods |
| CN109815593A (en) * | 2019-01-25 | 2019-05-28 | 周耘立 | A kind of analysis method of the distribution successive failure coupling mechanism of transmission |
| CN109815593B (en) * | 2019-01-25 | 2023-07-25 | 周耘立 | Analysis method for coupling propagation mechanism of successive faults of distribution network |
| CN110133445A (en) * | 2019-06-10 | 2019-08-16 | 集美大学 | A kind of submarine cable fault distance-finding method, terminal device and storage medium |
| CN111007427A (en) * | 2019-11-23 | 2020-04-14 | 清华大学 | Distribution line single-phase earth fault line selection method and computer readable storage medium |
| CN111007427B (en) * | 2019-11-23 | 2021-05-04 | 清华大学 | Distribution line single-phase earth fault line selection method and computer readable storage medium |
| US11543462B2 (en) | 2019-11-23 | 2023-01-03 | Tsinghua University | Single-phase-to-ground fault line selection method for distribution lines and computer readable storage medium |
| CN111426908A (en) * | 2020-03-31 | 2020-07-17 | 国电南瑞科技股份有限公司 | Single-phase earth fault protection method, device and system for small current earthing system |
| CN111426908B (en) * | 2020-03-31 | 2022-07-08 | 国电南瑞科技股份有限公司 | Single-phase earth fault protection method, device and system for small current earthing system |
| CN112881864A (en) * | 2021-01-25 | 2021-06-01 | 云南电网有限责任公司昆明供电局 | Fault traveling wave acquisition device and method for medium-voltage power distribution network |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107219440B (en) | 2018-11-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107219440B (en) | The localization method of single-ended radiation type distribution network singlephase earth fault | |
| Naidu et al. | Precise traveling wave-based transmission line fault location method using single-ended data | |
| Daisy et al. | Single phase fault location in electrical distribution feeder using hybrid method | |
| CN105137281A (en) | Hybrid line fault point positioning method based on single-end electric quantity and transient traveling wave comprehensive characteristic analysis | |
| CN108802564B (en) | distribution network T-shaped line parameter-independent fault location algorithm and system | |
| CN108120899A (en) | A kind of single-ended Section Location of one-phase earthing failure in electric distribution network | |
| US11114843B2 (en) | Method and system for protection in a mixed line | |
| CN103592575A (en) | Self-adaptation weighting data fusion fault distance measurement method based on multi-sensor system | |
| CN104391229A (en) | Transmission line fault fast phase selection method based on S conversion | |
| CN113659547B (en) | Power distribution network differential protection data synchronization method and system based on effective zero crossing point | |
| CN109490706A (en) | A kind of multi branch electric power lines road Fault Locating Method | |
| CN106908692B (en) | A kind of adaptive reclosing judgment method of transmission line one-phase earth fault | |
| CN103116115A (en) | Network-distribution single-phase failure line selection method based on virtual active power increment maximal principle | |
| CN107817420A (en) | Non- whole parallel lines on same tower double back transmission line non-synchronous data fault distance-finding method | |
| CN109342888B (en) | Single-ended power transmission line fault location method based on reclosing action analysis | |
| CN107037324B (en) | Fault location method free from transition resistance based on single-end electric quantity | |
| CN106253244A (en) | A kind of based on electric current from the sense of current longitudinal protection method of structure reference quantity | |
| Bak et al. | High Voltage AC underground cable systems for power transmission–A review of the Danish experience: Part 2 | |
| Zahran et al. | Improved ground distance protection for cascaded overhead-submarine cable transmission system | |
| CN103675521B (en) | Protection test method for same-tower double-circuit direct current line of high-voltage direct current transmission system | |
| CN103424627B (en) | The method of double-end measurement parallel grid line zero-sequence impedance | |
| Yun et al. | Fault location method for three-terminal lines in distribution network based on line voltage measured by μMPMU | |
| CN102147443B (en) | Single-end distance measuring method based on self-adaptive current | |
| JP2018531575A (en) | Improvement of power system or related | |
| JP2018531575A6 (en) | Improvement of power system or related |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Address after: 221116 No. 1 University Road, copper mountain, Jiangsu, Xuzhou Applicant after: China University of Mining & Technology Address before: 221116 Xuzhou University Road, Jiangsu, No. 1 Applicant before: China University of Mining & Technology |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |