实验研究

基于支持向量机与粗糙集的隔爆电动机故障诊断

马宪民1, 张兴1, 张永强1,2

(1.西安科技大学 电气与控制工程学院, 陕西 西安 710054;2.神华宁煤集团 矿山机械制造维修分公司, 宁夏 银川 750001)

摘要:针对煤矿井下隔爆电动机故障数据获取难且故障数据杂乱、非线性等问题,提出了一种基于支持向量机与粗糙集的隔爆电动机故障诊断方法。该方法采用小波包对隔爆电动机定子瞬时功率进行频谱分析,并提取故障特征量;利用粗糙集的约简特性消除故障特征量冗余数据,将约简后的故障特征量作为支持向量机的输入样本,实现隔爆电动机转子故障诊断和分类。仿真结果表明,该方法故障诊断结果准确率达到92.857 1%。

关键词:隔爆电动机; 故障诊断; 支持向量机; 粗糙集; 瞬时功率; 小波包

0 引言

在煤矿产业中,隔爆电动机作为重要的动力部件,在生产、传输、通风、排水、掘进、破碎等工作中频繁启动、制动,加上煤矿工作环境恶劣,因此,故障发生率相对较高,占井下机电故障的60%以上[1],且故障种类复杂多变。在隔爆电动机故障类型中,转子断条故障和转子偏心故障最常见。常用的转子故障检测方法有定子电流信号检测法、振动信号检测法、轴电压检测法、轴向漏磁检测法等[2-4]。定子电流信号检测法因具有非侵入式检测的特点而被普遍应用,但当隔爆电动机转子发生断条故障时,因转差率较小,定子电流信号产生的故障特征频率会被基波频率淹没,所以,采集定子电流时,也要对相应的电压进行采集,求取隔爆电动机的瞬时功率,并对瞬时功率故障信号进行分析。传统的信号分析法为傅里叶变换法,但快速傅里叶变换频谱分析法只对平稳信号处理有效,不能对非平稳信号进行完全分析。因此,本文采用具有多分辨率和时频局部化特点并可对非平稳、突变和非连续信号进行分解的小波包分析法提取隔爆电动机转子故障特征量[5]

煤矿井下环境恶劣,数据采集困难且采集的数据量有限,而采集的数据中可能存在冗余或不准确数据,这些数据不仅影响故障诊断速度,而且会降低诊断准确率,因此,本文首先采用具有约简特点的粗糙集(Rough Set, RS)理论对故障特征量进行约简,消除冗余;再将约简后的故障特征量作为支持向量机(Support Vector Machine, SVM)的输入样本,对隔爆电动机转子故障进行诊断和分类。基于SVM与RS的隔爆电动机故障诊断方法总体框图如图1所示。

图1 基于SVM与RS的隔爆电动机故障诊断方法总体框图

1 故障特征提取

1.1 瞬时功率信号频谱分析

当隔爆电动机转子发生断条故障时,定子电流会产生相应的变化,进而反映在转子的边频故障特征频率fb中,fb可以表示为[5-7]

fb=(1±2ks)f

(1)

式中:f为基频频率;k=1,2,…;s为隔爆电动机转差率。

当隔爆电动机转子发生偏心故障时,定子电流会产生相应的变化,进而反映在转子的故障特征频率fecc中,fecc可以表示为

(2)

式中:p为隔爆电动机极对数;m=1,2,…;fr为转子的旋转频率。

隔爆电动机转子发生断条故障时,取k=1,此时的定子电流特征量幅值变化最明显,因此,反映在转子上的瞬时功率表达式为

I1[cos((2+2s)wt-φ1)+cos(2swt-φ1)]+

I2[cos((2-2s)wt-φ2)+cos(2swt-φ2)]}

(3)

式中:U为基波电压幅值;I为基波电流幅值;w为基波电压角频率;φ0为转子基波电流分量对应的初相位;I1I2为转子断条故障后对应的转子电流故障特征分量的幅值;φ1φ2为转子电流故障特征分量对应的初相位值。

从式(3)可以看出,特征分量2sf完全可以作为判断断条故障的依据。

当隔爆电动机转子偏心时,取m=1,此时反映在转子上的瞬时功率表达式为

(4)

式中为转子断条故障后对应的故障特征转子电流分量的幅值;为故障特征转子电流分量对应的初相位值;wr为对应fr产生的角频率。

