一种基于小波阈值去噪的IFRA频响曲线去噪方法与流程

本发明涉及电力设备诊断与保护技术领域，尤其是指一种基于小波阈值去噪的ifra频响曲线去噪方法。

背景技术：

电力变压器是电力系统的核心设备之一，其状态影响着电力系统运行的安全性和可靠性。相关统计数据表明，绕组变形是大量变压器事故的重要起因之一。目前，绕组变形检测仅局限于离线状态，然而电网容量的扩增和智能电网的建设对绕组变形检测提出了在线运行的要求。ifra法具有对系统影响小，测试时间短，频带易拓展的特点，得到了众多研究人员的支持。然而目前ifra法还存在频率分辨率低等缺点。

频率响应分析法(frequencyresponseanalysis，fra)主要通过分析故障前后频率响应曲线谐振点的位置和幅值的变化对故障进行诊断。ifra是一种频率响应分析法，如果ifra频响曲线的频率分辨率低，即谱线间隔较大，有可能导致谐振点的遗漏，从而影响绕组变形故障诊断的准确性。

然而ifra频响曲线频率分辨率的提升必然导致其信噪比的降低。这是由于为了得到ifra频率响应曲线，首先要对离散的时域采样信号进行快速傅立叶变换以得到其频谱，受离散采样的限制，该频谱的谱线间为采样频率与采样总点数之商，在为了保障时域采样精度采样频率不可减小的前提下，若要减小谱线间隔，提升频率分辨率，只能增加采样总点数。然而，采样总点数的增加，必然导致信息冗余，这将降低ifra频率响应曲线的信噪比，使得其谐振点位置和幅值特征不清晰。

因此，在得到高频率分辨率ifra频响曲线后可根据其特性对其进行去噪处理，使得其频率分辨率的提升具有实际意义。ifra频响曲线通常具有低频段(1khz～1mhz)谐振点密集，信噪比高，高频段(≥1mhz)谐振点稀疏，信噪比低的特点。谐振点位置和幅值的变化是进行绕组变形诊断的重要依据，去噪时应尽可能避免对其造成影响。因此不能笼统地对整条ifra频响曲线进行去噪处理，本方法提出按照谐振点的分布将其分为多个频段，对ifra频响曲线的不同频段采用不同强度的去噪方法进行去噪。本方法提出采用调整小波阈值去噪过程中小波分解的层数对去噪强度进行控制。

技术实现要素：

本发明目的在于克服现有技术的不足与缺点，提出了一种基于小波阈值去噪的ifra频响曲线去噪方法，该方法在保留谐振点位置和幅值信息的前提下有效提升了ifra频响曲线的频率分辨率，提升了ifra绕组变形检测的精确度，具有实际应用意义。

为实现上述目的，本发明所提供的技术方案为：一种基于小波阈值去噪的ifra频响曲线去噪方法，包括以下步骤：

1)根据ifra频响曲线上谐振点分布密度将其分为若干频段；

2)找到ifra频响曲线上各个频段的最适分解层数，所述最适分解层数定义为：在对一条ifra频响曲线进行小波阈值去噪过程中，若采用某层数作为小波分解的层数，去噪后ifra频响曲线在某频段的噪声得到有效滤除，曲线毛刺得以消除且95％以上的谐振点不被滤除，则称该层数为该频段的最适分解层数；其中，一条ifra频响曲线上某频段的最适分解层数可能有多个，取其一即可，通过多次尝试找到各个频段的最适分解层数；

3)对含噪的ifra频响曲线整体进行多次小波阈值去噪，去噪的次数与步骤1)分的频段数相同，并依次采用步骤2)找到的各个频段的最适分解层数作为去噪过程中的小波分解层数进行小波阈值去噪，获取每次小波阈值去噪后能够得到有效去噪的频段曲线，即一整条含噪的ifra频响曲线需进行与步骤1)分的频段数相同数量的小波阈值去噪次数，且每次小波阈值去噪用的是步骤2)找到的所有频段当中某一频段的最适分解层数作为去噪过程中的小波分解层数；

4)将步骤3)得到的多条有效去噪的频段曲线依次拼接，即可得到所需的一整条完整的最终去噪后的ifra频响曲线。

