一种改进樽海鞘群优化支持向量机的变压器故障诊断方法

1.本发明涉及变压器故障诊断技术领域，具体涉及一种改进樽海鞘群优化支持向量机 的变压器故障诊断方法。


背景技术：

2.变压器是电力系统中关键的变电装置，其状况能够迅速联系到整个电网的稳定性。 当其发生故障时不能及时诊断，会使电力系统正常安全运行得不到保证，会影响对后续 设备的供电，带来不可预计的经济损失。因此如何确保故障诊断的高精准度意义重大.
3.当油浸式变压器突发内部故障，通常反映为油色谱在线监测气体异常升高，溶解气 体分析(dis-solved gas analysis，dga)已经成为变压器故障诊断最常用的方法，dga 主要用于分析溶于变压器油中的气体含量，依据不同气体含量判断变压器故障类型。近 年来，相关学者提出了采用dga方法得到的溶解气含量比值作为特征量，据此诊断变 压器的故障类型。目前，故障诊断大多依据工作人员现场经验，通过溶解气含量或比 值判断运行状态和故障类型，导致工作人员经验对故障诊断影响较大。油浸式变压器中 主要包括油裂解气体和绝缘纤维素裂解气体(氢气(h2)、总烃(th)、一氧化碳(co) 和二氧化碳(co2)，其中总烃为甲烷(ch4)、乙炔(c2h2)、乙烯(c2h4)及 乙烷(c2h6)含量之和)。相关学者研究表明采用dga气体含量作为特征量时，会 影响油浸式变压器故障诊断的准确性，主要是因为相同油浸式变压器故障类型，dga 气体含量变化范围较大。不同学者采用的dga气体含量比值作为特征量存在区别，且 未形成统一的标准。
4.随着人工智能、机器学习和数据挖掘等技术日渐成熟,各种机器学习方法被一一应用 到变压器故障诊断当中,如采用支持向量机(svm)等技术建立变压器故障诊断模型。 然而不同文献中svm诊断模型采用的特征量有明显的区别，若采用过多的特征量作为 故障诊断模型的输入特征量会引起输入特征量冗余及干扰故障诊断结果。因此利用主成 分分析法(pca)对输入特征量进行降维处理达到降低特征量冗余的目的再结合svm 模型的优点，进行变压器故障诊断，提升变压器故障诊断的准确率。
5.目前svm参数c和g的优化选取，国际上并没有公认统一的最好的方法，随着群 体智能优化算法的出现，樽海鞘群算法自2017年被提出来就被广泛使用，很大程度上 解决了svm在变压器故障诊断中存在的训练时收敛速度慢、诊断精度低等问题。但原 樽海鞘群算法存在仅有一个食物源导致的陷入局部最优的问题。


技术实现要素：

