CN102289682B - 基于集成学习Bagging算法的变压器故障诊断方法 - Google Patents

Abstract

本发明公开了变压器故障诊断技术领域中的一种基于集成学习Bagging算法的变压器故障诊断方法。本发明选用球向量机作为Bagging算法的弱学习算法，并使用油中溶解气体分析技术作为数据采集方法，通过数据归一化、类别数值化等数据处理方法得到适用于球向量机的样本集；集成学习Bagging算法通过反复调用弱学习算法完成对样本集的训练得到强学习机H；利用强学习机H作为变压器故障诊断模型，以完成待诊断记录的故障判别。本发明在改善变压器故障诊断精度方面，有较好的适应性和较低的诊断误差。

Description

基于集成学习Bagging算法的变压器故障诊断方法

技术领域

本发明属于变压器故障诊断技术领域，尤其涉及一种基于Bagging算法的变压器故障诊断方法。

背景技术

随着智能电网技术的飞速发展，传统的变压器故障诊断技术也发生巨大的变化，越来越多的智能化方法被引入到故障诊断领域。变压器作为一种常见电力设备，其结构相对复杂、故障种类繁多，如何快速、准确的诊断出变压器故障对于提高其运行维护水平和电力系统安全具有很强的理论及现实意义。目前变压器故障诊断方法主要是使用模式识别领域中的单分类器方法，将油中溶解气体分析技术DGA(Dissolved Gas Analysis)采集到的目标变压器中各个时间监测点的氢气H₂、一氧化碳CO、甲烷CH₄、乙烯C₂H₄、乙烷C₂H₆、乙炔C₂H₂和二氧化碳CO₂等七种气体的含量，及各个时间监测点对应的变压器运行状态(包括“正常、低温过热故障、中温过热故障、高温过热故障、局部放电故障、低能放电故障、低能放电兼过热故障、电弧放电故障、电弧放电兼过热故障”)作为数据样本进行训练，从而得到基于单分类器的变压器故障诊断模型，该过程主要包括三个部分，分别是构造样本集、分类器训练、故障诊断。

构造样本集是将变压器历史运行记录转化为模式识别算法可读的数据样本集，即将油中溶解气体分析技术DGA采集到的各个时间监测点的七种气体含量作为模式识别算法的输入向量，将各个时间监测点对应的变压器运行状态作为模式识别算法的输出向量。

分类器训练是将样本集输入到已选定的某种模式识别算法，通过对样本集的训练，得到精度满足要求的分类器，并将分类器作为变压器故障诊断模型。这是整个故障诊断过程最重要的一步。

故障诊断是将待诊断记录样本输入分类器，从而得到故障诊断结果。

目前变压器故障诊断方法，主要是利用决策树、神经网络、支持向量机等单分类器方法，而这些单分类器方法均为不稳定分类算法，所谓不稳定分类算法就是指训练样本集发生一个微小的变化，分类器的分类结果就会产生巨大变化。虽然经多年研究，通过各种群智能优化算法已使单分类器的稳定性有所提高，但其仍存在误差相对较大、易出现过拟合现象、计算过程复杂等难以克服的问题。为此本发明采用集成学习Bagging算法，并使用球向量机BVM作为Bagging算法的弱学习算法，完成变压器故障诊断，有效改善了诊断精度。

发明内容

针对上述背景技术中提到的基于单分类器的故障诊断技术因易出现过拟合现象而导致的诊断误差大、计算过程复杂等不足，本发明提出了一种基于集成学习Bagging算法的变压器故障诊断方法。

本发明的技术方案是，基于集成学习Bagging算法的变压器故障诊断方法，其特征是所述方法包括以下步骤：

步骤1：数据采集及数据预处理，得到初始历史数据样本集S；

步骤2：将初始历史数据样本集S构造成集成学习Bagging算法中弱学习算法可读的数据集S_sample，选定球向量机BVM作为弱学习算法；

步骤3：选定球向量机BVM作为集成学习Bagging算法中的弱学习算法并对该弱学习算法循环调用，完成样本集S_sample的训练，从而得到弱学习机序列，序列中包含各代弱学习机h_i，该序列即为强学习机H；

