CN107153155A - 一种电缆局部放电信号特征向量提取方法 - Google Patents
一种电缆局部放电信号特征向量提取方法 Download PDFInfo
- Publication number
- CN107153155A CN107153155A CN201710494855.8A CN201710494855A CN107153155A CN 107153155 A CN107153155 A CN 107153155A CN 201710494855 A CN201710494855 A CN 201710494855A CN 107153155 A CN107153155 A CN 107153155A
- Authority
- CN
- China
- Prior art keywords
- discharge signal
- imf
- local discharge
- cable local
- cable
- 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/12—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing
- G01R31/1227—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing of components, parts or materials
- G01R31/1263—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing of components, parts or materials of solid or fluid materials, e.g. insulation films, bulk material; of semiconductors or LV electronic components or parts; of cable, line or wire insulation
- G01R31/1272—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing of components, parts or materials of solid or fluid materials, e.g. insulation films, bulk material; of semiconductors or LV electronic components or parts; of cable, line or wire insulation of cable, line or wire insulation, e.g. using partial discharge measurements
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Relating To Insulation (AREA)
Abstract
本发明提供一种电缆局部放电信号特征向量提取方法。一种电缆局部放电信号特征向量提取方法,对已知来源的电缆局部放电信号分别进行小波包分解和经验模态分解,并分别构建相应的向量,最后将得到的两个向量合并组成电缆局部放电信号的特征向量。本发明提供的方法结合了小波包分解和经验模态分解的优势,提取的特征向量能够更好表征电缆局部放电信号。
Description
技术领域
本发明涉及电缆局部放电在线监测技术领域,更具体的涉及一种电缆局部放电信号特征向量提取方法。
背景技术
电缆局部放电在线监测中,检测到的局部放电信号可能来自电缆本体和电缆终端头,也可能来自与其相连的开关柜。由于不同来源的局部放电对设备危害不同,判断标准也不同,所以对局部放电信号来源进行识别有重要的现实意义。
在局部放电信号识别方面,信号特征提取和分类器选择是最关键部分。特征提取是局部放电信号识别第一步,特征提取的好坏直接影响到识别的效果。在特征提取方法中,小波分析法由于具有良好的局部时频特性,已广泛应用于局部放电信号的特征提取,但其容易收到噪声的影响,且基小波难于选择。而经验模式分解(empirical modedecomposition,EMD)方法基于信号的局部特征时间尺度,能把非平稳信号分解为有限的基本模态分量(intrinsic modefunction,IMF)之和,是一种自适应的信号处理方法,适合非线性和非平稳过程。
发明内容
本发明的目的是提供一种电缆局部放电信号特征向量提取方法,该方法结合了小波包分解和经验模式分解的优势,提取的特征向量能够更好表征局部放电信号。
为解决上述技术问题,本发明采用的技术方案是:一种电缆局部放电信号特征向量提取方法,其中,包括以下步骤:
(1)获取已知来源的电缆局部放电信号;
(2)对已知来源的各个电缆局部放电信号进行经验模态分解,得到各个固有模态函数分量,计算各个固有模态函数分量的对数能量熵,并构建相应的向量;
(3)对已知来源的各个电缆局部放电信号进行小波包分解,得到各个子频带上的小波包系数,计算各个子频带上小波包系数的对数能量熵,并构建相应的向量;
(4)将步骤(2)和步骤(3)得到的向量合并,组成电缆局部放电信号的特征向量。
本发明中,所述已知来源的电缆局部放电信号包括电缆本体局部放电信号、电缆终端头局部放电信号、开关柜的电晕放电信号和开关柜的表面放电信号。
进一步的,所述对已知来源的各个电缆局部放电信号进行经验模态分解,得到各个固有模态函数分量,计算各个固有模态函数分量的对数能量熵,并构建相应的向量的步骤为:
A.电缆局部放电信号进行经验模态分解步骤如下:
1)对于给定的一个电缆局部放电信号x(t),求出其所有局部极大值和局部极小值;
2)对这些极值点进行三次样条插值,得到由所有局部极大值点构成的上包络线和所有局部极小值点构成的下包络线,分别记为emax(t)和emin(t);并计算上、下包络线的均值m(t)=(emax(t)+emin(t))/2;
3)提取细节d(t)=x(t)-m(t),如果d(t)满足在整个数据序列中,极值点的数量与过零点的数量相等或者至多相差1个,且关于时间轴局部对称,则d(t)为第一个固有模态函数(IMF)分量;否则,以d(t)为输入重复步骤1)至步骤2),直到得到第一个IMF分量,记为IMF1(t);
4)记余项r1(t)=x(t)-IMF1(t),并作为新的新的待分析信号重复步骤1)至步骤3),以得到第二个IMF分量,记为IMF2(t),此时,余项r2(t)=r1(t)-IMF2(t);重复上述步骤,直到得到的余项rn(t)是一个单调信号或其值小于某个预先给定的阈值时,分解结束;
5)最后,得到n个IMF分量IMF1(t),IMF2(t)…IMFn(t)及余项rn(t),于是电缆局部放电信号x(t)可以表示为:
B.对于每一个IMF分量,按下式计算其对数能量熵:
式中,IMFi为第i个IMF分量,T为IMFi的长度,即电缆局部放电信号x(t)的长度,WEE(IMFi)为第i个IMF分量IMFi的对数能量熵;
C.根据计算得到的各个IMF分量对数能量熵,构建相应的向量λ1:
λ1=[WEE(IMF1),WEE(IMF2),.....,WEE(IMFn)]
进一步的,所述对已知来源的各个电缆局部放电信号进行小波包分解,得到各个子频带上的小波包系数;计算各个子频带上小波包系数的对数能量熵,并构建相应的向量的步骤为:
A.对电缆局部放电信号进行l层小波包分解,得到2l个子频带上的小波包系数,记为W1(k),W2(k)….WJ(k),J=2l;
B.对于每一个子频带上的小波包系数,按下式计算其对数能量熵:
式中,Wj为第j个子频带上的小波包系数,K为每个子频带上小波包系数的长度,WEE(Wj)为第j个子频带上小波包系数Wj的对数能量熵;
C.根据计算得到的各个子频带上小波包系数的对数能量熵,构建相应的向量λ2:
λ2=[WEE(W1),WEE(W2),.....,WEE(WJ)]
进一步的,所述将步骤(2)和步骤(3)得到的向量合并,组成电缆局部放电信号的特征向量的步骤为:
λ=[b1λ1,b2λ2]
式中,λ为电缆局部放电信号的特征向量;b1,b2为向量系数。
本发明与现有技术相比,其有益效果是:
本发明提供的一种电缆局部放电信号特征向量提取方法,对已知来源的电缆局部放电信号分别进行小波包分解和经验模态分解,并分别构建相应的向量,最后将得到的两个向量合并组成电缆局部放电信号的特征向量,该方法结合了小波包分解和经验模态分解的优势,提取的特征向量能够更好表征电缆局部放电信号。
附图说明
图1为本发明的一种电缆局部放电信号特征向量提取方法的流程示意图。
图2为电缆本体局部放电信号波形图。
图3为电缆终端头局部放电信号波形图。
图4为开关柜电晕放电信号波形图。
图5为开关柜表面放电信号波形图。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;为了更好说明本实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;对于本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。附图中描述位置关系仅用于示例性说明,不能理解为对本专利的限制。
如图1所示,一种电缆局部放电信号特征向量提取方法,其特征在于,包括以下步骤:
(1)获取已知来源的电缆局部放电信号。本实施例中,所述已知来源的局部放电信号包括电缆本体局部放电信号、电缆终端头局部放电信号、开关柜中的电晕放电信号和开关柜中的表面放电信号。如图2到图5所示的波形图,波形的采样频率为100MHz,每个波形的时域长度为1500个采样点。
(2)对已知来源的各个电缆局部放电信号进行经验模态分解,得到各个固有模态函数分量;计算各个固有模态函数分量的对数能量熵,并构建相应的向量。
具体地,包括如下步骤:
A.对电缆局部放电信号进行经验模态分解,步骤如下:
1)对于给定的一个电缆局部放电信号x(t),求出其所有局部极大值和局部极小值;
2)对这些极值点进行三次样条插值,得到由所有局部极大值点构成的上包络线和所有局部极小值点构成的下包络线,分别记为emax(t)和emin(t);并计算上、下包络线的均值m(t)=(emax(t)+emin(t))/2;
3)提取细节d(t)=x(t)-m(t),如果d(t)满足在整个数据序列中,极值点的数量与过零点的数量相等或者至多相差1个,且关于时间轴局部对称,则d(t)为第一个固有模态函数(IMF)分量;否则,以d(t)为输入重复步骤1)至步骤2),直到得到第一个IMF分量,记为IMF1(t);
4)记余项r1(t)=x(t)-IMF1(t),并作为新的新的待分析信号重复步骤1)至步骤3),以得到第二个IMF分量,记为IMF2(t),此时,余项r2(t)=r1(t)-IMF2(t);重复上述步骤,直到得到的余项rn(t)是一个单调信号或其值小于某个预先给定的阈值时,分解结束。
5)最后,得到n个IMF分量IMF1(t),IMF2(t)…IMFn(t)及余项rn(t),于是电缆局部放电信号x(t)可以表示为:
B.对于每一个IMF分量,按下式计算其对数能量熵:
式中,IMFi为第i个IMF分量,T为IMFi的长度,即电缆局部放电信号x(t)的长度,WEE(IMFi)为第i个IMF分量IMFi的对数能量熵;
C.根据计算得到的各个IMF分量对数能量熵,构建相应的向量λ1:
λ1=[WEE(IMF1),WEE(IMF2),.....,WEE(IMFn)]
(3)对已知来源的各个电缆局部放电信号进行小波包分解,得到各个子频带上的小波包系数;计算各个子频带上小波包系数的对数能量熵,并构建相应的向量。
具体地,包括如下步骤:
A.对电缆局部放电信号进行3层小波包分解,得到8个子频带上的小波包系数,记为W1(k),W2(k)….WJ(k),J=8;
B.对于每一个子频带上的小波包系数,按下式计算其对数能量熵:
式中,Wj为第j个子频带上的小波包系数,K为每个子频带上小波包系数的长度,WEE(Wj)为第j个子频带上小波包系数Wj的对数能量熵;
C.根据计算得到的各个子频带上小波包系数的对数能量熵,构建相应的向量λ2:
λ2=[WEE(W1),WEE(W2),.....,WEE(WJ)]
(4)将步骤(2)和步骤(3)得到的向量合并,组成电缆局部放电信号的特征向量。
具体地,包括如下步骤:
λ=[b1λ1,b2λ2]
式中,λ为电缆局部放电信号的特征向量;b1,b2为向量系数。
显然,本发明的上述实施例仅仅是为了清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。
Claims (5)
1.一种电缆局部放电信号特征向量提取方法,其特征在于,包括以下步骤:
(1)获取已知来源的电缆局部放电信号;
(2)对已知来源的各个电缆局部放电信号进行经验模态分解,得到各个固有模态函数分量,计算各个固有模态函数分量的对数能量熵,并构建相应的向量;
(3)对已知来源的各个电缆局部放电信号进行小波包分解,得到各个子频带上的小波包系数,计算各个子频带上小波包系数的对数能量熵,并构建相应的向量;
(4)将步骤(2)和步骤(3)得到的向量合并,组成电缆局部放电信号的特征向量。
2.根据权利要求1所述的一种电缆局部放电信号特征向量提取方法,其特征在于,所述已知来源的电缆局部放电信号包括电缆本体局部放电信号、电缆终端头局部放电信号、开关柜的电晕放电信号和开关柜的表面放电信号。
3.根据权利要求1所述的一种电缆局部放电信号特征向量提取方法,其特征在于,所述对已知来源的各个电缆局部放电信号进行经验模态分解,得到各个固有模态函数分量,计算各个固有模态函数分量的对数能量熵,并构建相应的向量的步骤为:
A.电缆局部放电信号进行经验模态分解步骤如下:
1)对于给定的一个电缆局部放电信号x(t),求出其所有局部极大值和局部极小值;
2)对这些极值点进行三次样条插值,得到由所有局部极大值点构成的上包络线和所有局部极小值点构成的下包络线,分别记为emax(t)和emin(t);并计算上、下包络线的均值m(t)=(emax(t)+emin(t))/2;
3)提取细节d(t)=x(t)-m(t),如果d(t)满足在整个数据序列中,极值点的数量与过零点的数量相等或者至多相差1个,且关于时间轴局部对称,则d(t)为第一个固有模态函数(IMF)分量;否则,以d(t)为输入重复步骤1)至步骤2),直到得到第一个IMF分量,记为IMF1(t);
4)记余项r1(t)=x(t)-IMF1(t),并作为新的新的待分析信号重复步骤1)至步骤3),以得到第二个IMF分量,记为IMF2(t),此时,余项r2(t)=r1(t)-IMF2(t);重复上述步骤,直到得到的余项rn(t)是一个单调信号或其值小于某个预先给定的阈值时,分解结束;
5)最后,得到n个IMF分量IMF1(t),IMF2(t)…IMFn(t)及余项rn(t),于是电缆局部放电信号x(t)可以表示为:
B.对于每一个IMF分量,按下式计算其对数能量熵:
式中,IMFi为第i个IMF分量,T为IMFi的长度,即电缆局部放电信号x(t)的长度,WEE(IMFi)为第i个IMF分量IMFi的对数能量熵;
C.根据计算得到的各个IMF分量对数能量熵,构建相应的向量λ1:
λ1=[WEE(IMF1),WEE(IMF2),.....,WEE(IMFn)] 。
4.根据权利要求1所述的一种电缆局部放电信号特征向量提取方法,其特征在于,所述对已知来源的各个电缆局部放电信号进行小波包分解,得到各个子频带上的小波包系数;计算各个子频带上小波包系数的对数能量熵,并构建相应的向量的步骤为:
A.对电缆局部放电信号进行l层小波包分解,得到2l个子频带上的小波包系数,记为W1(k),W2(k)….WJ(k),J=2l;
B.对于每一个子频带上的小波包系数,按下式计算其对数能量熵:
式中,Wj为第j个子频带上的小波包系数,K为每个子频带上小波包系数的长度,WEE(Wj)为第j个子频带上小波包系数Wj的对数能量熵;
C.根据计算得到的各个子频带上小波包系数的对数能量熵,构建相应的向量λ2:
λ2=[WEE(W1),WEE(W2),.....,WEE(WJ)] 。
5.根据权利要求1所述的一种电缆局部放电信号特征向量提取方法,其特征在于,所述将步骤(2)和步骤(3)得到的向量合并,组成电缆局部放电信号的特征向量的步骤为:
λ=[b1λ1,b2λ2]
式中,λ为电缆局部放电信号的特征向量;b1,b2为向量系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710494855.8A CN107153155A (zh) | 2017-06-26 | 2017-06-26 | 一种电缆局部放电信号特征向量提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710494855.8A CN107153155A (zh) | 2017-06-26 | 2017-06-26 | 一种电缆局部放电信号特征向量提取方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107153155A true CN107153155A (zh) | 2017-09-12 |
Family
ID=59796424
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710494855.8A Pending CN107153155A (zh) | 2017-06-26 | 2017-06-26 | 一种电缆局部放电信号特征向量提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107153155A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107991590A (zh) * | 2017-11-28 | 2018-05-04 | 广东电网有限责任公司珠海供电局 | 一种基于频域功率谱的电缆局部放电信号特征向量提取方法 |
CN109324247A (zh) * | 2018-11-13 | 2019-02-12 | 广东电网有限责任公司 | 一种电能质量扰动信号特征向量提取方法 |
CN109324250A (zh) * | 2018-11-29 | 2019-02-12 | 广东电网有限责任公司 | 一种电能质量扰动信号识别方法 |
CN109858490A (zh) * | 2018-12-21 | 2019-06-07 | 广东电网有限责任公司 | 一种电气设备红外图像特征向量提取方法 |
CN110175508A (zh) * | 2019-04-09 | 2019-08-27 | 杭州电子科技大学 | 一种应用于超声波局放检测的特征值提取方法 |
CN111337798A (zh) * | 2020-03-05 | 2020-06-26 | 国网四川电力送变电建设有限公司 | 一种特高压换流变压器绝缘监测及局部放电故障诊断方法 |
CN111239554B (zh) * | 2019-11-29 | 2021-04-13 | 深圳供电局有限公司 | 一种基于大数据的超声波局放检测分析模型 |
CN112905958A (zh) * | 2021-01-27 | 2021-06-04 | 南京国电南自电网自动化有限公司 | 基于测控装置的短时数据窗遥测数据状态辨识方法及系统 |
-
2017
- 2017-06-26 CN CN201710494855.8A patent/CN107153155A/zh active Pending
Non-Patent Citations (2)
Title |
---|
孙抗等: ""基于小波包和EMD的断路器机械故障诊断研究"", 《测控技术》 * |
张沫: "基于小波包分析的含噪序列主元特征提取和描述", 《中国优秀硕士学位论文全文数据库(信息科技辑)》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107991590A (zh) * | 2017-11-28 | 2018-05-04 | 广东电网有限责任公司珠海供电局 | 一种基于频域功率谱的电缆局部放电信号特征向量提取方法 |
CN109324247A (zh) * | 2018-11-13 | 2019-02-12 | 广东电网有限责任公司 | 一种电能质量扰动信号特征向量提取方法 |
CN109324247B (zh) * | 2018-11-13 | 2020-02-04 | 广东电网有限责任公司 | 一种电能质量扰动信号特征向量提取方法 |
CN109324250A (zh) * | 2018-11-29 | 2019-02-12 | 广东电网有限责任公司 | 一种电能质量扰动信号识别方法 |
CN109858490A (zh) * | 2018-12-21 | 2019-06-07 | 广东电网有限责任公司 | 一种电气设备红外图像特征向量提取方法 |
CN110175508A (zh) * | 2019-04-09 | 2019-08-27 | 杭州电子科技大学 | 一种应用于超声波局放检测的特征值提取方法 |
CN111239554B (zh) * | 2019-11-29 | 2021-04-13 | 深圳供电局有限公司 | 一种基于大数据的超声波局放检测分析模型 |
CN111337798A (zh) * | 2020-03-05 | 2020-06-26 | 国网四川电力送变电建设有限公司 | 一种特高压换流变压器绝缘监测及局部放电故障诊断方法 |
CN111337798B (zh) * | 2020-03-05 | 2022-04-08 | 国网四川电力送变电建设有限公司 | 一种特高压换流变压器绝缘监测及局部放电故障诊断方法 |
CN112905958A (zh) * | 2021-01-27 | 2021-06-04 | 南京国电南自电网自动化有限公司 | 基于测控装置的短时数据窗遥测数据状态辨识方法及系统 |
CN112905958B (zh) * | 2021-01-27 | 2024-04-19 | 南京国电南自电网自动化有限公司 | 基于测控装置的短时数据窗遥测数据状态辨识方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107153155A (zh) | 一种电缆局部放电信号特征向量提取方法 | |
CN107305774B (zh) | 语音检测方法和装置 | |
CN110175508A (zh) | 一种应用于超声波局放检测的特征值提取方法 | |
CN102314884B (zh) | 语音激活检测方法与装置 | |
KR20140059754A (ko) | 스펙트럼 모션 변환을 구현하는 사운드 신호 처리 시스템 및 방법 | |
CN109101910A (zh) | 一种基于噪声甄别的大地电磁信号去噪方法 | |
CN107274911A (zh) | 一种基于声音特征的相似度分析方法 | |
CN106646205A (zh) | 一种声振联合分析断路器故障的随机大扰动信号剔除算法 | |
CN107561420A (zh) | 一种基于经验模态分解的电缆局部放电信号特征向量提取方法 | |
Manjula et al. | Comparison of empirical mode decomposition and wavelet based classification of power quality events | |
CN105845149A (zh) | 声音信号中主音高的获取方法及系统 | |
CN109541455A (zh) | 一种基于s变换时频谱svd降噪的oltc冲击特性提取方法 | |
CN109598175A (zh) | 一种基于多小波基函数和超正交前向回归的时频分析方法 | |
CN105893976B (zh) | 一种行波信号的参数化识别方法 | |
CN107037338A (zh) | 一种gis振荡冲击耐压试验用缺陷类型识别方法 | |
CN106033669A (zh) | 语音识别方法及装置 | |
CN107210029A (zh) | 用于处理一连串信号以进行复调音符辨识的方法和装置 | |
CN204117590U (zh) | 语音采集降噪装置以及语音质量评价系统 | |
CN110367986A (zh) | 基于ewt的呼吸信号近似熵特征提取方法 | |
Lim et al. | Classification of underwater transient signals using MFCC feature vector | |
CN106291293A (zh) | 一种基于谱峭度与s变换的局放信号自适应去噪方法 | |
CN107430850A (zh) | 确定谐波信号的特征 | |
CN110648685B (zh) | 设备检测方法、装置、电子设备和可读存储介质 | |
CN110136741A (zh) | 一种基于多尺度上下文的单通道语音增强方法 | |
CN106340310B (zh) | 语音检测方法及装置 |
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: 20170912 |
|
RJ01 | Rejection of invention patent application after publication |