CN108732463A - A kind of improved wavefront extracting method - Google Patents
A kind of improved wavefront extracting method Download PDFInfo
- Publication number
- CN108732463A CN108732463A CN201810400009.XA CN201810400009A CN108732463A CN 108732463 A CN108732463 A CN 108732463A CN 201810400009 A CN201810400009 A CN 201810400009A CN 108732463 A CN108732463 A CN 108732463A
- Authority
- CN
- China
- Prior art keywords
- point
- wave
- waveform
- traveling
- virtual value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/085—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution lines, e.g. overhead
Abstract
The present invention relates to a kind of improved wavefront extracting methods, the element for traveling wave of being struck by lightning is divided into multi-section point by this method using unit vision window, virtual value is solved to various pieces and then obtains the wave head of single traveling wave, for eliminating influence of the sharp punching property noise to wave head recognition accuracy.Compared with prior art, the present invention has many advantages, such as that traveling-wave waveform signal, without particular/special requirement, method is easy.
Description
Technical field
The present invention relates to a kind of wavefront extracting methods, more particularly, to a kind of improved wavefront extracting method.
Background technology
With the growth of national economy, transmission line of electricity is increasing, in recent years, it is a kind of based on traveling wave positioning method obtain
Extensive use, this method key point are to find system monitoring to the wave head of the traveling wave corresponding time, pass through the time difference pair
Fault point is accurately positioned.Currently, both at home and abroad for there are many recognition methods of wavefront, such as derivative method, cross-correlation
The shortcomings that the methods of function method, morphological approach, Wavelet Transform, preceding several method is not to be suitable for non-linear non-stable letter
Number, and wavelet algorithm needs to select suitable wavelet basis, and the selection of wavelet basis has a great difficulty in practical applications, and this
The selection of kind wavelet basis is the whether suitable major criterion of signal, therefore Wavelet Transform has certain difficulty in practical applications
Degree.In addition, being commonly present certain electric discharge phenomena on transmission line of electricity, high frequency oscillation can be generated on the line, can be adulterated in measured waveform
High-frequency impulse interferes, and certain difficulty is brought to the extraction of transmission line travelling wave wave head.
Invention content
The purpose of the present invention aiming in current wavefront discrimination method both at home and abroad there are narrow application range, using tired
It is difficult to and present situation the deficiencies of anti-interference ability is poor, and a kind of improved wavefront extracting method is provided, this method can have
Effect identification wavefront, while the point punching property noise that can avoid accidentally occurring in traveling wave is influenced caused by looking for wave head, to traveling wave
For waveform signal without particular/special requirement, method is easy, is applicable to the identification of all wavefronts substantially.
The purpose of the present invention can be achieved through the following technical solutions：
The element for traveling wave of being struck by lightning is divided by a kind of improved wavefront extracting method, this method using unit vision window
Multi-section point solves virtual value to various pieces and then obtains the wave head of single traveling wave, knows to wave head for eliminating point punching property noise
The influence of other accuracy.
Preferably, the method specifically includes following steps：
1) it monitors in real time and extracts the traveling wave on transmission line of electricity；
2) by traveling-wave waveform according to the regular adding window of setting；
3) virtual value of traveling wave each point is calculated one by one；
4) it according to step 3), since first point of original waveform, calculates each point virtual value backward point by point, obtains traveling wave virtual value
Waveform；
5) the effective section of former traveling-wave waveform is divided；
6) extraction wave head sequentially finds maximum absolute value point in step 5) the virtual value section in each section, this
A little points are former wavefront.
Preferably, in the step 2) set rule as：The window of setting appropriate length, the length of Δ t, the value are big
It is small to set, for transmission line malfunction traveling wave, choose the when a length of 30-60 μ s of window.
Preferably, the computation rule of arbitrary i-th collection point virtual value is as follows in the step 3)：I-th collection point is made
For starting point, interception is from the starting point to the waveform segment of Δ t width below on the basis of former traveling-wave waveform, it is assumed that is intercepted
Segment points are n, and segment last point is original waveform jth point, then i-th virtual value I of original waveformrmsiFor：
Wherein xkFor the amplitude of the selected arbitrary kth point of waveform segment, k ∈ [i, j].
Preferably, the step 5) division principle is：Calculate virtual value waveforms amplitude Ip, and each maximum is sought successively
Point, if maximum point size > kIp, it is believed that the maximum point is effective, and k is a positive coefficient；At the time of with effective maximum point
On the basis of, each Δ t intercepts the effective section of original waveform forward and backward, it is assumed that certain maximum point moment is t0, then original waveform
[t0-Δt,t0+ Δ t] section is considered as effective section, all effective sections are sought successively.
Preferably, the k is taken as 0.5.
At this stage, the method for mainly using wavelet analysis in wavefront automatic identification field both at home and abroad, this method need
Suitable wavelet basis is first found, at part, the suitable wavelet basis of the searching is more difficult, or even needs oneself to create newly
Wavelet basis is difficult to carry out in the practical application of partial picture.Compared with prior art, the present invention is utilized adds frame to calculate to waveform
The principle of virtual value in unit frame can effectively give up the point punching property noise occurred once in a while in lightning stroke traveling wave, avoid it to system
In the influence of automatic identification wavefront, meanwhile, in conjunction with simple mathematical method, operation is simple, can be applied to various traveling wave waves
Shape, strong applicability.
Description of the drawings
Fig. 1 is the lightning fault current traveling wave schematic diagram with disturbing pulse；
Fig. 2 is travelling wave current virtual value waveform diagram.
Specific implementation mode
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation describes, it is clear that described embodiment is a part of the embodiment of the present invention, rather than whole embodiments.Based on this hair
Embodiment in bright, the every other reality that those of ordinary skill in the art are obtained without making creative work
Example is applied, the scope of protection of the invention should be all belonged to.
Fig. 1 is the lightning fault current traveling wave with the interference of high frequency spine of certain transmission line of electricity one actual measurement, wherein 1,2 be lightning stroke
The Main glow discharging pulse that moment is formed, i.e., main wave；3,4 back waves formed after substation or failure point reflection for main wave；5 are
A kind of common high-frequency impulse interference on transmission line of electricity.It can be seen from the figure that impulse disturbances amplitude has been more than traveling wave amplitude, use
It when conventional method carries out that wave head is looked for analyze, is easy to disturbing pulse 5 when decide wave or back wave, then passes through traveling wave and position
Analysis has obtained the positioning result of mistake so that fault location fails.
The improved wavefront extracting method of the present invention is implemented as follows.
(1) traveling wave on transmission line of electricity is monitored and extracts in real time, this traveling wave is made of N number of collection point, the width of each collection point
Value is followed successively by, and example waveform is as shown in Figure 1.
(2) by traveling-wave waveform according to following regular adding window：The window of appropriate length is set, the length of Δ t, which can
If for transmission line malfunction traveling wave, it is proposed that choose the when a length of 30-60 μ s of window.
(3) virtual value of traveling wave each point is calculated one by one.The computation rule of arbitrary i-th collection point virtual value is as follows：I-th is adopted
Collection point is used as starting point, and interception is from the starting point to the waveform segment of Δ t width below on the basis of former traveling-wave waveform, it is assumed that quilt
The segment points of interception are n, and segment last point is original waveform jth point, then i-th virtual value of original waveform is：
(4) (3) described rule is pressed, since first point of original waveform, calculates each point virtual value backward point by point, obtaining traveling wave has
Valid value waveform.Due to original waveform the last one when window width Δ t in collection point cannot be satisfied calculating when need to continue to choose backward
The requirement of the points of Δ t can not carry out virtual value calculating.From the point of view of practical experience, when being acquired relative to original waveform due to Δ t
Very little for length, and actual monitoring master involves back wave generally in the waveform middle front part of record, therefore when giving up final stage
Long Δ t waveform segments, will not lose original waveform important information.According to above-mentioned rule, node-by-node algorithm virtual value is carried out to Fig. 1,
Its virtual value waveform is obtained, as shown in Figure 2：
(5) the effective section of former traveling-wave waveform is divided.Division principle is as follows：Calculate virtual value waveforms amplitude Ip, and ask successively
Each maximum point is taken, if maximum point size > kIp, it is believed that the maximum point is effective, and k is a positive coefficient, and empirically value is desirable
It is 0.5.On the basis of at the time of effective maximum point, each Δ t intercepts the effective section of original waveform forward and backward, it is assumed that certain
The maximum point moment is t0, then original waveform [t0-Δt,t0+ Δ t] section is considered as effective section, it seeks successively all
Effective section.In exemplary plot 2,1,2,3 be maximum point, but only 1 and 2 meet effective maximum point requirement, with 1 and 2 stroke
The virtual value section divided just includes that original waveform master involves back wave part, and High-frequency Interference is then excluded.
(6) wave head is extracted.In (5) described virtual value section, maximum absolute value point is sequentially found in each section, these
Point is former wavefront.
The above description is merely a specific embodiment, but scope of protection of the present invention is not limited thereto, any
Those familiar with the art in the technical scope disclosed by the present invention, can readily occur in various equivalent modifications or replace
It changes, these modifications or substitutions should be covered by the protection scope of the present invention.Therefore, protection scope of the present invention should be with right
It is required that protection domain subject to.
Claims (6)
1. a kind of improved wavefront extracting method, which is characterized in that this method utilizes unit vision window by traveling wave of being struck by lightning
Element is divided into multi-section point, solves virtual value to various pieces and then obtains the wave head of single traveling wave, makes an uproar for eliminating point punching property
Influence of the sound to wave head recognition accuracy.
2. according to the method described in claim 1, it is characterized in that, the method specifically includes following steps：
1) it monitors in real time and extracts the traveling wave on transmission line of electricity；
2) by traveling-wave waveform according to the regular adding window of setting；
3) virtual value of traveling wave each point is calculated one by one；
4) it according to step 3), since first point of original waveform, calculates each point virtual value backward point by point, obtains traveling wave virtual value wave
Shape；
5) the effective section of former traveling-wave waveform is divided；
6) extraction wave head sequentially finds maximum absolute value point in step 5) the virtual value section in each section, these points
As former wavefront.
3. according to the method described in claim 2, it is characterized in that, in the step 2) set rule as：Setting is suitable
The window of length, the length of Δ t, which can set, and for transmission line malfunction traveling wave, choose the when a length of 30- of window
60μs。
4. according to the method described in claim 2, it is characterized in that, arbitrary i-th collection point virtual value in the step 3)
Computation rule is as follows：Using the i-th collection point as starting point, interception is from the starting point to Δ t below on the basis of former traveling-wave waveform
The waveform segment of width, it is assumed that the segment points being intercepted are n, and segment last point is original waveform jth point, then original waveform the
I point virtual values IrmsiFor：
Wherein xkFor the amplitude of the selected arbitrary kth point of waveform segment, k ∈ [i, j].
5. according to the method described in claim 2, it is characterized in that, the step 5) division principle is：Calculate virtual value wave
Shape amplitude Ip, and each maximum point is sought successively, if maximum point size > kIp, it is believed that the maximum point is effective, and k is a positive system
Number；On the basis of at the time of effective maximum point, each Δ t intercepts the effective section of original waveform forward and backward, it is assumed that certain pole
Big value point moment is t0, then original waveform [t0-Δt,t0+ Δ t] section is considered as effective section, it seeks all having successively
Imitate section.
6. according to the method described in claim 5, it is characterized in that, the k is taken as 0.5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810400009.XA CN108732463B (en) | 2018-04-28 | 2018-04-28 | Improved traveling wave head extraction method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810400009.XA CN108732463B (en) | 2018-04-28 | 2018-04-28 | Improved traveling wave head extraction method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108732463A true CN108732463A (en) | 2018-11-02 |
| CN108732463B CN108732463B (en) | 2020-10-09 |
Family
ID=63940077
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810400009.XA Active CN108732463B (en) | 2018-04-28 | 2018-04-28 | Improved traveling wave head extraction method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108732463B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109901016A (en) * | 2019-03-13 | 2019-06-18 | 清华四川能源互联网研究院 | A kind of both-end distance measuring method eliminating velocity of wave and influencing |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101291055A (en) * | 2008-06-18 | 2008-10-22 | 昆明理工大学 | Method for precisely marking arriving time of initial wave of fault generated traveling waves for electricity transmission line |
| CN104931855A (en) * | 2015-06-10 | 2015-09-23 | 深圳市索图科技有限公司 | Power transmission line fault-based traveling wave header identification and extraction device and method |
| CN105938173A (en) * | 2016-04-14 | 2016-09-14 | 云南电网有限责任公司电力科学研究院 | Precise range finding method for fault of ultrahigh-voltage line |
| CN106771877A (en) * | 2017-01-11 | 2017-05-31 | 北京衡天北斗科技有限公司 | The determination method and apparatus of the position of failure point of system with non effectively earth ed neutral |
-
2018
- 2018-04-28 CN CN201810400009.XA patent/CN108732463B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101291055A (en) * | 2008-06-18 | 2008-10-22 | 昆明理工大学 | Method for precisely marking arriving time of initial wave of fault generated traveling waves for electricity transmission line |
| CN104931855A (en) * | 2015-06-10 | 2015-09-23 | 深圳市索图科技有限公司 | Power transmission line fault-based traveling wave header identification and extraction device and method |
| CN105938173A (en) * | 2016-04-14 | 2016-09-14 | 云南电网有限责任公司电力科学研究院 | Precise range finding method for fault of ultrahigh-voltage line |
| CN106771877A (en) * | 2017-01-11 | 2017-05-31 | 北京衡天北斗科技有限公司 | The determination method and apparatus of the position of failure point of system with non effectively earth ed neutral |
Non-Patent Citations (3)
| Title |
|---|
| 吴仲朗 等: "基于离散求导法的铁路电力线路故障行波波头识别新算法", 《电力系统保护与控制》 * |
| 李传健 等: "基于卡尔曼滤波的行波波头检测算法研究", 《电网与清洁能源》 * |
| 董原竹 等: "基于希尔伯特 - 黄变换的行波波头识别算法研究", 《湖南农机》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109901016A (en) * | 2019-03-13 | 2019-06-18 | 清华四川能源互联网研究院 | A kind of both-end distance measuring method eliminating velocity of wave and influencing |
| CN109901016B (en) * | 2019-03-13 | 2020-11-20 | 清华四川能源互联网研究院 | Double-end ranging method for eliminating wave velocity influence |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108732463B (en) | 2020-10-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108508320B (en) | Arc grounding fault identification method based on harmonic energy and wave distortion feature | |
| CN108469560B (en) | Electromagnetic interference objective complexity evaluation method based on rapid S-transform time-frequency space model | |
| CN110852201B (en) | Pulse signal detection method based on multi-pulse envelope spectrum matching | |
| CN104730424B (en) | Cable local discharge localization method based on auto-correlation-Wavelet Modulus Maxima analysis | |
| CN107991706B (en) | Coal seam hydraulic fracturing microseismic signal combined noise reduction method based on wavelet packet multiple threshold and improved empirical mode decomposition | |
| CN102135558B (en) | Acquisition and hierarchical identification system of overvoltage data and hierarchical pattern identification method of overvoltage types | |
| CN111313932A (en) | Transmission line based on electromagnetic time reversal transfer function correlation and network fault positioning method | |
| CN103698669B (en) | A kind of high-tension electricity apparatus local discharge pulse separation method | |
| CN107086566B (en) | LMD interconnected electric power system low-frequency oscillation analysis method based on Wide-area Measurement Information | |
| CN106646205A (en) | Random big-disturbance signal removing algorithm for analyzing circuit breaker fault through sound and vibration combination | |
| CN102680860A (en) | Automatic fault-point locating method for traveling-wave based fault location of high-voltage electric power lines | |
| CN110068759A (en) | A kind of fault type preparation method and device | |
| CN108344922B (en) | Power transmission line direct lightning strike fault identification method based on phase classification and S transformation | |
| CN105067963A (en) | Distribution network fault early warning method and system based on transient waveforms | |
| CN107632225A (en) | A kind of small current system Earth design method | |
| CN110247420B (en) | Intelligent fault identification method for HVDC transmission line | |
| CN108199356B (en) | Wave-front information-based DC transmission line ultra-high speed protection method | |
| CN109270346A (en) | Polarization state waveform signal processing method at a kind of OPGW lightning stroke | |
| CN104635126A (en) | Local discharge single-pulse extraction method based on sliding window | |
| CN105893976B (en) | A kind of parametrization recognition methods of travelling wave signal | |
| CN110412467A (en) | A kind of lithium battery fault data screening technique of normalized mutual information criterion constraint | |
| CN109085477A (en) | Signal identification and localization method for power cable distribution partial discharge monitoring system | |
| CN109239551A (en) | Local discharge signal temporal signatures extracting method based on Chirplet transformation | |
| CN106885975A (en) | A kind of high-tension cable ageing testing method and device based on impulse response | |
| CN108732463A (en) | A kind of improved wavefront extracting 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 |