CN110133445A - A kind of submarine cable fault distance-finding method, terminal device and storage medium - Google Patents
A kind of submarine cable fault distance-finding method, terminal device and storage medium Download PDFInfo
- Publication number
- CN110133445A CN110133445A CN201910496512.4A CN201910496512A CN110133445A CN 110133445 A CN110133445 A CN 110133445A CN 201910496512 A CN201910496512 A CN 201910496512A CN 110133445 A CN110133445 A CN 110133445A
- Authority
- CN
- China
- Prior art keywords
- wave
- fault
- traveling wave
- submarine cable
- backward
- 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.)
- Pending
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02H—EMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
- H02H7/00—Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions
- H02H7/26—Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured
- H02H7/265—Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured making use of travelling wave theory
-
- 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/083—Locating faults in cables, transmission lines, or networks according to type of conductors in cables, e.g. underground
Abstract
The present invention relates to a kind of submarine cable fault distance-finding method, terminal device and storage mediums, in the method the following steps are included: S1: extracting the direct wave and backward-travelling wave when failure occurs in measurement point according to single end distance measurement method；S2: the time of the generation of the modulus maximum in direct wave and backward-travelling wave is calculated separately by wavelet transformation；S3: the position Lx of fault point is calculated.The present invention proposes the submarine cable fault distance-finding method combined with wavelet transformation and single end distance measurement, it is emulated in the submarine cable short trouble ranging model that MATLAB/Simulink is built, and fault traveling wave wave head information is identified using wavelet transformation and calculates modulus maximum, the one-end fault ranging precision based on wavelet transformation is verified under different fault distances and fault type, measurement result shows the lower measurement error of this method, therefore, measurement accuracy with higher.
Description
Technical field
The present invention relates to submarine cable ranging detection field more particularly to a kind of submarine cable fault distance-finding methods, terminal
Equipment and storage medium.
Background technique
As the exploitation on remote island increases and the continuous expansion of ocean development of resources scale, electrical energy demands over strait are sharply
Increase, the submarine cable for needing to be laid with is more and more, and increasing for mankind's Activities of Ocean causes submarine cable to be caught by drag anchor, fishery
The breakoff phenomenon of the factors such as fishing, ship dragging, bank base operation frequently occurs, and when submarine cable is damaged, will necessarily result in cable
Abnormal electrical power supply even cannot continue to power.Therefore, it is necessary to carry out ranging, positioning to fault point in time, and then repair failure and extensive
Multiple power supply.And the accuracy of ranging largely affects the time used in fault restoration.
Measuring distance of transmission line fault is divided into impedance method and traveling wave method by principle, due to impedance method ranging and transmission line of electricity
Parameter distribution is closely related, and the parameter that theoretically meet transmission line of electricity is uniformly distributed, but marine environment locating for submarine cable and
Landform is sufficiently complex, and parameter distribution otherness is larger, is not suitable for impedance method ranging.And common traveling wave method is due to subsea environment
Complicated reason, detection accuracy be not high.
Summary of the invention
In view of the above-mentioned problems, the present invention is intended to provide a kind of submarine cable fault distance-finding method, terminal device and storage are situated between
Matter improves detection accuracy for subsea environment.
Concrete scheme is as follows:
A kind of submarine cable fault distance-finding method, comprising the following steps:
S1: the direct wave and backward-travelling wave when failure occurs are extracted in measurement point according to single end distance measurement method；
S2: the time of the generation of the modulus maximum in direct wave and backward-travelling wave is calculated separately by wavelet transformation；
S3: the position Lx of fault point is calculated:
Wherein, t1For the time that the modulus maximum of direct wave occurs, t2For backward-travelling wave modulus maximum occur when
Between, v is the spread speed of traveling wave.
Further, the detailed process of step S1 are as follows:
S11: the three-phase voltage current value before the three-phase voltage current value after failure to be subtracted to failure obtains three-phase electricity piezoelectricity
The transient of stream；
S12: the transient of three-phase voltage current is converted to the Aerial mode component and zero _exit of voltage and current, and selects certain
A mold component is as travelling wave signal to be measured；
S13: the direct wave and backward-travelling wave of travelling wave signal to be measured are calculated.
Further, it is converted in step S12 by Clarke transform, specific conversion formula are as follows:
Wherein, ua、ub、ucThree-phase voltage traveling-wave component respectively on transmission line of electricity；uα、uβ、u0Respectively voltage traveling wave
α, β, 0 mold component；ia、ib、icFor the current traveling wave component on transmission line of electricity；iα、iβ、i0α, β, 0 mould point for current traveling wave
Amount.
Further, using α mold component as travelling wave signal to be measured in step S12.
Further, the calculation method of the direct wave of fault-signal and backward-travelling wave are as follows:
Wherein, S1α、S1β、S10Respectively α, β of direct wave, 0 mold component；S2α、S2β、S20The respectively α of backward-travelling wave,
β, 0 mold component；Zα、Zβ、Z0For the corresponding wave impedance of mold component traveling wave.
Further, the wavelet transformation is db4 small echo.
A kind of submarine cable fault localization terminal device, including processor, memory and storage are in the memory
And the computer program that can be run on the processor, the processor realize that the present invention is real when executing the computer program
The step of applying example above-mentioned method.
A kind of computer readable storage medium, the computer-readable recording medium storage have computer program, feature
The step of being, above-mentioned method of the embodiment of the present invention realized when the computer program is executed by processor.
The present invention uses technical solution as above, proposes the submarine cable fault localization combined with wavelet transformation and single end distance measurement
Method is emulated in the submarine cable short trouble ranging model that MATLAB/Simulink is built, and is known using wavelet transformation
Other fault traveling wave wave head information simultaneously calculates modulus maximum, and verifying is based on wavelet transformation under different fault distances and fault type
One-end fault ranging precision, measurement result shows the lower measurement error of this method, therefore, measurement accuracy with higher.
Detailed description of the invention
Fig. 1 show the fault traveling wave production principle figure in the embodiment of the present invention one.
Fig. 2 show the Single Terminal Traveling Wave Fault Location schematic diagram in the embodiment.
Fig. 3 show the both-end method travelling wave ranging schematic diagram in the embodiment.
Fig. 4 show the voltage traveling wave catadioptric figure in the embodiment.
Fig. 5 show the current traveling wave catadioptric figure in the embodiment.
Fig. 6 show the Morlet wavelet function in the embodiment.
Fig. 7 show the waveform diagram of the db4 small echo in the embodiment.
Fig. 8 show the analogue system figure in the embodiment.
Fig. 9 show the A phase ground short circuit false voltage waveform diagram in the embodiment.
Figure 10 show the A phase ground short circuit fault current waveform figure in the embodiment.
Figure 11 show the A phase ground fault current traveling wave waveform diagram in the embodiment.
The alternate ground short circuit voltage oscillogram of the AB that Figure 12 is shown in the embodiment.
The alternate earth short circuit current waveform diagram of the AB that Figure 13 is shown in the embodiment.
The alternate earth short circuit current traveling wave figure of the AB that Figure 14 is shown in the embodiment.
Figure 15 show three relative ground circuit voltage oscillograms in the embodiment.
Figure 16 show the three-phase short circuit current to earth waveform diagram in the embodiment.
Figure 17 show the three-phase short circuit current to earth traveling wave figure in the embodiment.
Figure 18 show the fault point travelling wave analysis figure in the embodiment.
Figure 19 show another fault point travelling wave analysis figure in the embodiment.
Figure 20 show the wavelet transformation of the end the singlephase earth fault M electric current direct wave in the embodiment.
Figure 21 show the wavelet transformation of the end the singlephase earth fault M electric current backward-travelling wave in the embodiment.
Specific embodiment
To further illustrate that each embodiment, the present invention are provided with attached drawing.These attached drawings are that the invention discloses one of content
Point, mainly to illustrate embodiment, and the associated description of specification can be cooperated to explain the operation principles of embodiment.Cooperation ginseng
These contents are examined, those of ordinary skill in the art will be understood that other possible embodiments and advantages of the present invention.
Now in conjunction with the drawings and specific embodiments, the present invention is further described.
Embodiment one:
Measuring distance of transmission line fault is divided into impedance method and traveling wave method by principle, due to impedance method ranging and transmission line of electricity
Parameter distribution is closely related, and the parameter that theoretically meet transmission line of electricity is uniformly distributed, but marine environment locating for submarine cable and
Landform is sufficiently complex, and parameter distribution otherness is larger, is not suitable for impedance method ranging.And travelling wave ranging reliability and precision are not
It is influenced by route distribution parameter unevenness, is suitable for submarine cable fault localization.
Hot spot using the Research on fault locating of traveling wave method is one-end fault ranging and Two-terminal Fault Location, and both-end failure is surveyed
Although away from the influence for being avoided that traveling wave catadioptric in communication process causes measurement accuracy, its equipment cost and computation complexity
It is higher than single-ended method ranging.
Since the superiority that wavelet transformation is showed when extracting local signal Characteristics of Mutation is widely used in electric system
The fields such as perturbation analysis and electromagnetic transient analysis can effectively extract local signal feature.Therefore it proposes in the embodiment based on small
Submarine cable fault distance-finding method of the wave analysis in conjunction with single-ended method.
1, fault travelling wave ranging is theoretical
(1) basic conception of traveling wave
Traveling wave is the transient voltage generated when transmission line of electricity breaks down, current signal.Fault Location With Traveling Wave is to apply at present
Fault distance-finding method the most extensive.When certain point A of submarine cable route breaks down, can be analyzed using principle of stacking,
It is the failure additivity figure that fault component figure can be equivalent to Fig. 1 (c) if Fig. 1 (a) is normal load component figure and Fig. 1 (b).
It is equivalent to by the visible fault component of Fig. 1 (b) when system potential is zero, in fault point, A goes out to add one and this normal duty state
Under equal in magnitude, contrary voltage.Under the action of this voltage, the fault traveling wave phase bus both ends that fault point generates are passed
It broadcasts.
By on single lossless sea cable transmission line of electricity voltage u and electric current i on the line position x and time t be change
Several partial differential equation indicate, then have:
Wherein, L and C is respectively the inductance and capacitor of route unit length.
By formula (1) and formula (2) respectively to x, t differential, the wave equation of characterization transmission line travelling wave, fluctuation side can be obtained
Journey expression formula is as follows:
Then it is up to Lang Beier (D ' Alembert) solution
Wherein,For x positive direction propagate the positive traveling wave of failure,For along the anti-row of the failure of x opposite direction
Wave,It is the spread speed of fault traveling wave,The wave impedance of fault traveling wave.
In the submarine cable transmission line of electricity of three-phase, due to there is coupling between each phase, lead to the traveling-wave component of every phase not
It can be mutually indepedent.Therefore, it is necessary to carry out phase-model transformation to traveling-wave component, convert the phase component that three mutually not only stand to mutually solely
Then vertical phase component recycles modulus traveling-waves to realize the corresponding function of traveling wave.The method of phase-model transformation is mainly Clarke
(Clarke) transformation and the transformation of triumphant human relations Bell (Karenbauer) are realized, since triumphant human relations Bell transformation can not achieve according to single amount
Judge failure phase mould, in the embodiment realize that phase-model transformation, transform method are as follows in Clarke transform method:
Wherein, ua、ub、ucThree-phase voltage traveling-wave component respectively on transmission line of electricity；uα、uβ、u0Respectively voltage traveling wave
α, β, 0 mold component；ia、ib、icFor the current traveling wave component on transmission line of electricity；iα、iβ、i0α, β, 0 mould point for current traveling wave
Amount.
Therefore, the mold component of direction traveling wave may be expressed as:
Wherein, S1α、S1β、S10Respectively α, β of direct wave, 0 mold component；S2α、S2β、S20The respectively α of backward-travelling wave,
β, 0 mold component；Zα、Zβ、Z0For the corresponding wave impedance of mold component traveling wave.
(2) travelling wave ranging principle
The travelling wave ranging method of overhead line has just been proposed early in the 1930s, i.e., event is occurred by measurement transmission line of electricity
The wavefront generated when barrier or reflection, refracted wave wave head reach the time of measurement point to calculate between fault point and measurement end
Distance.Travelling wave ranging scheme is attributed to tetra- seed type of A, B, C, D: A type is to reach bus back reflection according to the traveling wave of fault point to arrive
Fault point, then fault distance is calculated by the time difference of arrival bus after failure point reflection；Type B and c-type include pulse or signal
Generator is artificially to apply high frequency or direct current signal after failure, is made according to Principles of Radar that wherein B is both-end method, and c-type is single
Hold-carrying；D type is that two side bus of initial row wave direction generated using fault point moves to the time difference of two side ends bus to determine
The distance between fault point and measurement end.A type and D type have preferable be applicable in for the instantaneity and permanent fault of route
Property, and Type B and c-type are only applicable to permanent fault.Since transmission line of electricity is mostly high voltage power transmission feature, artificial poor signal is applied
Add, and excessively expense is prohibitively expensive, Type B and c-type are poorly suited for the fault localization of the transmission line of electricity of high pressure.
Single Terminal Traveling Wave Fault Location be using traveling wave for the first time reach measurement end and its reflect back into side both ends from fault point when
Between difference calculate the distance between fault point and measurement end, belong to A type travelling wave ranging method.As shown in Fig. 2, in the F of transmission line of electricity
Point breaks down, and fault traveling wave is t at the time of reaching the end bus M for the first time1, fault traveling wave is through the reflection of the end bus M again by failure
Second of point reflection is t at the time of reaching the end bus M2, the overall length of transmission line of electricity is L, and the spread speed of traveling wave (is in theoretical value
The light velocity) it is v, the absolute moment that failure occurs is t0.Therefore, the equation group that can obtain one-end fault travelling wave ranging is as follows:
Above-mentioned equation group is solved to obtain:
The key of one-end fault travelling wave ranging is accurately to find out fault traveling wave to reach side both ends for the first time with it by reflecting
The time difference of measurement end is returned to, accurately can extract traveling-wave component characteristic information is the principal element for influencing range accuracy.
Both-end travelling wave ranging is to reach the time of measuring difference at route both ends for the first time using fault traveling wave to realize that failure is surveyed
Away from as shown in Figure 3.Transmission line of electricity is pointed out in F and is broken down, and the time that fault traveling wave reaches the end bus M is t1, fault traveling wave arrives
Time up to bus N-terminal is t2, the overall length of transmission line of electricity is L, and the absolute moment that failure occurs is t0, both-end fault traveling wave can be obtained
Equation group it is as follows:
Solution obtains the distance that traveling wave reaches the end bus M:
Therefore, the time that detection transmission line malfunction traveling wave reaches both ends bus can also measure fault point and measurement end
Fault distance, with the development of GPS synchronised clock technology, the precision of fault traveling wave both-end distance measuring is also continuously improved, but both-end is surveyed
Away from technical requirements and cost it is more much higher than single end distance measurement.
(3) reflection of traveling wave and refracting feature are analyzed
Voltage traveling wave and current traveling wave meet the relationship of impedance Z, but submarine cable in uniform transmission line of electricity upload sowing time
In complicated laying environment, including the rough topography such as coral reef and clay causes line stress to be unevenly distributed, and transmits electricity
The distribution parameter of route changes, i.e. the ratio wave impedance of voltage traveling wave and current traveling wave during moving ahead becomes
Change, the reflection and refraction of traveling wave will occur in the changed place of wave impedance.
If Fig. 4 and Fig. 5 are respectively the catadioptric figure of voltage, current traveling wave, when the voltage and current row that voltage magnitude is U and I
Wave impedance value is Z1Transmission line of electricity is propagated forward, and impedance changes at F point, and the current traveling wave of voltage will occur instead here
It penetrates and reflects, be Z in wave impedance1Transmission line of electricity on the voltage and current traveling wave that reflects be respectively uf2And if2, in wave resistance
Resist for Z2Transmission line of electricity on refraction voltage and current traveling wave be respectively uq2And iq2.At the F point that wave impedance mutates, electricity
Stream and voltage traveling wave have continuity, so that
Meet following relationship between voltage, current traveling wave and wave impedance:
Voltage, current traveling wave transmission process in introduce voltage traveling wave folding, the folding of reflection coefficient and current traveling wave, reflection
Coefficient is anti-other are as follows: αu, βu, αi, βi。
It can be obtained by formula (15) and formula (16): uq2=αuU, uf1=βuU；iq2=αiI, if1=βiI.Voltage can be calculated as a result,
Traveling wave catadioptric coefficient and current traveling wave catadioptric coefficient:
The variation of the catadioptric coefficient of transmission line malfunction traveling wave is as caused by surge impedance of a line.Therefore, submarine cable
The environment of laying and the normalization of construction will affect the variation of wave impedance, and then have an impact to fault travelling wave ranging precision.
2, wavelet transformation
(1) wavelet transformation and its application
There are Heisenberg's uncertainties for signal analysis due to Fourier transformation and Short Time Fourier Transform to non-stationary, i.e.,
Frequency domain and time-domain analysis cannot be met simultaneously.Although Short Time Fourier Transform improves the part of Fourier transformation with adding window method
Property analysis it is insufficient, but narrow width of window function also will affect the accuracy that signal is analyzed.When window function is wider, frequency can be improved
Resolution ratio, but temporal resolution can be reduced；When window function is relatively narrow, temporal resolution can be improved, but frequency resolution can reduce.
Then French petroleum signal processing engineer J.Morle proposed wavelet transformation in 1974, it overcomes Fourier transformation
The shortcomings that.Small echo is that effective expanded range is limited and telescopic damped oscillation function, be illustrated in figure 6 Morlet small echo when
Domain and frequency domain figure.The essence of wavelet transformation is the local similarity for calculating analysis signal and wavelet shapes, i.e. analysis jump signal
Characteristic information at mutation, this is also the important process of measuring distance of transmission line fault.
(2) Fundamentals of Mathematics of wavelet transformation
If function ψ ∈ L2(R)∩L1(R), andFamily's function is obtained by flexible and translation by ψ
Claim ψ (t)a,bFor analysis wavelet or continuous wavelet, ψ is referred to as wavelet or morther wavelet.Wherein a is contraction-expansion factor, b
For shift factor.By conditionIt can obtain,
For given wavelet ψ, in space L2(R) continuous wavelet transform of signal f (t) is defined as:
Wf(a, b)=∫Rf(t)ψa,b(t)dt (19)
If wavelet meets following condition:
Then there is the inverse transformation of continuous wavelet, inverse transformation formula are as follows:
Since the data obtained in practical operation are mostly discrete value, therefore, it is necessary to by continuous wavelet discretization.When continuous
When wavelet discrete, contraction-expansion factor a and shift factor b also must be discretizations, and the discretization for studying wavelet transformation reforms into
Parameter discretization, and in time variable t is not comprised in, therefore, stretch because there are following discrete forms by a and shift factor b:
Continuous wavelet can be obtained are as follows:
It is defined as follows after discretization:
DWf(m, n)=∫Rf(t)ψm,n(t)dt (26)
By above-mentioned demonstration, wavelet transform DWfWhen (m, n) restores original signal, reconstruction formula are as follows:
(3) singular points detection and the modulus maximum to small echo are theoretical
When detecting the singular point of transmission line malfunction traveling wave, need to be smoothed travelling wave signal, and calculate signal
First derivative or second dervative, and then using first derivative maximum and second dervative zero point detect fault-signal it is unusual
Point.The smooth function θ (t) that wavelet transformation utilizes is defined as follows:
Smooth function meet it is secondary can be micro-, first derivative and second dervative are asked to above-mentioned smooth function:
Wavelet transformation of the function f (t) at scale S and variable t is that " * " indicates convolution:
Calculable functions f (t) is in the single order derived function and second order derived function after scale S smoothing processing as a result:
When signal smoothing treated second dervativeWhen there is zero, show that first derivative obtains extreme value；When signal is flat
Sliding treated first derivativeWhen taking extreme value, it is the strongest to indicate that signal changes at this.
As wavelet transformation zero crossing T existing for the first derivative under scale S of function f (t)0The referred to as office of wavelet transformation
Portion's extreme value, in T0A certain field δ in all in the presence of other any points wavelet transformation absolute value all than in T0The absolute value at place
It is small, then T0For the modulus maximum point of wavelet transformation, the value of the point is modulus maximum.
(4) selection of wavelet basis
In the short trouble ranging emulation of submarine cable, fault message is collectively constituted by high frequency section and low frequency part,
Therefore to consider that morther wavelet can characterize the Characteristics of Mutation of place's signal when choosing wavelet basis, and can reflect that high frequency is believed simultaneously
Breath and low-frequency information, and Daubechies (dbN) small echo is extracted with good abrupt local information characteristics, additionally it is possible to accurate weight
The local feature information of structure signal.Db4 small echo is selected to identify for submarine cable fault traveling wave wave head in the embodiment.
Daubechies (dbN) small echo is proposed by the wavelet analysis scholar Inrid Daubechies of France, small
DbN is abbreviated as in wave analysis, wherein N is the order of small echo, which does not have specific expression formula, wavelet function ψ and scale letter
The effective length of number φ is 2N-1, the as Haar small echo when N takes 1.Fig. 7 is the time-frequency waveform diagram of db4.Small echo dbN small echo tool
There are following 4 features:
It (1) is finite support in time domain, i.e. ψ (t) is finite length；
(2) in frequency domain ψ (ω), there is N rank zero point at ω=0；
(3) scaling function φ (t) is low-pass filter function, and finite length, supporting domain is in the range of t=0~(2N-1)；
(4) ψ (t) and its integer are displaced orthonomalization, i.e. ∫ ψ (t) ψ (t-k) dt=δk, wherein k is integer displacement.
3, fault localization simulation analysis
Submarine cable short trouble is built using the Powersystem in MTALAB/Simulink platform in the embodiment
Ranging simulation model respectively imitates single-phase grounding fault, two-phase phase fault and three-phase ground short trouble
Very, analogue system figure is as shown in figure 8, line construction parameter are as follows: R0=0.1148 Ω/km, R1=0.02083 Ω/km；L0=
0.0022886H/km, L1=0.0008984H/km；C0=0.00523uF/km, C1=0.01294uF/km.Sample frequency is
500kHz, corresponding line mould speed are v1=293516km/s, zero mould speed v2=298170km/s.It is small using single-ended measurement and db4
Wave carries out 6 layers of decomposition to travelling wave signal α mold component.
Now with A phase ground short circuit failure, AB phase phase fault and ABC tri-, relatively shorted to earth failure is carried out respectively
The step of emulating and write Program extraction current traveling wave using MATLAB, extracting the direct wave and backward-travelling wave when failure occurs
It is as follows:
(1) transient for extracting three-phase voltage current subtracts event with the three-phase voltage current value in a period of time after failure
Three-phase voltage current value before barrier in corresponding a period of time, just obtains the transient of three-phase voltage current；
(2) transient of three-phase voltage current is subjected to Clarke (Clark) transformation using formula (7) and formula (8) and obtains electricity
The Aerial mode component and zero _exit of current voltage；
(3) it recycles formula (9) and formula (10) to calculate the direct wave and backward-travelling wave of fault-signal, considers in the embodiment
The influence that the wave head of traveling wave is extracted in the reflected refraction of fault-signal is reduced, α mold component is selected to believe as the traveling wave of fault localization
Number.
The MATLAB principal function for extracting the fault traveling wave of submarine cable singlephase earth fault is as follows:
clc
clear
load xingbo.mat；% is loaded into .mat file
M=n'；
Ua=m (3501:2500,2)-m (1501:2900,2)；
Ia=m (3501:4900,5)-m (1501:2900,5)；
Ub=m (3501:4900,3)-m (1501:2900,3)；
Ib=m (3501:4900,6)-m (1501:2900,6)；
Uc=m (3501:4900,4)-m (1501:2900,4)；
Ib=m (3501:4900,7)-m (1501:2900,7)；
Q=1/3* [2-1-1；0sqrt(3)sqrt(3)；1 1 1]；
Um1=Q (1 :) * [ua ub uc] '；
Im1=Q (1 :) * [ia ib ic] '；% carries out Clarke and converts to obtain the modulus of voltage, electric current
Lm1=0.0894e-3；
Cm1=12.94e-9；
Zcm1=sqrt (Lm1/Cm1)；% seeks wave impedance
Uf=(um1+im1*Zcm1)；
Ur=(um1-im1*Zcm1)；
IF1=uf/Zcm；
IF2=ur/Zcm；
T1=0:10:4990；
T=t1'；
plot(t,IF1,'r',t,IF2,'b--')；
xlabel('t/us')；Ylabel (' electric current/A')；
Legend (' direct wave ', ' backward-travelling wave ', ' location', ' northwest')；
Meet hair by the voltage waveform and current waveform of current-voltage measurement module and oscilloscope module fault inspecting
The voltage and current feature of raw single-phase grounding fault, two-phase phase fault and three relative ground circuit failures, to illustrate to emulate
The correctness of model.
It carries out carrying out inflection point detection in fault traveling wave using wavelet transformation tool, extracts direct wave and backward-travelling wave
The wave head time and modulus maximum is found out using MATLAB algorithm.By taking A phase ground fault as an example, fault point is observed from M measurement end
Traveling wave, such as Figure 18,1 be the initial traveling wave of failure；2 be fault point back wave, and corresponding traveling wave path is F-M-F-M；3 be failure
Point refracted wave；4 be opposite end bus refracted wave；5 be the back wave of opposite end N, and corresponding traveling wave path is F-N-F-M.Such as Figure 19,1'
For the transmitted wave of the initial traveling wave in fault point；The transmitted wave of the fault point 2' back wave, corresponding traveling wave path are F-M-F-M；3' is
The projection wave of measurement end M；4' is the back wave for measuring the end M；5' is the transmitted wave of opposite end N, and corresponding traveling wave path is F-N-F-
M。
Such as the forward and reverse traveling wave of voltage failure traveling wave of Figure 20 and Figure 21 (abscissa is time/ms, and ordinate is electric current/mA)
Wavelet transformation in the failure catastrophe point of traveling wave can be observed in scale d1.Now connected with the A of single-phase grounding fault
For ground short circuit failure, occurs for setting failure for 0.15ms the moment, and fault point is 10.2km, traveling wave speed choosing at a distance from the end M
With above-mentioned line mould speed, hair is calculated at the time of generation by wavelet transformation modulus maximum and using single-ended measure equation (12)
Raw single-phase grounding fault point is at a distance from the end M are as follows:
Measuring relative error is 0.69%.It is verifying different faults distance and different faults type to utilization wavelet transformation
The influence of one-end fault ranging, it is shown in Fig. 8 on the basis of by setting different faults distance and type approval this method it is accurate
Degree.Such as the single-ended measurement range error of wavelet transformation that table 1 is different faults distance, table 2 is the wavelet transformation of different faults type
Single-ended measurement range error.By the one-end fault ranging method based on wavelet transformation in different faults distance and different faults type
Simulation result shows that the measurement error of this method is held in the range of 1%, measurement accuracy with higher.
Table 1
Table 2
4, conclusion
The theory of travelling wave of transmission line of electricity is introduced in the embodiment, the two methods in traveling wave fault location: one-end fault is surveyed
Away from Two-terminal Fault Location range measurement principle and fault localization Affecting Factors of Accuracy analysis, and analyze compare two methods superiority and inferiority
Property, double-end distance measurement method is by the complicated limitation of at high cost and calculating, and single-ended method is at low cost, calculating process is simple.But due in single-ended event
It is influenced by fault traveling wave wave head comprising complicated high and low frequency information and catadioptric when hindering ranging, traditional Fourier transformation side
Method can not accurately identify and extract fault characteristic information, and wavelet transformation has good local feature information extraction with high accuracy
Feature, the method for proposing in the embodiment to be combined with wavelet transformation and single end distance measurement.In the seabed that MATLAB/Simulink is built
Cable short trouble ranging model emulation simultaneously identifies fault traveling wave wave head information using wavelet transformation and calculates modulus maximum, not
The one-end fault ranging precision based on wavelet transformation is verified under same fault distance and fault type, measurement result shows this method
Lower measurement error.
Embodiment two:
The present invention also provides a kind of submarine cable fault localization terminal device, including memory, processor and it is stored in
In the memory and the computer program that can run on the processor, when the processor executes the computer program
Realize the step in the above method embodiment of the embodiment of the present invention one.
Further, as an executable scheme, the submarine cable fault localization terminal device can be desktop
Computer, notebook, palm PC and cloud server etc. calculate equipment.The submarine cable fault localization terminal device can wrap
It includes, but is not limited only to, processor, memory.It will be understood by those skilled in the art that above-mentioned submarine cable fault localization terminal is set
Standby composed structure is only the example of submarine cable fault localization terminal device, is not constituted whole to submarine cable fault localization
The restriction of end equipment may include perhaps combining certain components or different components than above-mentioned more or fewer components,
Such as the submarine cable fault localization terminal device can also be including input-output equipment, network access equipment, bus etc., this
Inventive embodiments do not limit this.
Further, as an executable scheme, alleged processor can be central processing unit (Central
Processing Unit, CPU), it can also be other general processors, digital signal processor (Digital Signal
Processor, DSP), it is specific integrated circuit (Application Specific Integrated Circuit, ASIC), existing
At programmable gate array (Field-Programmable Gate Array, FPGA) or other programmable logic device, discrete
Door or transistor logic, discrete hardware components etc..General processor can be microprocessor or the processor can also
To be any conventional processor etc., the processor is the control centre of the submarine cable fault localization terminal device, benefit
With the various pieces of various interfaces and the entire submarine cable fault localization terminal device of connection.
The memory can be used for storing the computer program and/or module, and the processor is by operation or executes
Computer program in the memory and/or module are stored, and calls the data being stored in memory, described in realization
The various functions of submarine cable fault localization terminal device.The memory can mainly include storing program area and storing data
Area, wherein storing program area can application program needed for storage program area, at least one function；Storage data area can store
Created data etc. are used according to program.In addition, memory may include high-speed random access memory, can also include
Nonvolatile memory, such as hard disk, memory, plug-in type hard disk, intelligent memory card (Smart Media Card, SMC), safety
Digital (Secure Digital, SD) card, flash card (Flash Card), at least one disk memory, flush memory device or
Other volatile solid-state parts.
The present invention also provides a kind of computer readable storage medium, the computer-readable recording medium storage has computer
Program, when the computer program is executed by processor the step of the realization above method of the embodiment of the present invention.
If the integrated module/unit of the submarine cable fault localization terminal device is real in the form of SFU software functional unit
Now and when sold or used as an independent product, it can store in a computer readable storage medium.Based in this way
Understanding, the present invention realize above-described embodiment method in all or part of the process, can also be instructed by computer program
Relevant hardware is completed, and the computer program can be stored in a computer readable storage medium, the computer program
When being executed by processor, it can be achieved that the step of above-mentioned each embodiment of the method.Wherein, the computer program includes computer
Program code, the computer program code can be source code form, object identification code form, executable file or certain centres
Form etc..The computer-readable medium may include: can carry the computer program code any entity or device,
Recording medium, USB flash disk, mobile hard disk, magnetic disk, CD, computer storage, read-only memory (ROM, ROM, Read-Only
Memory), random access memory (RAM, Random Access Memory) and software distribution medium etc..
Although specifically showing and describing the present invention in conjunction with preferred embodiment, those skilled in the art should be bright
It is white, it is not departing from the spirit and scope of the present invention defined by the appended claims, it in the form and details can be right
The present invention makes a variety of changes, and is protection scope of the present invention.
Claims (8)
1. a kind of submarine cable fault distance-finding method, it is characterised in that: the following steps are included:
S1: the direct wave and backward-travelling wave when failure occurs are extracted in measurement point according to single end distance measurement method；
S2: the time of the generation of the modulus maximum in direct wave and backward-travelling wave is calculated separately by wavelet transformation；
S3: the position Lx of fault point is calculated:
Wherein, t1For the time that the modulus maximum of direct wave occurs, t2For the time that the modulus maximum of backward-travelling wave occurs, v is
The spread speed of traveling wave.
2. according to the method described in claim 1, it is characterized by: the detailed process of step S1 are as follows:
S11: the three-phase voltage current value before the three-phase voltage current value after failure to be subtracted to failure obtains three-phase voltage current
Transient；
S12: the transient of three-phase voltage current is converted to the Aerial mode component and zero _exit of voltage and current, and selects some mould
Component is as travelling wave signal to be measured；
S13: the direct wave and backward-travelling wave of travelling wave signal to be measured are calculated.
3. according to the method described in claim 1, it is characterized by: being converted in step S12 by Clarke transform, specifically
Conversion formula are as follows:
Wherein, ua、ub、ucThree-phase voltage traveling-wave component respectively on transmission line of electricity；uα、uβ、u0The respectively α of voltage traveling wave,
β, 0 mold component；ia、ib、icFor the current traveling wave component on transmission line of electricity；iα、iβ、i0α, β, 0 mold component for current traveling wave.
4. according to the method described in claim 2, it is characterized by: using α mold component as travelling wave signal to be measured in step S12.
5. according to the method described in claim 3, it is characterized by: the direct wave of fault-signal and the calculating side of backward-travelling wave
Method are as follows:
Wherein, S1α、S1β、S10Respectively α, β of direct wave, 0 mold component；S2α、S2β、S20Respectively α, β of backward-travelling wave, 0
Mold component；Zα、Zβ、Z0For the corresponding wave impedance of mold component traveling wave.
6. according to the method described in claim 1, it is characterized by: the wavelet transformation is db4 small echo.
7. a kind of submarine cable fault localization terminal device, it is characterised in that: including processor, memory and be stored in described
The computer program run in memory and on the processor, the processor are realized such as when executing the computer program
In claim 1~6 the step of any the method.
8. a kind of computer readable storage medium, the computer-readable recording medium storage has computer program, and feature exists
In when the computer program is executed by processor in realization such as claim 1~6 the step of any the method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910496512.4A CN110133445A (en) | 2019-06-10 | 2019-06-10 | A kind of submarine cable fault distance-finding method, terminal device and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910496512.4A CN110133445A (en) | 2019-06-10 | 2019-06-10 | A kind of submarine cable fault distance-finding method, terminal device and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110133445A true CN110133445A (en) | 2019-08-16 |
Family
ID=67580835
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910496512.4A Pending CN110133445A (en) | 2019-06-10 | 2019-06-10 | A kind of submarine cable fault distance-finding method, terminal device and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110133445A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111948493A (en) * | 2020-08-21 | 2020-11-17 | 兰州理工大学 | MMC-HVDC direct current transmission line fault positioning method |
| CN113030648A (en) * | 2021-03-19 | 2021-06-25 | 国网河北省电力有限公司邢台供电分公司 | Method and device for determining position of fault point of power cable and terminal equipment |
| CN115267417A (en) * | 2022-06-15 | 2022-11-01 | 北京妙微科技有限公司 | Accurate positioning method for transmission line fault and transmission line traveling wave measuring device |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0572249A (en) * | 1991-09-10 | 1993-03-23 | Furukawa Electric Co Ltd:The | Method for locating fault point of transmission line |
| CN102135588A (en) * | 2011-02-18 | 2011-07-27 | 昆明理工大学 | Intelligent ranging method for radiation network fault by utilizing S transformation spectrum characteristics |
| CN104569744A (en) * | 2014-11-26 | 2015-04-29 | 国家电网公司 | Comprehensive single-end fault positioning method applicable to power distribution network lines |
| 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 |
| CN105738760A (en) * | 2014-12-12 | 2016-07-06 | 国家电网公司 | Frequency domain method and traveling wave method-combined high-resistance fault location method |
| CN107219440A (en) * | 2017-04-11 | 2017-09-29 | 中国矿业大学 | The localization method of single-ended radiation type distribution network singlephase earth fault |
| CN109725236A (en) * | 2019-02-26 | 2019-05-07 | 四川轻化工大学 | A kind of bus-bar fault area determination method based on power frequency quantity Polarity comparision |
-
2019
- 2019-06-10 CN CN201910496512.4A patent/CN110133445A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0572249A (en) * | 1991-09-10 | 1993-03-23 | Furukawa Electric Co Ltd:The | Method for locating fault point of transmission line |
| CN102135588A (en) * | 2011-02-18 | 2011-07-27 | 昆明理工大学 | Intelligent ranging method for radiation network fault by utilizing S transformation spectrum characteristics |
| CN104569744A (en) * | 2014-11-26 | 2015-04-29 | 国家电网公司 | Comprehensive single-end fault positioning method applicable to power distribution network lines |
| CN105738760A (en) * | 2014-12-12 | 2016-07-06 | 国家电网公司 | Frequency domain method and traveling wave method-combined high-resistance fault location method |
| 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 |
| CN107219440A (en) * | 2017-04-11 | 2017-09-29 | 中国矿业大学 | The localization method of single-ended radiation type distribution network singlephase earth fault |
| CN109725236A (en) * | 2019-02-26 | 2019-05-07 | 四川轻化工大学 | A kind of bus-bar fault area determination method based on power frequency quantity Polarity comparision |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111948493A (en) * | 2020-08-21 | 2020-11-17 | 兰州理工大学 | MMC-HVDC direct current transmission line fault positioning method |
| CN113030648A (en) * | 2021-03-19 | 2021-06-25 | 国网河北省电力有限公司邢台供电分公司 | Method and device for determining position of fault point of power cable and terminal equipment |
| CN115267417A (en) * | 2022-06-15 | 2022-11-01 | 北京妙微科技有限公司 | Accurate positioning method for transmission line fault and transmission line traveling wave measuring device |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bahmanyar et al. | A comparison framework for distribution system outage and fault location methods | |
| Pourahmadi-Nakhli et al. | Path characteristic frequency-based fault locating in radial distribution systems using wavelets and neural networks | |
| CN110133445A (en) | A kind of submarine cable fault distance-finding method, terminal device and storage medium | |
| Goudarzi et al. | Improved fault location algorithm for radial distribution systems with discrete and continuous wavelet analysis | |
| CN103293449B (en) | Method for removing single-terminal traveling wave fault location dead area of high-voltage power grid in coal mine | |
| Deng et al. | Travelling-wave-based fault location algorithm for hybrid transmission lines using three-dimensional absolute grey incidence degree | |
| Shu et al. | Single-ended fault location for hybrid feeders based on characteristic distribution of traveling wave along a line | |
| CN105004969B (en) | The localization method of overhead line one-phase earthing failure in electric distribution network | |
| Naik et al. | Adaptive differential relay coordination for PV DC microgrid using a new kernel based time-frequency transform | |
| Xie et al. | A novel fault location method for hybrid lines based on traveling wave | |
| CN105353269B (en) | A kind of online fault distance-finding method of high-tension cable | |
| Ngaopitakkul et al. | Combination of discrete wavelet transform and probabilistic neural network algorithm for detecting fault location on transmission system | |
| Nan et al. | Asynchronous fault location scheme for half‐wavelength transmission lines based on propagation characteristics of voltage travelling waves | |
| CN108562820A (en) | Single-phase ground fault line selecting method of small-electric current grounding system and device | |
| CN106154117A (en) | A kind of combined method of distributed grid-connected isolated island detection | |
| Khalili et al. | Traveling wave-based protection for SVC connected transmission lines using game theory | |
| Chen et al. | A novel method for SLG fault location in power distribution system using time lag of travelling wave components | |
| Zhu et al. | Multiple random forests based intelligent location of single-phase grounding fault in power lines of DFIG-based wind farm | |
| Gaur et al. | New ground fault location method for three-terminal transmission line using unsynchronized current measurements | |
| Wu et al. | A new principle for initial traveling wave active power differential busbar protection | |
| Ngaopitakkul et al. | Identification of fault locations in underground distribution system using discrete wavelet transform | |
| Abasi et al. | Location of double-circuit grounded cross-country faults in GUPFC-compensated transmission lines based on current and voltage phasors analysis | |
| Lin et al. | Fault location for three-ended ring-topology power system using minimum GPS-based measurements and CVT/CT sensing | |
| Hamidi et al. | A learning-based framework for locating faults on power grid lines based on distortion of traveling waves | |
| CN112782532A (en) | Power distribution network fault location method based on traveling wave signal generated by circuit breaker closing |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190816 |