CN114372693A - 一种基于云模型和改进ds证据理论的变压器故障诊断方法 - Google Patents

Abstract

本发明公开了一种基于云模型和改进DS证据理论的变压器故障诊断方法，包括：步骤1、获取变压器油中溶解气故障样本，建立变压器的故障样本集；步骤2、采用故障样本集建立各类型故障的故障标准云模型；步骤3、计算待诊断样本与故障标准云模型间的最终隶属度矩阵；步骤4、根据最终隶属度矩阵确定不同故障下的基本概率分配；步骤5、利用Pignistic概率距离构建相似度并对基本概率分配进行修正，引入平均支持度加权优化证据融合规则，并使用改进D‑S证据理论对基本概率分配进行融合，以得到最终融合结果；步骤6、在最终融合结果中，选取基本概率最大值对应的变压器状态作为诊断结果。本发明有效地提高了变压器智能诊断方法的准确率。

Description

一种基于云模型和改进DS证据理论的变压器故障诊断方法

技术领域

本发明涉及电气设备技术领域，尤其涉及一种基于云模型和改进DS证据理论的变压器故障诊断方法。

背景技术

在油浸式变压器故障诊断中，油中溶解气分析法是目前常用的方法之一，基于油中溶解气分析法的传统诊断方法，如IEC比值法、改良三比值法、Duval三角法等，上述方法使用简单、诊断快速，经常运用于工程现场，但是存在编码不完备、诊断边界过于绝对化的问题，导致其识别率不高。

随着人工智能在电网中广泛应用，学者们将智能诊断方法引入变压器故障诊断领域。由于变压器故障与其特征之间对应关系存在模糊性的问题，所以常规智能诊断方法应用在变压器故障诊断上效果不明显。支持向量机、相关向量机和神经网络等不能表示除诊断结果以为的故障信息，因此上述方法在表达变压器各个故障的模糊性和不确定性存在不足。

发明内容

本发明所要解决的技术问题在于，提供一种基于云模型和改进DS证据理论的变压器故障诊断方法，以解决由于变压器故障与其特征之间对应关系存在模糊性，而导致变压器智能诊断方法准确率不高的问题。

为了解决上述技术问题，本发明提供了一种基于云模型和改进DS证据理论的变压器故障诊断方法，包括：

步骤1、获取变压器油中溶解气故障样本，建立变压器的故障样本集D；

步骤2、采用所述故障样本集D建立各类型故障的故障标准云模型；

步骤3、计算待诊断样本与所述故障标准云模型间的最终隶属度矩阵；

步骤4、根据所述最终隶属度矩阵确定不同故障下的基本概率分配；

步骤5、利用Pignistic概率距离构建相似度Sim并对所述基本概率分配进行修正，引入平均支持度加权优化证据融合规则，并使用改进D-S证据理论对所述基本概率分配进行融合，以得到最终融合结果；

步骤6、在所述最终融合结果中，选取基本概率最大值对应的变压器状态作为诊断结果。

在一些实施方式中，所述故障样本集D包含有多种变压器故障类型的油色谱样本，所述油色谱样本由多种特征气体的含量特征构成。

在一些实施方式中，所述步骤2包括：

步骤21、针对所述故障样本集D，对原始数据进行归一化处理；

步骤22、通过分析所述故障样本集D计算变压器各状态下，各故障特征值的期望E_x、标准差E_n和离散程度H_e，进而构建所述故障标准云模型。

在一些实施方式中，所述步骤3包括：

步骤31、根据故障特征值的期望E_x、标准差E_n和离散程度H_e，计算所述待诊断样本的第i个故障特征下第j类故障的隶属度μ_ij；

步骤32、所述待诊断样本的各特征与故障标准云模型间的隶属度值构建隶属度矩阵μ_ij＝[μ_i1,μ_i2,…,μ_in]，其中n为变压器故障类型总数；

步骤33、计算第i个故障特征对应的不确定性隶属度φ_i＝1-max(μ_ij)；计算所述待诊断样本的所有特征下各类故障的隶属度和对应的不确定隶属度，构建最终隶属度矩阵μ_m×(n+1)。

在一些实施方式中，所述步骤4包括：

步骤41、将所述最终隶属度矩阵每一行进行归一化处理，得到所述待诊断样本的所述基本概率分配T_m×(n+1)。

在一些实施方式中，所述步骤5包括：

步骤51、基于Pignistic概率函数构建相似度Sim，并将所述基本概率分配T_m×(n+1)与对应的所述相似度Sim相乘，以对所述基本概率分配进行修正；

