局部放电在线监测中基于小波变换的阈值消噪算法研究

发布于:2007/6/19 9:20:31 | 2053 次阅读

 

    摘  要:小波变换的正交多分辨率分析可将信号在不同尺度下进行分解。通过对各尺度上的分解系数进行合适的阈值处理,可有效地消除局部放电在线监测中的窄带干扰,提取脉冲信号。通过对仿真信号和现场信号的处理,结果表明该方法行之有效。
   
关键词:局部放电;小波变换;阈值;消噪

引言
    局部放电在线监测是一种广泛用于检测电力设备绝缘性能的重要手段,但由于现场监测过程中存在大量的干扰信号,如何消除各种干扰的影响,特别是载波通信等缓变窄带干扰,是提高局部放电在线监测灵敏度的一个关键技术问题。在信号处理领域中,传统的傅氏(Fourier)变换是滤除噪声、提高信噪比(SNR)的常用方法,但这种方法也会带来一定的负面效应,当用低通滤波器对信号滤波时,在滤除噪声的同时也对信号的边沿细节有一定的平滑作用,从而使信号发生失真。
    小波变换是继Fourier变换以后发展起来的一种新的信号处理方法,不同于Fourier变换,小波变换的基函数是紧支撑的,这就使得小波变换可以描述信号的时频局部化信息。它可以在多个分辨率下对信号进行分解,给出不同尺度下的信号细节信息。信号和噪声一般都分布在不同的频带上,对应于小波变换各个不同尺度。本文通过对小波变换的不同尺度进行研究,找出了一种行之有效的滤除噪声保存信号的方法。
2  基于小波变换的多分辨率分析
    小波是函数空间L2(R)中满足下述条件的一个函数Ψ(x)

对于任意的实数对(a,b),其中a为非零实数,称如下形式的函数为由小波母函数Ψ(x)生成的依赖于参数(a,b)的连续小波函数

式中 a为尺度参数;b为时间参数。
    1988年S.Mallat在信号多分辨率分析(Multi-Resolution Analysis)空间分解概念的基础上,将二进尺度小波变换与多分辨率分析结合起来,提出了小波变换的快速算法¾¾Mallat算法,从空间概念上形象地说明了小波的多分辨率特性[1]。多分辨率分析是将信号在L2(R)的两个正交子空间上逐级分解,将每级输入分解为低频概貌和高频细节两部分,从而使输出采样率减半。利用多分辨率分析及尺度方程和小波方程的系数,可得到信号小波变换的Mallat算法,即对于任意信号f(t)L2(R),有如下关系式

式中aj(k),dj(k)分别为尺度j上信号的尺度变换系数和小波变换系数,即信号的概貌和细节部分,h0(k),h1(k)是满足二尺度差分方程[2] 的两个滤波器系数。
    正交多分辨率分析是将整个信号的频带不断以二进方式划分,而实际测量信号的分辨率是有限的,这样就得到了由低到高地分布在不同频带的信号成分,即信号的高低频成分,从而可将信号的有用频率成分和属于干扰的频率成分区分开来,再经过阈值处理和信号重构,可得到有用信号。
3  基于小波变换的消噪算法
3.1  消噪模型
    自Witikin首先提出利用信号频带分解的尺度空间系数清除噪声的思想以来,已发展出了许多基于小波分析的噪声滤波方法[3],其中Donoho[4]提出的含噪声的一维信号模型表示为

式中 di为含噪声信号,f(ti)为有用信号,zi为噪声信号,σ为噪声水平。
   
消噪目的是使得信号的均方差达到,即

    式(7)是一个zi-N(0,1)的白噪声统计模型,但此模型同样适用于非白噪声,可将其用于非白噪声的消噪处理。
