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
- point
- fault
- correlation
- interpolation
- waveform
- 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
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000004364 calculation method Methods 0.000 claims abstract description 24
- 238000005070 sampling Methods 0.000 abstract description 17
- 238000004891 communication Methods 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 238000003780 insertion Methods 0.000 description 2
- 230000037431 insertion Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 1
Images
Landscapes
- Locating Faults (AREA)
Abstract
一种通过插值卷积计算精确定位电力电缆故障位置的方法, 属于电学通信技术领域。选取实际运行故障时的波形如电流波、电压波、电磁波等波形数据x;用Matlab画出数据二维图;对波形数据进行卷积计算,初步定位故障位置;对波形数据进行插值,再次进行卷积计算,精确定位故障位置;波形相关性公式:
,其中,Rxy是X、Y两个序列相关性的值,Xn与Yn是序列中对应的每一个数值,m与n代表序列中的第几个数,Y*代表数列中数值的共轭,E代表期望值,再次插值之后的相关性最大点就是精确定位的故障点。这样,通过本方法,在不提高硬件条件的情况下更加精确定位故障点位置。本发明在不提高硬件采样率的情况下,能够提高实际故障定位精度。A method for accurately locating a power cable fault location through interpolation and convolution calculations belongs to the technical field of electrical communication. Select waveform data x such as current wave, voltage wave, electromagnetic wave, etc. during the actual operation fault; use Matlab to draw a two-dimensional graph of the data; perform convolution calculation on the waveform data, and initially locate the fault location; interpolate the waveform data, and then perform Convolution calculation, precise fault location; waveform correlation formula:
, where R xy is the value of the correlation between X and Y sequences, X n and Y n are each corresponding value in the sequence, m and n represent the number in the sequence, and Y * represents the number of values in the sequence Conjugate, E represents the expected value, and the point with the greatest correlation after interpolation is the precisely located fault point. In this way, through this method, the location of the fault point can be located more precisely without improving the hardware condition. The present invention can improve the accuracy of actual fault location without increasing the hardware sampling rate.Description
技术领域 technical field
一种通过插值卷积计算精确定位电力设备故障位置的方法,尤其涉及定位如GIS、电力电缆、架空线等分布式设备放电性故障的方法,属于电学通信技术领域。The invention relates to a method for accurately locating fault locations of power equipment through interpolation and convolution calculation, in particular to a method for locating discharge faults of distributed equipment such as GIS, power cables, and overhead lines, and belongs to the technical field of electrical communication.
背景技术 Background technique
电力设备,尤其如GIS、电力电缆、架空线等分布式设备,都会一定概率发生放电性故障,如电力电缆电力设备等就会产生放电性故障。Power equipment, especially distributed equipment such as GIS, power cables, overhead lines, etc., will have a certain probability of discharge faults, such as power cables and power equipment, etc. will produce discharge faults.
电力电缆故障测距在原理上可分为两大类:行波法和阻抗法。现有的技术由于行波采样率受到限制,使得故障定位在一个较大的范围。本发明主要是基于行波法,提出一种利用卷积插值,提高检测精度的方法。而通过本技术可以大大缩小故障范围,为快速寻找故障点提供较大帮助。Power cable fault location can be divided into two categories in principle: traveling wave method and impedance method. Due to the limitation of traveling wave sampling rate in the existing technology, the fault location is in a larger range. The present invention is mainly based on the traveling wave method, and proposes a method for improving detection precision by using convolution interpolation. However, this technology can greatly reduce the scope of faults and provide great help for quickly finding fault points.
行波法测距利用行波在测量点到故障点之间往返一次的时间,经过简单运算即可得到距离。行波法具体分为两类:第一类是利用电压行波信号的方法。其方法是使电缆故障点在直流高压或脉冲高压信号的作用下击穿,然后通过观察放电脉冲在观察点与故障点之间往返一次的时间测距。第二类是采用电流行波信号的测距方法。这类方法与前一类方法的区别在于,前者通过一线性电流耦合器测量电缆故障击穿时产生的电流信号,传感器耦合出的电流波形也较容易分辨。我们采用的方法是电流行波进行故障测距。其测距精度基本不受线路的故障位置、故障类型、线路长度、接地电阻等因素的影响它依据的公式主要为:Traveling wave distance measurement uses the time for traveling waves to go back and forth between the measurement point and the fault point, and the distance can be obtained after simple calculation. The traveling wave method is specifically divided into two categories: the first category is a method using a voltage traveling wave signal. The method is to make the cable fault point break down under the action of DC high voltage or pulse high voltage signal, and then measure the distance by observing the time when the discharge pulse goes back and forth between the observation point and the fault point. The second type is the ranging method using current traveling wave signal. The difference between this method and the former method is that the former uses a linear current coupler to measure the current signal generated when the cable fault breaks down, and the current waveform coupled by the sensor is also easier to distinguish. The method we adopt is current traveling wave for fault location measurement. Its ranging accuracy is basically not affected by factors such as line fault location, fault type, line length, grounding resistance, etc. The formula it is based on is mainly:
公式中:l1和l2分别为故障点到两端的距离;t为行波到达线路两端的时间。从公式中我们可以看到,找到故障点的关键是找到电流波到达的时间差(即t1-t2)。而到达电缆两端的波形其实是同一故障点发出的电流波,寻找到达两端的两个波形的相关性最大点,就是找到同一波形的自相关最大点。确定波形相关性的公式为。公式中:Rxy是X、Y两个序列相关性的值,Xn与Yn是序列中对应的每一个数值,m与n代表序列中的第几个数,Y*代表数列中数值的共轭,E代表期望值。In the formula: l 1 and l 2 are the distances from the fault point to the two ends respectively; t is the time for the traveling wave to reach the two ends of the line. We can see from the formula that the key to finding the fault point is to find the time difference between the arrival of the current wave (ie t 1 -t 2 ). The waveforms arriving at both ends of the cable are actually current waves from the same fault point. Finding the maximum correlation point of the two waveforms reaching both ends is to find the maximum autocorrelation point of the same waveform. The formula for determining waveform correlation is . In the formula: R xy is the value of the correlation between the two sequences of X and Y, X n and Y n are each corresponding value in the sequence, m and n represent the number in the sequence, Y * represents the value of the value in the sequence Conjugate, E stands for expected value.
从公式中可知,找到的自相关最大点其实就是两端波形到达相差的点数(若采样率为1Mhz,电流波传播速度为3×108m/s,则每两点间的距离L=V×T=V/f=300m)。根据相差点数以及实际两点间代表的距离,既可求出故障点位置。It can be seen from the formula that the maximum autocorrelation point found is actually the number of points where the waveforms at both ends differ (if the sampling rate is 1Mhz, and the current wave propagation speed is 3×108m/s, then the distance between each two points L=V×T =V/f=300m). According to the number of difference points and the distance represented by the actual two points, the location of the fault point can be obtained.
发明内容 Contents of the invention
本发明通过对信号采样后,进行插值与卷积计算,可以提高故障定位精度,为故障准确定位提供很大的帮助。The invention performs interpolation and convolution calculation after the signal is sampled, which can improve the accuracy of fault location and provide great help for accurate fault location.
一种通过插值卷积计算精确定位电力电缆故障位置的方法,该方法步骤如下:A method for accurately locating a fault location of a power cable through interpolation convolution calculation, the steps of the method are as follows:
1、选取实际运行故障时的波形(电流波、电压波、电磁波等波形都可以)数据x;用Matlab画出数据二维图;1. Select the waveform (current wave, voltage wave, electromagnetic wave, etc.) data x of the actual operation failure; use Matlab to draw a two-dimensional graph of the data;
2、对波形数据进行卷积计算,初步定位故障位置;2. Carry out convolution calculation on the waveform data, and initially locate the fault location;
3、对波形数据进行插值,再次进行卷积计算,精确定位故障位置。3. Interpolate the waveform data, perform convolution calculation again, and accurately locate the fault location.
到达电缆两端的波形其实是同一故障点发出的电流波,寻找到达两端的两个波形的相关性最大点,就是找到同一波形的自相关最大点。确定波形相关性的公式为,其中:Rxy是X、Y两个序列相关性的值,Xn与Yn是序列中对应的每一个数值,m与n代表序列中的第几个数,Y*代表数列中数值的共轭,E代表期望值。从公式中可知,找到的自相关最大点其实就是两端波形到达相差的点数(若采样率为1Mhz,电流波传播速度为3×108m/s,则每两点间的距离L=V×T=V/f=300m)。根据相差点数以及实际两点间代表的距离,既可求出故障点位置。The waveforms arriving at both ends of the cable are actually current waves from the same fault point. Finding the maximum correlation point of the two waveforms reaching both ends is to find the maximum autocorrelation point of the same waveform. The formula for determining waveform correlation is , where: R xy is the value of the correlation between the two sequences of X and Y, X n and Y n are each corresponding value in the sequence, m and n represent the number in the sequence, Y * represents the number of values in the sequence Conjugate, E stands for expected value. It can be seen from the formula that the maximum autocorrelation point found is actually the number of points where the waveforms at both ends differ (if the sampling rate is 1Mhz, and the current wave propagation speed is 3×108m/s, then the distance between each two points L=V×T =V/f=300m). According to the number of difference points and the distance represented by the actual two points, the location of the fault point can be obtained.
再次插值之后的相关性最大点就是精确定位的故障点,这样,通过本方法,提高采样率,在不提高硬件条件的情况下更加精确定位故障点位置。The point with the greatest correlation after interpolation is the precisely located fault point. In this way, the sampling rate can be increased through this method, and the fault point can be located more accurately without improving the hardware condition.
本发明在不提高硬件采样率的情况下,能够提高实际故障定位精度。The present invention can improve the accuracy of actual fault location without increasing the hardware sampling rate.
附图说明 Description of drawings
图1是本发明的流程图。Fig. 1 is a flow chart of the present invention.
图2是电流波形数据时域图。Figure 2 is a time domain diagram of the current waveform data.
图3是波形数据卷积计算结果图。Fig. 3 is a diagram of the calculation results of waveform data convolution.
图4是第一组数据插值后数据图。Fig. 4 is a data diagram after interpolation of the first group of data.
图5插值后两组数据相关性计算结果图。Figure 5. The calculation results of the correlation between the two groups of data after interpolation.
具体实施方式 Detailed ways
下面结合具体实施方式进一步说明本发明。The present invention will be further described below in combination with specific embodiments.
图1是本发明的流程图。一种通过插值卷积计算精确定位电力设备位置的方法,该方法步骤如下:Fig. 1 is a flow chart of the present invention. A method for accurately locating the position of electric equipment through interpolation and convolution calculation, the steps of the method are as follows:
选取实际运行时的电流波数据x(采样率可以根据设备实际情况选择,若设备采样率越高则测量精度越高。本实施例下面列举出的数据采样率为100Ms/s,一组数据40个点)。用Matlab画出第一组数据二维图,如图2。图2是电流波形数据时域图。图2为我们对故障现场连续波形进行的采样,采样数据一共40个点。横坐标为点数,纵坐标为每一点的电流幅值。Select the current wave data x during actual operation (the sampling rate can be selected according to the actual situation of the equipment. If the sampling rate of the equipment is higher, the measurement accuracy will be higher. The data sampling rate listed below in this embodiment is 100Ms/s, and a set of
(2)对两组数据进行卷积计算,求出相关性最大点,如图3(图3是运用matlab中的xcorr函数对第一组数据与第二组数据进行相关性计算后的计算结果图。横坐标是两组数据相对位置,0点表示两组数据完全对应,1点表示两组数据错开一个点;纵坐标为相关性计算值),卷积值最大点在第40个点,即两组数据在第40个点相关性最大。(2) Carry out convolution calculations on two sets of data to find the maximum correlation point, as shown in Figure 3 (Figure 3 is the calculation result after using the xcorr function in matlab to calculate the correlation between the first set of data and the second set of data Fig. The abscissa is the relative position of the two sets of data, 0 points indicate that the two sets of data are completely corresponding, 1 point indicates that the two sets of data are staggered by one point; the ordinate is the correlation calculation value), the maximum convolution value is at the 40th point, That is, the two sets of data have the greatest correlation at the 40th point.
(3)对两组数据进行插值(每两个点之间插入4个点,本实施例为方便只插入4个点,实际使用时插入点数根据实际精度需要情况选择。插入点数越多则测量精度越高。),插值后的第一组数据如图4,图4是第一组数据插值后数据图(横坐标为点数,纵坐标为每一点的电流幅值),(3) Interpolate two sets of data (insert 4 points between each two points, this embodiment only inserts 4 points for convenience, the number of insertion points is selected according to the actual accuracy requirements in actual use. The more the number of insertion points, the more the measurement The higher the accuracy.), the first set of data after interpolation is shown in Figure 4, and Figure 4 is the data map after the first set of data interpolation (the abscissa is the number of points, and the ordinate is the current amplitude of each point),
第一组First group
第二组Second Group
(4)对插值后的两组数据再次进行卷积计算,得到结果如图5,图5插值后两组数据相关性计算结果图,在第201点卷积值最大。(4) Carry out convolution calculation again on the two sets of data after interpolation, and the result is shown in Figure 5. Figure 5 shows the correlation calculation results of the two sets of data after interpolation, and the convolution value is the largest at the 201st point.
从计算结果可知,通过合理的采样率,可以准确表达波形之间的相关性,同时又大大降低了采样率。It can be seen from the calculation results that the correlation between waveforms can be accurately expressed through a reasonable sampling rate, and the sampling rate is greatly reduced at the same time.
这样通过软件方法提高采样率,可以使我们的定位在不提高硬件条件的情况下更加精确。In this way, increasing the sampling rate through software can make our positioning more accurate without improving hardware conditions.
在第二步中,我们的采样率是100Ms/s,若故障点的位置在第40个点和第41个点之间,这样故障范围是3m。在重新插值卷积计算之后,故障点在第200个点与201个,故障范围是0.6m。这样就把采样率提高了5倍In the second step, our sampling rate is 100Ms/s, if the location of the fault point is between the 40th point and the 41st point, the fault range is 3m. After the re-interpolation convolution calculation, the fault point is at the 200th point and the 201st point, and the fault range is 0.6m. This increases the sampling rate by a factor of 5
本发明在不提高硬件采样率的情况下,通过计算方法,提高实际故障定位精度。The present invention improves the accuracy of actual fault location through a calculation method without increasing the hardware sampling rate.
Claims (2)
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 | |
Lin et al. | Travelling wave time–frequency characteristic-based fault location method for transmission lines | |
CN101718833B (en) | Transmission line fault single-ended location method based on natural frequency extraction of traveling waves | |
CN105705956B (en) | The method of one-end fault positioning in HVDC transmission circuit | |
CN102255291B (en) | A post-measurement simulation method for longitudinal protection of AC transmission lines based on Bergeron model | |
CN108896874B (en) | A transmission line fault location method with short lines connected at the end | |
CN102545177A (en) | Bergeron-model-based simulation-after-test method for fault phase selection of alternating current transmission line | |
CN103278747A (en) | High-tension transmission line single-ended traveling wave fault distance detection method combined with time-frequency characteristics | |
CN103245878B (en) | A kind of transmission line one-phase earth fault method of single end distance measurement | |
CN109643890A (en) | Method and system for being positioned to the failure in mixed type transmission line of electricity | |
CN105021950A (en) | Power line fault location method, device and system | |
CN102590704A (en) | After-test simulation method for internal and external failure recognition of double-circuit transmission line region based on Bergeron model | |
CN102998593B (en) | Method for accurately positioning electrical cable faults through interpolation convolution calculating | |
CN103941152B (en) | A kind of cable mixed line fault distance-finding method of k-NN algorithm based on waveform similarity | |
CN105403812B (en) | It is a kind of based on decomposing along fault traveling wave and the triangle loop grid one-end fault ranging method of distance calibration | |
CN109564256B (en) | Traveling wave based method for locating faults in transmission lines and apparatus therefor | |
CN106841918A (en) | A kind of cable low resistance faults localization method combined using single both-end | |
Subramani et al. | Fault investigation methods on power transmission line: A comparative study | |
CN103245887A (en) | Method utilizing lumped parameters to realize line interphase fault single-terminal location | |
CN103245890B (en) | Single-end distance measurement method for line single-phase ground faults against the influence of transition resistance and load current | |
CN105182186B (en) | A kind of radiation network Fault branch identification method based on voltage's distribiuting along the line and traveling wave information all standing | |
CN103217628B (en) | A kind of electric transmission line phase fault method of single end distance measurement | |
CN103267929B (en) | Voltage drop imaginary part characteristic is utilized to realize circuit inter-phase fault single-end ranging | |
CN103245876A (en) | Method for realizing line interphase fault single-ended distance measurement by means of real part characteristics of voltage drop | |
CN103217630B (en) | Method of achieving single-phase ground fault single-end distance measurement of line by means of voltage drop real part characteristics |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150415 Termination date: 20191115 |