诱导有序加权证据推理的变压器故障诊断方法及系统

技术领域

本发明属于变压器故障诊断领域，更具体地，涉及一种将诱导有序加权证据推理的多源信息融合理论运用于变压器FRA的故障诊断方法及系统。

背景技术

电力变压器是电力系统最重要的设备之一，其结构复杂，在实际运行中面临着各种类型的危险。其中绕组机械故障造成的损坏是最常见的电力变压器故障原因。具体来说，短路电流与电磁力的作用可能会导致机械轴向/径向位移变形。机械故障检测中使用最广泛的方法是频率响应分析(Frequence Response Analysis，FRA)，其思路是获取变压器的传递函数，并与该设备的标准值进行比对来辨识故障。然而，进行比较的过程与结果解释都要取决于专家经验的判断，目前尚无统一的故障辨别标准。专家经验与研究专长的差异可能导致分析结果的不准确，这是基于FRA方法的共同缺点。常见方法是通过研究计算出扫频响应曲线的统计指标，并设置一定范围的设备健康阈值，但该方法的适应性较差，设定好的指标只针对特定设备有效。目前的研究还尝试采用机器学习的分类算法实现扫频响应曲线的故障辨识，但该方法所需数据量较大，在部分实践场景中可能较难实现。

发明内容

针对现有技术的以上缺陷或改进需求，本发明提出了一种诱导有序加权证据推理的变压器故障诊断方法及系统，其目的是为提高电力设备故障诊断方法的智能化程度及诊断准确率，并能够适用于数据信息不完全、故障样本缺乏及样本来源可靠性多样化的情形。

为实现上述目的，按照本发明的一个方面，提供了一种诱导有序加权证据推理的变压器故障诊断方法，包括：

(1)加载变压器FRA的典型数据样本，将诊断标签设置为识别框架；

(2)加载待诊断设备的检测数据；

(3)计算待诊断设备的检测数据曲线与识别框架内所有特征数据曲线的基本可信度分配，构建信度决策矩阵；

(4)根据待诊断设备的检测数据样本来源，计算诱导有序加权平均算子及其诱导向量；

(5)计算指标权重向量；

(6)通过诱导有序加权的证据理论融合所有证据，计算综合评价信度，从而确定诊断结果。

在一些可选的实施方案中，步骤(1)中加载的变压器FRA的典型数据样本包括以下类型：实测样本：该设备的历史检测数据或是厂商检测数据，该数据样本类型是最准确的参考数据；同型号设备样本：其他同型号设备的检测数据，该数据是较准确的参考数据；精确仿真样本：针对该设备建立详细仿真模型后得到各类标签的样本；快速仿真样本：考虑到硬件、时效性的各种限制，针对该设备建立简化模型，并通过仿真获取各类标签的样本。

在一些可选的实施方案中，在步骤(1)中，将诊断标签设置为识别框架的具体方法为：对于诊断的第一阶段，首先将待检测数据划分为健康、不健康两种诊断类型，则识别框架为Θ＝{健康，不健康}，或是划分为不同的健康等级；若诊断结果为不健康状态，并且需要进一步检测故障类型，此时进入诊断的第二阶段，故障类型为F

在一些可选的实施方案中，在步骤(3)中，计算待诊断设备的检测数据曲线与识别框架内所有特征数据曲线的基本可信度分配的具体方法为：

将待诊断设备的检测数据曲线

通过曲线相似性算法计算待诊断设备的检测数据曲线的各段与识别框架内所有特征数据曲线的各段的距离P

在一些可选的实施方案中，在步骤(3)中，构建信度决策矩阵的具体方法为：

定义分布式评估向量S为：S(X

在一些可选的实施方案中，在步骤(4)中，计算诱导有序加权平均算子及其诱导向量的具体方法为：

确定待诊断设备的诱导变量u

定义有序加权平均算子F满足：

在一些可选的实施方案中，在步骤(5)中，计算指标权重向量的具体方法为：

采用平均分配，针对属性X

或者，采用客观权重分配方法计算各自的权重

在一些可选的实施方案中，步骤(6)的具体实现方法为：

将不同属性X

m

合成可信度分配：

归一化处理，得到综合诊断结果：

则对于识别框架的各个状态Θ＝{H

按照本发明的另一方面，提供了一种诱导有序加权证据推理的变压器故障诊断系统，包括：

预处理模块，用于加载变压器FRA的典型数据样本，将诊断标签设置为识别框架；

待检测数据获取模块，用于加载待诊断设备的检测数据；

信度构建模块，用于计算待诊断设备的检测数据曲线与识别框架内所有特征数据曲线的基本可信度分配，构建信度决策矩阵；

诱导有序加权计算模块，用于根据待诊断设备的检测数据样本来源，计算诱导有序加权平均算子及其诱导向量；

权重计算模块，用于计算指标权重向量；

诊断模块，用于通过诱导有序加权的证据理论融合所有证据，计算综合评价信度，从而确定诊断结果。

按照本发明的另一方面，提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述任一项所述方法的步骤。

总体而言，通过本发明所构思的以上技术方案与现有技术相比，能够取得下列有益效果：

本发明将多源信息融合理论运用到扫频响应分析的故障诊断情形中，结合了指标分析法的客观性与人工智能证据理论的自适应性这两者的优点。通过诱导有序加权(Induced Ordered Weighted，IOW)因子，考虑不同场景信息，从而实现权重的优化，灵活的改变诊断的重点，从而提高诊断结果的针对性。此外，证据推理理论可以处理不完全信息，在实际诊断过程中常存在着的数据缺失或样本类型不全问题，该方法通过引入“不确定”状态H

附图说明

图1是本发明实施例提供的一种基于诱导有序加权证据推理的变压器故障诊断方法的流程示意图；

图2是本发明实施例提供的一种曲线分段与距离计算方法示意图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。

证据理论方法是一种可以表达不确定性的推理理论，通常用于融合多源不确定信息，例如多传感器、多位专家的意见等，从而提高单一来源信息的准确性。本发明将多源信息融合理论运用到扫频响应分析的故障诊断情形，结合了指标分析法的客观性与人工智能推理理论的自适应性这两者优点，且能够适用包含不确定信息的场景。这是因为在实际诊断过程中往往存在着数据的缺失或样本类型不全的问题。此外，本发明通过引入诱导有序加权因子，可以灵活的改变诊断的着重点，提高诊断结果的针对性。

如图1所示是本发明实施例提供的一种诱导有序加权证据推理的变压器故障诊断方法的流程示意图，在图1所示的方法中包括以下步骤：

S1：加载变压器FRA的典型数据样本，将诊断标签设置为识别框架；

S2：加载待诊断设备的检测数据；

S3：计算待诊断设备的检测数据曲线与识别框架内所有特征数据曲线的基本可信度分配(Basic Probability Assignment，BPA)，构建信度决策矩阵；

S4：根据待诊断设备的检测数据样本来源(可构建指标：设备来源、时间、修正因子)，计算诱导有序加权平均算子及其诱导向量；

S5：计算指标权重向量；

S6：通过诱导有序加权的证据理论融合所有证据，计算综合评价信度，从而确定诊断结果；

S7：若还需进一步判断，则更改识别框架，返回步骤S1。

进一步地，步骤S1中加载的典型数据样本包括以下类型：a)实测样本：该设备的历史检测数据或是厂商检测数据，该数据样本类型是最准确的参考数据；b)同型号设备样本：其他同型号设备的检测数据，该数据是较准确的参考数据；c)精确仿真样本：针对该设备建立详细仿真模型后得到各类标签的样本；d)快速仿真样本：考虑到硬件、时效性等限制，针对该设备建立简化模型，并通过仿真获取各类标签的样本。