步骤52、基于平均支持度加权方法对D-S证据理论Dempster融合规则进行优化；

步骤53、使用优化后的融合规则对修正后的所述基本概率分配进行信息融合，得到所述最终融合结果。

在一些实施方式中，所述步骤6包括：

步骤61、采用最大值规则来确定最终融合结果，即是选择最大基本概率分配函数值作为决策结果。

在一些实施方式中，所述多种特征气体为氢气(H₂)、甲烷(CH₄)、乙烷(C₂H₆)、乙烯(C₂H₄)和乙炔(C₂H₂)，所述多种变压器故障类型包括正状态(N)、低温过热(T1)、中温过热(T2)、高温过热(T3)、局部放电(PD)、低能放电(D1)和高能放电(D2)。

在一些实施方式中，所述步骤21中采用百分比归一化方法，计算公式如下所示：

式中，X_i表示第i个特征气体占所有特征气体的百分比；x_i表示第i个特征气体含量；THC表示总烃气体含量；H₂表示氢气气体含量。

在一些实施方式中，所述步骤22中各故障特征值的期望E_x、标准差E_n和离散程度H_e，具体计算方法如下：

设x_ijk为某一变压器故障状态下的油中溶解气数据，表示在第i个特征的第j类故障，有k个样本，则第i个特征的第j类故障下期望E_x，计算公式如下所示：

则第i个特征的第j类故障下标准差E_n，计算公式如下所示：

离散程度H_e取值0.01。

在一些实施方式中，所述步骤31中隶属度μ_ij计算方法，如下所示：

式中，E_nij′是E_n和H_e产生的一个正态随机数。

在一些实施方式中，所述步骤51中相似度Sim计算方法如下：

取多种变压器故障类型为辨识框架U，取m表示识别框架U中的基本概率分配函数，Pignistic概率函数BetPm(ω)是从U到[0,1]的映射关联函数，计算公式如下所示：

式中，|A|表示子集A的基数；

基于Pignistic概率函数BetPm(ω)计算事件A的两证据体m₁和m₂之间的差值矩阵Dif_Bet，计算公式如下所示：

基本概率分配函数值在[0,1]上，且Dif_Bet可表示证据体m₁和m₂之间差异性，所以定义两证据体m₁和m₂之间相似度Sim，计算公式如下所示：

通过基本概率分配T_m×(n+1)与相应相似度Sim相乘，完成对所述基本概率分配的修正。

在一些实施方式中，所述步骤52中优化证据理论融合规则，计算方法如下：

采用证据之间的平均支持度加权值进行融合规则优化，融合规则如下所示：

式中，f(A)为平均支持度加权值，计算公式所示：

f(A)＝K×q(A)

式中，K为Dempster规则中的证据冲突，计算公式如下所示：

本发明与现有技术相比具有明显的优点和有益效果。借由上述技术方案，本发明一种基于云模型和改进DS证据理论的变压器故障诊断方法可达到相当的技术进步性及实用性，并具有产业上的广泛利用价值，其至少具有下列优点：

本发明的一种基于云模型和改进D-S证据理论的变压器故障诊断方法，利用云模型表示变压器故障类型与其对应故障特征之间的关系，再使用改进D-S证据理论得出决策结果。云模型能有效解决存在模糊性和随机性的问题；证据理论在处理不确定信息时具有明显优势，两者结合在解决变压器故障类型与其对应故障特征之间的模糊性问题具有明显效果。同时，该方法相比支持向量机和神经网络等方法不需要预设初始参数，解决了无法寻找最优初始参数的问题。

上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，而可依照说明书的内容予以实施，并且为了让本发明的上述和其他目的、特征和优点能够更明显易懂，以下特举较佳实施例，并配合附图，详细说明如下。

附图说明

图1示出了本发明一实施例的基于云模型和改进DS证据理论的变压器故障诊断方法的流程示意图。

具体实施方式

为更进一步阐述本发明为达成预定发明目的所采取的技术手段及功效,以下结合附图及较佳实施例，对依据本发明提出的一种基于云模型和改进DS证据理论的变压器故障诊断方法的具体实施方式及其功效，详细说明如后。

本发明实施例提供了一种基于云模型和改进DS证据理论的变压器故障诊断方法，请参阅图1所示：本发明为一种基于云模型和改进D-S证据理论的变压器故障诊断方法，诊断方法包括以下步骤：

步骤1、收集并整理变压器油中溶解气故障样本，建立变压器的故障样本集D；

步骤2、利用故障样本集D建立各类型故障的标准云模型；

步骤3、计算待诊断样本与故障标准云模型间的最终隶属度矩阵；

步骤4、根据最终隶属度矩阵确定不同故障下的基本概率分配；

步骤5、利用Pignistic概率距离构建相似度Sim并对基本概率分配进行修正，再引入平均支持度加权优化证据融合规则，并通过改进D-S证据理论对基本概率分配进行融合，以得到最终融合结果；

步骤6、在最终融合结果中，选取基本概率最大值对应的变压器状态作为诊断结果。

在一实施例中，步骤1中故障样本集D包括多种变压器故障类型的油色谱样本，油色谱样本由多种特征气体含量组成。

优选地，变压器故障类型为7种，特征气体为5种。其中，5种特征气体为氢气(H₂)、甲烷(CH₄)、乙烷(C₂H₆)、乙烯(C₂H₄)和乙炔(C₂H₂)，7种变压器故障类型包括正状态(N)、低温过热(T1)、中温过热(T2)、高温过热(T3)、局部放电(PD)、低能放电(D1)和高能放电(D2)。

在一实施例中，步骤2包括以下步骤：

步骤21、针对所建立的故障样本集D，对原始数据进行归一化处理，其目的是消除样本之间量纲不一致的影响。

进一步地，步骤21中采用的百分比归一方法，计算公式如下所示：

式中，X_i表示第i个特征气体占所有特征气体的百分比；x_i表示第i个特征气体含量；THC表示总烃气体含量；H₂表示氢气气体含量。

步骤22、通过分析故障样本集D计算变压器各状态下，各故障特征值的期望E_x、标准差E_n和离散程度H_e，进而构建起故障标准云模型。

进一步地，步骤22中各故障特征值的期望E_x、标准差E_n和离散程度H_e，具体计算方法如下：

设x_ijk为某一变压器故障状态下的油中溶解气数据，表示在第i个特征的第j类故障，有k个样本。则第i个特征的第j类故障下期望E_x，计算公式如下所示：

则第i个特征的第j类故障下标准差E_n，计算公式如下：

根据相关学术期刊离散程度H_e通常取值0.01。

在一实施例中，步骤3包括以下步骤：

步骤31、根据故障特征值的期望E_x、标准差E_n和离散程度H_e，计算待诊断样本的第i个故障特征下第j类故障的隶属度μ_ij。

进一步地，其中，步骤31中隶属度μ_ij计算方法如下所示：

式中，E_nij′是E_n和H_e产生的一个正态随机数。

步骤32、待诊断样本的各特征与故障标准云间的隶属度值构建隶属度矩阵μ_ij＝[μ_i1,μ_i2,…,μ_in]，其中n为变压器故障类型总数；

步骤33、计算第i个故障特征对应的不确定性隶属度φ_i＝1-max(μ_ij)，其中不确定性隶属度表示诊断过程中的数据误差、外界干扰等原因，引起的不确定性。计算待诊断样本所有特征下各类故障的隶属度和对应的不确定隶属度后，最终隶属度矩阵μ_m×(n+1)，如下所示：

式中，m为待诊断样本的特征总数。最终隶属度矩阵μ_m×(n+1)中每行分别表示变压器的各个故障特征，每列表示各个故障状态。

在一实施例中，步骤4包括以下步骤：

步骤41、将隶属度矩阵每一行进行归一化处理，得到待诊断样本的基本概率分配T_m×(n+1)，如下所示：

在一实施例中，步骤5包括以下步骤：

步骤51、基于Pignistic概率函数构建相似度Sim，并使用基本概率分配T_m×(n+1)与对应相似度Sim相乘，完成对基本概率分配修正。

进一步地，步骤51中相似度Sim计算方法如下：

取7种变压器故障类型为辨识框架U，取m表示识别框架U中的基本概率分配函数，Pignistic概率函数BetPm(ω)是从U到[0,1]的映射关联函数，计算公式如下所示：

式中，|A|表示子集A的基数。

基于Pignistic概率函数BetPm(ω)计算事件A的两证据体m₁和m₂之间的差值矩阵Dif_Bet，计算公式如下所示：

因为基本概率分配函数值在[0,1]上，且Dif_Bet可表示证据体m₁和m₂之间差异性，所以定义两证据体m₁和m₂之间相似度Sim，计算公式如下所示：

通过基本概率分配T_m×(n+1)与相应相似度Sim相乘，完成对基本概率分配修正。

