基于电机电流分析的万能式断路器机械故障诊断*

摘 要：电机电流信号常用于分析电动机本身的故障问题，但对其应用于与电机相连机构的故障分析的研究较少。提出一种基于储能电机电流分析的万能式断路器操作机构故障诊断方法。首先采用Hilbert幅值解调法和改进的小波包阈值法相结合获取交流电流信号的包络线，以解决随机噪声干扰造成的所提取包络线粗糙的问题；然后通过包络线提取电流信号的时间量、电流量以及峭度作为不同故障状态电流波形的特征参数；最后融合模糊聚类和量子粒子群优化的相关向量机实现对断路器正常状态、传动齿轮卡涩、储能弹簧卡涩以及脱落的4种状态的辨识。构建了基于电流分析的万能式断路器故障诊断系统，在不同工况下进行了验证，结果表明该方法能有效提取操作机构储能相关部件的故障特征，实现了对操作机构储能相关部件的故障诊断。

关键词：万能式断路器；储能电机电流；量子粒子群；模糊聚类；相关向量机

1 引 言

万能式断路器是低压大电流断路器，广泛应用在低压配电系统中总进线处或重要设备处。操作机构作为断路器重要的组成部分，其动态特性是决定断路器动作可靠性的重要因素[1-3]。操作机构一旦发生故障，轻则储能不畅、分合闸速度下降，导致断路器动作性能下降；重则出现误动或者拒动故障，以至于造成严重事故。因此，有必要针对万能式断路器操作机构的运行状态进行监测以及故障诊断。

当前，大部分研究人员主要针对高压断路器故障诊断技术进行研究[2]，并且由于振动与声音检测方法无须与断路器电气连接，可实现对断路器状态的非侵入式监测，所以基于振动与声音的监测是断路器机械状态的主流监测手段[3-7]。虽然振声诊断方法的检测效果较为理想，但其弊端也较显著。振动信号对传感器的安装位置很敏感，且易受机械谐振的影响；声音信号在实际的工作环境中，易受工业噪声的影响，诊断效果不佳；为了克服振声监测的这些缺点，将储能电机电流信号监测作为断路器机械故障的无创诊断方法，其优势在于电流信号容易获取，不会对断路器的正常运行产生干扰，而且摆脱了只能从机械角度进行断路器故障诊断的束缚。

电机电流信号分析(motor current signal analysis, MCSA)法是一种利用电机定子电流信号的特征来分析相应故障的方法，已成功应用于检测电机本体故障，如对电机定子、转子和气隙磁通不对称等故障进行诊断[8-11]。虽然MCSA法在检测电机所驱动的设备上不及电机本体故障方面的应用广泛，但由于传动系统故障时产生的波动转矩可很好地反映在电机定子电流中[12-13]，该检测方法逐渐成为研究热点。文献[2]推导并验证了GW6H-252 型隔离开关的分合闸电机电流与操作机构输出力矩之间存在二次函数关系，测试验证了隔离开关机械故障情况下的电机电流信号存在畸变现象；文献[14]将电机电流信号的调制信号双谱波峰的变化作为评判故障状态的依据，实现对电机所驱动的齿轮箱故障的监测；文献[15]通过分析主轴电机电流信号的时频特性来提取相关特征，实现对机床切削加工中刀齿失效的判断。这些以电机电流作为信号来源的检测方法所提取的特征主要为时频域的特征，而对于故障识别的方法，随着人工智能的发展，常采用神经网络(neural network, NN)[5]、支持向量机(support vector machine, SVM)[4,7]等。NN具有一定的抗噪声和泛化能力，但是分类器训练需要较多样本，且存在局部收敛问题。近年来，SVM为解决小样本、高维数、非线性等问题提供了新的思路，得到广泛应用[3-4]，但该算法有规则化系数确定困难、预测结果不具有统计意义、核函数受Mercer条件限制等固有局限。相关向量机(relevance vector machine, RVM)可以有效弥补上述缺陷，还能以概率形式输出结果，更加具备实用性，采用RVM作为模式分类器，可有效提高故障识别率[16]。虽然RVM仅需选择核函数参数，但其分类结果对核参数依然敏感，并且根据经验知识确定的二叉树拓扑结构没有反映出实际故障样本的特性，不利于故障类别的正确识别。

针对上述问题，本文将MCSA法进一步推广应用到低压万能式断路器操作机构故障诊断，同时优化故障诊断模型。对于储能电机电流的交流形式，在降噪的电流信号基础上，利用Hilbert幅值解调法提取信号的包络线，并采用改进的小波包阈值法对包络线降噪，从而精确获取平滑的包络线，以提取电流量、时间节点以及峭度特征，并采用融合模糊聚类和量子粒子群优化(quantum particle swarm optimization, QPSO)的相关向量机构建完全二叉树分类器，实现对万能式断路器操作机构故障的分类。

2 电流信号预处理

2.1 储能电流分析

在电动操作机构中，电机通过传动齿轮拉伸储能弹簧做功，储存分合闸所需能量。储能电机的动态方程为[17]：

(1)

(2)

式中：C为电机常数，具体为，ua为电枢电压，Ra为电枢电阻，ia为电枢绕组中的电流，La为电枢电感，w为转子角速度，J为电机转子及负载的转动惯量，TL为负载力矩，Rw为旋转阻力系数，p为电机的极对数，N为电枢绕组有效导体数，Φ为每极磁通，a为支路对数。

由式(1)和(2)的电压和转矩方程分析，储能电机的电流信号不仅与自身参数有关，负载特性的改变也能反映到电流信号中。以DW15-1600万能式断路器为研究对象，其储能电机为磁阻同步电动机，额定电压为AC380 V，起动转矩不小于8 N·m ，转速为750 r/min。采集其电流的故障检测实验系统如图1所示。

实验系统中工控机通过PCL720+发出断路器动作指令，PCL-818HD完成数据采集功能，其中数据采集卡的采样频率设置为20 kHz，而工控机中的测试软件采用LabVIEW与MATLAB混合编程实现，解决复杂数据的在线处理问题。考虑到工业现场的用电环境，由霍尔电流传感器所测得的储能电机电流信号，需经中值滤波算法平滑脉冲噪声，并结合集合经验模态分解(ensemble empirical mode decomposition, EEMD)阈值去噪法剔除掉电流信号中存在的白噪声干扰[18]，即按照通用阈值函数对模态分量作软阈值处理。图2所示为正常状态下所采集电机电流信号的去噪效果。

由图2的放大窗口可以更好地对比去噪效果，可看出中值滤波和EEMD阈值去噪的联合保证了实际工程中噪声的滤除要求。操作机构正常与故障状态下，消除噪声后的储能电流波形如图3所示。

通过分析图3所示电机电流信号，可知，断路器操作机构正常状态下，储能电机的瞬时启动电流为5.654 A，电流从电机刚启动时的冲击状态到回落的启动时间约为0.2 s，断路器的储能动作历时为1.3 s，储能弹簧到位后，

限位开关动作，切断储能电机电源。断路器每次储能的持续时间应保持相对稳定，若时间特性异常，则表明断路器储能异常。断路器储能弹簧卡涩时的储能持续时间虽改变不明显，但电机平稳运行时电流信号的幅值有所上升；而储能弹簧脱落时，电机的启动电流虽变化较小，但运行总时间和平稳运行幅值有所下降；传动齿轮卡涩状态下的电流信号不仅在电机平稳运行的后期幅值有所上升，而且储能的持续时间有明显的变化。

2.2 电流信号包络线的提取

储能过程中，随着储能弹簧的拉伸，负载增大，在电机平稳运行的过程中，电流信号出现了多个峰值。为更好地提取电流信号的特征，需对信号波形的包络线进行提取。常用的包络线提取方法有检波滤波法、全波整流法、复调制法、希尔伯特变换法。在实际应用中，包络提取算法一般需要与滤波器配合使用。由于Hilbert幅值解调法简单有效且多用于工程应用，故文中选用该方法提取包络。

Hilbert幅值解调法的原理简述如下。

设储能电流信号为x(t)，其希尔伯特变换可定义为[19]：

⊗

(3)

式中：f为与时间t区别的时间位移。以x(t)为实部，(t)为虚部构造一个解析信号：

(4)

(5)

式中(t)/x(t)}为信号相位，而即为电流信号x(t)的幅值包络。由Hilbert变换所提的包络线可更好保留原始信号的局部特性。

以断路器正常状态下的储能电流为例，其Hilbert包络线如图4所示。