其中，当存在多个同样标签的数据时，一般选择可靠性最高的数据。

进一步地，步骤S1中将诊断标签设置为识别框架的具体方法为：对于诊断的第一阶段，首先将待检测数据划分为健康、不健康两种诊断类型，则识别框架为Θ＝{健康，不健康}，或是划分为不同的健康等级；若诊断结果为“不健康”状态，并且需要进一步检测故障类型，则需要返回步骤S1，此时即诊断的第二阶段，故障类型为F

为了统一描述，将识别框架中的所有元素记为Θ＝{H

其中，在实际运用中，可能不需要进行故障类型的判别，只需要进行故障定位，则此时的故障位置检测过程属于诊断的第二阶段。

进一步地，步骤S3中计算基本可信度分配的步骤为：

S3.1：首先将待诊断设备的检测数据曲线

其中，对于同样隶属于识别框架中相同元素的多条曲线，则选取数据来源可信度较高的曲线作为研究对象。

S3.2：通过曲线相似性算法(例如欧氏距离、动态时间规整、弗雷歇距离或多线位置距离等)，计算待诊断设备的检测数据曲线的各段与识别框架内所有特征数据曲线的各段的距离P

S3.3：然后将各曲线距离转换为信度β

其中，该距离-信度转换公式可以有多种方式，例如：①

例如，假设对于某一段曲线X

进一步地，步骤S3中信度决策矩阵的构建方法为：

定义分布式评估向量S为：S(X

其中，被评价对象A

进一步地，步骤S4中计算诱导有序加权平均算子的方法为：

S4.1：首先确定待诊断设备的诱导变量u

其中，诱导变量u

S4.2：计算有序加权平均(Ordered weighted averaging,OWA)算子，定义OWA算子F满足：

进一步地，步骤S5中指标权重向量的计算方法为：一般先平均分配，即针对属性X

进一步地，步骤S6通过证据融合计算诊断结果的方法为：

S6.1：将不同属性X

m

S6.2：合成可信度分配：

S6.3：归一化处理，得到综合诊断结果：

则对于识别框架的各个状态Θ＝{H

进一步地，步骤S7通过更改识别框架，返回步骤S1进一步诊断定位，从而实现降低复杂度的目的。通常根据实际需求，按照健康状况、故障类型、故障位置的层次依次辨识。如果不需要进行故障类型判断，则按照健康状况、故障位置进行辨识。

以下结合具体实例对本发明进行详细说明。本发明不仅适用于变压器绕组扫频响应分析的故障诊断中，还可以推广至其他设备诊断领域。

流程图如图1所示，根据步骤S1，首先加载变压器FRA的测试指纹数据

表1 加载数据概况

上表中的数据

然后根据步骤S3，计算基本可信度分配，构建信度决策矩阵。首先将待诊断曲线

其中x

将曲线距离转换为信度β

然后计算信度决策矩阵：

表2 变压器FRA故障诊断信度决策矩阵

注：“*”表示不确定量，可能由于数据缺失、损坏等造成的不确定性结果，其概率平均分配。

在本实施例中，需要重点关注绕组短路这一故障类型(H

然后根据步骤S5，计算权重因子。对于不同的故障类型，其初始权重相同：ω’＝(ω’

接下来执行步骤S6，通过诱导有序加权的证据理论融合所有证据，计算综合评价信度，得到：β

本申请还提供一种计算机可读存储介质，其上存储有计算机程序，程序被处理器执行时实现方法实施例中的诱导有序加权证据推理的变压器故障诊断方法。

需要指出，根据实施的需要，可将本申请中描述的各个步骤/部件拆分为更多步骤/部件，也可将两个或多个步骤/部件或者步骤/部件的部分操作组合成新的步骤/部件，以实现本发明的目的。

本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。