CN102998593A - Method for accurately positioning electrical cable faults through interpolation convolution calculating - Google Patents
Method for accurately positioning electrical cable faults through interpolation convolution calculating Download PDFInfo
- Publication number
- CN102998593A CN102998593A CN2012104613180A CN201210461318A CN102998593A CN 102998593 A CN102998593 A CN 102998593A CN 2012104613180 A CN2012104613180 A CN 2012104613180A CN 201210461318 A CN201210461318 A CN 201210461318A CN 102998593 A CN102998593 A CN 102998593A
- Authority
- CN
- China
- Prior art keywords
- interpolation
- point
- waveform
- wave
- sequence
- 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
Images
Landscapes
- Locating Faults (AREA)
Abstract
The invention relates to a method for accurately positioning electrical cable faults through interpolation convolution calculating and belongs to the technical field of electric communication. The method includes: selecting wave form data x of wave forms such as current waves, voltage waves and electromagnetic waves during operation faults; drawing two-dimensional diagram with Matlab; performing convolution calculating to the wave form data to position faults primarily; performing interpolation to wave form data, and performing convolution calculating again to position faults accurately. The related wave form formula is as followed, wherein X and Y in Rxy are two values of sequence relevance, Xn and Yn are corresponding numerical values in the sequence, m and n refers to ordinal numbers in the sequence, Y* refers to conjugate of numbers in the sequence, and E refers to expectation value. The point with most relevance after the second interpolation is the fault point accurately positioned. Faults can be more accurately positioned on the premise of not improving hardware conditions. Actual fault positioning precision can be improved under the condition of not improving hardware sampling rate.
Description
Technical field
A kind ofly accurately locate the method for electrical equipment fault position by the interpolation convolutional calculation, relate in particular to the location such as the method for GIS, power cable, pole line distributed equipment discharging fault, belong to the electricity communication technical field.
Background technology
Power equipment, especially such as GIS, power cable, pole line distributed equipment, all probability generation discharging fault necessarily will produce discharging fault such as power cable power equipment etc.
Power cable fault location can be divided into two large classes on principle: traveling wave method and impedance method.Existing technology is owing to row ripple sampling rate is restricted, so that localization of fault is a larger scope.The present invention mainly is based on traveling wave method, proposes a kind of convolution interpolation of utilizing, and improves the method for accuracy of detection.And can greatly dwindle fault coverage by present technique, for the fast searching trouble spot provides larger help.
Fault Location With Traveling Wave utilization row ripple can obtain distance in the time of measurement point to round trip between the trouble spot through simple operation.Traveling wave method specifically is divided into two classes: the first kind is to utilize the method for voltage traveling wave signal.Its method is that Method of Cable Trouble Point is punctured under the effect of high direct voltage or high voltage pulse signal, then by observing the time range finding of discharge pulse round trip between observation point and trouble spot.Equations of The Second Kind is the distance-finding method that adopts the current traveling wave signal.The difference of these class methods and last class methods is that the former measures the current signal that produces when cable fault punctures by a linear current coupling mechanism, and the current waveform that sensors coupled goes out also is easier to differentiate.The method that we adopt is that current traveling wave carries out fault localization.Its distance accuracy is not mainly by the formula that affects its foundation of the factors such as the abort situation, fault type, line length, stake resistance of circuit substantially:
In the formula: l
1And l
2Be respectively the trouble spot to the distance at two ends; T arrives the time at circuit two ends for the row ripple.We can see from formula, and the key that finds the trouble spot is that the mistiming of finding current wave to arrive (is t
1-t
2).And the waveform that arrives the cable two ends is the current wave that send same trouble spot in fact, seeks the correlativity maximum point of two waveforms that arrive two ends, finds exactly the auto-correlation maximum point of same waveform.The formula of determining the waveform correlativity is
In the formula: R
XyThe value of X, two serial correlations of Y, X
nWith Y
nBe each corresponding in sequence numerical value, m and n represent which number in the sequence, Y
*Represent the conjugation of numerical value in the ordered series of numbers, E represents expectation value.
As can be known, the auto-correlation maximum point that finds is exactly that the two ends waveforms arrives counting of differing (if sampling rate is 1Mhz, the current wave velocity of propagation is 3 * 108m/s, then distance L=the V of every point-to-point transmission * T=V/f=300m) in fact from formula.Count and the distance of actual point-to-point transmission representative according to differing, both can obtain position of failure point.
Summary of the invention
The present invention carries out interpolation and convolutional calculation by to behind the signal sampling, can improve the localization of fault precision, for fault accurately provides very large help in the location.
A kind ofly accurately locate the method for power cable fault position by the interpolation convolutional calculation, the method step is as follows:
Waveform when 1, choosing the actual motion fault (waveforms such as current wave, voltage wave, electromagnetic wave can) data x; Draw the data X-Y scheme with Matlab;
2, Wave data is carried out convolutional calculation, the Primary Location abort situation;
3, Wave data is carried out interpolation, again carry out convolutional calculation, accurately locate abort situation.
The waveform that arrives the cable two ends is the current wave that send same trouble spot in fact, seeks the correlativity maximum point of two waveforms that arrive two ends, finds exactly the auto-correlation maximum point of same waveform.The formula of determining the waveform correlativity is
, wherein: R
XyThe value of X, two serial correlations of Y, X
nWith Y
nBe each corresponding in sequence numerical value, m and n represent which number in the sequence, Y
*Represent the conjugation of numerical value in the ordered series of numbers, E represents expectation value.As can be known, the auto-correlation maximum point that finds is exactly that the two ends waveforms arrives counting of differing (if sampling rate is 1Mhz, the current wave velocity of propagation is 3 * 108m/s, then distance L=the V of every point-to-point transmission * T=V/f=300m) in fact from formula.Count and the distance of actual point-to-point transmission representative according to differing, both can obtain position of failure point.
Correlativity maximum point after the interpolation is exactly pinpoint trouble spot again, like this, by this method, improves sampling rate, more accurately locates position of failure point in the situation that does not improve hardware condition.
The present invention can improve the physical fault bearing accuracy in the situation that does not improve the hardware sampling rate.
Description of drawings
Fig. 1 is process flow diagram of the present invention.
Fig. 2 is current waveform data time-domain diagrams.
Fig. 3 is as a result figure of Wave data convolutional calculation.
Fig. 4 is data plot behind first group of data interpolating.
Two groups of data dependence result of calculation figure after Fig. 5 interpolation.
Embodiment
Further specify the present invention below in conjunction with embodiment.
Fig. 1 is process flow diagram of the present invention.A kind ofly accurately locate the method for power equipment position by the interpolation convolutional calculation, the method step is as follows:
Current wave data x(sampling rate when choosing actual motion can be selected according to the equipment actual conditions, if the higher then measuring accuracy of equipment sampling rate is higher.The data sampling rate that present embodiment lists below is 100Ms/s, 40 points of one group of data).Draw first group of data X-Y scheme with Matlab, such as Fig. 2.Fig. 2 is current waveform data time-domain diagrams.The sampling that Fig. 2 carries out the fault in-situ continuous wave for us, sampled data is 40 points altogether.Horizontal ordinate is for counting, and ordinate is the current amplitude of every bit.
(2) two groups of data are carried out convolutional calculation, obtain the correlativity maximum point, (Fig. 3 uses the xcorr function among the matlab that first group of data and second group of data are carried out result of calculation figure after the correlation calculations such as Fig. 3.Horizontal ordinate is two groups of data relative positions, and two groups of data of 0 expression are fully corresponding, two groups of data of 1 expression point that staggers; Ordinate is the correlation calculations value), convolution value maximum point is at the 40th point, and namely two groups of data are maximum at the 40th spot correlation.
(3) (insert 4 points between per two points, present embodiment is for conveniently only inserting 4 points, and several factually border, insertion point precision needs situation to select during actual the use two groups of data to be carried out interpolation.Insertion is counted, and measuring accuracy is higher more at most.), first group of data after the interpolation such as Fig. 4, Fig. 4 is data plot behind first group of data interpolating (horizontal ordinate is for counting, and ordinate is the current amplitude of every bit),
First group
Second group
(4) two groups of data after the interpolation are carried out convolutional calculation again, obtain result such as Fig. 5, two groups of data dependence result of calculation figure after Fig. 5 interpolation, maximum the 201st convolution value.
From result of calculation as can be known, by rational sampling rate, can accurately express the correlativity between the waveform, greatly reduce again sampling rate simultaneously.
Improve sampling rate by software approach like this, it is more accurate that our being positioned at do not improved in the situation of hardware condition.
In second step, our sampling rate is 100Ms/s, if the position of trouble spot between the 40th point and the 41st point, fault coverage is 3m like this.After interpolation convolutional calculation again, the trouble spot is the 200th point and 201, and fault coverage is 0.6m.So just sampling rate 5 times have been improved
The present invention by computing method, improves the physical fault bearing accuracy in the situation that does not improve the hardware sampling rate.
Claims (2)
1. accurately locate the method for power cable fault position by the interpolation convolutional calculation for one kind, it is characterized in that the method step is as follows:
Wave data x when (1) choosing the actual motion fault; Draw the data X-Y scheme with Matlab;
(2) Wave data is carried out convolutional calculation, the Primary Location abort situation;
(3) Wave data is carried out interpolation, again carry out convolutional calculation, accurately locate abort situation;
The waveform that arrives the cable two ends is the current wave that send same trouble spot, seeks the correlativity maximum point of two waveforms that arrive two ends, finds exactly the auto-correlation maximum point of same waveform, determines that the formula of waveform correlativity is
Wherein, R
XyThe value of X, two serial correlations of Y, X
nWith Y
nBe each corresponding in sequence numerical value, m and n represent which number in the sequence, Y
*Represent the conjugation of numerical value in the ordered series of numbers, E represents expectation value; As can be known, the auto-correlation maximum point that finds is exactly that the two ends waveform arrives counting of differing in fact from formula; Count and the distance of actual point-to-point transmission representative according to differing, both can obtain position of failure point;
Correlativity maximum point after the interpolation is exactly pinpoint trouble spot again.
2. according to claim 1ly a kind ofly accurately locate the method for power cable fault position by the interpolation convolutional calculation, it is characterized in that described waveform is current wave, voltage wave or electromagnetic wave waveform.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210461318.0A CN102998593B (en) | 2012-11-15 | 2012-11-15 | Method for accurately positioning electrical cable faults through interpolation convolution calculating |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210461318.0A CN102998593B (en) | 2012-11-15 | 2012-11-15 | Method for accurately positioning electrical cable faults through interpolation convolution calculating |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102998593A true CN102998593A (en) | 2013-03-27 |
| CN102998593B CN102998593B (en) | 2015-04-15 |
Family
ID=47927412
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210461318.0A Expired - Fee Related CN102998593B (en) | 2012-11-15 | 2012-11-15 | Method for accurately positioning electrical cable faults through interpolation convolution calculating |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102998593B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106841923A (en) * | 2017-02-27 | 2017-06-13 | 广州思泰信息技术有限公司 | Distribution network line fault localization method based on difference Convolution Analysis method |
| CN109564258A (en) * | 2016-08-17 | 2019-04-02 | 通用电器技术有限公司 | About or being located in improvement in the failure in power transmission pipeline |
| CN110865286A (en) * | 2019-11-25 | 2020-03-06 | 武汉三相电力科技有限公司 | Discharge phase calculation method |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102004209A (en) * | 2010-12-03 | 2011-04-06 | 丹东欣泰电气股份有限公司 | Distribution network cable fault on-line distance measuring equipment and distance measuring method |
| US20110227582A1 (en) * | 2010-03-19 | 2011-09-22 | Honeywell International Inc. | Methods for detecting a hidden peak in wire fault location applications - improving the distance range resolution |
-
2012
- 2012-11-15 CN CN201210461318.0A patent/CN102998593B/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20110227582A1 (en) * | 2010-03-19 | 2011-09-22 | Honeywell International Inc. | Methods for detecting a hidden peak in wire fault location applications - improving the distance range resolution |
| CN102004209A (en) * | 2010-12-03 | 2011-04-06 | 丹东欣泰电气股份有限公司 | Distribution network cable fault on-line distance measuring equipment and distance measuring method |
Non-Patent Citations (3)
| Title |
|---|
| 王星海 等: "基于小波变换和自相关分析的电力电缆故障测距", 《继电器》, vol. 33, no. 12, 16 June 2005 (2005-06-16), pages 32 - 35 * |
| 范毅 等: "小波变换在电缆故障定位中的应用", 《高电压技术》, vol. 26, no. 4, 30 August 2000 (2000-08-30) * |
| 许珉 等: "电力电缆故障低压脉冲自动测距方法", 《继电器》, vol. 35, no. 7, 1 April 2007 (2007-04-01), pages 37 - 40 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109564258A (en) * | 2016-08-17 | 2019-04-02 | 通用电器技术有限公司 | About or being located in improvement in the failure in power transmission pipeline |
| US11474140B2 (en) | 2016-08-17 | 2022-10-18 | General Electric Technology Gmbh | Relating to locating faults in power transmission conduits |
| CN106841923A (en) * | 2017-02-27 | 2017-06-13 | 广州思泰信息技术有限公司 | Distribution network line fault localization method based on difference Convolution Analysis method |
| CN110865286A (en) * | 2019-11-25 | 2020-03-06 | 武汉三相电力科技有限公司 | Discharge phase calculation method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102998593B (en) | 2015-04-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103176107B (en) | High-voltage direct-current power transmission line hybrid fault ranging method | |
| CN108107319A (en) | A kind of multiterminal flexible direct current electric network fault localization method and system | |
| CN102426393B (en) | Electric prospecting method | |
| CN101718833B (en) | Method of single end distance measurement of power transmission line malfunction based on traveling wave inherent frequency extraction | |
| CN101666848B (en) | Fault location method based on travelling wave | |
| CN107632236A (en) | A kind of single outgoing-feeder line one-end fault ranging method based on the identification of opposite end bus back wave | |
| CN100480734C (en) | High resolution static frequency domain removed magnetotelluric method | |
| CN105607131B (en) | Coding source electromagnetic depth method obtains the method and device of the earth APPARENT RESISTIVITY spectrum | |
| CN103777115B (en) | Electric transmission line single-terminal positioning method based on fault transient state and steady-state signal wave velocity difference | |
| CN103207413B (en) | Electrical prospecting device and system | |
| CN103941153B (en) | A kind of outgoing lines radiation network fault distance-finding method of k-NN algorithm based on waveform similarity | |
| CN103941159B (en) | Mixed line fault location method | |
| CN104931855B (en) | The device and method for being identified and being extracted based on powerline fault traveling wave head | |
| CN202330736U (en) | Electric-method exploration device | |
| CN101489238B (en) | Time difference measuring method, system and apparatus | |
| CN103837799B (en) | A kind of frequency domain method of voltage DC ground electrode circuit fault based on R-L model range finding | |
| CN105116295A (en) | Direct distribution overhead line fault range finding method based on traveling wave abrupt change distance calibration | |
| CN101833058A (en) | Earth fault distance measurement method based on traveling wave modulus time difference | |
| CN103760460B (en) | A kind of shape high voltage DC earthing pole line time-domain fault distance-finding method based on Bei Jielong model | |
| CN103235237A (en) | High-voltage direct-current grounding electrode line high-resistance fault positioning method | |
| CN105759178A (en) | Double-end range finding method for single-phase earth fault of overhead-cable mixed line | |
| CN110927510A (en) | Frequency domain method for power transmission line double-end traveling wave fault location | |
| CN104977504A (en) | Cable fault on-line detecting and positioning device | |
| CN102998593B (en) | Method for accurately positioning electrical cable faults through interpolation convolution calculating | |
| CN102095929B (en) | Method for rapidly measuring frequency of alternating-current signals |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150415 Termination date: 20191115 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |