本发明涉及gis绝缘缺陷故障诊断领域,具体涉及一种电力设备绝缘缺陷模式识别方法。
背景技术:
封闭式组合电器(gasinsulatedswitchgear,gis)是一种以sf6气体作为绝缘介质的气体绝缘金属封闭开关设备。与传统的敞开式变电站相比,gis具有占地面积小、运行可靠性高、安全性强、维修周期长等显著优点。因此自20世纪60年代实用化以来,gis在国内外电力系统中得到广泛应用。
由于gis在电网中应用的广泛性及其重要程度,它的运行情况与整个电网能否正常安全工作息息相关。对gis局部放电进行在线监测则可以在掌握其绝缘情况的同时可以避免停运带来的负面影响,并且在不停运的状态下检测更能够表征gis当前的绝缘状态,对于保证整个电力系统的安全稳定运行具有重要意义。超声检测是目前常用的一种局部放电检测方法,但是基于超声检测局放信号对绝缘缺陷进行识别的方法准确率较低,无法实现gis绝缘缺陷类别的精确诊断。
技术实现要素:
针对现有技术中的不足,本发明的目的是为克服现有技术的不足,提供了一种gis局放故障诊断方法,该方法基于mfcc特征图谱和omp稀疏重构,实现gis绝缘缺陷的智能识别,提高gis局部放电检测系统的智能化水平。
为实现上述目的,本发明的技术方案如下:
一种电力设备绝缘缺陷模式识别方法,其包括:
步骤1:采集gis局部放电的超声信号;
步骤2:从超声信号中提取mfcc特征参数,并构造mfcc二维图谱;
步骤3:将样本按比例划分为训练集与测试集,使用训练集构建过完备字典,并根据建立的过完备字典采用omp法对测试集进行稀疏重构;
步骤4:依据非零项的累积误差判断放电类型。
如上所述的电力设备绝缘缺陷模式识别方法,进一步地,构造mfcc二维图谱具体包括:
步骤2-1:对信号分帧加窗:取帧长为一个工频周期的长度即20ms,重叠率为50%,并为信号加汉明窗以防止在后续过程中发生频谱泄露;
步骤2-2:对加窗后信号进行fft变换;
步骤2-3:使用mel滤波器组过滤信号,mel滤波器组的频率响应如下:
式中,f(m)为第m个三角滤波器的中心频率,共有m个三角滤波器;
步骤2-4:对mel频谱进行对数能量处理,得到对数频谱;
式中x(k)为分帧信号的fft,hm(k)为mel滤波器频率响应,n为fft长度;
步骤2-5:将对数频谱进行离散余弦变换,得到m个mfcc特征参数;
步骤2-6:得到信号的mfcc特征参数后,依据分帧顺序排列构成特征矩阵,最终得到信号的二维mfcc特征图谱。
如上所述的电力设备绝缘缺陷模式识别方法,进一步地,omp法对测试集进行稀疏重构具体包括:
omp的输入量有特征量矩阵x,传感矩阵φ,采样向量y,稀疏度k;输出量为x的k-稀疏的逼近
步骤3-1:初始化参数,残差r0=y,索引集λ0=φ,迭代次数t=1,最大迭代次数m;
步骤3-2:通过计算残差r和传感矩阵的列
步骤3-3:由λt更新索引集λt=λt-1∪{λt},记录找到的传感矩阵中的重建原子集合
步骤3-4:由最小二乘法得
步骤3-5:更新残差
步骤3-6:循环步骤3-1到步骤3-5,直到判断迭代次数是否满足t≥m,若满足,则停止迭代,获得非零项累计误差即最终残差值,残差值最小的类型即为待识别样本的绝缘缺陷类型;若不满足,则返回执行步骤3-1。
本发明与现有技术相比,其有益效果在于:本发明首先提取超声局放信号的mfcc声纹特征参数,并构造形成mfcc特征图谱,其次使用样本构建过完备字典,并基于已建立的过完备字典采用omp法对待测样本进行稀疏重构,最后依据非零项的累积误差判断放电类型,实现gis绝缘缺陷的智能识别,提高gis局部放电检测系统的智能化水平。。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例中所需要使用的附图进行简单的介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为基于omp稀疏重构的缺陷识别流程图;
图2为不同绝缘缺陷下gis局放超声信号,其中,图2(a)为金属微粒,图2(b)为悬浮电位;图2(c)为沿面放电;图2(d)尖板放电;
图3为不同绝缘缺陷下的mfcc特征图谱,其中,图3(a)为金属微粒,图3(b)为悬浮电位;图3(c)为沿面放电;图3(d)尖板放电。
具体实施方式
下面结合附图和具体实施方式对本发明的内容做进一步详细说明。
实施例:
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
本发明基于gis局放超声信号的特点,提出一种基于mfcc特征图谱和omp稀疏重构的gis绝缘缺陷诊断识别方法。首先提取超声局放信号的mfcc声纹特征参数,并构造形成mfcc特征图谱,其次使用样本构建过完备字典,并基于已建立的过完备字典采用omp法对待测样本进行稀疏重构,最后依据非零项的累积误差判断放电类型。通过对大量试验数据测试,验证了本文的识别方法在gis绝缘缺陷识别诊断的应用效果。
参见图1,图1为基于omp稀疏重构的缺陷识别流程图,该方法可以包括以下步骤:
步骤1:使用超声传感器采集gis局部放电的超声信号。
步骤2:从超声信号中提取mfcc特征参数,并构造mfcc二维图谱。
步骤3:将样本按比例划分为训练集与测试集,使用训练集构建过完备字典,并根据建立的过完备字典采用omp法对测试集进行稀疏重构。
步骤4:依据非零项的累积误差判断放电类型。
梅尔频率倒谱系数(melfrequencycoefficient,mfcc)是将实际频率转换到mel频率域中以获得声信号的声纹特性。梅尔频率与实际频率的关系为
式中,fmek为梅尔频率,f为实际频率。
建立mfcc特征图谱的具体流程如下:
步骤2-1:对信号分帧加窗:取帧长为一个工频周期的长度即20ms,重叠率为50%,并为信号加汉明窗以防止在后续过程中发生频谱泄露;
步骤2-2:对加窗后信号进行fft变换;
步骤2-3:使用mel滤波器组过滤信号,mel滤波器组的频率响应如下:
式中,f(m)为第m个三角滤波器的中心频率,共有m个三角滤波器;
步骤2-4:对mel频谱进行对数能量处理,得到对数频谱;
式中x(k)为分帧信号的fft,hm(k)为mel滤波器频率响应,n为fft长度;
步骤2-5:将对数频谱进行离散余弦变换,得到m个mfcc特征参数;
步骤2-6:得到信号的mfcc特征参数后,依据分帧顺序排列构成特征矩阵,最终得到信号的二维mfcc特征图谱。
正交匹配追踪法(omp)的基本思想是在每次迭代的过程中,从全息矩阵中选出与测量信号相关度(内积)最大的一列,将该列从全息矩阵中移除而加入到扩展矩阵中,然后利用最小二乘法原理求出使残差最小的估计,再不断从全息矩阵中减去相关列重复上述过程,直至达到规定的迭代次数或达到稀疏度要求。
omp的输入量有特征量矩阵x,传感矩阵φ,采样向量y,稀疏度k;输出量为x的k-稀疏的逼近
步骤3-1:初始化参数,残差r0=y,索引集λ0=φ,迭代次数t=1,最大迭代次数m;
步骤3-2:通过计算残差r和传感矩阵的列
步骤3-3:由λt更新索引集λt=λt-1∪{λt},记录找到的传感矩阵中的重建原子集合
步骤3-4:由最小二乘法得
步骤3-5:更新残差
步骤3-6:循环步骤3-1到步骤3-5,直到判断迭代次数是否满足t≥m,若满足,则停止迭代,获得非零项累计误差即最终残差值,残差值最小的类型即为待识别样本的绝缘缺陷类型;若不满足,则返回执行步骤3-1。
为验证本专利中所提出方法的有效性,实验室中搭建四种gis绝缘缺陷模型,分别是金属微粒放电、悬浮放电、尖板放电和沿面放电,并采集其超声局放信号,采集所得的原始信号如图2所示。
之后根据上述识别方法对采集所得超声信号进行分析诊断。提取每组信号的mfcc特征参数,形成mfcc特征矩阵后图像化该矩阵,所得mfcc特征图谱如图3所示。
随机选择70%样本为训练组,30%样本为测试组。使用训练集构建过完备字典,对测试集使用omp法在过完备字典上进行稀疏重构,依据最小残差判断种放电类型,最终识别结果统计如下表1所示,平均识别准确率达到91.67%。
表1omp稀疏重构绝缘缺陷识别结果
为体现所提方法的优越性,同时建立svm模型与ann模型对mfcc特征图谱进行识别,不同识别方法的识别准确率统计如下表2:
表2不同模型的识别结果
对比表2中不同识别模型的识别准确率结果可知,omp稀疏重构模型对于mfcc特征图谱识别准确率高于svm与ann,验证了本专利中所提方法的有效性。
本发明基于gis局放超声信号的特点,提出一种基于mfcc特征图谱和omp稀疏重构的gis绝缘缺陷诊断识别方法。首先提取超声局放信号的mfcc声纹特征参数,并构造形成mfcc特征图谱,其次使用样本构建过完备字典,并基于已建立的过完备字典采用omp法对待测样本进行稀疏重构,最后依据非零项的累积误差判断放电类型。通过对大量试验数据测试,验证了本文的识别方法在gis绝缘缺陷识别诊断的应用效果。
上述实施例只是为了说明本发明的技术构思及特点,其目的是在于让本领域内的普通技术人员能够了解本发明的内容并据以实施,并不能以此限制本发明的保护范围。凡是根据本发明内容的实质所做出的等效的变化或修饰,都应涵盖在本发明的保护范围内。