CN107064728B - The single-ended holographic frequency domain Fault Locating Method of ultra-high-tension power transmission line - Google Patents
The single-ended holographic frequency domain Fault Locating Method of ultra-high-tension power transmission line Download PDFInfo
- Publication number
- CN107064728B CN107064728B CN201610851981.XA CN201610851981A CN107064728B CN 107064728 B CN107064728 B CN 107064728B CN 201610851981 A CN201610851981 A CN 201610851981A CN 107064728 B CN107064728 B CN 107064728B
- Authority
- CN
- China
- Prior art keywords
- fault
- measurement point
- component
- phase
- voltage
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000005540 biological transmission Effects 0.000 title claims abstract description 15
- 230000007704 transition Effects 0.000 claims abstract description 13
- 238000005259 measurement Methods 0.000 claims description 83
- 238000005070 sampling Methods 0.000 claims description 23
- 238000001228 spectrum Methods 0.000 claims description 16
- 238000004422 calculation algorithm Methods 0.000 claims description 14
- 230000005611 electricity Effects 0.000 claims description 7
- 238000009434 installation Methods 0.000 claims description 3
- 206010034133 Pathogen resistance Diseases 0.000 claims 1
- 230000009897 systematic effect Effects 0.000 abstract description 8
- 230000001052 transient effect Effects 0.000 description 7
- 230000004807 localization Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 235000007926 Craterellus fallax Nutrition 0.000 description 1
- 240000007175 Datura inoxia Species 0.000 description 1
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Locating Faults (AREA)
- Medicines That Contain Protein Lipid Enzymes And Other Medicines (AREA)
Abstract
The single-ended holographic frequency domain Fault Locating Method of ultra-high-tension power transmission line, this method is from fault component network equation, if peer-to-peer system impedance, transition resistance and fault distance are undetermined parameter, the network equation of frequency domain under lumped parameter model is given, solution obtains system parameter and location information；The method of the present invention eliminates Systematic Errors existing for one-terminal data ranging from principle, improves range accuracy.
Description
Technical field
The present invention relates to power system transmission line technical field of relay protection, and in particular to ultra-high-tension power transmission line is single-ended complete
Cease frequency domain Fault Locating Method.
Background technique
UHV transmission network delivery capacity is huge, and line fault reparation plays the operational reliability of power grid extremely great
Effect.The quick reparation of line fault is decided by accurate fault location, and China's land resource is poor, and water power is sent outlet outside and walked
Corridor is with a varied topography, while line walking equipment is simple and crude, more makes the accuracy ever more important of fault location.
Summarize all kinds of algorithms of existing single-end electrical quantity, common issue be due to peer-to-peer system impedance and transition resistance not
Range equation caused by knowing is unsatisfactory for definite condition, it has to by assuming that supplementary condition realization fault location, and various hypothesis
It can not be suitable for the ever-changing operating status of power grid, so as to cause the generation of Systematic Errors.Here it is existing single-ended events
It is not accurate enough to hinder distance measuring method distance measurement result, is not able to satisfy the basic reason of engineer application requirement.
In recent years, intelligent substation was built is in full swing in China, and electronic mutual inductor is used widely.It has height
It the characteristics of acquisition rate, high linearity and high bandwidth, can electric signal of the true progress of disease.Due to electronic mutual inductor
Excellent Transfer characteristic, so that relay protection device is possible to using true transient fault signal, extracting more horn of plenty can
The fault message leaned on.If transient information itself can provide more location conditions, the deficiency of location condition can be supplemented, this
Sample just need not be made any it is assumed that and directly giving accurate positioning result by single-end information to system.
It to sum up analyzes, under the conditions of intelligent substation and electronic mutual inductor, application failure transient information, research has more
High-performance, the more accurate novel power transmission line fault range measurement principle of ranging are of great significance.
The technical problem to be solved in the invention has:
1, how fault transient information is utilized, more location conditions is supplemented, to avoid the systematicness of one-end fault ranging
Error reaches accurate ranging in principle.
2, the method for solving of the frequency domain range equation based on parameter identification thought.
3, the computing technique of signal spectrum.
Summary of the invention
In order to solve the above-mentioned technical problem, the purpose of the present invention is to provide a kind of single-ended holographic frequency domains of ultra-high-tension power transmission line
Fault Locating Method supplements more location conditions, to avoid the Systematic Errors of one-end fault ranging, reaches quasi- in principle
True ranging.
In order to achieve the above object, the present invention adopts the following technical scheme:
The single-ended holographic frequency domain Fault Locating Method of ultra-high-tension power transmission line, includes the following steps:
A, B, C three-phase voltage u of step 1, acquisition protection installation place current transformer and voltage transformerA(k),uB(k),
uC(k) and three-phase current iA(k),iB(k),iC(k)；
Step 2, to A, B, C three-phase voltage and electric current collected, utilize following formula calculate three-phase voltage and electric current therefore
Hinder component Δ uA(k),ΔuB(k),ΔuC(k), Δ iA(k),ΔiB(k),ΔiC(k) and the sampled value 3i of zero-sequence component0(k):
3i0(k)=Δ iA(k)+ΔiB(k)+ΔiC(k)
In above formula, k is sampling instant, and N is every cycle sampling number；A, B, C are calculated using complete cycle Fourier algorithm simultaneously
The real and imaginary parts of three-phase voltage and electric current cosine phasor；
Step 3 defines constant:
z1=r1+jωl1
z0=r0+jωl0
In formula, ω0- power frequency angular frequency；- m is divider ratio, can be equal to 2,3 or 4 etc.；r1,l1,r0,l0- line
Road unit length positive sequence, zero sequence resistance and inductance；Kr、Kl- zero-utility theory；
Step 4 calculates inductance and interior resistance parameter R in each mold component equivalence of this end systemMi,LMi(i=0,1,2), takes two
The three-phase current and three-phase voltage sampled value of a above different moments k obtains 0 modulus of Current Voltage, 1 mould using Clarke transform
Amount, 2 modulus instantaneous values are brought into system of linear equations
In, this equation group is solved using linear least-squares algorithm and obtains RMi,LMi(i=0,1,2)；
In formula, Δ iMiEach modulus fault component current instantaneous value of the side-M；ΔuMiEach modulus of measurement point M when-failure
Fault component voltage；DT is sampling interval duration；
Step 5 determines fault type
(1) it if singlephase earth fault, enables,
K2=-j ω UMAHC-ω0UMAHS
K3=j ω (IMAHCr1+IMAHSω0l1)+ω0(IMAHSr1-IMAHCω0l1)
In formula, UMAHC, UMAHSThe real part itself and the phase of imaginary part of-measurement point M failure phase normal duty voltage cosine phasor
Anti- number；IMAHC, IMAHSThe real part itself and the opposite number of imaginary part of-measurement point M failure phase normal duty electric current cosine phasor；
Calculate measurement point faulted phase current fault component Δ IMA(ω), failure phase busbar voltage fault component Δ UMA(ω)
And zero mould electric current frequency point ω spectrum component IM0(ω):
In formula, Δ IMAC、ΔIMASThe real part itself and the opposite number of imaginary part of-measurement point M failure phase fault electric current phasor；
ΔiMA(k)-measurement point M failure phase fault electric current component sampled value；ΔUMAC、ΔUMAS- measurement point M failure phase busbar voltage
The opposite number of real part of failure phasor itself and imaginary part；ΔuMA(k) sampling of-measurement point M failure phase busbar voltage fault component
Value；ΔIM0C、ΔIM0SThe real part itself and the opposite number of imaginary part of zero mould electric current of-measurement point M；ΔiMA(k)-zero mould of measurement point M
The sampled value of electric current；The sampling number of a cycle of N-；DT-sampling interval duration；K-is sampling instant；M is divider ratio；
A0=K1ΔUMA(ω)z0D-K2z0D
A1=-K1ΔUMA(ω)z0-K1az0D+K2z0-K3z0D
A2=-3K1IM0(ω)(RM0+jωLM0+z0D)
A3=K1ΔUMA(ω)-K2
A4=j ω A3
A5=-K1a-K3
A6=j ω A5
A7=-3K1IM0(ω)
A8=j ω A7
A9=K1z0a+K3z0
A=r1(ΔIMA(ω)+KrIM0(ω))+jωl1(ΔIMA(ω)+KlIM0(ω))
In formula, D-total track length；RM0、LM0- measurement point M side system zero sequence resistance and inductance；
(2) it if phase-to phase fault, enables,
K2=-j ω UMBCHC-ω0UMBCHS
K3=j ω (IMBCHCr1+IMBCHSω0l1)+ω0(IMBCHSr1-IMBCHCω0l1)
In formula, UMBCHC, UMBCHSThe real part of 2 modulus normal duty network load voltage cosine phasor of-measurement point M itself
With the opposite number of imaginary part；IMBCHC, IMBCHSThe real part itself and imaginary part of 2 modulus normal duty electric current cosine phasor of-measurement point M
Opposite number；
Calculate 2 modulus fault current component Δ I of measurement point MMBCThe spectrum component of (ω), false voltage component in frequency point ω
ΔUMBC(ω):
In formula, Δ IMBCC、ΔIMBCSThe real part itself and the opposite number of imaginary part of 2 modulus fault current phasor of-measurement point M；
ΔiMBC(k) -2 modulus fault current component of measurement point M sampled value；ΔUMBCC、ΔUMBCS2 modulus failure electricity of-measurement point M
Press the real part itself and the opposite number of imaginary part of phasor；ΔuMBC(k) -2 modulus false voltage component of measurement point M sampled value；
The sampling number of a cycle of N-；DT-sampling interval duration；
It enables,
A0=K1ΔUMBC(ω)z1D-K2z1D
A1=-K1ΔUMBC(ω)z1-K1az1D+K2z1-K3z1D
A3=K1ΔUMBC(ω)-K2
A4=-j ω A3
A5=-K1a-K3
A6=j ω A5
A8=j ω A7
A9=K1z1a+K3z1
A=z1ΔIMBC(ω)
In formula, D-total track length；RM2、LM2- measurement point M side system negative sequence resistance and inductance；
(3) it if three-phase fault, enables,
K2=-j ω UMAHC-ω0UMAHS
K3=j ω (IMAHCr1+IMAHSω0l1)+ω0(IMAHSr1-IMAHCω0l1)
In formula, UMAHC, UMAHSThe real part itself and the phase of imaginary part of-measurement point M failure phase normal duty voltage cosine phasor
Anti- number；IMAHC, IMAHSThe real part itself and the opposite number of imaginary part of-measurement point M failure phase normal duty electric current cosine phasor；
Calculate measurement point faulted phase current fault component Δ IMA(ω), failure phase busbar voltage fault component are frequency point ω's
Spectrum component Δ UMA(ω):
In formula, Δ IMAC、ΔIMASThe real part itself and the opposite number of imaginary part of-measurement point M failure phase fault electric current phasor；
ΔiMA(k)-measurement point M failure phase fault electric current component sampled value；ΔUMAC、ΔUMAS- measurement point M failure phase busbar voltage
The opposite number of real part of failure phasor itself and imaginary part；ΔuMA(k) sampling of-measurement point M failure phase busbar voltage fault component
Value；The sampling number of a cycle of N-；DT-sampling interval duration；
It enables,
A0=K1ΔUMA(ω)z1D-K2z1D
A1=-K1ΔUMA(ω)z1-K1az1D+K2z1-K3z1D
A2=-K1ΔIMA(ω)(RM1+jωLM1+z1D)
A3=K1ΔUMA(ω)-K2
A4=-j ω A3
A5=-K1a-K3
A6=j ω A5
A7=-K1ΔIMA(ω)
A8=j ω A7
A9=K1z1a+K3z1
A=z1ΔIMA(ω)
In formula, D-total track length；RM1、LM1- measurement point M side system positive sequence resistance and inductance
Step 6, range equation are as follows:
A0+A1d+A2RF+A3RN+A4LN+A5RNd+A6LNd+A7RFRN+A8RFLN+A9d2=0
In formula, distance of the d-fault point to measurement point M；RN,LNThe resistance and inductance parameters of-opposite side power supply；RF- failure
Point transition resistance；A in above-mentioned equationiType is calculated by the definition of step 5 according to different faults；
Under any frequency, fault current and the corresponding spectrum component of voltage strictly meet range equation, when ω is base
Wave frequency rate ω0When, 2 reality, imaginary part nonlinear equations are obtained according to range equation, when ω is other optional frequencies, according to survey
2 reality, imaginary part nonlinear equations can be also obtained away from equation, meet definite condition, using Newton iteration method or least square
The estimation technique solves equation group, obtains the estimated value of fault distance, transition resistance, opposite end power supply equivalent inductance and resistance parameter.
Compared to the prior art, the present invention has following advantage:
The prior art is all made of power frequency component, and unknown variable number is more than equation number, is unsatisfactory for definite condition, therefore usual root
Do appropriate according to system features it is assumed that eliminating Partial Variable.Processing certainly will cause systematic error in this way, because of real system
It is multifarious, it is specific to assume that all scenario is suitble to.This method is also used in signal other than using power frequency component
Other frequency components, therefore do not need to do system any hypothesis and can solve fault distance, there is no the systems of principle
Error.
Detailed description of the invention
Fig. 1 (a)-Fig. 1 (c) is transmission line fault network, and wherein Fig. 1 (a) is failure whole network, and Fig. 1 (b) is normal negative
Lotus network, Fig. 1 (c) are fault component network.
Specific embodiment
The present invention is described in further detail below in conjunction with the accompanying drawings.
The present invention is that a kind of ultra-high-tension power transmission line fault localization based on parameter recognition principle is new under the conditions of intelligent substation
Method, main purpose is to solve the problems, such as traditional one-terminal data fault distance-finding method there are Systematic Errors, in intelligent substation
Under the conditions of using electronic mutual inductor provide abundant transient information, with parameter identify principle carry out transmission line malfunction survey
Away from.
The basic ideas of parameter identification distance measuring method proposed by the invention are to obtain fault network equation, fault point side
On the basis of boundary's conditional equation, the current in the fault point and opposite end electric current, voltage in equation are eliminated, is obtained comprising home terminal current, electricity
Pressure, opposite side system impedance, fault point transition resistance and fault distance network equation, then utilize local terminal sample rate current, voltage,
The estimation for realizing unknown parameter in equation, to realize fault location.
R-L model is widely used in single-end electrical quantity Fault Location Algorithm, ignores the influence of distribution capacity.Under normal conditions,
For low pressure overhead transmission line or high pressure short-term road, approximation is carried out using R-L model, can achieve the requirement of engineer application.
Basic principle is introduced by two-shipper single loop line system model first, shown in model such as Fig. 1 (a)-Fig. 1 (c), it is assumed that system model is
Single-phase, route uses R-L model.
For Linear Network, failure whole network can be decomposed into normal duty network and fault component using principle of stacking
Network.Obtain measurement point M electric current, voltage failure component and fault point superimposed current, the network equation that voltage meets
It is calculated according to route two sides to the equal condition of fault point voltage, obtains route both ends electric current, voltage meets as off line
Network equation
For fault point, there are following failure boundary conditional equations
ΔiM+ΔiN=iFFormula (3)
In formula: Δ uMThe side-M voltage failure component instantaneous value；ΔiM,ΔiNThe side-M and N side line road current failure component wink
Duration；uFThe instantaneous value in-fault component network superimposed voltage source；iF- current in the fault point instantaneous value；D-fault point is to measurement point
The distance of M；D-total track length；R, l-route unit length resistance and inductance parameters；RM,LM,RN,LN- two sides power supply
Resistance and inductance parameters；RF- fault point transition resistance.
Fourier transform is carried out to formula (1) to formula (3) two sides, it is as follows to obtain corresponding frequency domain equation
ΔUM(ω)=rd Δ IM(ω)+jωldΔIM(ω)+RFIF(ω)+UF(ω) formula (4)
(RM+dr+jω(LM+dl))ΔIM(ω)=(RN+(D-d)r+jωLN+jω(D-d)l)ΔIN(ω) formula (5)
ΔIM(ω)+ΔIN(ω)=IF(ω) formula (6)
Simultaneous obtains eliminating the network equation of route opposite end electric current, voltage and current in the fault point, has following general type
f(ΔUM(ω),ΔIM(ω),d,RM,LM,RN,LN,RF,UF(ω))=0 formula (7)
In formula: the corresponding angular frequency of the investigated frequency point f of ω-；ΔUM(ω), Δ IMThe side (ω)-M voltage, current failure point
Measure the spectrum component in frequency point ω；ΔUN(ω), Δ INThe frequency spectrum of the side (ω)-N voltage, current failure component in frequency point ω divides
Amount；UFSpectrum component of the fault point superimposed voltage source in frequency point ω in (ω)-fault component network；IF(ω)-current in the fault point
In the spectrum component of frequency point ω.
Equation (7) is nonlinear equation.Wherein Δ UM(ω),ΔIM(ω) can calculate frequency spectrum by this side measurement data and obtain
It arrives；RM,LMFault component network can be solved by local terminal sampled data, be obtained by parameter identification method, is considered as known
Amount.
According to principle of stacking, uFTime-domain expression it is as follows
uF(t)=- uH(t) u (t)=- (UHCcosω0t+UHSsinω0T) u (t) formula (8)
Wherein, uHIt (t) is the voltage of fault point under normal operating conditions, u (t) represents unit-step function.UF(ω's)
Analytical expression is
With voltage phasor before failure at measurement point MAnd electric current phasorRelationship be
To obtain
It can be seen that each variable in formula (7), removes RN,LN,RF, other than d, remaining variables can be before the failure by measurement point M
Measurement data afterwards is calculated, or is expressed as expression formula (such as U of local measurement data and fault distanceF(ω)), therefore can
With the range equation that abbreviation is containing only following known variables
f(d,RN,LN,RF)|ω=0 formula (12)
After failure, it can be established with the Current Voltage frequency spectrum that fault transient process sampled data is calculated in different frequent points
2 or 2 or more range equations corresponding with formula (12).Since formula (12) are complex number equation, thus available 4 or 4
The equation group of above real number equation composition, solves this Nonlinear System of Equations by least square optimization, available
Fault distance d, while can also be in the hope of unknown parameter R in R-L modelN,LN,RFValue.
Actual route is three-phase system, and location algorithm of the invention is introduced for singlephase earth fault occurs.Still with
Fig. 1 (a)-Fig. 1 (c) Lai Jinhang algorithmic derivation, but the figure should regard three-phase system as.Assuming that singlephase earth fault occurs for fault point,
And assume that fault type is A phase ground fault.Then for fault component network, electric current, the electricity of available measurement point M failure phase
The network equation for pressing fault component to meet with fault point superimposed current, voltage is as follows
In formula: rs,rmResistance and the mutual resistance certainly of-route unit length；ls,lmThe self-induction and mutual inductance of-route unit length；
ΔuMA(t)-measurement point M fault component instantaneous voltage；ΔiMA(t),ΔiMB(t),ΔiMC(t)-measurement point M three-phase current
Fault component instantaneous value；uFA(t)-fault component network fault point superimposed voltage instantaneous value；iFA(t)-fault point fault current
Instantaneous value.
Notice r1=rs-rm, l1=ls-lmWith Δ iMA+ΔiMB+ΔiMC=3 Δ iM0=3iM0And define zero sequence compensation system
NumberFormula (13) is rewritten as
For zero lay wire network, calculate that the equal condition of voltage to fault point obtains by route two sides protection installation place
For A phase ground fault, current in the fault point meets such as downstream condition
iFA(t)=3 (iM0(t)+iN0(t)) formula (16)
Fourier transform is carried out, it is as follows that corresponding frequency domain equation can be obtained
ΔUMA(ω)=UFA(ω)+IFA(ω)RF+r1d(ΔIMA(ω)+krIM0(ω))+jωl1d(ΔIMA(ω)+
klIM0(ω)) formula (17)
(RM0+r0d+jω(LM0+l0d))IM0(ω)=(RN0+r0D-r0d+jω(LN0+l0D-l0d))IN0(ω) formula (18)
IFA(ω)=3 (IM0(ω)+IN0(ω)) formula (19)
In formula: the π of the ω=2 f-corresponding angular frequency of frequency f；ΔUMA(ω),ΔIMA(ω),IM0(ω)-measurement point bus
The spectrum component of voltage failure component, failure phase A phase current fault component and zero mould electric current in frequency point ω；IFA(ω),IN0
The spectrum component of (ω)-fault point fault current and the side N zero-sequence current in frequency point ω；r1,l1,r0,l0- route unit length
Positive sequence, zero sequence resistance and inductance；RN0=RNs+2RNm,LN0=LNs+2LNm- N side system zero sequence resistance and inductance；RNs,RNm,LNs,
LNm- N side system is each mutually from resistance, self-induction and alternate mutual resistance, mutual inductance.
The load voltage electric current phasor that the entire spectrum and measurement point M of fault component network fault point superimposed voltage operate normally
Between relationship be
It is above-mentioned
UFAC=-UMAHC+IMAHCr1d+IMAHSω0l1d
UFAS=-UMAHS+IMAHSr1d-IMAHCω0l1D formula (21)
In formula: UFAC, UFASThe real part itself and the phase of imaginary part of-fault component network fault point superimposed voltage cosine phasor
Anti- number；UMAHC, UMAHSThe real part itself and the opposite number of imaginary part of-measurement point M failure phase normal duty voltage cosine phasor；
IMAHC, IMAHSThe real part itself and the opposite number of imaginary part of-measurement point M failure phase normal duty electric current cosine phasor.
Define constant
K2=-j ω UMAHC-ω0UMAHS
K3=j ω (IMAHCr1+IMAHSω0l1)+ω0(IMAHSr1-IMAHCω0l1)
It is collated to obtain following nonlinear equation
A0+A1d+A2RF+A3RN0+A4LN0+A5RN0d+A6LN0d+A7RFRN0+A8RFLN0+A9d2=0 formula (22)
It is above-mentioned
A0=K1ΔUMA(ω)z0D-K2z0D
A1=-K1ΔUMA(ω)z0-K1az0D+K2z0-K3z0D
A2=-3K1IM0(ω)(RM0+jωLM0+z0D)
A3=K1ΔUMA(ω)-K2, A4=j ω A3
A5=-K1a-K3, A6=j ω A5
A7=-3K1IM0(ω), A8=j ω A7
A9=K1z0a+K3z0, z0=r0+jωl0
A=r1(ΔIMA(ω)+KrIM0(ω))+jωl1(ΔIMA(ω)+KlIM0(ω))
Under any frequency, fault current and the corresponding spectrum component of voltage strictly meet formula (22), certainly for base
Wave, each harmonic and Fractional Frequency are equally set up.When ω is fundamental frequency ω0When, can be obtained by stable state network reality, imaginary part 2 it is non-
Linear equation.When ω is any other frequency, such as ω=ω0/ 2 or ω0/ 4, from transient network equation also can be obtained it is real,
2 nonlinear equations of imaginary part can be used Newton iteration method or least squares estimate solve the party to meet definite condition
Journey group obtains the estimated value of the parameters such as fault distance.For two-phase short-circuit fault and three phase short circuit fault, location algorithm is pushed away
It is similar to lead process.
For the correctness and validity of verification method, the emulation of fault localization is carried out to Beijing-Tianjin-Tangshan 500kV transmission line of electricity, is
Unite wiring such as Fig. 1 (a)-Fig. 1 (c), sample frequency 500kHz, and simulation system parameters are as follows.
Line length: D=300km
Line parameter circuit value: r1=0.02083 Ω/km；l1=0.8984mH/km；
r0=0.1148 Ω/km；l0=2.2886mH/km；
M side system parameter: RM1=1.0515 Ω；LM1=0.13743H；
RM0=0.6 Ω；LM0=0.0926H；
N side system parameter；RN1=26 Ω；LN1=0.14298H；
RN0=20 Ω；LN0=0.11927H；The advanced N side system of M side system
Range error:
Table 1 simulate route different location occur it is single-phase through 200 Ω transition resistance failure when distance measurement result.It can from table 1
To find out, under the practical operation situation deviation algorithm of the system assumed condition almost the same to two sides system impedance angle, this hair
Bright range accuracy under various fault ' conditions is all higher, can satisfy requirement of engineering.
Distance measurement result (R of 1 different faults of table underF=200 Ω)
Table 2 simulates in the line away from generation singlephase earth fault, sheet when transition resistance takes different value at measurement point 200km
The distance measurement result of invention.From table 2 it can be seen that using modal identification algorithm carry out ranging, distance measurement result almost with transition resistance
It is unrelated, and range accuracy is also to meet engine request.
The asynchronous emulation distance measurement result (d=200km) of 2 fault point transition resistance of table
Table 3 gives when different types of faults occur for route different location, uses the emulation of parameter identification location algorithm
As a result, wherein transition resistance is used uniformly 50 Ω.From table 3 it can be seen that parameter identifies location algorithm in different faults type feelings
Accurate distance measurement result can be provided under condition.
Emulation distance measurement result (RF=50 Ω) when 3 different faults type of table
To sum up, traditional power frequency single end distance measurement algorithm can only list two equations, and range equation owes fixed, all directly or
Indirectly to variable RN,LN,RFCarried out simplify it is assumed that thus there are Systematic Errors.The present invention is by RN,LN,RFAlso as to
Identification parameter, this is different with traditional power frequency single end distance measurement algorithm, fundamentally avoids the systematic error of single end distance measurement.
Simulation result also demonstrates the present invention can overcome systematic error from principle, and distance measurement result is not by peer-to-peer system impedance and transition
The influence of resistance adapts to the requirement applied in intelligent substation.
Claims (1)
1. the single-ended holographic frequency domain Fault Locating Method of ultra-high-tension power transmission line, characterized by the following steps:
A, B, C three-phase voltage u of step 1, acquisition protection installation place current transformer and voltage transformerA(k),uB(k),uC(k)
With three-phase current iA(k),iB(k),iC(k)；
Step 2, to A, B, C three-phase voltage and electric current collected, utilize following formula to calculate three-phase voltage and failure of the current point
Measure Δ uA(k),ΔuB(k),ΔuC(k), Δ iA(k),ΔiB(k),ΔiC(k) and the sampled value 3i of zero-sequence component0(k):
3i0(k)=Δ iA(k)+ΔiB(k)+ΔiC(k)
In above formula, k is sampling instant, and N is every cycle sampling number；A, B, C three-phase are calculated using complete cycle Fourier algorithm simultaneously
The real and imaginary parts of voltage and current cosine phasor；
Step 3 defines constant:
z1=r1+jωl1
z0=r0+jωl0
In formula, ω0- power frequency angular frequency；- m is divider ratio, is equal to 2,3 or 4；r1- route unit length positive sequence electricity
Resistance, l1- route unit length positive sequence inductance, r0- route unit length zero sequence resistance, l0- route unit length zero sequence inductance；
Kr、Kl- zero-utility theory；
Step 4 calculates inductance and interior resistance parameter in each mold component equivalence of this end system, including measurement point M positive sequence resistance and electricity
Feel RM1,LM1；Measurement point M negative sequence resistance and inductance RM2,LM2；Measurement point M zero sequence resistance and inductance RM0,LM0；Take two or more not
The three-phase current of k and three-phase voltage sampled value in the same time obtain 0 modulus of Current Voltage, 1 modulus, 2 modulus using Clarke transform
Instantaneous value is brought into system of linear equationsIn, wherein i=0,
1,2, this equation group is solved using linear least-squares algorithm and obtains RMi,LMi；
In formula, Δ iMiEach modulus fault component current instantaneous value of-measurement point M；ΔuMiEach modulus of measurement point M when-failure
Fault component instantaneous voltage；DT is sampling interval duration；
Step 5 determines fault type
(1) it if singlephase earth fault, enables,
K2=-j ω UMAHC-ω0UMAHS
K3=j ω (IMAHCr1+IMAHSω0l1)+ω0(IMAHSr1-IMAHCω0l1)
In formula, UMAHC, UMAHSThe real part itself and the opposite number of imaginary part of-measurement point M failure phase normal duty voltage cosine phasor；
IMAHC, IMAHSThe real part itself and the opposite number of imaginary part of-measurement point M failure phase normal duty electric current cosine phasor；
Calculate measurement point faulted phase current fault component Δ IMA(ω), failure phase busbar voltage fault component Δ UMA(ω) and zero mould
Spectrum component I of the electric current in frequency point ωM0(ω):
In formula, Δ IMAC、ΔIMASThe real part itself and the opposite number of imaginary part of-measurement point M failure phase fault electric current phasor；ΔiMA
(k)-measurement point M failure phase fault electric current component sampled value；ΔUMAC、ΔUMAS- measurement point M failure phase busbar voltage failure
The opposite number of real part of phasor itself and imaginary part；ΔuMA(k)-measurement point M failure phase busbar voltage fault component sampled value；
IM0C、IM0SThe real part itself and the opposite number of imaginary part of zero mould electric current of-measurement point M；iM0(k) sampling of-zero mould electric current of measurement point M
Value；The sampling number of the every cycle of N-；DT-sampling interval duration；K-is sampling instant；M is divider ratio；
A0=K1ΔUMA(ω)z0D-K2z0D
A1=-K1ΔUMA(ω)z0-K1az0D+K2z0-K3z0D
A2=-3K1IM0(ω)(RM0+jωLM0+z0D)
A3=K1ΔUMA(ω)-K2
A4=j ω A3
A5=-K1a-K3
A6=j ω A5
A7=-3K1IM0(ω)
A8=j ω A7
A9=K1z0a+K3z0
A=r1(ΔIMA(ω)+KrIM0(ω))+jωl1(ΔIMA(ω)+KlIM0(ω))
In formula, D-total track length；
(2) it if phase-to phase fault, enables,
K2=-j ω UMBCHC-ω0UMBCHS
K3=j ω (IMBCHCr1+IMBCHSω0l1)+ω0(IMBCHSr1-IMBCHCω0l1)
In formula, UMBCHC, UMBCHSThe real part of the 2 modulus normal duty voltage cosine phasors of-measurement point M itself is opposite with imaginary part
Number；IMBCHC, IMBCHSThe real part itself and the opposite number of imaginary part of the 2 modulus normal duty electric current cosine phasors of-measurement point M；
Calculate the 2 modulus fault current component Δ I of measurement point MMBCThe spectrum component Δ of (ω), false voltage component in frequency point ω
UMBC(ω):
In formula, Δ IMBCC、ΔIMBCSThe real part itself and the opposite number of imaginary part of the 2 modulus fault current phasors of-measurement point M；Δ
iMBC(k) sampled value of 2 modulus fault current component of-measurement point M；ΔUMBCC、ΔUMBCSThe 2 modulus failure electricity of-measurement point M
Press the real part itself and the opposite number of imaginary part of phasor；ΔuMBC(k) sampled value of 2 modulus false voltage component of-measurement point M；
The sampling number of a cycle of N-；DT-sampling interval duration；
It enables,
A0=K1ΔUMBC(ω)z1D-K2z1D
A1=-K1ΔUMBC(ω)z1-K1az1D+K2z1-K3z1D
A3=K1ΔUMBC(ω)-K2
A4=-j ω A3
A5=-K1a-K3
A6=j ω A5
A8=j ω A7
A9=K1z1a+K3z1
A=z1ΔIMBC(ω)
In formula, D-total track length；
(3) it if three-phase fault, enables,
K2=-j ω UMAHC-ω0UMAHS
K3=j ω (IMAHCr1+IMAHSω0l1)+ω0(IMAHSr1-IMAHCω0l1)
In formula, UMAHC, UMAHSThe real part itself and the opposite number of imaginary part of-measurement point M failure phase normal duty voltage cosine phasor；
IMAHC, IMAHSThe real part itself and the opposite number of imaginary part of-measurement point M failure phase normal duty electric current cosine phasor；
Calculate measurement point faulted phase current fault component Δ IMAThe frequency spectrum of (ω), failure phase busbar voltage fault component in frequency point ω
Component Δ UMA(ω):
In formula, Δ IMAC、ΔIMASThe real part itself and the opposite number of imaginary part of-measurement point M failure phase fault electric current phasor；ΔiMA
(k)-measurement point M failure phase fault electric current component sampled value；ΔUMAC、ΔUMAS- measurement point M failure phase busbar voltage failure
The opposite number of real part of phasor itself and imaginary part；ΔuMA(k)-measurement point M failure phase busbar voltage fault component sampled value；
The sampling number of a cycle of N-；DT-sampling interval duration；
It enables,
A0=K1ΔUMA(ω)z1D-K2z1D
A1=-K1ΔUMA(ω)z1-K1az1D+K2z1-K3z1D
A2=-K1ΔIMA(ω)(RM1+jωLM1+z1D)
A3=K1ΔUMA(ω)-K2
A4=-j ω A3
A5=-K1a-K3
A6=j ω A5
A7=-K1ΔIMA(ω)
A8=j ω A7
A9=K1z1a+K3z1
A=z1ΔIMA(ω)
In formula, D-total track length；
Step 6, range equation are as follows:
A0+A1d+A2RF+A3RN+A4LN+A5RNd+A6LNd+A7RFRN+A8RFLN+A9d2=0
In formula, distance of the d-fault point to measurement point M；RN,LNThe resistance and inductance parameters of-opposite side power supply；RF- fault point mistake
Cross resistance；A in above-mentioned range equationiType is calculated by the definition of step 5 according to different faults；
Under any frequency, fault current and the corresponding spectrum component of voltage strictly meet range equation, when ω is power frequency angle
Frequencies omega0When, 2 reality, imaginary part nonlinear equations are obtained according to range equation, when ω is other optional frequencies, according to ranging
Equation can also obtain 2 reality, imaginary part nonlinear equations, meet definite condition, estimated using Newton iteration method or least square
Meter method solves equation group, obtains the estimated value of fault distance, transition resistance, the resistance of opposite side power supply and inductance parameters.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610851981.XA CN107064728B (en) | 2016-09-26 | 2016-09-26 | The single-ended holographic frequency domain Fault Locating Method of ultra-high-tension power transmission line |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610851981.XA CN107064728B (en) | 2016-09-26 | 2016-09-26 | The single-ended holographic frequency domain Fault Locating Method of ultra-high-tension power transmission line |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107064728A CN107064728A (en) | 2017-08-18 |
| CN107064728B true CN107064728B (en) | 2019-10-25 |
Family
ID=59617070
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610851981.XA Active CN107064728B (en) | 2016-09-26 | 2016-09-26 | The single-ended holographic frequency domain Fault Locating Method of ultra-high-tension power transmission line |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107064728B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107979076B (en) * | 2017-11-30 | 2022-06-10 | 中国电力科学研究院有限公司 | Power transmission line protection method and system based on transient component frequency spectrum characteristics |
| CN109884471B (en) * | 2019-03-22 | 2020-12-11 | 张安斌 | Phasor product active method for judging fault phase and selecting line |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3244975A (en) * | 1961-12-20 | 1966-04-05 | Siemens Ag | Apparatus for determining disturbances along electrical lines utilizing the impulse reflection method including adjustable equalizing means and stroboscopic sampling |
| JPS59135377A (en) * | 1983-01-24 | 1984-08-03 | Fuji Electric Co Ltd | Method for evaluating grounding fault point of three- phase power transmission distribution line |
| CN101975910A (en) * | 2010-09-07 | 2011-02-16 | 昆明理工大学 | Intelligent fault classification and location method for ultra-high voltage direct current transmission line |
| CN102353875A (en) * | 2011-09-30 | 2012-02-15 | 山东理工大学 | Fault distance measuring method for combined travelling wave of power transmission line |
| CN103018636A (en) * | 2012-12-14 | 2013-04-03 | 昆明理工大学 | Cable single-ended travelling wave distance measurement method using fault signature frequency bands and TT conversion |
| CN103091603A (en) * | 2013-01-14 | 2013-05-08 | 华北电力大学 | Breakdown intelligent classification and positioning method of electric transmission line |
| CN103207351A (en) * | 2013-03-12 | 2013-07-17 | 西安工程大学 | Transmission line fault locating method based on reclosure |
| CN103245893A (en) * | 2013-04-10 | 2013-08-14 | 昆明理工大学 | Radial distribution network layered and distributed ANN fault location method based on natural frequency |
| CN103777114A (en) * | 2014-01-26 | 2014-05-07 | 浙江大学 | Method for recognizing single-phase permanent fault of single-ended electric transmission line with paralleling reactor |
| EP2741390A1 (en) * | 2012-12-06 | 2014-06-11 | Schneider Electric Industries SAS | Directional detection of a fault, in particular in a network with compensated or insulated neutral |
| CN105223468A (en) * | 2015-09-29 | 2016-01-06 | 湖南大学 | Based on the transmission line of electricity one-end fault ranging method of mapping function |
-
2016
- 2016-09-26 CN CN201610851981.XA patent/CN107064728B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3244975A (en) * | 1961-12-20 | 1966-04-05 | Siemens Ag | Apparatus for determining disturbances along electrical lines utilizing the impulse reflection method including adjustable equalizing means and stroboscopic sampling |
| JPS59135377A (en) * | 1983-01-24 | 1984-08-03 | Fuji Electric Co Ltd | Method for evaluating grounding fault point of three- phase power transmission distribution line |
| CN101975910A (en) * | 2010-09-07 | 2011-02-16 | 昆明理工大学 | Intelligent fault classification and location method for ultra-high voltage direct current transmission line |
| CN102353875A (en) * | 2011-09-30 | 2012-02-15 | 山东理工大学 | Fault distance measuring method for combined travelling wave of power transmission line |
| EP2741390A1 (en) * | 2012-12-06 | 2014-06-11 | Schneider Electric Industries SAS | Directional detection of a fault, in particular in a network with compensated or insulated neutral |
| CN103018636A (en) * | 2012-12-14 | 2013-04-03 | 昆明理工大学 | Cable single-ended travelling wave distance measurement method using fault signature frequency bands and TT conversion |
| CN103091603A (en) * | 2013-01-14 | 2013-05-08 | 华北电力大学 | Breakdown intelligent classification and positioning method of electric transmission line |
| CN103207351A (en) * | 2013-03-12 | 2013-07-17 | 西安工程大学 | Transmission line fault locating method based on reclosure |
| CN103245893A (en) * | 2013-04-10 | 2013-08-14 | 昆明理工大学 | Radial distribution network layered and distributed ANN fault location method based on natural frequency |
| CN103777114A (en) * | 2014-01-26 | 2014-05-07 | 浙江大学 | Method for recognizing single-phase permanent fault of single-ended electric transmission line with paralleling reactor |
| CN105223468A (en) * | 2015-09-29 | 2016-01-06 | 湖南大学 | Based on the transmission line of electricity one-end fault ranging method of mapping function |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107064728A (en) | 2017-08-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103792465B (en) | A kind of method of the range finding of the one-phase earthing failure in electric distribution network based on residual voltage | |
| CN103941147B (en) | Utilize the distribution network cable single-phase ground fault distance measuring method of transient state principal frequency component | |
| CN103399209B (en) | A kind of extra-high voltage bipolar direct current transmission line power frequency parameters measuring method | |
| CN108802564B (en) | distribution network T-shaped line parameter-independent fault location algorithm and system | |
| CN102135571B (en) | Anti-interference measurement method for zero sequence impedance of super-high-voltage/ultrahigh-voltage multi-loop power transmission line | |
| CN109541392A (en) | A kind of one-end fault ranging method suitable for flexible HVDC transmission system | |
| CN102445638B (en) | Time-domain fault location method of multibranch direct current line in multiterminal direct current power transmission system | |
| CN106680636B (en) | A kind of double back does not hold power transmission line zero-sequence parameter accurate measurement method altogether | |
| CN109490706A (en) | A kind of multi branch electric power lines road Fault Locating Method | |
| CN106501612B (en) | Based on the double-back portions of non uniform transmission line with tower route zero sequence parameter measurement method | |
| CN109142974A (en) | A kind of alternating current circuit one-end fault ranging method suitable for flexible direct current feed-in | |
| CN107064728B (en) | The single-ended holographic frequency domain Fault Locating Method of ultra-high-tension power transmission line | |
| Zhang et al. | Voltage-sag-profiles-based fault location in high-speed railway distribution system | |
| CN103424627B (en) | The method of double-end measurement parallel grid line zero-sequence impedance | |
| CN106405337B (en) | Inverse distributed power accesses the Fault Locating Method of power distribution network | |
| Dutta et al. | Assessment of impedance based fault locator for AC micro-grid | |
| CN104865452B (en) | Not transposed transmission line electric capacity anti-jamming measurement methods based on harmonic component | |
| Li et al. | Augmented state estimation method for fault location based on on-line parameter identification of PMU measurement data | |
| Qianqian et al. | A new smart distribution grid fault self-healing system based on traveling-wave | |
| CN106786498B (en) | Master station-transformer substation data collaborative identification method and device | |
| Zhao et al. | A precise fault location method independent of line parameters based on positive and negative sequence fault components | |
| Sun et al. | Generalized impedance-based fault distance calculation mehtod for power distribution systems | |
| CN113009275A (en) | Double-end fault location method for flexible direct-current access alternating-current hybrid line | |
| Wang et al. | A fault location method for multi-branch distribution lines | |
| Lin et al. | RL model algorithm-based fault location for DC distribution networks |
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 |