掌桥专利:专业的专利平台
掌桥专利
首页

一种基于CS-MUSIC矿井提升机故障源定位分析方法

文献发布时间:2023-06-19 16:06:26



技术领域

本发明涉及音频信号分析及定位领域,具体是一种基于CS-MUSIC矿井提升机故障源定位分析方法。

背景技术

矿井提升机在矿井的运输工程中承担着极其重要的作用,提升机的安全和稳定的性能影响着井下工作人员的安全以及煤矿的经济效益。目前,矿井提升机监控设备不完善,发生故障时维修人员很难即时排查出故障源头,大大影响了矿井工作的效率。当矿井提升机发生故障时,其声音也会随着提升机的状态发生改变,因此通过分析提升机发出的音频特征来检测设备是否出现故障。而现有的基于声音检测矿井提升机故障主要侧重于预测识别故障类型,这种方式的缺点是发生故障时要进一步寻找故障的源头,浪费大量人力物力。

为此,本发明提出了一种基于CS-MUSIC矿井提升机故障源定位分析方法解决上述问题。

梅尔频率倒谱系数(Mel-Frequency Cepstral Coefficient,MFCC)是声音信号频域分析方法的一种,不限制其输入的音频信号类型,具有更好的鲁邦性,抗干扰性强。MUSIC算法(Multiple Signal classification)是一种只要已知天线阵的布阵形式,不管阵元是否是等间距分布,构建谱函数求得波峰估计值,就会得到高分辨率的定位结果的算法发。布谷鸟算法(Cuckoo Search,CS)通过模拟布谷鸟的寄生育雏的方式,解决最优解相关的问题,具有强大的搜寻能力,需要的参数少等优点。

发明内容

本发明的目的是提出一种基于CS-MUSIC矿井提升机故障源定位分析方法,对矿井提升机的音频信号进行实时监测和识别,通过对音频信号的预处理,梅尔频率倒谱系数提取以及故障识别和故障定位,实现对矿井提升机的监测故障诊断和故障定位等功能,可快速有效的定位矿井提升机故障发生的具体位置,从而保证提升机的安全稳定的运行。

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

一种基于CS-MUSIC矿井提升机故障源定位分析方法,包括以下步骤:

步骤1、利用五元十字形麦克风阵列采集并存储矿井提升机完整运行过程的音频信号,对采集到的音频信号进行预处理,其中:包含预加重、分帧和加窗;

步骤2、对预处理后的音频信号进行MFCC特征参数提取,包括:FFT变换、Mel滤波器滤波和取对数DCT;

步骤3、对识别出故障的音频用MUSIC算法进行声音定位;

步骤4、用CS算法优化步骤3得到的定位分布峰值结果;

进一步的,步骤1中的预加重、分帧和加窗具体包括:

步骤1-1、预加重:通过一阶FIR高通滤波器以此来增强矿井提升机音频信号的高频部分,以此来减小噪声对音频信号的影响。结果为:

z'(t)=z(t)-hz(t-1) (1)

其中:z(t)表示第t个采样点出的音频信号,h为为加强系数,h=0.97;

步骤1-2、分帧加窗:先对z'(t)先进行分帧,以256点为一帧,步长为512,对音频信号每一帧采用汉明加窗,变为更平滑的z(m);

进一步的,步骤2中的短时傅里叶变换、Mel滤波器滤波和离散余弦变换具体包括:

步骤2-1、FFT。对每帧提升机音频信号z(m)进行快速傅立叶变换,时域音频信号z

Z(i,j)=FFT[z

步骤2-2、Mel滤波器滤波。设计Mel滤波器为24阶让每一个都转换为Mel刻度下的值,其频率公式表达为:

E(i,j)=[Z(i,j)]

f

其中:i为第i帧;j为频域中的第j条谱线,E(i,j)代表每一帧的能量。f

步骤2-3、取对数DCT。对经过转换后能量值进行对数变换,再进行DCT变换,得到MFCC参数。具体为:

其中:

进一步的步骤3中识别出故障的音频用MUSIC算法进行声音定位具体为:

步骤3-1、、将采集的五路声音信号进行二次分帧,长度为L,步长为分帧长度的一半,然后对分帧后的信号做FFT变换。公式如下:

其中:k=0,1,…L-1,

可以得到音频信号频域模型为:

Z(f

其中:

步骤3-2、求协方差矩阵的估计值。

其中:N为采样个数。

步骤3-3、对R

U

从而得到MUSIC的谱估计函数为:

其中:

通过θ在观察扇面

进一步的步骤4中用CS算法优化步骤3得到的定位分布峰值结果具体为:

步骤4-1、初始化参数:设置迭代次数M,巢穴数量n,目标参数即提升机故障信号的方位γ,并在搜索可行范围内随机初始化巢穴位置,

步骤4-2、计算每个巢穴的适应度值

所述适应度值计算公式为:

步骤4-3、第一次更新巢穴位置

第一次更新巢穴位置公式为:

其中:α为步长控制因子,

步骤4-4、判断步骤4-3所述新巢穴是否为“劣质”巢穴,如果不是,则新巢穴

其中:

步骤4-5、计算步骤4-4所述新巢穴

步骤4-6、判断是否当前迭代t≥最大迭代次数T,若否,则返回到步骤4-3,若是,则输出适应度值最小的巢穴位置,即为目标参数γ的最优解。

优选的步骤4-4中所述“劣质”巢穴通过如下方法判断:若当前适应度值

本发明现有技术相比,其有益效果体现在:

1.通过应用提升机运行时的声学信号的对提升机的故障进行分析,声音采集设备成本较低、操作简单且实时性高。

2.通过MFCC特征提取音频信号的信息,该算法抗干扰性能强,具有更好的鲁棒性。

3.通过应用CS算法对MUSIC定位结果寻优,大大提高了定位结果的精准性和准确性。

附图说明

图1为本发明的一种基于CS-MUSIC矿井提升机故障源定位分析方法的流程示意图。

图2为本发明中音频信号的原始信号波形图。

图3为本发明中音频信号预处理后的示意图。

图4为本发明中MUSIC定位部分算法流程图。

图5为本发明中五元十字形阵列建立在直角坐标系下的结构图。

具体实施方式

以下通过具体实施例来对本发明做进一步解释说明。

如图1所示,本发明的一种基于CS-MUSIC矿井提升机故障源定位分析方法,具体包括如下步骤:

步骤1、利用五元十字形麦克风阵列采集并存储矿井提升机完整运行过程的音频信号,对采集到的原始音频信号z(t)如图2所示进行预处理,其中:包含预加重、分帧和加窗,预处理后的图像如图3所示。

步骤1-1、预加重:通过一阶FIR高通滤波器以此来增强矿井提升机音频信号的高频部分,以此来减小噪声对音频信号的影响。结果为:

z'(t)=z(t)-hz(t-1) (1)

其中:z(t)表示第t个采样点出的音频信号,h为为加强系数,h=0.97;

步骤1-2、分帧加窗:先对z'(t)先进行分帧,以256点为一帧,步长为512,对音频信号每一帧采用汉明加窗,变为更平滑的z(m)

步骤2、对预处理后的音频信号进行MFCC特征参数提取,包括:FFT变换、Mel滤波器滤波和取对数DCT;

步骤2-1、FFT。对每帧提升机音频信号z(m)进行快速傅立叶变换,时域音频信号z

Z(i,j)=FFT[z

步骤2-2、Mel滤波器滤波。设计Mel滤波器为24阶让每一个都转换为Mel刻度下的值,其频率公式表达为:

E(i,j)=[Z(i,j)]

f

其中,i为第i帧;j为频域中的第j条谱线,E(i,j)代表每一帧的能量。f

步骤2-3、取对数DCT。对经过转换后能量值进行对数变换,再进行DCT变换,得到MFCC参数。具体为:

其中:

步骤3、对识别出故障的音频用MUSIC算法进行声音定位部分流程如图4所示。

步骤3-1、将采集的五路声音信号进行二次分帧,长度为L,步长为分帧长度的一半,然后对分帧后的信号做FFT变换。公式如下:

其中:k=0,1,…L-1,

可以得到音频信号频域模型为:

Z(f

其中:

步骤3-2、求协方差矩阵的估计值。

其中:N为采样个数。

步骤3-3、对R

U

从而得到MUSIC的谱估计函数为:

其中:

通过θ在观察扇面

步骤4中用CS算法优化步骤3得到的定位分布峰值结果。

步骤4-1、初始化参数:设置迭代次数M,巢穴数量n,目标参数即提升机故障信号的方位γ,并在搜索可行范围内随机初始化巢穴位置,

步骤4-2、计算每个巢穴的适应度值

所述适应度值计算公式为:

步骤4-3、第一次更新巢穴位置

第一次更新巢穴位置公式为:

其中:α为步长控制因子,

步骤4-4、判断步骤4-3所述新巢穴是否为“劣质”巢穴,如果不是,则新巢穴

其中:

步骤4-5、计算步骤4-4所述新巢穴

步骤4-6、判断是否当前迭代t≥最大迭代次数T,若否,则返回到步骤4-3,若是,则输出适应度值最小的巢穴位置,即为目标参数γ的最优解。

优选的步骤4-4中所述“劣质”巢穴通过如下方法判断:若当前适应度值

需要说明的是,以上公开的仅为本发明的具体实例,根据本发明提供的思想,本领域的技术人员能思及的变化,都应落入本发明的保护范围内。

技术分类

06120114704879