CN109375065A - Travelling wave identification method and positioning device based on three-dimensional grey absolute correlation degree - Google Patents
Travelling wave identification method and positioning device based on three-dimensional grey absolute correlation degree Download PDFInfo
- Publication number
- CN109375065A CN109375065A CN201811516653.XA CN201811516653A CN109375065A CN 109375065 A CN109375065 A CN 109375065A CN 201811516653 A CN201811516653 A CN 201811516653A CN 109375065 A CN109375065 A CN 109375065A
- Authority
- CN
- China
- Prior art keywords
- frequency
- traveling wave
- wave
- pulverised
- curved surface
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Locating Faults (AREA)
Abstract
The invention discloses a kind of travelling wave identification method based on three-dimensional grey absolute correlation degree, it the steps include: the travelling wave signal for obtaining physical fault, draw out for the first when-frequency traveling wave Full wave shape；Obtain the second when-frequency traveling wave Full wave shape；Calculate the three-dimensional grey absolute correlation degree of the first when-frequency traveling wave Full wave shape and the second when-frequency traveling wave Full wave shape；According to three-dimensional grey absolute correlation degree, it identifies and the highest second when-frequency traveling wave Full wave shape of the first when-frequency traveling wave Full wave shape degree of association, by the quantitative analysis of different faults point traveling wave Full wave shape similitude and otherness, effective identification of different faults point traveling wave Full wave shape is realized.The present invention also provides the traveling wave positioning devices based on three-dimensional grey absolute correlation degree, including acquisition module, acquisition module, calculation of relationship degree module and identification module, realize the accurate positioning of physical fault.Methods and apparatus of the present invention is not limited by factors such as sample rate, transition resistances, and the practical application for traveling wave positioning opens a kind of new approaches.
Description
Technical field
Present invention relates generally to Relay Protection Technology in Power System fields, are related to one kind and are absolutely associated with based on three-dimensional grey
The travelling wave identification method and positioning device of degree.
Background technique
Fault traveling wave protection and localization method positioning accuracy theoretically with higher and reliability, obtain in power transmission network
It is widely applied.Fault traveling wave is a broadband step signal, it includes the time and frequency zones such as T/F-amplitude-polarity
Traveling wave panoramic information can characterize the abundant fault message including network topology structure, position of failure point, fault point parameter etc..Mesh
Before, fault traveling wave protection and localization method are broadly divided into following two: based on frequency domain information and based on time-domain information；Based on frequency domain
The fault traveling wave protection of information and localization method are based only upon certain two frequency range travelling wave signals (50~100kHz high frequency and 1~10kHz
Low frequency signal), for long distance transmission line end or opposite end bus-bar fault, since high frequency traveling wave signal decaying is serious, pole
Protection erroneous judgement even failure is easily led to, there are dead zones with localization method for the fault traveling wave protection based on frequency domain information；Based on time domain
Information fault traveling wave protection and localization method need the initial wavefront of accurate recognition and its subsequent catadioptric wave head polarity,
It is amplitude or arrival time, very high to sample rate requirement, 1MHz is generally taken, it is even higher, when failure initial phase angle is smaller (close to 0 °)
Or when high resistive fault, fault traveling wave wave head weak output signal, it is difficult to accurately identify.It can be seen that fault traveling wave in the prior art
Protection and influence of the localization method vulnerable to factors such as sample frequency, failure initial phase angle, transition resistances, depend critically upon wavefront
Accurately identify, had some limitations in practical application.
Summary of the invention
The present invention is directed to propose a kind of travelling wave identification method and positioning device based on three-dimensional grey absolute correlation degree, with solution
It is certainly influenced present in fault traveling wave in the prior art protection and locating scheme vulnerable to factors such as sample frequencys and applying has
The technical issues of limitation.
In order to solve the above technical problems, the invention adopts the following technical scheme:
The first aspect of the embodiment of the present invention provides a kind of traveling wave identification side based on three-dimensional grey absolute correlation degree
Method, comprising steps of
The initial traveling wave for obtaining currently practical failure reaches travelling wave signal after test point in preset duration, draws out current
The corresponding first when-frequency traveling wave Full wave shape of physical fault；
Obtain the second when-total evidence of frequency traveling wave Full wave shape；
The three-dimensional grey for successively calculating the described first when-frequency traveling wave Full wave shape and the second when-frequency traveling wave Full wave shape is exhausted
To the degree of association；
According to the three-dimensional grey absolute correlation degree, identify and the described first when-frequency traveling wave Full wave shape degree of association highest
The second when-frequency traveling wave Full wave shape.
Optionally, the total evidence of the second when-frequency traveling wave Full wave shape is obtained, including constructs fault waveform database in advance, from institute
It states and obtains the described second when-total evidence of frequency traveling wave Full wave shape in fault waveform database；
Building fault waveform database in advance, comprising:
The initial travelling wave signal for collecting history physical fault and/or simulated failure in advance, draws out corresponding fault traveling wave
When-frequency waveform, as the second when-frequency traveling wave Full wave shape；
When constructing the second of the described second when-frequency traveling wave Full wave shape-frequency can spectrum matrix；
When using described second-frequency can spectrum matrix as the second behavioural matrix, according to the element in second behavioural matrix
Initial line annhilation operation is calculated, the second initial line pulverised curved surface of second behavioural matrix is obtained；
Fault waveform database is constructed, when for storing described second-frequency traveling wave Full wave shape and the second initial line pulverised
Curved surface.
Optionally, the three-dimensional grey of the described first when-frequency traveling wave Full wave shape and the second when-frequency traveling wave Full wave shape is calculated
Absolute Correlation Analysis, comprising:
When constructing the first of the described first when-frequency traveling wave Full wave shape-frequency can spectrum matrix；
When using described first-frequency can spectrum matrix as the first behavioural matrix, according to the element in first behavioural matrix
Initial line annhilation operation is calculated, the first initial line pulverised curved surface of first behavioural matrix is obtained；
The second initial line pulverised curved surface is obtained from the fault waveform database；
Calculate the first initial line pulverised curved surface and the second initial line pulverised curved surface three-dimensional grey absolute correlation degree.
Optionally, the corresponding first when-frequency traveling wave Full wave shape of physical fault is drawn out, comprising:
Phase-model transformation is carried out to the travelling wave signal, extracts Aerial mode component, by continuous wavelet transform, is obtained containing sometimes
Between, the first when-frequency traveling wave Full wave shape including frequency, amplitude and polarity information.
Optionally, construct the described first when-frequency traveling wave Full wave shape first when-frequency energy spectrum matrix, comprising:
Described first when-frequency traveling wave Full wave shape is decomposed into M frequency range on frequency domain, and to continuous small obtained by each frequency range
Wave system number is divided into N number of period, obtains M × N number of time-frequency fritter, calculates the energy e of each time-frequency fritterji, when obtaining first-frequency
It can spectrum matrix EM×N。
Optionally, construct the described first when-frequency traveling wave Full wave shape first when-frequency energy spectrum matrix, comprising:
Described first when-frequency traveling wave Full wave shape is decomposed into M frequency range on frequency domain, and to continuous small obtained by each frequency range
Wave system number is divided into N number of period, obtains M × N number of time-frequency fritter；
Calculate the energy e of each time-frequency fritterji, frequency range j is now defined in the energy of the time-frequency fritter of period i are as follows:
In formula, i1And ipStarting and end sampled point for i-th of period, cj(k) in frequency band j, sampled point is corresponding
Continuous wavelet coefficient；Wherein, 1≤j≤M, 1≤i≤N；
Successively calculate the energy e for each the time-frequency fritter being divided intoji, when can obtain first-frequency can spectrum matrix EM×N:
When optionally, using described first-frequency can spectrum matrix as the first behavioural matrix, according in first behavioural matrix
Element calculate initial line annhilation operation, obtain the initial line pulverised curved surface of first behavioural matrix, comprising:
When by described first-frequency can spectrum matrix EM×NAs the first behavioural matrix, according to the member in first behavioural matrix
Element calculates initial line annhilation operation, obtains the initial line pulverised curved surface of first behavioural matrix:
EM×ND=(ejid)M×N=eji-ej1
Wherein D is known as the initial line annhilation operation of behavioural matrix, EM×ND is known as EM×NInitial line pulverised curved surface expression matrix
Form.
Optionally, the first initial line pulverised curved surface is calculated absolutely to be associated with the second initial line pulverised curved surface three-dimensional grey
Degree includes:
Calculate first behavioural matrix and the corresponding first initial line pulverised curved surface of second behavioural matrix and
The three-dimensional grey absolute correlation degree ε of two initial line pulverised curved surfacespq:
In formula,
Wherein, XpFor the first behavioural matrix, XqFor the second behavioural matrix,Indicate the first initial line pulverised curved surface,It indicates
Second initial line pulverised curved surface, | sp| and | sq| be the first initial line pulverised curved surface and the second initial line pulverised curved surface respectively with coordinate plane
The cylinder volume surrounded, | sp-sq| it is the cylinder that the first initial line pulverised curved surface and the second initial line pulverised curved surface surround
Volume,WithIt is the first initial line pulverised curved surface respectivelyWith the second initial line pulverised curved surfaceAt two-dimensional space midpoint (i, j)
The behavior value at place, WithIndicate meaning to should refer toWithεpqIt takes
Value is [0,1], and i indicates that the i-th column in behavioural matrix, j indicate jth row in behavioural matrix, 1≤j≤M, 1≤i≤N.
Optionally, the initial traveling wave for obtaining physical fault reaches travelling wave signal after test point in preset duration, comprising:
Travelling wave signal after the initial traveling wave arrival test point of acquisition physical fault in 1ms time window respectively.
The second aspect of the embodiment of the present invention provides a kind of traveling wave positioning dress based on three-dimensional grey absolute correlation degree
It sets, including acquisition module, acquisition module, calculation of relationship degree module and identification module；
Acquisition module, the initial traveling wave for obtaining currently practical failure reach the traveling wave letter after test point in preset duration
Number, draw out the corresponding first when-frequency traveling wave Full wave shape of currently practical failure；
Module is obtained, for obtaining the total evidence of the second when-frequency traveling wave Full wave shape；
Calculation of relationship degree module, when for calculating described first-frequency traveling wave Full wave shape and the second when-frequency traveling wave all-wave
The three-dimensional grey absolute correlation degree of shape；
Identification module, for identifying and the described first when-frequency traveling wave all-wave according to the three-dimensional grey absolute correlation degree
It is fixed according to the failure of the described second when-frequency traveling wave Full wave shape position to occur for the highest second when-frequency traveling wave Full wave shape of the shape degree of association
The generation position of the currently practical failure in position.
Compared with prior art, the present invention can reach following technical effect:
Travelling wave identification method proposed by the present invention based on three-dimensional grey absolute correlation degree obtains after physical fault generation
The physical fault travelling wave signal for taking preset duration, draws out the corresponding first when-frequency traveling wave Full wave shape of physical fault, and with second
When-frequency traveling wave Full wave shape does three-dimensional grey absolute correlation degree and calculates, the second when-frequency traveling wave Full wave shape is pre-set multiple moulds
The corresponding Full wave shape of travelling wave signal that quasi- fault point generates, or the physical fault travelling wave signal that history occurs is corresponding complete
Waveform.It is calculated by the three-dimensional grey absolute correlation degree, obtains being associated with angle value accordingly, and then obtain and currently practical failure
When first-the highest simulated failure waveform of the frequency traveling wave Full wave shape degree of association or history physical fault waveform, it realizes to current reality
The identification of border failure Full wave shape；And the generation position of simulated failure and history physical fault can be pre-recorded, currently practical failure
Specific location on the line is then determined as the hair of the highest simulated failure of the Full wave shape degree of association or history physical fault therewith
Raw position or near, and then realize fault location.The recognition methods provided by the invention is based on time and frequency zone traveling wave Full wave shape,
By three-dimensional grey absolute correlation degree, the traveling wave Full wave shape in the case of different faults point is effectively identified, and to its similitude and difference
The opposite sex carries out quantitative analysis, realizes the identification of fault traveling wave；Also, importantly, the method for waveform identification is without detection
Initial wavefront arrival time, amplitude and polarity, without recognizing the 2nd backward-travelling wave, without detecting the transmission of different modulus
Velocity of wave is not limited by sample rate, especially in the Weak faults such as failure initial phase angle smaller (close to 0 °) or high resistive fault,
Waveform identification degree still with higher substantially increases the reliability of traveling-wave protection and positioning theoretically without waveform recognition dead zone.
Detailed description of the invention
It in order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, below will be to institute in embodiment
Attached drawing to be used is needed to be briefly described, it should be apparent that, the accompanying drawings in the following description is only some implementations of the invention
Example, for those of ordinary skill in the art, without any creative labor, can also be according to these attached drawings
Obtain other attached drawings.
Fig. 1 is that the present invention is based on the signals of the process of one embodiment of the travelling wave identification method of three-dimensional grey absolute correlation degree
Figure；
Fig. 2 is that the present invention is based on the system knots of a simulated example of the travelling wave identification method of three-dimensional grey absolute correlation degree
Structure schematic diagram；
Fig. 3 be one embodiment of the invention in when-frequency traveling wave Full wave shape schematic diagram.
Specific embodiment
The following is a clear and complete description of the technical scheme in the embodiments of the invention, it is clear that described embodiment
Only a part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, the common skill in this field
Art personnel every other embodiment obtained without making creative work belongs to the model that the present invention protects
It encloses.
To keep technical solution of the present invention more obvious and easy to understand, the present invention is made With reference to embodiment further
Detailed description.
Embodiment 1
The embodiment of the present invention 1 provides a kind of travelling wave identification method based on three-dimensional grey absolute correlation degree, referring to Fig. 1 institute
Show comprising step:
Step S101, the initial traveling wave for obtaining currently practical failure reach travelling wave signal after test point in preset duration,
Draw out the corresponding first when-frequency traveling wave Full wave shape of currently practical failure.
Since at the time of the initial traveling wave of failure reaches test point, the waveform of certain time window T is intercepted, continuous wavelet is utilized
Transformation can be drawn out: in time window T, when fault traveling wave-frequency waveform, define waveform when being-frequency traveling wave Full wave shape.In this hair
In bright embodiment, when-frequency traveling wave Full wave shape include the first when-frequency traveling wave Full wave shape and the second when-frequency traveling wave Full wave shape.It is current real
Failure corresponding Full wave shape in border is the first when-frequency traveling wave Full wave shape.
Step S102 obtains the total evidence of the second when-frequency traveling wave Full wave shape.
When second-frequency traveling wave Full wave shape is totally according to including the second when-frequency traveling wave Full wave shape and/or for calculating three-dimensional grey
The necessary data (such as when second-the corresponding initial line pulverised curved surface of frequency traveling wave Full wave shape) of Absolute Correlation Analysis.Wherein second when-frequency
Traveling wave Full wave shape is the Full wave shape of the travelling wave signal of simulated failure and the travelling wave signal of history physical fault.
Since at the time of the initial traveling wave of simulated failure or history physical fault reaches test point, certain time is intercepted
The waveform of window T can be drawn out using continuous wavelet transform: in time window T, when fault traveling wave-frequency waveform, which is
When two-frequency traveling wave Full wave shape.
When second-frequency traveling wave Full wave shape is totally according to further including specific position that simulated failure or history physical fault occur
It sets.
It specifically, can be in every line of the transmission line of electricity local topology network in protection, orientation range for simulated failure
A simulated failure point is arranged at interval of certain distance length (such as 0.5-1.5km) in road, and simulation generates failure and acquisition should
The initial traveling wave of simulated failure reaches the travelling wave signal in test point preset duration, draw out traveling wave Full wave shape and record to get
When corresponding to simulated failure second-frequency traveling wave Full wave shape；The initial traveling wave for obtaining each physical fault in history reaches detection
Travelling wave signal after point in preset duration is drawn out traveling wave Full wave shape and is recorded to get history physical fault corresponding second is arrived
When-frequency traveling wave Full wave shape.
Step S103, the three-dimensional grey for successively calculating the first when-frequency traveling wave Full wave shape and the second when-frequency traveling wave Full wave shape are exhausted
To the degree of association.
Step S104 is identified and the first when-frequency traveling wave Full wave shape degree of association highest according to three-dimensional grey absolute correlation degree
The second when-frequency traveling wave Full wave shape.
The embodiment of the present invention is directed to time and frequency zone traveling wave Full wave shape, proposes a kind of based on three-dimensional grey absolute correlation degree
Method for waveform identification, by interception certain time window in time and frequency zone fault traveling wave panorama waveform, can show network topology structure,
Abundant primary fault information including position of failure point, fault point parameter；Due to different faults point in area and internal fault external fault point,
Traveling wave Full wave shape has similitude and otherness, by three-dimensional grey absolute correlation degree to this similitude and the otherness amount of progress
Change analysis, realizes effective identification of different faults point traveling wave Full wave shape.The method for waveform identification is not limited by sample rate, is not necessarily to
The initial wavefront of accurate recognition and its subsequent reflection wave head realize reliably accurately identifying for Weak fault travelling wave signal, are
The practical application of protection and localization method based on traveling wave information opens a kind of new approaches.
Embodiment 2
One embodiment of travelling wave identification method of the offer of the embodiment of the present invention 2 based on three-dimensional grey absolute correlation degree.
In this embodiment, the method comprising the steps of:
S201 constructs fault waveform database in advance.
1: 1 true type traveling wave experiment porch is built, fault traveling wave detection is built in dynamic model experiment room and is surveyed with fault location dynamic model
Test system constructs power grid traveling-wave waveform database by simulation analysis, laboratory test and dynamic simulation test: including phase-to phase fault base
Quasi- database PP_faultWith ground fault benchmark database PG_fault；Simulated failure data are acquired at the scene, or collect historical failure
Data establish fault waveform database as complete and closing to reality as possible.That is, collecting history physical fault in advance
And/or the travelling wave signal of simulated failure, when drawing out the fault traveling wave of preset duration-frequency waveform, it is complete as the second when-frequency traveling wave
Waveform, building include the database of the total evidence of the second when-frequency traveling wave Full wave shape.In the present embodiment, the second when-frequency traveling wave all-wave
The definition of the total evidence of shape is repeated no more with embodiment 1.
Preferably, in this embodiment, it when constructing fault waveform database, also completes three-dimensional grey absolute correlation degree and calculates
In the second when-frequency traveling wave Full wave shape basic calculation, when including being previously-completed second-the second initial line pulverised of frequency traveling wave Full wave shape
The calculating of curved surface.
S202, building the second when-frequency traveling wave Full wave shape second when-frequency energy spectrum matrix.
S203, when using second-frequency can spectrum matrix as the second behavioural matrix, calculated according to the element in the second behavioural matrix
Initial line annhilation operation obtains the second initial line pulverised curved surface of the second behavioural matrix, stores to fault waveform database.
S204, the initial traveling wave for obtaining currently practical failure reach travelling wave signal after test point in 1ms time window, draw
It has the corresponding first when-frequency traveling wave Full wave shape of currently practical failure.
Preferably, phase-model transformation is carried out to travelling wave signal, extracts Aerial mode component and is contained by continuous wavelet transform
The first when-frequency traveling wave Full wave shape including time, frequency, amplitude and polarity information.When second-drafting of frequency traveling wave Full wave shape
Journey is the same as the first when-frequency traveling wave Full wave shape.
In view of traveling wave Full wave shape should be reflected with reflection process in the multiple of wave impedance discontinuity point (through scene comprising traveling wave
Test discovery, traveling wave by 5 times or more folding, after reflection, waveform attenuating is serious, so the present invention consider traveling wave surge folding,
Order of reflection is less than or equal to 5 times), traveling wave Full wave shape just contains network topology structure abundant and position of failure point information；As far as possible
The influence of perimeter circuit back wave is avoided, to reduce or even eliminate influence of the network topology structure variation to algorithm；When consideration
Between window and sample rate cooperation.Comprehensively consider factors above, in embodiments of the present invention, time window value 1ms.
S205, building the first when-frequency traveling wave Full wave shape first when-frequency energy spectrum matrix；
S206, when using first-frequency can spectrum matrix as the first behavioural matrix, calculated according to the element in the first behavioural matrix
Initial line annhilation operation obtains the first initial line pulverised curved surface of the first behavioural matrix；
S207 obtains the second initial line pulverised curved surface from fault waveform database；
S208 calculates the first initial line pulverised curved surface and the second initial line pulverised curved surface three-dimensional grey absolute correlation degree.
S209 is identified and the first when-frequency traveling wave Full wave shape degree of association highest according to three-dimensional grey absolute correlation degree
When two-frequency traveling wave Full wave shape.
It is the positioning result as currently practical failure that position, which occurs, for the corresponding failure of the highest travelling wave signal of the degree of association.
When constructing fault traveling wave database, while it being previously-completed the basic calculation of the second when-frequency traveling wave Full wave shape, it can be
After physical fault occurs, reduce the calculation amount that Trouble Match is carried out by three-dimensional grey absolute correlation degree, so reduce identification and
Positioning time, more rapidly, efficiently realize identify and position.
Embodiment 3
The embodiment of the present invention 3 provides an a kind of preferred reality of travelling wave identification method based on three-dimensional grey absolute correlation degree
Apply example.The preferred embodiment comprising steps of
A simulated failure point is arranged at interval of 1km on every route of transmission line of electricity network in advance, respectively more in S301
A simulated failure point simulated failure generates, the 1ms after the initial traveling wave that test point acquires each simulated failure reaches the test point
Travelling wave signal in time is drawn the time and frequency zone waveform of 1ms internal fault traveling wave using continuous wavelet transform, obtains multiple second
When-frequency traveling wave Full wave shape.
S302: when the initial traveling wave of physical fault reaches test point, when interception a length of 1ms traveling-wave waveform, the company of utilization
Continuous wavelet transformation draws the time and frequency zone waveform of 1ms internal fault traveling wave, i.e. the first when-frequency traveling wave Full wave shape.
Preferably, in the step, since the initial traveling wave of failure reaches test point, interception 1ms line line wave component passes through
Continuous wavelet transform is reflected the traveling wave Full wave shape of signal time and frequency zone feature simultaneously.Traveling wave wave initial for different faults point
Head, at the time of identical criterion should being used to judge that the wavefront reaches test point, until timing 1ms, intercepts wave in the time window
Shape information, using continuous wavelet transform obtain the waveform when-frequency waveform.Since traveling wave is multiple at wave impedance discontinuity point
Refraction and reflection process, traveling wave by 5 times or more folding, after reflection, decaying is serious, and the embodiment of the present invention will intercept duration value
1ms, to avoid the influence of surrounding network back wave as far as possible.
S303: when by above-mentioned 1ms the first-frequency traveling wave Full wave shape is decomposed into M frequency range on frequency domain, and to each frequency range institute
It obtains continuous wavelet coefficient and is divided into N number of period, obtain M × N number of time-frequency fritter, calculate the energy e of each time-frequency fritterji, obtain
When reflecting the first of original signal-frequency energy spectrum matrix EM×N。
Preferably as a kind of embodiment, step S303 be may particularly include:
The traveling wave Full wave shape signal that will be truncated is decomposed into M frequency range on frequency domain, and to continuous small obtained by each frequency range
Wave system number is divided into N number of period, obtains M × N number of time-frequency fritter；
Calculate the energy e of each time-frequency fritterji, frequency range j is now defined in the energy of the time-frequency fritter of period i are as follows:
In formula, iNAnd i1Starting and end sampled point for i-th of period, cj(k) in frequency band j, sampled point is corresponding
Continuous wavelet coefficient.
Successively calculate the energy e for each the time-frequency fritter being divided intojiAfterwards can be obtained first when-frequency energy spectrum matrix
EM×N:
S304: using for reference gray system theory, will above-mentioned first when-frequency energy spectrum matrix EM×NAs the first behavioural matrix, to this
First behavior matrix element calculates initial line annhilation operation, obtains the first initial line pulverised curved surface of first behavioural matrix, EM×ND=
(ejid)M×N=eji-ej1, wherein D is known as the initial line annhilation operation of behavioural matrix, EM×ND is EM×NInitial line pulverised curved surface square
Battle array expression-form, is denoted asIn the step of subsequent calculating three-dimensional grey absolute correlation degree, to retouch
It states conveniently, by the first initial line pulverised curved surface and the second initial line pulverised curved surface respectively with symbolWithIt indicates.
Calculate behavioural matrix initial line pulverised curved surface physical significance are as follows: by the cycling of elements in behavioural matrix be dimensionless
Corresponding behavior curved surface is reset to same initial line by data, performance in space.
Building simulated failure second when-frequency spectrum matrix and can calculate initial line annhilation operation and obtain the second initial line pulverised song
The step of face, with step S303 and 304, repeats no more.
S305: three-dimensional grey absolute correlation degree ε is calculatedpq, pass through εpqTo the phase for the traveling wave Full wave shape that different faults point generates
Quantitative analysis is carried out like property and otherness, to realize effective identification of different faults point traveling wave Full wave shape.Calculate three-dimensional grey
Absolute Correlation Analysis εpqFor
In formula,
The physical significance of three-dimensional grey absolute correlation degree are as follows: measure the difference of folded volume between two initial line pulverised curved surfaces
Degree.Wherein εpqValue is [0,1], and two behavioural matrixes correspond to that curve form is closer, then the degree of association is bigger, i.e. εpqMore connect
Nearly 1.
Wherein, XpFor the first behavioural matrix, XqFor the second behavioural matrix,Indicate the first initial line pulverised curved surface,It indicates
Second initial line pulverised curved surface, | sp| and | sq| be the first initial line pulverised curved surface and the second initial line pulverised curved surface respectively with coordinate plane
The cylinder volume surrounded, | sp-sq| it is the cylinder that the first initial line pulverised curved surface and the second initial line pulverised curved surface surround
Volume,WithIt is the first initial line pulverised curved surface respectivelyWith the second initial line pulverised curved surfaceAt two-dimensional space midpoint (i, j)
Behavior value, WithIndicate meaning to should refer toWithεpqValue is
[0,1], i indicate that the i-th column in behavioural matrix, j indicate jth row in behavioural matrix, 1≤j≤M, 1≤i≤N.
Specifically,Indicate the first initial line pulverised curved surfaceBehavior value at two-dimensional space point (i+1, j),
Indicate the first initial line pulverised curved surfaceBehavior value at two-dimensional space point (i, j+1),Indicate that the first initial line pulverised is bent
FaceBehavior value at two-dimensional space point (i+1, j+1)；Indicate the second initial line pulverised curved surfaceIn two-dimensional space point (i
+ 1, j) at behavior value,Indicate the second initial line pulverised curved surfaceBehavior value at two-dimensional space point (i, j+1),Indicate the second initial line pulverised curved surfaceBehavior value at two-dimensional space point (i+1, j+1).
Those skilled in the art will be appreciated that,WithOnly indicate that symbol should substitute into the first initial line zero when reality calculates
Change curved surface or the corresponding matrix E of the second initial line pulverised curved surface institute0 M×NIn element value accordingly calculated.
As shown in Fig. 2, being IEEE14 node standard test system figure, Tu2Zhong substation 4 is equipped with traveling wave acquisition device.
Based on test macro shown in Fig. 2, emulation experiment is carried out to travelling wave identification method of the invention.
Shown in Fig. 2 1-14 substation and its between transmission line of electricity.(the F in Fig. 2 assuming that singlephase earth fault1) hair
The raw head end in route 5-4 (route i.e. between substation 5 and substation 4), away from No. 5 substation 21.4km, failure initial phase angle
The singlephase earth fault for being 35 Ω for 90 °, transition resistance.And one simulated failure point (figure is set every 1km on every route
Simulated failure point is not shown in 2), wherein route 5-4 gives simulated failure piont mark since reference end 5, be denoted as Nn (n=1,2,
3 ...)；The physical fault travelling wave signal that No. 5 substation's test points detect is obtained, Aerial mode component is extracted by phase-model transformation,
Continuous wavelet transform is being utilized, traveling wave Full wave shape is being obtained, as shown in figure 3, calculating simulation fault point and actual fault point traveling wave are complete
The three-dimensional grey absolute correlation degree ε of waveformpq, wherein interception waveform time window is 1ms, sample rate is set as 0.2MHz, and will count
According to filling the following table 1.
Singlephase earth fault traveling wave Full wave shape recognition result occurs for 1 route 5-4 head end of table
Faulty line external analog fault point traveling wave Full wave shape recognition result
(note: non-fault line simulated failure point and the three-dimensional grey absolute correlation degree of actual fault point traveling wave Full wave shape take
The average value of this route each point)
Table 1 statistics indicate that: for different faults point in area, each secondary traveling wave surge transmission path is identical, reflects and reflected
Cheng Xiangtong, arrival timing is consistent, and traveling wave Full wave shape similarity is high, but due to the transmission distance of different location fault traveling wave to test point
From having differences, cause waveforms amplitude attenuation degree different, when fault point is opposite closer apart, similarity is higher, three-dimensional grey
Absolute Correlation Analysis numerical value more levels off to 1, and when fault point relative distance increases, similarity is reduced, three-dimensional grey absolute correlation degree number
Value reduces；For internal fault external fault, each secondary traveling wave surge transmission path is different, and refraction is different with reflection process, reaches timing not
Together, the traveling wave Full wave shape waveform difference opposite sex is very big, and three-dimensional grey absolute correlation degree reduces rapidly, close to 0.
It (sets fault point to verify different transition resistances, initial phase angle and sample frequency under same failure and single-phase connect occurs
Earth fault) traveling wave Full wave shape recognition result influence, be respectively established and emulated, the data obtained such as table 2, table 3,4 institute of table
Show:
The traveling wave Full wave shape recognition result of the different transition resistances of table 2
The fault traveling wave Full wave shape recognition result of the different initial phase angles of table 3
The fault traveling wave Full wave shape recognition result of the different sample rates of table 4
Table 2, table 3, table 4 statistics indicate that: three-dimensional grey absolute correlation degree traveling wave Full wave shape recognition methods proposed by the present invention
Substantially it is not influenced by transition resistance, failure initial phase angle and sample rate, there is practical engineering application value.
Embodiment 4
The embodiment of the present invention also provides a kind of traveling wave positioning device based on three-dimensional grey absolute correlation degree, including acquisition mould
Block, calculation of relationship degree module and identification module.
Acquisition module, the initial traveling wave for obtaining physical fault reach the travelling wave signal after test point in preset duration,
The corresponding first when-frequency traveling wave Full wave shape of physical fault is drawn out, and obtains the second when-frequency traveling wave Full wave shape.
Calculation of relationship degree module, for calculating the three-dimensional of the first when-frequency traveling wave Full wave shape and the second when-frequency traveling wave Full wave shape
Grey absolute correlation degree.
Identification module, for identifying and the first when-frequency traveling wave Full wave shape degree of association according to three-dimensional grey absolute correlation degree
When highest second-frequency traveling wave Full wave shape, in turn, the highest according to the degree of association second when-corresponding simulation of frequency traveling wave Full wave shape
Failure point determines actual fault point position.
Preferably, calculation of relationship degree module, when for constructing the first of the first when-frequency traveling wave Full wave shape-frequency energy spectrum matrix；
When using first-frequency can spectrum matrix as the first behavioural matrix, initial line annhilation operation is calculated according to the element in the first behavioural matrix,
Obtain the first initial line pulverised curved surface of the first behavioural matrix；And the second when-frequency energy for constructing the second when-frequency traveling wave Full wave shape
Spectrum matrix；When using second-frequency can spectrum matrix as the second behavioural matrix, initial line zero is calculated according to the element in the second behavioural matrix
Change operator, obtains the second initial line pulverised curved surface of the second behavioural matrix；The first initial line pulverised curved surface and described second are calculated later
Initial line pulverised curved surface three-dimensional grey absolute correlation degree.
Preferably, acquisition module extracts Aerial mode component, passes through continuous wavelet for carrying out phase-model transformation to travelling wave signal
Transformation, obtains containing the first when-frequency traveling wave Full wave shape including having time, frequency, amplitude and polarity information.
To sum up, the travelling wave identification method proposed by the present invention based on three-dimensional grey absolute correlation degree, by collected traveling wave
When signal draws traveling wave using continuous wavelet transform-frequency schemes, and when calculating-frequency energy spectrum matrix, row can be showed visual and clearly
Wave signal shows rich including network topology structure, position of failure point, fault point parameter in the time-frequency characteristics information of transmission process
Rich primary fault information carries out fault location according to the recognition methods, compensates in traditional traveling wave detector method and only analyze traveling wave
Time domain specification and can not judge the frequency distribution of traveling wave and the drawbacks of according to frequency variation characteristic, it is special also to compensate for only analysis traveling wave frequency domain
Property and the drawbacks of wavefront arrival time and timing can not be obtained；
Also, the travelling wave identification method proposed through the invention is, it can be achieved that the power grid based on time and frequency zone traveling wave Full wave shape is protected
Shield and fault location, are not limited by sample rate, are not necessarily to the initial wavefront of accurate recognition and its subsequent reflection wave head, theoretically
Without waveform recognition dead zone, for based on traveling wave information protection and localization method provide a kind of new technical concept, for reality
The functionization of existing fault traveling wave protection and positioning has important theory and realistic meaning.
The above is only the preferred embodiment of the present invention, protection scope of the present invention is not limited merely to above-described embodiment,
All technical solutions belonged under thinking of the present invention all belong to the scope of protection of the present invention.It should be pointed out that for the art
For those of ordinary skill, several improvements and modifications without departing from the principles of the present invention should be regarded as protection of the invention
Range.
Claims (10)
1. the travelling wave identification method based on three-dimensional grey absolute correlation degree, which is characterized in that comprising steps of
The initial traveling wave for obtaining currently practical failure reaches travelling wave signal after test point in preset duration, draws out currently practical
The corresponding first when-frequency traveling wave Full wave shape of failure；
Obtain the second when-total evidence of frequency traveling wave Full wave shape；
The three-dimensional grey for successively calculating the described first when-frequency traveling wave Full wave shape and the second when-frequency traveling wave Full wave shape is absolutely closed
Connection degree；
According to the three-dimensional grey absolute correlation degree, identify and the described first when-frequency traveling wave Full wave shape degree of association highest the
When two-frequency traveling wave Full wave shape.
2. the travelling wave identification method according to claim 1 based on three-dimensional grey absolute correlation degree, which is characterized in that described
Step obtains the total evidence of the second when-frequency traveling wave Full wave shape, including constructs fault waveform database in advance, from the fault wave figurate number
According to obtaining the described second when-total evidence of frequency traveling wave Full wave shape in library；
The preparatory building fault waveform database, comprising:
The initial travelling wave signal for collecting history physical fault and/or simulated failure in advance, when drawing out corresponding fault traveling wave-frequency
Waveform, as the second when-frequency traveling wave Full wave shape；
When constructing the second of the described second when-frequency traveling wave Full wave shape-frequency can spectrum matrix；
When using described second-frequency can spectrum matrix as the second behavioural matrix, calculated according to the element in second behavioural matrix
Initial line annhilation operation obtains the second initial line pulverised curved surface of second behavioural matrix；
Fault waveform database is constructed, when for storing described second-frequency traveling wave Full wave shape and the second initial line pulverised curved surface.
3. the travelling wave identification method according to claim 2 based on three-dimensional grey absolute correlation degree, which is characterized in that described
Step calculates the three-dimensional grey absolute correlation degree of the described first when-frequency traveling wave Full wave shape and the second when-frequency traveling wave Full wave shape,
Include:
When constructing the first of the described first when-frequency traveling wave Full wave shape-frequency can spectrum matrix；
When using described first-frequency can spectrum matrix as the first behavioural matrix, calculated according to the element in first behavioural matrix
Initial line annhilation operation obtains the first initial line pulverised curved surface of first behavioural matrix；
The second initial line pulverised curved surface is obtained from the fault waveform database；
Calculate the first initial line pulverised curved surface and the second initial line pulverised curved surface three-dimensional grey absolute correlation degree.
4. the travelling wave identification method according to claim 1-3 based on three-dimensional grey absolute correlation degree, feature
It is, the step draws out the corresponding first when-frequency traveling wave Full wave shape of physical fault, comprising:
Phase-model transformation is carried out to the travelling wave signal, extracts Aerial mode component, by continuous wavelet transform, is obtained containing having time, frequency
The first when-frequency traveling wave Full wave shape including rate, amplitude and polarity information.
5. the travelling wave identification method according to claim 4 based on three-dimensional grey absolute correlation degree, which is characterized in that described
When step constructs the first of the described first when-frequency traveling wave Full wave shape-frequency can spectrum matrix, comprising:
Described first when-frequency traveling wave Full wave shape is decomposed into M frequency range on frequency domain, and to continuous wavelet system obtained by each frequency range
Number is divided into N number of period, obtains M × N number of time-frequency fritter, calculates the energy e of each time-frequency fritterji, obtain the first when-frequency power spectrum
Matrix EM×N。
6. the travelling wave identification method according to claim 5 based on three-dimensional grey absolute correlation degree, which is characterized in that described
When step constructs the first of the described first when-frequency traveling wave Full wave shape-frequency can spectrum matrix, comprising:
Described first when-frequency traveling wave Full wave shape is decomposed into M frequency range on frequency domain, and to continuous wavelet system obtained by each frequency range
Number is divided into N number of period, obtains M × N number of time-frequency fritter；
Calculate the energy e of each time-frequency fritterji, frequency range j is now defined in the energy of the time-frequency fritter of period i are as follows:
In formula, i1And ipStarting and end sampled point for i-th of period, cj(k) in frequency band j, sampled point is corresponding continuous
Wavelet coefficient；Wherein, 1≤j≤M, 1≤i≤N；
Successively calculate the energy e for each the time-frequency fritter being divided intoji, when can obtain first-frequency can spectrum matrix EM×N:
7. the travelling wave identification method according to claim 6 based on three-dimensional grey absolute correlation degree, which is characterized in that described
When step is using described first-frequency can spectrum matrix calculated as the first behavioural matrix, according to the element in first behavioural matrix
Initial line annhilation operation obtains the initial line pulverised curved surface of first behavioural matrix, comprising:
When by described first-frequency can spectrum matrix EM×NAs the first behavioural matrix, according to the element meter in first behavioural matrix
Initial line annhilation operation is calculated, the initial line pulverised curved surface of first behavioural matrix is obtained:
EM×ND=(ejid)M×N=eji-ej1
E0 M×N=EM×ND=(e0 ji)M×N
Wherein D is known as the initial line annhilation operation of behavioural matrix, EM×ND is known as EM×NInitial line pulverised curved surface expression matrix form.
8. the travelling wave identification method according to claim 3 based on three-dimensional grey absolute correlation degree, which is characterized in that described
Step calculates the first initial line pulverised curved surface and the second initial line pulverised curved surface three-dimensional grey absolute correlation degree includes:
First behavioural matrix and the corresponding first initial line pulverised curved surface of second behavioural matrix and second is calculated to begin
The three-dimensional grey absolute correlation degree ε of side pulverised curved surfacepq:
In formula,
Wherein, XpFor the first behavioural matrix, XqFor the second behavioural matrix,Indicate the first initial line pulverised curved surface,Indicate second
Initial line pulverised curved surface, | sp| and | sq| it is that the first initial line pulverised curved surface and the second initial line pulverised curved surface are surrounded with coordinate plane respectively
Cylinder volume, | sp-sq| for the body for the cylinder that the first initial line pulverised curved surface and the second initial line pulverised curved surface surround
Product,WithIt is the first initial line pulverised curved surface respectivelyWith the second initial line pulverised curved surfaceAt two-dimensional space midpoint (i, j)
Behavior value, WithIndicate meaning to should refer toWithεpqValue is
[0,1], i indicate that the i-th column in behavioural matrix, j indicate jth row in behavioural matrix, 1≤j≤M, 1≤i≤N.
9. according to claim 1-3, travelling wave identification method based on three-dimensional grey absolute correlation degree described in any one of 5-7,
It is characterized in that, the travelling wave signal after the initial traveling wave arrival test point of the step acquisition physical fault in preset duration, packet
It includes:
Travelling wave signal after the initial traveling wave arrival test point of acquisition physical fault in 1ms time window respectively.
10. the traveling wave positioning device based on three-dimensional grey absolute correlation degree, which is characterized in that including acquisition module, obtain module,
Calculation of relationship degree module and identification module；
The acquisition module, the initial traveling wave for obtaining currently practical failure reach the traveling wave letter after test point in preset duration
Number, draw out the corresponding first when-frequency traveling wave Full wave shape of currently practical failure；
The acquisition module, for obtaining the total evidence of the second when-frequency traveling wave Full wave shape；
The calculation of relationship degree module, when for calculating described first-frequency traveling wave Full wave shape and the second when-frequency traveling wave all-wave
The three-dimensional grey absolute correlation degree of shape；
The identification module, for identifying and the described first when-frequency traveling wave all-wave according to the three-dimensional grey absolute correlation degree
It is fixed according to the failure of the described second when-frequency traveling wave Full wave shape position to occur for the highest second when-frequency traveling wave Full wave shape of the shape degree of association
The generation position of the currently practical failure in position.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811516653.XA CN109375065A (en) | 2018-12-12 | 2018-12-12 | Travelling wave identification method and positioning device based on three-dimensional grey absolute correlation degree |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811516653.XA CN109375065A (en) | 2018-12-12 | 2018-12-12 | Travelling wave identification method and positioning device based on three-dimensional grey absolute correlation degree |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109375065A true CN109375065A (en) | 2019-02-22 |
Family
ID=65373387
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811516653.XA Pending CN109375065A (en) | 2018-12-12 | 2018-12-12 | Travelling wave identification method and positioning device based on three-dimensional grey absolute correlation degree |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109375065A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110703128A (en) * | 2019-04-23 | 2020-01-17 | 武汉理工大学 | Single-phase earth fault phase selection method for non-effectively-grounded power distribution network based on grey correlation analysis |
| CN113447756A (en) * | 2021-05-24 | 2021-09-28 | 南方电网科学研究院有限责任公司 | Fault positioning method, device, equipment and medium for power distribution network line |
| CN114184886A (en) * | 2021-11-24 | 2022-03-15 | 昆明理工大学 | Method for quantizing complexity of fault traveling wave of power transmission line |
| CN114355112A (en) * | 2022-03-22 | 2022-04-15 | 智联新能电力科技有限公司 | Transient fault and defect discharge fault simulation test platform and data analysis method |
| CN117330890A (en) * | 2023-09-18 | 2024-01-02 | 浙江德清迪生电力科技有限公司 | Power transmission line fault diagnosis system and method |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101526433A (en) * | 2009-04-20 | 2009-09-09 | 江苏省交通科学研究院股份有限公司 | Method for automatically diagnosing faults of monitoring system |
| CN101975910A (en) * | 2010-09-07 | 2011-02-16 | 昆明理工大学 | Intelligent fault classification and location method for ultra-high voltage direct current transmission line |
| CN102680860A (en) * | 2012-06-08 | 2012-09-19 | 东华理工大学 | Automatic fault-point locating method for traveling-wave based fault location of high-voltage electric power lines |
| CN103487722A (en) * | 2013-08-02 | 2014-01-01 | 深圳市智远能科技有限公司 | Distributed cable state detection system based on TDR |
| CN103809061A (en) * | 2012-11-06 | 2014-05-21 | 西安元朔科技有限公司 | Technique for analyzing power transmission line fault type based on grey correlation analysis |
| KR101420567B1 (en) * | 2013-03-21 | 2014-07-17 | 울산대학교 산학협력단 | Method for Fault Classification of Induction Motors |
| CN104155582A (en) * | 2014-08-29 | 2014-11-19 | 上海交通大学 | Distribution network line fault section positioning method based on full-waveform information |
| CN104898024A (en) * | 2015-06-10 | 2015-09-09 | 国网上海市电力公司 | Transformer station grounding network fault diagnosis method based on synthesis of surface potential and magnetic induction strength |
| CN106405651A (en) * | 2016-11-14 | 2017-02-15 | 中国石油化工股份有限公司 | Logging-matching-based method for constructing full-waveform inversion initial model |
| CN107102236A (en) * | 2017-03-24 | 2017-08-29 | 国网浙江余姚市供电公司 | A kind of fault line selection method for single-phase-to-ground fault based on waveform correlation analysis after failure |
| CN107255743A (en) * | 2017-05-03 | 2017-10-17 | 昆明理工大学 | A kind of extra high voltage direct current transmission line lightning fault recognition methods based on power spectrum similarity |
| CN108429245A (en) * | 2018-03-23 | 2018-08-21 | 上海电力学院 | Wide area backup protection method based on sequence electric current grey correlation and Multi-information acquisition |
-
2018
- 2018-12-12 CN CN201811516653.XA patent/CN109375065A/en active Pending
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101526433A (en) * | 2009-04-20 | 2009-09-09 | 江苏省交通科学研究院股份有限公司 | Method for automatically diagnosing faults of monitoring system |
| CN101975910A (en) * | 2010-09-07 | 2011-02-16 | 昆明理工大学 | Intelligent fault classification and location method for ultra-high voltage direct current transmission line |
| CN102680860A (en) * | 2012-06-08 | 2012-09-19 | 东华理工大学 | Automatic fault-point locating method for traveling-wave based fault location of high-voltage electric power lines |
| CN103809061A (en) * | 2012-11-06 | 2014-05-21 | 西安元朔科技有限公司 | Technique for analyzing power transmission line fault type based on grey correlation analysis |
| KR101420567B1 (en) * | 2013-03-21 | 2014-07-17 | 울산대학교 산학협력단 | Method for Fault Classification of Induction Motors |
| CN103487722A (en) * | 2013-08-02 | 2014-01-01 | 深圳市智远能科技有限公司 | Distributed cable state detection system based on TDR |
| CN104155582A (en) * | 2014-08-29 | 2014-11-19 | 上海交通大学 | Distribution network line fault section positioning method based on full-waveform information |
| CN104898024A (en) * | 2015-06-10 | 2015-09-09 | 国网上海市电力公司 | Transformer station grounding network fault diagnosis method based on synthesis of surface potential and magnetic induction strength |
| CN106405651A (en) * | 2016-11-14 | 2017-02-15 | 中国石油化工股份有限公司 | Logging-matching-based method for constructing full-waveform inversion initial model |
| CN107102236A (en) * | 2017-03-24 | 2017-08-29 | 国网浙江余姚市供电公司 | A kind of fault line selection method for single-phase-to-ground fault based on waveform correlation analysis after failure |
| CN107255743A (en) * | 2017-05-03 | 2017-10-17 | 昆明理工大学 | A kind of extra high voltage direct current transmission line lightning fault recognition methods based on power spectrum similarity |
| CN108429245A (en) * | 2018-03-23 | 2018-08-21 | 上海电力学院 | Wide area backup protection method based on sequence electric current grey correlation and Multi-information acquisition |
Non-Patent Citations (5)
| Title |
|---|
| REZA JALILZADEH HAMIDI: "Traveling-Wave Detection Technique Using Short-Time Matrix Pencil Method", 《IEEE TRANSACTIONS ON POWER DELIVERY》 * |
| 刘思峰: "《灰色***理论及其应用》", 30 June 2014, 科学出版社 * |
| 徐松林: "基于小波能量谱和灰色综合关联度的HVDC换相失败故障诊断", 《电力***保护与控制》 * |
| 邓丰: "基于全波形信息的混联线路单端行波定位方法", 《电工技术学报》 * |
| 邓丰: "基于波形唯一和时–频特征匹配的单端行波保护和故障定位方法", 《中国电机工程学报》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110703128A (en) * | 2019-04-23 | 2020-01-17 | 武汉理工大学 | Single-phase earth fault phase selection method for non-effectively-grounded power distribution network based on grey correlation analysis |
| CN113447756A (en) * | 2021-05-24 | 2021-09-28 | 南方电网科学研究院有限责任公司 | Fault positioning method, device, equipment and medium for power distribution network line |
| CN114184886A (en) * | 2021-11-24 | 2022-03-15 | 昆明理工大学 | Method for quantizing complexity of fault traveling wave of power transmission line |
| CN114184886B (en) * | 2021-11-24 | 2024-06-07 | 昆明理工大学 | Method for quantifying complexity of fault traveling wave of power transmission line |
| CN114355112A (en) * | 2022-03-22 | 2022-04-15 | 智联新能电力科技有限公司 | Transient fault and defect discharge fault simulation test platform and data analysis method |
| CN117330890A (en) * | 2023-09-18 | 2024-01-02 | 浙江德清迪生电力科技有限公司 | Power transmission line fault diagnosis system and method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109375065A (en) | Travelling wave identification method and positioning device based on three-dimensional grey absolute correlation degree | |
| Shu et al. | Single-ended fault location for hybrid feeders based on characteristic distribution of traveling wave along a line | |
| CN105137293B (en) | The localization method of trouble spot in a kind of power distribution network joint line | |
| CN108508320B (en) | Arc grounding fault identification method based on harmonic energy and wave distortion feature | |
| CN110596539B (en) | Power distribution network fault positioning method based on transient waveform correlation | |
| CN103217626B (en) | A kind of single-ended traveling wave fault location method utilizing positive-negative polarity wave head timing intervals | |
| CN103323741B (en) | A kind of D molded line cable mixed line fault section compared based on false voltage initial row wave amplitude for strong fault sentences method for distinguishing | |
| CN103777083B (en) | Based on the capacitive equipment dielectric loss on-line monitoring method of Kalman's frequency-tracking | |
| CN106646121B (en) | A kind of discrimination method of distribution network failure wavefront | |
| CN105445624A (en) | Cable fault positioning method according to combination of wavelet transformation and curve fitting | |
| Deng et al. | A method for distribution network line selection and fault location based on a hierarchical fault monitoring and control system | |
| CN104833851B (en) | Harmonious Waves in Power Systems method of estimation based on distributed related Kalman filtering | |
| Ma et al. | A Method of line fault location based on traveling wave theory | |
| CN106093705A (en) | A kind of computational methods of one-phase earthing failure in electric distribution network wavefront | |
| CN103995181A (en) | Method for analyzing electric energy quality harmonic waves of digital substation | |
| CN106018958A (en) | Current voltage frequency tracking method of secondary side in isolated grid frequency conversion system | |
| CN110231554A (en) | A kind of 4 acoustic fix ranging methods of GIL/GIS breakdown fault | |
| CN107179473B (en) | A kind of power transmission line fault locating method | |
| CN112083263A (en) | Multi-source data fusion high-voltage cable line fault online positioning device and method based on traveling wave principle | |
| CN115201635A (en) | Multi-branch distribution network fault location method and device based on traveling waves | |
| CN103941152A (en) | Cable hybrid line fault distance measuring method for k-NN algorithm based on waveform similarity | |
| CN111007315B (en) | Circuit branch identification method based on current harmonic spectrum | |
| El Din et al. | A wavelet-based fault location technique for aged power cables | |
| Xu et al. | A novel single-ended traveling wave fault location method based on reflected wave-head of adjacent bus | |
| Shu et al. | Single-ended fault location for direct distribution overhead feeders based on characteristic distribution of traveling waves along the line |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190222 |
|
| RJ01 | Rejection of invention patent application after publication |