3.2  消噪方法
    Donoho[4]提出的基于小波变换的消噪方法可分以下三个步骤: 选择小波基函数和小波分解的尺度,对含噪声的信号进行小波分解; 选择阈值,对各个尺度的小波系数选择合适的阈值来进行处理; 将经过处理后的各尺度系数进行重构,得到消噪后的信号。
    由此可见,小波消噪就是用阈值对小波分解系数的量化处理,其中最重要的环节就是如何选取阈值和如何进行阈值的量化。在阈值选取过程中,既要保证能有效去除噪声,又要保证不能将有用信号当作噪声消除掉,而且还要尽可能保持有用信号的大小不变。基于Donoho[4]的工作出现了源于统计学原理的阈值选取算法,如Stein的无偏似然估计原理,即对一个给定的阈值t,得到它的似然估计,再将非似然t最小化,就得到了所选的阈值,这是一种软件阈值估计器。
    (1)小波基函数的选择
    小波变换不同于傅氏变换,它对信号进行变换时可采用不同的基函数,而且对于特定信号而言采用的基函数不同,其分析结果也会相差很大。小波基主要有下列几个特征:紧支性、衰减性、正则性,对称性、消失距。这些特征关系到如何选择合适的小波基,以便高效地分析信号。为了分析局部放电这种突变信号,在选择小波基时主要考虑满足给定区间的紧支性和足够的消失距,这样能有效地消除噪声,提取突变信号。Daubechies系列小波基是典型的具有紧支光滑的正交小波基,其它几大类(双正交Biorthogonal小波基系列,Coiflets小波基系列,Symlets小波基系列)都由Daubechies系列小波基推广、引伸得到[5]。因此本文选择Daubechies系列小波基作为分析局部放电信号的基函数。
    (2)小波分解尺度的选择
    在正交多分辨率分析中,小波分解的尺度越大,信号的分辨率越高。但在用小波变换对信号进行消噪处理过程中,如果小波分解的尺度选择过大,并不能明显地提高消噪效果,反而会增大计算量,但小波分解尺度如果选择过小,就不能有效地消除噪声。如果采样频率如为fs,由香农采样定理知采样信号的频率fmax=fs/2,窄带干扰的频率范围为f1-f2,且令


    由于信号的正交小波变换是将信号的频率按二进制划分成一系列子频带,消噪时只要尽量将干扰信号频率分解到各个子频带中,就可以在各尺度上进行有效的阈值处理。由式(8)可见,J值是分解窄带干扰信号的尺度,因此在对信号进行正交小波变换时,可将J值作为参考值来选取合适的分解尺度。
    (3)阈值的处理方法
    小波消噪通常是一种基于软阈值或硬阈值处理的消噪方法。在根据信号的小波系数选取合适的阈值后,可选取合适的阈值处理方法来应用于小波系数。
    由图1(b)可知硬阈值是将信号的与阈值进行比较,小于或等于阈值的点变为0,大于阈值的点保持不变,见图1(b)。硬阈值函数可表示为[6]


    软阈值是将信号的与阈值进行比较,小于或等于阈值的点变为0,大于阈值的点变为该点值与阈值的差值,见图1(c)。软阈值函数可表示为[6]

    由此可见,软阈值是对硬阈值的扩展。一般情况下,采用软阈值处理方法可获得比较平滑的消噪信号,然而对于从大量的窄带干扰等平衡信号中提取出具有突变特性的局部放电信号而言,采用硬阈值处理方法更为合适,见图2。比较图1(b)与图1(c)可知,软阈值处理方法会使处理后的脉冲峰值减小(有时甚至到0),从而无法判断出是否发生局部放电。而硬阈值处理方法可比较好地保持脉冲峰值不变,从而能正确判断出是否发生了局部放电。

4 实验分析
4.1 小波消噪对模拟信号的处理结果
    仿真的局部放电信号为衰减振荡信号,幅值为2,仿真的背景干扰为现场常见的一些窄带干扰信号。由于局部放电信号属于突变信号,所以这里选择Daubechies小波基作为小波变换的基函数,小波分解尺度为6,采用基于Stein的无偏以然估计原理的阈值估计器、硬阈值处理方法选择阈值来对模拟信号进行消噪处理,处理结果见图3。

    由图可见,经过消噪处理后局部放电信号可完全被提取出来。由此可见,在选择合适的小波基函数和阈值处理方式的情况下,可有效地消除窄带干扰等平衡信号。