步骤4：将强学习机H作为故障诊断模型，并将当前待故障诊断的数据样本输入到强学习机H，强学习机H利用其各个弱学习机h_i做初步故障诊断，进而以投票的方式判定当前待测数据样本的故障类别，得票数多的故障类别即为强学习机H的最终诊断结果；

所述步骤1包括以下步骤：

步骤1.1：数据采集：采集目标变压器各个时间监测点的气体含量及运行状态；

步骤1.2：数据预处理：设定时间监测点的气体含量及运行状态的特征属性，构成初始历史数据样本集S。

所述步骤1.1中数据采集的方法为油中溶解气体分析技术DGA。

所述步骤1.1中气体包括氢气H₂、一氧化碳CO、甲烷CH₄、乙烯C₂H₄、乙烷C₂H₆、乙炔C₂H₂和二氧化碳CO₂。

所述步骤2包括以下步骤：

步骤2.1：数据归一化：将初始历史数据样本集S中气体的含量值按照各自的取值范围全部归一化到[0，1]区间；

步骤2.2：数值化处理：将变压器运行的运行状态设定为数值型类别标号；

步骤2.3：将初始历史数据样本集S中的气体含量作为故障诊断模型的输入向量；将初始历史数据样本集S中的变压器的运行状态作为故障诊断模型的输出向量，故障诊断模型的输入向量和输出向量构成了数据集S_sample。

所述步骤3具体包括下列步骤：

步骤3.1：对集成学习Bagging算法初始化，设定Bagging算法最大迭代次数t，并设定弱学习算法的训练参数；

步骤3.2：以指定概率从数据集S_sample中有放回地选取样训练本子集S_i，i∈[1，…t]，作为弱学习算法的训练样本子集S_i；

步骤3.3：将训练样本子集S_i输入到弱学习算法训练，得到对应的弱学习机h_i；

步骤3.4：检查当前集成学习Bagging算法是否达到算法的最大迭代次数t，若已达到，则执行步骤3.5；否则，返回步骤3.2；

步骤3.5：输出弱学习机序列，即强学习机H。

所述步骤3.1中弱学习算法的训练参数包括核函数类型、核函数参数和惩罚因子。

所述步骤3.2中的指定概率为50％。

本发明利用集成学习Bagging算法并行生成多个弱学习机来完成对目标问题的求解。采用本发明对变压器故障进行诊断，不但克服了原有基于单分类器故障诊断技术的缺陷，而且提高了诊断的准确率。

附图说明

图1为本发明方法的流程图；

图2为变压器故障诊断样本集的生成过程的流程图；

图3为集成学习Bagging算法训练弱学习机的流程图；

图4为强学习机诊断变压器故障过程的示意图。

具体实施方式

下面结合附图，详细说明本发明方法的实施过程。应该强调的是，下述说明仅仅是示例性的，而不是为了限制本发明的范围及其应用。

变压器运行数据采集是故障诊断的前提，只有准确、充分的数据样本才能保证故障诊断方法的有效性。为此，本发明采用油中溶解气体分析技术DGA(Dissolved Gas Ahalysis)收集目标变压器中各个时间监测点的氢气H₂、一氧化碳CO、甲烷CH₄、乙烯C₂H₄、乙烷C₂H₆、乙炔C₂H₂和二氧化碳CO₂等七种气体的含量，并记录各个时间监测点的变压器运行状态“正常、低温过热故障、中温过热故障、高温过热故障、局部放电故障、低能放电故障、低能放电兼过热故障、电弧放电故障、电弧放电兼过热故障”，通过数据归一化处理及数值化处理构造出变压器运行数据集。进而将数据集按照Bagging算法中弱学习算法要求改造成样本集，并将改造好的样本集输入集成学习Bagging算法，由集成学习Bagging算法调用弱学习算法进行样本集训练以得到各代弱学习机，再将各代弱学习机组合成强学习机。最后利用强学习机完成变压器故障诊断。这样从局部到整体，使集成学习Bagging算法能诊断更一般的变压器故障问题。

图1是本发明方法的流程图，可以看出，本发明方法包括如下步骤：

