基于三比值特征量的小波PSO‑SVM变压器故障诊断方法与流程

本发明属于电气设备故障诊断方法
技术领域：
，具体涉及一种基于三比值特征量的小波pso-svm变压器故障诊断方法。
背景技术：
：油浸式变压器是电网的核心设备，承担着电压变换和电能传输等重任。变压器部分一旦发生故障可能导致电网严重故障，进而造成大停电等巨大损失。因此，如何准确判断变压器的运行状态，评估变压器的故障情况，降低变压器的故障风险，是电力企业急需解决的关键问题，也是保障电力系统稳定运行的重中之重。溶解气体分析法(dga)作为应用最广泛、最可靠的缺陷类型诊断方法之一，可有效区分变压器内部存在的缺陷类型。其原理主要是：变圧器在运行过程中，主要受到热应力及电应力的影响，变压器中的油纸绝缘系统逐渐受到破坏，进而分解产生co、co2及低分子烃类等相关产物，这些特征气体溶解于油中。因为不同性质故障所产生的气体含量和种类有所差异，所以能够根据变压器油中溶解气体的特征来判断故障类型。常用的故障表征气体包括：h2、ch4、c2h2、c2h4、c2h6、co和co2。dga油中溶解气体分析法不仅可有效诊断变压器内部的潜伏性故障，且方便在线监测。但是现有的dga诊断方法涉及特征量较多，且一些理化特征量(如油中溶解co、co2和糠醛含量等)作为评估变压器故障诊断的参考依据时受到现场换滤油、温度、适度等因素影响，有效性难以保证。在实际工程应用中多采用气体体积分数相对比值的方法进行诊断，如我国现行的dl/t722-2000《导则》中推荐的即是三比值法。比值法计算简单，判断法则直观,对操作人员要求不高，适合现场工作人员使用。因此引入三比值法对采集到的dga特征量进行处理和归一变换，以保证其有效性。但比值法存在编码边界过于绝对，编码不全等问题。目前变压器故障诊断的研究主要集中在基于已有的dga特征量(即dga比值，如roger比值、iec比值等)，采用支持向量机(svm)等技术建立变压器故障诊断模型。然而，不同文献中svm诊断模型采用的dga特征量有明显的区别，且部分dga气体比值不具备变压器故障诊断的能力，采用这些气体比值反而会影响故障诊断的效果。因此，亟需综合应用dga三比值特征量和svm模型的优点，进行变压器故障诊断，提升变压器故障诊断的准确率。技术实现要素：本发明的目的在于针对现有技术中的上述问题，提供一种基于三比值特征量的小波pso-svm变压器故障诊断方法。为实现上述发明目的，本发明采用了如下技术方案：基于三比值特征量的小波pso-svm变压器故障诊断方法，包括如下步骤：(1)采集变压器故障样本数据，选取变压器油中dga特征量，对dga特征量的比值进行分析；(2)对所述步骤(1)分析后的dga特征量的比值进行归一化预处理，得到归一化后的溶解气体分析三比值特征量数据；(3)构建非线性多分类的支持向量机模型，并求解优化函数；(4)采用小波核函数作为所述步骤(3)中支持向量机模型的核函数，通过所述步骤(3)的优化函数和小波核函数，计算出非线性多分类时支持向量机模型的分类决策函数；(5)构建粒子群优化算法以优化支持向量机核函数的参数；(6)采用所述步骤(5)的粒子群优化算法对所述步骤(4)中的核函数进行参数优化，构建小波pso-svm故障诊断模型；(7)利用dga特征量、待定惩罚因子和核参数的支持向量机模型形成故障诊断的目标函数，计算出目标函数的最优解，即为故障诊断的最优函数适应度；(8)计算所述步骤(7)中目标函数最优解所对应的支持向量机模型惩罚因子和核参数，即为最优的支持向量机模型惩罚因子和核参数；(9)将所述步骤(8)中最优的支持向量机模型惩罚因子和核参数代入所述步骤(7)中的故障诊断目标函数，构建小波核函数支持向量机诊断模型，采用该诊断模型对电力变压器油中溶解气体进行故障诊断，判断变压器的运行状态。上述步骤(1)中溶解气体分析特征量的比值包括c2h2/c2h4、ch4/h2和c2h4/c2h6三组气体的浓度比。上述步骤(2)中归一化预处理表达式为：式中，xsn为归一化计算后的溶解气体分析三比值特征量数据，xn为归一化前的溶解气体分析特征量比值，xnmax为归一化前溶解气体分析特征量比值的最大值，xnmin为归一化前溶解气体分析特征量比值的最小值。上述步骤(3)中的优化函数为：式中，ω为超平面的法向量，c为惩罚因子，ξi为松弛变量，l为松弛变量个数；同时满足以下约束条件：其中，设{(x1,y1),(x2,y2),...,(xn,yn)}的样本个数为n，样本xi(i＝1,2,...,n)对应的类别yi∈{-1,1}，为非线性映射，b为偏差量。上述步骤(4)中小波核函数的表达式为：式中，n为粒子群算法的种群数，a′为尺度因子。上述步骤(4)中的分类决策函数为：式中，ω为超平面的法向量，b为偏差量，jk表示多分类，为非线性映射，即核参数σ。上述步骤(5)中粒子群优化算法的数学表达式为：vsd(t+1)＝vsd(t)+c1(t)r1(t)(psd(t)-xsd(t))+c2(t)r2(t)(pgd(t)-xsd(t))；xsd(t+1)＝xsd(t)+vsd(t+1)；式中：t为进化代数，随机变量r1(t)、r2(t)为服从(0，1)区间的均匀分布，c1(t)、c2(t)为加速常数；xsd(t)为t代第s个粒子的位置；psd为t代群体内所有粒子搜索到的最优位置，vsd(t)为t代第s个粒子的速度，pgd(t)为t代群体内所有粒子搜索到的极值。上述步骤(1)中的dga特征量包括总烃含量、总烃产气率、油击穿电压和油介质损耗的特征量。相比于现有技术，本发明的优势在于：本发明所提供的基于三比值特征量的小波pso-svm变压器故障诊断方法，对所选变压器油中溶解气体(dga)特征量的比值进行分析，由归一化预处理，得到dga三比值特征量数据，实现了三比值特征量的有效性和简洁性；采用粒子群优化算法优化支持向量机的核函数参数，构建小波核函数支持向量机诊断模型，采用该模型对电力变压器的油中溶解气体进行诊断，从而判断变压器的运行状态并分析出相应的故障类型；本发明采用dga比值作为特征量进行故障诊断，并采用粒子群优化算法进行核函数的参数优化，提高了故障诊断的准确率；本发明将二分类svm拓展为多分类svm，能够对变压器的多种故障类型进行一次性的诊断识别，简单快捷，便于工程中的实际运用。本发明为判断变压器的运行状态、评估变压器的故障情况、降低变压器的故障风险、保障电力系统稳定运行提供了新思路。附图说明图1是本发明基于三比值特征量的小波pso-svm变压器故障诊断方法的流程图。图2是本发明油浸式变压器内dga气体分类图。图3是本发明最优函数的适应度迭代图。图4是本发明的小波核函数支持向量机诊断模型的训练准确率。图5是本发明的小波核函数支持向量机诊断模型的测试准确率。具体实施方式以下结合实施例及其附图对本发明技术方案进行清楚、完整地描述。所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。如图1所示，基于三比值特征量的小波pso-svm变压器故障诊断方法，包括如下步骤：(1)采集变压器故障样本数据，结合《国家电网公司状态评价导则》和老化特征量分析部分结果进行变压器油中dga特征量的选取，所选取的dga特征量包括总烃含量、总烃产气率、油击穿电压、油介质损耗等与变压器老化相关的特征量。采用三比值法对溶解气体分析特征量的比值进行分析：油浸式变压器受热应力及电应力的影响，会导致变压器油和绝缘纸(板)分解，产物包括h2、ch4、c2h2、c2h4、c2h6、co和co2等特征气体，如图2所示(油浸式变压器内dga气体分类)，该七种特征气体和总烃th(ch4、c2h2、c2h4、c2h6之和)采用三比值法对溶解气体分析特征量的比值进行分析；根据变压器在故障下分解产生气体组分含量的相对浓度与温度的相互依赖关系，从七种特征气体中选择两种溶解度和扩散系数相近似的气体组分组成三对比值数据进行分析，这三对比值数据分别为c2h2/c2h4、ch4/h2和c2h4/c2h6。(2)对步骤(1)分析后的溶解气体分析特征量比值进行归一化预处理，得到归一化后的溶解气体分析三比值特征量数据；归一化预处理的表达式如下式：式中，xsn为归一化计算后的溶解气体分析三比值特征量数据，xn为归一化前的溶解气体分析特征量比值，xnmax为归一化前溶解气体分析特征量比值的最大值，xnmin为归一化前溶解气体分析特征量比值的最小值。(3)构建非线性多分类的支持向量机模型，并求解优化函数minφ(ω,ξ)；式中，ω为超平面的法向量，c惩罚因子，ξi为松弛变量，l为松弛变量个数；同时满足以下约束条件：其中，设{(x1,y1),(x2,y2),...,(xn,yn)}的样本个数为n，样本xi(i＝1,2,...,n)对应的类别yi∈{-1,1}，为非线性映射，b为偏差量。本文中的t均指代矩阵转置符号。(4)对于式(2)和式(3)的qp问题(quadraticprogrammingproblem)，构建拉格朗日函数为：式中，αi和βi均为拉格朗日乘子，且αi>0和βi>0，有将式(5)代入拉格朗日函数式(4)中，得到如下对偶优化形式：由式(5)中c与αi的关系转换，得出计算分类问题的决策函数为：式中，k(x,x′)为小波核函数，采用小波核函数作为支持向量机模型的核函数：式中，n为粒子群算法的种群数，a′为尺度因子。采用oao(one-agains-onedecomposition)方法将二分类支撑向量机拓展为多分类支持向量机，即求解式(6)的优化问题，得出非线性多分类时支持向量机模型的分类决策函数为：式中，ω为超平面的法向量，b为偏差量，jk表示多分类，为非线性映射，即核参数σ。(5)构建粒子群优化算法：粒子x＝{x1,...,xs,...,xs}以速度vs＝[vs1,vs2,...,vsd]t(s＝1,2,...,s)在解空间(设在d的搜索空间)内飞行，且该粒子速度可按照其历史行为逐渐向其最优搜索位置动态调整(设第s个粒子搜索到的最优位置ps＝[ps1,ps2,...,psd]t)，其中psd为d维空间中，第s个粒子搜索到的最优解，vsd为d维空间中，第s个粒子搜索到的最优解)，粒子群优化算法的数学表达式为：vsd(t+1)＝vsd(t)+c1(t)r1(t)(psd(t)-xsd(t))+c2(t)r2(t)(pgd(t)-xsd(t))(10)xsd(t+1)＝xsd(t)+vsd(t+1)(11)式中：t为进化代数，随机变量r1(t)、r2(t)为服从(0，1)区间的均匀分布，c1(t)、c2(t)为加速常数；xsd(t)为t代第s个粒子的位置；psd为t代群体内所有粒子搜索到的最优位置，vsd(t)为t代第s个粒子的速度，pgd(t)为t代群体内所有粒子搜索到的极值。(6)支持向量机的惩罚因子c和核参数σ决定支持向量机的最优分类曲面即故障诊断准确率，因此，采用步骤(5)的粒子群优化算法对步骤(4)中的核函数进行参数优化。构建小波pso-svm故障诊断模型：粒子群优化算法首先随机生成初始例子x1,...,xi,...,xs，组成种群x，其粒子适应度值由f(xi1,x2,...,xin)来表示；小波pso-svm故障诊断模型表达式式中，an为运用svm算法分类时第n个验证集中的样本个数(n表示参加验证样本总个数，在算法中表现为最后一个)，为运用svm算法分类时该验证集中被正确分类的个数，e为交叉验证cv的折数，本实施例取值e＝5；f(xi1,x2,...,xin)中粒子适应度是评价粒子优劣的标准，粒子适应度越大越好，f(xi1,x2,...,xin)表现为f(c,σ)，f(c,σ)的值利用训练集和待定惩罚因子c、核参数σ不断迭代产生，当它达到最优时，对应了最佳的惩罚因子c、核参数σ，从而代入式(9)得出小波pso-svm故障诊断模型。(7)计算故障诊断的最优函数适应度：利用dga特征量、待定惩罚因子c、核参数σ的支持向量机形成故障诊断的目标函数，经过迭代、交叉等处理，计算出目标函数的最优解，即为f(xi1,x2,...,xin)中粒子适应度最大值。(8)计算步骤(7)中目标函数的最优解所对应的支持向量机模型惩罚因子c和核参数σ，即为最优的支持向量机模型惩罚因子c和核参数σ。(9)将步骤(8)中最优的支持向量机惩罚因子c和核参数σ代入步骤(7)中的故障诊断目标函数，得到小波核函数支持向量机诊断模型，并用测试集进行计算验证；采用该小波核函数支持向量机诊断模型对电力变压器油中溶解气体进行故障诊断，判断变压器的运行状态。实施例：本实例采用118组iectc10故障数据进行算法测试。将变压器故障划分为：低能放电(l-d)、高能放电(h-d)、中低温过热(l-t)、高温过热(h-t)、正常状态(n-c)五种状态，得到118组iectc10变压器故障样本数据如表1所示。表1变压器故障样本数据故障类型le-dhe-dlm-th-tn-ciectc10故障样本2345101426对搜集的118组dga数据进行三比值计算和归一化预处理得到归一化后的dga特征量(三比值特征量数据)，并通过不同的气体比值得到不同的dga特征量(iec三比值法数据，包括c2h2/c2h4，ch4/h2和c2h4/c2h6比值)。本发明提出的小波核函数支持向量机诊断模型中，惩罚因子c和核参数σ的搜索区间分别设为[1，103]和[10–3，10]，cv的折数e为9，最大迭代次数为10，初始生成的粒子群数为5，应用小波核函数支持向量机诊断模型对三比值特征数据进行故障诊断。对分析后的dga特征量比值进行归一化预处理，得到归一化后的溶解气体分析三比值特征量数据:构建非线性多分类的支持向量机模型：满足以下约束条件：构建拉格朗日函数为：式中，αi和βi均为拉格朗日乘子，且αi>0和βi>0，有由对偶优化得：得出决策函数为：其中小波核函数为：粒子群函数优化得到计算函数f2(c,σ)：由此，通过函数粒子适应度的评价优劣性，求得svm惩罚因子c和核参数σ，为非线性映射，即核参数σ。采用oao(one-agains-onedecomposition)方法将二分类支撑向量机拓展为多分类支持向量机，得出非线性多分类时支持向量机模型的分类决策函数为：pso每一次迭代的平均结果和最优结果如图3至图5所示。参数c和σ粒子群优化算法对svm故障诊断模型进行参数优化的适应度收敛曲线，由适应度值得到对应的最优参数。从图3中可以看到，适应度曲线在前期内收敛较快，而后经短暂波动，最终趋近一致收敛，即实现了参数的优化。通过图4得出的svm惩罚因子c＝100(即图4、图5中的参数bestc)和核参数σ＝77.7949(即图4、图5中的参数bestg)对实际的测试集进行上述步骤的计算，其结果由图5所示，其准确率为82.61％，这表明了该变压器故障诊断模型的准确性和有效性。得出的特征值的训练最大、最小、平均准确率，优选dga比值的故障诊断准确率如表2所示。另选取23组训练集中的10组以供校验，变压器故障诊断的准确率比较如表3所示，从表3可以看出，本发明的pso-svm模型的故障诊断准确率相对于其他方法更高。10组训练集校验数据如表4所示。表2表3表4以上所揭露的仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本
技术领域：
的技术人员在本发明揭露的技术范围内，可轻易想到变化或变型，都应涵盖在本发明的保护范围之内。当前第1页12