4.2 现场数据的处理结果
    现场信号为某变电所220kV三相变压器的地线电流,检测系统通频带为500kHz,采样频率为2MHz。图4(a)(b)分别为原始信号波形及其频谱分析结果,由图可见周期信号形成尖的谱线。选取“db2”作为小波变换的基函数,小波分解的尺度为6,同样采用基于Stein的无偏似然估计原理选取阈值,对信号进行硬阈值处理,最终得到消噪后的信号,见图4(c)(d)。由消噪后所得的脉冲信号、剩余的噪声信号及它们的频谱可见,基于小波变换的消噪方法可以比较有效地消除窄带信号和白噪等平稳信号对局部放电信号的干扰,从而为更好地提取和识别局部放电信号打下基础。

5  结论
    小波消噪方法的有效性主要取决于小波基函数、小波分解尺度和阈值处理函数等几个因素。从对仿真信号和现场信号的小波消噪结果来看,基于小波变换和通过对信号的各层分解系数进行阈值处理的消噪技术,可以对含有噪声的局部放电信号进行有效的消噪处理。

参考文献

[1]   Mallat SSingularity detection and processing with wavelets[J]IEEE Trans on IT199238(2)617-643
[2]   冉启文.小波变换与分数傅里叶变换理论及应用[M].哈尔滨:哈尔滨工业大学出版社,2001
[3]   仇智华,艾德才.基于子波分析的信号噪声分离方法研究[J].天津大学学报,200033(2)262-264
[4]   Donoho D LDe-noise by soft-thresholding[J]IEEE Trans. on Information Theory199541(2)613-627
[5]   谭善文,秦树人,汤宝平.小波基时频特性及其在分析突变信号中的应用[J].重庆大学学报,200124(2)12-16
[6]   Shim ISoraghan J JSiew W HDetection of PD utilizing digital signal porcessing methods part 3open-loop noise reduction[J]IEEE Electrical Insulation Magazine200117(1)6-12

参与讨论
后参与讨论

//评论区

推荐阅读

智能网联汽车国际标准法规协调专家组(HEAG)召开工作会议

近年来智能网联汽车快速发展,新技术不断涌现,与相关产业融合度持续提升,正在推动全球汽车产业发生深刻变革。为应对此种形势,欧、美、日等汽车工业发达国家和地区都加大了智能网联汽车的国际标准法规协调的参与力度,在联合国世界车辆法规论坛(UN/WP.29)和国际标准化组织(ISO)层面,智能网联汽车相关国际标准法规协调活动正快速推进。 为更有效地支撑上述组织的国际标准法规协调活动,2017年全国汽车标准

0215jiejie | 发布于:2022-12-01 0评论 0赞

苹果推出搭载M2芯片的新款iPad Pro 799美元起售

据苹果官网,苹果推出搭载M2芯片的新款iPadPro。 11英寸wifi版起售价为799美元,wifi+蜂窝网络版起售价为999美元;12.9英寸wifi版起售价为1099美元,wifi+蜂窝网络版起售价为1299美元。

0215jiejie | 发布于:2022-10-19 0评论 0赞

新能源汽车领衔 “中国智造”加速登陆欧洲市场

全球五大车展之一巴黎车展时隔四年再度启幕。在这场被视为“全球汽车行业风向标”的盛会上,国内外汽车品牌云集,长城汽车、比亚迪等再次领衔中国汽车出海。 长城汽车欧洲区域总裁孟祥军表示:“欧洲是长城汽车最重要的海外市场之一,巴黎车展是长城汽车向欧洲市场展示GWM品牌和产品的最佳机会。长城汽车正在研究汽车行业碳排放的整个生命周期,到2025年,将推出50多款新能源产品,全力支持可再生能源使用,为全球用户

0215jiejie | 发布于:2022-10-19 0评论 0赞

严监管时代来临,电子烟“通配”大战走向何方?

针对通配烟弹厂商的一系列诉讼的结果,将对生产通配烟弹的品牌未来在电子烟行业的发展产生深远影响。 10月1日,《电子烟强制性国家标准》正式实施,中国电子烟监管全面生效。而在电子烟行业进入规范化、法治化阶段前夕,一场围绕着通配烟弹的争论在行业里发酵。 “通配”是电子烟从业者约定俗成的概念。换弹式电子烟由烟杆和烟弹组成,“通配”烟弹指的是非品牌商生产、可与品牌烟杆匹配使用的烟弹。多位业内人士表示,被