步骤1：数据采集及数据预处理，得到初始历史数据样本集S；

步骤2：将初始历史数据样本集S构造成集成学习Bagging算法中弱学习算法可读的数据集S_sample；

步骤3：选定球向量机BVM作为集成学习Bagging算法中的弱学习算法并对该弱学习算法循环调用，完成样本集S_sample的训练，从而得到一个弱学习机序列h＝{h₁，h₂，…h_t}，t为集成学习Bagging算法设定的迭代次数，弱学习机序列中的每一个h_i代表第i代弱学习机，i∈[1，…t]，弱学习机序列h＝{h₁，h₂，…h_t}即为强学习机H，j∈[1，+∞]；

步骤4：将强学习机H作为故障诊断模型，并将当前待故障诊断的数据样本输入到强学习机H，强学习机H利用其中各个弱学习机h_i做初步故障诊断，进而以投票的方式判定当前待测数据样本的故障类别，得票数多的故障类别为强学习机H最终诊断结果；

设已知某台变压器一段连续时间的运行状态记录，按照步骤1的方法进行数据采集和数据预处理：记录中包含的信息有时间监测点、“氢气H₂、一氧化碳CO、甲烷CH₄、乙烯C₂H₄、乙烷C₂H₆、乙炔C₂H₂和二氧化碳CO₂”七种气体含量的监测数据、以及变压器的运行状态“正常、低温过热故障、中温过热故障、高温过热故障、局部放电故障、低能放电故障、低能放电兼过热故障、电弧放电故障、电弧放电兼过热故障”，现按照图2的方法将设备记录改造成集成学习Bagging算法中弱学习算法可读的样本集。

首先，样本集中每条样本均包含9个特征属性，其中第一条为该样本的时间监测点，第2至第8条特征属性分别对应时间监测点和该时间监测点的氢(H₂)、甲烷(CH₄)、乙烷(C₂H₆)、乙烯(C₂H₄)、乙炔(C₂H₂)、一氧化碳(CO)、二氧化碳(CO₂)含量，第9条特征属性对应变压器当时的运行状态；

其次，执行步骤2.1：数据归一化，初始历史数据样本集S中氢(H₂)、甲烷(CH₄)、乙烷(C₂H₆)、乙烯(C₂H₄)、乙炔(C₂H₂)、一氧化碳(CO)、二氧化碳(CO₂)等七种气体的含量值按照各自的取值范围全部归一化到[0，1]的区间。比如H₂的含量最高值为100毫升/小时，最低值为0.1毫升/小时，则按下面公式进行归一化，其中x为未归一化之前的含量，为归一化之后的含量，由此可将所有含量值都约束在[0，1]的值域，

$\hat{x}=\frac{x-0.1}{100-0.1}$

再次，执行步骤2.2：数值化处理，变压器运行状态分为“正常、低温过热故障、中温过热故障、高温过热故障、局部放电故障、低能放电故障、低能放电兼过热故障、电弧放电故障、电弧放电兼过热故障”等九种状态，因而基于集成学习Bagging算法的变压器故障诊断模型共有九种诊断结果，分别对应变压器的九种运行状态，将九种运行状态设定为“-4、-3、-2、-1、0、+1、+2、+3、+4”等九个数值型类别标号。

最后，执行步骤2.3：将初始历史数据样本集S中氢(H₂)、甲烷(CH₄)、乙烷(C₂H₆)、乙烯(C₂H₄)、乙炔(C₂H₂)、一氧化碳(CO)、二氧化碳(CO₂)等七种气体的含量值作为故障诊断模型的输入向量，即集成学习Bagging算法中弱学习算法BVM的输入向量，将初始历史数据样本集S中变压器的运行状态作为故障诊断模型的输出向量，即集成学习Bagging算法中弱学习算法BVM的输出向量，至此已将初始历史数据样本集S构造成集成学习Bagging算法中弱学习算法可读的数据集S_sample，数据集S_sample中包含故障诊断模型的输入向量和输出向量。