本发明与现有技术相比，具有如下优点与有益效果：

1、本发明首次解决了ifra频率响应曲线频率分辨率提升会使其信噪比降低，从而影响频率分辨率提升的实际效用的矛盾，使得ifra检测频率精确度得以有效提升。

2、本发明首次实现了在保留谐振点位置和幅值信息的前提下提升了ifra频响曲线频率分辨率。

3、本发明操作步骤简单、快捷，成本低，在不需更新ifra检测硬件的基础上提升了检测的精确度。

4、本发明有利于将ifra检测法推广到工业实践中，有利于变压器绕组变形在线检测平台的搭建。

附图说明

图1为本发明的ifra频响曲线去噪方法流程图。

图2为去噪前ifra频率响应曲线。

图3为去噪后ifra频率响应曲线。

具体实施方式

下面结合具体实施例对本发明作进一步说明。

参见图1所示，本实施例所提供的基于小波阈值去噪的ifra频响曲线去噪方法，包括以下步骤：

1)将含噪ifra频响曲线分为若干频段。参见图2所示，含噪的ifra频响曲线在1khz～2.9mhz的范围内，谐振点分布密度为高，在2.9mhz～5mhz范围内谐振点分布密度为中，而在5mhz～10mhz范围内，谐振点分布密度为低。因此，将ifra频响曲线分为3个频段，频段1为1khz～2.9mhz，频段2为2.9mhz～5mhz，频段3为5mhz～10mhz。

2)找到各个频段的最适分解层数。在matlab中编写小波阈值去噪程序，载入含噪ifra频响曲线，小波分解层数初始值设置为3。针对频段1，在小波分解层数初始值设置为3时，去噪后频响曲线已经满足噪声得到有效滤除，曲线毛刺得以消除且95％以上的谐振点不被滤除的要求，因此设置频段1最适分解层数l1为3。改变小波分解层数，找到频段2和频段3的最适分解层数l2、l3分别为5、7。

3)对含噪的ifra频响曲线整体进行3次小波阈值去噪(其实去噪的次数要与步骤1)分的频段数相同)，然后依次采用步骤2)找到的3个频段的最适分解层数作为去噪过程中的小波分解层数进行小波阈值去噪(即该一整条含噪的ifra频响曲线每次小波阈值去噪用的是步骤2)找到的3个频段当中某一频段的最适分解层数作为去噪过程中的小波分解层数，且进行3次小波阈值去噪时，步骤2)找到的3个频段的最适分解层数都会用上)，获取每次小波阈值去噪后能够得到有效去噪的频段曲线(每次去噪实质得到的也仅有一段频响曲线)。

4)将步骤3)得到的3条有效去噪的频段曲线依次拼接，即可得到所需的一整条完整的最终去噪后的ifra频响曲线，具体参见图3所示。

以上所述实施例只为本发明之较佳实施例，并非以此限制本发明的实施范围，故凡依本发明之形状、原理所作的变化，均应涵盖在本发明的保护范围内。

技术特征：

技术总结
本发明公开了一种基于小波阈值去噪的IFRA频响曲线去噪方法，包括步骤：1)将IFRA频响曲线分为若干频段；2)找到IFRA频响曲线上各个频段的最适分解层数；3)对含噪的IFRA频响曲线整体进行多次小波阈值去噪，去噪的次数与步骤1)分的频段数相同，并依次采用步骤2)找到的各个频段的最适分解层数作为去噪过程中的小波分解层数进行小波阈值去噪，获取每次小波阈值去噪后能够得到有效去噪的频段曲线；4)将步骤3)得到的多条有效去噪的频段曲线依次拼接，即可得到所需的一整条完整的最终去噪后的IFRA频响曲线。本发明在保留谐振点位置和幅值信息的前提下有效提升了IFRA频响曲线的频率分辨率，提升了IFRA绕组变形检测的精确度。

技术研发人员：唐文虎;黄晶晶;周文婷;周九江;辛妍丽
受保护的技术使用者：华南理工大学
技术研发日：2017.05.05
技术公布日：2017.11.03