CN109188189A - Based on arc characteristic ultra-high/extra-high voltage transmission line of electricity permanent fault recognition methods - Google Patents

Based on arc characteristic ultra-high/extra-high voltage transmission line of electricity permanent fault recognition methods Download PDF

Info

Publication number
CN109188189A
CN109188189A CN201810757496.5A CN201810757496A CN109188189A CN 109188189 A CN109188189 A CN 109188189A CN 201810757496 A CN201810757496 A CN 201810757496A CN 109188189 A CN109188189 A CN 109188189A
Authority
CN
China
Prior art keywords
frequency
voltage
fault
value
transmission line
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810757496.5A
Other languages
Chinese (zh)
Other versions
CN109188189B (en
Inventor
李博通
赵盛
李斌
姚斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Electric Power Research Institute of State Grid Henan Electric Power Co Ltd
Original Assignee
Tianjin University
Electric Power Research Institute of State Grid Henan Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tianjin University, Electric Power Research Institute of State Grid Henan Electric Power Co Ltd filed Critical Tianjin University
Priority to CN201810757496.5A priority Critical patent/CN109188189B/en
Publication of CN109188189A publication Critical patent/CN109188189A/en
Application granted granted Critical
Publication of CN109188189B publication Critical patent/CN109188189B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/081Locating faults in cables, transmission lines, or networks according to type of conductors
    • G01R31/085Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution lines, e.g. overhead
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/088Aspects of digital computing

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 relates to power transmission and transforming equipments, for the accurate automatic identification for realizing power transmission and transforming equipment failure.The present invention is based on arc characteristic ultra-high/extra-high voltage transmission line of electricity permanent fault recognition methods, steps are as follows: step A: when singlephase earth fault occurs for system, after circuit breaker trip, extract the three-phase voltage and electric current at transmission line of electricity both ends measurement point (m, n), the modulus voltage and modulus electric current of route both ends measurement point are obtained through phase-model transformation, using Clarke transform;Step B: the value step C of fault point arc current feature expression is calculated using Two-Terminal Electrical Quantities: to ux(t) EMD decomposition is carried out, and finds out the marginal spectrum energy value of each frequency component;Step D: the ratio between high frequency spectrum energy and fundamental frequency spectrum energy P are found outR, with setting valve PsetCompare, works as PRLess than PsetWhen, it is determined as permanent fault, otherwise is transient fault.Present invention is mainly applied to the accurate automatic identification occasions of power transmission and transforming equipment failure.

Description

Based on arc characteristic ultra-high/extra-high voltage transmission line of electricity permanent fault recognition methods
Technical field
The present invention relates to power transmission and transforming equipments, specifically, are related to based on the permanently event of arc characteristic ultra-high/extra-high voltage transmission line of electricity Hinder recognition methods.
Background technique
The operating experience of electric system shows that overhead transmission line failure is all " instantaneity " mostly, and singlephase earth fault Account for about 70%~80%, therefore single phase autoreclosing is widely applied in line protection.It is single-phase automatic Reclosing is overlapped after transmission line of electricity tripping regardless of whether failure disappears, if overlapping to transient fault, failure phase Route will restore electricity;If overlapping to permanent fault, tripping again and being no longer overlapped.Such reclosing mode cannot judge Nature of trouble is not only unable to the normal power supply of recovery system when coinciding with permanent fault, and bring secondary pulse may be used also The stable operation of system be able to can be caused to seriously affect, therefore, accurate Judging fault property is after system failure tripping to realize Adaptive reclosing is of great significance.
When singlephase earth fault occurs for system, it is instantaneous/two kinds of permanent fault in the case of arc characteristic have differences. When permanent fault, electric arc is very short there are the time.When transient fault, the arc duration is longer and exists in the secondary arc phase "-blow-out of restriking " phenomenon.It can be seen that including fault message abundant in fault electric arc.However most of fault identification so far Method is all mutually to restore the characteristic expansion of voltage around tripping, although such method can accurately determine nature of trouble, It is that could start decision process after needing the complete blow-out of electric arc, it is possible to lead to system cut-off overlong time and cause biggish Economic loss.For traditional transmission line of electricity, system be can only obtain electrical at both-end bus after failure occurs How amount using route both ends electrical data is calculated the arc information of fault point to accurate Judging fault property, With important researching value.
Summary of the invention
In order to overcome the deficiencies of the prior art, the present invention is directed to propose ultra-high/extra-high voltage transmission line of electricity permanent fault identification side Method realizes the accurate automatic identification of power transmission and transforming equipment failure.For this purpose, the technical solution adopted by the present invention is that, it is based on arc characteristic Ultra-high/extra-high voltage transmission line of electricity permanent fault recognition methods, steps are as follows:
Step A: it after system generation singlephase earth fault, circuit breaker trip, extracts transmission line of electricity both ends measurement point (m, n) The three-phase voltage and electric current at place obtain the modulus voltage and modulus electric current of route both ends measurement point through phase-model transformation, using gram Clarke transformation, phase mode transformation matrix T are specific as follows:
Voltage, current transformation process are as follows:
In formula, the three-phase voltage and electric current at the both ends mn are respectively umABC、unABC、imABC, inABC, the both ends mn modulus voltage and Modulus electric current is respectively um123、un123、im123、in123
Step B: the value of fault point arc current feature expression is calculated using Two-Terminal Electrical Quantities
Rectify the calculation formula to voltage traveling wave amount and the end n backward voltage traveling wave amount using m are as follows:
[um(t)+zc(t)*im(t)]*amn(t)
-[un(t)-zc(t)*in(t)]
=zc(t)*ik(t)*ank(t)
Rectify the calculation formula to voltage traveling wave amount using the end m backward voltage traveling wave amount and n are as follows:
[um(t)-zc(t)*im(t)]
-[un(t)+zc(t)*in(t)]*amn(t)
=zc(t)*ik(t)*amk(t)
In formula, wherein amkIt (t) is the time-domain expression of the propagation constant at fault point k to m side bus, ankIt (t) is failure The time-domain expression of propagation constant at point k to n side bus, amnIt (t) is propagation constant time domain of the m side bus at n side bus Expression formula.zcIt (t) is the time-domain expression of surge impedance of a line, um(t) the modulus voltage to be obtained according to the end route m phase voltage value Value, un(t) the modulus voltage value to be obtained according to the end route n phase voltage value.ik(t) be fault point arc current value.
U obtained in selecting step Am1、un1、im1、in1Ground modulus value substitutes into calculation formula, obtains the calculating of feature expression Value ux(t);
Step C: to ux(t) EMD decomposition is carried out, and finds out the marginal spectrum energy value of each frequency component.
Step D: the ratio between high frequency spectrum energy and fundamental frequency spectrum energy P are found outR, with setting valve PsetCompare, works as PRLess than PsetWhen, It is determined as permanent fault, otherwise is transient fault.
Further, step C: to ux(t) EMD decomposition is carried out, and finds out the marginal spectrum energy value of each frequency component;
EMD decomposable process is as follows:
(1) signal u is found outx(t) all maximum point and minimum point use cubic spline letter to obtained extreme point Number is fitted to the upper and lower envelope of former data sequence, calculates the mean value of envelope up and down, is denoted as m1(t), former data sequence ux(t) Subtract m1(t) the new data sequence h for removing low frequency component is obtained1(t), i.e.,
ux(t)-m1(t)=h1(t) (10)
Repeat above-mentioned treatment process k times, until h1(t) the definition requirement for meeting IMF drives off the 1st in this way IMF component c1(t), c1It (t) is ux(t) component of highest frequency in;
(2) by c1(t) from ux(t) it is separated in get the difference signal r of high fdrequency component is removed to one1(t), that is, have
r1(t)=ux(t)-c1(t) (11)
By r1(t) it is used as initial data, repeats step (1) (2), obtains second IMF component c2(t), n times are repeated, until rn(t) become monotonic function or amplitude is small to stopping decomposing when can ignore, obtain n IMF component and a residual components rn (t), it may be assumed that
Wherein, rnIt (t) is survival function, the average tendency of representation signal, and each IMF component c1(t)~cn(t) respectively The ingredient of different frequency sections from high to low is contained, after carrying out Hilbert transformation to each IMF component, is obtained
Wherein ai(t) c is indicatedi(t) amplitude, ωi(t) it indicates instantaneous frequency, can further be defined with Hilbert spectrum Marginal spectrum is
In an example, sample frequency 10000Hz, highest frequency is 5000Hz in marginal spectrum, if frequency 44Hz with Under component be low-frequency range, 45Hz~55Hz be fundamental frequency section, 56Hz-5000Hz is high band, indicates it using marginal spectral amplitude ratio The energy value of respective frequencies, if the energy value for the component that frequency is i is Pi, low frequency energy, fundamental frequency energy, high-frequency energy is respectively Pd,PjAnd Pg, calculation formula is such as shown in (15):
Further, step D: finding out the ratio between high frequency spectrum energy and fundamental frequency spectrum energy, with setting valve PsetCompare:
According to P obtained in step CjAnd PgCalculate the ratio between high frequency spectrum energy and fundamental frequency spectrum energy, calculation formula PR= Pg/Pj
In an example, PsetIt is 0.5.As the P being calculatedRWhen less than 0.5, failure be it is permanent, work as PRIt is greater than When 0.5, failure is instantaneity.
The features of the present invention and beneficial effect are:
1. compared to based on the fault recognition method for restoring voltage characteristic, this method needs not wait for electric arc and extinguishes it completely Just start entire determination flow afterwards, is conducive to the quick discrimination of nature of trouble.
2. compared to the fault recognition method of failure phase terminal voltage during directly analysis electric arc arcing, this method passes through both-end Electrical quantity has calculated the characteristic information of fault point and has distinguished nature of trouble with this, and criterion principle is tighter, determines result more Add reliable.
3. this method can be adaptive compared to the conventional method for the frequency content extracted in signal using Fourier decomposition Decomposition non-stationary signal, obtained each frequency content is more accurate, and criterion is with a high credibility.
Detailed description of the invention:
Fig. 1 is extra high voltage network system schematic.
System schematic after the tripping of Fig. 2 transmission line one-phase earth fault.
U when Fig. 3 is transient faultx(t) marginal spectrum.
U when Fig. 4 is permanent faultx(t) marginal spectrum.
Fig. 5 is the transmission line of electricity permanent fault recognition methods flow chart based on fault electric arc characteristic.
Specific embodiment
To solve the problems of the prior art, the invention proposes a kind of ultra-high/extra-high voltage based on fault electric arc characteristic is defeated Electric line permanent fault recognition methods.This method is directed to ultrahigh-voltage alternating-current transmission line of electricity, using transmission line travelling wave principle, structure The feature expression of fault point arc current is made.This feature expression is carried out using the electrical data at transmission line of electricity both ends It calculates, and spectrum energy analysis is carried out to calculated result.When nature of trouble difference, the analysis result of marginal spectrum energy is different, root According to different as a result, can be with accurate judgement nature of trouble.
A kind of permanent fault recognition methods based on fault electric arc characteristic the following steps are included:
Step A: it after system generation singlephase earth fault, circuit breaker trip, extracts transmission line of electricity both ends measurement point (m, n) The three-phase voltage and electric current at place obtain the modulus voltage and modulus electric current of route both ends measurement point through phase-model transformation, using gram Clarke transformation, if phase mode transformation matrix are T, transformation matrix is as follows
Voltage, current transformation process are as follows:
In formula, the three-phase voltage and electric current at the both ends mn are respectively umABC、unABC、imABC, inABC, the both ends mn modulus voltage and Modulus electric current is respectively um123、un123、im123、in123
Step B:
The value of fault point arc current feature expression is calculated using Two-Terminal Electrical Quantities.
Rectify the calculation formula to voltage traveling wave amount and the end n backward voltage traveling wave amount using m are as follows:
[um(t)+zc(t)*im(t)]*amn(t)
-[un(t)-zc(t)*in(t)]
=zc(t)*ik(t)*ank(t)
Rectify the calculation formula to voltage traveling wave amount using the end m backward voltage traveling wave amount and n are as follows:
[um(t)-zc(t)*im(t)]
-[un(t)+zc(t)*in(t)]*amn(t)
=zc(t)*ik(t)*amk(t)
In formula, wherein amkIt (t) is the time-domain expression of the propagation constant at fault point k to m side bus, ankIt (t) is failure The time-domain expression of propagation constant at point k to n side bus, amnIt (t) is propagation constant time domain of the m side bus at n side bus Expression formula.zcIt (t) is the time-domain expression of surge impedance of a line, um(t) the modulus voltage to be obtained according to the end route m phase voltage value Value, un(t) the modulus voltage value to be obtained according to the end route n phase voltage value.ik(t) be fault point arc current value.
U obtained in selecting step Am1、un1、im1、in1(mould) value substitutes into calculation formula, obtains the meter of feature expression Calculation value ux(t)。
Step C: to ux(t) EMD decomposition is carried out, and finds out the marginal spectrum energy value of each frequency component.
Step D: the ratio between high frequency spectrum energy and fundamental frequency spectrum energy P are found outR, with setting valve PsetCompare, works as PRLess than PsetWhen, It is determined as permanent fault, otherwise is transient fault.
Invention is further explained in the following with reference to the drawings and specific embodiments.
Step A: in extra high voltage network system shown in Fig. 1, after circuit breaker trip, transmission line of electricity both ends is extracted and are surveyed The three-phase voltage and electric current at point (m, n) are measured, obtains the modulus voltage and modulus electricity of route both ends measurement point through phase-model transformation Stream, using Clarke transform, transformation matrix is as follows
Voltage, current transformation process are as follows:
In formula, the three-phase voltage and electric current at the both ends mn are respectively umABC、unABC、imABC, inABC, the both ends mn modulus voltage and Modulus electric current is respectively um123、un123、im123、in123
Step B: calculating the value of faulty channel current characteristic expression formula using both-end traveling wave amount, and calculating process is as follows:
1. calculating characteristic impedance ZcWith the frequency-domain function of propagation coefficient A
Z '=R '+j ω L '
In formula, R is the resistance of transmission line of electricity unit length, and L ' is transmission line of electricity unit length inductance.
2. frequency domain value is transformed to time-domain value (rational function fitting process) using rational function fitting process:
To the Z of frequency domainc(ω) and A (ω) function is fitted, and fitting function scope discontinuity matching line segment slope is enabled to have The point of variation is the pole and zero of rational function, then
In formula, s=j ω and m > n;zn(n=1,2,3 ...) is zero point, pnFor pole, and all zero points are with pole Negative real number;K is coefficient;τminIndicate the shortest time that wave is propagated on the transmission line;Zc,approx(s) it is characterized having for impedance Manage approximation to function formula;AapproxIt (s) is the rational-function approximation of propagation coefficient.
Further obtain:
In formula, zc(t) it is characterized the time-domain value of impedance;A (t) is the time-domain value of propagation coefficient;knFor rational function expansion Coefficient;
3. calculating the value of fault point arc current feature expression
Such as Fig. 2, if transmission line of electricity overall length is l, fault point k and the end m bus distance are lmk, it is l with the end n bus distancenk。 There is following relationship to the end m direct wave:
uk(t)+zc(t)*imk(t)=[um(t)+zc(t)*im(t)]*amk(t) (2)
Wherein amkIt (t) is the time-domain expression of the propagation constant at fault point k to m side bus.zcIt (t) is surge impedance of a line Time-domain expression.
Similarly, there is following relationship to the end n backward-travelling wave:
uk(t)-zc(t)*ink(t)=[un(t)-zc(t)*in(t)]/*ank(t) (3)
Wherein ankIt (t) is the time-domain expression of the propagation constant at fault point k to n side bus ,/* indicates deconvolution.
To the arc current i of fault pointk(t), have
ik(t)=imk(t)+ink(t) (4)
Simultaneous (2) (3) (4) obtains
[um(t)+zc(t)*im(t)]*amk(t)
-[un(t)-zc(t)*in(t)]/*ank(t)
=zc(t)*ik(t) (5)
By (5) formula equal sign two sides while convolution ank(t) it obtains:
[um(t)+zc(t)*im(t)]*amn(t)
-[un(t)-zc(t)*in(t)]
=zc(t)*ik(t)*ank(t) (6)
Similarly, rectified using the end m backward voltage traveling wave and n to voltage traveling wave, can be obtained:
[um(t)-zc(t)*im(t)]
-[un(t)+zc(t)*in(t)]*amn(t)
=zc(t)*ik(t)*amk(t) (7)
Wherein amkIt (t) is the time-domain expression of the propagation constant at fault point k to m side bus, ank(t) it is arrived for fault point k The time-domain expression of propagation constant at n side bus, amnIt (t) is propagation constant time domain expression of the m side bus at n side bus Formula.zcIt (t) is the time-domain expression of surge impedance of a line.
By u obtained in step Am1、un1、im1、in1Value substitutes into (6) (7) formula, obtains the calculated value u of feature expressionx (t)。
Step C: to ux(t) EMD decomposition is carried out, and finds out the marginal spectrum energy value of each frequency component;
EMD decomposable process is as follows:
(1) signal u is found outx(t) all maximum point and minimum point use cubic spline letter to obtained extreme point Number is fitted to the upper and lower envelope of former data sequence.The mean value for calculating envelope up and down, is denoted as m1(t), former data sequence ux(t) Subtract m1(t) the new data sequence h for removing low frequency component can be obtained1(t).I.e.
ux(t)-m1(t)=h1(t) (10)
In general h1(t) it is unsatisfactory for the condition of IMF sequence, needs to repeat it above-mentioned treatment process k times thus, Until h1(t) the definition requirement for meeting IMF drives off the 1st IMF component c in this way1(t), it is ux(t) highest frequency in Component.
(2) by c1(t) from ux(t) it is separated in get the difference signal r of high fdrequency component is removed to one1(t), that is, have
r1(t)=ux(t)-c1(t) (11)
By r1(t) it is used as initial data, repeats step (1) (2), obtains second IMF component c2(t), n times are repeated, until rn(t) become monotonic function or amplitude is small to stopping decomposing when can ignore, obtain n IMF component and a residual components rn (t), it may be assumed that
Wherein, rnIt (t) is survival function, the average tendency of representation signal.And each IMF component c1(t)~cn(t) respectively The ingredient of different frequency sections from high to low is contained, after carrying out Hilbert transformation to each IMF component, is obtained
Wherein ai(t) c is indicatedi(t) amplitude, ωi(t) it indicates instantaneous frequency, can further be defined with Hilbert spectrum Marginal spectrum is
Sample frequency is 10000Hz in the present invention, and highest frequency is 5000Hz in marginal spectrum.If below point of frequency 44Hz Amount is low-frequency range, and 45Hz~55Hz is fundamental frequency section, and 56Hz~5000Hz is high band.Transient fault with two kinds of permanent fault In the case of marginal spectrum it is as shown in Figure 3, Figure 4.The energy value of its respective frequencies is indicated using marginal spectral amplitude ratio, if frequency is i's The energy value of component is Pi, low frequency energy, fundamental frequency energy, high-frequency energy is respectively Pd,PjAnd Pg, calculation formula is such as shown in (15):
Step D: finding out the ratio between high frequency spectrum energy and fundamental frequency spectrum energy, with setting valve PsetCompare:
According to P obtained in step CjAnd PgCalculate the ratio between high frequency spectrum energy and fundamental frequency spectrum energy, calculation formula PR= Pg/Pj
When transient fault occurs, arc characteristic is more apparent, comprising a large amount of high in arc current after fault phase trip Frequency component, ux(t) the close even more than fundamental frequency limit spectrum energy of the high frequency limit spectrum energy for including in.When generation permanent fault When, electric arc extinguishes rapidly, and the major frequency components of arc current are still fundamental frequency cycles signal in Grounding after fault phase trip, ux(t) high-frequency energy contained in can be ignored, at this time PRTheoretical value be zero.In view of the complexity of actual condition, answer There are enough allowances, therefore set PsetIt is 0.5.As the P being calculatedRWhen less than 0.5, failure be it is permanent, work as PRGreater than 0.5 When, failure is instantaneity.

Claims (5)

1. one kind is based on the permanent fault recognition methods of arc characteristic ultra-high/extra-high voltage transmission line of electricity, characterized in that steps are as follows:
Step A: it after system generation singlephase earth fault, circuit breaker trip, extracts at transmission line of electricity both ends measurement point (m, n) Three-phase voltage and electric current obtain the modulus voltage and modulus electric current of route both ends measurement point through phase-model transformation, using Clarke Transformation, phase mode transformation matrix T are specific as follows:
Voltage, current transformation process are as follows:
In formula, the three-phase voltage and electric current at the both ends mn are respectively umABC、unABC、imABC, inABC, the both ends mn modulus voltage and modulus are electric Stream is respectively um123、un123、im123、in123
Step B: the value of fault point arc current feature expression is calculated using Two-Terminal Electrical Quantities
Rectify the calculation formula to voltage traveling wave amount and the end n backward voltage traveling wave amount using m are as follows:
[um(t)+zc(t)*im(t)]*amn(t)
-[un(t)-zc(t)*in(t)]
=zc(t)*ik(t)*ank(t)
Rectify the calculation formula to voltage traveling wave amount using the end m backward voltage traveling wave amount and n are as follows:
[um(t)-zc(t)*im(t)]
-[un(t)+zc(t)*in(t)]*amn(t)
=zc(t)*ik(t)*amk(t)
In formula, wherein amkIt (t) is the time-domain expression of the propagation constant at fault point k to m side bus, ank(t) it is arrived for fault point k The time-domain expression of propagation constant at n side bus, amnIt (t) is propagation constant time domain expression of the m side bus at n side bus Formula.zcIt (t) is the time-domain expression of surge impedance of a line, umIt (t) is the modulus voltage value obtained according to the end route m phase voltage value, un (t) the modulus voltage value to be obtained according to the end route n phase voltage value.ik(t) be fault point arc current value.
U obtained in selecting step Am1、un1、im1、in1Ground modulus value substitutes into calculation formula, obtains the calculated value u of feature expressionx (t);
Step C: to ux(t) EMD decomposition is carried out, and finds out the marginal spectrum energy value of each frequency component.
Step D: the ratio between high frequency spectrum energy and fundamental frequency spectrum energy P are found outR, with setting valve PsetCompare, works as PRLess than PsetWhen, determine For permanent fault, on the contrary is transient fault.
2. being based on the permanent fault recognition methods of arc characteristic ultra-high/extra-high voltage transmission line of electricity, feature as described in claim 1 It is, further, step C: to ux(t) EMD decomposition is carried out, and finds out the marginal spectrum energy value of each frequency component;EMD was decomposed Journey is as follows:
(1) signal u is found outx(t) all maximum point and minimum point, it is quasi- using cubic spline function to obtained extreme point The upper and lower envelope for synthesizing former data sequence, calculates the mean value of envelope up and down, is denoted as m1(t), former data sequence ux(t) it subtracts m1(t) the new data sequence h for removing low frequency component is obtained1(t), i.e.,
ux(t)-m1(t)=h1(t) (10)
Repeat above-mentioned treatment process k times, until h1(t) the definition requirement for meeting IMF drives off the 1st IMF points in this way Measure c1(t), c1It (t) is ux(t) component of highest frequency in;
(2) by c1(t) from ux(t) it is separated in get the difference signal r of high fdrequency component is removed to one1(t), that is, have
r1(t)=ux(t)-c1(t) (11)
By r1(t) it is used as initial data, repeats step (1) (2), obtains second IMF component c2(t), n times are repeated, until rn (t) become monotonic function or amplitude is small to stopping decomposing when can ignore, obtain n IMF component and a residual components rn (t), it may be assumed that
Wherein, rnIt (t) is survival function, the average tendency of representation signal, and each IMF component c1(t)~cn(t) it separately includes The ingredient of different frequency sections from high to low obtains after carrying out Hilbert transformation to each IMF component
Wherein ai(t) c is indicatedi(t) amplitude, ωi(t) it indicates instantaneous frequency, can further define limit with Hilbert spectrum Spectrum is
3. being based on the permanent fault recognition methods of arc characteristic ultra-high/extra-high voltage transmission line of electricity, feature as described in claim 1 It is that in an example, sample frequency 10000Hz, highest frequency is 5000Hz in marginal spectrum, if below point of frequency 44Hz Amount is low-frequency range, and 45Hz~55Hz is fundamental frequency section, and 56Hz-5000Hz is high band, indicates that it corresponds to frequency using marginal spectral amplitude ratio The energy value of rate, if the energy value for the component that frequency is i is Pi, low frequency energy, fundamental frequency energy, high-frequency energy is respectively Pd,PjWith Pg, calculation formula is such as shown in (15):
4. being based on the permanent fault recognition methods of arc characteristic ultra-high/extra-high voltage transmission line of electricity, feature as claimed in claim 3 It is that further, step D: finding out the ratio between high frequency spectrum energy and fundamental frequency spectrum energy, with setting valve PsetCompare:
According to P obtained in step CjAnd PgCalculate the ratio between high frequency spectrum energy and fundamental frequency spectrum energy, calculation formula PR=Pg/ Pj
5. being based on the permanent fault recognition methods of arc characteristic ultra-high/extra-high voltage transmission line of electricity, feature as claimed in claim 3 It is, in an example, PsetIt is 0.5.As the P being calculatedRWhen less than 0.5, failure be it is permanent, work as PRWhen greater than 0.5, Failure is instantaneity.
CN201810757496.5A 2018-07-11 2018-07-11 Method for identifying permanent fault of ultra/ultra-high voltage transmission line based on arc characteristic Active CN109188189B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810757496.5A CN109188189B (en) 2018-07-11 2018-07-11 Method for identifying permanent fault of ultra/ultra-high voltage transmission line based on arc characteristic

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810757496.5A CN109188189B (en) 2018-07-11 2018-07-11 Method for identifying permanent fault of ultra/ultra-high voltage transmission line based on arc characteristic

Publications (2)

Publication Number Publication Date
CN109188189A true CN109188189A (en) 2019-01-11
CN109188189B CN109188189B (en) 2020-04-28

Family

ID=64935984

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810757496.5A Active CN109188189B (en) 2018-07-11 2018-07-11 Method for identifying permanent fault of ultra/ultra-high voltage transmission line based on arc characteristic

Country Status (1)

Country Link
CN (1) CN109188189B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112067944A (en) * 2020-09-09 2020-12-11 国网江苏省电力有限公司盐城供电分公司 Single-phase grounding fault detection control method for power transmission and distribution line
CN113608084A (en) * 2021-08-12 2021-11-05 山东理工大学 Series fault arc detection method

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0894698A (en) * 1994-09-27 1996-04-12 Tohoku Electric Power Co Inc Method and system for standardizing intermittent arc ground fault section in high voltage power distribution system with non-grounded neutral point
US20050134285A1 (en) * 2003-12-22 2005-06-23 Ross Mernyk Arc fault detector
CN101242095A (en) * 2008-02-03 2008-08-13 西安西瑞保护控制设备有限责任公司 A permanent failure judgement method for power transmission line with parallel reactor
CN102062831A (en) * 2010-10-29 2011-05-18 昆明理工大学 Single-phase permanent fault recognition method for extra-high voltage AC transmission line
CN102135588A (en) * 2011-02-18 2011-07-27 昆明理工大学 Intelligent ranging method for radiation network fault by utilizing S transformation spectrum characteristics
CN103116116A (en) * 2013-02-25 2013-05-22 浙江省电力公司电力科学研究院 Judging method and device of over line stoppage property of same pole double-circuit line with paralleling reactor
CN103529332A (en) * 2013-10-28 2014-01-22 昆明理工大学 Ultra-high voltage direct current transmission line lightning stroke interference recognition method based on voltage relevancy and wavelet transformation transient state energy distribution characteristics
CN103809078A (en) * 2014-01-26 2014-05-21 浙江大学 Single-phase permanent fault identification method for electric transmission line with paralleling reactors at two ends of electric transmission line
US8866487B2 (en) * 2010-09-10 2014-10-21 Southern States, Llc Directional fault sectionalizing system
CN107219440A (en) * 2017-04-11 2017-09-29 中国矿业大学 The localization method of single-ended radiation type distribution network singlephase earth fault

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0894698A (en) * 1994-09-27 1996-04-12 Tohoku Electric Power Co Inc Method and system for standardizing intermittent arc ground fault section in high voltage power distribution system with non-grounded neutral point
US20050134285A1 (en) * 2003-12-22 2005-06-23 Ross Mernyk Arc fault detector
CN101242095A (en) * 2008-02-03 2008-08-13 西安西瑞保护控制设备有限责任公司 A permanent failure judgement method for power transmission line with parallel reactor
US8866487B2 (en) * 2010-09-10 2014-10-21 Southern States, Llc Directional fault sectionalizing system
CN102062831A (en) * 2010-10-29 2011-05-18 昆明理工大学 Single-phase permanent fault recognition method for extra-high voltage AC transmission line
CN102135588A (en) * 2011-02-18 2011-07-27 昆明理工大学 Intelligent ranging method for radiation network fault by utilizing S transformation spectrum characteristics
CN103116116A (en) * 2013-02-25 2013-05-22 浙江省电力公司电力科学研究院 Judging method and device of over line stoppage property of same pole double-circuit line with paralleling reactor
CN103529332A (en) * 2013-10-28 2014-01-22 昆明理工大学 Ultra-high voltage direct current transmission line lightning stroke interference recognition method based on voltage relevancy and wavelet transformation transient state energy distribution characteristics
CN103809078A (en) * 2014-01-26 2014-05-21 浙江大学 Single-phase permanent fault identification method for electric transmission line with paralleling reactors at two ends of electric transmission line
CN107219440A (en) * 2017-04-11 2017-09-29 中国矿业大学 The localization method of single-ended radiation type distribution network singlephase earth fault

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
张云柯等: "带并联电抗器的超高压电缆-架空混合线路三相永久性故障识别方法", 《电力自动化设备》 *
李博通等: "带并联电抗器的超高压输电线路自适应重合闸新判据", 《电力***自动化》 *
钟佑明等: "希尔伯特-黄变换中边际谱的研究", 《***工程与电子技术》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112067944A (en) * 2020-09-09 2020-12-11 国网江苏省电力有限公司盐城供电分公司 Single-phase grounding fault detection control method for power transmission and distribution line
CN113608084A (en) * 2021-08-12 2021-11-05 山东理工大学 Series fault arc detection method

Also Published As

Publication number Publication date
CN109188189B (en) 2020-04-28

Similar Documents

Publication Publication Date Title
RU2489724C1 (en) Method of identifying type of short-circuiting in electric power line
CN109283428B (en) Zero-sequence component high-order differential transformation-based single-phase grounding transient protection method at feeder outlet
CN103226176B (en) A kind of wire selection method for power distribution network single phase earthing failure
CN109923423B (en) Method, system and device for fault detection
CN109283430B (en) Power distribution network fault location method based on voltage distribution principle
CN107543998B (en) Direct-current side fault positioning system and method for multi-terminal flexible direct-current power transmission system
CN109888744B (en) Protection method for high-voltage direct-current transmission line
CN104237731A (en) Single-phase earth fault line selection method on basis of EEMD (ensemble empirical mode decomposition) and energy method for resonance earthed power distribution network
CN103323741A (en) Method for judging D-shaped cable hybrid circuit fault section directing at strong fault and based on fault voltage initial row amplitude value comparison
CN109188189A (en) Based on arc characteristic ultra-high/extra-high voltage transmission line of electricity permanent fault recognition methods
Gajare et al. Synchrophasor-based intelligent autoreclosing scheme for series compensated transmission lines
CN110456219A (en) A kind of appraisal procedure of distribution line short voltage dip degree
CN109061397A (en) Line fault area identification method
CN105606955A (en) Numerical differentiation and empirical mode decomposition-based fault line distinguishing method
CN103311909B (en) Positive sequence Sudden Changing Rate and zero-sequence component is utilized to realize line single phase grounding failure voltage protection method
CN102435896A (en) Intermittent grounding fault rapid identification method of ship medium-voltage power system
CN110244181A (en) Marine wind electric field collection electric line sea cable electric fault localization method and its device
CN102035178B (en) Method for transformer excitation surge current latching
CN105811389B (en) The high anti-neutral point MOA of 750kV circuits and dielectric level determination method for parameter
CN104698312A (en) System and method for electrically detecting DC arrester of DC converter station
CN107179476A (en) A kind of Distribution Network Failure distance-finding method
CN112649694B (en) Method for judging single-phase grounding fault of low-current grounding system
Elangovan Fault noise based approach to phase selection using wavelets based feature extraction
CN109188191B (en) Rapid fault phase selection method for half-wavelength alternating-current transmission line based on power mutation
CN109245057B (en) Time domain full-waveform protection device for power transmission line and interphase mutation direction judgment method

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