2.3 包络分析中噪声的剔除

由于电流信号易受白噪声的干扰，Hilbert幅值解调法所提取的信号包络很粗糙[19]。含噪信号经小波包变换后，噪声能量散布于整个小波域，而有用信号能量集中于一些大的小波分解系数。因此，可对小波系数阈值量化处理，重构信号即可达到有用信号和噪声分离的目的。

具体的基于改进阈值的小波包降噪步骤如下：

1)选取合适的小波基和分解层数对含噪电流信号进行小波包分解；

2)选取合适的变化函数对小波分解系数进行阈值量化处理，得到相应的小波系数；

3)将处理后的小波系数反变换，重构去噪后的信号。

其中，步骤2)的阈值量化处理过程中，利用Donoho等提出的通用阈值公式获取阈值λ后，为克服直接选用软、硬阈值函数所带来的伪吉布斯、边缘模糊等失真现象[20]，采用软硬阈值的折中方法，阈值函数为[21]：

(6)

式中，j,i为经阈值函数处理后的小波系数,Wj,i为处理之前的小波系数。β为调节系数，取值范围为0≤β≤1，当β为0和1时，式(6)分别为软、硬阈值函数。

同样以断路器正常状态为例，选取小波基db25(Daubechies小波)、分解层数为3以及调节系数β为0.3，利用改进的小波包阈值法去除包络分析中的噪声，经阈值处理的小波系数重构结果如图5所示。

由图5可见，采用改进的小波包阈值算法，可有效滤除噪声干扰，获得清晰的包络轮廓，提高了包络提取的精度。

3 特征提取及有效性分析

以万能式断路器DW15-1600为研究对象，人为进行故障模拟。采用硬铁丝固定不同长度的储能弹簧，将木楔置于传动齿轮之间，从而模拟了断路器不同程度的储能弹簧锈蚀与传动齿轮卡涩状态，以模拟不同工况下的断路器机械故障。由于弹簧脱落发生概率较小且较为极端，故仅模拟了一根弹簧脱落情况。

建议中国石油企业选择君主制国家的油气项目开展合作。阿联酋、卡塔尔在风险防控上相对容易，可作为投资首选。阿曼资源有限，且现有油田均掌握在西方公司手中，需要等待时机。巴林的海上非常规资源勘探取得突破，是该地区潜在的投资机会。

在断路器的储能过程中，当储能弹簧与传动齿轮发生卡涩故障时，电机电流的幅值和时间会发生相应的变化。因此，本文通过电流信号的包络线来提取电流的特征参数：时间量和电流量，即提取电机的启动电流，电机平稳运行时电流包络的最值及其所对应的时刻，电机运行总时间以及包络线的峭度。

其中，电流包络的峭度为描述包络线尖峰度的参数，对电流信号中的冲击成分十分敏感，能更好地反映不同工况下操作机构的故障情况。

(7)

式中，x为电流信号的瞬时幅值，为电流均值，p(x)为概率密度，σ为信号标准差。

针对操作机构的传动齿轮和储能弹簧这两个关键部件，主要模拟试验了操作机构的卡涩故障。图6所示为正常、传动齿轮卡涩以及弹簧卡涩状态及脱落下典型信号的时间与幅值特征。

本文利用断路器储能电流包络的幅值、时间节点以及峭度信息作为分析操作机构故障的特征量，来构造断路器操作机构故障诊断模型所需向量序列，表1所示为部分实验数据的特征序列，分别将正常、传动齿轮卡涩、储能弹簧卡涩、弹簧脱落状态记为类1～4，其中I1表示电机启动时电流最大值，t1表示电机运行的总时间，I2、I3、t2、t3分别表示电机平稳运行时的最小值和最大值及其对应的时刻，K为包络线峭度。

表1 部分实验数据特征
Table 1 Characteristics of partial experimental data

3.2 特征量提取的有效性分析

对于万能式断路器储能电流特征量提取的有效性分析是开展断路器操作机构故障诊断的基础。而模糊聚类作为非监督学习方法，其主要思想就是将特征集中的每一个特征向量，按照同一类特征向量所具有的相似性分隔成不同的类。因此，通过对正常状态、传动齿轮卡涩、储能弹簧卡涩及弹簧脱落等4种状态下万能式断路器储能电流特征量的模糊聚类识别，可验证所提取特征的有效性。

