优化变分模态分解的分布式光纤声传感系统的信号处理方法

1.本发明涉及信号源监测领域，具体为一种优化变分模态分解的分布式光纤声传感系统的信号处理方法。


背景技术：

2.决定冶金和煤矿行业持续高效生产的关键技术是保障技术路线中各个环节尽可能的有效运行。原料输送胶带机以输送能力大、输送距离长、结构简单、维护方便、通用性强等优点，成为散装材料运输中的首选设备。但是原料输送胶带机操作、维护作业具有极大风险，且一个简单的故障或者某一小部分出现异常警报，都会引起后续工作暂停甚至瘫痪等问题。修复某些故障可能需要很长时间，需要工作人员在故障诊断、确保安全、首次操作决策、拆解、材料运输、维修和装配等领域进行大量参与。然而，胶带机托辊早期故障特征信号微弱、设备运行环境恶劣、背景噪声强等因素导致监测系统无法检测到这些故障。因此，强背景噪声条件下原料输送胶带机的托辊故障、胶带撕裂等诊断与运行趋势预测尤为必要。
3.原料输送胶带机的实时早期故障监测的传统手段是采用点式的电学加速度或振动传感器，但是由于原料输送胶带机传输距离长，现场工况环境恶劣，传统的点式电学传感器往往无法实现长距离连续监测和故障的精准定位。相比于传统点式电学传感器，基于相敏型光时域反射仪的分布式光纤声传感系统在原料输送胶带机的故障监测定位应用领域具有长距离分布式检测，环境适应性强等突出优势，但是由于监测现场背景噪声强等因素的影响，使得常规的信号时频分析方法无法实现对故障信息的有效检测和准确定位。


技术实现要素：