从式(4)可以看出,在这些特征分量中,fr可以作为判断转子偏心故障的依据。

1.2 小波包分析法

多分辨分析是根据不相同的尺度因子i(iZ)把Hilbert空间L2(R)分解成所有子空间Wi(Wi是小波函数Ψ(t)的小波子空间)的正交和,。用一个新的子空间统一表征小波子空间Wi与尺度子空间Vi,则Hilbert空间中的正交分解Vi+1=Wi+Vi,可用分解统一表征为[8-10]

(5)

1.3 基于瞬时功率的故障诊断应用

给定型号为YBSS-200G的破碎机三相异步隔爆电动机,转差率s为0.02,极对数p为2,电压为3 300 V,额定电流为43 A,额定转速为1 483 r/min。对隔爆电动机的转子运行状态进行频谱分析,如图2所示。

图2可以反映电动机转子故障特征频率的变化,但是没有将反映故障的特征频率位置具体化。因此,为了更加清楚地反映电动机转子故障特征频率,对隔爆电动机的转子瞬时功率采用小波基sym3进行4层分解,对分解后的信号系数[4,0],[4,2]重构,最后求得隔爆电动机转子断1根并经过小波包处理后的频谱,如图3所示。

(a) 转子正常时的瞬时功率采样信号

(b) 转子偏心时的瞬时功率采样信号

(c) 转子断条和偏心时的瞬时功率采样信号

(d) 转子断条为1根时的瞬时功率采样信号

(e) 转子断条为2根时的瞬时功率采样信号

图2 电动机转子状态瞬时功率频谱

(a) 转子正常时的瞬时功率采样信号

(b) 转子偏心时的瞬时功率采样信号

(c) 转子断条和偏心时的瞬时功率采样信号

(d) 转子断条为1根时的瞬时功率采样信号

(e) 转子断条为2根时的瞬时功率采样信号

图3 经小波包重构后的电动机转子状态瞬时功率频谱

从图3可以很清晰地看出,在转子发生断条时,分别在2sf,4sf处出现故障特征量波动。转子发生偏心故障时瞬时功率分别在fr,2f-fr处出现故障特征量波动。

对图3进行分析并提取不同频段的隔爆电动机数据样本,这些数据样本中可能存在不准确、冗余的故障样本,这些样本不仅影响诊断速度,更会降低诊断准确性。因此,在对数据样本进行故障诊断分类之前,先利用RS理论的约简特性对数据样本进行预处理,在不影响数据分类结果的前提下,消除冗余条件属性,剔除重复、缺陷样本。

2 RS约简与SVM分类原理

2.1 RS约简基本原理

RS理论的知识表达系统即信息系统,可表达为四元组S=(H,B,R,l),其中(H,B)组成一个决策表,H为对象集合,H={y1,y2,…,ym};B为属性集合,B={b1,b2,…,bn}=CDC为条件属性,D为决策属性,且D={d};R为属性的值域;l表示一种映射:H×BR,反映对象集合的一种属性值。在一个信息系统中,每个条件属性之间通常存在一定程度的依赖关系,这种依赖关系可表示为[11-12]

(6)

式中posP(Q)为QP的正域。

RS理论虽然可以根据约简后数据的潜在规则对新的数据进行推理分类,但这种规则更多地依赖知识库的逻辑推理,规则之间不存在可寻的关系,因此,这种诊断方法速度低、效率不高。引入RS算法后,会使样本数据大量减少,所以对于具有局限性的专家系统及存在经验风险最小化原则的人工神经网络而言,由于数据样本太少,误判率就会大大增加。因此,本文采用基于结构风险最小化原则,并可对小样本、非线性数据进行较好诊断分类的SVM算法。

2.2 SVM分类原理

假定有训练样本F={(xi,yi)},xRd,yi=1,2,…,F能被超平面w·x+b=0正确分成2类,且使样本集满足式(7)[13]

yi[(w·xi)+b]-1≥0

(7)

式中:w为权值;b为阈值。

在Lagrange乘子αi满足αi>0的约束条件下,构建最优超平面的目标函数αi,即

(8)

对于线性不可分问题[14-15],在线性问题的基础上引入一个松弛变量ξi(ξi>0),当求解的αi为最优解时可得最优分类函数为

(9)

3 故障诊断分类实例分析