模糊聚类的基本思路是将每一个特征向量按照一定的隶属度，划分为某一聚类的中心[22]。设X={x1,x2,…,xn}为特征集，特征集中特征向量的长度为p，n为特征向量个数。将X分为c类，则各特征向量将以某一程度隶属于不同的c个类别，以uij表示X的第j个特征向量属于第i类的隶属度，得到Uc×n隶属度矩阵。其中，Uc×n={uij}，uij∈[0,1]，且uij还需满足一下约束条件：

(8)

进而有式(9)目标函数的成立。

(9)

式中:V为聚类中心，vi为第i类的聚类中心，m为模糊加权指数，常取m=2，dij(xj,vi)表示特征向量xi到聚类中心vi的欧式距离。利用迭代法不断修正U和V中的值，最终得到最优聚类中心V={vi}和隶属度矩阵Uc×n={uij}。

因此，可将表1所示的每行特征序列作为特征向量，即形成12×7的特征向量矩阵X，其中X中特征向量的顺序与表中的特征序列一致。模糊聚类的输入为X，聚类数目c=4，模糊加权指数m=2，迭代中止因子为10-5，最大迭代次数为50。此时得到的4类特征数据的聚类中心V和隶属度矩阵U如式(10)、(11)所示。

(10)

U=

(11)

隶属度矩阵U的行号与所属类别对应，列号与特征向量的顺序编号对应，即U中的第1～12列与表1中的第1～12行的顺序一一对应。U中的每个数据表示特征向量对应于所在行号类别的隶属度，故隶属度矩阵每一列中最大数据所在的行即为特征向量所对应的类别。由矩阵U可知，第1～3列的最大值同在第一行，故表1的第1～3个特征数据属于同一类别。同理，第4～6、7～9、10～12列的最大值分别处在不同的行，且同一类的特征向量的隶属度最大值出现在相同的行，与实际的状态类别一致，验证了所提取的特征数据有效，能够正确地区分操作机构不同的机械状态。

4 相关向量机优化及故障识别

4.1 完全二叉树相关向量机分类器的建立

相关向量机是一种基于贝叶斯定理和马尔可夫性质的稀疏概率模型[23]，对于给定数据的特征序列为，类别标签序列为，(i=1,2,…,n)。RVM的分类函数模型定义为：

(12)

式中：x为输入的特征序列集，n是序列个数，ωi是模型对应的权重，K(x,xi)是核函数。将S函数(logistic sigmoid)σ(y)=1/(1+e-y)引入分类模型中，假设)服从贝努利分布，这里的分类问题不包含噪声变量σ2，得到似然估计概率：

(13)

在贝叶斯框架下，通过极大似然法可以获得权值，为了克服过学习的缺陷，这里的ω服从高斯先验分布，即)为每个权值定义了约束参数以实现光滑模型：

(14)

式中：η=(η0,η1,…,ηn)T是n+1维超参数，对每个权值引入超参数，从而产生稀疏概率模型。

单一相关向量机仅能实现数据的二分类，而要解决多分类问题还需建立分类器的二叉树模型。由于二叉树的拓扑结构影响着分类模型的性能，为获得最优完全二叉树，同样可利用模糊聚类来反映各类样本间的相似性和关联性。通过模糊聚类使相似的类聚集在一起，并根据隶属度矩阵来建立二叉树结构，可最大限度地减少误差累积，提高识别准确率。

按照上节中的模糊聚类分析过程，此时仅需将聚类数目c改为2，其余参数均保持不变，可得出所需的隶属度矩阵Utree如式(15)所示。

Utree=

(15)

分析以上隶属度矩阵，可得出特征集X中的1～3和10～12所属的类别更相近，4～9所属的类别更类似，即正常状态与弹簧脱落相近，传动齿轮卡涩与储能弹簧卡涩相近，从而得出图7所示的完全二叉树。

4.2 QPSO优化RVM分类模型

基于RVM多分类模型包含3个分类器，故有3个核函数参数需要设置。由于核参数人为经验的选取易导致诊断结果不精确，所以本文利用量子粒子群优化RVM的训练过程，获得核参数的最优值。

