往复式压缩机拉缸故障多参数统计监测诊断方法

本资料包含pdf文件1个，下载需要1积分

石化行业属于易燃、易爆及有毒物质开发与利用的高风险行业，企业生产属于典型的流程工业。往复式压缩机是石化企业中常见的关键设备，虽然其设计制造水平已经趋于成熟，但是由于其易损零件多、结构复杂，而被业内公认为可靠性较差的设备。

往复式压缩机常见的故障类型有：气阀漏气或损坏、活塞环及支撑环磨损、活塞杆沉降、十字头磨损或松动、轴承损坏等。由于其易损件多、结构复杂等-系列收稿日期：2012-07-19基金项 目：国家重点基础研究发展计划(973”计划)资助项 目(2012CB026000)作者简介：党 露(1988-)，女，陕西富平人，主要从事信号处理与故障诊断方面的研究通信联系人：江志农，男，研究员，博士生导师.E-mail：lunadl###163.con· 10 · 机 电 工 程 第30卷特点，对往复式压缩机进行在线监测及故障诊断尤为重要。而这些功能的实现需要针对其振动、位移、温度等多种类型信号，通过安装在不同测点、不同类型的传感器(例如加速度传感器、压力传感器、温度传感器及电涡流传感器)，以获取往复式压缩机的运行状态信息n。]。

化工 、石油行业的统计数据表明：目前在国内化工行业企业中在运行的100 kW以上的往复式压缩机已经超过20 000台，1000 kW以上的大型往复式压缩机估计在6 000台左右；近20年来，化工行业发生往复式压缩机重大事故平均每年32起，约占行业内重大事故总量的25%。鉴于往复式压缩机具有应用量大、面广、结构复杂、易损件多等特点，对往复式压缩机的故障诊断技术研究-直以来都得到了国内外学者的广泛关注。在国内，有些专家通过对往复式压缩机缸体振动信号进行分析，从而对压缩机工况进行判断。在国外，有些专家借助或依靠先进的传感器技术和计算机信号处理技术判断设备运行状态。获取反映往复式压缩机工作状态的有效特征参数对其故障诊断起着举足轻重的作用。凭借经验或设想去试凑特征参数，再对其进行实验验证”的方法距离实际应用还存在着较大的差距。

近年来，人工智能技术和专家系统被越来越多地应用于工程实际中，文献[4]将PCA方法应用于 自适应过程监控 ，文献[5]讨论了采用PCA实现间歇过程的故障诊断。

本研究将主成分分析法(PCA)应用于往复式压缩机特征参数的选择中，以实现对压缩机拉缸故障的早期预警。

1 拉缸故障及其严重性拉缸故障是往复压缩机中的-类严重故障。当活塞环、导向环异常磨损或断裂后，活塞体与气缸内壁或者气缸衬套严重磨损，导致衬套表面甚至气缸内壁出现严重拉痕，即为拉缸。拉缸故障-旦发生，轻则造成活塞部件、缸体衬套严重磨损 ，重则导致活塞杆断裂、撞缸 、易燃易爆气体外泄，引起着火、爆炸等恶性事故。拉缸故障发生后需要更换活塞部件、气缸衬套甚至气缸体，经济损失通常在数十万元至上百万元 。而-旦造成往复压缩机爆炸等恶性事故，经济损失和社会危害性则难以估量 。

2 拉缸故障的检测方法往复式压缩机拉缸故障的状态监测手段通常包括以下两种：(1)在填料函外安装电涡流传感器，检测活塞杆沉降量；(2)在缸体上安装加速度传感器，检测缸体的冲击振动。

本研究采用第2种方法，即基于加速度传感器信号的方法。理由是安装方便、测量准确，测量值为缸体相对振动加速度。相比较而言，电涡流传感器安装则较为复杂，初始安装间隙调节非常重要，测量值为相对振动位移。除此之外，加速度传感器的价格较低，经济性好 I。

3 拉缸故障早期特征提取方法的研究对于这些关键设备的在线监测和故障诊断，振动是最重要的信息来源及诊断依据，蕴含了丰富的机械设备异常或故障信息，直接反映机械设备运行状态。

通过对装置测点的振动信号进行信号采集及处理分析，本研究对机械设备或结构的运行状态和可能出现的故障提出预测、估计、判断，并为设备的更换和修复乃至生产过程提供决策依据。

在设备的实际运行过程中，计算机监控系统记录和存储了大量的生产过程数据，但是由于设备故障早期的特征信息很微弱，信噪比很小，有用的特征信号往往淹没在其他部件运行中引起的振动信号和大量的随机振动信号中。因此，如何提取故障早期的特征信息已成为工程中实际面临的重要问题。

3.1 特征参数的计算特征参数在机械设备的状态检测及故障诊断中起着举足轻重的作用。若研究者所选取的特征参数能够准确、灵敏地反映设备运行的状态特征 ，那么对机械故障的自动诊断就成为可能n ]。

常用的特征参数有：(1)有效值Xm幅(RMS Value)。有效值是振动信号振幅的均方根值，定义为： 式中： -采集信号，Ⅳ-采样点数。

有效值 也称均方根值 ，用来反映信号的能量大小，特别适用于具有随机性质的振动测量。有效值反映-个波形的整体总能量，但无法反映短时脉冲振动波形的幅值。

(2)峰值 X (Peak Value)。其定义为： ∑ (2)～ ： 1通常计算峰值的方法是把 的 个采样点分成若干段，在每-段中找出-个绝对值最大的采样点作为该段峰值，然后对所得峰值取平均数。本研究将采第 1期 党 露，等：往复式压缩机拉缸故障多参数统计监测诊断方法 ·11·样数据均分为l0段，求出每段峰值，计算式为：X.Amax-Ami.(p1，2，..，10P ) (3) ' J j/再将算得的10个峰值取平均数。

峰值可以反映设备某-局部故障点的冲击力大校冲击力越大，峰值越高，在冲击性振动方面，峰值比有效值更有效。

(3)波形指标。有效值与均值之比称为波形指标，其定义式为：SX /x (4)式中： -信号的均值。

(4)脉冲指标。峰值与均值之比称为脉冲指标，其定义式为：Ix x 5、)(5) 因子(K Factor)。峰值与有效值的乘积称为 因子，其定义式为：KX ×X (6)由上式可知 K因子是有量纲的参数，随信号的峰值和有效值的变化而变化。

(6)峰值因子C(Crest Factor)。峰值与有效值之比称为峰值因子，无量纲参数，其定义式为：C ok/X (7)峰值因子 C不受振动信号绝对幅值的影响，不随传感器、放大器等仪表灵敏度变化而变化 ，-旦设备出现损伤，产生冲击信号，峰值因子明显增大，但是当故障不断扩展，峰值因子又逐渐减小，直至恢复到无故障水平。因此，研究者根据峰值因子，可进行有效的早期预警，并预测故障的发展趋势。

(7)峭度 (Kurtosis)。峭度离散化的公式为：善( ) (8)峭度 K 是无量纲参数 ，反应了振动信号分布特性。峭度的绝对值越大，说明设备偏离其正常状态，故障越严重。

(9) - 7，- 歪度 5 是无量纲参数，对轴承的转速及负载变化不敏感。与峭度相比，其值变化范围较校由于不同的特征参数对于不同的机械故障的灵敏度各不相同，如何选择能够明显反映故障特征的特征参数对有效反映机械设备的运行状态非常关键。

本研究采用了主成分分析法(PCA)，对多种特征参数进行分析，并找出关键的特征参数。

3.2 主成分分析(PCA)算法本研究利用主成分分析(PCA)算法，对原始数据进行降维处理。具体算法如下：(1)对原始 数据分 4段滤波 ，并求 其有效值RMS1、RMS2、砌 和RMS4。

(3)共得到8个特征参数X ( l，2，..，8)，特征参数 置维度为：NX Xi .，X ， ( l，2.，N)(4)建立样本矩阵 ：M ] (10)(5)求样本矩阵肼 的协方差矩阵： M×M (11)(6)求协方差矩阵 的特征值及相应的特征向量。

(7)将协方差矩阵 肘 的特征值按照其绝对值大续行降序排列，同时调整其对应的特征向量的序列 ，得到矩阵 P。

(8)计算得到特征矩阵：KP×M (12)(9)特征矩阵 的前两行即为对多特征参数处理后最能代表故障征兆的2个特征参数 。

在实际运用中，本研究首先对往复式压缩机正常运行下所采集到的大量数据进行上述运算处理，得到正常工况下的值域；设定合理的阈值后，再对未知工况下所采集到的数据进行实时处理；超出阈值后报警，从而实现了对往复式压缩机拉缸故障的早期预警。

4 实际数据验证本研究所选用的数据来源于中石油某石化企业2010年发生的-起拉缸故障(往复式压缩机活塞与气缸内壁磨损图如图1所示 )，该机组安装有 BH5000R往复压缩机在线监测系统 ，机组已安装检测十字头冲图1 往复式压缩机活塞与气缸内壁磨损图机 电 工 程 第30卷击的加速度传感器。

本研究将 8个特征参数两两组合 ，分别作为轴、y轴的值，得到XY图(如图2～7所示)。

本研究将数据进行PCA处理后，提取影响系数最大的两组数据作为新得到的特征参数，并反映在 XY图中(如图7所示)。

将特征参数两两组合的效果图进行比较和总结，得到结果如表 1所示。