步骤52、基于平均支持度加权方法对D-S证据理论Dempster融合规则进行优化。

进一步地，步骤52中优化证据理论融合规则，计算方法如下：

由于经典证据理论的融合规则是Dempster规则，在高冲突的情况下会导致出现有悖常理的决策结果，因此为得到合理的决策结果，本方法采用证据之间的平均支持度加权值进行融合规则优化，融合规则如下所示：

式中，f(A)为平均支持度加权值，计算公式所示：

f(A)＝K×q(A)

式中K为Dempster规则中的证据冲突，计算公式如下所示：

在一实施例中，步骤6包括以下步骤：

步骤61、采用最大值规则来确定最终融合结果，即是选择最大基本概率分配函数值作为决策结果。

通过上述步骤，即可完成对变压器故障样本的诊断，为检修人员提供可靠的诊断信息，减少变压器检修的盲目性。

在一具体实施例中，以变压器低温过热故障(T1)的油中溶解气样本进行说明。

步骤10、收集并整理变压器故障样本集D共1400条，其中7种变压器故障类型每类200条样本。实际待测样本X＝[141.19，199.39，40.38，21.96，0]，其中X的数值依次代表H₂、CH₄、C₂H₆、C₂H₄和C₂H₂的故障气体含量，单位是ppm；已知其故障类型为T1。

步骤20、基于已建立的变压器故障样本集合D，随机选取1260条样本(每类样本180条)作为训练集，实际待测样本X为测试集。根据训练集建立各类型故障的标准云参数，即计算各故障特征值的期望E_x、标准差E_n和离散程度H_e，其中离散程度H_e＝0.01，不同变压器状态下的标准云参数，如表1所示。

表1 不同变压器状态下的标准云参数

步骤30、根据待诊断样本S计算所有特征下各类故障的隶属度(包括不确定度隶属度)，如表2所示：

表2 各个故障特征的隶属度矩阵

步骤40、将隶属度矩阵(即最终隶属度矩阵)每一行进行归一化处理，得到待诊断样本的基本概率分配，如表3所示：

表3 各个故障特征的基本概率分配

步骤50、基于Pignistic概率函数构建相似度Sim并对所建立的基本概率分配进行修正。其中，各个故障特征的相似度Sim，如表4所示；基本概率分配T_m×(n+1)与对应相似度Sim相乘后，得到修正后各个故障特征的基本概率分配，如表5所示。

表4 各个故障特征的相似度Sim

表5 修正后各个故障特征的基本概率分配

步骤60、最终融合结果如表6所示，选取基本概率最大值0.219对应的变压器状态T1作为诊断结果，与该待诊断样本固有标签一致。

表6 样本S的故障诊断融合结果

本发明的一种基于云模型和改进D-S证据理论的变压器故障诊断方法，利用云模型表示变压器故障类型与其对应故障特征之间的关系，再使用改进D-S证据理论得出决策结果。云模型能有效解决存在模糊性和随机性的问题；证据理论在处理不确定信息时具有明显优势，两者结合在解决变压器故障类型与其对应故障特征之间的模糊性问题具有明显效果。同时，该方法相比支持向量机和神经网络等方法不需要预设初始参数，解决了无法寻找最优初始参数的问题。

以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明,任何熟悉本专业的技术人员，在不脱离本发明技术方案范围内，当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例，但凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属于本发明技术方案的范围内。

Claims (12)

1.一种基于云模型和改进DS证据理论的变压器故障诊断方法，其特征在于，包括：

步骤1、获取变压器油中溶解气故障样本，建立变压器的故障样本集；

步骤2、采用所述故障样本集建立各类型故障的故障标准云模型；

步骤3、计算待诊断样本与所述故障标准云模型间的最终隶属度矩阵；

步骤4、根据所述最终隶属度矩阵确定不同故障下的基本概率分配；

步骤5、利用Pignistic概率距离构建相似度并对所述基本概率分配进行修正，引入平均支持度加权优化证据融合规则，并使用改进D-S证据理论对所述基本概率分配进行融合，以得到最终融合结果；

步骤6、在所述最终融合结果中，选取基本概率最大值对应的变压器状态作为诊断结果。

2.根据权利要求1所述的基于云模型和改进DS证据理论的变压器故障诊断方法，其特征在于，所述故障样本集包含有多种变压器故障类型的油色谱样本，所述油色谱样本由多种特征气体的含量特征构成。

