专利名称：一种设备振动信号的处理方法
技术领域：
本发明涉及检测设备故障隐患的方法，具体涉及利用设备振动信号检测设备故障隐患的 方法，特别是对测得的设备振动信号进行处理的方法。
背景技术：
为了能提前检测到设备的故障隐患，人们发明了以机械设备振动信号为状态参量的设备 故障隐患检测方法。
目前，利用设备振动信号检测设备故障隐患的方法， 一般是基于平稳过程的传统信号处 理方法，如信号滤波法、时域平均法、自适应消噪法等，这些处理方法的采用使故障的早期 诊断取得了许多进展。但上述处理方法也存在缺陷如信号滤波法一般不能同时得到滤波通 带以外的信息；时域平均法要求信号平稳并严格按周期采样，这在实际应用中受到运行工况 非平稳的限制，难以实施；自适应消噪法要求信号和噪声平稳且相互独立，在现场信号和噪 声不仅是非平稳，而且存在不同程度的相关性，因此，目前的设备振动信号处理方法对于变 工况设备振动信号处理不理想，影响诊断设备故障隐患的准确性。
发明内容
本发明所要解决的技术问题是提供一种设备振动信号的处理方法，该方法对变工况设 备振动信号进行处理后，可获得利于判断的波形图，为准确诊断设备故障隐患创造条件。 本发明解决上述技术问题所采用的技术方案是 一种设备振动信号的处理方法，它包括以下步骤 步骤(1)、测得设备振动信号的步骤； 歩骤(2)、对设备振动信号进行处理的步骤，具体为
步骤(2.1)、用第二代小波包技术将设备振动信号分解为四个以上独立的频带； 步骤(2.2)、分别计算分解后的各频带能量；
步骤(2.3)、选取能量最大的频带，对其对应的第二代小波包分解信号做自回归谱处理， 得到自回归谱图。
上述方案中，步骤(2.2)具体为:
分别计算分解后的各频带能量，归一化的各个频带的能量，画出能量特征直方图。 与现有技术相比，本发明具有以下优点
1、采用第二代小波包技术，除了具有传统小波分析的优点外，还克服了传统小波分析的 不足。第二代小波方法相对于传统小波算法而言，是一种更为快速、有效的小波变换实现方 法，它不依赖于Fourier变换，完全在时域完成了对双正交小波滤波器的构造。这种构造方法在结构化设计和自适应构造方面的突出优点弥补了传统频域构造方法的不足。由第二代小 波变换得到的第二代小波包技术能够把平稳的或非平稳的信号映射到一个由小波伸縮而成的 一组基函数上，信息量完整无缺，在通频范围内的得到分布在不同频带内的分解序列，具有 对非平稳信号进行局部化分析的功能，弥补了现有处理方法的缺陷，也弥补了小波分析只对 信号低频进行细致刻画，对信号的高频的分析不够细致的缺陷。考虑到第二代小波包分解经 过一层分解后数据的个数减半的特点，本发明结合自回归谱分析用于提取故障信息。自回归 谱分析较傅里叶分析得到的频率成分更精确而且在处理小量数据时更能体现其优越性。
2、自回归谱处理在短数据的情况下较FFT能取得更好的效果，能有效的提取信号中的频 率信息，获得利于判断的波形图，为准确诊断设备故障隐患创造条件。
图1为测得的设备振动原始信号时域波形 图2为三层小波包分解示意图 图3、 4为第二代小波分析的分解过程示意图 图5为归一化后各频带能量的直方图 图6 、 7为第二代小波分析的重构过程示意图 图8为第二代小波包三层分解第一频段自回归谱图
具体实施例方式
本发明设备振动信号的处理方法实施例，本实施例用于对变工况低速重载变速箱的振动
信号进行处理，本实施例包括以下步骤
步骤(1)、测得设备振动信号的步骤，具体为
将SRI50声发射固定设置在变速箱的外壳上，声发射的位置尽量靠近变速箱主要轴承，
声发射的输出端通过A/D转换输入计算机，采样频率,为1000Hz，计算机得到的振动信号时 域波形如图1所示，该时域波形无法直接用于诊断设备故障隐患。 步骤(2)、对设备振动信号进行处理的步骤，具体为
步骤(2.1)、用第二代小波包技术将设备振动信号分解为八个独立的频带；分解示意图 如图2所示，它为3层小波包分解。
第二代小波包分解过程由剖分、预测和更新组成。如图3、 4所示，设数据序列 S = e Z}，基于提升模式的第二代小波变换的分解过程如下，其中尸(，)为预测器，C/(，)
为更新器。
A、剖分，将数据序列^O0AeZ)分为奇样本序列^(^)和偶样本序<formula>formula see original document page 4</formula>B、 预领(j，设户O)为预测器，用、(/t)预测&(/t)，定义预测偏差为细节信号d(A:)。 则细节信号序列O = W(",/t e Z}
C、 更新，设t/(，)为更新器，在细节信号d(A:)的基础上更新^(W，其结果定义为逼近信 号c("。
c(A;)、(it) + C/[牟)]A:sZ 则逼近信号序列为C = {c(",A: e Z}
本实施例利用第二代小波包技术将信号无冗余、无疏漏、正交地分解到独立的频带。本 实施例信号的采样频率,为1000Hz，则其分析频率500Hz，对信号进行3层第二代小波包分 解，分解得到八个频带依次是
、 [62.5Hz， 125Hz]、 [125Hz， 187. 5Hz] 、 [187. 5Hz， 250Hz]、 [250Hz， 312.5Hz]、 [312.5Hz， 375Hz]、 [375Hz， 437. 5Hz]、 [437. 5Hz， 500Hz]。
步骤(2.2)、分别计算分解后再重构的各频带能量；具体为分别计算分解后的各频带 能量(假设^(")为3层第二代小波包分解后[O， 62.5Hz]频带的信号，贝IJ[O， 62.5Hz]频带的 能量计算为
e为能量，n为^0i)的数据长度。同理其它七个频带依次为5 、、 53( )、 (n)、 a(")。
归一化的各个频带的能量，画出能量特征直方图(如图5所示)，从直方图中看出第一频 带[O， 62.5Hz]所占百分比最大。
步骤(2.3)、选取能量最大的频带——
频带，对其对应的第二代小波包分 解信号做自回归谱处理，得到自回归谱图——
频带的自回归谱图(如图8所示)。
如图6、 7所示，第二代小波重构过程为分解过程的逆过程，由恢复更新、恢复预测和合 并组成。
se(A;) = c(A:)-
5。(" = , +化(/1;)] A; e Z
由奇偶样本序列合并构成重构信号s 。 假设预测器长度w二2，更新器长度#=4。 自回归谱处理的原理
设通过第二代小波包分解后的某一频段序列为{&}， k=l， 2，…，N,，其M阶自回归模型AR(M)为
S仏一,+ "
其中A， ^，……，&为自回归模型系数A为白噪声输入且^ NID (0，《)。N,为分界 了/层后每个频道内的数据长度。
自回归模型的阶数由经验公式AIC指标定阶，AK> ^式中[2
N表示采样点数，M表示模型阶数，/^表示第M阶AR模型预报误差， 当取得AIC值最小时的阶数M作为AR模型的最佳阶数。自回归系数^通过Marple算法 求出。则第二代小波包分解得到的第/层第/序列的M阶自回归普表达式如下
式中，采样频率厶、起始频率,,与第二代小波包分解层数有关。原始数据采样频率为/w ， 数据长度为N,则第二代小波包分解第/层时可得到位于不同频段的2'个序列，它们的采样频
率为厶=2—每个序列的长度为N,=2W*N，起始频率按小波包理论为
A = _/L* (i-l) /2。
从自回归谱图(如图8所示)中看到了 4.364HZ频率成分存在。这个频率成分与轴承的 轴频4. 469Hz非常接近，对此我们可以得到信息是变速箱存在故障隐患，建议停产检修变 速箱轴承的轴。
权利要求
1、一种设备振动信号的处理方法，其特征在于它包括以下步骤步骤(1)、测得设备振动信号的步骤；步骤(2)、对设备振动信号进行处理的步骤，具体为步骤(2.1)、用第二代小波包技术将设备振动信号分解为四个以上独立的频带；步骤(2.2)、分别计算分解后的各频带能量；步骤(2.3)、选取能量最大的频带，对其对应的第二代小波包分解信号做自回归谱处理，得到自回归谱图；步骤(3)、根据自回归谱图判断设备故障隐患的步骤。
2、 如权利要求1所述的处理方法，其特征在于步骤(2.2)具体为-分别计算分解后的各频带能量，归一化的各个频带的能量，画出能量特征直方图。
全文摘要
本发明涉及一种设备振动信号的处理方法，它包括以下步骤步骤(1)、测得设备振动信号的步骤；步骤(2)、对设备振动信号进行处理的步骤，具体为步骤(2.1)、用第二代小波包技术将设备振动信号分解为四个以上独立的频带；步骤(2.2)、分别计算分解后的各频带能量；步骤(2.3)、选取能量最大的频带，对其对应的第二代小波包分解信号做自回归谱处理，得到自回归谱图；步骤(3)、根据自回归谱图判断设备故障隐患的步骤。本发明处理方法对变工况设备振动信号进行处理后，可获得利于判断的波形图，为准确诊断设备故障隐患创造条件。
文档编号G01H17/00GK101634589SQ20091006367
公开日2010年1月27日 申请日期2009年8月21日 优先权日2009年8月21日
发明者辉 叶, 汤文亮, 辉 罗, 胡邦喜, 胥永刚, 苏善斌, 邓崎琳, 邹江华, 马海龙, 高立新 申请人:武汉钢铁(集团)公司;北京工业大学