CN105467270B - 基于频谱相似性评价的单端行波测距反射波辨识算法 - Google Patents
基于频谱相似性评价的单端行波测距反射波辨识算法 Download PDFInfo
- Publication number
- CN105467270B CN105467270B CN201510811828.XA CN201510811828A CN105467270B CN 105467270 B CN105467270 B CN 105467270B CN 201510811828 A CN201510811828 A CN 201510811828A CN 105467270 B CN105467270 B CN 105467270B
- Authority
- CN
- China
- Prior art keywords
- wave
- traveling wave
- wave head
- frequency spectrum
- traveling
- 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.)
- Active
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 25
- 238000011156 evaluation Methods 0.000 title claims abstract description 13
- 238000005259 measurement Methods 0.000 claims abstract description 20
- 238000001514 detection method Methods 0.000 claims description 16
- 238000004088 simulation Methods 0.000 claims description 7
- 238000009499 grossing Methods 0.000 claims description 5
- 230000001960 triggered effect Effects 0.000 claims description 5
- 238000009795 derivation Methods 0.000 claims description 4
- 230000003595 spectral effect Effects 0.000 claims 2
- 238000000034 method Methods 0.000 abstract description 9
- 230000008901 benefit Effects 0.000 abstract description 6
- 230000005540 biological transmission Effects 0.000 abstract description 5
- 230000001052 transient effect Effects 0.000 abstract description 4
- 230000004807 localization Effects 0.000 abstract 1
- 230000007257 malfunction Effects 0.000 abstract 1
- 238000010183 spectrum analysis Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Locating Faults (AREA)
Abstract
本发明提出了一种基于频谱相似性评价的单端行波测距反射波辨识算法，根据行波的高频特性，截取适当的行波特征时间窗，构造了充分考虑高频分布特性的波头频谱相似函数，并运用频谱分析方法来分析故障暂态行波，并通过相似频谱的优选方法，进一步地提高了波头识别的工作效率与准确率。通过本发明提出的方法寻找与初始行波频谱相近的行波来确定故障第一个反射行波波头到达时刻，能快速准确完成故障测距，并且避免了直接人工确定波头，使波头到达时刻更易确定。本发明属于高压线路中的继电保护领域，具有能提高输电线路故障单端测距的准确率的优点。
Description
技术领域
本发明涉及输电技术领域，特别是涉及一种基于频谱相似性评价的单端行波测距反射波辨识算法。
背景技术
输电线路在发生线路故障时，突发的电压电流变化会以电磁波的形式沿着导线向两侧传播，其不但有低频暂态分量，而且还有包含着丰富故障信息的高频暂态行波，因此这些高频暂态行波可以用来进行故障和接地特征分析。当输电线路发生故障时，可以利用这些故障行波信息来确定故障位置，从而使故障能够得到快速处理，因此准确的测距方法对迅速排除故障和及时恢复供电具有十分重要的意义。
现代行波故障测距常见的基本原理有A型、D型、E型、F型四种，其中A型为单端测距，只需要在线路一端安装装置，具有节省大量成本的优点。但同时也存在波形分析困难，可靠性差，存在定位死区等问题。
发明内容
为解决克服A型单端测距的问题，本发明提供一种基于频谱相似性评价的单端行波测距反射波辨识算法，以有效提高单端行波测距时对于反射波的识别率，从而大大提高单端行波测距的成功率。
为了解决上述技术问题，本发明的技术方案如下：
基于频谱相似性评价的单端行波测距反射波辨识算法，包括以下步骤：
(1)在线路发生故障时触发行波检测装置后，记录和保存故障行波数据与时钟数据，并对所述数据进行平滑处理；
(2)对平滑后的数据进行求导，根据数据的极值点得出波头峰值；
(3)取所述波头峰值所在区域的一段行波作为该行波的特征时间窗；
(4)对所述特征时间窗的行波段进行傅立叶变换，得到所述行波检测装置接收到的各个波头信号对应的频域分布；
(5)由所述频域分布算出每个波头信号中每一个频率幅值占总幅值的比值，进而构建各个波头信号与初始波头信号的频谱相似度函数；
(6)利用所述频谱相似度函数，判断各个波头信号与初始波头信号的相似程度，找出线路故障的第一次反射信号；
(7)利用单端行波测距原理，求得故障点到测站的距离。
优选地，所述的波头峰值所在区域的一段行波是指包含了波头峰值且幅值达到波头峰值的0.3倍的行波段。
进一步地，所述步骤(5)包括：
利用依次计算各个波头信号中每一个频率幅值占总幅值的比值，其中，Xm(k)表示行波检测装置接收到的第m个波头信号的第k个数据点的频率幅值，表示行波检测装置接收到的第m个波头信号的总幅值，m＝0,1,2,3,...,表示行波检测装置接收到的波头信号的序号，k＝0,1,...N-1，N为特征时间窗的行波段上采集到的数据点的数量；
设不同波头信号与初始波头信号的相似程度为zm，令
其中λ为大于1的仿真参数；
根据仿真结果得出λ的合适取值，并代入上式获得各个波头信号与初始波头信号的频谱相似度函数。
进一步地，所述步骤(6)包括：
算出各个波头信号对应的相似程度zm；
判断行波检测装置接收到的波头信号的数量是否大于等于预设数量，若是，则从幅值大于预设阈值的波头信号中选取相似程度zm最大的波头信号为所述的第一次反射信号，若否，则直接选取相似程度zm最大的波头信号为所述的第一次反射信号。
本发明的优点是：通过对原始波形进行平滑后，利用快速傅立叶变换算法，计算得到了平滑波形的频域分布，通过对比波形之间的频域分布，能够更准确地寻找出第一次反射波形，从而能够大大提高单端测距的准确率，既很好地克服了原有的单端测距的缺点，同时能够完好地保留其原有的优点。
附图说明
图1是本发明辨识算法的流程图；
图2是获取行波特征时间窗的示意图；
图3是图2所示的行波所对应的频域分布；
图4是在一次线路故障中寻找第一次反射波的示意图。
具体实施方式
本发明提出了一种基于连续频谱相似性评价的单端行波测距反射波辨识算法，当配置了单端行波测距装置的线路上发生故障时，会产生故障行波，故障行波会在线路上经过多次的折射与反射，而单端行波测距装置也会接收到包括多次的折射、反射以及干扰在内多个波头。行波检测装置被触发记录保存故障行波数据与时钟数据。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
如图1至4所示，基于频谱相似性评价的单端行波测距反射波辨识算法，包括以下步骤：
(1)在线路发生故障时触发行波检测装置后，记录和保存故障行波数据与时钟数据，并对所述数据进行平滑处理；
(2)对平滑后的数据进行求导，根据数据的极值点得出波头峰值；
(3)取所述波头峰值所在区域的一段行波作为该行波的特征时间窗；
(4)对所述特征时间窗的行波段进行傅立叶变换，得到所述行波检测装置接收到的各个波头信号对应的频域分布；
(5)由所述频域分布算出每个波头信号中每一个频率幅值占总幅值的比值，进而构建各个波头信号与初始波头信号的频谱相似度函数；
(6)利用所述频谱相似度函数，判断各个波头信号与初始波头信号的相似程度，找出线路故障的第一次反射信号；
(7)利用单端行波测距原理，求得故障点到测站的距离。
优选地，以上所述步骤可以按下述具体实施方式展开。
1)线路发生故障时，产生的电压故障行波由故障点向两端传播，行波检测装置被触发记录保存故障行波数据与时钟数据。装置首先需要利用五点滑动平均法对采集到的数据进行平滑处理，以消除因噪音而产生的波形毛刺。即
其中，x′(n)为平滑后的行波数据，x(n)为平滑前的行波数据，N为装置对于这次故障行波所记录到的总数据点数。
2)利用平滑后的数据x′(n)对故障行波进行求导。求导公式如下：
其中，t为装置记录的采样间隔。
然后，根据f′(n)的数值变化，找出x′(n)的第一个极值点的大致时刻，并认为该时刻为波头到达时刻。
例如：若有f′(n)＝a，f′(n+1)＝-b，a,b≥0，则认为波头到达时刻t为
显然，这里计算的波头到达时刻t的准确度与采样间隔t有关，t越小，t的准确度越高，记初始波头到达时刻为t1。
3)取极值点前后两点中较大者作为行波的波头峰值xmax，并取峰值附近一段行波数据作为该行波的行波特征时间窗(如附图2所示)。按照仿真结果及运行经验，该行波特征时间窗一般取幅值达到0.3xmax的行波段。当然，以包含波头峰值的其它行波段作为行波特征时间窗也是可以的，本领域的技术人员可以根据实际情况选取，但根据运维和仿真的经验，选取幅值达到0.3xmax的行波段更能集中体现行波的特征，数据处理量也比较小。
4)将以上得到的波头特征时间窗进行傅立叶变换，有
其中，N为特征时间窗的行波段上采集到的数据点的总数。
以上公式可通过快速傅立叶变换(FFT)计算，得出装置接收到的每个波头的频率分布Xm(k)，m为接收到的行波信号的序号，m＝0,1,2,3,...。如附图3所示。其中X0(k)为线路发生故障后装置接收到的第一个行波信号，即初始波头信号，转换而来的频域分布。
5)将以上各次波头频域分布进行如下处理：
先利用下式计算X0(k)上每一个频率幅值占总幅值的比值，即
同理，可算出装置收到的每个波头信号的每一个频率幅值占总幅值的比值。
为衡量不同波头信号与初始波头信号的相似程度，优选地，可以使用如下方法：
设不同波头信号与初始波头信号的相似程度为zm，令
按照仿真结果，λ应大于1，一般可取1.02～1.07。
6)比较各个波头信号对应的相似程度zm值的大小。
若装置接收到的波头数据除初次波头外，总个数少于3个，则认为zm值最大的波头信号与原始波头信号最为相似。
若总个数不少于3个，则需要对此进行一次判定，由于第一次的反射波的幅值相对其他反射波的幅值，应该是相对较大的，因此应过滤掉幅值较小的一部分波头信号，仅考虑幅值较大的一部分相似程度zm值，认为其中相似程度zm值最大的波头信号与原始波头信号最为相似。即
{xmmax}＝{xMmax|xMmax>M0.5} (9)
其中，M0.5为原波头峰值集合{xMmax}的中位数，{xmmax}为新生成的波头峰值集合。
取此zm值对应的波头信号作为故障行波信号的第一次反射信号。并记录下该波头信号到达测站的精确时间，记为t2。
如附图4所示，在一次线路故障中，行波监测装置接收到原始故障行波及4个反射波，因接收到的波头超过3个，因此要进行判定，其中，x3max与x4max较小，因此不可能是我们需要的找的第一次反射波，将其过滤后，只需要分别对比x1max和x2max与初始行波的相似度即可。
7)利用单端行波测距原理可求得故障点到测站的距离，即
其中，L表示故障点到测站的距离；t1和t2分别表示故障产生电气分量的初始行波从第一次传到测量点的时间和从故障点反射回到测量点的时间；v为行波的传播速度。
本发明利用了故障行波特征时间窗的识别方法，很好地找出了每个行波的特征时间窗，并能够很好地适应各种类型的波形。同时，利用了频谱相似性评价方法，对不同的频段进行了不同的加权，增加了高频段的权重，使其更能够反映不同频谱之间的相似性。本发明还利用了相似频谱优选方法，对幅值过低的波头进行了过滤，大大减少了工作量及误差，大大提高了装置的工作效率。总之，通过本发明提出的方法寻找与初始行波频谱相近的行波来确定故障第一个反射行波波头到达时刻，能快速准确完成故障测距，其避免了直接人工确定波头，使波头到达时刻更易确定，具有能提高输电线路故障单端测距的准确率的优点。
上列详细说明是针对本发明可行实施例的具体说明，该实施例并非用以限制本发明的专利范围，凡未脱离本发明所为的等效实施或变更，均应包含于本案的专利范围中。
Claims (4)
1.基于频谱相似性评价的单端行波测距反射波辨识算法，其特征在于，包括以下步骤：
(1)在线路发生故障时触发行波检测装置后，记录和保存故障行波数据与时钟数据，并对所述数据进行平滑处理；
(2)对平滑后的数据进行求导，根据数据的极值点得出波头峰值；
(3)取所述波头峰值所在区域的一段行波作为该行波的特征时间窗；
(4)对所述特征时间窗的行波段进行傅立叶变换，得到所述行波检测装置接收到的各个波头信号对应的频域分布；
(5)由所述频域分布算出每个波头信号中每一个频率幅值占总幅值的比值，进而构建各个波头信号与初始波头信号的频谱相似度函数；
(6)利用所述频谱相似度函数，判断各个波头信号与初始波头信号的相似程度，找出线路故障的第一次反射信号；
(7)利用单端行波测距原理，求得故障点到测站的距离。
2.根据权利要求1所述的基于频谱相似性评价的单端行波测距反射波辨识算法，其特征在于，所述的波头峰值所在区域的一段行波是指包含了波头峰值且幅值达到波头峰值的0.3倍的行波段。
3.根据权利要求1或2所述的基于频谱相似性评价的单端行波测距反射波辨识算法，其特征在于，所述步骤(5)包括：
利用
依次计算各个波头信号中每一个频率幅值占总幅值的比值，其中，Xm(k)表示行波检测装置接收到的第m个波头信号的第k个数据点的频率幅值，表示行波检测装置接收到的第m个波头信号的总幅值，m＝0,1,2,3,...,表示行波检测装置接收到的波头信号的序号，N为特征时间窗的行波段上采集到的数据点的数量；
设不同波头信号与初始波头信号的相似程度为zm，令
<mrow>
<msub>
<mi>z</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msup>
<mi>&lambda;</mi>
<mi>k</mi>
</msup>
<mo>|</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>a</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
</mfrac>
<mo>,</mo>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中λ为大于1的仿真参数；
根据仿真结果得出λ的合适取值，并代入上式获得各个波头信号与初始波头信号的频谱相似度函数。
4.根据权利要求3所述的基于频谱相似性评价的单端行波测距反射波辨识算法，其特征在于，所述步骤(6)包括：
算出各个波头信号对应的相似程度zm；
判断行波检测装置接收到的波头信号的数量是否大于等于预设数量，若是，则从幅值大于预设阈值的波头信号中选取相似程度zm最大的波头信号为所述的第一次反射信号，若否，则直接选取相似程度zm最大的波头信号为所述的第一次反射信号。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510811828.XA CN105467270B (zh) | 2015-11-20 | 2015-11-20 | 基于频谱相似性评价的单端行波测距反射波辨识算法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510811828.XA CN105467270B (zh) | 2015-11-20 | 2015-11-20 | 基于频谱相似性评价的单端行波测距反射波辨识算法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105467270A CN105467270A (zh) | 2016-04-06 |
| CN105467270B true CN105467270B (zh) | 2018-05-11 |
Family
ID=55605186
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510811828.XA Active CN105467270B (zh) | 2015-11-20 | 2015-11-20 | 基于频谱相似性评价的单端行波测距反射波辨识算法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105467270B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112083271B (zh) * | 2020-08-18 | 2021-10-22 | 昆明理工大学 | 一种基于护层电流行波时频复合分析的10kV电缆在线测距方法 |
| CN113687192B (zh) * | 2021-10-26 | 2022-03-01 | 武汉华瑞伏安电力科技有限公司 | 输电线路放电信号采集及定位方法 |
| CN118275825B (zh) * | 2024-06-03 | 2024-09-10 | 东方电子股份有限公司 | 一种用于配电线路行波测距的数据处理方法及系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4095073B2 (ja) * | 2005-03-04 | 2008-06-04 | 関西電力株式会社 | 送電線の故障点標定方法および送電線の故障点標定装置および故障点標定プログラム |
| CN101666848A (zh) * | 2009-10-15 | 2010-03-10 | 北京四方继保自动化股份有限公司 | 基于行波的故障定位方法 |
| CN103217626A (zh) * | 2013-03-26 | 2013-07-24 | 昆明理工大学 | 一种利用正负极性波头时序间隔的单端行波故障测距方法 |
| CN104820165A (zh) * | 2015-04-20 | 2015-08-05 | 华南理工大学 | 一种基于单极信息的同塔双回直流输电线路故障选线方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2780998B1 (en) * | 2011-11-17 | 2016-05-11 | ABB Research Ltd. | Fault protection in mixed high-voltage transmission lines |
-
2015
- 2015-11-20 CN CN201510811828.XA patent/CN105467270B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4095073B2 (ja) * | 2005-03-04 | 2008-06-04 | 関西電力株式会社 | 送電線の故障点標定方法および送電線の故障点標定装置および故障点標定プログラム |
| CN101666848A (zh) * | 2009-10-15 | 2010-03-10 | 北京四方继保自动化股份有限公司 | 基于行波的故障定位方法 |
| CN103217626A (zh) * | 2013-03-26 | 2013-07-24 | 昆明理工大学 | 一种利用正负极性波头时序间隔的单端行波故障测距方法 |
| CN104820165A (zh) * | 2015-04-20 | 2015-08-05 | 华南理工大学 | 一种基于单极信息的同塔双回直流输电线路故障选线方法 |
Non-Patent Citations (1)
| Title |
|---|
| 基于导数法的故障行波波头识别改进算法;徐伟宗;《电网技术》;20100105;第34卷(第1期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105467270A (zh) | 2016-04-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108469560B (zh) | 一种基于快速s变换时频空间模型的电磁干扰客观复杂度评估方法 | |
| CN102680860B (zh) | 一种高压电力线路行波测距用故障点自动定位方法 | |
| CN104597376B (zh) | 一种考虑实测波速的高压直流输电线路故障测距方法 | |
| CN109375060B (zh) | 一种配电网故障波形相似度计算方法 | |
| CN103323741B (zh) | 一种针对强故障的基于故障电压初始行波幅值比较的d型线缆混合线路故障区段判别的方法 | |
| CN105467270B (zh) | 基于频谱相似性评价的单端行波测距反射波辨识算法 | |
| CN104155650A (zh) | 一种基于熵权值法点迹质量评估的目标跟踪方法 | |
| CN102323518A (zh) | 一种基于谱峭度的局部放电信号识别方法 | |
| CN106646205A (zh) | 一种声振联合分析断路器故障的随机大扰动信号剔除算法 | |
| CN107329046B (zh) | 基于模量分析的直流架空线雷击识别方法 | |
| CN106771598B (zh) | 一种自适应谱峭度信号处理方法 | |
| CN105223481A (zh) | 基于差值能量函数的局部放电特高频信号起始时刻确定方法 | |
| CN108344922B (zh) | 一种基于相类区分和s变换的输电线路直击雷故障识别方法 | |
| CN114217164B (zh) | 一种基于放电波形智能识别的电缆故障测距方法及系统 | |
| CN107179473B (zh) | 一种输电线路故障定位方法 | |
| WO2016004687A1 (zh) | 超高频局放信号初始时刻判别方法 | |
| CN108199356B (zh) | 基于波前信息的直流输电线路超高速保护方法 | |
| CN102621452A (zh) | 一种基于信号距离和π型线路模型的纵联保护方法 | |
| CN103490830B (zh) | 基于物联网电力测温设备的去噪声射频频谱峰值获取方法 | |
| CN105893976B (zh) | 一种行波信号的参数化识别方法 | |
| CN109991520A (zh) | 一种电缆振荡波局放检测系统用波速计算新方法 | |
| CN110095691B (zh) | 基于全波形主频分量的初始行波波头提取方法和装置 | |
| CN109596354B (zh) | 基于自适应共振频带识别的带通滤波方法 | |
| CN104965155B (zh) | 一种基于滤波特性的波头查找方法 | |
| CN110048741A (zh) | 一种基于短时分数阶傅里叶变换的跳频信号的参数估计方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |