矿用排水泵噪声多时域特征联合分析方法与流程

本发明涉及噪声分析技术领域，具体是矿用排水泵噪声多时域特征联合分析方法。

背景技术：

当下，常见的电机故障及健康诊断多数的是利用振动、电流等参数进行电机的故障诊断。此种方法目前已经比较成熟，应用非常多，可以查验绝大多数的电机健康状况。但是，采用振动及电流等参数进行故障诊断时，大多是在频域或时频域下进行的，此类方法判定过程较为复杂，不够直观，并且此类参数的采集受设备约束较高，某些场景信号的采集并不是很方便。

而近年来，利用音频信号进行电机故障诊断也逐渐进行人们的视野。音频信号更直接地反应的是振源的状态，可以更加直观地反应出电机的问题所在。同时，很多工况下对振动等信号的采集要求非常高，音频信号是非接触采集，不影响设备的工作，而且音频采集装置更加简单，采集速度更快，在电机健康诊断中可以更早地发现问题。

通过音频信号对水泵电机健康状况分析是十分可行的，不仅可以解决不同工业场景的环境束缚，同时，利用音频信号的多时域特征进行故障的联合分析，简化了分析过程，也具有较高的准确度。

技术实现要素：

本发明的目的是提出矿用排水泵噪声多时域特征联合分析方法，该方法利用水泵电机工作时的音频信号来进行电机的健康状态分析，通过对原始音频信号的滤波、预加重及分帧处理得到平滑连接的音频数据，再利用音频数据的峭度、短时能量、短时平均过零率三个时域特征对水泵电机的健康状态进行综合分析，将水泵电机的故障扼杀在萌芽之中，保证水泵电机的正常运行。

本发明为实现发明目的采用如下技术方案：

矿用排水泵噪声多时域特征联合分析方法，包括以下步骤：

通过低通滤波器对水泵电机的音频信号进行滤波处理，消除原始音频信号中的混叠失真和高频噪声；

对滤波后的音频数据进行预加重和分帧处理，减少音频数据高频段的信号衰减，将长段音频数据分成若干小段数据；

计算各段音频数据的峭度，根据计算出的峭度值判定水泵电机此刻的健康状况；

当通过峭度值无法对水泵电机健康状况进行判定时，计算音频数据的短时能量、短时平均过零率，利用多时域特征参量联合判定水泵电机的健康状况。

进一步的，低通滤波器采用gammatone滤波器，其时域表达式为：

式中，a是滤波器的增益，bi是滤波器的衰减因子，fi是滤波器的中心频率，φi是偏移相位，u(t)是阶跃函数。

进一步的，音频数据的峭度计算，具体步骤如下:

定义音频数据的n阶中心统计矩，具体为:

mn＝e[(x-μ)ⁿ],

式中，x为音频信号，e表示期望，μ为均值

对信号n阶中心统计矩来说，一阶中心矩是均值，二阶是方差，三阶是偏度，而四阶为峭度。因此，峭度定义如下：

式中，σ为标准差；

根据计算结果进行水泵电机的健康诊断，当峭度指标的绝对值大于8时，判定为电机存在故障，当峭度指标的绝对值介于3和8之间时，进行短时能量及短时过零率的计算，进一步判定水泵电机的健康状况，当峭度指标的绝对值小于3时，判定为水泵电机健康。

进一步的，计算音频数据的短时能量，公式如下：

yi(n)＝w(n)·x((i-1)·inc n),1≤n≤l,1≤i≤fn

式中，yi(n)是一帧的数值，w(n)是分帧函数，l是帧移长度，fn是总帧数

第i帧音频数据的yi(n)的短时能量e(i)为:

计算音频数据的短时平均过零率，公式如下：

结合计算出的短时能量和短时平均过零率对水泵电机的健康进行二次诊断。当短时能量、短时平均过零率超出正常范围时，则判定此时水泵电机处于故障状态，若短时能量或短时平均过零率仅有一项不在正常范围内，则判定此时的水泵电机处于异常状态。

本发明现有技术相比，其有益效果体现在：

1.利用音频信号进行水泵电机的健康诊断，较传统的故障诊断方式，更加直接地反应了故障时的问题所在，具有较高的诊断精度。

2.通过音频信号的多种时域特征联合判断水泵电机的健康状态，精度高、准确率高，可以更加全面的反映问题。同时，直接通过音频信号的时域频谱来进行健康诊断，比频域、时频域诊断过程更为简便。

3.对音频信号的处理加入了滤波、预加重及分帧，可以降低噪声的干扰，减少音频信号高频部分的损失，获得更加平滑、更精确的短时音频信号，方便处理。

附图说明

图1为本发明的工作流程示意图

图2为本发明的分帧示例图

图3为本发明的三种峭度值分布曲线图

其中，图2中，n表示音频数据帧长。

具体实施方式

下面结合附图来说明本发明专利的技术方案。应当理解，此处描述的具体实施方式仅仅是用以解释本次发明专利，并不限定于本发明。

所述的矿用排水泵噪声多时域特征联合分析方法，通过gammatone滤波器对水泵电机运行时的音频信号进行滤波，消除原始音频信号中的噪声，在将滤波后的音频数据进行预加重，减弱音频信号中高频部分的衰减，再通过分帧，将预加重后的音频数据分成短段数据来处理，增强信号的平稳性，最后利用时域中音频数据的多种特征进行联合判断，判断水泵电机的健康状况。

参见图1，图2，图3，其具体实施过程如下：

矿用排水泵噪声多时域特征联合分析方法，包括以下步骤：

(一)利用所述gammatone滤波器对水泵电机的音频信号进行滤波处理，将原始音频信号中的混叠失真和高频噪声进行滤除，减少干扰；

所述低通滤波器采用gammatone滤波器，其时域表达式为：

式中，a是滤波器的增益，bi是滤波器的衰减因子，fi是滤波器的中心频率，φi是偏移相位，u(t)是阶跃函数。

(二)对滤波后的音频数据进行预加重和分帧处理，减少音频数据中高频部分的衰减，同时，将长段的音频数据分成若干短段数据，使音频数据更加平滑，更容易处理；

1.音频数据预加重，具体步骤如下：

假设n时刻音频采样值为x(n)，预加重处理后的结果为：

yn(n)＝x(n)-αx(n-1),0.9＜α＜1

式中，α为预加重系数。

2.音频数据的分帧处理，具体步骤如下：

分帧步骤可根据音频数据的大小，合理选择分帧数，进行分帧处理，具体公式如下：

式中，w(n)是分帧函数，f(x)是音频数据总长，n是分帧数。

(三)利用处理好的音频数据，计算音频数据的峭度，根据计算结果进行水泵电机的健康诊断；

1.音频数据的峭度计算，具体步骤如下:

定义音频数据的n阶中心统计矩，具体为:

mn＝e[(x-μ)ⁿ]，

式中，x为音频信号，e表示期望，μ为均值

对信号n阶中心统计矩来说，一阶中心矩是均值，二阶是方差，三阶是偏度，而四阶为峭度。因此，峭度定义如下：

式中，σ为标准差。

2.利用计算出的峭度的绝对值进行水泵电机的健康诊断，当峭度指标的绝对值大于8时，判定为电机存在故障；当峭度指标的绝对值介于3和8之间时，进行短时能量及短时过零率的计算，进一步判定水泵电机的健康状况；当峭度指标的绝对值小于3时，判定为水泵电机健康。

(四)当峭度值无法直接判断水泵电机的健康状况时，计算音频数据的短时能量、短时平均过零率等时域重要特征参量，对水泵电机的健康状况进行联合判定；

1.计算音频数据的短时能量，公式如下：

yi(n)＝w(n)·x((i-1)·inc n),1≤n≤l,1≤i≤fn

式中，yi(n)是一帧的数值，w(n)是分帧函数，l是帧移长度，fn是总帧数

第i帧音频数据的yi(n)的短时能量e(i)为:

2.进一步计算音频数据的短时平均过零率，公式如下：

3.根据计算的短时能量和短时平均过零率再次判断水泵电机的健康状况。当短时能量、短时平均过零率超出正常范围时，则判定此时水泵电机处于故障状态，若短时能量或短时平均过零率仅有一项不在正常范围内，则判定此时的水泵电机处于异常状态。

以上所述的具体实施方式，对本发明的目的、技术方案和有益效果做了进一步的详细说明，所属领域技术人员应该明白，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。