4.针对上述胶带机轴承早期损伤状态进行及时诊断和实时定位难题，本发明提供了一种优化变分模态分解的分布式光纤声传感系统的信号处理方法。本发明从强背景噪声环境中采集声传感系统发射脉冲并获取光纤长度的声振信号，通过变分模态分解处理所获取的声振信号并对存在故障全部段数实时报警。本发明能够很好地解决强背景噪声环境中的微弱故障特征信号提取难题，能够在低信噪比的原始声振信号中检测出光纤沿线全部存在故障段数并实时定位报警，实现基于分布式光纤声传感对多点故障的实时判别和定位。
5.一种优化变分模态分解的分布式光纤声传感系统的信号处理方法，包括以下步骤：
6.步骤s1：分布式光纤声传感系统对传感光纤发射探测光脉冲，对传感光纤沿线的声振信息进行实时定量监测，并获取原始声振信号s(p,q)，其中p为获取脉冲个数，q为光纤沿线声振信号获取段数。
7.步骤s2：对步骤s1获取的原始声振信号s(p,q)按照光纤沿线长度分段s(p,i)(1≤i≤q)进行优化变分模态分解处理，并记录存在故障全部段数；
8.步骤s21：通过对采集数据的变分模态分解分析，根据变分模态分解结果确定变分模态分解的最佳分解层数k与惩罚因子α，无需通过算法对每一分段数据分别优化分解参
数，使用测试数据进行测试并成功提取故障信号，从而达到算法简化的目的；
9.步骤s22：通过希尔伯特(hilbert)变换计算每个数组s(p,i)的解析信号，各解析信号分别与预估中心频率相乘，将每个模态的频谱调制到相应的基频带得到的各个模态解调信号；
10.步骤s23：计算步骤s22得到的各个模态解调信号的梯度的平方l2范数，根据约束变分模型估计出各个模态解调信号的带宽，受约束的变分模型表达式为：
[0011][0012]
其中，δ(t)表示狄利雷克分布函数，*表示卷积，{μk}＝{μ1,μ2,
…
,μk}表示k个模态解调信号，{ωk}＝{ω1,ω2,
…
,ωk}代表每个模态解调信号的中心频率，s为原始信号，t为时间；
[0013]
步骤s24：引入α和拉格朗日(lagrange)乘法算子λ，将约束性变分问题变为非约束性变分问题，α参数的作用是保证重构信号的精度，扩展的lagrange表达式为：
[0014][0015]
式中，λ(t)是严格约束条件；
[0016]
步骤s25：利用交替方向乘子算法求解，通过交替更新和λ
n+1
来求解扩展lagrange表达式的“鞍点”最小值，最优解为各本征模态分量(imfs)及中心频率相应的变量，更新表达式分别为：
[0017][0018][0019][0020]
其中，∧表示傅里叶变换，n为迭代次数，τ为保真系数。
[0021]
在迭代求解过程中各imfs的频率及带宽不断更新，直至满足迭代停止条件结束整个循环，最终根据实际信号的频域特征完成分割，并通过傅里叶逆变换将转为时域imfs。
[0022]
步骤s26：判断每个数组s(p,i)中是否存在故障并记录存在故障段数，循环直至判断全部光纤沿线段数q时结束。
[0023]
步骤s3：输出步骤s2记录的存在故障全部段数，并实时判别和定位。
[0024]
本发明旨在通过优化变分模态分解算法：在保证处理分布式光纤声传感系统采集的多通道声信息处理实时性的同时对环境噪声进行抑制和对目标信号进行识别，解决现有
基于分布式光纤声波传感系统中无法有效提取强背景噪声下皮带机轴承早期损伤振动信号中的故障特征频率，难以对胶带机轴承早期损伤状态进行及时诊断和实时定位问题。
[0025]
综上所述，本发明有效提高了处理后声振信号的信噪比，降低系统误报与漏报振动事件的发生，对光纤沿线实现基于分布式光纤声波传感系统在强背景噪声环境中多点故障特征信号的实时判别和精准定位。
附图说明
[0026]
图1是本发明信号处理方法流程图；
[0027]
图2是本发明数据变分模态分解处理示意图；
[0028]
图3是安静环境中和强背景噪声环境中测试点a的模拟故障信号强度图；
[0029]
图4是实施例强背景噪声环境中测试点a模拟故障特征频率为86hz的处理结果示意图；
[0030]
图5是强背景噪声环境中测试点a模拟故障特征频率为71hz的处理结果示意图；
[0031]
图6是强背景噪声环境中测试点a和测试点a、b的定位报警示意图。
具体实施方式
[0032]
下面结合附图和实施例对本发明做进一步的详细说明。
[0033]
如图1所示，一种优化变分模态分解的分布式光纤声传感系统的信号处理方法，具体实施步骤为：
[0034]
步骤s1：分布式光纤声传感系器对传感光纤发射探测光脉冲，对传感光纤沿线的声振信息进行实时定量监测，并获取原始声振信号s(p,q)，其中p＝5120，q＝800；
[0035]
步骤s2：在本实施例中，对步骤s1获取的原始声振信号s(5120,800)按照光纤沿线长度分段s(5120,i)(1≤i≤800)进行优化变分模态分解处理，并记录存在故障全部段数，对原始信号进行优化变分模态分解处理如图2所示，具体包括以下步骤：
[0036]
步骤s21：通过对采集数据的变分模态分解分析，确定变分模态分解的分解层数k＝3与惩罚因子α＝1000；
[0037]
步骤s22：通过hilbert变换计算每个数组s(5120,i)的解析信号，解析信号分别与预估中心频率将每个模态的频谱调制到相应的基频带得到的各个模态解调信号；
[0038]
步骤s23：计算以上解调信号的梯度的平方l2范数，根据约束变分模型估计出各个模态解调信号的带宽，受约束的变分模型表达式如式(s23)所示；
[0039][0040]
其中，δ(t)表示狄利雷克分布函数，*表示卷积，{μk}＝{μ1,μ2,
…
,μk}表示k个模态解调信号，{ωk}＝{ω1,ω2,
…
,ωk}代表每个模态解调信号的中心频率，s为原始信号；
[0041]
步骤s24：引入α和lagrange乘法算子，将约束性变分问题变为非约束性变分问题，α参数的作用是保证重构信号的精度，扩展的lagrange表达式如式(s24)所示：
[0042][0043]
式中λ(t)则是严格约束条件，λ为lagrange算子；
[0044]
步骤s25：利用交替方向乘子算法求解，通过交替更新和λ
n+1
来求解扩展lagrange表达式的“鞍点”最小值，最优解为各imfs及中心频率相应的变量，更新表达式如下：
[0045][0046][0047][0048]
其中，∧表示傅里叶变换，n为迭代次数，τ为保真系数。
[0049]
在迭代求解过程中各imfs的频率及带宽不断更新，直至满足迭代停止条件结束整个循环，最终根据实际信号的频域特征完成自适应分割，并通过傅里叶逆变换将转为时域imfs。
[0050]
步骤s26：判断每个数组s(5120,i)中是否存在故障并记录存在故障段数，循环直至判断全部光纤沿线段数q＝800时结束。
[0051]
步骤s3：输出步骤s2记录的存在故障全部段数(即各个故障段数对应的i)，并实时判别和定位。
[0052]
在本实施例中：
[0053]
通过步骤s1采集数据得到安静环境下和强背景噪声条件下的测试点a(i＝78)敲击模拟故障信号如图3所示，其中图3(a)和(b)分别为安静环境下和强背景噪声条件下的光纤沿线测试点a附近的敲击信号强度图，由图3(a)可知安静环境条件下是能够判断故障信号(如图中点划线框选所示)存在的，而由图3(b)可知强背景噪声下是无法判断是否存在故障信号的。
[0054]
通过步骤s2中的变分模态分解降噪处理后得到测试点a强背景噪声环境中的模拟故障特征频率信号，模拟故障频率86hz处理结果如图4所示，如图可见本发明方法成功提取了该测试点的故障特征频率。进一步在测试点a强背景噪声环境中的模拟故障频率71hz，并经上述方法步骤s2中降噪处理后得到结果如图5所示，如图可见该方法成功提取了该测试点的故障特征频率。
[0055]
通过上述方法采集和降噪处理后得到强背景噪声环境中在测试点a和同时在测试点a、b(i＝228)的模拟故障信号定位报警结果分别如图6(a)和(b)所示，如图可见测试在测试点a和同时在测试点a、b分别发生报警。
[0056]
通过以上实施例可见，本发明对采集的声振信号s(p,q)按照光纤沿线长度分段s(p,i)进行变分模态分解处理并记录存在故障的各段数，利用优化变分模态分解对强背景
噪声环境中光纤沿线声振信号的合理分解，提高了声振信号的信噪比，降低系统误报与漏报振动事件的发生，对光纤沿线基于分布式光纤声波传感系统在强背景噪声环境中多点故障特征信号的实时判别和精准定位；并最终输出光纤长度存在故障报警段数。
[0057]
本发明通过优化变分模态分解算法在保证处理分布式光纤声传感系统采集的多通道声信息处理实时性的同时对环境噪声进行抑制和对目标信号进行识别，解决了现有基于分布式光纤声波传感系统中无法有效提取强背景噪声下皮带机轴承早期损伤振动信号中的故障特征频率，难以对胶带机轴承早期损伤状态进行及时诊断和实时定位问题。