CN110376522B - 一种数据融合的深度学习网络的电机故障诊断方法 - Google Patents

一种数据融合的深度学习网络的电机故障诊断方法 Download PDF

Info

Publication number
CN110376522B
CN110376522B CN201910827717.6A CN201910827717A CN110376522B CN 110376522 B CN110376522 B CN 110376522B CN 201910827717 A CN201910827717 A CN 201910827717A CN 110376522 B CN110376522 B CN 110376522B
Authority
CN
China
Prior art keywords
network
data
motor
vector
formula
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910827717.6A
Other languages
English (en)
Other versions
CN110376522A (zh
Inventor
王博涛
齐亚舟
张忠德
王孔照
朱耿超
沈传文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NINGXIA NORTHWEST HORSE ELECTRIC MANUFACTURING CO LTD
Xian Jiaotong University
Original Assignee
NINGXIA NORTHWEST HORSE ELECTRIC MANUFACTURING CO LTD
Xian Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by NINGXIA NORTHWEST HORSE ELECTRIC MANUFACTURING CO LTD, Xian Jiaotong University filed Critical NINGXIA NORTHWEST HORSE ELECTRIC MANUFACTURING CO LTD
Priority to CN201910827717.6A priority Critical patent/CN110376522B/zh
Publication of CN110376522A publication Critical patent/CN110376522A/zh
Application granted granted Critical
Publication of CN110376522B publication Critical patent/CN110376522B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/34Testing dynamo-electric machines
    • G01R31/343Testing dynamo-electric machines in operation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Control Of Electric Motors In General (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

一种数据融合的深度学习网络的电机故障诊断方法;该神经网络中包括数据压缩网络、特征提取网络和分类网络;其确定与训练方法为:1、采集电机的A,B两相电流信号和电机端部轴承的振动信号。进行数据标准化,通过希尔伯特‑黄变换得到频谱序列,建立神经网络的数据集;2、建立深度神经网络,确定网络结构并初始化参数;3、将训练集输入神经网络,分别计算不同神经网络的损耗函数,并利用损耗值更新神经网络参数;4、将测试集的数据输入神经网络,计算准确率,重复步骤3直至准确率满足要求;训练结束后,神经网络可以将输入的电流、振动数据映射至特征平面,分类网络根据其所在区域对应的故障状态可预测电机是否发生故障。

Description

一种数据融合的深度学习网络的电机故障诊断方法
技术领域
本发明涉及电气设备的故障诊断技术领域,具体涉及一种数据融合的深度学习网络的电机故障诊断方法,通过改变深度学习神经网络的结构,利用自动编码器进行数据融合,实现利用电流与振动数据实现综合判断电机故障。
背景技术
电机是现代工业、生产和生活中最常用的驱动设备,其使用范围十分广泛。在工业生产长期服役的情况下,电机可能发生轴承疲劳磨损、定子匝间短路、转子断条等故障。通过电机状态监测和故障诊断能够及时检修,保证电机的可靠运行。
随着近年深度学习的发展与硬件计算速度的提升,很多基于深度神经网络的故障诊断模型和算法被提出并应用于工业生产中。相对于传统的建模分析法,基于数据的故障诊断技术不会只考虑故障的主要因素,分析更加全面,能够通过已有的全部数据进行训练,学习信号中的规律。
通过数据融合技术,可利用不同类型的传感器测得的数据进行综合判断,提高故障辨识的准确率与可靠性。但由于神经网络通常采用多层全连接结构,特征融合方法不易实现。
发明内容
为提高电机故障监测的可靠性,本发明的目的在于提出一种数据融合的深度学习网络的电机故障诊断方法,利用深度神经网络学习数据特征,融合电机电流、振动数据并综合判断电机故障。
为达到上述目的,本发明所采用的技术方案是:
一种数据融合的深度学习网络的电机故障诊断方法,包括如下步骤:
步骤1:采集电机的A,B两相电流信号iA(t)、iB(t)和电机端部轴承的振动信号v(t),进行数据标准化,通过希尔伯特-黄变换得到A,B两相电流信号与振动信号的频谱序列iA(z)、iB(z)、v(z);建立神经网络的数据集:
对采集到的电流信号、振动信号进行标准化,使其均值为0,方差为1;对于A相电流式(1):
Figure BDA0002189645880000021
式中iAmean、iAvar为A相电流的平均值、方差;
Figure BDA0002189645880000022
为A相电流标准化后的结果;n为序列长度;
同样,根据式(2)、式(3)得到B相电流和振动数据标准化后的结果
Figure BDA0002189645880000023
v*(t);
Figure BDA0002189645880000024
式中iBmean、iBvar为B相电流的平均值、方差;
Figure BDA0002189645880000031
为B相电流标准化后的结果;
Figure BDA0002189645880000032
式中vmean、vvar为振动信号的平均值、方差;v*为振动信号标准化后的结果;
对标准化后的结果进行HHT变换:对于
Figure BDA0002189645880000033
首先对标准化后的结果进行经验模态分解并计算幅频特性函数:
a)
Figure BDA0002189645880000034
的极大值点经三次样条插值得到的上包络线记为Smax(t),极小值点经三次样条插值得到的下包络线记为Smin(t),则上、下包络线的中线m(t)表示为式(4):
Figure BDA0002189645880000035
b)根据式(5)计算
Figure BDA0002189645880000036
的本征模态函数h1(t):
Figure BDA0002189645880000037
c)若m(t)极大值、极小值个数均大于1,则重复a)、b),对m(t)进行进一步的经验模态分解,到高次的本征模态函数h2(t)、h3(t)、...hp(t),直至极大值或极小值个数为1,其中p为本征模态函数的个数;
d)将所有本征模态函数的幅频特性相加,得到原始信号的频域信号,如下式(6):
Figure BDA0002189645880000041
式中
Figure BDA0002189645880000042
表示Hilbert变换,Re(·)表示取实部,hk为第k个本征模态函数;
上述经验模态分解与幅频特性函数计算过程中,将
Figure BDA0002189645880000043
分别替换为
Figure BDA0002189645880000044
v*(t),则能够分别得到B相电流和振动数据的频域信号iB(z)、v(z);
将iA(z)、iB(z)连接为i(z),作为电流输入信号,信号长度为Ni;v(z)为振动输入信号,长度为Nv;一组电流、振动信号共同组成一个神经网络的输入;
若电机状态类型数为s,电机处于第j种状态时,状态标签向量表示为下式(7)
Figure BDA0002189645880000045
将多个时间段的电流、振动信号进行预处理,得到多个神经网络输入,建立数据集,并根据输入数据对应的电机状态通过式(7)确定训练标签;随机选取一半作为训练集,另外一半作为测试集;
步骤2:建立深度神经网络,确定网络结构并初始化参数:
深度神经网络包括数据压缩网络、特征提取网络和特征分类网络三个部分;
数据压缩网络采用等效自动编码器架构,对于电流数据压缩网络,输入、输出神经元个数为Ni,经过压缩的特征向量层的神经元个数为τNi,其中τ为常系数,取值范围为0.1~0.5;对于振动数据压缩网络,输入、输出神经元个数为Nv,经过压缩的特征向量层的神经元个数为τNv;从输入层到特征向量层共3~4层,神经元个数逐层递减,解码器与编码器对称;
特征提取网络以两个数据压缩网络的压缩结果为输入,神经元个数为τ(Ni+Nv),采用自动编码器结构,最中间的隐藏层神经元个数为2,编码部分层数为5,神经元数目逐层减少,解码器与编码器对称;
特征分类网络是一个Soft-max分类器,输入为二维特征潜向量[f1,f2]T;分类先通过一个两层的神经网络,将特征潜向量映射至每种故障的置信值vi,i=1,2,...,n,计算公式如式(8):
v=[v1,v2,…,vs]T=g(Ws×2·[f1,f2]T+bs) (8)
式中Ws×2、bs分别为神经网络的权重矩阵、偏置向量,g(·)为激活函数;
由公式(9)根据置信值计算电机发生每种故障的概率;
Figure BDA0002189645880000051
式中pi为被测电机发生第i种故障的概率;
建立好对应结构的神经网络后,生成均值为0、方差为1的随机数,将神经网络参数初始化;
步骤3:将步骤1建立的训练集输入神经网络,分别计算不同神经网络的损耗函数,并利用损耗值更新神经网络参数:
将训练集中的电流数据输入至电流数据压缩网络,该网络输出层为
Figure BDA0002189645880000052
特征向量层的为电流数据的特征向量fi(n),电流数据压缩网络的损失函数表示为式(10)
Figure BDA0002189645880000053
式中Li为电流数据压缩网络的损耗值,||·||2为向量的模;i为电流数据压缩网络输入的电流数据;
Figure BDA0002189645880000061
为电流数据压缩网络输出层输出的向量;
将训练集中的振动数据输入至振动数据压缩网络,该网络输出层为
Figure BDA0002189645880000062
特征向量层的为振动数据的特征向量fv(n),振动数据压缩网络的损失函数表示为式(11)
Figure BDA0002189645880000063
式中Lv为振动数据压缩网络的损耗值;v为振动数据压缩网络输入的振动数据;
Figure BDA0002189645880000064
为振动数据压缩网络输出层输出的向量;
根据损失函数Li,Lv,利用梯度下降法更新电流数据压缩网络与振动数据压缩网络的参数;
将fi(n)、fv(n)合并为fin(n),输入至特征提取网络,得到输出层fout(n)和二维特征潜向量[f1,f2]T;特征提取网络损失函数表示为式(12),并根据特征压缩网络的损耗值Lf更新该网络的参数;
Lf=||fin(n)-fout(n)||2 (12)
将二维的特征潜向量输入特征分类网络,得到最终预测不同故障的概率向量p=[pi,p2,…,pn]T,则概率最大的状态为预测的电机状态类型如式(13),并计算特征分类网络的损耗函数Lp如式(14),根据损耗函数更新分类网络的参数;
Figure BDA0002189645880000065
Lp=||q-l||2 (14)
式中q为故障预测向量,qk是q的k个元素,是对第k种故障的预测结果,若qk=1表示输入信号是第i种故障,若qk=0表示输入信号不是第i种故障;
步骤4:将测试集的数据输入神经网络:电流数据输入电流数据压缩网络,提取电流数据特征;振动数据输入振动数据压缩网络,提取振动数据特征;将该两网络的输出合并,作为特征提取网络的输入,提取输入数据的特征潜向量;由特征分类网络计算故障概率向量,并由式(13)计算故障预测向量q,判断电机状态;
对于一组数据,其对应电机的真实状态的标签向量l由式(7)计算得到,若q=l,则神经网络正确判断了电机的状态,否则判断错误;计算准确率a如式(15)
Figure BDA0002189645880000071
式中m为测试集样本的总数,mc为故障判断正确的个数;
提前设定常数amin作为准确率阈值低限,若a<amin,则重复步骤3,直至准确率满足a≥amin时训练结束;
训练结束后,通过数据压缩网络和特征提取网络,将电机的电流、振动数据映射为二维的特征潜向量;若将二维向量表示在平面直角坐标系中,该直角坐标系即能够作为电机状态的特征平面;特征分类网络就是根据特征潜向量在特征平面的分布,预测电机是否发生故障;将不同输入数据的特征潜向量用离散的点表示在特征平面,即实现故障预测的可视化。
和现有技术相比较,本发明具备如下优点:
本发明采用深度神经网络进行信号特征的提取与分类,避免了复杂的建模过程,直接对数据进行分类,可以减小建模过程考虑不全面、忽略次要因素而造成的误差;采用了数据融合的神经网络结构,一般的神经网络为了考虑不同类型的数据,只能将他们连接为一个向量并从输入层输入,这样的提取特征值计算成本高、特征向量准确性差、效率低,本发明在步骤2中建立了数据压缩网络和特征融合网络,先通过数据压缩网络对不同类型的数据进行降维,再将压缩后的数据输入特征压缩网络提取用于分类的特征潜向量,这样减小了计算量、提高了数据融合的效率;传统的基于神经网络的故障检测只能简单的给出对监测对象是否发生故障,难以直观的观察,而本发明中采用的特征分类网络最中间一层的向量维数为2,因此可以将特征潜向量的训练结果表示在平面直角坐标系中,并将平面分为不同的区域,对应不同的故障类型,在特征潜向量跃变到不同区域时即表示监测对象发生故障,实现故障判断的可视化。
附图说明
图1为振动信号EMD分解为本征模态函数示例,其中,图1a到图1h分别为第一本征模态函数h1(t)到第八本征模态函数h1(t)。
图2为故障诊断网络总体结构。
图3为数据压缩网络结构,其中,图3a为电流数据压缩网络,图3b为振动数据压缩网络。
图4为特征潜向量提取网络。
图5为Soft-max分类器结构。
图6为训练集和测试集的特征潜向量在特征平面分布,其中,图6a为训练集的特征潜向量在特征平面分布,图6b为测试集的特征潜向量在特征平面分布。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细说明。
根据以上对数据融合的深度学习网络的电机故障诊断方法的详细介绍,下面以同步电机轴承故障的开源数据集为例说明本发明的具体实施方式。
开源数据由帕德博恩大学轴承中心提供。该数据集对33个不同状态的电机进行测试,包括正常、轴承内环故障、轴承外环故障三种状态。不同电机故障的严重程度、形成方式不同。数据集中测量了电机在不同转速、负载下的两相电流数据与振动数据,可用于故障判断。按照表1选定数据集,选取部分人工制造的故障电机作为神经网络的训练集,选取部分加速老化的故障电机作为测试集,验证本发明方案的故障辨识能力。
表1数据集选择
Figure BDA0002189645880000091
步骤1:取不同轴承编号对应的电机的数据,截取同一时刻对应的两相电流数据号iA(t)、iB(t)和电轴承振动信号v(t);进行数据标准化,通过希尔伯特-黄变换得到A,B两相电流信号与振动信号的频谱序列iA(z)、iB(z)、v(z);建立神经网络的数据集:
截取电流数据号和电轴承振动信号,每次截取数据2000点的数据;对采集到的电流信号、振动信号进行标准化,使其均值为0,方差为1;对于A相电流式(1):
Figure BDA0002189645880000101
式中n=2000;
同样,根据式(2)、式(3)得到B相电流和振动数据标准化后的结果
Figure BDA0002189645880000102
v*(t);
Figure BDA0002189645880000103
Figure BDA0002189645880000104
对标准化后的数据进行HHT变换:对于
Figure BDA0002189645880000105
首先对标准化后的结果进行经验模态分解并计算幅频特性函数:
a)
Figure BDA0002189645880000106
的极大值点经三次样条插值得到的上包络线记为Smax(t),极小值点经三次样条插值得到的下包络线记为Smin(t),则上、下包络线的中线m(t)表示为式(4):
Figure BDA0002189645880000111
b)根据式(5)计算
Figure BDA0002189645880000112
的本征模态函数h1(t):
Figure BDA0002189645880000113
c)若m(t)极大值、极小值个数均大于1,则重复a)、b),对m(t)进行进一步的经验模态分解,到高次的本征模态函数h2(t)、h3(t)、...hp(t),直至极大值或极小值个数为1,其中p为本征模态函数的个数;一个振动数据样本的本征模态分解结果如图1所示,图1a到图1h分别为第一本征模态函数h1(t)到第八本征模态函数h8(t);
d)将所有本征模态函数的幅频特性相加,得到原始信号的频域信号,如下式(6):
Figure BDA0002189645880000114
上述经验模态分解与幅频特性函数计算过程中,将
Figure BDA0002189645880000115
分别替换为
Figure BDA0002189645880000116
v*(t),则能够分别得到B相电流和振动数据的频域信号iB(z)、v(z);
其中iA(z)、iB(z)分别截取前400点的低频的有效成分,连接为i(z),作为电流输入信号,信号长度为800;截取振动信号前1000点的低频的有效成分作为v(z)为振动输入信号,长度为1000;一组电流、振动信号共同组成一个神经网络的输入,信号长度为1800;
开源数据集中电机状态类型数为3,若电机处于第j种状态时,状态向量标签可以表示为下式(7)
Figure BDA0002189645880000121
其中s=3;
因此正常电机、内环故障电机、外环故障电机输入信号的标签分别为[1,0,0]T、[0,1,0]T、[0,0,1]T;选择与输入信号对应的标签,作为一组数据;
选取不同时间片段的信号进行上述预处理操作,使得每个编号的轴承都有1600组数据,用于训练、测试神经网络,建立规模为8000的训练数据集和规模为8000的训练数据集;
步骤2:建立深度神经网络,确定网络结构并初始化参数:
深度神经网络总体结构如图2,包括数据压缩网络、特征提取网络和特征分类网络三个部分;
数据压缩网络采用等效自动编码器架构,其结构如图3;如图3a,对于电流数据压缩网络,输入、输出神经元个数为1000,经过压缩的特征向量层的神经元个数为750;如图3b,对于振动数据压缩网络,输入、输出神经元个数为800,经过压缩的特征向量层的神经元个数为400;
数据压缩网络从输入层到特征向量层,神经元个数逐层递减;电流数据压缩网络从输入层到特征向量层,每层神经元个数为1000、900、750,解码器与编码器对称;振动数据压缩网络从输入层到特征向量层,每层神经元个数为800、700、600,解码器与编码器对称;
特征提取网络结构如图4,以两个数据压缩网络的压缩结果为输入,神经元个数为1350,采用自动编码器结构,最中间的隐藏层神经元个数为2,编码部分层数为5,神经元数目逐层减少,每层神经元个数为1350、700、500、350、2,解码器与编码器对称;
分类网络是一个Soft-max分类器,其结构如图5,输入为二维特征潜向量[f1,f2]T;分类先通过一个两层的神经网络,将特征潜向量映射至每种故障的置信值vi,i=1,2,...,n,计算公式如式(8):
v=[v1,v2,…,vs]T=g(Ws×2·[f1,f2]T+bs) (8)
由公式(9)根据置信值计算电机发生每种故障的概率;
Figure BDA0002189645880000131
建立好对应结构的神经网络后,生成均值为0、方差为1的随机数,将神经网络参数初始化;
步骤3:将步骤1建立的训练集输入神经网络,分别计算不同神经网络的损耗函数,并利用损耗值更新神经网络参数:
将训练集中的电流数据输入至电流数据压缩网络,该网络输出层为
Figure BDA0002189645880000132
特征向量层的为电流数据的特征向量fi(n),根据式(10)计算电流数据压缩网络的损失函数;
Figure BDA0002189645880000133
将训练集中的振动数据输入至振动数据压缩网络,该网络输出层为
Figure BDA0002189645880000134
特征向量层的为振动数据的特征向量fv(n),根据式(11)计算振动数据压缩网络的损失函数;
Figure BDA0002189645880000135
根据损失函数Li,Lv,利用梯度下降法更新电流数据压缩网络与振动数据压缩网络的参数;
将fi(n)、fv(n)合并为fin(n),输入至特征提取网络,得到输出层fout(n)和二维的特征潜向量[f1,f2]T;特征提取网络损失函数可表示为式(9),并根据特征压缩网络的损耗值Lf更新该网络的参数;
Lf=||fin(n)-fout(n)||2 (9)
将二维的特征潜向量输入特征分类网络,得到最终预测不同故障的概率向量p=[pi,p2,…,pn]T,则概率最大的状态为预测的电机状态类型如式(13),并计算特征分类网络的损耗函数Lp如式(14),根据损耗函数更新分类网络的参数;
Figure BDA0002189645880000141
Lp=||q-l||2 (14)
步骤4:对于一组数据,其对应电机的真实状态的标签向量l由式(7)计算得到,若q=l,则神经网络正确判断了电机的状态,否则判断错误;计算准确率a如式(15)
Figure BDA0002189645880000142
提前设定常数amin=0.8作为准确率阈值低限,若a<amin,则重复步骤3,直至准确率满足a≥amin时训练结束;
训练结束后,通过数据压缩网络和特征提取网络,将电机的电流、振动数据映射为二维的特征潜向量;若将二维向量表示在平面直角坐标系中,该直角坐标系即能够作为电机状态的特征平面;特征分类网络就是根据特征潜向量在特征平面的分布,预测电机是否发生故障;将不同输入数据的特征潜向量用离散的点表示在特征平面,即实现故障预测的可视化;
图6中图6a为训练结束后训练集特征潜向量的分布,同种灰度的点代表同一种故障状态的向量;经过训练集对神经网络的优化,平面直角坐标系被不同的特征向量点分割成不同的区域,对应不同的故障状态;测试集被输入神经网络后,也被映射到该特征平面;根据特征潜向量所在的区域,根据区域对应的故障状态可以预测输入信号包含的故障类型;图6b为测试集特征潜向量的分布;可见其分布与测试集类似,同一故障类型的特征潜向量所在的区域与训练集的近乎相同;根据预测的故障类型和实际标签对比,可计算出该故障诊断系统的准确率达到84%。
可见,本发明提出的数据融合的深度学习网络的电机故障诊断方法,可以根据已有的数据集训练、优化神经网络参数,使其将输入信号映射到特征平面,将其分割成不同的故障区域;当待测电机信号输入网络,根据其潜向量所在的区域,可以由电流和机械振动信号综合预测其故障类型。