3.根据权利要求1所述的基于云模型和改进DS证据理论的变压器故障诊断方法，其特征在于，所述步骤2包括：

步骤21、针对所述故障样本集，对原始数据进行归一化处理；

步骤22、通过分析所述故障样本集计算变压器各状态下，各故障特征值的期望、标准差和离散程度，进而构建所述故障标准云模型。

4.根据权利要求3所述的基于云模型和改进DS证据理论的变压器故障诊断方法，其特征在于，所述步骤3包括：

步骤31、根据故障特征值的期望、标准差和离散程度，计算所述待诊断样本的第i个故障特征下第j类故障的隶属度；

步骤32、所述待诊断样本的各特征与故障标准云模型间的隶属度值构建隶属度矩阵μ_ij＝[μ_i1,μ_i2,…,μ_in]，其中n为变压器故障类型总数；

步骤33、计算第i个故障特征对应的不确定性隶属度φ_i＝1-max(μ_ij)；计算所述待诊断样本的所有特征下各类故障的隶属度和对应的不确定隶属度，构建最终隶属度矩阵。

5.根据权利要求4所述的基于云模型和改进DS证据理论的变压器故障诊断方法，其特征在于，所述步骤4包括：

步骤41、将所述最终隶属度矩阵每一行进行归一化处理，得到所述待诊断样本的所述基本概率分配。

6.根据权利要求5所述的基于云模型和改进DS证据理论的变压器故障诊断方法，其特征在于，所述步骤5包括：

步骤51、基于Pignistic概率函数构建相似度，并将所述基本概率分配与对应的所述相似度相乘，以对所述基本概率分配进行修正；

步骤52、基于平均支持度加权方法对D-S证据理论Dempster融合规则进行优化；

步骤53、使用优化后的融合规则对修正后的所述基本概率分配进行信息融合，得到所述最终融合结果。

7.根据权利要求2-6中任一项所述的基于云模型和改进DS证据理论的变压器故障诊断方法，其特征在于，所述多种特征气体为氢气、甲烷、乙烷、乙烯和乙炔，所述多种变压器故障类型包括正状态、低温过热、中温过热、高温过热、局部放电、低能放电和高能放电。

8.根据权利要求3所述的基于云模型和改进DS证据理论的变压器故障诊断方法，其特征在于，所述步骤21中采用百分比归一化方法，计算公式如下所示：

式中，X_i表示第i个特征气体占所有特征气体的百分比；x_i表示第i个特征气体含量；THC表示总烃气体含量；H₂表示氢气气体含量。

9.根据权利要求3所述的基于云模型和改进DS证据理论的变压器故障诊断方法，其特征在于，

所述步骤22中各故障特征值的期望、标准差和离散程度，具体计算方法如下：

设x_ijk为某一变压器故障状态下的油中溶解气数据，表示在第i个特征的第j类故障，有k个样本，则第i个特征的第j类故障下的期望，计算公式如下所示：

则第i个特征的第j类故障下的标准差，计算公式如下所示：

离散程度取值0.01。

10.根据权利要求4所述的基于云模型和改进DS证据理论的变压器故障诊断方法，其特征在于，

所述步骤31中隶属度计算方法，如下所示：

式中，E_nij′是标准差和离散程度产生的一个正态随机数。

11.根据权利要求6所述的基于云模型和改进DS证据理论的变压器故障诊断方法，其特征在于，所述步骤51中相似度Sim计算方法如下：

取多种变压器故障类型为辨识框架，取m表示识别框架中的基本概率分配函数，Pignistic概率函数BetPm(ω)是从辨识框架到[0,1]的映射关联函数，计算公式如下所示：

式中，|A|表示子集的基数；

基于Pignistic概率函数BetPm(ω)计算事件的两证据体m₁和m₂之间的差值矩阵Dif_Bet，计算公式如下所示：

基本概率分配函数值在[0,1]上，且Dif_Bet可表示证据体m₁和m₂之间差异性，所以定义两证据体m₁和m₂之间相似度，计算公式如下所示：

通过基本概率分配与相应相似度相乘，完成对所述基本概率分配的修正。

12.根据权利要求6所述的基于云模型和改进DS证据理论的变压器故障诊断方法，其特征在于，所述步骤52中优化证据理论融合规则，计算方法如下：

采用证据之间的平均支持度加权值进行融合规则优化，融合规则如下所示：

式中，f(A)为平均支持度加权值，计算公式所示：

f(A)＝K×q(A)

式中，K为Dempster规则中的证据冲突，计算公式如下所示：

Images