CN108693444A - A kind of failure point of power transmission line localization method based on electromagnetic signal time reversal - Google Patents
A kind of failure point of power transmission line localization method based on electromagnetic signal time reversal Download PDFInfo
- Publication number
- CN108693444A CN108693444A CN201810498192.1A CN201810498192A CN108693444A CN 108693444 A CN108693444 A CN 108693444A CN 201810498192 A CN201810498192 A CN 201810498192A CN 108693444 A CN108693444 A CN 108693444A
- Authority
- CN
- China
- Prior art keywords
- transmission line
- real
- matrix
- norm
- fault
- 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/085—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution lines, e.g. overhead
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/088—Aspects of digital computing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/11—Locating faults in cables, transmission lines, or networks using pulse reflection methods
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Locating Faults (AREA)
Abstract
The present invention is a kind of failure point of power transmission line localization method based on electromagnetic signal time reversal, this method is based on electromagnetism time reversal principle, single observation point is arranged in any one terminal in distribution network, measures and records one section of fault high frequency transient signal having in limit;Measuring signal is subjected to time backward, is injected into distribution network Simulation Calculation by observation point, the norm value at different conjecture positions in power distribution network is calculated, by finding peak value maximum and energy root maximum come failure judgement position.This method is different from traditional judgement traveling wave method and impedance method, the distribution network of a complex topology structure can be covered by an observation point, it is not influenced by fault impedance size, while there is stronger anti-noise ability, there is stronger practical value in transmission open acess.
Description
Technical field
The invention belongs to relay protection of power system and distribution network failure positioning field, and failure is based on more particularly to one kind
High frequency transient signal analyzes fault diagnosis and localization method.
Background technology
Transmission line of electricity is the element of power grid, and due to having a very wide distribution, and ultra-high-tension power transmission line is exposed to for a long time
In the severe natural environment such as wind, rain, thunder, electricity, inevitably break down.These failures directly affect transmission line of electricity peace
Full reliability service even can cause large area blackout when serious.Therefore, failure should investigate abort situation as early as possible after occurring,
Process accident, fast recovery of power supply, to improving, safe operation of power system is horizontal, operational reliability is particularly significant for this.Circuit is sent out
Trouble-shooting point is a key technology for ensureing power network safety operation after raw failure, and the positioning failure of precise and high efficiency can
Mitigate line walking burden, save a large amount of manpower and materials, and restoring electricity for circuit can be accelerated, reduces the comprehensive warp caused by power failure
Ji loss, therefore have important practical significance to transmission line of electricity accurate fault location.
Traditional transmission open acess technology is broadly divided into impedance method and traveling wave method two types at present, the former algorithm
Simply, it is easy to accomplish, disadvantage is that positioning result is affected by fault resstance, can not be suitable for multiple terminals Distributed power
Net.The main stream approach that traveling wave method is positioned as current power transmission line fault needs to estimate the spread speed of traveling wave in the line, profit
Abort situation is analyzed and estimated to the transient state travelling wave signal generated with transmission line malfunction.Traveling wave method in recent years at home and abroad from
It is theoretical to be obtained for faster development to application, but there is also some still unsolved technical barriers, such as high resistance earthing fault
When transient state travelling wave weak output signal lead to not extract pulse characteristics, the arrival time of main bang and catadioptric signal is difficult to know
Not etc..Localization method can be divided into single-ended mensuration, double-end measurement method and multiple-end measurement method again from observation point quantity.For more
Terminal distribution power network, there are multi-load, the complex topology structure of more bifurcateds, distributed power generation, single-ended mensuration is difficult to
It is applicable in, double-end measurement method (or multiple-end measurement method) needs to carry out high-precision time synchronization and data to the measurement data of each terminal
Telecommunication, this undoubtedly increases the complexity and cost of monitoring system.The above method improves traveling wave method to a certain extent
Reliability, but the validity of fault location still depends critically upon the extraction effect of wavefront.
Invention content
The purpose of this invention is to provide a kind of electrical power distribution network fault location method based on single-ended measurement, is not necessarily to time synchronization
Communication apparatus is adaptable to more bifurcateds, the fault diagnosis of multiple terminals distribution power network and quickly positioning.
Technical scheme is as follows:
A kind of failure point of power transmission line localization method based on electromagnetic signal time reversal, includes the following steps:
[1]According to the circuit types parameter of transmission line of electricity, wire topologies, transmission line conductive parameters, it is calculated defeated
The characterisitic parameter of electric line;
[2]Several failures conjecture position x is set in the model of power transmission system of foundationG;
[3]By the electric current, voltage transformer of relay protection system on transmission line of electricity, monitoring obtains faulty transmission circuit
The time domain Full wave shape of voltage u (t), electric current i (t), t ∈ [tf,tf+T], wherein tfIt is to measure initial time, when T is signal sampling
Between length, and the waveform in finite time length is acquired record;
[4]To step[3]The acquired electromagnetic data of acquisition carries out time backward inverting, obtains the time in acquisition duration T
Backward signal;
ur(t)=u (T-t), ir(t)=i (T-t);
[5]Based on step[1]The characterisitic parameter of transmission line of electricity establish transmission line simulation model, backward signal is input to
In transmission line simulation model, fault observation sites position x is calculatedGThe current or voltage at place responds, and according to norm theory meter
The peak value norm responded in duration 2T is calculated, the fault observation sites position corresponding to the maximum of peak value norm is obtained
xp,realFault observation sites position x corresponding to maximum with energy norme,real, and judge to obtain the position where fault point
It sets;
The mathematic(al) representation of wherein peak value Number Norm Criterion is:
The mathematic(al) representation of wherein energy norm criterion is:
X in formulaGTo guess the position of fault point;I(xG, t) and it is conjecture position xGT changes the short circuit current at place at any time;
Sup is xGLocate the peak value of short circuit current absolute value;It indicates when norm takes maximum, corresponding independent variable xGValue.
[6]Fault point judges:
[6.1]According to step[2]In x is calculated when distance is d between adjacent fault pointp,realWith xe,real, realize just
Walk fault location:
[6.2]Work as xp,realWith xe,realAs a result in the case of consistent or close, in xp,realOr xe,realNear zone takes step
[2]In the distance between adjacent fault point be d/5~d/20, repeat step[1]~[5], obtain the exact position of fault point
Realize that fault point is accurately positioned;
[6.3]Work as xp,realWith xe,realWhen result difference is larger, in xe,realNear zone takes step[2]In adjacent failure
The distance between point is d/5~d/20, repeats step[1]~[5], it is accurately positioned, obtains the exact position of fault point
Realize that fault point is accurately positioned.
Further, step[1]The characterisitic parameter calculation formula of middle transmission line of electricity is:
In formula:Z(ω)'Indicate unit length transmission line impedance matrix;Y(ω)'Indicate unit length transmission-line admittance square
Battle array;L'Indicate unit length transmission line inductance matrix;Z'wIndicate unit length transmission line resistor matrix;Z'gIndicate unit length
Transmission line the earth impedance matrix;C'Indicate unit length transmission line capacitor matrix;Y'gWith indicating unit length transmission line area mobility
Matrix.
Further, step[5]The analytic equation of middle transmission line simulation model is:
Wherein, V, I are node voltage, current phasor, Y respectivelyCIt is transmission line character admittance matrix, S is cable network
Scattering parameter hypermatrix, T are pipeline incidence matrix, and Γ is to propagate hypermatrix, W(s)It is synthesis driving source super vector , [1]Indicate with
The consistent unit matrix of S and Γ matrix dimensionalities.
Further, step[6.2]Middle xp,realNear zone refers to x on transmission line of electricityp,real- d~xp,realThe position of+d
Place;xe,realNear zone refers to x on transmission line of electricitye,real- d~xe,realAt the position of+d.
Compared with prior art, the advantageous effects that the present invention has are as follows:
1. the present invention is based on the data of single observation station, without carrying out multiple observation station data communications and time synchronization, contract
It while small hardware cost, can eliminate due to the introduced position error of time synchronization error, while in application by failure
The influence of noise introduced when impedance and measurement is smaller, improves positional accuracy.
2. present invention utilizes the adapted local cosine transform principle of time reversal technology, it is temporary for failure that it is different from conventional method
The discriminatory analysis of state signal traveling wave details, is integrally injected, Inversion Calculation is sentenced by norm based on the sampled signal in finite time
According to realizing fault location.Avoid the erroneous judgement of traveling wave arrival time and wave head feature the problems such as being difficult to extract.In multiple terminals point
Cloth distribution network, make somebody a mere figurehead-application backgrounds such as buried mixing distribution line (inhomogeneous medium) under with unique advantage.
3. the present invention realizes the accurate positionin to transmission line malfunction using energy norm and the double criterions of peak value norm, fully
Consider the noise immunity of the stability of root criterion and peak value Number Norm Criterion in energy norm.Using first in abort situation calculating
, length resolution lower coarse positioning larger into row span d obtains the approximate region of fault point, then smaller into row span d, long
The higher fine positioning of resolution ratio is spent, has greatly saved computing cost in this way, accelerates Inversion Calculation speed, and measurement result meets
Long transmission line, high-precision positioning requirements.The mean value of energy norm locator value and peak value norm locator value is used in calculating so that
Positioning result takes into account stronger anti-noise ability, focusing and stability requirement simultaneously.
Description of the drawings
Fig. 1 is technical solution of the present invention flow chart;
Fig. 2 is test club tower dimensional drawing of the present invention;
Fig. 3 experimental configuration schematic diagrames of the present invention;
Fig. 4 is the faulted phase current waveform for not considering cross talk conditions;
Fig. 5 is the source-series injection transmission line model schematic diagram of inverting voltage;
Fig. 6 is that transmission line model schematic diagram is injected in inverting current source parallel connection;
Fig. 7 is not consider the norm value of cross talk conditions along circuit distribution map;
Fig. 8 is the three-phase current waveform for considering cross talk conditions;
Fig. 9 is to consider that the norm value of cross talk conditions is distributed along circuit.
Specific implementation mode
In recent years, time reversal (Time Reversal, TR) is as a kind of new skill based on signal time domain backout
Art has the characteristics that time and space double focus, in acoustics, biomedicine, radar imagery, communication, electromagnetic radiation positioning
Etc. multiple fields be rapidly developed.Consider under lossless case, describing two-conductor line, (or ideal big ground is singly led
Line circuit), telegraph equation that we can be following describes the communication process of wave
Wherein v (x, t) and i (x, t) is the voltage and current on transmission line of electricity, L&apos respectively;And C'It is transmission line unit length
Inductance and capacitance.
To above-mentioned equation application time inversion (t →-t), can obtain
Wherein, since the direction of motion of charge changes, while t takes negative sign, reply electric current i (t) takes negative.It can be with
Find out, in lossless transmission line of electricity, telegraph equation meets time reversal invariance.In the case where the time taking inverse, electromagnetic field
Variation and communication process are referred to as electromagnetism time reversal (Electromagnetic Time Reversal, EMTR), and EMTR is from reason
It can be used for the positioning of transmission line fault localization by above having confirmed, since EMTR has space-time focusing in non-uniform dielectric,
The technology has under multiple terminals distribution power network, the application backgrounds such as aerial-buried mixing distribution line (inhomogeneous medium)
There is unique advantage.It is different from discriminatory analysis of the conventional method for fault transient signals traveling wave details, avoids traveling wave arrival
The problems such as moment erroneous judgement and wave head feature are difficult to extract.
Norm is commonly used to the feature of quantitative description electromagnetic pulse waveform, such as peak value, energy, rectification integrated value, norm carry
A kind of unified mathematical form is supplied, to state the sensitive parameter of electromagnetic pulse.The present invention is based on electromagnetism time reversal technologies to carry
The novel fault localization method for fault high frequency transient signal gone out has been invented in conjunction with norm theory for being based on 2- norms
Energy criterion and peak value criterion based on ∞-norm, can be realized to Distributed power net in the case of single observation point
It is accurately positioned.In the method, we are only concerned about the 2- norms and ∞-norm of transient signal, expression formula and corresponding physics
Meaning is as shown in table 1.
The mathematic(al) representation and its physical significance of 1 p- norms of table
The mechanism that failure point of power transmission line positioning is carried out based on norm theory is as follows:
One, the peak value criterion based on ∞-norm
By the radiation field propagated in strong scattering, reflecting medium due to multiple reflection and scattering effect, when passing through
Between after inverting, electromagnetic wave can be converged in a certain particular moment at the position of radiation source, i.e., when electromagnetism time reversal process has
Empty focusing, and the electromagnetic field converged can realize the maximum of amplitude at source position.This property is equally applicable to transmit
Lineation opinion, the communication process of electromagnetic wave on the transmission line are equivalent to external electromagnetic field and are propagated in waveguiding structure in a tem mode, can
The propagation property of electromagnetic field three dimensions is applied in one-dimensional scene with regarding as.When short trouble has occurred in transmission line, therefore
The high frequency transient pulse signal that barrier generates can be by a step voltage source V over the groundFTo indicate.Consider that is maked somebody a mere figurehead a lossless biography
Defeated line, reference ground are perfact conductor, and shorted to earth failure, the transient state electromagnetic pulse that fault overvoltage generates occur at position
It is propagated to both ends along cable, electromagnetic wave encounters at short fault location or cable terminations can reflect, transmission line head end electric current
Time domain analytic solutions are:
Here, Z1Indicate the load of cable head end, ρ1Indicate the reflectance factor of cable head end, TdtIndicate fault traveling wave
Travel to the time delay of cable head end.Transient current is made of the fault-current signal and subsequent reflection signal reached for the first time.
When transient signal is carried out time backward by us, and refill in telegram in reply road, then the short circuit electricity of short-circuit point
Flowing expression formula is:
Here TrtIndicate that traveling wave is propagated to the time delay for guessing abort situation by cable head end in refutation process.When conjecture position
When setting equal to true fault position, Trt=Tdt, head reaches signal and follow-up signal forms in-phase stacking, the expression formula of short circuit current
It can be written as:
It follows that when the position of conjecture is consistent with abort situation, short circuit current signal that Inversion Calculation obtains due to
The in-phase stacking of pulse forms the great effect of amplitude, this is that the space-time focussing property of time reversal technology is determined.
We can be judged using this property by finding the maximum of peak value of short at different conjecture positions
Abort situation, the peak value criterion (∞-norm) for obtaining the fault location based on time reversal technology are:
In formula:xG--- the position of conjecture fault point;I(xG, t) --- conjecture position xGT becomes the short circuit current at place at any time
Change;Sup --- indicate xGLocate the peak value of short circuit current absolute value;--- it indicates when norm takes maximum, it is corresponding from change
Measure xGValue.
Two, the energy criterion based on 2- norms
For Lossless transport line model, it can be regarded as Two-port netwerk model, we use h (t) and hr(t) failure is indicated
With the system shock response during time reversal.Since Lossless transport line model meets reciprocity, i.e.,:
H (t)=hr(t)
Assuming that from the two-port network input signal e (t) that one end constitutes transmission line, output signal isIts
InIndicate convolution.Output function is done into time reversal operation, and transmission back gauze network is refilled from output end, in input terminal
Obtained response be:
Here, Rhh(t) cross-correlation function is indicated.The frequency-domain expression of g (t) can be written as:
G (f)=E*(f)·H*(f) H (f)=E*(f)·|H(f)|2
According to Holder inequality,
And if only if X (f)=KY*(f) when, equal sign is set up.
According to Pa Saiwaer (Parseval) theorem,
When guessing that abort situation is equal with physical fault position, Hr(f)=H (f), equal sign are set up.That is, when guessing
Location set place transmission function it is equal with the transmission function of physical fault position when, the output signal energy after time reversal
Reach maximum.We can come by comparing the 2- norms (energy root) of the short circuit current at different conjecture positions as a result,
Trouble-shooting position.Energy criterion (2- Number Norm Criterions) is written as:
X in formulaG--- the position of conjecture fault point;I(xG, t) --- conjecture position xGThe short circuit current at place t at any time
Variation;E(xG, t) --- conjecture position xGT changes the short circuit current energy at place at any time;--- it indicates to take most when norm
When big, corresponding independent variable xGValue.
The present invention is based on the failure point of power transmission line localization methods of electromagnetic signal time reversal, specifically comprise the following steps:
[1]According to the circuit types parameter of transmission line of electricity, wire topologies, transmission line conductive parameters, it is calculated defeated
The characterisitic parameter of electric line;Wherein circuit types parameter includes overhead transmission line or buried cable, phase line quantity, cable dimensions, leads
Body overhead height, conductor spacing etc.;Wire topologies include wiring path, prong numbers, bifurcation site, leg length;It passes
Defeated line conductive parameters include the electric parameter of conductor and the earth, are used for the foundation of simulation model.
It says in general sense, the characterisitic parameter calculation formula of multi-conductor transmission lines is:
In formula:Z(ω)'Indicate unit length transmission line impedance matrix;Y(ω)'Indicate unit length transmission-line admittance square
Battle array;L'--- unit length transmission line inductance matrix;Z'w--- unit length transmission line conductor impedance matrix;Z'g--- unit
Length transmission line the earth impedance matrix;C'--- unit length transmission line capacitor matrix;Y'g--- unit length transmission line ground
Admittance matrix.
For multi-conductor transmission lines, the inductance matrix of unit length includes the mutual inductance between the self-induction of conducting wire and conducting wire, is calculated
Formula is:
In formula:μ0--- the magnetic capacity in free space;hi--- i-th transmission line height;ri--- i-th transmission
The radius of line;hj--- jth root transmission line height;dij--- the spacing of i-th transmission line and jth root transmission line.Lii'It is inductance
The diagonal entry of matrix represents the self-induction of conducting wire;Lij'It is the off diagonal element of inductance matrix, represents the mutual inductance between conducting wire,
They together constitute inductance matrix.Impedance matrix is similar with its.
The capacitance matrix of unit length can be calculated from inductance matrix:
C'=μ0ε0L'-1
Unit length transmission line resistor matrix generally only considers that conductor DC resistance, calculation formula are:
In formula, σw--- it is transmission line conductors conductivity.
Consider that the earth is lost, the expression formula of the earth impedance is:
In formula, γgFor the earth propagation constant, expression formula is:
In formula, σg--- the conductivity of the earth, εg--- the earth relative dielectric constant.
The expression formula of the earth admittance is:
Yg'=γg 2Zg'-1
Transmission line characteristic impedance expression formula is:
Transmission line character admittance expression formula is:
Yc=Zc'-1
[2]Several failures conjecture position x is set in the model of power transmission system of foundationG, adjacent failure conjecture position it
Between distance be d;
[3]By electric current, the voltage transformer in relay protection system, voltage u (t), the electricity on faulty transmission road are monitored
Flow the time domain Full wave shape of i (t), t ∈ [tf,tf+T], wherein tfIt is to measure initial time, T is signal sampling time span, and will
Waveform in finite time length is acquired record;
[4]To step[3]The acquired electromagnetic data of acquisition carries out time backward inverting, obtains the backward in acquisition duration T
Signal ur(t)=u (T-t), ir(t)=i (T-t);
[5]Based on step[1]The characterisitic parameter of transmission line of electricity establish transmission line simulation model, inverting signal is input to
In transmission line simulation model, failure conjecture position x is calculatedGThe current or voltage at place responds, and according to norm theoretical calculation
The energy norm and peak value norm responded in duration 2T.
The mode of transmission line modeling can be time domain modeling and Modeling In Frequency Domain.Time domain modeling mainly utilizes EMTP-RV etc.
Business software calls software model to generate circuit model, the time-domain-simulation embedded by software by inputting transmission line parameter
Algorithm completes simulation calculation.Modeling In Frequency Domain is first established single based on BLT frequency domains equation method (Baum-Liu-Tesche equation)
Bit length equivalence linear electrical parameter calculates under frequency domain voltage, current-responsive on transmission line by the form of analytic equation.
The expression formula of wherein energy norm is:
The expression formula of peak value norm is:
[5.1]Building of Simulation Model
According to BLT Equation Theories, transmission circuit network can be abstractively with a series of pipeline (tube), node
(junction) etc. physical quantitys characterize, and the voltage, electric current on each node can write out its solution by BLT frequency domains equation
Analyse expression formula:
Wherein, V, I are node voltage, current phasor, Y respectivelyCIt is transmission line character admittance matrix, S is cable network
Scattering parameter hypermatrix, T are pipeline incidence matrix, and Γ is to propagate hypermatrix, W(s)It is synthesis driving source super vector.Wherein transmit
Line Characteristic mobility matrix YCAccording to step[1]In unit length linear electrical parameter calculate and obtain.
[5.2]Calculate scattering parameter hypermatrix
The connection of node both ends is all pipeline, cannot be indicated with reflectance factor matrix, will scatter determining for hypermatrix at this time
Adopted formula is combined the scattering hypermatrix for solving desired node with the Kirchoff s voltage current law of node.According to voltage and current and conjunction
At the relational expression of voltage wave, in conjunction with Kirchoff s voltage current law, the reflection of wave is converted into voltage-current relationship, you can
To the scattering hypermatrix expression formula of desired node:
CV, CIFor the coefficient matrix of the Kirchoff s voltage current law of node, YCIt is characteristic admittance matrix, ZCIt is characteristic
Impedance matrix.
[5.3]Calculate incidence matrix T
Wherein, Wu, WvIndicate the wave vector propagated in transmission network.
[5.4]It calculates and propagates hypermatrix Γ
Internet communication hypermatrix is made of propogator matrix of the wave in each pipeline, and the propogator matrix of each pipeline is by cable
Mode propagation constant γ and duct length l determine.It propagates hypermatrix and characterizes phase-amplitude during the wave in whole network is propagated
Situation of change, expression formula are:
[5.5]Calculate driving source super vector
As lump driving source, expression formula is the transient state pulse signal that transmission line failure generates:
The Internet communication hypermatrix being calculated, network scattering hypermatrix, driving source super vector are substituted into BLT hypermatrix
Equation acquires each resultant voltage wave.According to the relationship of resultant voltage wave and each node voltage current-responsive, parsed according to BLT frequency domains
Equation solution obtains the voltage and current response of each node of series transmission lines.
[5.6]The mathematic(al) representation that energy norm criterion is calculated is:
The expression formula of peak value criterion is:
[6]Fault point judges:
[6.1]According to step[2]In x is calculated when distance is d between adjacent fault pointp,realWith xe,real, realize just
Walk fault location:
[6.2]Work as xp,realWith xe,realAs a result in the case of consistent or close, in xp,realOr xe,realNear zone takes step
[2]In the distance between adjacent fault point be d/5~d/20, repeat step[1]~[5], obtain the exact position of fault point
Wherein xp,realNear zone refers to x on transmission line of electricityp,real- d~xp,realAt the position of+d;xe,realNear zone
Refer to x on transmission line of electricitye,real- d~xe,realAt the position of+d.
[6.3]Work as xp,realWith xe,realWhen result difference is larger, in xe,realNear zone takes step[2]In adjacent failure
The distance between point is d/5~d/20, usually selects d/10, repeats step[1]~[5], it is accurately positioned, obtains fault point
Exact position
This allows for the criterion that energy norm uses root, thus has stronger stability, therefore with energy norm
Fine positioning is carried out subject to Primary Location result and in its vicinity, to realize that fault point is accurately positioned.
With reference to field experiment practical measuring examples, the present invention is described in further detail.
(1) experimental configuration
This experiment is based on the double loops 10kV power distribution network, and overhead line parameter and shaft tower scene pictorial diagram are as shown in Figure 2,3, passes
Defeated line design parameter is as shown in table 2.
2 experimental line major parameter of table
Experimental line overall length 677m, totally 11 shaft towers, number are 22 to 32.Using the double exponential type high-tension pulses of a nanosecond
It rushes source and a pulse signal is injected for simulated failure overvoltage pulse signal, pulse source output terminal to circuit at No. 23 shaft towers
Failure phase is connect, ground terminal is connected with shaft tower the earth.The clock is only used for simulating generating short trouble overvoltage in this experiment,
The not necessarily condition in the application of practical this method.Observation point is arranged on No. 22 shaft towers of circuit head end for we, end load
With 1500pF thin film capacitors come analogue transformer busbar entrance capacitance, this is chosen according to the characterisitic parameter of real transformer
Alternative model, transmission line of electricity is directly accessed in transformer during physical fault judges.Experimental configuration schematic diagram is as shown in Figure 3.With
Electric current loop measures the fault current on conductor, and measuring apparatus and its major parameter are as shown in table 3.
3 measuring apparatus of table and its major parameter
(2) fault transient waveform is measured
We consider two kinds of experiment measurement schemes, are respectively free of the single-phase measurement of crosstalk factor and include crosstalk factor
Three-phase measures.
A) the single-phase measurement and its interpretation of result of crosstalk factor are not considered
For single-phase measurement, capacitive load only is connected in failure phase end, other non-faulting are mutually open circuit, in circuit head end
The electric current time domain waveform that place measures is as shown in Figure 4.It is larger in noise, denoising is carried out to waveform, utilizes formula
(1) time domain backward is carried out.
We are swashed the measuring signal after inverting with the parsing equation method modeling of BLT frequency domains at line port as lump
The source of encouraging is injected into circuit model, as shown in Figure 5.If measuring signal is voltage signal, according to Thevenin's theorem, desired voltage
Source should connect in access transmission line circuit, wherein V1(T-t) voltage signal after backward, Z are indicated1, Z2Indicate the negative of first and last end
It carries, IF(t) short circuit current at conjecture position is indicated.
The electric current after backward as ideal current source and is coupled according to Nortons theorem if measuring signal is current signal
Enter, then it is in parallel with fault impedance to inject source, as shown in Figure 6.Wherein, I1(T-t) current signal after backward, Z are indicated1, Z2It indicates
The load at first and last end, IF(t) short circuit current at conjecture position is indicated.
And the short circuit current norm value at different conjecture positions is calculated, distribution results are as shown in Figure 7 along the line.From result figure
In as can be seen that the energy and peak value of short circuit current at No. 23 shaft towers are maximum.
B) consider the three-phase measurement and its interpretation of result of crosstalk factor
The three-phase considered under cross talk conditions is measured, capacitive load is all connected at the first and last end of three-phase line, in circuit
Head end measures the current signal in triple line simultaneously, and time domain waveform is as shown in Figure 8.By Inversion Calculation, obtained failure phase
On norm along be distributed it is as shown in Figure 9.It can be seen that energy and peak value maximum still match with abort situation.By
It is far smaller than norm value in failure phase in the short circuit current norm value in non-faulting phase, only depicts the norm of failure phase here
Distribution.It follows that by the statistics to norm value, we can also be by this method come failure judgement phase.
Studies have shown that ∞-Number Norm Criterion based on peak value has for the stronger anti-noise of the noise signal introduced in measurement
Ability, and the amplitude at each moment of inverting signal is carried out due to introducing integral and calculating based on the 2- norms of energy root
Cumulative summation, inverting output result more they tend to stablize.A large amount of verification experimental verification shows that both methods judges to obtain the one of result
Cause property is fine.Two kinds of criterions should be used in combination in practical applications to estimate abort situation, taken under higher distance resolution
The arithmetic mean of instantaneous value of the two positioning result may make result both to have stronger anti-noise ability, while also have stronger stabilization
Property.
Claims (5)
1. a kind of failure point of power transmission line localization method based on electromagnetic signal time reversal, which is characterized in that including walking as follows
Suddenly:
[1]According to the circuit types parameter of transmission line of electricity, wire topologies, transmission line conductive parameters, power transmission line is calculated
The characterisitic parameter on road;
[2]Several failures conjecture position x is set in the model of power transmission system of foundationG, between adjacent failure conjecture position away from
From for d;
[3]By the electric current, voltage transformer of relay protection system on transmission line of electricity, monitoring obtains the voltage of faulty transmission circuit
The time domain Full wave shape of u (t), electric current i (t), t ∈ [tf,tf+T], wherein tfIt is to measure initial time, T is that the signal sampling time is long
Degree, and the waveform in finite time length is acquired record;
[4]To step[3]The acquired electromagnetic data of acquisition carries out time backward inverting, obtains the time backward in acquisition duration T
Signal data;
ur(t)=u (T-t), ir(t)=i (T-t)
[5]Based on step[1]The characterisitic parameter of transmission line of electricity establish transmission line simulation model, inverting signal is input to transmission
In line simulation model, failure conjecture position x is calculatedGThe current or voltage at place responds, and is being held according to norm theoretical calculation
The peak value norm responded in continuous time 2T, obtains the fault observation sites position x corresponding to the maximum of peak value normp,realAnd energy
Measure the fault observation sites position x corresponding to the maximum of norme,real, and judge to obtain the position where fault point;
The mathematic(al) representation of wherein peak value Number Norm Criterion is:
The mathematic(al) representation of wherein energy norm criterion is:
X in formulaGTo guess the position of fault point;I(xG, t) and it is conjecture position xGThe short circuit current at place;Sup indicates xGPlace's short circuit
The peak value of current absolute value;It indicates when norm takes maximum, corresponding independent variable xGValue.
[6]Fault point judges:
[6.1]According to step[2]In x is calculated when distance is d between adjacent fault pointp,realWith xe,real, realize preliminary event
Barrier positioning;
[6.2]Work as xp,realWith xe,realAs a result in the case of consistent or close, in xp,realOr xe,realNear zone takes step[2]
In the distance between adjacent fault point be d/5~d/20, repeat step[1]~[5], obtain the exact position of fault point
Realize that fault point is accurately positioned;
[6.3]Work as xp,realWith xe,realWhen result difference is larger, in xe,realNear zone takes step[2]In adjacent fault point it
Between distance be d/5~d/20, repeat step[1]~[5], it is accurately positioned, obtains the exact position of fault point
Realize that fault point is accurately positioned.
2. the failure point of power transmission line localization method according to claim 1 based on electromagnetic signal time reversal, step[1]
The characterisitic parameter calculation formula of middle transmission line of electricity is:
In formula:Z(ω)'Indicate unit length transmission line impedance matrix;Y(ω)'Indicate unit length transmission-line admittance matrix;L'
Indicate unit length transmission line inductance matrix;Z'wIndicate unit length transmission line resistor matrix;Z'gIndicate unit length transmission line
The earth impedance matrix;C'Indicate unit length transmission line capacitor matrix;Y'gIndicate unit length transmission line ground admittance matrix.
3. the failure point of power transmission line localization method according to claim 1 based on electromagnetic signal time reversal, feature
It is:Step[5]The analytic equation of middle transmission line simulation model is:
Wherein, V, I are node voltage, current phasor, Y respectivelyCIt is transmission line character admittance matrix, S is the scattering ginseng of cable network
Number hypermatrix, T are pipeline incidence matrix, and Γ is to propagate hypermatrix, W(s)It is synthesis driving source super vector , [1]It indicates and S and Γ
The consistent unit matrix of matrix dimensionality.
4. the failure point of power transmission line localization method according to claim 1 based on electromagnetic signal time reversal, feature
It is:Step[6.2]Middle xp,realNear zone refers to x on transmission line of electricityp,real- d~xp,realAt the position of+d;xe,realNear
Region refers to x on transmission line of electricitye,real- d~xe,realAt the position of+d.
5. the failure point of power transmission line localization method according to claim 1 based on electromagnetic signal time reversal, feature
It is:Step[6.3]Middle xe,realNear zone refers to x on transmission line of electricitye,real- d~xe,realAt the position of+d.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810498192.1A CN108693444B (en) | 2018-05-23 | 2018-05-23 | A kind of failure point of power transmission line localization method based on electromagnetic signal time reversal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810498192.1A CN108693444B (en) | 2018-05-23 | 2018-05-23 | A kind of failure point of power transmission line localization method based on electromagnetic signal time reversal |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108693444A true CN108693444A (en) | 2018-10-23 |
| CN108693444B CN108693444B (en) | 2019-08-06 |
Family
ID=63847615
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810498192.1A Active CN108693444B (en) | 2018-05-23 | 2018-05-23 | A kind of failure point of power transmission line localization method based on electromagnetic signal time reversal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108693444B (en) |
Cited By (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110531216A (en) * | 2019-07-15 | 2019-12-03 | 重庆大学 | A kind of overhead line and cable mixed power transmission line fault section method of discrimination |
| CN111157845A (en) * | 2020-01-07 | 2020-05-15 | 国网辽宁省电力有限公司朝阳供电公司 | Time domain inversion technology-based fault positioning method suitable for power distribution network |
| CN111313932A (en) * | 2020-02-07 | 2020-06-19 | 西安交通大学 | Transmission line based on electromagnetic time reversal transfer function correlation and network fault positioning method |
| CN111458597A (en) * | 2020-02-06 | 2020-07-28 | 云南电网有限责任公司电力科学研究院 | Fault positioning method based on lightning overvoltage phase transmission characteristic |
| CN111596165A (en) * | 2020-04-13 | 2020-08-28 | 上海紫通信息科技有限公司 | Power distribution network fault simulation method based on steady-state transient smooth splicing output |
| EP3726681A1 (en) * | 2019-04-19 | 2020-10-21 | Supergrid Institute | Transient based method for identifying faults in a high / medium voltage electric power transmission system, fault identification module and power transmission system |
| CN112083282A (en) * | 2020-09-04 | 2020-12-15 | 昆明理工大学 | New energy station outgoing line fault location method based on electromagnetic time reversal method |
| CN112444887A (en) * | 2019-09-03 | 2021-03-05 | 中国科学院物理研究所 | Electromagnetic detection information inversion device and method |
| CN112906262A (en) * | 2021-01-27 | 2021-06-04 | 哈尔滨工程大学 | Method for predicting overhead power line crosstalk caused by external transient electromagnetic interference |
| CN112946369A (en) * | 2021-01-25 | 2021-06-11 | 西安理工大学 | Cable aging identification method using power line communication equipment |
| CN113030244A (en) * | 2021-03-22 | 2021-06-25 | 国网江苏省电力有限公司无锡供电分公司 | Power transmission line tower corrosion defect magnetic leakage detection signal inversion imaging method and system |
| CN113030651A (en) * | 2021-03-31 | 2021-06-25 | 国网上海市电力公司 | Method for positioning impact surge current |
| CN113219298A (en) * | 2021-03-24 | 2021-08-06 | 昆明理工大学 | Complex alternating current power grid fault current traveling wave numerical simulation method |
| CN113311294A (en) * | 2021-07-29 | 2021-08-27 | 清华大学 | Branch line short circuit fault positioning method based on time domain inversion |
| CN113311287A (en) * | 2021-04-28 | 2021-08-27 | 上海电力大学 | Fault location method applied to direct-current power distribution network controlled by two-end VSC |
| CN113466612A (en) * | 2021-06-02 | 2021-10-01 | 湖南省计量检测研究院 | Cable fault identification method |
| CN113625200A (en) * | 2021-08-16 | 2021-11-09 | 哈尔滨工业大学 | In-situ detection method and detection device for performance degradation of coaxial connector |
| CN114325244A (en) * | 2022-01-13 | 2022-04-12 | 清华大学 | Line short-circuit fault positioning method and system based on direct signal convolution |
| CN114689997A (en) * | 2022-05-06 | 2022-07-01 | 国网江苏省电力有限公司宿迁供电分公司 | Distribution network cable fault identification and positioning method |
| CN114844776A (en) * | 2022-04-22 | 2022-08-02 | 电子科技大学 | Electromagnetic imaging positioning method and device based on time reversal time domain characteristics |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103499770A (en) * | 2013-09-29 | 2014-01-08 | 华北电力大学（保定） | Grid fault positioning method based on wide area voltage |
| 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 |
| CN107228988A (en) * | 2017-06-28 | 2017-10-03 | 西安交通大学 | Electromagnetic pulse multiport effective matrix method based on coupling current signals hypernorm |
-
2018
- 2018-05-23 CN CN201810498192.1A patent/CN108693444B/en active Active
Patent Citations (3)
| 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 |
| CN103499770A (en) * | 2013-09-29 | 2014-01-08 | 华北电力大学（保定） | Grid fault positioning method based on wide area voltage |
| CN107228988A (en) * | 2017-06-28 | 2017-10-03 | 西安交通大学 | Electromagnetic pulse multiport effective matrix method based on coupling current signals hypernorm |
Cited By (36)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3726681A1 (en) * | 2019-04-19 | 2020-10-21 | Supergrid Institute | Transient based method for identifying faults in a high / medium voltage electric power transmission system, fault identification module and power transmission system |
| WO2020212533A1 (en) * | 2019-04-19 | 2020-10-22 | Supergrid Institute | Transient based method for identifying faults in a high / medium voltage electric power transmission system, related fault identification module and high / medium voltage electric power transmission and / or distribution system equipped with such fault identification module |
| CN110531216B (en) * | 2019-07-15 | 2024-02-23 | 重庆大学 | Overhead line and cable hybrid transmission line fault section distinguishing method |
| CN110531216A (en) * | 2019-07-15 | 2019-12-03 | 重庆大学 | A kind of overhead line and cable mixed power transmission line fault section method of discrimination |
| CN112444887A (en) * | 2019-09-03 | 2021-03-05 | 中国科学院物理研究所 | Electromagnetic detection information inversion device and method |
| CN112444887B (en) * | 2019-09-03 | 2022-03-11 | 中国科学院物理研究所 | Electromagnetic detection information inversion device and method |
| CN111157845B (en) * | 2020-01-07 | 2021-10-26 | 国网辽宁省电力有限公司朝阳供电公司 | Time domain inversion technology-based fault positioning method suitable for power distribution network |
| CN111157845A (en) * | 2020-01-07 | 2020-05-15 | 国网辽宁省电力有限公司朝阳供电公司 | Time domain inversion technology-based fault positioning method suitable for power distribution network |
| CN111458597A (en) * | 2020-02-06 | 2020-07-28 | 云南电网有限责任公司电力科学研究院 | Fault positioning method based on lightning overvoltage phase transmission characteristic |
| CN111458597B (en) * | 2020-02-06 | 2022-06-10 | 云南电网有限责任公司电力科学研究院 | Fault positioning method based on lightning overvoltage phase transmission characteristic |
| CN111313932A (en) * | 2020-02-07 | 2020-06-19 | 西安交通大学 | Transmission line based on electromagnetic time reversal transfer function correlation and network fault positioning method |
| CN111313932B (en) * | 2020-02-07 | 2021-07-13 | 西安交通大学 | Transmission line based on electromagnetic time reversal transfer function correlation and network fault positioning method |
| CN111596165A (en) * | 2020-04-13 | 2020-08-28 | 上海紫通信息科技有限公司 | Power distribution network fault simulation method based on steady-state transient smooth splicing output |
| CN112083282A (en) * | 2020-09-04 | 2020-12-15 | 昆明理工大学 | New energy station outgoing line fault location method based on electromagnetic time reversal method |
| CN112083282B (en) * | 2020-09-04 | 2021-11-09 | 昆明理工大学 | New energy station outgoing line fault location method based on electromagnetic time reversal method |
| CN112946369A (en) * | 2021-01-25 | 2021-06-11 | 西安理工大学 | Cable aging identification method using power line communication equipment |
| CN112946369B (en) * | 2021-01-25 | 2024-04-26 | 西安理工大学 | Cable aging identification method using power line communication device |
| CN112906262A (en) * | 2021-01-27 | 2021-06-04 | 哈尔滨工程大学 | Method for predicting overhead power line crosstalk caused by external transient electromagnetic interference |
| CN113030244A (en) * | 2021-03-22 | 2021-06-25 | 国网江苏省电力有限公司无锡供电分公司 | Power transmission line tower corrosion defect magnetic leakage detection signal inversion imaging method and system |
| CN113030244B (en) * | 2021-03-22 | 2024-01-23 | 国网江苏省电力有限公司无锡供电分公司 | Inversion imaging method and system for transmission line tower corrosion defect magnetic flux leakage detection signal |
| CN113219298B (en) * | 2021-03-24 | 2022-10-11 | 昆明理工大学 | Fault current traveling wave numerical simulation method for complex alternating current power grid |
| CN113219298A (en) * | 2021-03-24 | 2021-08-06 | 昆明理工大学 | Complex alternating current power grid fault current traveling wave numerical simulation method |
| CN113030651B (en) * | 2021-03-31 | 2023-03-14 | 国网上海市电力公司 | Method for positioning impact surge current |
| CN113030651A (en) * | 2021-03-31 | 2021-06-25 | 国网上海市电力公司 | Method for positioning impact surge current |
| CN113311287A (en) * | 2021-04-28 | 2021-08-27 | 上海电力大学 | Fault location method applied to direct-current power distribution network controlled by two-end VSC |
| CN113466612B (en) * | 2021-06-02 | 2022-07-01 | 湖南省计量检测研究院 | Cable fault identification method |
| CN113466612A (en) * | 2021-06-02 | 2021-10-01 | 湖南省计量检测研究院 | Cable fault identification method |
| CN113311294A (en) * | 2021-07-29 | 2021-08-27 | 清华大学 | Branch line short circuit fault positioning method based on time domain inversion |
| CN113625200A (en) * | 2021-08-16 | 2021-11-09 | 哈尔滨工业大学 | In-situ detection method and detection device for performance degradation of coaxial connector |
| CN113625200B (en) * | 2021-08-16 | 2024-02-13 | 哈尔滨工业大学 | In-situ detection method and detection device for performance degradation of coaxial connector |
| CN114325244A (en) * | 2022-01-13 | 2022-04-12 | 清华大学 | Line short-circuit fault positioning method and system based on direct signal convolution |
| CN114325244B (en) * | 2022-01-13 | 2023-10-17 | 清华大学 | Line short-circuit fault positioning method and system based on direct signal convolution |
| CN114844776A (en) * | 2022-04-22 | 2022-08-02 | 电子科技大学 | Electromagnetic imaging positioning method and device based on time reversal time domain characteristics |
| CN114844776B (en) * | 2022-04-22 | 2023-04-18 | 电子科技大学 | Electromagnetic imaging positioning method and device based on time reversal time domain characteristics |
| CN114689997B (en) * | 2022-05-06 | 2023-09-26 | 国网江苏省电力有限公司宿迁供电分公司 | Distribution network cable fault identification and positioning method |
| CN114689997A (en) * | 2022-05-06 | 2022-07-01 | 国网江苏省电力有限公司宿迁供电分公司 | Distribution network cable fault identification and positioning method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108693444B (en) | 2019-08-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108693444B (en) | A kind of failure point of power transmission line localization method based on electromagnetic signal time reversal | |
| Meliopoulos et al. | Transient analysis of grounding systems | |
| Papalexopoulos et al. | Frequency dependent characteristics of grounding systems | |
| Mohanty et al. | An accurate noniterative fault-location technique for low-voltage DC microgrid | |
| He et al. | Natural frequency-based line fault location in HVDC lines | |
| Paolone et al. | Lightning electromagnetic field coupling to overhead lines: Theory, numerical simulations, and experimental validation | |
| Zhang et al. | EMTR‐based fault location for DC line in VSC‐MTDC system using high‐frequency currents | |
| Formisano et al. | Assessment of induced voltages in common and differential‐mode for a PV module due to nearby lightning strikes | |
| Dobakhshari | Noniterative parameter-free fault location on untransposed single-circuit transmission lines | |
| CN106054023A (en) | Method for estimating system impedances at two sides in single-ended range finding of power transmission line | |
| Yu et al. | A novel traveling wave fault location method for transmission network based on time linear dependence | |
| Razzaghi et al. | Electromagnetic time reversal applied to fault location in power networks | |
| Ametani | Lightning surge analysis by EMTP and numerical electromagnetic analysis method | |
| Brignone et al. | EM fields generated by a scale model helical antenna and its use in validating a code for lightning-induced voltage calculation | |
| Maslowski et al. | Experimental investigation and numerical modeling of surge currents in lightning protection system of a small residential structure | |
| Nie et al. | Convolution based time domain fault location method for lines in MMC-HVDC grids with distributed and frequency dependent line model | |
| Noda et al. | Error in propagation velocity due to staircase approximation of an inclined thin wire in FDTD surge simulation | |
| Barakou et al. | Online transient measurements of EHV cable system and model validation | |
| El Din et al. | A wavelet-based fault location technique for aged power cables | |
| He et al. | On the spatial resolution of fault-location techniques based on full-fault transients | |
| Wang et al. | Research on traveling wave fault technology based on ground potential | |
| Elkalashy et al. | Earth fault distance estimation using travelling waves provided with triacs‐based reclosing in distribution networks | |
| Anderson et al. | Algorithms for calculating of lightning induced voltages on distribution lines | |
| Zou et al. | Transient simulation model for a lightning protection system using the approach of a coupled transmission line network | |
| Pokharel et al. | Applications of time-domain numerical electromagnetic code to lightning surge analysis |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |