CN105223482B - 局部放电特高频信号波形的小波分解二值去噪方法 - Google Patents
局部放电特高频信号波形的小波分解二值去噪方法 Download PDFInfo
- Publication number
- CN105223482B CN105223482B CN201510697143.7A CN201510697143A CN105223482B CN 105223482 B CN105223482 B CN 105223482B CN 201510697143 A CN201510697143 A CN 201510697143A CN 105223482 B CN105223482 B CN 105223482B
- Authority
- CN
- China
- Prior art keywords
- mrow
- scale
- signal
- wavelet
- partial discharge
- 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.)
- Expired - Fee Related
Links
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 77
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000001514 detection method Methods 0.000 claims abstract description 9
- 230000000877 morphologic effect Effects 0.000 claims abstract description 5
- 238000005070 sampling Methods 0.000 claims description 14
- 230000001186 cumulative effect Effects 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 230000007797 corrosion Effects 0.000 claims description 6
- 238000005260 corrosion Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 4
- 238000007619 statistical method Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 2
- 230000003628 erosive effect Effects 0.000 claims description 2
- 238000010606 normalization Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 6
- 230000006872 improvement Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000009413 insulation Methods 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000010339 dilation Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Testing Relating To Insulation (AREA)
Abstract
本发明局部放电特高频信号波形的小波分解二值去噪方法,包括:1)在被检测的电力设备上安装局部放电特高频传感器,利用检测设备采集局部放电特高频信号;2)根据选取的母小波与分解层数,对局部放电信号进行多尺度小波分解,得到各尺度小波系数与信号波形;3)以减小局部放电信号波形畸变为目标,对小波分解系数进行处理:计算各尺度信号累积能量函数的数学形态学梯度,并计算数学形态学梯度的陡峭度;将陡峭度与某一阈值进行比较,若陡峭度大于该阈值认为该尺度包含局部放电信号,该尺度小波系数保持不变;否则该尺度仅包含噪声信号,将其小波系数置为零;4)根据处理后的小波系数进行小波重构,得到去噪后的局部放电特高频信号。
Description
技术领域:
本发明属于电力设备绝缘状态评估技术领域,具体涉及一种局部放电特高频信号波形的小波分解二值去噪方法。
背景技术:
局部放电检测作为发现电力设备内隐藏缺陷及绝缘状态评估的有效手段,在变压器及气体组合绝缘电器等众多电力设备中获得广泛应用。而特高频局部放电检测方法具有灵敏度高、抗噪性好及可实现局部放电源定位等优点,逐渐成为检测局部放电的重要方法。
国内外众多学者开展了采用局部放电特高频信号时间差进行局部放电源定位的研究。利用多个特高频传感器接收信号的时间差,建立时间差定位方程组,通过求解该方程组获得局部放电源在电力设备内的位置,方便后期检修时及时发现局部放电源的所在位置,提高检修效率。但变电站现场噪声源众多,检测到的特高频信号包含较大幅值的噪声,使计算的特高频信号时间差精度降低,因此需要对特高频信号进行去噪处理。小波分解去噪方法具有优异的去噪性能,逐渐被众多学者及工程人员采用。然而,传统小波分解去噪方法(小波阈值去噪)以信噪比最高为目标,会造成局部放电信号波形的畸变,劣化特高频信号时差精度,使得去噪后的精度在许多情况下反而降低。
发明内容:
本发明的目的在于解决传统小波分解阈值去噪方法会畸变特高频信号波形的问题,提供了一种以减小局部放电特高频信号波形畸变为目标的小波分解二值去噪方法。
为了实现上述目的,本发明采用如下技术方案来实现的:
局部放电特高频信号波形的小波分解二值去噪方法,包括如下步骤:
1)在被检测的电力设备上安装局部放电特高频传感器,利用检测设备采集局部放电特高频信号;
2)选取与局部放电特高频信号波形相似的母小波,设置分解层数,对局部放电信号进行小波多尺度分解,得到各尺度小波系数与信号波形;
3)以减小局部放电特高频信号波形畸变为目标,对小波分解系数进行处理,具体方法为:计算各尺度信号累积能量函数的数学形态学梯度,并计算数学形态学梯度的陡峭度Ku;将各尺度信号的陡峭度与某一阈值Kt进行比较,若陡峭度大于该阈值,该尺度包含局部放电信号,否则仅包含噪声信号;将仅包含噪声的尺度的小波系数置为零,包含局部放电信号的尺度的小波系数保持不变;
4)根据处理后的小波系数进行小波重构,得到去噪后的局部放电特高频信号。
本发明进一步的改进在于,步骤2)所述的选取与局部放电特高频信号波形相似的母小波,为db系列与sym系列母小波。
本发明进一步的改进在于,步骤2)所述的设置分解层数为5~20层。
本发明进一步的改进在于,步骤2)所述的对局部放电信号进行小波多尺度分解,具体为采用Mallat多尺度算法进行小波分解。
本发明进一步的改进在于,步骤3)所述的计算各尺度信号累积能量函数的数学形态学梯度,并计算数学形态学梯度的陡峭度,包括以下步骤:
(a)若第i个尺度的信号为Ui(t),求取该信号的累积能量函数:
其中,E(tk)为tk时刻的累积能量函数值,N为信号的采样点数,tj为第j个采样点对应的时刻;为了去除信号幅度的影响,将累积能量函数除以总累积能量进行归一化;
(b)求取累积能量函数的数学形态学梯度,数学形态学的基本运算为腐蚀与膨胀,若累积能量函数为E,则E关于结构元素g的膨胀E⊕g和腐蚀EΘg分别定义为:
式中,DE、Dg分别为累积能量函数E及结构元素g的定义域,定义域为DE=1,2,…,N;结构元素采用扁平型,其定义域Dg=-SEL/2,…,0,…,SEL/2,SEL为结构元素的长度,取为50个采样点,n与m为采样点索引;累积能量函数的数学形态学梯度为膨胀与腐蚀运算之差:
mg(n)=E⊕g(n)-EΘg(n) (3)
(c)计算累积能量函数数学形态学梯度的陡峭度Ku,计算公式为:
其中,
其中,tk表示第k个采样点对应的时刻,mgk为对应时刻的数学形态学梯度数值;将tk作为随机变量,mgk作为对应的概率值,则pk为将数学形态学梯度数值作为概率计算得到概率密度值,μ为tk的1阶原点矩,σ2为tk的2阶中心矩。
本发明进一步的改进在于,步骤3)所述的阈值Kt,通过对大量特高频信号的统计分析得到,具体步骤为:
(a)在实验室内及现场电力设备内采集不少于20组局部放电特高频信号;
(b)根据选取的母小波与分解层数,对局部放电信号进行小波多尺度分解,得到各尺度小波系数与信号波形;对于各尺度信号,若信号波形呈现脉冲型,则该尺度包含局部放电信号;若信号幅值随时间变化不大,则该尺度仅包含噪声信号;
(c)根据陡峭度计算方法,计算所有信号各尺度的陡峭度;
(d)根据步骤(b)中判断各尺度是否包含局部放电信号,及步骤(c)计算的各尺度的陡峭度,分别计算仅包含噪声的尺度与包含局部放电信号的尺度的陡峭度范围,选取某一可将两者范围有效区分的阈值,作为步骤3)中的比较阈值Kt。
与现有技术相比,本发明具有以下的有益效果:
(1)本发明提出了局部放电特高频信号波形的小波分解二值去噪方法,基本思路为:将仅包含噪声的尺度的小波系数置为零,同时包含局放信号与噪声的尺度小波系数保持不变,去噪后波形畸变很小。
(2)将传统小波分解去噪应用于电力设备内部局部放电源的定位时,特高频信号被严重畸变,使后续确定信号起始时刻精度降低,反而劣化了定位精度;本发明提出的小波分解二值去噪方法可有效降低噪声水平,且对局部放电特高频信号波形畸变很小,较大程度提高特高频信号起始时刻的精度,提高局部放电源定位的准确度。
(3)将本发明提出的小波分解二值去噪方法应用于局部放电特高频信号特征提取时,使后续提取的特征参量受噪声影响程度降低,且不会改变去噪前后特征参量数值。
(4)本发明采用累积能量函数数学形态学梯度的陡峭度判断各尺度是否包含局部放电信号,该方法可用于局部放电信号与噪声信号的识别,进而去除无效的噪声信号。
附图说明:
图1为本发明小波分解二值去噪方法流程图。
图2为本发明检测得到的特高频信号波形图;其中,图2(a)为时间尺度在-500ns-500ns内的波形图,图2(b)为时间尺度放大后在-10ns-50ns之间的波形图。
图3为本发明小波分解的各尺度信号及其数学形态学梯度波形图;其中,图3(a)为小波多尺度分解后各尺度的信号波形图,图3(b)为各尺度信号的累积能量函数数学形态学梯度波形图。
图4为本发明仅包含噪声尺度与包含局部放电信号尺度的陡峭度统计图。
图5为本发明采用小波分解去噪方法前后的波形对比图;其中,图5(a)为原始采集的信号、加入信噪比为-3dB的白噪声后信号、传统小波分解去噪后信号及小波分解二值去噪后信号波形对比图,图5(b)为原始采集的信号、传统小波分解去噪后信号及小波分解二值去噪后信号波形在-5ns-10ns时间内对比图。
图6为本发明利用累积能量拐点法确定的原始信号、传统小波分解去噪后及小波二值分解去噪后时间差值对比图。
具体实施方式:
下面结合附图和实施例对本发明进一步详细说明。
如图1所示,本发明一种局部放电特高频信号波形的小波分解二值去噪方法,包括如下步骤:
(1)在被检测的电力设备上安装局部放电特高频传感器,利用检测设备采集局部放电特高频信号;
作为一种实施例,利用一套气体绝缘组合电器局部放电检测系统,采用带宽为500MHz~3GHz的圆盘形内置传感器检测局部放电信号,利用高速采集数字示波器以10GS/s的采样率检测该特高频信号,检测得到特高频信号在-500ns-500ns内波形如图2(a)所示,在时间尺度上放大后在-10ns-50ns之间的波形如图2(b)所示。
(2)选取与局部放电特高频信号波形相似的母小波,设置分解层数,对局部放电信号进行小波多尺度分解,得到各尺度小波系数与信号波形;
本步骤中所述的选取与局部放电特高频信号波形相似的母小波,主要为db系列与sym系列母小波;所述的设置分解层数,主要为5~20层。
本步骤中所述的对局部放电信号进行小波多尺度分解,具体为采用Mallat多尺度算法进行小波分解。
作为一种实施例,将采集到的特高频信号加入信噪比为-3dB的白噪声,选取db6母小波作为小波分解的基小波,采用Mallat多尺度算法对加入噪声后的特高频信号进行6层小波分解,图3(a)为小波分解得到的各尺度的波形。
(3)以减小局部放电特高频信号波形畸变为目标,对小波分解系数进行处理,具体方法为:计算各尺度信号累积能量函数的数学形态学梯度,并计算数学形态学梯度的陡峭度Ku;将各尺度信号的陡峭度与某一阈值Kt进行比较,若陡峭度大于该阈值,该尺度包含局部放电信号,否则仅包含噪声信号;将仅包含噪声的尺度的小波系数置为零,包含局部放电信号的尺度的小波系数保持不变;
本步骤中所述的计算各尺度信号累积能量函数的数学形态学梯度,并计算数学形态学梯度的陡峭度,包括以下步骤:
(a)若第i个尺度的信号为Ui(t),求取该信号的累积能量函数:
其中,E(tk)为tk时刻的累积能量函数值,N为信号的采样点数,tj为第j个采样点对应的时刻。为了去除信号幅度的影响,将累积能量函数除以总累积能量进行归一化。
(b)求取累积能量函数的数学形态学梯度,数学形态学的基本运算为腐蚀与膨胀,若累积能量函数为E,则E关于结构元素g的膨胀E⊕g和腐蚀EΘg分别定义为:
式中,DE、Dg分别为累积能量函数E及结构元素g的定义域,定义域为DE=1,2,…,N;结构元素采用扁平型,其定义域Dg=-SEL/2,…,0,…,SEL/2,SEL为结构元素的长度,取为50个采样点,n与m为采样点索引。累积能量函数的数学形态学梯度按下式计算:
mg(n)=E⊕g(n)-EΘg(n) (3)
(c)计算累积能量函数数学形态学梯度的陡峭度,计算公式为:
其中,
其中,tk表示第k个采样点对应的时刻,mgk为对应时刻的数学形态学梯度数值;将tk作为随机变量,mgk作为对应的概率值,则pk为将数学形态学梯度数值作为概率计算得到概率密度值,μ为tk的1阶原点矩,σ2为tk的2阶中心矩。
作为一种实施例,计算图3(a)小波分解各尺度信号的数学形态学梯度,数学形态学梯度波形如图3(b)所示,计算得到d1、d2、d3、d4、d5、d6及a6尺度的陡峭度分别为-1.2、-1.2、-0.54、5.49、3.59、-1.15及-1.26。根据各尺度波形,尺度d1、d2、d6及a6波形在整个时间内幅值基本相同,故这些尺度仅包含噪声;尺度d3、d4及d5均包含明显的脉冲信号,说明这些尺度包含局部放电信号。尺度d3、d4及d5包含局部放电信号,可以看出其陡峭度明显大于仅包含噪声的d1、d2、d6及a6尺度。
本步骤中所述的阈值Kt,需要通过对大量特高频信号的统计分析得到,具体步骤为:
(a)在实验室内及现场电力设备内采集不少于20组局部放电特高频信号;
(b)根据选取的母小波与分解层数,对局部放电信号进行小波多尺度分解,得到各尺度小波系数与信号波形;对于各尺度信号,若信号波形呈现脉冲型,则该尺度包含局部放电信号;若信号幅值随时间变化不大,则该尺度仅包含噪声信号;
(c)根据陡峭度计算方法,计算所有信号各尺度的陡峭度;
(d)根据步骤(b)中判断各尺度是否包含局部放电信号,及步骤(c)计算的各尺度的陡峭度,分别计算仅包含噪声的尺度与包含局部放电信号的尺度的陡峭度范围,选取某一可将两者范围有效区分的阈值,作为步骤(3)中的比较阈值Kt。
作为一种实施例,本发明采用步骤(1)中实验布置测量了20组特高频信号,并根据步骤(2)对所有信号进行小波分解,针对小波分解后各尺度信号是否含有脉冲判断该尺度是否含有局部放电信号,根据步骤(3)计算了所有信号小波分解后数学形态学梯度的陡峭度,如图4所示。可以看出包含局部放电信号尺度的数学形态学梯度陡峭度范围为-0.9~6,仅包含噪声的尺度的数学形态学梯度陡峭度范围为-1.3~-0.7,两者数值范围存在相交区域-0.9~-0.7,但该区域内样本点数目已经很少,说明将阈值Kt取为-0.8可基本将包含局部放电尺度与噪声尺度区分开。由于某些尺度包含较弱的局部放电信号,其Ku与仅含噪声尺度比较接近,因此很难得到一个可以完全区分包含局部放电尺度与仅包含噪声尺度的阈值。如果需要保留更多的局部放电信息,可适当调低Kt值;若想更多地去除噪声,可适当提高Kt值。
本步骤中所述的将陡峭度与某一阈值进行比较,若陡峭度大于该阈值认为该尺度包含局部放电信号,否则仅包含噪声信号;将仅包含噪声的尺度的小波系数置为零,包含局部放电信号的尺度小波系数保持不变。
作为一种实施例,以前述步骤确定的-0.8作为阈值,判断图1信号各尺度是否包含局部放电信号,根据前述计算,d1、d2、d3、d4、d5、d6及a6尺度的陡峭度分别为-1.2、-1.2、-0.54、5.49、3.59、-1.15及-1.26。尺度d3、d4及d5陡峭度大于阈值-0.8,说明其包含局部放电信号,这些尺度小波系数保持不变;d1、d2、d6及a6尺度的陡峭度小于-0.8,这些尺度仅包含噪声,其小波系数置为零。
(4)根据处理后的小波系数进行小波重构,得到去噪后的局部放电特高频信号。
作为一种实施例,对图2(a)所示信号加入信噪比为-3dB的白噪声,并分别采用传统小波分解阈值去噪与小波分解二值去噪方法对其进行去噪。图5(a)为检测到的特高频信号、加入信噪比为-3dB白噪声后的特高频信号、传统小波分解阈值去噪后的特高频信号及小波分解二值去噪后的特高频信号波形。可以看出,传统小波分解阈值去噪去掉了全部噪声;小波分解二值去噪将原始峰-峰值0.012V的噪声降低为0.005V,去除了大部分噪声。图5(b)为原始采集的信号、传统小波分解去噪后信号及小波分解二值去噪后信号波形在-5ns-10ns时间内对比图,分析了去噪前后特高频信号波形的变化及去噪对信号起始时刻的影响,可以看出小波分解二值去噪对信号波形畸变很小,对信号起始时刻影响较小;传统小波分解阈值去噪后信号波形畸变严重,造成信号起始时刻误差增大。
为了进一步说明本发明提出的小波分解二值去噪算法的有益效果,本发明同时利用两个不同位置的特高频传感器测量局部放电信号,根据局部放电源与两个传感器位置的距离差及传播速度,可计算出信号真实时间差。本发明通过检测得到的特高频信号确定时间差,具体步骤为:利用小波分解算法进行对特高频信号进行去噪,利用累积能量拐点法确定信号之间的时间差。图6为利用某次测量原始信号、传统小波分解去噪后信号及小波二值分解去噪后信号确定的时间差值,可以看出传统小波分解去噪后确定的时间差与真实时间差-2.55ns相差极大,而小波二值分解去噪后确定的时间差与真实时间差-2.55ns最接近,且受噪声水平的影响较小。
本领域的技术人员容易理解,以上所述仅作为本发明的实施案例,并不用以限制本发明,凡在本发明的精神和原则之内所做的任何修改、替换或变更,均应包含在本发明的保护范围之内。
Claims (4)
1.局部放电特高频信号波形的小波分解二值去噪方法,其特征在于,包括如下步骤:
1)在被检测的电力设备上安装局部放电特高频传感器,利用检测设备采集局部放电特高频信号;
2)选取与局部放电特高频信号波形相似的母小波,设置分解层数,对局部放电特高频信号进行小波多尺度分解,得到各尺度小波系数与信号波形;
3)以减小局部放电特高频信号波形畸变为目标,对小波分解系数进行处理,具体方法为:计算各尺度信号累积能量函数的数学形态学梯度,并计算数学形态学梯度的陡峭度Ku;将各尺度信号的陡峭度与某一阈值Kt进行比较,若陡峭度大于该阈值,该尺度包含局部放电特高频信号,否则仅包含噪声信号;将仅包含噪声的尺度的小波系数置为零,包含局部放电特高频信号的尺度的小波系数保持不变;
其中,所述的计算各尺度信号累积能量函数的数学形态学梯度,并计算数学形态学梯度的陡峭度,包括以下步骤:
(a)若第i个尺度的信号为Ui(t),求取该信号的归一化累积能量函数:
<mrow>
<mi>E</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>/</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>*</mo>
<mn>100</mn>
<mo>,</mo>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>N</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,E(tk)为tk时刻的归一化累积能量函数值,N为信号的采样点数,tj为第j个采样点对应的时刻;为了去除信号幅度的影响,将累积能量函数除以总累积能量进行归一化;
(b)求取累积能量函数的数学形态学梯度,数学形态学的基本运算为腐蚀与膨胀,若累积能量函数为E,则E关于结构元素g的膨胀和腐蚀EΘg分别定义为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>E</mi>
<mi>&Theta;</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>max</mi>
<mo>{</mo>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>-</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>&Element;</mo>
<msub>
<mi>D</mi>
<mi>E</mi>
</msub>
<mo>,</mo>
<mi>m</mi>
<mo>&Element;</mo>
<msub>
<mi>D</mi>
<mi>g</mi>
</msub>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>E</mi>
<mi>&Theta;</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>min</mi>
<mo>{</mo>
<mi>E</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>)</mo>
</mrow>
<mo>&Element;</mo>
<msub>
<mi>D</mi>
<mi>E</mi>
</msub>
<mo>,</mo>
<mi>m</mi>
<mo>&Element;</mo>
<msub>
<mi>D</mi>
<mi>g</mi>
</msub>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,DE、Dg分别为累积能量函数E及结构元素g的定义域,定义域为DE=1,2,…,N;结构元素采用扁平型,其定义域Dg=-SEL/2,…,0,…,SEL/2,SEL为结构元素的长度,取为50个采样点,n与m为采样点索引;累积能量函数的数学形态学梯度为膨胀与腐蚀运算之差:
<mrow>
<mi>m</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>E</mi>
<mo>&CirclePlus;</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>E</mi>
<mi>&Theta;</mi>
<mi>g</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
(c)计算累积能量函数数学形态学梯度的陡峭度Ku,计算公式为:
<mrow>
<mi>K</mi>
<mi>u</mi>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<mi>&mu;</mi>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
<msub>
<mi>p</mi>
<mi>k</mi>
</msub>
</mrow>
<msup>
<mi>&sigma;</mi>
<mn>4</mn>
</msup>
</mfrac>
<mo>-</mo>
<mn>3</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,
其中,tk表示第k个采样点对应的时刻,mgk为对应时刻的数学形态学梯度数值;将tk作为随机变量,mgk作为对应的概率值,则pk为将数学形态学梯度数值作为概率计算得到概率密度值,μ为tk的1阶原点矩,σ2为tk的2阶中心矩;
所述的阈值Kt,通过对大量特高频信号的统计分析得到,具体步骤为:
(a)在实验室内及现场电力设备内采集不少于20组局部放电特高频信号;
(b)根据选取的母小波与分解层数,对局部放电特高频信号进行小波多尺度分解,得到各尺度小波系数与信号波形;对于各尺度信号,若信号波形呈现脉冲型,则该尺度包含局部放电特高频信号;若信号幅值随时间变化不大,则该尺度仅包含噪声信号;
(c)根据陡峭度计算方法,计算所有信号各尺度的陡峭度;
(d)根据步骤(b)中判断各尺度是否包含局部放电特高频信号,及步骤(c)计算的各尺度的陡峭度,分别计算仅包含噪声的尺度与包含局部放电特高频信号的尺度的陡峭度范围,选取某一可将两者范围有效区分的阈值,作为步骤3)中的比较阈值Kt;
4)根据处理后的小波系数进行小波重构,得到去噪后的局部放电特高频信号。
2.根据权利要求1所述的局部放电特高频信号波形的小波分解二值去噪方法,其特征在于,步骤2)所述的选取与局部放电特高频信号波形相似的母小波,为db系列与sym系列母小波。
3.根据权利要求1所述的局部放电特高频信号波形的小波分解二值去噪方法,其特征在于,步骤2)所述的设置分解层数为5~20层。
4.根据权利要求1所述的局部放电特高频信号波形的小波分解二值去噪方法,其特征在于,步骤2)所述的对局部放电特高频信号进行小波多尺度分解,具体为采用Mallat多尺度算法进行小波分解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510697143.7A CN105223482B (zh) | 2015-10-22 | 2015-10-22 | 局部放电特高频信号波形的小波分解二值去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510697143.7A CN105223482B (zh) | 2015-10-22 | 2015-10-22 | 局部放电特高频信号波形的小波分解二值去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105223482A CN105223482A (zh) | 2016-01-06 |
CN105223482B true CN105223482B (zh) | 2018-01-19 |
Family
ID=54992546
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510697143.7A Expired - Fee Related CN105223482B (zh) | 2015-10-22 | 2015-10-22 | 局部放电特高频信号波形的小波分解二值去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105223482B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106950475A (zh) * | 2017-03-23 | 2017-07-14 | 广东电网有限责任公司珠海供电局 | 一种基于小波变换的局部放电信号提取方法及装置 |
CN108804832B (zh) * | 2018-06-14 | 2020-06-16 | 东南大学 | 一种基于vmd的间隔阈值局部放电去噪方法 |
CN109145762B (zh) * | 2018-07-27 | 2022-02-01 | 西南石油大学 | 一种基于数学形态学和分形理论的电缆附件局部放电识别方法 |
CN109541405B (zh) * | 2018-10-23 | 2020-10-16 | 深圳供电局有限公司 | 一种电力电缆局部放电定位方法 |
CN109932624B (zh) * | 2019-04-01 | 2021-07-02 | 珠海华网科技有限责任公司 | 一种基于高斯尺度空间的电缆局放周期窄带干扰去噪方法 |
CN110244202B (zh) * | 2019-06-21 | 2021-07-20 | 三峡大学 | 基于同步压缩小波变换域变压器局部放电去噪方法 |
CN111007369B (zh) * | 2019-12-27 | 2021-10-19 | 广东电网有限责任公司电力科学研究院 | 一种超高频电磁波信号到达时间差计算方法和装置 |
CN111535999B (zh) * | 2020-05-22 | 2021-08-24 | 三一重能有限公司 | 风机脱落物监测方法、装置、系统和存储介质 |
CN111707908B (zh) * | 2020-07-29 | 2022-05-17 | 中国科学技术大学先进技术研究院 | 多负载回路串联故障电弧检测方法、装置以及存储介质 |
CN112464876A (zh) * | 2020-12-12 | 2021-03-09 | 南方电网数字电网研究院有限公司 | 电力设备的故障诊断方法、装置、计算机设备和存储介质 |
CN113189457A (zh) * | 2021-04-26 | 2021-07-30 | 天津大学 | 一种基于局部放电原始时域波形的多尺度特征提取技术 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102590716A (zh) * | 2012-02-11 | 2012-07-18 | 广东电网公司电力科学研究院 | 一种基于小波分析的电力电缆局部放电信号去噪方法 |
EP2518520A2 (en) * | 2011-04-29 | 2012-10-31 | General Electric Company | Device and method for detecting and locating defects in underground cables |
CN103543390A (zh) * | 2013-09-25 | 2014-01-29 | 国家电网公司 | 电力变压器局部放电超高频信号的去噪方法 |
CN103558517A (zh) * | 2013-10-29 | 2014-02-05 | 云南电网公司楚雄供电局 | 基于能量变化定位高压设备放电源的方法 |
CN103576060A (zh) * | 2013-10-11 | 2014-02-12 | 华南理工大学 | 基于小波自适应阈值的局部放电信号去噪方法 |
CN104198901A (zh) * | 2014-08-13 | 2014-12-10 | 广东电网公司电力科学研究院 | 变电站局部放电信号的定位方法和系统 |
CN104614578A (zh) * | 2015-02-12 | 2015-05-13 | 国家电网公司 | 一种波头时间的确定方法及系统 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0504600D0 (en) * | 2005-03-04 | 2005-04-13 | Univ Strathclyde | Detecting partial discharge in high voltage cables |
KR20070094680A (ko) * | 2006-03-18 | 2007-09-21 | 한국전기연구원 | 차동 기법을 이용한 외부 잡음 제거 방식을 채용한발전기의 운전중의 부분 방전 측정방법 |
JP2007292700A (ja) * | 2006-04-27 | 2007-11-08 | Toshiba Corp | 静止誘導機器の部分放電位置特定方法 |
KR101803559B1 (ko) * | 2010-12-01 | 2017-12-01 | 한양대학교 에리카산학협력단 | 가스절연기기의 부분방전 검출방법 및 검출시스템 |
KR20140067653A (ko) * | 2012-11-27 | 2014-06-05 | 한국전력공사 | 부분방전 감지장치의 진단을 위한 부분방전 분석장치 및 방법 |
-
2015
- 2015-10-22 CN CN201510697143.7A patent/CN105223482B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2518520A2 (en) * | 2011-04-29 | 2012-10-31 | General Electric Company | Device and method for detecting and locating defects in underground cables |
CN102590716A (zh) * | 2012-02-11 | 2012-07-18 | 广东电网公司电力科学研究院 | 一种基于小波分析的电力电缆局部放电信号去噪方法 |
CN103543390A (zh) * | 2013-09-25 | 2014-01-29 | 国家电网公司 | 电力变压器局部放电超高频信号的去噪方法 |
CN103576060A (zh) * | 2013-10-11 | 2014-02-12 | 华南理工大学 | 基于小波自适应阈值的局部放电信号去噪方法 |
CN103558517A (zh) * | 2013-10-29 | 2014-02-05 | 云南电网公司楚雄供电局 | 基于能量变化定位高压设备放电源的方法 |
CN104198901A (zh) * | 2014-08-13 | 2014-12-10 | 广东电网公司电力科学研究院 | 变电站局部放电信号的定位方法和系统 |
CN104614578A (zh) * | 2015-02-12 | 2015-05-13 | 国家电网公司 | 一种波头时间的确定方法及系统 |
Non-Patent Citations (9)
Title |
---|
A combined CWT-DWT method using model-based design simulator for partial discharges online detection;R. Candela 等;《Annual Report Conference on Electrical Insulation and Dielectric Phenomena》;20090831;429-432 * |
Time Delay Estimation for UHF Signals in PD location of Transformers;L. Yang 等;《Annual Report Conference on Electrical Insulation and Dielectric Phenomena》;20041213;414-417 * |
UHF传感器GIS局部放电检测特性的实验研究;邵先军 等;《高压电器》;20150116;第56卷(第1期);46-55 * |
变压器局部放电定位技术及新兴UHF方法的关键问题;唐志国 等;《南方电网技术》;20080228;第2卷(第1期);36-40 * |
变压器绕组对局部放电特高频定位方法的影响;郑书生 等;《高电压技术》;20130228;第39卷(第2期);324-329 * |
基于特高频法的GIS局部放电故障诊断研究;王彩雄;《中国博士学位论文全文数据库工程科技Ⅱ辑》;20150715(第7期);C042-133 * |
局部放电在线监测中小波阑值去噪法的最优阑值自适应选择;李剑 等;《电网技术》;20060430;第30卷(第8期);25-30 * |
应用特高频检测方法处理GIS局部放电故障;廖天明 等;《华东电力》;20110430;第39卷(第4期);0652-0654 * |
特高频方法用于套管局部放电检测的可行性研究;肖天为 等;《电网技术》;20150630;第39卷(第6期);1719-1724 * |
Also Published As
Publication number | Publication date |
---|---|
CN105223482A (zh) | 2016-01-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105223482B (zh) | 局部放电特高频信号波形的小波分解二值去噪方法 | |
CN105223481B (zh) | 基于差值能量函数的局部放电特高频信号起始时刻确定方法 | |
CN109557429B (zh) | 基于改进小波阈值去噪的gis局部放电故障检测方法 | |
CN105785324B (zh) | 基于mgcstft的线性调频信号参数估计方法 | |
CN109633633B (zh) | 一种基于分段分类增强处理的生命信号增强方法 | |
CN106856594A (zh) | 基于rssi的室内定位方法与系统 | |
CN113325277A (zh) | 一种局部放电处理方法 | |
CN109212503B (zh) | 一种基于ar谱分形相关特征的海上目标检测方法及装置 | |
CN114785379A (zh) | 一种水声janus信号参数估计方法及系统 | |
CN112986922B (zh) | 基于顺序统计滤波和二元检测的被动雷达信号检测方法 | |
CN111505598A (zh) | 一种基于frft域三特征联合检测装置及方法 | |
CN113589253A (zh) | 一种基于伪时域的小波变换算法对微弱回波信号进行检测的方法 | |
CN110287853B (zh) | 一种基于小波分解的暂态信号去噪方法 | |
CN110632563A (zh) | 一种基于短时傅里叶变换的脉内频率编码信号参数测量方法 | |
CN105636061B (zh) | 一种在一般性衰落信道上统一精确的能量检测方法 | |
CN106330361B (zh) | 基于图连通性分析的弱信号检测方法 | |
CN106771928A (zh) | 一种局部放电脉冲初至时刻在线拾取方法 | |
CN110595529A (zh) | 一种强背景噪声下微弱周期信号的快速检测方法 | |
CN107315169B (zh) | 基于二阶统计量相似度的杂波协方差矩阵估计方法 | |
JP6258574B2 (ja) | パッシブソーナー装置、方位集中処理方法、及び、パッシブソーナー信号処理プログラム | |
CN114966550A (zh) | 一种基于经验小波变换的自适应海杂波抑制方法 | |
CN103441975B (zh) | 一种基于功率谱的二相编码信号参数估值方法 | |
CN105738878A (zh) | 一种基于频偏分类的精密测量雷达测频方法 | |
CN109793511A (zh) | 基于深度学习技术的心电信号噪声检测算法 | |
CN110941908A (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 | ||
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: 20180119 |