基于小波包的uhfpd信号特征提取

基于小波包的uhfpd信号特征提取

1uhfpd信号的分解最常见的电气故障特点是，在完全穿透或闪络前，设备中产生局部放电（pd）。采用超高频(UltraHighFrequency,UHF)检测GIS内PD可以及时发现设备内部缺陷。对局部放电信号进行有效模式识别,可以进一步准确地了解和掌握GIS内缺陷类型性质和特征,对于保证GIS安全可靠运行,掌握GIS绝缘状况及指导GIS检修工作有着十分重要的意义。特征提取品质的优劣极大地影响着分类器的设计和性能,是模式识别最重要的核心问题之一。在PD模式识别中,由于图像或者波形所获得的数据量相当大,对放电模式直接识别是很困难的。为了有效地实现UHFPD分类识别,就要对原始数据进行变换,得到最能反映分类本质的特征。小波包变换(WaveletPacketTransform,WPT)在满足海森堡测不准原理的前提下,自由选择不同时间点、频率点上的时频分辨率,可以更精练地在时频联合分析相平面上提取非稳定信号的特征。20世纪90年代末,BernhardScholkopf小组在主分量分析的基础上,将核主分量分析(KernelPrincipalComponentAnalysis,KPCA)推广到非线性领域,它通过核函数将输入空间的非线性问题转化到特征空间成为普通的特征值问题,比主成分分析具有更优越的降维效果。本文从WPT对UHFPD信号分解过程入手,根据已建立的GIS内4种典型缺陷UHFPD数学模型,分别采用熵最小原则选取最优小波包基,综合选用了UHFPD信号在各频带投影序列的能量、在各个尺度下的模极大值和绝对平均值,构造出能完整描述UHFPD信号的特征空间;利用KPCA法将高维特征空间降到低维特征空间,解决了维数危机,消除了类内散度矩阵的奇异性,并最大限度地保持原有信号的特性。研究结果表明,基于熵最小原则的最优小波包变换和KPCA法获得的信号特征均能够较好的应用于UHFPD信号模式识别。2uhfpd的物理模型作者根据GIS设备绝缘缺陷放电形式和特点,设计了4种GIS模拟装置内局部放电物理模型:分别是:高压导体金属突出物缺陷,简称N(needle)类放电模型缺陷;自由金属微粒缺陷,简称P(particle)类放电模型;绝缘子表面固定金属微粒缺陷,简称M(metal)类放电模型;气隙缺陷,简称G(gap)类放电模型。将上述4种人工物理缺陷模型分别置于研制的GIS模拟装置中,并充以0.5MPa的SF6与N2的混合气体(体积比4∶1),通过研制的内置超高频传感器和高速数字采集系统(模拟带宽1GHz、采样率为20GS/s、存储深度48MB)获取UHFPD信号的数据和波形。实测放电波形如图1所示。3基于wpt的高频局部压电位特性的空间结构3.1双尺度方程法小波包分解是一种更加精细的离散小波变换,对频带进行多层次划分,对每一节点系数采用近似值矢量分离法分解成两部分,从而得到完整的二进树。双尺度方程可以定义为式中,h(k),g(k)为共轭滤波器组;W1,0(t)为尺度函数;W1,1(t)为小波函数。每一节点(j,n)的小波包系数(WaveletPacketCoefficients,WPC)由时间序列x(t)计算得到WPC{Cj,n(k)}包含了重构信号的所有信息,其平方值代表了小波包树特定子带中信号的能量。3.2shrenon熵和上as-pb的编码及实验结果本文采用db系列小波包对UHFPD信号进行分解,分解层数4层,最优小波基的选择采用熵最小原则,且花费函数为Shannon熵。对于离散的时间序列x=(x1,x2,,xN),Shannon熵定义如下:当xi=0,则有(xi)2log(xi)2=0,M(x)反映了时间序列能量集中的程度,其最优小波包基快速搜索算法如下:(1)对已进行归一化处理了的UHFPD信号的数学模型,采用小波包函数进行4层小波包分解,并计算出各节点的Shannon熵值。(2)对最低层的各节点作标记,而对其他各层节点不作标记。(3)将每层从左至右每两个节点构成一组,由下往上进行相邻两层的下一层每组的两个节点之和与上一层对应节点的熵值比较。(4)若下一层的两个节点的Shannon熵和大于等于上一层对应节点的Shannon熵值,则对上一层的该节点作标记,而去掉下一层的该两个节点的标记;否则,以下一层该两个节点的Shannon熵之和取代上一层的该节点的Shannon熵值,但对标记不作处理。(5)从下往上对相邻两层按上述过程逐层进行比较,直至所有层均处理完毕。被作上标记的所有节点所对应的小波包函数即是搜索到的最优小波包基。针对GIS中4种常见缺陷研制的物理模型,作者采用5阶高斯函数构建了相应的数学模型,其数学表达式为式中,参数ai表示波峰的高度;bi表示波峰所在位置的横坐标x的值;ci反映了波峰的陡度;n表示放电脉冲极值的个数;对应参数值参见文献。作者选用了db系列小波构成的小波包并采用Shannon熵对GIS内4种典型缺陷对UHFPD数学模型进行4层分解后得到了最优小波包基。图2分别表示了以db4小波包为例经分解后得到的最优小波包基。从图中可以看出,对于不同的信号,即使采用同一小波包进行分解,所得到的最优小波包基也可能是不同的。UHFPD信号经过快速最优小波包基选取后,使信号的能量更加集中,在“最佳”频率分片上得到一系列的小波包系数,这些系数完备地描述了信号的特征,因而可以用作分类的特征子集。信号经m层小波包分解后,原UHFPD信号的能量被分解到2m个正交频带上,信号在各频带上的能量总和与原信号的能量一致,每个频带内的UHFPD信号表征原信号在该频率范围内的放电信息。因此,可将UHFPD信号在各频带投影序列的能量或与能量对应的值作为特征矢量。另外,Mallat证明了当各个尺度下的模极大值已知时,可利用投影迭代法恢复各尺度的子波变换值,从而完全实现原信号的重建。这表明利用各个尺度下的极大值作特征具有信息不丢失的优点,因而有望取得较高的识别能力。此外,每个尺度下小波包系数的统计特征也是工程上常用的特征之一。在本文的研究中,综合选用了UHFPD信号在各频带投影序列的能量、在各个尺度下的模极大值和绝对平均值完整地构造了其特征空间。但是该特征空间的特征维数仍然非常多,如果都作为识别的特征量,势必严重降低分类器的性能,而且很不适合实时应用的场合,因此必须进行降维处理,以除去那些对分类没有太多贡献的特征,同时可以用较低维数的向量作为分类器的输入向量,从而简化分类器的设计,并提高其识别能力。4kpca对超级高频局部放炮信号的性能进行了二维处理，并确定了识别结果4.1kpca基本思想PCA是一种基于目标统计特性的最佳正交变换,具有提取模式中具有最大描述特征的能力,然而PCA是一种线性算法,只考虑数据中的二阶统计特性,当UHFPD信号特征存在着大量非线性关系时,已经不能满足要求。因此,本文采用PCA与核学习方法的有机融合而形成的KPCA不仅特别适合处理非线性问题,而且能够提供更多的信息。KPCA的基本思想是通过一个非线性映射,将输入数据映射到一个高维特征空间,再在特征空间上进行线性主成分分析。设输入数据xk被映射为Φ(xk)假设它们已被中心化,即满足式中,xk(k=1,,N)是N个输入训练样本;Φ(xk)是变换后的训练样本。映射后训练样本的协方差矩阵C为求解下列特征方程:根据再生核理论,特征向量v一定位于由Φ(x1),…,Φ(xN)张成的空间内,即v可以由Φ(x1),…,Φ(xN)的线性组合表示:式中,α1,…,αi为常数。定义一个N×N的矩阵K称K为核矩阵。将式(7)、式(9)和式(10)代入式(8)得这样求解式(8)的特征向量v的问题就转化为求特征方程(11)的特征向量α。由式(10)可知核矩阵K是对称、半正定方阵,且它的特征值将是非负的。通过对特征方程(11)的求解,可得到一组非零特征值λj以及对应的满足归一化条件(12)的特征矢量αj(j=1,…,N)。根据式(9),可得到特征空间上的投影主分量vj=(j=1,…,N)。设x是一检验样本,则其在vj上的投影为这样就可得到原始信号的投影特征矢量。4.2多尺度归一化建模技术采用第2节所述试验得到的UHFPD信号为样本,该样本由4种放电缺陷模型在不同条件下测得,每类同种条件下均有100个有效放电样本,取其中50个样本作为训练样本,其余50个为测试样本。因此,N类样本数为3200,P类样本数为1200,M类样本数为1600,G类样本数为1600。将归一化处理后训练样本按照第3节介绍的方法经最优小波包变换后,得到一系列的小波包系数,计算出UHFPD信号在各频带投影序列的能量、在各个尺度下的模极大值和绝对平均值,构造出其特征空间分别用x1,…,xN来表示。然后利用第4.1节的KPCA法从特征空间中抽取核主成分特征,采用高斯核函数,其表达式为式中,参数σ=0.3×N,N为原始输入空间的维数。图3为4种GIS典型缺陷归一化的最大特征值递减示意图。图中n代表第n个特征量。通过计算可得前20个特征值对应的累计方差贡献率达到了87.5%以上,这样既减少了计算机的存储量,又大大缩短计算时间。然后将训练样本向该特征向量构成的空间投影,这些投影值构成原测试样本经KPCA提取后的新特征向量,作为分类器的输入进行模式识别。4.3uhfpd信号分解识别率按照上述方法和步骤采用db系列最优小波包4层分解对UHFPD信号提取的特征向量进行放电模式识别,分类器选用了径向基函数神经网络,训练方法为正交最小二乘法,结果如下表所示。由上表看出,识别率均在75%以上,其中采用db3、db4最优小波包对UHFPD信号分解识别率最高,均在85%以上。这是因为:当支集长度不消失矩太低时,不利于信号能量的集中;当支集长度太长时则会产生边界问题。另外,N类和M类缺陷识别率高于P类和G类缺陷的识别率。这是由于对于自由金属微粒缺陷(P类)存在的状态比较复杂,随着外加电压的变化,自由金属微粒缺陷有可能处于静止、起跳、舞动等不同运动状态,因此测量数据具有一定的分散性;对于气隙缺陷(G类)来说,主要是由于产生气隙缺陷的绝缘子在多次放电后,有数次击穿并伴有轻微炭化现象,因此导致对该缺陷所采集的波形发生较大的变化,从而导致识别率下降。5uhfpd信号特征(1)根据UHFPD的特点,提出一种适合UHFPD模式识别特征提取的新方法。在GIS内典型缺陷UHFPD数学模型上,分别采用熵最小原则选取最优小波包基,综合选用了UHFPD信号在各频带投影序列的能量分部特征、在各个尺度下的模极大值和绝对平均值,构造出能完整描述UHFPD信号的特征空间,并用KPCA法从特征空间中成功提取了GIS内UHFPD信号特征量。(2)采用Shannon熵最小原则法