选取型号为YBSS-200G的破碎机隔爆电动机进行实例分析。设E1表示隔爆电动机转子断条为1根,E2表示转子断条为2根,E3表示转子偏心故障,E4表示转子断条为1根及偏心故障,选取≤0.04f,0.08f,0.49f,1.51f,1.92f,1.98f,2.04f,≥2.08f这些频段或频率的故障数据对原始故障数据进行离散化。隔爆电动机转子出现故障用1表示,无故障用0表示,可得故障决策表,见表1。

表1 故障决策表

样本≤0.04f(C1)0.08f(C2)0.49f(C3)1.51f(C4)1.92f(C5)1.98f(C6)2.04f(C7)≥2.08f(C8)故障类型Y110000110E1Y210000111E1Y300000110E1Y400000111E1Y510000110E1Y600000111E1Y711001111E2Y811001110E2Y901001111E2Y1001001110E2Y1111001111E2Y1201001111E2Y1300110001E3Y1400110000E3Y1500110001E3Y1600110000E3Y1710110110E4Y1810110111E4Y1900110110E4Y2000110111E4Y2110110110E4Y2200110111E4

表1中,U={Y1,Y2,Y3,Y4,Y5,Y6,Y7,Y8,Y9,Y10,Y11,Y12,Y13,Y14,Y15,Y16,Y17,Y18,Y19,Y20,Y21,Y22};条件属性C={C1,C2,C3,C4,C5,C6,C7,C8};决策属性D={E1,E2,E3,E4}。U关于C的等价关系:

U/C={{Y1,Y5},{Y2},{Y3},{Y4,Y6},{Y7,Y11},{Y8},{Y9,Y12},{Y10},{Y13,Y15},{Y14,Y16},{Y17,Y21},{Y18},{Y19},{Y20,Y22}}

(10)

从式(10)可以看出,在同一个决策属性E1中,样本Y1Y5重复,样本Y4Y6重复,消除样本Y5Y6;同样的,在同一个决策属性E2中,消除样本Y11Y12;在同一个决策属性E3中,消除样本Y15Y16;在同一个决策属性E4中,消除样本Y21Y22

在条件属性中,C2C5C3C4C6C7对于各个决策属性的表现状态相同,所以去除条件属性C5C4C7。约简后的故障决策表见表2。

将约简后的故障决策表作为SVM训练样本,对SVM进行训练,经过多次训练后,选取径向基核函数k(x,y)作为SVM训练模型,误差惩罚因子g取值为10。

表2 约简后的故障决策表

样本C1C2C3C6C8故障类型Y110010E1Y210011E1Y300010E1Y400011E1Y711011E2Y811010E2Y901011E2Y1001010E2Y1300101E3Y1400100E3Y1710110E4Y1810111E4Y1900110E4Y2000111E4

(11)

式中:x,y为训练输入;d为径向基核函数的宽度,取为0.5。

根据SVM的训练模型结果,选取44组测试样本故障数据,按RS约简方法对测试数据进行约简,得到最优决策表,将最优决策表作为SVM的测试样本集进行故障预测,预测结果如图4所示。其中类别标签1,2,3,4分别代表4种故障类型E1E4。从图4可看出,只有2组样本的分类与实际不符,故障预测准确度达到92.857 1%。

图4 经RS处理后的SVM仿真结果

将未经过RS处理的44组测试数据进行SVM仿真,结果如图5所示。从图5可看出,有6组样本与实际不相符,故障预测准确度为86.363 6%。与图4的诊断结果相比,经过RS处理后的故障诊断结果准确率更高。

图5 未经RS处理的SVM仿真结果

4 结语

采用小波包对隔爆电动机转子瞬时功率进行频谱仿真,可避免故障边频被基波淹没的情形,提取更准确的故障特征值。应用RS的约简特性消除隔爆电动机故障数据中的冗余信息,并采用具有小样本特点的SVM实现数据的快速提取,利用“一对多”的多故障诊断方法对电动机故障进行诊断,实现了故障的快速分类。仿真结果表明,基于SVM与RS的隔爆电动机故障诊断方法简单有效,诊断结果准确率达到92.857 1%。

参考文献:

[1] 马祥军.矿井电机故障分析及预防[J].中国设备工程,2008(10):23-24.