将数据集S_sample输入到集成学习Bagging算法，以完成变压器故障诊断模型的训练。与其他集成学习方法不同的是：集成学习Bagging算法通过并行生成多个弱学习机完成对目标问题的求解。对于固定的初始样本集，集成学习Bagging算法采用有放回的方式每次以随机概率抽取相同个数的样本组成样本子集，并输入弱学习算法训练，从而得到弱学习机序列，该序列即为强学习机。最终分类判别时，根据各个弱学习机的判别结果投票决定待分类样本的类别归属。集成学习Bagging算法可以有效地提高泛化能力，因为其每次样本子集生成过程均为有放回的随机抽取Bootstrap Aggregating方法，因而各个弱学习机之间不存在依赖关系，保证了集成学习Bagging算法的可靠性。

集成学习Bagging算法的详细训练过程如图3，主要包括下列步骤：

步骤3.1：对集成学习Bagging算法初始化，选用球向量机BVM作为弱学习算法，并设定球向量机BVM的训练参数，如核函数类型、核函数参数、惩罚因子等；

步骤3.2：以指定概率从数据集S_sample中有放回的选取样本子集S_i，i∈[1，…t]，t为Bagging算法最大迭代次数，样本子集S_i作为球向量机BVM的训练样本集，具体选取规则采用Bootstrap Aggregating方法，即某个样本可能在不同子训练集中反复出现，也可能从不出现；

步骤3.3：将样本子集S_i输入到球向量机BVM中做训练，得到对应的弱学习机h_i，i∈[1，t]，t为算法最大迭代次数；

步骤3.4：检查当前集成学习Bagging算法的执行步骤是否达到迭代次数t，若已达到，则执行步骤3.5；否则，返回步骤3.2；

步骤3.5：输出弱学习机序列h＝{h₁，h₂，…h_t}，即最终强学习机H。

特别注意的是在步骤3.1中的训练过程是利用球向量机BVM完成对样本数据的训练。球向量机BVM改进于支持向量机SVM(Support Vector Machine)。球向量机BVM的改进在于利用最小包含球算法MEB(Minimum Enclosing Ball)代替了支持向量机SVM中的凸二次规划，从而大大节省了运算时间。最小包含球算法MEB算法理论基础来自于“计算几何(Computational Geometry)”。球向量机BVM利用最小包含球算法MEB算法求解原n维欧氏空间Rⁿ中目标问题Φ，其过程如下：

(1)将原n维欧氏空间Rⁿ中的目标问题Φ映射到Hilbert空间，并在Hilbert空间中构造对偶问题Φ′。

(2)根据对偶问题Φ′的样本集S构造初始球。

(3)迭代求解样本集S的核子集S_e，即完成对偶问题Φ′到最小闭包球MEB问题的转化。设c、r分别为初始球的重心和半径，使用B(c，r)表示一个重为c，半径为r的球，r∈[0，R]，r的上限为R，当r增加到R时，此时的球即为MEB球。再设误差阈值δ＞0，球B(c，(1+δ)r)视为MEB(S)的(1+δ)近似球。则核子集S_c可定义为：若真子集S_c以因子(1+δ)扩展的最小闭包球MEB包含了所有S中的样本点，即：

其中B(c，R)＝MEB(S_c)，则真子集S_c称为S的核子集。

(4)在中心约束条件下，求解最小闭包球MEB问题，即求解原n维欧氏空间Rⁿ的目标问题Φ。

完成变压器故障诊断模型构造后，将待诊断变压器记录输入到强学习机H，以得到诊断结果，具体包括步骤4，如图4：

步骤4：将强学习机H作为故障诊断模型，并将当前待故障诊断的数据样本输入到强学习机H，强学习机H利用其中各个弱学习机h_i做初步故障诊断，进而以投票的方式判定当前待测数据样本的故障类别，即得票数多的故障类别为强学习机H最终诊断结果(“弱学习机h_i投票过程”为步骤44)；

“弱学习机投票过程”阐述，因为强学习机为若干弱学习机的组合，而对于同一条待测样本各弱学习机的故障诊断结果不尽相同，因此需投票决定该待测样本的最终诊断结果。下面举例说明：

