基于压缩感知的密封端面声发射信号提取技术研究.pdf

年 月第 卷 第 期润滑与密封 ：文献引用：张尔卿，陈平，周燕基于压缩感知的密封端面声发射信号提取技术研究润滑与密封，（）：，（）：基金项目：国家重点研发计划项目（）收稿日期：；修回日期：作者简介：张尔卿（），男，博士，高级工程师，研究方向为密封摩擦机制及状态监测技术。：。基于压缩感知的密封端面声发射信号提取技术研究张尔卿 陈 平 周 燕（中密控股股份有限公司 四川成都）摘要：机械密封动、静环高速转动摩擦时会产生一定频段的声发射（）信号，根据声发射信号可以判断密封的运行状态。为提高 信号数据传输效率，提出一种基于压缩感知的密封环声发射数据降维方法：利用稀疏矩阵将采集的原始 数据进行稀疏表示，再通过观测矩阵将稀疏表示后的 数据进行压缩采样以达到降维的目的，数据传输完成后利用 等信号重构算法将压缩后的信号恢复。将 原始数据与重构后的数据进行对比，验证了该方法的可靠性。通过实验对比证明，在压缩比为 ，训练样本集 值为 时，选择 作为稀疏矩阵，伯努利随机矩阵作为观测矩阵，采用 算法的信号重构效果最佳，其均方根误差为 。关键词：机械密封；声发射信号；压缩感知；数据传输；信号重构中图分类号：（，）：，（），：；机械密封是一种旋转轴封装置，常用于压缩机、反应釜等装置。其工作原理是通过密封端面摩擦副动、静环的贴合、旋转、补偿来实现阻漏的过程，因此密封端面的接触状态决定着机械密封的性能。机械密封在运转状态时，因端面摩擦会产生一定频段的声频，且不同摩擦副材料、不同运行转速或不同端面状态下，摩擦所产生的声音频率不同，因此可以通过摩擦副声发射信号数据判断摩擦副的运行状态，从而监测机械密封是否失效。声发射（，）检测是一种无损检测方法，其信号属于高频剧变信号，对数据采样频率要求较高，一般采集的声发射数据数量巨大且存在数据冗余，非常不利于数据的传输，需要通过数据降维再恢复的方法提高数据传输效率。压缩感知（，）是一种寻找欠定线性系统稀疏解的技术，由 等在 年提出，常被应用于信号处理中，可进行信号压缩、剔除冗余信息和降维降噪等，也用于重构稀疏和信号恢复等。本文作者将压缩感知技术应用于机械密封端面声发射信号传输过程，利用压缩感知相关算法对原始 信号进行压缩，再通过重构算法将压缩后的 信号进行恢复，对比原始 信号与重构信号的相似度，以探究机械密封端面声发射信号的压缩传输效果。压缩感知原理 压缩感知压缩感知的理论主要分为稀疏表示、信号压缩、信号重构 个部分，其中信号稀疏表示是压缩感知理论成立的前提条件，信号压缩是数据传输过程的关键步骤，信号重构是压缩感知的最终目的。压缩感知算法流程图如图 所示。图 压缩感知算法流程 信号稀疏表示信号的稀疏表示是指经过稀疏变换后的信号中非零元素尽可能减少，如果一个 的信号里，非零元素不超过 个，那么该信号可称为 稀疏信号，因此信号通过稀疏变换后是离散信号。在实际采集信号过程中，往往信号都不是以稀疏的形式存在的，为了便于信号的处理，提高计算效率，则需要对原始信号进行稀疏变换。对于采集的原始信号 ，为 长度的列向量，稀疏变换基 ，为 维矩阵，其中的任一元素 ，为 长度的列向量，则信号 可以表示为（）其中 ，当 中有 个非零元素时，且 时，则 为 稀疏信号，同时 也称为原始信号 的稀疏表示。从公式（）中可以看出，信号的稀疏表示的 值与 密切相关，稀疏变换基 也称为稀疏矩阵或稀疏字典，压缩感知理论在不同领域中研究最广泛的就是如何找到合适的字典。目前稀疏字典主要分为两类，固定字典和学习字典。其中比较典型的固定字典有傅里叶字典、离散余弦字典等。快速傅里叶变化（）的核心就是将时域信号转换为频域信号，因此在频域较少的信号中，对信号的稀疏效果是非常好的。离散余弦变换（）是离散傅里叶变换（）的一种特殊形式，省略了 中的复数运算，通过构造偶函数的方式来实现，相对于，具有较好的频域聚集度，可以将不重要的频域信号去掉。奇异值分解（）是一种比较主流的学习算法，其综合了 算法和奇异值分解（）算法的思想，利用 算法对误差项进行分解，将分解得到误差最小项作为修正系数对字典元素进行更新，经过迭代得到最优解。信号压缩信号稀疏表示后虽然可以简化计算，但信号的维度并没有降低，稀疏系数 仍然是 维数据，对于数据传输效率并没有得到优化。这里引入一个与稀疏矩阵 不相关的观测矩阵，观测矩阵 ，为 维矩阵，其中的任一元素，为 长度的列向量，且要求。则压缩后的信号 可以表示为（）即：（）其中 ，。由式（）可以明显看出，通过观测矩阵 与原始信号 的乘积，得到了压缩后的信号，信号数量从 个降低到了 个，由于，起到了数据降维的作用。当然寻找合适的观测矩阵 会直接影响信号重构的效果。观测矩阵也分为 种类型：确定矩阵和随机矩阵。确定矩阵主要有部分哈达玛矩阵等，随机矩阵主要有高斯随机矩阵、伯努利随机矩阵等。信号重构信号压缩完成后即可将数据进行传输发送至服务终端，最后就是将接收到的数据进行信号重构，信号重构相当于对公式（）中，已知 和 的值，求解线性方程组得到 的值，但由于，使得方程组中未知量大于方程组的个数，因此该解析式是一个欠定方程组。这里可以引入一个感知矩阵，使得：（）即：（）从式（）和式（）可以看出，这里方程将求 年第 期张尔卿等：基于压缩感知的密封端面声发射信号提取技术研究 值转变为求，然后通过间接的方式求得原始信号，因为 是稀疏信号，极大地减少了未知量，优化了问题求解过程。因此将信号重构问题转化为 （）根据式（）求得 值后，代入式（）中，即可重构出信号，将原始信号 与重构信号 进行对比，用均方根误差（）来判定 与 相似度，再根据信号压缩比（）进行综合评判，找到适合机械密封声发射信号传输过程的最优压缩感知算法。均方根误差如式（）所示。（）（）压缩比如式（）所示。（）贪婪算法是一种主流的信号重构方法，其基本思想是利用迭代的方式重建信号的支撑集，通过某种贪婪准则，求得待估信号构成的元素，贪婪算法可以显著降低计算的复杂程度。正交匹配追踪（）是一种常用的贪婪算法，它主要有识别、增强、残差更新 个阶段，通过反复迭代计算信号的支撑集，不断更新残差系数，最终获得重构信号。实验数据采集 实验设备实验设备采用国内某知名企业机械密封运转实验台，机械密封采用双端面机械密封，动、静环材料选用石墨硬质合金组合，润滑介质采用油润滑，声发射探头安装在机械密封外压盖开孔处，要求尽可能贴近动、静环处，机械密封声发射信号采集实验台如图 所示。图 机械密封声发射信号采集实验台 实验参数及过程声发射数据采集前需要预先拟定实验参数。实验设置最大静压为 ，开机转速为 。设置采样频率为 ，每 保存一次信号采集文件。机械密封声发射采集传输实验流程图如图 所示。图 机械密封声发射采集传输实验流程 实验数据分析文中的核心旨在研究机械密封声发射信号传输效率和数据恢复效果，因此实验数据分析最重要的环节在信号压缩和信号重构 个阶段。信号压缩阶段信号压缩的关键是寻找合适的稀疏矩阵 和观测矩阵，稀疏矩阵 会影响稀疏系数 的稀疏程度，从而间接影响信号重构的效果。稀疏矩阵 和观测矩阵 需要保证 有点等距性质，性质的等价条件是保持稀疏矩阵 和观测矩阵 不相关性。提取实验采集时长 的原信号数据，其声发射信号的长度为，如图 所示。图 声发射原始信号波形 在实现压缩感知算法过程中，信号的稀疏性对重构的效果至关重要，波形从时域到频域等空间域的转换，可以对实验装置的特有声频突出显示，同时可以排除不存在的声频和淡化周围环境噪声，信号从时域到频域的转化也是对信号稀疏的过程。从图 中明显可以看出，声发射信号是由不同频率的声波叠加而成，这里利用 和 对原始声发射信号进行稀疏表示，效果如图 所示，其中图 （）为快速傅里叶变换后的效果，图 （）为离散余弦变换后的效果。润滑与密封第 卷图 声发射信号稀疏表示效果 ：（）；（）对比图 和图 可以发现，稀疏表示后的信号大部分值都明显趋近于，而声发射频域等空间域信息被明显突出。从图（）、（）中可以看出，不同的稀疏算法，对信号的稀疏分布也不一致，也会直接影响信号重构的效果。通常在实际信号应用中，声发射的频域是非常复杂的，可能会出现大量声频的情况，使得信号在 或 变换后达不到稀疏的效果或稀疏效果很差，使得重构信号的误差很大。而学习字典可以在不同工况环境条件下，自适应地去寻找信号的稀疏表示，得到在有限训练样本条件下，输出最好的算法模型用于信号重构。在压缩感知算法中，不同的压缩比对信号的稀疏系数也会有很大的影响，文中主要研究压缩比为 之间声发射信号的重构效果，以 算法为例，分别选取压缩比 为 和 时的信号稀疏系数效果示意图如图 所示。对比图 （）和图 （），在压缩比不同时，稀疏信号的频率和幅值看不出明显的差别。而对比图（）和图 （）可以发现，训练样本集 值越大时，其信号稀疏程度更大。除此之外，与 和 不同的是，稀疏表示后的信号数量会减少，这与选择的训练样本集 值有关，当训练样本 时，稀疏表示后的信号数量为 个；而当训练样本 时，稀疏表示后的信号数量为 个。图 的声发射信号稀疏信号效果 ：（），；（），；（），；（），压缩比的大小同时也决定着观测矩阵的大小，压缩信号的长度与观测矩阵的行数相等，因此观测矩阵的选择与信号重构效果也有直接关系。文中主要运用部分哈达玛矩阵、高斯随机矩阵和伯努利随机矩阵作为观测矩阵，渐变云图可以很清楚地显示各种矩阵特征情况，分别选取矩阵中 像素大小图像进行 年第 期张尔卿等：基于压缩感知的密封端面声发射信号提取技术研究 观察，其图像效果如图 所示。图 （）为部分哈达玛矩阵，该矩阵是由 组成的正交矩阵。图 （）为高斯随机矩阵，该矩阵数值的随机概率符合正态分布的特点，从图像中也可以清楚地看出大部分的数值接近于。图（）为伯努利随机矩阵，该矩阵也是由 组成，但与部分哈达玛矩阵不同的是，和 的概率均为 且随机分布，从图 （）和图 （）的纹理分布也可明显看出。以上 种观测矩阵对信号进行压缩处理时的运算效率以伯努利随机矩阵最快，其次为部分哈达玛矩阵，最后是高斯随机矩阵，但它们处理数据的速度都在秒级。任意一种稀疏矩阵和这 种观测矩阵组合，以处理 一个通道的声发射数据为例需要的最长时间依次为：伯努利随机矩阵需要 ，部分哈达玛矩阵需要 ，高斯随机矩阵需要 。图 观测矩阵渐变云图效果 ：（）；（）；（）信号重构阶段信号重构的目标是将压缩后的信号恢复成原始信号，其原理是通过重构算法，将压缩信号 转化求解得到稀疏系数，再利用已知的稀疏矩阵 与稀疏系数 相乘，便可求得重构信号。采用不同稀疏矩阵、观测矩阵 以及不同压缩比时，信号重构的效果是不一致的，信号重构效果的好坏可以用均方根误差（）来进行判定。如图 所示。对比图 （）、（）、（）可以看出，采用部分哈达玛矩阵作为观测矩阵重构的均方根误差值整体偏大，高斯随机矩阵和伯努利随机矩阵的重构效果非常接近。由于伯努利随机矩阵是由 和 组成的，相对于高斯随机矩阵而言，降低了计算难度，提高了计算效率，因此在效果相近的情况下优先选择伯努利随机矩阵。图 在不同压缩比下的重构均方根误差 ：（）；（）；（）通过 算法求得的稀疏矩阵计算均方根误差值整体较小，特别在压缩比为 左右时，均方根误差值为（），信号压缩量较大且信号重构效果较好，此处 选取的训练样本 值为，该算法用于重构信号的效果也与选取 值的大小有关。因此针对 算法，选择伯努利随机矩阵，在不同压缩比、不同训练 值情况下，重构均方根误差效果曲面图如图 所示。润滑与密封第 卷图 重构效果 由图 可以看出，随着训练样本和压缩比的增大，其均方根误差值明显减小，信号重构性能提升。但对于文中压缩感知的根本目的是为了解决声发射信号传输问题，若训练样本较多时，会增大数据运算量，降低程序运行效率。而压缩比较大时，会影响信号压缩效果，信号传输效率也无法提高。从图 中可以发现，三维曲面梯度下降最快点在压缩比 、训练样本集数量 的位置，因此在该范围内信号重构性能最佳。这里取 、，稀疏矩阵采用 算法、观测矩阵采用伯努利随机矩阵，重构算法采用 算法，声发射信号原始采集信号与重构信号对比效果如图 所示。从图 中可以看出，重构信号曲线（红色曲线）对原始信号曲线（蓝色曲线）的拟合效果是非常好的，基本还原了信号的波形、轮廓及峰谷形状，均方根误差为 。图 声发射信号重构效果对比 针对机械密封声发射信号，采用不同稀疏矩阵和观测矩阵，基于 算法进行重构，对于不同压缩比，其压缩感知重构性能参数如表 所示。从表 中也可以明显看出，和 的重构效果最好，考虑到伯努利随机矩阵的运算效率是高于高斯随机矩阵的，因此对于声发射信号重构采用 算法是最优的，再考虑到其信号压缩和信号传输的综合性能，因此选择训练样本集 ，压缩比 的重构效果最优，其均方根误差值为（），小于压缩比 时的平均基准值 ，满足工业现场应用基本要求。表 机械密封声发射信号压缩感知重构性能参数 稀疏矩阵观测矩阵均方根误差 注：表示部分哈达玛矩阵，表示高斯随机矩阵，表示伯努利随机矩阵，其中 算法的训练样本集。结论提出一种基于压缩感知的密封环声发射数据降维方法，即利用稀疏矩阵将采集的原始 数据进行稀疏表示，再通过观测矩阵将稀疏表示后的 数据进行压缩采样以达到降维的目的，数据传输完成后利用 等信号重构算法将压缩后的信号恢复。在某些特定工况条件下，对采集原始声发射信号进行不同算法处理，对比不同算法对信号数据压缩和重构效果，探究适用于机械密封声发射信号的最优算法，以应用于信号传输。主要结论如下：年第 期张尔卿等：基于压缩感知的密封端面声发射信号提取技术研究 （）对于稀疏矩阵的研究表明，在频率较少的信号波形中，使用 正交基和 正交基的信号稀疏表示效果较好。但对于复杂的机械密封声发射数据而言，使用 的重构效果明显优于 正交基和 正交基。（）对于观测矩阵的研究表明，使用部分哈达玛矩阵的整体重构误差值较大，高斯随机矩阵和伯努利随机矩阵的重构效果非常接近，但在算法计算量和运行效率上，伯努利随机矩阵是优于高斯随机矩阵的。（）综合考虑信号压缩率和重构效果，当选择训练样本集，压缩比 时，选择 作为稀疏矩阵，伯努利随机矩阵作为观测矩阵，采用 算法的信号重构效果最佳，其均方根误差为（）。参考文献 张尔卿机械密封端面状态监测及寿命预测关键技术研究成都：西南交通大学，：，：尹源，黄伟峰，刘向锋，等机械密封智能化的技术基础和发展趋势机械工程学报，（）：，（）：，：高志，林尤滨，黄伟峰，等干气密封启动过程中的声发射信号特征清华大学学报（自然科学版），（）：，（），（）：，：，：王聪基于压缩感知的机械故障诊断方法研究昆明：昆明理工大学，：，李泽芃，乔百杰，文璧，等基于压缩感知的多级风扇周向声模态重构航空动力学报，（）：，（）：，（）：，：王番，魏超，刘智，等基于 稀疏压缩感知域内遥感图像融合测绘科学技术学报，（）：，（）：，：李帅永，毛维培，程振华，等基于 和 字典学习的供水管道泄漏振动信号压缩感知方法仪器仪表学报，（）：，（）：，：，：李沐阳，苑宇基于压缩感知的滚动轴承故障信号提取技术研究大连交通大学学报，（）：，（）：，：魏子然基于压缩感知的信号重建技术及其成像系统的研究成都：电子科技大学，：，：，：润滑与密封第 卷