Claims (1)

1.一种数据融合的深度学习网络的电机故障诊断方法,其特征在于:包括如下步骤:
步骤1:采集电机的A,B两相电流信号iA(t)、iB(t)和电机端部轴承的振动信号v(t),进行数据标准化,通过希尔伯特-黄变换得到A,B两相电流信号与振动信号的频谱序列iA(z)、iB(z)、v(z);建立神经网络的数据集:
对采集到的电流信号、振动信号进行标准化,使其均值为0,方差为1;对于A相电流式(1):
Figure FDA0002946175800000011
式中iAmean、iAvar为A相电流的平均值、方差;
Figure FDA0002946175800000012
为A相电流标准化后的结果;n为序列长度;
同样,根据式(2)、式(3)得到B相电流和振动数据标准化后的结果
Figure FDA0002946175800000013
v*(t);
Figure FDA0002946175800000014
式中iBmean、iBvar为B相电流的平均值、方差;
Figure FDA0002946175800000015
为B相电流标准化后的结果;
Figure FDA0002946175800000021
式中vmean、vvar为振动信号的平均值、方差;v*为振动信号标准化后的结果;
对标准化后的结果进行HHT变换:对于
Figure FDA0002946175800000022
首先对标准化后的结果进行经验模态分解并计算幅频特性函数:
a)
Figure FDA0002946175800000023
的极大值点经三次样条插值得到的上包络线记为Smax(t),极小值点经三次样条插值得到的下包络线记为Smin(t),则上、下包络线的中线m(t)表示为式(4):
Figure FDA0002946175800000024
b)根据式(5)计算
Figure FDA0002946175800000025
的本征模态函数h1(t):
Figure FDA0002946175800000026
c)若m(t)极大值、极小值个数均大于1,则重复a)、b),对m(t)进行进一步的经验模态分解,到高次的本征模态函数h2(t)、h3(t)、...、hp(t),直至极大值或极小值个数为1,其中p为本征模态函数的个数;
d)将所有本征模态函数的幅频特性相加,得到原始信号的频域信号,如下式(6):
Figure FDA0002946175800000027
式中
Figure FDA0002946175800000028
表示Hilbert变换,Re(·)表示取实部,hk为第k个本征模态函数;
上述经验模态分解与幅频特性函数计算过程中,将
Figure FDA0002946175800000031
分别替换为
Figure FDA0002946175800000032
v*(t),则能够分别得到B相电流和振动数据的频域信号iB(z)、v(z);
将iA(z)、iB(z)连接为i(z),作为电流输入信号,信号长度为Ni;v(z)为振动输入信号,长度为Nv;一组电流、振动信号共同组成一个神经网络的输入;
若电机状态类型数为s,电机处于第j种状态时,状态标签向量表示为下式(7)
Figure FDA0002946175800000033
将多个时间段的电流、振动信号进行预处理,得到多个神经网络输入,建立数据集,并根据输入数据对应的电机状态通过式(7)确定训练标签;随机选取一半作为训练集,另外一半作为测试集;
步骤2:建立深度神经网络,确定网络结构并初始化参数:
深度神经网络包括数据压缩网络、特征提取网络和特征分类网络三个部分;
数据压缩网络采用等效自动编码器架构,对于电流数据压缩网络,输入、输出神经元个数为Ni,经过压缩的特征向量层的神经元个数为τNi,其中τ为常系数,取值范围为0.1~0.5;对于振动数据压缩网络,输入、输出神经元个数为Nv,经过压缩的特征向量层的神经元个数为τNv;从输入层到特征向量层共3~4层,神经元个数逐层递减,解码器与编码器对称;
特征提取网络以两个数据压缩网络的压缩结果为输入,神经元个数为τ(Ni+Nv),采用自动编码器结构,最中间的隐藏层神经元个数为2,编码部分层数为5,神经元数目逐层减少,解码器与编码器对称;
特征分类网络是一个Soft-max分类器,输入为二维特征潜向量[f1,f2]T;分类先通过一个两层的神经网络,将特征潜向量映射至每种故障的置信值vi,i=1,2,...,n,计算公式如式(8):
v=[v1,v2,…,vs]T=g(Ws×2·[f1,f2]T+bs) (8)
式中Ws×2、bs分别为神经网络的权重矩阵、偏置向量,g(·)为激活函数;
由公式(9)根据置信值计算电机发生每种故障的概率;
Figure FDA0002946175800000041
式中pi为被测电机发生第i种故障的概率;
建立好对应结构的神经网络后,生成均值为0、方差为1的随机数,将神经网络参数初始化;
步骤3:将步骤1建立的训练集输入神经网络,分别计算不同神经网络的损耗函数,并利用损耗值更新神经网络参数:
将训练集中的电流数据输入至电流数据压缩网络形成网络输出层,该网络输出层为
Figure FDA0002946175800000042
特征向量层的为电流数据的特征向量fi(n),电流数据压缩网络的损失函数表示为式(10)
Figure FDA0002946175800000043
式中Li为电流数据压缩网络的损耗值,||·||2为向量的模;i为电流数据压缩网络输入的电流数据;
Figure FDA0002946175800000044
为电流数据压缩网络输出层输出的向量;
将训练集中的振动数据输入至振动数据压缩网络形成网络输出层,该网络输出层为
Figure FDA0002946175800000051
特征向量层的为振动数据的特征向量fv(n),振动数据压缩网络的损失函数表示为式(11)
Figure FDA0002946175800000052
式中Lv为振动数据压缩网络的损耗值;v为振动数据压缩网络输入的振动数据;
Figure FDA0002946175800000053
为振动数据压缩网络输出层输出的向量;
根据损失函数Li,Lv,利用梯度下降法更新电流数据压缩网络与振动数据压缩网络的参数;
将fi(n)、fv(n)合并为fin(n),输入至特征提取网络,得到输出层fout(n)和二维特征潜向量[f1,f2]T;特征提取网络损失函数表示为式(12),并根据特征压缩网络的损耗值Lf更新该网络的参数;
Lf=||fin(n)-fout(n)||2 (12)
将二维的特征潜向量输入特征分类网络,得到最终预测不同故障的概率向量p=[pi,p2,…,pn]T,则概率最大的状态为预测的电机状态类型如式(13),并计算特征分类网络的损耗函数Lp如式(14),根据损耗函数更新分类网络的参数;
Figure FDA0002946175800000054
Lp=||q-l||2 (14)
式中q为故障预测向量,qk是q的k个元素,是对第k种故障的预测结果,若qk=1表示输入信号是第i种故障,若qk=0表示输入信号不是第i种故障;
步骤4:将测试集的数据输入神经网络:电流数据输入电流数据压缩网络,提取电流数据特征;振动数据输入振动数据压缩网络,提取振动数据特征;将该两网络的输出合并,作为特征提取网络的输入,提取输入数据的特征潜向量;由特征分类网络计算故障概率向量,并由式(13)计算故障预测向量q,判断电机状态;
对于一组数据,其对应电机的真实状态的标签向量l由式(7)计算得到,若q=l,则神经网络正确判断了电机的状态,否则判断错误;计算准确率a如式(15)
Figure FDA0002946175800000061
式中m为测试集样本的总数,mc为故障判断正确的个数;
提前设定常数amin作为准确率阈值低限,若a<amin,则重复步骤3,直至准确率满足a≥amin时训练结束;
训练结束后,通过数据压缩网络和特征提取网络,将电机的电流、振动数据映射为二维的特征潜向量;若将二维向量表示在平面直角坐标系中,该直角坐标系即能够作为电机状态的特征平面;特征分类网络就是根据特征潜向量在特征平面的分布,预测电机是否发生故障;将不同输入数据的特征潜向量用离散的点表示在特征平面,即实现故障预测的可视化。
CN201910827717.6A 2019-09-03 2019-09-03 一种数据融合的深度学习网络的电机故障诊断方法 Active CN110376522B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910827717.6A CN110376522B (zh) 2019-09-03 2019-09-03 一种数据融合的深度学习网络的电机故障诊断方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910827717.6A CN110376522B (zh) 2019-09-03 2019-09-03 一种数据融合的深度学习网络的电机故障诊断方法