： 嘲 - 18重305045柏 353Ou 252015105由表 1可知：(1)有效值与峭度、歪度相组合的效果较好；(2)峭度、歪度、峰值因子两两组合，效果都不理想，不能判断是否有早期故障；(3) 因子与其他特征值组合的效果普遍较好。

PCA处理后的结果图如图8所示，该图可明显区分正常工况和早期拉缸故障。

以上比较分析结果说明，对数据进行主成分分析L I勘 盛 L。

f- UP0 2 5 5 7.5 l0 12.5 l5 17 5 20 22.5 25 0 2.5 5 7.5 10 12.5 15 17.5 20 22.5 25RMS1/(m s 1 q 27 52522 52017 515邑12 5107.552.5 1/(m s --403020lO时 0- 10- 2000- l0 2·5 5 7·5 10 12-5 15 17·5 20 22·5 25RMS2/(m s，l - l/(m )t- r 0 2.5 5 7 5 l0 l2 5 15 17.5 20 22.5 25RM'Sl/(m s，窖 - 0 2.5 5 7 5 10 12.5 15 l7.5 20 22.5 2540302010时 0- 1O20- 30RMS1/(m·fIRMSl/(m )图2 特征参数组合图之-导RMS1/fn10 2.5 5 7.5 10 12.5 1 5 17·5 20 22·5 25RMS2/(ml -- - r27 52522 520l7 5I5量12.5107 552.50RMS2/(ml-啊0 2.5 5 7.5 10 12 5 l5 17.5 20 22.5 25RMS2/(m s). III0 2 5 5 7.5 10 12.5 15 17 5 20 22.5 25 0P-M S2/(m-5 7.5 10 l2.5RMS3/(m s )- 0 2.5 5 7.5 10 12.5 15 17·5 20 22-5 25 ORMS 2，(m )图3 特征参数组合图之二5 7.5 10 l2.5且 3/(m s图4 特征参数组合图之三， 量 ∞i÷ ～如 蚰 ∞ ∞ ：2 mu m 9 8 7 6 5 4 3 2 0 ∞- 日lJ，ng 奄 如 加 m 5 O柏 ∞ 加 ：2 m 0 堪 j m 6 4 2 0- 旦寻 ， 由毪譬8 6 4 2 O 8 6 4 2 0 08耋 8 6 4 2 O 8 6 4 2 O 由氇 釜第1期 党 露，等：往复式压缩机拉缸故障多参数统计监测诊断方法柏 30201O时 0- l0-20.30.40且 1/(m·027·52522.52017 515量12.5107 552.5054 543.53U2.521.5105lL l rP-. 4/(m s- ，IPk - O 2.5 5 7 5 10RMS3/(m·lI柏 3020lO时 0· lO- 20· 30枷 0 2 3 4 5RMS'4，'(m ：IIr0 2 3 4 5 6RMS钏m 。

图5 特征参数组合图之四3020lO时 0- 1O-20-30·曙壤-，.- 蔓： 。

j I 霹 龟 ：- 毋 r IIO 2 4 6 8 10 12 14 16 18 20Kurtoxv'0 2 4 6 8 10 12 J4 16 l8 20Kuno。

图6 特征参数组合图之五5045403530ru25201Sl05. I t·- - ： .。兰粥 酆 -1 f.f - F O。。 J r瓢1 lJ -40 -30 -20 -10 0 10 20 30 405% wB譬l ，27.52522 52017 5蛩 1512.5l07 52.505.O5%mqr 0 10 20 3O 40 50CttFactor图7 特征参数组合图之六蘸警 lh 1，-30.0 -200 -10-0 0.0 10.0 200 3n0， 8 A l图8经PCA处理过的结果r(， 目挪 ： 2 。 宝 。

T. 星 N a星 8 6 4 2 O 8 6 4 2 O 喜g遗 揶 ∞ m" ：3。

r 暑· 14 · 机 电 工 程 第30卷(PCA)够有效反映设备的运行状态，并能很好地实现拉缸故障的早期预警，同时解决了对特征参数的选择问题。

5 结束语笔者研究了对往复式压缩机拉缸故障的早期预警方法，得到了如下结果：(1)不同的特征参数表征着不同的故障类型，对于现场采集到的原始数据，选取适当的特征参数显得尤为重要；虽然有些特征参数两两组合也能很好地区分往复式压缩机的设备工况，但是也存在个别不能区分设备工况的情况 ，总体状态不稳定，易给压缩机早期预警埋下隐患。

(2)经PCA处理之后，针对压缩机预警状态稳定，能很好地区分设备工况。工程实际应用过程中，可通过对压缩机缸体加速度信号的提取，实现对往复式压缩机拉缸故障的早期预警。