设有某条待测样本，以及“正常、低温过热故障、中温过热故障、高温过热故障、局部放电故障、低能放电故障、低能放电兼过热故障、电弧放电故障、电弧放电兼过热故障”等九种变压器运行状态，十五个弱学习机h₁，h₂，Λ，h₁₅。十五个弱学习机的故障诊断结果分别为11个“低温过热故障”、3个“低能放电兼过热故障”、1个“低能放电故障”，则该待测样本的最终诊断结果为“低温过热故障”。

若某条待测样本的故障诊断结果出现无法确定的情况，如诊断“正常”的弱学习机数和诊断“低温过热故障”的弱学习机数相同，则需对该待测样本重新进行故障诊断。经研究，Bagging算法中出现上述情况的可能性会随着弱学习机个数的增加而逐渐减小。

经过上述4个步骤的训练学习之后，形成基于集成学习Bagging算法的变压器故障诊断模型，从而实现对待故障诊断设备记录的精确诊断。

本发明在诊断变压器故障方面，相比传统单分类器方法，有较好的诊断精度，提高了变压器故障诊断算法的实用性。

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。

Claims (2)

1.基于集成学习Bagging算法的变压器故障诊断方法，其特征是所述方法包括以下步骤：

步骤1：数据采集及数据预处理，得到初始历史数据样本集                                                

；该步骤具体由步骤1.1-步骤1.2组成：

步骤1.1：数据采集：采集目标变压器各个时间监测点的气体含量及运行状态；

步骤1.2：数据预处理：设定时间监测点的气体含量及运行状态的特征属性，构成初始历史数据样本集

；

步骤2：将初始历史数据样本集构造成集成学习Bagging算法中弱学习算法可读的数据集

；该步骤具体由步骤2.1-步骤2.3组成：

步骤2.1：数据归一化：将初始历史数据样本集

中气体含量按照各自的取值范围全部归一化到[0,1]区间；

步骤2.2：数值化处理：将变压器运行的运行状态设定为数值型类别标号；

步骤2.3：将初始历史数据样本集

中的气体含量作为故障诊断模型的输入向量；将初始历史数据样本集

中的变压器的运行状态作为故障诊断模型的输出向量，故障诊断模型的输入向量和输出向量构成了数据集

；

步骤3：选定球向量机BVM作为集成学习Bagging算法中的弱学习算法并对该弱学习算法循环调用，完成数据集

的训练，从而得到弱学习机序列，序列中包含各代弱学习机

，该序列即为强学习机

；该步骤具体由步骤3.1-步骤3.5组成：

步骤3.1：对集成学习Bagging算法初始化，设定Bagging算法最大迭代次数，并设定弱学习算法的训练参数；

步骤3.2：以指定概率从数据集

中有放回地选取训练样本子集

，

，作为弱学习算法的训练样本子集

；

步骤3.3：将训练样本子集

输入到弱学习算法训练，得到对应的弱学习机；

步骤3.4：检查当前集成学习Bagging算法是否达到算法的最大迭代次数

，若已达到，则执行步骤3.5；否则，返回步骤3.2；

步骤3.5：输出弱学习机序列，即强学习机

；

步骤4：将强学习机

作为故障诊断模型，并将当前待故障诊断的数据样本输入到强学习机

，强学习机

利用其各代弱学习机

做初步故障诊断，进而以投票的方式判定当前待故障诊断的数据样本的故障类别，得票数多的故障类别即为强学习机

的最终诊断结果。

2.根据权利要求1所述基于集成学习Bagging算法的变压器故障诊断方法，其特征是所述步骤1.1中数据采集的方法为油中溶解气体分析技术DGA。

3. 根据权利要求1所述基于集成学习Bagging算法的变压器故障诊断方法，其特征是所述步骤1.1中气体包括氢气H₂、一氧化碳CO、甲烷CH₄、乙烯C₂H₄、乙烷C₂H₆、乙炔C₂H₂和二氧化碳CO₂。

4. 根据权利要求1所述基于集成学习Bagging算法的变压器故障诊断方法，其特征是所述步骤3.1中弱学习算法的训练参数包括核函数类型、核函数参数和惩罚因子。

5. 根据权利要求1所述基于集成学习Bagging算法的变压器故障诊断方法，其特征是所述步骤3.2中的指定概率为50%。

Images