一种优化iceemdan的管道泄漏声波信号检测方法
技术领域
1.本发明属于管道泄漏次声波信号的处理领域，具体涉及一种优化iceemdan的管道泄漏声波信号检测方法。
背景技术：2.管道运输因其具有经济型、便捷性，在工业上得到了众多应用。长距离输送管道大多会铺设在环境比较恶劣的地方，管道容易受到腐蚀和破坏等，这就导致管道极易出现泄漏。管道泄漏时，管道内的介质会对管道周边的环境造成不可挽回的破坏，还会造成企业的经济损失。因此，对长距离管道中的泄漏信号进行采集分析尤为重要。
3.次声波具有良好的抗衰减性，能够进行远距离的传播。负压波法具有较高的灵敏度和不易受噪声影响，但是负压波法无法反映管道微小的泄漏，而次声波信号可以反映管道的微小泄漏。次声波法在传播时，噪声会对其产生一定的影响。所以，需要对次声波信号进行去噪。
4.对次声波进行去噪的方法有很多，ceemdan是在emd的基础上改进而来的。它借用了eemd的加入高斯白噪声和多次叠加的思想，较好的解决了eemd存在的模态混叠现象。但是ceemdan分解后，信号中仍然存在一定的噪声残余。iceemdan是对ceemdan算法的优化改进得到的。iceemdan不是直接添加高斯白噪声，而是选择白噪声被emd分解后的第k阶imf。这就解决了ceemdan中残留噪声较多的问题。
技术实现要素：5.本发明的目的是提供一种优化iceemdan的管道泄漏声波信号检测方法，这种优化iceemdan的管道泄漏声波信号检测方法能够对管道中的泄漏声波信号进行检测和处理，提高对管道的泄漏检测能力，并准确判断管道的工况。
6.本发明提供了一种优化iceemdan的管道泄漏声波信号检测方法，包括：步骤1，采用声压传感器获取管道中的次声波信号；步骤2，用iceemdan分解采集的次声波信号，残余值不能在分解时停止分解，得到多个个本征模态函数imf；步骤3，计算各imf的mpe值，根据mpe值对imf进行筛选；步骤4，对信号进行互谱分析，得到信号的互谱相位谱和相干函数，通过互谱相位谱和相干函数确定信号的特征频带，根据特征频带选取符合要求的imf；步骤5，将满足mpe值要求并且频率位于特征频带范围内的imf用于管道泄漏声波信号的重构，得到去噪后的声波信号。
7.上述方案中步骤3的具体方法为：对采集的信号进行多尺度的粗粒化并进行重构，得到：
(1) (2)对粗粒化并重构后的信号进行升序排列，计算信号的排列熵并对其进行归一化处理：(3)(4)(5)根据预置的mpe数值来选择imf，mpe小于预设值的imf中含有更多的有用信息，mpe大于预设值的imf中，含有的噪声更多。
8.上述方案中步骤4的具体方法为：原始信号用来表示，管道两端信号分别用x1和x2，它们之间的互相关函数为： (1)信号之间的相互依赖程度可以用互谱表示，两个信号之间的互谱为：
ꢀꢀꢀ
(2)一般情况下，会采用welch平均周期图标法来对互谱进行估计表示： (3)在互谱相位谱中，若相位随着频率呈现线性变化，那么信号之间具有较大的相关性。两个信号的想干函数可以用下式表示：(4)在相干函数中，当两个信号的频率越大时，表明两个声波信号的相关程度越高。
9.本发明具有下列有益效果：本发明采用基于互谱分析和mpe的iceemdan算法对传感器采集的次声波信号进行去噪处理，iceemdan将次声波信号分解为多个本征模态函数imf，对各个imf的mpe数值进行
计算，并根据预设值对imf进行筛选。对管道两端传感器采集的次声波信号进行互谱分析，得到两个信号之间的互谱相位谱和相干函数，根据互谱相位谱和相干函数来确定特征频带，进而筛选含有有用信息多的imf。最后，将mpe值小于预设值并且位于特征频带内的imf用于次声波信号的重构。重构后的次声波信号包含的有用信息多，能够准确反映管道的状态。
附图说明
10.以上步骤大概的介绍了本发明的技术手段，为了能更清楚的了解本发明的技术手段，以下为本发明实施例所需要的配图。
11.图1为本发明提供的基于互谱分析和mpe的iceemdan的管道泄漏声波信号检测方法的流程图；图2为本发明所提供的传感器采集的次声波信号去噪前的曲线图；图3为本发明所提供的传感器采集的次声波信号去噪后的曲线图；图4为本发明所提供的得到的各个imf的mpe；图5为互谱分析得到的互谱相位谱图；图6为互谱分析得到的相干函数图。
具体实施方式
12.以下实施例只是用来说明本发明，但不用来限制本发明。
13.如图1所示，本实施例提供了基于互谱分析和mpe的iceemdan的管道泄漏声波信号检测方法，包括：步骤1，采用声压传感器获取管道中的次声波信号；次声波信号均在管道实验场地中的管道中进行采集，管道的长度为3000m，管径为dn100，管道两端的压力分别为2.0kpa和0.43kpa，管道的首端和末端分别安装一个声压传感器，采样频率为1000hz，在管道的200m、1000m和1800m处各有一个泄漏阀门，用于模拟泄漏的产生。在所有阀门关闭时，采集管道内的声波信号，即正常工况下的声波信号。打开管道的泄漏阀门模拟泄漏，去噪之前的声波信号曲线如图2所示；步骤2，采用iceemdan对管道两端传感器所采集到的次声波信号进行分解，直到残余不能再分解时停止分解；本案例中对信号进行iceemdan分解的具体步骤如下：定义原始信号为x, em(.)为emd分解后的第m阶imf，z(.)是次声波信号的局部均值，w(i)是高斯白噪声。
14.向原始数据中添加一组高斯白噪声，构成序列，并得到第一个残余值。
15.第一阶imf为：。
16.将剩余的信号作为新的源信号，继续添加高斯白噪声，根据局部均值分解得到。
17.计算第k阶残余值。
18.步骤3，计算各imf的mpe值，根据mpe值对imf进行筛选；本案例中各imf的mpe值求解步骤如下：对采集的信号进行多尺度的粗粒化并进行重构，得到； (1) (2)对粗粒化并重构后的信号进行升序排列，计算信号的排列熵并对其进行归一化处理；
ꢀꢀꢀ
(3)
ꢀꢀꢀ
(4)
ꢀꢀꢀ
(5)根据预置的mpe数值来选择imf，mpe小于预设值的imf中含有更多的有用信息，mpe大于预设值的imf中，含有的噪声更多，各个imf的mpe值如图4所示。
19.步骤4，对信号进行互谱分析，得到信号的互谱相位谱和相干函数，通过互谱相位谱和相干函数确定信号的特征频带，根据特征频带选取符合要求的imf；本案例中互谱分析的步骤如下：原始信号用来表示，管道两端信号分别用x1和x2，它们之间的互相关函数为：
ꢀꢀ
(1)信号之间的相互依赖程度可以用互谱表示，两个信号之间的互谱为： (2)一般情况下，会采用welch平均周期图标法来对互谱进行估计表示：
ꢀꢀꢀ
(3)在互谱相位谱中，若相位随着频率呈现线性变化，那么信号之间具有较大的相关性。两个信号的想干函数可以用下式表示：
ꢀꢀꢀ
(4)在相干函数中，当两个信号的频率越大时，表明两个声波信号的相关程度越高互谱相位谱和相干函数如图5、图6所示。
20.步骤5，将满足mpe值要求并且频率位于特征频带范围内的imf用于管道泄漏声波信号的重构，得到去噪后的声波信号，去噪后的次声波曲线如图3所示。
技术特征：1.一种优化iceemdan的管道泄漏声波信号检测方法包括以下步骤：步骤1，采用声压传感器获取管道中的次声波信号；步骤2，用iceemdan分解采集的次声波信号，残余值不能在分解时停止分解，得到多个个本征模态函数imf；采用iceemdan方法对采集的次声波信号进行分解的步骤如下：（1）定义原始信号为x, e
m
(.)为emd分解后的第m阶imf，z(.)是次声波信号的局部均值，w
(i)
是高斯白噪声；（2）向原始数据中添加一组高斯白噪声，构成序列，并得到第一个残余值；（3）第一阶imf为：；（4）将剩余的信号作为新的源信号，继续添加高斯白噪声，根据局部均值分解得到；（5）计算第k阶残余值；步骤3，计算各imf的mpe值，根据mpe值对imf进行筛选；步骤4，对信号进行互谱分析，得到信号的互谱相位谱和相干函数，通过互谱相位谱和相干函数确定信号的特征频带，根据特征频带选取符合要求的imf；步骤5，将满足mpe值要求并且频率位于特征频带范围内的imf用于管道泄漏声波信号的重构，得到去噪后的声波信号。2.根据权利要求1所述的一种优化iceemdan的管道泄漏声波信号检测方法，其特征在于，步骤3的具体方法为：对采集的信号进行多尺度的粗粒化并进行重构，得到：(1) (2)对粗粒化并重构后的信号进行升序排列，计算信号的排列熵并对其进行归一化处理： (3) (4)
ꢀ
(5)根据预置的mpe数值来选择imf，mpe小于预设值的imf中含有更多的有用信息，mpe大于预设值的imf中，含有的噪声更多。3.根据权利要求1所述的一种优化iceemdan的管道泄漏声波信号检测方法，其特征在于，步骤4的具体方法为：原始信号用来表示，管道两端信号分别用x1和x2，它们之间的互相关函数为： (1)信号之间的相互依赖程度可以用互谱表示，两个信号之间的互谱为： (2)一般情况下，会采用welch平均周期图标法来对互谱进行估计表示： (3)在互谱相位谱中，若相位随着频率呈现线性变化，那么信号之间具有较大的相关性；两个信号的想干函数可以用下式表示： (4)在相干函数中，当两个信号的频率越大时，表明两个声波信号的相关程度越高。
技术总结本发明提供一种优化ICEEMDAN的管道泄漏声波信号检测方法，其涉及次声波信号的采集与去噪处理领域。方法包括：一、利用安装在管道上的声压传感器采集管道中的次声波信号；二、采用ICEEMDAN算法对采集得到的声波信号进行分解，直到残余信号不能再分解；三、分别计算各个模态函数IMF的多尺度排列熵MPE，根据得到MPE值选择符合要求的IMF；四、对信号进行互谱分析，得到信号的互谱相位谱和相干函数，通过互谱相位谱和相干函数确定信号的特征频带，根据特征频带选取符合要求的IMF；五、将满足MPE值要求并且频率位于特征频带范围内的IMF用于管道泄漏声波信号的重构，得到去噪后的声波信号。号。号。
技术研发人员：朱永强 郎宪明
受保护的技术使用者：辽宁石油化工大学
技术研发日：2022.05.04
技术公布日：2022/7/4