QPSO相较于标准粒子群优化算法，由于QPSO算法中的粒子在量子空间移动时没有确定的轨迹，可自由在可行解空间中寻优，从而保证了更好的全局收敛[24]，QPSO对RVM分类模型的寻优实现步骤如下：

1)确定RVM多分类模型的拓扑结构；

2)在解空间中初始化一组粒子，其中种群规模为15，粒子的维数为3，最大迭代次数为50，并赋其为个体最好位置；

3)根据适应度函数即K-折交叉验证识别率，其中折数K设置为10，计算粒子当前位置的适应度值，并根据适应度值更新粒子的个体最好位置和全局最好位置；

4)根据粒子的进化方程计算得到粒子的新位置[25]；

5)迭代更新满足结束条件时退出，此时的粒子位置即为RVM的最优核参数。

基于上述研究方案，本文在万能式断路器DW15-1600上进行了验证。万能式断路器的每种机械状态分别做30组，共120组实验并提取电流包络的特征量。为避免训练数据的偏斜问题，随机选取每种机械状态中的20组用于分类模型的训练及参数优化。此时，RVM1、RVM2、RVM3的核参数依次为5.64、2.21、3.57。采用剩余的样本，即剩余40组特征用于分类器的测试，分类情况如表2所示。

表2 部分实验数据特征
Table 2 Characteristics of partial experimental data

测试结果中分类正确的样本有37个，总体的正确率为92.5%，可有效实现对万能式断路器操作机构的故障诊断。其中，易分错的故障类型为弹簧卡涩故障，其原因在于弹簧卡涩受到不同工况的影响，从而出现储能弹簧的轻微卡涩更倾向于正常状态，而所模拟的储能弹簧严重卡涩更倾向于齿轮卡涩造成的故障状况，但整体上识别模型能够对常见工况下的弹簧卡涩 故障正确分类。

4.3 不同识别方法诊断结果比较

为评估本文故障诊断模型的识别效果，分别在不同的二叉树结构、核参数设定条件下进行测试。其中，首先按照图7的二叉树依次构建表3所示的两种节点训练样本不同的完全二叉树层次结构；其次，构建了一种次序二叉树，以更好对比说明。

表3 二叉树的层次结构
Table 3 Hierarchy construction of binary tree

利用表3所建立的不同二叉树拓扑结构，分别测试基于不同二叉树的RVM以及融合QPSO的RVM识别模型，并与QPSO优化的SVM识别方法进行比较。为评估诊断模型的有效性，训练样本以及测试样本均采用同一数据集，同一特征参数进行诊断，其中每个节点的RVM人为经验设定的核参数δ=2.0，各诊断模型的诊断结果如表4所示。

表4 不同诊断方法的正确率对比
Table 4 Accuracy rate comparison of different diagnosis methods (%)

由表4可以得出，整体上基于RVM的诊断模型的识别效果优于基于SVM的诊断模型。具体分析在相同的模型条件下，二叉树层次结构的不同影响着诊断结果，即完全二叉树的识别模型整体上优于次序二叉树模型的诊断结果，而由各类的聚类中心采用分裂层次法生成的完全二叉树模型识别率最好。原因在于：次序二叉树的树型结构深度最大，二叉树根节点处的分类错误会导致整个分类器较高的错误积累；而模糊聚类算法使相似的类聚集在一起，由此所建立的完全二叉树结构在诊断时保证了上层节点分类器诊断的准确性，减少了误差累积，从而提高诊断的准确率。而在相同的二叉树结构下，引入QPSO优化的RVM算法优于单一人为经验设定参数的RVM算法的诊断效果，QPSO能很好地解决全局优化问题并且收敛速度快，可以有效地搜索全局最优核参数值，保证了诊断模型的准确性与可靠性。

5 结 论

储能电机的电流信号中蕴含了操作机构的机械状态信息，利用电流信号的包络线特征能有效表征万能式断路器的故障状态。提出一种基于储能电机电流分析的低压断路器操作机构故障诊断新方法。利用中值滤波与EEMD阈值去噪法相结合，有效去除工业现场的噪声干扰，利用模糊聚类算法，验证了电流包络特征量的有效性，并为后续分类器的二叉树建立提供依据。对不同工况下的断路器故障进行了实验，故障识别率较为理想，验证了该方法的有效性。该模型具有较好的实用价值，为实现断路器操作机构非侵入式的故障诊断提供了一种新的参考。