Publications (2)

Publication Number Publication Date
CN110376522A CN110376522A (zh) 2019-10-25
CN110376522B true CN110376522B (zh) 2021-07-16

Family

ID=68261343

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910827717.6A Active CN110376522B (zh) 2019-09-03 2019-09-03 一种数据融合的深度学习网络的电机故障诊断方法

Country Status (1)

Country Link
CN (1) CN110376522B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112903001A (zh) * 2019-12-03 2021-06-04 财团法人纺织产业综合研究所 织物定型机的操作方法
CN111504676B (zh) * 2020-04-23 2021-03-30 中国石油大学(北京) 基于多源监控数据融合的设备故障诊断方法、装置及系统
CN111582463B (zh) * 2020-06-08 2024-02-09 佛山金华信智能科技有限公司 伺服电机故障识别及模型训练方法、装置、介质及终端
CN114200302A (zh) * 2020-09-18 2022-03-18 株洲中车时代电气股份有限公司 一种电机定子线圈匝间短路故障在线诊断方法和系统
CN112199295B (zh) * 2020-10-29 2023-05-05 北京信息科技大学 一种基于频谱的深度神经网络缺陷定位方法及系统
CN112839349B (zh) * 2020-12-25 2022-10-25 朗坤智慧科技股份有限公司 基于5g技术的角域平稳信号消噪方法、装置和网络服务端
CN112834959A (zh) * 2021-01-04 2021-05-25 福州大学 基于高频特征提取的直流配电系统故障检测方法
CN114332621B (zh) * 2021-12-30 2023-06-13 华智生物技术有限公司 一种基于多模型特征融合的病虫害识别方法及系统
CN114509158B (zh) * 2022-01-04 2023-05-30 东南大学 一种声振融合的叶片裂纹故障检测方法及应用
CN114781552B (zh) * 2022-06-17 2022-09-02 深圳硅山技术有限公司 电机性能测试方法、装置、设备及存储介质
CN116738354B (zh) * 2023-08-15 2023-12-08 国网江西省电力有限公司信息通信分公司 一种电力物联网终端行为异常检测方法及系统

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3182114A1 (en) * 2015-12-15 2017-06-21 General Electric Company Partial discharge monitoring of electrical machines using acoustic emission sensors and electrical sensors

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2476246A (en) * 2009-12-15 2011-06-22 Univ Portsmouth Diagnosing an operation mode of a machine
WO2015047121A1 (en) * 2013-09-25 2015-04-02 Siemens Aktiengesellschaft Method and apparatus for embedded current signature analysis and remote condition monitoring for industrial machinery
CN104166095B (zh) * 2014-08-29 2017-02-15 东南大学 一种基于双边型直线电机的故障信息融合诊断方法
CN206515445U (zh) * 2016-11-22 2017-09-22 中国石油化工股份有限公司天然气榆济管道分公司 电机故障诊断系统
CN109873581A (zh) * 2017-12-04 2019-06-11 北京自动化控制设备研究所 一种基于神经网络的永磁同步电机直接速度控制方法
CN108931724A (zh) * 2018-07-30 2018-12-04 袁小芳 一种伺服电机故障诊断方法
CN110119397B (zh) * 2019-04-18 2023-06-30 东南大学 一种同时实现数据异常检测及数据压缩的深度学习方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3182114A1 (en) * 2015-12-15 2017-06-21 General Electric Company Partial discharge monitoring of electrical machines using acoustic emission sensors and electrical sensors