[2] 杨存祥,朱琛,解豪杰.基于RPROP神经网络算法的异步电动机故障诊断[J].电力自动化设备,2012,32(1):80-83.

[3] 孙丽玲,许伯强,李志远.基于旋转不变信号参数估计技术与模式搜索算法的异步电动机转子故障检测新方法[J].机械工程学报,2012,48(13):89-95.

[4] 阳同光,蒋新华,付强.瞬时功率频谱分析在牵引电机转子故障诊断的应用研究[J].电机与控制学报,2012,16(10):95-99.

[5] 杜占峰.基于小波包分析的电动机转子断条故障检测[D].北京:北方工业大学,2011.

[6] 郑玉芳,马松龄.基于电流分析法的电动机故障诊断虚拟仪器系统的研制[J].计算机测量与控制,2010,18(3):512-514.

[7] 杨存祥,娄艺,罗丹雨,等.基于定子电流的笼型异步电动机转子断条故障检测[J].微电机,2008,41(9):49-51.

[8] 张敬春,谷爱昱,莫慧芳.基于小波能量谱分析的电机故障诊断[J].电力系统及其自动化学报,2006,18(3):55-58.

[9] 于志伟,苏宝库,曾鸣.小波包分析技术在大型电机转子故障诊断系统中的应用[J].中国电机工程学报,2005,25(22):158-162.

[10] 李川,肖蕙蕙.小波包算法在异步电机故障检测中的应用[J].重庆工学院学报(自然科学版),2007,21(7):107-110.

[11] 侯媛彬.基于粗糙集理论的运输机故障诊断与分类[J].长安大学学报(自然科学版),2004,24(1):104-107.

[12] 臧宏志,俞晓冬.基于粗糙集与支持向量机的变压器故障诊断[J].变压器,2008,45(8):61-65.

[13] 王旭红,何怡刚.基于小波包分析和支持向量机的异步电机转子断条故障诊断[J].微特电机,2010,38(5):34-36.

[14] 张柯,陆剑.小波包分析和最小二乘支持向量机的电机故障诊断[J].微型电脑应用,2015,31(6):6-8.

[15] 林茂,李孝全,张兴.小波包分解法与粗糙集神经网络在电机故障诊断中的应用[J].大电机技术,2012(5): 31-34.

Fault diagnosis of explosion proof motor based on SVM and RS

MA Xianmin1, ZHANG Xing1, ZHANG Yongqiang1,2

(1.School of Electrical and Control Engineering, Xi'an University of Science and Technology, Xi'an 710054, China; 2.Branch of Mining Machinery Manufacturing Maintenance Co., Ltd., Shenhua Ningxia Coal Industry Group, Yinchuan 750001, China)

Abstract:In view of problems that fault data acquisition of explosion proof motor in underground coal mine was difficult and fault data was clutter and nonlinear, a fault diagnosis method of explosion proof motor based on SVM and RS was proposed. Wavelet packet is used to analyze instantaneous power of stator of explosion proof motor, and extract fault feature. Redundant data of fault feature is eliminated using reduction feature of rough set, and the reduced feature is used as input sample of support vector machine, to realize rotor fault diagnosis and classification of explosion proof motor. The simulation results show that fault diagnosis accuracy of the proposed method reaches 92.857 1%.

Key words:explosion proof motor; fault diagnosis; SVM; RS; instantaneous power; wavelet packet

文章编号:1671-251X(2017)02-0035-06

DOI:10.13272/j.issn.1671-251x.2017.02.008

收稿日期:2016-08-30;

修回日期:2017-01-05;责任编辑:胡娴。

基金项目:国家自然科学基金项目(51277149);陕西省教育厅科研计划项目(14JK1472)。

作者简介:马宪民(1954-),男,陕西杨凌人,教授,博士研究生导师,主要研究方向为智能仪器与检测,E-mail:maxianminphd@sina.com。通信作者:张兴(1991-),女,河北衡水人,硕士研究生,研究方向为检测技术与自动化装置,E-mail:1515228151@qq.com。

中图分类号:TD614

文献标志码:A

网络出版:时间:2017-01-22 10:30

网络出版地址:http://www.cnki.net/kcms/detail/32.1627.TP.20170122.1030.008.html

马宪民,张兴,张永强.基于支持向量机与粗糙集的隔爆电动机故障诊断[J].工矿自动化,2017,43(2):35-40.