0215jiejie | 发布于:2022-10-19 0评论 0赞

Bourns 全新大功率分流电阻器

采用金属感应引脚,专用于大电流应用中进行精确测量 全新分流电阻器专为电池管理系统、大电流工业控制和电动汽车充电站 提供高可靠性、高成本效益的解决方案 美国柏恩Bourns全球知名电子组件领导制造供货商,宣布新增12款CSM2F系列功率分流电阻器,扩展其产品组合。全新系列采用铆接通孔金属传感引脚,可满足大电流应用中对电压测试点精确定位日益增长的需求。最新型Bourns?CSM2F系列分流电阻器

0215jiejie | 发布于:2022-10-18 0评论 0赞

请尊重元宇宙“这个筐”

元宇宙是个筐,啥都往里装,但区别在于有的像聚宝盆,有的像垃圾桶。国庆假期刚结束,中青宝“90后”董事长李逸伦便亲自上阵,玩起了元宇宙婚礼。靠着老板首秀和代言,中青宝顺势推出“MetaLove元囍”App,正式进军元宇宙婚礼赛道。 就产品而言,如同其他元宇宙产品,李逸伦的元宇宙婚礼“新奇与吐槽齐飞”:有人说是有趣的尝试,有人则认为像QQ炫舞结婚系统。要知道,QQ炫舞是一款推出了十余年的老游戏。

0215jiejie | 发布于:2022-10-13 0评论 0赞

边缘计算:突围商业模式痛点

截至8月末,中国5G基站总数达210.2万个,中国5G发展已经进入下半场。随着5G加速融入千行百业,互动直播、vCDN、安防监控等场景率先大规模落地,车联网、云游戏、工业互联网、智慧园区、智慧物流等场景也快速走向成熟,这些更大流量、更低时延、更高性能的场景涌现,对边缘计算的刚性需求势必爆发。 GrandViewResearch预测,即使在新型冠状病毒肺炎疫情肆虐全球的背景下,边缘计算和5G网络市

0215jiejie | 发布于:2022-10-13 0评论 0赞

商务部回应美商务部升级半导体等领域对华出口管制并调整出口管制“未经验证清单”

商务部新闻发言人10日就美商务部升级半导体等领域对华出口管制并调整出口管制“未经验证清单”应询答记者问。 有记者问:近日,美国商务部在半导体制造和先进计算等领域对华升级出口管制措施。同时,在将9家中国实体移出“未经验证清单”过程中,又将31家中国实体列入,请问中方对此有何回应? 对此,商务部新闻发言人回应称,中方注意到相关情况。首先,通过中美双方前一阶段共同努力,9家中国实体zui终

0215jiejie | 发布于:2022-10-13 0评论 0赞

TCL华星官宣与奔驰合作:推出全球首款横贯A柱的车载显示屏

今年1月,奔驰带来了VISIONEQSS概念车,其中控台采用了一块完全无缝的47.5英寸曲面显示屏,横贯整个A柱,令人印象深刻。今天,TCL华星正式官宣与奔驰达成合作,并认领了VISIONEQSS上这块全球首款横贯整个A柱曲面的车载显示屏。 根据TCL介绍,这款显示屏采用了完全无缝的超薄一体化设计,将仪表盘、中控与副驾娱乐显示融为一体,并能够与3D实时导航系统相辅相成。 同时,这块显示屏还采用

0215jiejie | 发布于:2022-10-12 0评论 0赞

半导体板块暴跌 谁最受伤

国庆假期后首日开盘,上证综指时隔5个月再次失守3000点,与此同时,半导体板块也再度走低,其中,北方华创、雅克科技等个股跌停。10月11日早盘期间,半导体板块持续下挫,北方华创、雅克科技再度跌停。截至下午收盘,北方华创、雅克科技维持跌停状态,华海清科、拓荆科技-U、盛美上海、清溢光电、海光信息的跌幅则超10%。同日,半导体板块中的119只个股中超五成呈现下跌趋势。 在半导体板块遭遇下挫的同时,北

0215jiejie | 发布于:2022-10-12 0评论 0赞