6.因此，本发明要解决的技术问题在于克服现有技术中的缺陷，从而提供一种改进樽 海鞘群优化支持向量机的变压器故障诊断方法。
7.本发明提供了一种改进樽海鞘群优化支持向量机的变压器故障诊断方法，包括：
8.s1：采集变压器故障样本，获取数据；所述数据包括变压器油中dga气体以及dga 气体含量比值；
9.s2：对所述数据进行归一化预处理，得到优选特征量数据；对优选特征量数据进行 特征优选；
10.s3：采用径向基函数作为核函数，基于径向基函数和目标函数得到分类决策函数，； 根据所述分类决策函数和核函数构建支持向量机模型；
11.s4：根据改进的樽海鞘群算法优化所述支持向量机模型，得到优化后的惩罚因子以 及优化后的核参数；
12.改进的樽海鞘群算法包括：
13.步骤1：种群随机初始化；
14.步骤2：计算适应度的值；
15.步骤3：将樽海鞘群分为领导者和追随者，对领导者的位置以及追随者的位置进行 更新；
16.步骤4：领导者组成领导者群体，计算领导者群体的加权平均位置；并对追随者的 舒适度进行排序，确定适应度前三的追随者；领导者的位置在更新时，随机选取所述适 应度前三的追随者中一个追随者更新食物源；根据更新后的食物源再次更新领导者的位 置；
17.步骤5：判断终止条件，终止条件包括迭代次数是否达到预设的迭代次数或适应度 的值不再提高；若满足终止条件，则输出再次更新后领导者的位置；若不满足则返回步 骤3；
18.s5：基于优化后的支持向量机模型、特征优选后的优选特征量数据、优化后的惩罚 因子以及优化后的核参数构建变压器诊断模型，根据变压器诊断模型对变压器进行故障 诊断。
19.优选的，步骤3中，对领导者的位置进行更新的计算公式为：
[0020][0021][0022]
其中，表示领导者的位置；fj表示更新前的食物源；c1为收敛因子c2，c3为[0,1] 内随机产生的数字；
[0023]
对追随者的位置进行更新的计算公式为：
[0024][0025]
其中，为更新之前第r个追随者在第j维的位置；为更新之前第r-1个追随者 在第j维的位置。
[0026]
优选的，步骤4中，更新后的食物源记为：
[0027][0028]
[0029][0030]
其中，分别为适应度前三的追随者；表示领导者群体 的加权平均位置；n表示樽海鞘群种群数量，ωr表示领导者种群中按适应度值降序排列 的权重系数；
[0031]
根据更新后的食物源再次更新领导者的位置，计算公式为：
[0032][0033]
其中，表示领导者的位置；fj表示更新前的食物源；c1为收敛因子c2，c3为[0,1] 内随机产生的数字；
[0034]
再次更新后领导者的位置用于优化所述支持向量机模型。
[0035]
优选的，s5中，将特征优选后的所述优选特征量数据作为优化后的所述支持向量机 模型的输入，并将优化后的惩罚因子以及优化后的核参数输入至所述支持向量机模型中， 构建变压器诊断模型。
[0036]
优选的，对所述数据进行归一化预处理的表达式为：
[0037][0038]
其中，x
sn
表示优选特征量数据，xn表示dga气体含量比值，x
nmax
是归一化处理前 的数据的最大值，x
nmin
是归一化处理前的数据的最小值。
[0039]
优选的，s2中，采用主成分分析法对优选特征量数据进行特征优选。
[0040]
优选的，s3中，目标函数表达式为：
[0041][0042]
其中，ω表示超平面的法向量，ξi为松弛变量，c为惩罚因子，l为缓和变量个数； i表示第i个样本；
[0043]
目标函数满足约束条件，则得到分类决策函数；
[0044]
约束条件的表达式为：
[0045][0046]
其中，xi表示变压器故障样本特征量；yi表示变压器故障样本类别；ξi为松弛变量；表示非线性映射函数；b表示偏差量；ω
t
表示转置后的超平面的法向量，l为缓和 变量个数；
[0047]
分类决策函数的表达式为：
[0048]
[0049]
其中，ω为超平面的法向量，b为偏差量，jk表示多分类，为非线性映射，即 核函数σ；t表示转置。
[0050]
优选的，径向基函数的表达式为：
[0051][0052]
其中，x表示核函数中任意一点；xi表示核函数中心点；x、xi均为变压器故障样本 特征量；i表示第i个样本；σ表示核参数。
[0053]
本发明技术方案，具有如下优点：改进的樽海鞘群算法避免了原樽海鞘群算法中仅 有一个食物源导致的陷入局部最优的不足。另外，基于获取的数据，并且通过改进后的 樽海鞘群算法优化支持向量机；通过优化后的支持向量机构建变压器诊断模型；采用该 模型对变压器的油中溶解气体进行诊断，从而判断变压器中的运行状态并分析出相应的 故障类型。
附图说明
[0054]
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实 施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的 附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的 前提下，还可以根据这些附图获得其他的附图。
[0055]
图1为本发明实施中故障诊断方法的流程图；
[0056]
图2为本发明实施中改进的樽海鞘群算法的流程图。
具体实施方式
[0057]
下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施 例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通 技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的 范围。
[0058]
此外，下面所描述的本发明不同实施方式中所涉及的技术特征只要彼此之间未构成 冲突就可以相互结合。
[0059]
如图1所示，本实施例提供了一种改进樽海鞘群优化支持向量机的变压器故障诊断 方法，包括：
[0060]
s1：采集变压器故障样本，获取数据；所述数据包括变压器油中dga气体以及dga 气体含量比值；
[0061]
dga气体包括ch4、c2h2、c2h4、c2h6、h2共5个。dga气体含量比值备选特征 量有ch4/c2h6、ch4/h2、c2h
2/
ch4、c2h
2/
c2h4、c2h
2/
c2h6、c2h
2/
h2、c2h4/ch4、 c2h4/c2h6、c2h4/h2、c2h6/h2共10个。
[0062]
表1为变压器故障样本表；
[0063][0064]
油浸式变压器故障时，绝缘油会发生热裂解，产生h2、ch4、c2h2、c2h4和c2h6五种主 要气体并溶解于绝缘油中，不同故障类型下溶解的气体浓度和各气体之间的比例关系各 不相同。为了能完整的反映变压器故障类型特征，提高变压器故障诊断精度，选用五种 气体油中溶解相对浓度、两两气体间比例关系作为特征数据，具体特征数据如表2；
[0065]
表2为油中溶解气体比值表；
[0066][0067]
s2：对所述数据进行归一化预处理，得到优选特征量数据；对优选特征量数据进行 特征优选；
[0068]
具体的，对表2中的数据进行归一化预处理，得到优选特征量数据；采用主成分分 析法对优选特征量数据进行特征优选(利用主成分分析法(pca)对dga气体含量比 值进行降维分析，达到减少冗余信息的目的)；
[0069]
对数据进行归一化预处理的表达式为：
[0070][0071]
其中，x
sn
表示优选特征量数据，xn表示dga气体含量比值，x
nmax
是归一化处理前 的数据的最大值，x
nmin
是归一化处理前的数据的最小值。
[0072]
利用主成分分析法对优选特征量数据的比值进行特征优选，pca的综合特征量是原 始特征量经线性组合得到的，当综合特征量的累计贡献率能够达到80％以上时，则认为 综合特征变量就能够代替原始特征量，表征该工况特征。对前k个综合特征量的累计贡 献率可表示为：
[0073][0074]
其中，v＝1,2，
┈
，p；λv为第v个综合特征量的特征值；η为前k个综合特征量的 累计贡献率；λh为第h个综合特征量的特征值。
[0075]
s3：采用径向基函数作为核函数，基于径向基函数和目标函数得到分类决策函数，； 根据所述分类决策函数和核函数构建支持向量机模型；
[0076]
目标函数表达式为：
[0077][0078]
其中，ω表示超平面的法向量，ξi为松弛变量，c为惩罚因子，l为缓和变量个数； i表示第i个样本；
[0079]
目标函数满足约束条件，则得到分类决策函数；
[0080]
约束条件的表达式为：
[0081][0082]
其中，xi表示变压器故障样本特征量；yi表示变压器故障样本类别，yi∈{-1,1}；ξi为松 弛变量；表示非线性映射函数；b表示偏差量；ω
t
表示转置后的超平面的法向量；
[0083]
对于以上两式的二次规划问题(quadraiic programming problem，qp)，可借助拉格朗 日函数来表示：
[0084][0085]
其中，αi和βi为拉格朗日乘子，且αi》0和βi》0，于是得到
[0086][0087]
将上述代入拉格朗日函数中，得出的对偶优化形式为：
[0088][0089][0090]
由c与αi的关系转换，可得出计算分类问题的决策函数为，分类决策函数的表达式 为：
[0091][0092]
其中，k(x，xi)表示径向基函数；αi为拉格朗日乘子，b为偏差量；yi表示变压器 故障样本类别；i表示第i个样本；l表示缓和变量的数量。
[0093]
径向基函数的表达式为：
[0094]
[0095]
其中，x表示核函数中任意一点；xi表示核函数中心点；x、xi均为变压器故障样本 特征量；i表示第i个样本；σ表示核参数。
[0096]
采用oao(one-agains-onedecomposiiion)方法将二分类支持向量机拓展为多分类 支持向量机，得出非线性多分类时支持向量机模型的分类决策函数为：
[0097][0098]
其中，ω为超平面的法向量，b为偏差量，jk表示多分类，为非线性映射，即 核函数σ；t表示转置。
[0099]
s4：根据改进的樽海鞘群算法优化所述支持向量机模型，得到优化后的惩罚因子以 及优化后的核参数；
[0100]
如图2所示，改进的樽海鞘群算法包括：
[0101]
步骤1：种群随机初始化；根据如下计算公式对种群进行初始化：
[0102]
xn×d＝rand(n，d)
·
(ub-lb)+lb
[0103]
其中，搜索空间为n
×
d维，n为种群规模，d为观测量维数，ub为搜索上界，lb 为搜索下界。
[0104]
步骤2：计算适应度的值；计算公式为：
[0105][0106]
其中，i为支持向量机模型将训练集样本分类正确的数量，f为支持向量机模型将 训练集样本分类错误的数量。
[0107]
步骤3：将樽海鞘群分为领导者和追随者，对领导者的位置以及追随者的位置进行 更新；
[0108]
对领导者的位置进行更新的计算公式为：
[0109][0110][0111]
其中，表示领导者的位置；fj表示更新前的食物源；c1为收敛因子c2，c3为[0,1] 内随机产生的数字；
[0112]
对追随者的位置进行更新的计算公式为：
[0113][0114]
其中，为更新之前第r个追随者在第j维的位置；为更新之前第r-1个追随者 在第j维的位置。
[0115]
步骤4：领导者组成领导者群体，计算领导者群体的加权平均位置；并对追随者的 舒适度进行排序，确定适应度前三的追随者；领导者的位置在更新时，随机选取所述适 应度前三的追随者中一个追随者更新食物源；根据更新后的食物源再次更新领导者的位 置；
[0116]
更新后的食物源记为：
[0117][0118][0119][0120]
其中，分别为适应度前三的追随者；表示领导者群体 的加权平均位置；n表示樽海鞘群种群数量，ωr表示领导者种群中按适应度值降序排列 的权重系数；
[0121]
根据更新后的食物源再次更新领导者的位置，计算公式为：
[0122][0123]
其中，表示领导者的位置；fj表示更新前的食物源；c1为收敛因子c2，c3为[0,1] 内随机产生的数字；
[0124]
再次更新后领导者的位置用于优化所述支持向量机模型。
[0125]
步骤5：判断终止条件，终止条件包括迭代次数是否达到预设的迭代次数或适应度 的值不再提高；若满足终止条件，则输出再次更新后领导者的位置；若不满足则返回步 骤3。
[0126]
s5：基于优化后的支持向量机模型、特征优选后的优选特征量数据、优化后的惩罚 因子以及优化后的核参数构建变压器诊断模型，根据变压器诊断模型对变压器进行故障 诊断；
[0127]
具体的，将特征优选后的所述优选特征量数据作为优化后的所述支持向量机模型的 输入，并将优化后的惩罚因子以及优化后的核参数输入至所述支持向量机模型中，构建 变压器诊断模型。
[0128]
依照传统方法选取的输入特征作为dga特征量会陷入局部最优问题，无法全面放映 变压器的故障分类问题，在实际应用中极有可能造成故障诊断错误。若选用全部气体特 征量作为故障诊断模型输入特征量，会引起输入特征量冗余及干扰故障诊断结果。因此 基于主成分分析的dga气体含量比值特征量选择方法，当综合特征量的累计贡献率能够 达到80％以上时，则认为综合特征变量就能够代替原始特征量，表征该工况特征，达到 降低信息损失和降低特征量冗余的目的。
[0129]
本发明采用改进的海鞘群优化算法对支持向量机的核函数和惩罚因子进行了优化， 一定程度上避免了原算法中仅有一个食物源导致的陷入局部最优的不足。并将二分类 svm拓展为多分类svm，能够对变压器的故障特征量进行优选，提高了变压器故障诊断 的精确率以及诊断速度，减少了变压器故障诊断的时间和成本。
[0130]
显然，上述实施例仅仅是为清楚地说明所作的举例，而并非对实施方式的限定。对 于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化 或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引伸出的显而易见的变 化或变动仍处于本发明创造的保护范围之中。