CN110927543A - 一种电力设备局部放电超声信号时差估计方法 - Google Patents
一种电力设备局部放电超声信号时差估计方法 Download PDFInfo
- Publication number
- CN110927543A CN110927543A CN201911305793.7A CN201911305793A CN110927543A CN 110927543 A CN110927543 A CN 110927543A CN 201911305793 A CN201911305793 A CN 201911305793A CN 110927543 A CN110927543 A CN 110927543A
- Authority
- CN
- China
- Prior art keywords
- time difference
- signal
- partial discharge
- imf
- noise
- 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
Images
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/1209—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 using acoustic measurements
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Testing Relating To Insulation (AREA)
Abstract
本发明公开了一种电力设备局部放电超声信号时差估计方法,涉及局部放电检测领域,解决了传统局放定位信噪比低、估计精度差的弊端,其技术方案要点是输入两路原信号,经过CEEMDAN算法将每一路原信号分解为多个不同频率成分的固有模态分量IMF;通过相关系数法对各IMF进行自相关计算以辨识噪声分量与局部放电超声信号的分界点;根据计算辨识的分界点设定噪声阈值,并对小于噪声阈值的IMF进行小波阈值去噪;将去噪后获得IMF与分解获得的大于噪声阈值的IMF进行重构,得到去噪局部放电超声信号;通过相关系数波形比较法对两路去噪局部放电超声信号进行时差估计,本发明的一种电力设备局部放电超声信号时差估计方法,能够提高信噪比,从而提高时差估计精度。
Description
技术领域
本发明涉及局部放电检测,特别涉及一种电力设备局部放电超声信号时差估计方法。
背景技术
超声波定位技术作为一种成本低、不受电磁干扰、非破坏性的故障定位技术被逐步应用于局部放电检测领域,而对超声信号的时差估计是定位的关键。现有的应用于局放领域的时延估计方法有相关分析法、自适应时差估计法和相位谱估计法等。
目前产生局放的电力设备现场环境十分复杂,在检测中往往会遇到各种干扰,严重影响检测的灵敏度和可靠性。变电站现场检测的干扰信号通常包含白噪声、随机脉冲干扰和窄带干扰。
现有的基于超声局放定位的时延估计方法大多是峰值检测算法,但实际PD超声信号为窄带信号,相关峰会因为带宽窄而在时间域上被延展,使得峰值不明显进而导致其估计精度对噪声十分敏感且精度较低。
发明内容
本发明的目的是提供一种电力设备局部放电超声信号时差估计方法,能够提高信噪比,从而提高时差估计精度。
本发明的上述技术目的是通过以下技术方案得以实现的:
一种电力设备局部放电超声信号时差估计方法,包括有以下步骤:
输入两路原信号,经过CEEMDAN算法将每一路原信号都分解为多个不同频率成分的固有模态分量IMF;
通过相关系数法对各IMF进行自相关计算以辨识噪声分量与局部放电超声信号的分界点;
根据计算辨识的分界点设定噪声阈值,并对小于噪声阈值的IMF进行小波阈值去噪;
将去噪后获得IMF与分解获得的大于噪声阈值的IMF进行重构,得到去噪局部放电超声信号;
通过相关系数波形比较法对两路去噪局部放电超声信号进行时差估计。
作为优选,CEEMDAN算法如下:
对每一路原信号均进行加零均值单位方差白噪声处理得到第一信号;
对第一信号进行N次EMD分解,得到N个第一阶分解分量和剩余残差;
对N个第一阶分解分量进行求集成均值,得到第一阶分量,并计算得到第一残差;
同样对计算得到的第一残差再次进行加零均值单位方差白噪声处理,以进行N次EMD分解计算,重复m次以对计算得到的第m残差进行处理,直至最终获得的第m残差信号无法继续分解,得到第m阶分量与最终残差;得到由第一阶分量至第m阶分量与最终残差累加形式的原信号。
作为优选,CEEMDAN算法具体步骤如下:
在原信号s(t)中加入正负成对且幅值为a的零均值单位方差白噪声,得到第一信号s'(t)=s(t)+(-1)qani(t);
对第一信号s'(t)进行N次EMD分解,得到N个第一阶分解分量IMF1 i(t)和剩余残差r1 i(t),即s'(t)=IMF1 i(t)+r1 i(t);
作为优选,通过相关系数法对各固有模态分量IMF做自相关计算具体计算公式为:
作为优选,小波阈值去噪具体步骤为:
对需要去噪的固有模态分量IMF进行小波分解,并求取对应的小波系数dj,k;
设置噪声阈值λ,对固有模态分量IMF分解的小波中小于噪声阈值的小波系数对应的小波判断为由噪声产生并对应去除,重构处理后的小波系数dj,k与近似系数,还原去噪的IMF。
作为优选,噪声阈值λ的选取计算如下:
作为优选,根据相关系数波形比较法对两路原信号的时差估算具体:
两路原信号模型为
式中,D为两个信号的时差,n1(t)和n2(t)为相互独立的高斯白噪声;
由相关函数性质,可得到
将两路原信号的时差转移到至相关函数中;
根据估计值进行表示:
N为信号采样点个数,比较估计值的波形以获得时差。
作为优选,通过比较估计值的波形以获得时差具体为:
A3、重复步骤A2,对得到的多个估计值进行直方图统计,选取重复次数最多的数值作为最终得结果。
综上所述,本发明具有以下有益效果:
利用相关系数波形比较算法对PD超声窄带信号进行估计,避开传统峰值检测算法对于窄带信号峰值在时间域被展平而导致的峰值不明显估计精度差的问题,在时差估计前利用CEEMDAN阈值去噪对带噪信号进行去噪,提升实际采集PD超声信号的信噪比,从而提高时差估计精度,提高了相关系数法在电力现场复杂噪声情况下的鲁棒性和稳定性。
附图说明
图1为CEEMDAN阈值去噪算法流程图;
图2为基于CEEMDAN改进的波形比较法流程图。
具体实施方式
以下结合附图对本发明作进一步详细说明。
根据一个或多个实施例,公开的一种电力设备局部放电超声信号时差估计方法,如图1及图2所示,具体步骤如下:。
输入两路原信号,经过CEEMDAN算法将每一路原信号都分解为多个不同频率成分的固有模态分量IMF;
通过相关系数法对各IMF进行自相关计算以辨识噪声分量与局部放电超声信号的分界点;
根据计算辨识的分界点设定噪声阈值,并对小于噪声阈值的IMF进行小波阈值去噪;
将去噪后获得IMF与分解获得的大于噪声阈值的IMF进行重构,得到去噪局部放电超声信号;
通过相关系数波形比较法对两路去噪局部放电超声信号进行时差估计。
如图1所示,原信号s(t)由局部放电超声信号x(t)以及白噪声n(t)构成的混合信号经过CEEMDAN分解成多个固有模态分量IMF,利用相关系数法,确定噪声占主导的IMF并对其进行小波阈值去噪,随后将去噪后的IMF和剩余的IMF进行重构,得到去噪后的局部放电超声信号即PD超声信号。最后利用相关系数波形比较法得出时差估计值整体时差估计算法图如图2所示。
具体的,设两路信号模型为:
式中x(t)和x(t-D)为源信号,D为两个信号的时差。n1(t)和n2(t)为相互独立的高斯白噪声。
CEEMDAN是由EMD原理的基础上改进而来。其计算步骤如下:
(1)在原信号s(t)中加入正负成对的幅值为a的零均值单位方差白噪声,得到:
s'(t)=s(t)+(-1)qani(t)
s'(t)=IMF1 i(t)+r1 i(t)
得到第一残差:
同理最终第二残差:
相关法即对各IMF作自相关,自相关公式为:
小波阈值去噪步骤为去噪信号进行小波分解,并求取对应的小波系数dj,k,设置噪声阈值λ,对小于噪声阈值的系数认为是由噪声产生并将其去除,最后重构处理后的小波系数dj,k与近似系数,还原去噪信号。选取阈值函数,其表达式为:
相关系数波形比较法具体步骤为
N为信号采样点个数,则比较估计值波形即可获得时差。避免了窄带信号在时间域的延展性对峰值检测所造成的的影响。具体计算步骤为:
3)重复步骤2),对得到的多个估计值进行直方图统计,选取重复次数最多的数值作为最终得结果。
利用相关系数波形比较算法对PD超声窄带信号进行估计,解决了传统峰值检测时延估计算法如相关法等对窄带信号估计精度低的问题,在时差估计前利用CEEMDAN阈值去噪对带噪信号进行去噪,提升实际采集PD超声信号的信噪比,提高了相关系数法在电力现场复杂噪声情况下的鲁棒性和稳定性。
在传统局放定位采取的峰值检测时延算法未充分考虑PD超声信号窄带信号对于其估计精度的影响,采用对带宽不敏感的相关系数波形比较法对其进行时差估计,同时利用CEEMDAN阈值去噪对采集信号进行预处理,解决传统波形比较法在低信噪比条件下估计精度差的问题。采用的CEEMDAN阈值去噪具有完备性,能够很好地保留信号中的时差信息,而通过相关系数波形比较将信号中的时差信息转移到自相关和互相关函数中,进一步提高信噪比,从而提高时差估计精度,为后续的定位打下良好的基础。
本发明可以广泛适用于各种基于到达时差的电力设备局放源超声定位的时差估计,算法计算量较小,在低信噪比环境下鲁棒性更强,易于配合相关空间定位算法实现在线定位工作。
本具体实施例仅仅是对本发明的解释,其并不是对本发明的限制,本领域技术人员在阅读完本说明书后可以根据需要对本实施例做出没有创造性贡献的修改,但只要在本发明的权利要求范围内都受到专利法的保护。
Claims (8)
1.一种电力设备局部放电超声信号时差估计方法,其特征是,包括有以下步骤:
输入两路原信号,经过CEEMDAN算法将每一路原信号都分解为多个不同频率成分的固有模态分量IMF;
通过相关系数法对各IMF进行自相关计算以辨识噪声分量与局部放电超声信号的分界点;
根据计算辨识的分界点设定噪声阈值,并对小于噪声阈值的IMF进行小波阈值去噪;
将去噪后获得IMF与分解获得的大于噪声阈值的IMF进行重构,得到去噪局部放电超声信号;
通过相关系数波形比较法对两路去噪局部放电超声信号进行时差估计。
2.根据权利要求1所述的电力设备局部放电超声信号时差估计方法,其特征是,CEEMDAN算法如下:
对每一路原信号均进行加零均值单位方差白噪声处理得到第一信号;
对第一信号进行N次EMD分解,得到N个第一阶分解分量和剩余残差;
对N个第一阶分解分量进行求集成均值,得到第一阶分量,并计算得到第一残差;
同样对计算得到的第一残差再次进行加零均值单位方差白噪声处理,以进行N次EMD分解计算,重复m次以对计算得到的第m残差进行处理,直至最终获得的第m残差信号无法继续分解,得到第m阶分量与最终残差;得到由第一阶分量至第m阶分量与最终残差累加形式的原信号。
3.根据权利要求2所述的电力设备局部放电超声信号时差估计方法,其特征是,CEEMDAN算法具体步骤如下:
在原信号s(t)中加入正负成对且幅值为a的零均值单位方差白噪声,得到第一信号s'(t)=s(t)+(-1)qani(t);
对第一信号s'(t)进行N次EMD分解,得到N个第一阶分解分量IMF1 i(t)和剩余残差r1 i(t),即s'(t)=IMF1 i(t)+r1 i(t);
5.根据权利要求4所述的电力设备局部放电超声信号时差估计方法,其特征是,小波阈值去噪具体步骤为:
对需要去噪的固有模态分量IMF进行小波分解,并求取对应的小波系数dj,k;
设置噪声阈值λ,对固有模态分量IMF分解的小波中小于噪声阈值的小波系数对应的小波判断为由噪声产生并对应去除,重构处理后的小波系数dj,k与近似系数,还原去噪的IMF。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911305793.7A CN110927543A (zh) | 2019-12-18 | 2019-12-18 | 一种电力设备局部放电超声信号时差估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911305793.7A CN110927543A (zh) | 2019-12-18 | 2019-12-18 | 一种电力设备局部放电超声信号时差估计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110927543A true CN110927543A (zh) | 2020-03-27 |
Family
ID=69863014
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911305793.7A Pending CN110927543A (zh) | 2019-12-18 | 2019-12-18 | 一种电力设备局部放电超声信号时差估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110927543A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112327108A (zh) * | 2020-10-09 | 2021-02-05 | 国网冀北电力有限公司检修分公司 | 一种罐式断路器局部放电超声信号去噪及时间差识别方法 |
CN112731292A (zh) * | 2020-12-25 | 2021-04-30 | 吉林大学 | 局部imf能量加权的低空飞行目标信号时延估计方法 |
CN112924831A (zh) * | 2021-03-05 | 2021-06-08 | 国网山东省电力公司电力科学研究院 | 一种特高频局部放电定位时延估计方法 |
CN113925482A (zh) * | 2021-10-26 | 2022-01-14 | 西安易朴通讯技术有限公司 | 心率计算方法、可穿戴电子设备和存储介质 |
CN114264724A (zh) * | 2021-12-23 | 2022-04-01 | 大连海洋大学 | 基于ceemdan联合小波包阈值的超声信号处理方法 |
CN114355284A (zh) * | 2022-01-04 | 2022-04-15 | 电子科技大学 | 一种利用频谱主分量的时差估计方法 |
DE102021206719A1 (de) | 2021-06-29 | 2022-12-29 | Siemens Aktiengesellschaft | Ortung eines Lichtbogens |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104360251A (zh) * | 2014-12-02 | 2015-02-18 | 广西电网公司电力科学研究院 | 一种变压器局部放电的超声波信号时延估计方法 |
CN104614647A (zh) * | 2015-01-23 | 2015-05-13 | 云南电网有限责任公司电力科学研究院 | 一种复小波变换局放定位试验方法及装置 |
KR20180032962A (ko) * | 2016-09-23 | 2018-04-02 | 한국전력공사 | 부분방전 진단 장치 및 방법 |
CN107966206A (zh) * | 2017-11-07 | 2018-04-27 | 华北电力大学(保定) | 一种光纤超声阵列传感器的阵列信号去噪方法 |
CN110470956A (zh) * | 2019-08-05 | 2019-11-19 | 上海电机学院 | 一种电力设备局部放电超声定位方法 |
-
2019
- 2019-12-18 CN CN201911305793.7A patent/CN110927543A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104360251A (zh) * | 2014-12-02 | 2015-02-18 | 广西电网公司电力科学研究院 | 一种变压器局部放电的超声波信号时延估计方法 |
CN104614647A (zh) * | 2015-01-23 | 2015-05-13 | 云南电网有限责任公司电力科学研究院 | 一种复小波变换局放定位试验方法及装置 |
KR20180032962A (ko) * | 2016-09-23 | 2018-04-02 | 한국전력공사 | 부분방전 진단 장치 및 방법 |
CN107966206A (zh) * | 2017-11-07 | 2018-04-27 | 华北电力大学(保定) | 一种光纤超声阵列传感器的阵列信号去噪方法 |
CN110470956A (zh) * | 2019-08-05 | 2019-11-19 | 上海电机学院 | 一种电力设备局部放电超声定位方法 |
Non-Patent Citations (2)
Title |
---|
张建文等: ""基于CEEMDAN与小波自适应阈值的去噪新方法"", 《电测与仪表》 * |
胡锦熙 等: ""Arrival Time Difference Estimation of Ultrasonic Signals from Partial Discharge in Electric Power Equipments"", 《JOURNAL OF ELECTRICAL AND ELECTRONIC ENGINEERING》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112327108A (zh) * | 2020-10-09 | 2021-02-05 | 国网冀北电力有限公司检修分公司 | 一种罐式断路器局部放电超声信号去噪及时间差识别方法 |
CN112731292A (zh) * | 2020-12-25 | 2021-04-30 | 吉林大学 | 局部imf能量加权的低空飞行目标信号时延估计方法 |
CN112731292B (zh) * | 2020-12-25 | 2024-04-05 | 吉林大学 | 局部imf能量加权的低空飞行目标信号时延估计方法 |
CN112924831A (zh) * | 2021-03-05 | 2021-06-08 | 国网山东省电力公司电力科学研究院 | 一种特高频局部放电定位时延估计方法 |
DE102021206719A1 (de) | 2021-06-29 | 2022-12-29 | Siemens Aktiengesellschaft | Ortung eines Lichtbogens |
CN113925482A (zh) * | 2021-10-26 | 2022-01-14 | 西安易朴通讯技术有限公司 | 心率计算方法、可穿戴电子设备和存储介质 |
CN113925482B (zh) * | 2021-10-26 | 2024-05-31 | 西安易朴通讯技术有限公司 | 心率计算方法、可穿戴电子设备和存储介质 |
CN114264724A (zh) * | 2021-12-23 | 2022-04-01 | 大连海洋大学 | 基于ceemdan联合小波包阈值的超声信号处理方法 |
CN114355284A (zh) * | 2022-01-04 | 2022-04-15 | 电子科技大学 | 一种利用频谱主分量的时差估计方法 |
CN114355284B (zh) * | 2022-01-04 | 2023-05-05 | 电子科技大学 | 一种利用频谱主分量的时差估计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110927543A (zh) | 一种电力设备局部放电超声信号时差估计方法 | |
CN109243483B (zh) | 一种含噪频域卷积盲源分离方法 | |
CN106597408B (zh) | 基于时频分析和瞬时频率曲线拟合的高阶pps信号参数估计方法 | |
CN109932624B (zh) | 一种基于高斯尺度空间的电缆局放周期窄带干扰去噪方法 | |
CN117390434A (zh) | 一种基于多尺度分解与重构的瞬态电磁信号时延估计方法 | |
CN109923430B (zh) | 用于进行相位差展开的装置及方法 | |
CN112766224B (zh) | 畸变信号中提取真实信号的方法、装置、设备和存储介质 | |
Djurovic | Estimation of the sinusoidal signal frequency based on the marginal median DFT | |
CN106980722B (zh) | 一种脉冲响应中谐波成分的检测和去除方法 | |
JP2016500847A (ja) | デジタルプロセッサベースの複素音響共鳴デジタル音声分析システム | |
CN108020761B (zh) | 一种局部放电去噪方法 | |
CN113822363B (zh) | 一种基于二次有理核函数的卷积神经网络时延估计方法 | |
CN109272054B (zh) | 一种基于独立性的振动信号去噪方法及系统 | |
CN114813129A (zh) | 基于wpe与emd的滚动轴承声信号故障诊断方法 | |
Sircar et al. | Parametric modeling of speech by complex AM and FM signals | |
Pace et al. | Nyquist folding analog-to-information receiver: Autonomous information recovery using quadrature mirror filtering | |
CN112929053B (zh) | 一种跳频信号特征提取与参数估计方法 | |
Sumithra et al. | Wavelet based speech signal de-noising using hybrid thresholding | |
Abdullah et al. | MODELLING OF PARTIAL DISCHARGE ANALYSIS SYSTEM USING WAVELET TRANSFORM DENOISING TECHNIQUE IN LABVIEW ENVIRONMENT | |
CN116699337A (zh) | 一种基于时延估计的城市地下电缆局部放电定位方法 | |
Jiang et al. | Data Assessment Method Based on the Statistical Analysis | |
CN115526204A (zh) | 一种混合振荡信号参数辨识方法及系统 | |
Zezula et al. | Audio signal watermarking in mclt domain with the aid of 2d pattern | |
CN117607624A (zh) | 基于声信号的配用电设备故障检测方法 | |
CN114781460A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200327 |