Also Published As

Publication number Publication date
CN110376522A (zh) 2019-10-25

Similar Documents

Publication Publication Date Title
CN110376522B (zh) 一种数据融合的深度学习网络的电机故障诊断方法
CN111914883B (zh) 一种基于深度融合网络的主轴轴承状态评估方法及装置
CN110006645B (zh) 一种多源融合的高压断路器机械故障诊断方法
CN111562108A (zh) 一种基于cnn和fcmc的滚动轴承智能故障诊断方法
CN111504676A (zh) 基于多源监控数据融合的设备故障诊断方法、装置及系统
CN113505655A (zh) 面向数字孪生系统的轴承故障智能诊断方法
CN111678699B (zh) 一种面向滚动轴承早期故障监测与诊断方法及系统
CN109948194B (zh) 一种高压断路器机械缺陷集成学习诊断方法
CN112926728B (zh) 一种永磁同步电机小样本匝间短路故障诊断方法
CN112305388B (zh) 一种发电机定子绕组绝缘局部放电故障在线监测诊断方法
CN114386537A (zh) 基于CatBoost的锂电池故障诊断方法、装置及电子设备
CN112949402A (zh) 极小故障样本量下行星齿轮箱故障诊断方法
CN112082639A (zh) 发动机状态的诊断方法及其诊断建模方法
CN113569990A (zh) 面向强噪声干扰环境的演艺装备故障诊断模型构建方法
CN114993669B (zh) 多传感器信息融合的传动系统故障诊断方法及系统
Dong Combining unsupervised and supervised learning for asset class failure prediction in power systems
CN115905832A (zh) 一种基于振动信号聚类识别的变压器故障诊断方法
CN114048787B (zh) 一种基于Attention CNN模型的轴承故障实时智能诊断方法与系统
CN115238785A (zh) 基于图像融合和集成网络的旋转机械故障诊断方法和系统
CN113092044B (zh) 一种基于加权水平可见图的旋转机械故障诊断方法
CN117076935B (zh) 数字孪生辅助的机械故障数据轻量级生成方法及系统
CN115791174B (zh) 一种滚动轴承异常诊断方法、系统、电子设备及存储介质
CN109784777B (zh) 基于时序信息片段云相似度度量的电网设备状态评估方法
CN116842379A (zh) 一种基于DRSN-CS和BiGRU+MLP模型的机械轴承剩余使用寿命预测方法
CN116956215A (zh) 一种传动系统故障诊断方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant