基于迁移学习和多源信息融合的电机故障诊断分析方法与流程

1.本发明涉及电机故障诊断技术领域，尤其涉及基于迁移学习和多源信息融合的电机三相电流信号诊断分析方法。


背景技术：

2.异步电机作为工农业生产活动的重要动能之一，具有结构简单易拆卸、制造材料价格低廉、可靠性强、持久耐用等一系列优点，广泛应用于采矿设备与冶金业、机床、电厂、风机等工业领域中。电机的稳定可靠运行不仅与单个机械设备的稳定、长时间运行相关，而且影响着整个生产系统的稳定性，异步电机一旦发生故障，轻则设备运行中断、工作效率降低、能源消耗增加，重则整个系统运行崩溃，产生巨大的经济损失，造成人员伤亡。
3.现阶段异步电机安全管理主要依靠定期维修，通过对电机各个部位的检查修护来判断电机的使用情况，容易产生维护盲目、资源浪费等情况。随着泛在电力物联网的提出与推进，通过先进的传感技术可监测异步电机各个部位实时运行状态，因此基于采集信号，实现电机在线早期故障诊断成为了重要研究方向。传统电机故障诊断往往依靠阈值判断的方法，但是异步电机通常在复杂的环境中运行，此时故障早期特征往往会淹没在噪声中，以至于故障发现太晚，造成电机故障修复不及时。现阶段工业设备整体趋向于无人化、智能化，人工智能算法飞速发展，故障诊断技术相应随之进行更新，无需专家经验的智能诊断算法具有诊断效果好、鲁棒性强等优点，是电机故障诊断领域的应用热点。
4.纵观电机故障诊断技术的发展，大致可分为四个类别：基于信号、基于机械理论、基于模型与基于仿真。发展最初，基于机械理论的故障诊断算法尤其流行，主要在理论上探究电机各个故障的表现形式，但是无法建立准确描述故障机理的数学模型，且模型无法体现电机实际运行状态，存在明显缺陷。由于信号检测技术、信号处理与人工智能算法的发展，基于信号与基于传统机器学习模型的故障诊断算法相应推进，在故障诊断领域中具有诊断效率高、速度快、精度好等显著优点，目前的重难点在于需人工选取特征、需大量数据训练网络等。基于仿真的故障诊断算法主要是采用流行的数字仿真软件构建近似于实际电机的有效物理模型，输出电流、振动甚至声信号等响应，直观反映出电机运行时的物理变化。随着大数据、智能化时代的来临，基于深度学习算法的智能化故障诊断技术渐渐取代基于传统智能算法的诊断技术，深度学习作为新时代研究重点，在电机故障领域得到了广泛的探索与尝试。
5.每种诊断算法都有其优越性与局限性，现有智能诊断算法在跨域故障诊断中准确率不高，且大多数算法需要大量训练数据。实际工程中，由于设备运行过程中基本处于健康状态，故障发生一般在极短时间内，导致故障数据不足、样本不均衡等问题。因此，本发明针对目标任务中样本较少的问题，研究基于迁移学习和多源信息融合的电机故障诊断分析方法，提高诊断算法在设备不同运行状态下与训练样本不足时的诊断准确率。


技术实现要素：

6.为解决上述技术问题，本发明提供一种基于迁移学习和多源信息融合的电机故障诊断分析方法，针对直接采用基于resnet的深度学习诊断模型在变工况下诊断准确率较低的问题，引入深度迁移学习理论，构建了基于特征迁移的变工况电机故障诊断模型；采用jmmd距离度量方式进行特征空间分布适配，用于目标域测试数据的故障分类；研究多源信息决策层融合方法，构建了基于多源信息融合的电机故障诊断模型；将原始电流信号与经eemd分解后的振动信号采用特征迁移模型进行分类预测，最终采用lance改进d-s证据融合理论对结果进行融合并输出最终的结果。
7.为实现上述目的，本发明是以如下技术方案实现的：一种基于迁移学习和多源信息融合的电机故障诊断分析方法，包括以下步骤：
8.s1.信号处理：对输入三相电流信号和两轴振动信号进行eemd分解，构建多维数据样本集；
9.s2.模型训练：基于特征迁移模型提取域不变深度特征并训练模型；
10.s3.决策融合：将模型的softmax输出的概率作为证据体的mass函数值，通过可信度计算修正mass函数值，进行证据融合；
11.s4.故障分类：取概率最大作为模型预测结果，输出模型诊断准确率。
12.优选的，对输入三相电流信号和两轴振动信号进行eemd分解，构建多维数据样本集的步骤包括：
13.(1)将原始信号分为源域数据、目标域训练集与目标域测试集；
14.(2)三相电流信号按电流输入策略1进行处理，转换为矩阵构建样本集；
15.(3)振动信号采用eemd分解后，选取其中4个有效imf分量，并转换为矩阵构建样本集。
16.优选的，所述划分原始信号的原则为，进行同工况实验时，源域数据与目标域数据来自同一转速；变工况实验时，源域数据与目标域数据来自来自不同转速；采用多通道深度学习网络对三相电流信号进行处理，将三相电流作为一个整体，送入3通道深度网络中提取特征，在卷积之后将特征叠加。对于同工况实验，由于源域数据与目标域数据来自于同一转速，特征分布相近，计算出jmmd损失接近于0。
17.优选的，步骤(3)具体过程如下：
18.31.给定输入信号x(t)，对其叠加一组高斯白噪声，得到信号xi(t)；
19.xi(t)＝x(t)+ni(t)
20.式中，xi(t)为叠加白噪声的第i次分解的信号，ni(t)为第i次叠加的白噪声(i＝1,2,3,
…
,m)；
21.32.对xi(t)进行emd分解，得到各阶imf分量；
[0022][0023]
式中，c
ij
(t)为第i次分解得到的第j个imf分量，r
ij
(t)为第i次分解的剩余分量。
[0024]
33.将m次分解得到对应imf分量进行求和以抵消白噪声，得到最终的各个imf分量；
[0025][0026]
34.计算各阶imf分量cj(t)与原信号x(t)的相关系数；
[0027][0028]
35.选取相关系数大于阈值的imf分量作为最终的振动信号输入，过滤高阶冗余分量。
[0029]
优选的，特征迁移模型采用联合最大平均偏差jmmd计算多层网络的联合分布距离，解决由网络结构引起的域位移问题，jmmd损失函数计算如下：
[0030][0031]
式中，h
l
表示第l层rkhs(再生希尔伯特空间)；|l|为对应集合层数；为张量积在rkhs中的特征映射；z
sl
与z
tl
分别为第l层中源域、目标域的激活。
[0032]
优选的，对于独立同分布的单一事件集合ω＝{a1，a2，...，a
t
}，其中m为概率分配函数，则ω事件满足如下约束：
[0033][0034]
这里，m(ai)表示事件ai的准确信任度。
[0035]
对于两个证据m1与m2，可以基于兰氏lance距离度量进行信任度计算，采用融合规则得到最终的融合结果，证据间的兰氏距离为：
[0036][0037]
对于多个证据，可计算兰氏距离矩阵：
[0038][0039]
式中，n为证据个数，矩阵为主对角线为0的对称矩阵。
[0040]
优选的，采用度量证据间的距离表示其相似性，并利用可信度指标对信任度函数进行优化，定义证据m1与m2的相似性计算式如下：
[0041]s12
＝s(m1，m2)＝1-d
12
＝1-d(m1，m2)，0≤s
12
≤1
[0042]
对于n个证据的问题，可以得到相似度矩阵为：
[0043][0044]
式中，s
ij
表示证据mi与mj之间的相似性，值越大，表示证据间具有一致性程度越高，这里定义证据mi的可信度指标reli如下：
[0045][0046]
利用信任度指标reli对信任度函数进一步修正，得到m
avg
(ai)，修正过程如下式：
[0047][0048]
与现有技术相比，本发明具有如下优点和技术效果：
[0049]
在样本预处理阶段，本发明对原始电流信号直接进行数据转换，对x、y方向振动信号经eemd分解后的有效imf分量进行数据转换，分别构建电流与振动样本集。在深度特征迁移模型训练阶段，本发明将训练好的网络用于测试集的预测结果输出，将softmax层的概率函数作为预测结果输出，得到电流信号与振动信号的3个证据体。本发明分别对每一个样本的概率输出依据d-s证据融合理论进行结果融合，并采用兰氏距离计算证据的可信度用于修正mass函数，最终得到的融合结果即为采用本发明得到的诊断准确率。
[0050]
在同工况下进行故障分类识别时，经过证据融合后的诊断结果相对于单一信号有所提升，多工况的平均诊断准确率为98％，误差率减少至2％；在变工况下进行多源信息融合故障分类识别时，相对于电流信号采用参数迁移诊断时诊断准确率有所提升，各工况下诊断准确率相较于传统方法提高3％以上。
附图说明
[0051]
图1为本发明流程示意图；
[0052]
图2为本发明实施例的振动信号采集示意图；
[0053]
图3为本发明实施例的振动信号eemd分解图；
[0054]
图4为本发明实施例的基于特征迁移模型的电机故障诊断框图；
[0055]
图5为本发明实施例的基于特征迁移模型的电机故障诊断模型运行结果；
[0056]
图6为本发明实施例的特征迁移模型a
→
b工况下训练集与测试集诊断结果的混淆矩阵；
[0057]
图7为本发明实施例的同工况下异类信息融合前后的诊断准确率对比；
[0058]
图8为本发明实施例的变工况下异类信息融合前后的诊断准确率对比。
具体实施方式
[0059]
需要说明的是，在不冲突的情况下，本技术中的实施例及实施例中的特征可以相
互组合。下面将参考附图并结合实施例来详细说明本技术。
[0060]
需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
[0061]
实施例
[0062]
如图1-3所示，本实施例提供了基于迁移学习和多源信息融合的电机故障诊断分析方法，包括以下步骤：
[0063]
s1.基于不同工况下的电机x轴和y轴方向振动信号，进行eemd分解，并经过相关性筛选得到原始振动的前四阶imf分量；eemd分解及筛选的具体步骤为：
[0064]
(1)给定输入信号x(t)，对其叠加一组高斯白噪声，得到信号xi(t)；
[0065]
xi(t)＝x(t)+ni(t)
[0066]
式中，xi(t)为叠加白噪声的第i次分解的信号，ni(t)为第i次叠加的白噪声(i＝1,2,3,
…
,m)；
[0067]
(2)对xi(t)进行emd分解，得到各阶imf分量；
[0068][0069]
式中，c
ij
(t)为第i次分解得到的第j个imf分量，r
ij
(t)为第i次分解的剩余分量。
[0070]
(3)将m次分解得到对应imf分量进行求和以抵消白噪声，得到最终的各个imf分量；
[0071][0072]
(4)计算各阶imf分量cj(t)与原信号x(t)的相关系数；
[0073][0074]
(5)选取相关系数大于阈值的imf分量作为最终的振动信号输入，过滤高阶冗余分量。
[0075]
s2.将相关系数较大的前四个imf分量送入4通道输入深度学习网络，三相电流直接送入3通道深度学习网络；
[0076]
目标工况的三相电流数据作为目标域数据，分为带标签的微调样本与无标签的目标域测试样本，其中源域数据与目标域测试数据均采用对应转速下三相电流的100个样本，微调样本采用带标签目标域数据30个，每个样本均构造为3*64*128大小矩阵输入网络中。
[0077]
s3.输入信号经resnet网络后提取出深度特征，最小化jmmd损失和resnet模型交叉熵损失；初始化resnet深度网络的参数与权值，对于源域数据与目标域训练数据同时输入到参数共享的源resnet模型与目标resnet模型中，通过计算特征空间的jmmd损失与源resnet模型的交叉熵损失，最小化网络整体损失函数更新网络参数；同时输出目标域测试数据的分类结果。其中，网络损失函数整体朝着降低的方向进行，而在迭代过程中，通过更
新网络参数，训练集验证结果与目标域测试数据结果应越来越高，直至诊断准确率达到顶峰并稳定下来。
[0078]
s4.将softmax输出的概率作为每个证据体的mass函数值，并通过可信度计算修正mass函数值，采用d-s进行证据融合；
[0079]
原始三相电流数据送入深度学习网络，在第一个卷积层后对三项电流进行特征融合，首先对特征迁移模型进行训练，并用于测试集输出电流信号诊断结果。对于振动信号，将x和y轴方向振动信号处理后分为源域数据与目标域数据，在第一个卷积层后将4个imf分量的特征进行融合，并采用训练集训练特征迁移模型，用于预测测试集的分类结果。采用lance改进d-s理论进行预测结果融合，融合计算过程如表1所示。
[0080]
表1
[0081][0082][0083]
s5.通过证据融合后的结果，将概率最大所对应的样本标签作为模型的预测结果，并输出模型的诊断准确率。
[0084]
更进一步地，在本实施例中，还利用本发明所提出的算法模型进行电机试验台进行模拟故障分析。
[0085]
(1)mfs-rds电机故障实验台
[0086]
使用mfs-rds作为电机故障仿真平台，采集正常、定子匝间短路、转子断条、气隙偏心、不平衡与轴承故障电机的稳定运行下的定子三相电流，如图2所示。以8192个连续采样点为一个电流信号样本，采集1200、1500、1800rpm下的每种状态下的定子三相电流样本各200个，随机选择100个样本构建训练集，剩余100个作为测试集，其中，用于样本类别说明如表3和表2所示。
[0087]
表2
[0088]
电机状态训练集测试集类别正常1001000转子断条故障1001001气隙偏心故障1001002定子匝间短路故障1001003
供电不平衡故障1001004轴承故障1001005
[0089]
(2)基于特征迁移的变工况故障诊断
[0090]
为研究基于特征迁移的电机故障诊断框架，引入jmmd，实现异步电机的变工况故障诊断研究，算法结构如图4。设置a
→
b、a
→
c、b
→
c三种工况用于研究特征迁移诊断模型在跨转速诊断下的可用性。并引入两个对比试验：
①
引入mmd、mk-mmd，对比不同的距离度量方法对于诊断准确率的影响；
②
设置目标域训练集分别为5、10、20、30个，分析目标域训练样本个数对于故障诊断准确率的影响。
[0091]
采用特征迁移模型分别在a
→
b、b
→
c、a
→
c三种工况下对电机故障进行诊断，输出诊断准确率准确率，同时计算诊断模型的micro-f1值、macro-f1值如表3所示：
[0092]
表3
[0093]a→
ba
→
cb
→
c准确率/％92.17％88％95.17％micro-f1值0.92170.880.9517macro-f1值0.92230.88180.9518
[0094]
由表可知，网络在b
→
c工况下达到了95.17％，而在其他工况下故障诊断正确率有所降低，a
→
c诊断准确率最低，与参数模型结论相同，故障诊断准确率受转速跨度影响。由于所采取的类别较为均衡，两种f1指标相差不大。
[0095]
以a
→
b为例，将1500rpm下的源域数据与1800rpm下的目标域练训数据分别送入参数共享的源模型与目标模型中进行训练，并采用1800rpm下的目标域测试数据进行分类预测，得到诊断结果如图5所示。由图中可以看出，jmmd损失函数在50次迭代之前迅速降低并达到稳定状态，这个过程主要为寻找特征映射空间以提取各个状态下的域不变特征；在此迭代之后，网络整体损失函数开始降低，这个过程表示更新网络参数以提高源域数据的诊断结果，同时目标域测试结果也随之上升。由诊断准确率曲线可以看出，准确率随着迭代次数增加缓慢上升，训练到951次时，源域训练准确率达到了97.83％，而目标域测试集的分类效果达到了91.83％，经过2000次的迭代过程，最终训练集与测试集的诊断准确率为97.17％、92.17％。
[0096]
训练集与测试集运行的混淆矩阵如图6所示，可以看出采用特征迁移模型对a
→
b进行跨域诊断，可有效区分电机正常状态、断条故障与匝间短路故障状态，对于其他状态，在训练集的状态区分结果要好于测试集，这是由于网络在训练阶段计算了源域的交叉熵损失。由诊断结果可看出，特征迁移后在偏心故障、绕组不平衡故障与轴承故障均出现了误诊，在总数100个样本中其误诊样本个数均在15个左右，偏心故障误诊个数稍多，达到了18个。
[0097]
(3)基于多源信息融合的电机故障诊断
[0098]
基于lance改进d-s证据理论，本发明对于电流信号与x、y两个方向的振动数据，分别利用本发明所提出的特征迁移模型提取域不变深度特征并输出分类结果，作为证据体m1、m2与m3，并采用信息融合理论输出最终结果。为验证本发明所选取多源信息融合诊断算法的优越性，分别对于所采集到的1500rpm、1800rpm、2100rpm下的电流与振动信号，设置同工况与变工况的实验对比，并对模型诊断结果进行分析。
[0099]
通过计算兰氏距离矩阵、相似度矩阵、证据的可行度对原始mass函数值进行修正，其中对于每个证据的可信度分别为0.3168、0.3452、0.3381，表示对于示例样本，融合过程中更加信任m2，其次为m3，但总体可信度相差不大。经过第一次融合，结果得到此样本属于标签4的概率为0.3525；第二次融合后，概率上升为0.4163，即经过3个证据的mass值融合后，可以判定为此样本属于绕组不平衡故障。
[0100]
对于示例样本的证据融合可知，即使两个证据源给出了错误判定且有效证据源的正确判定确信度不高，但是经过d-s证据融合，诊断结果依然被判定为正确结论，同时提升了样本所属类别的确信度。通过证据融合，可对三个信号源的诊断输出概率进行融合，并输出最终的模型分类结果，分别对比不同信号的输出结果与融合后的结果，对比如表4所示。表中为各状态下诊断正确的样本个数。
[0101]
表4
[0102]
电机状态类型三相电流x振动信号y振动信号融合后结果01009587971100100100100210010010010031009499994728499965999998100平均准确率95.17％95.33％97.17％98.67％
[0103]
将所研究的信息融合算法用于电机变工况故障诊断研究，考虑变工况下深度特征迁移模型诊断效果要好于深度参数迁移模型，本发明分别将电流信号与处理后的振动信号输入至深度特征迁移模型进行状态识别，设置三种工况，分别输出工况1：1500rpm
→
1800rpm、工况2：1500rpm
→
2100rpm、工况3：1800rpm
→
2100rpm下的故障识别结果，则融合前后的故障诊断结果如图7和图8所示。在变工况下，采用振动信号的诊断准确率相对于电流信号准确率较高，可以得知在跨转速诊断下，电流信号由于基频幅值较大，进而影响诊断准确率；而对于振动信号，没有此局限性，使得其在跨转速故障诊断效果表现较好。
[0104]
以上，仅为本技术较佳的具体实施方式，但本技术的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本技术揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本技术的保护范围之内。因此，本技术的保护范围应该以权利要求的保护范围为准。