本发明属于管道泄漏检测领域,涉及一种基于自适应多元变分模态分解的流体管道泄漏声学定位方法,主要应用于天然气、石油、供水、热力管道等的泄漏检测定位。
背景技术:
随着我国的天然气、水、石油等管道网络不断完善,其中管道网络规模不断上升,因管道材料的缺陷、老化、腐蚀,地质沉降损害以及不及时的维修等问题导致泄漏率递增,引起大量的资源浪费和大面积的环境污染,甚至导致爆炸、中毒等恶性事故,对人们的生命和财产安全造成巨大的威胁。因此,为了有效减小或避免管道泄漏造成的危害,针对流体管道泄漏进行准确检测定位和及时的维修是十分必要的。
基于管道泄漏声振动信号的流体管道泄漏检测定位方法因其检测便捷方便、精度高而得到广泛应用,但泄漏声振动信号易受周围环境噪声的干扰导致信噪比较低,从而引起泄漏定位误差较大。流体管道泄漏声振动信号是一个多模态、非平稳的信号,针对非平稳的信号的分解和模态信号的提取一直是研究的热点方向,比如经验模态分解(empiricalmodaldecomposition,emd)、小波变换(wavelettransform,wt)、奇异值分解(singularvaluedecomposition,svd)、局域均值分解(localmeandecomposition,lmd)等,变分模态分解(variationalmodedecomposition,vmd)是dragomiretskiy等人在2014年提取出一种新型的信号分解处理方法(k.dragomiretskiyandd.zosso,“variationalmodedecomposition”,ieeetransactionsonsignalprocessing,vol.62,no.3,pp.531-544,2014),相比以上分解算法,vmd可消除模态混叠,具有良好的鲁棒性。公开号为cn108644618a的专利公布了一种基于vmd分量相对熵分析的管道泄漏定位方法,此方法将泄漏振动信号经过vmd分解得到模态分量,分别与原信号和无泄漏振动信号进行相对熵值比较,从而通过一定的阀值提取模态分量,能够去除干扰信号与虚拟分量,提升一定信噪比,但是vmd分解的模态层数确定不能适应不同状态泄漏的振动信号,其次不能实现自适应模态分量提取,同时由于环境噪声的复杂性,无泄漏振动信号与泄漏振动信号所处的时间不一样导致相对熵比较的效果较差。公开号为cn110454687a的专利公布了一种基于改进vmd的管道多点泄漏定位方法,该方法利用总体局域均值分解将泄漏振动信号分解为瞬时频率具有物理意义的乘积函数(productfunction,pf),根据相关系数提取出有效pf分量,对提取的pf分量进行信号重构,然后进行变分模态分解,其vmd的模态分解数是提取pf分量的个数,最后计算各模态分量的多尺度熵值并分解,筛选泄漏模态分量进行信号重构后进行互相关定位。该方法利用总体局域均值分解和pf分量的相关系数来确定vmd分解的模态分解数,实现了vmd的自适应分解,其次通过两次信号的分解与重构,能够消除了原始泄漏信号中的低相关成分和噪声的影响,但是在pf分量和模态分量的提取时不能实现信号的适应性提取。rehman等人在2019年提出改进的多元变分模态分解(multivariatevariationalmodedecomposition,mvmd)算法(n.u.rehmanandh.aftab,"multivariatevariationalmodedecomposition,"inieeetransactionsonsignalprocessing,vol.67,no.23,pp.6039-6052,1dec.1,2019.),该算法针对于多元信号之间存在有联合或共同的频率分量,同时继承vmd消除模态混叠和端点效应的优点,具有较好的鲁棒性和稳定性,但在进行管道泄漏振动信号分解时无法自适应确定模态分解层数。
技术实现要素:
有鉴于此,本发明的目的在于提供一种基于自适应多元变分模态分解的流体管道泄漏定位方法,根据泄漏信号自身特性自适应实现自适应消除噪声和自适应模态分量提取,从而提高信噪比,有效减少泄漏定位误差。
为达到上述目的,本发明提供如下技术方案:
一种基于自适应多元变分模态分解的流体管道泄漏定位方法,具体包括以下步骤:
s1:漏信号自适应分解:采用自适应多元变分模态分解(adaptivemultivariatevariationalmodedecomposition,amvmd)算法,根据两路泄漏信号的频率曲线相加后提取频率曲线的驼峰数,自适应确定多元变分模态的模态分解个数;
s2:模态分量自适应提取:根据两路泄漏信号的模态分量的信号清晰度,利用最大微分法自适应提取模态分量,从而提高泄漏信号信噪比;根据测量管道的长度和泄漏传播信号的速度确定泄漏信号的延时范围,选择提取两路泄漏信号模态分量的时间延迟在延时范围的模态分量作为有效模态分量,从而重构有效的泄漏信号;
s3:泄漏定位:利用两路有效泄漏信号的互相关延时估计确定时间延迟,结合泄漏信号传播速度和两传感器的距离确定泄漏点位置。
进一步,步骤s1中,所述amvmd算法具体包括:x(t)是包含c个维度的数据,x(t)=[x1(t),x2(t),…,xc(t),…,xc(t)],c∈1,2,…,c;
1)对x(t)的c个维度的数据相加得到x′(t):
将x′(t)进行傅里叶变换得到频谱曲线序列p,提取p的驼峰数k,驼峰数指的是频谱曲线纵坡在一定距离内较为急剧变化的凸峰个数,驼峰数确定x(t)的每个维度有共同的k个模态分量;则x(t)表示为:
其中,
2)为了得到每个维度的模态分量函数uk,c(t)的单边频谱,对uk,c(t)进行hilbert变换得到对应相关的解析信号,取其复数的实部;
3)混合各模态解析信号与预估中心频率
其中,{ωk}:={ω1,…,ωk}表示各个模态分量的中心频率,*为卷积;
4)通过公式(4)的范数l2平方来估计出每个维度的模态分量的模态信号带宽,对应的约束变分问题模型为:
其中,
5)引入二次惩罚因子α和拉格朗日乘法算子λc(t)求出式(5)的最优解,将约束性变分问题变为非约束性变分问题,即:
6)求公式(6)约束变分模型的最优解,将x(t)的每个维度的数据都分解k个模态分量函数imf。
进一步,步骤s1中,采用amvmd算法对泄漏信号进行模态分解的具体步骤为:
s11:对泄漏信号x1(t)的c个维度的数据相加得到x′(t),将x′(t)进行傅里叶变换的频谱曲线序列p,提取p的驼峰数k;
s12:初始化
s13:令n=n+1,对于所有ω≥0使得,根据公式(7)~(9)更新
其中,
s14:对于所有ω≥0,使得
s15:重复步骤s13、s14直到满足约束条件式(11),精度要求满足ε>0条件时结束循环;
s16:最终分解得到{uk}即为泄漏信号x1(t)的模态分量
s17:同理采用amvmd算法分解泄漏信号x2(t)得到模态分量
进一步,所述步骤s2具体包括以下步骤:
s21:对管道泄漏振动信号x1(t)的模态分量
则两泄漏信号的模态分量的时间延迟dk为:
为表示分解后的各个模态信号包含的泄漏信息量,定义参数pk为:
其中,t为模态信号的时域的长度,γ为指数因子;
对公式(14)中pk进行归一化,得到信号清晰度
s22:利用最大微分法提取x1(t)和x2(t)的模态分量
s23:结合模态分量互相时延估计dk,根据两传感器之间的距离和泄漏信号传播的速度计算泄漏信号的时延的范围,选取
更进一步,所述步骤s22中,最大微分法的提取步骤包括:
1)对信号清晰度
2)计算mk相邻的差值erri=mi+1-mi,1≤i≤k-1;
3)寻找err的最大值errmax,提取errmax前面对应x1(t)和x2(t)的模态分量得到
进一步,所述步骤s3具体包括以下步骤:
s31:对重构有效泄漏信号x′1(t)和x′2(t)进行互相关时延估计如下:
则两泄漏信号的时间延迟d为:
泄漏信号模态信号声速模型分别为:
流体主导模态声速为:
管壁主导模态声速为:
其中,vf为声波在流体介质中的传播速度,vl为声波在管壁材料中的传播速度,b为流体介质的体积弹性模量,a为管壁半径,h为管壁厚度,δ=ω2ρha2,ω为角频率,ρ为管壁材料的密度,ζ为管壁材料的泊松比,e为管壁材料的弹性模量;
s32:将重构后的有效泄漏信号相干函数的峰值频率作为中心频率,代入速度公式(18)(19),确定管壁主导波和流体主导波的传播速度;
流体管道泄漏流体主导波和管壁主导波在管壁径向引起的位移为w1和w2,其中二者之间的比例关系为:
其中,bf表示管内流体体积弹性模型,β=(2bfa/eh)(1-ν2),将确定的流体主导波和管壁主导波传播速度和材料参数代入公式(18)确定二者引起的管壁径向位移的比例,从而判断泄漏主导模态类型,从将主导模态的传播速度作为管道泄漏信号的传播速度v;
s33:将管道泄漏信号的传播速度v和两传感器之间的距离l代入泄漏定位公式(21),确定泄漏点位置;
其中,d为传感器1距离泄漏点的长度,l为两传感器之间的距离,该距离可通过现场测量直接获得。
本发明的有益效果在于:本发明利用流体管道泄漏振动信号的频率曲线特性提取频率曲线驼峰数,进而自适应确定多元变分模态分解算法的模态分解层数,实现泄漏振动信号自适应分解;其次,根据模态分量的信号清晰度和对应两路的模态分量互相关时间延迟范围,利用最大微分法和时延范围选取泄漏振动信号的有效模态分量,本发明技术能够根据泄漏信号自身特性自适应实现自适应消除噪声和自适应模态分量提取,从而提高信噪比,有效减少泄漏定位误差。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。
附图说明
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作优选的详细描述,其中:
图1为本发明基于自适应多元变分模态分解的流体管道泄漏检测定位方法流程图。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需要说明的是,以下实施例中所提供的图示仅以示意方式说明本发明的基本构想,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
请参阅图1,为一种基于自适应多元变分模态分解的流体管道泄漏声学定位方法,当流体管道发生泄漏时,泄漏源会释放压力形成振动信号沿着流体管道传播,此信号可被分解成多种不同类型的模态信号,因而具有多模态性。流体管道可以看作由管内流体、管壁及管道外部介质三层不同属性的介质构成的圆柱形壳体复合结构,低频泄漏振动信号在这种圆柱形壳体复合结构中传输时只激励两种基本的模态类型:管壁主导模态和流体主导模态,其中管壁主导模态波主要在管壁中传输,而流体主导模态波则为流体中传输的纵波。由于受泄漏源的自身的特性、管道材质、传输路径、管道周围环境噪声等多种复杂因素的影响,此信号是一个复杂的非平稳信号,同时包含大量环境噪声。当流体管道泄漏时产生振动信号沿着管道向两端传输,被两端加速度传感器拾取后形成两路管道泄漏信号分别为泄漏信号1为x1(t)和泄漏信号2为x2(t),然后泄漏信号进行模态分解、提取,重构、以及对管道泄漏定位,具体步骤如下:
1、泄漏信号自适应分解
本发明改进多元变分模态分解算法,在此基础上提出了amvmd算法,利用amvmd算法对管道泄漏信号x1(t)和x2(t)进行模态分解,惩罚参数α设为默认值2000,则x1(t)和x2(t)分解的结果如下;管道泄漏振动信号x1(t)的模态分量为
amvmd算法的分解原理如下:
x(t)是包含c个维度的数据,x(t)=[x1(t),x2(t),…,xc(t),…,xc(t)],c∈1,2,…,c。
1)对x(t)的c个维度的数据相加得到x′(t),
将x′(t)进行傅里叶变换得到频谱曲线序列p,提取p的驼峰数k,驼峰数指的是频谱曲线纵坡在一定距离内较为急剧变化的凸峰个数,驼峰数确定x(t)的每个维度有共同的k个模态分量。则x(t)可以表示为
其中,
2)为了得到每个维度的模态分量函数uk,c(t)的单边频谱,对其进行hilbert变换得到对应相关的解析信号(取其复数的实部):
3)混合各模态解析信号与预估中心频率
4)通过公式(4)的范数l2平方来估计出每个维度的模态分量的模态信号带宽,对应的约束变分问题模型如下:
其中,
5)引入二次惩罚因子α和拉格朗日乘法算子λc(t)求出式(5)的最优解,将约束性变分问题变为非约束性变分问题,即:
6)求公式(6)约束变分模型的最优解,将x(t)的每个维度的数据都分解k个模态分量函数imf。则利用amvmd算法对泄漏信号进行模态分解的具体实施过程如下:
(1)对泄漏信号x1(t)的c个维度的数据相加得到x′(t),将x′(t)进行傅里叶变换的频谱曲线序列p,提取p的驼峰数k。
(2)初始化
(3)令n=n+1,对于所有ω≥0使得,根据公式(7)、(8)、(9)更新
其中,
(4)对于所有ω≥0,使得
(5)重复步骤(3),(4)直到满足约束条件式(11),精度要求满足ε>0条件时结束循环。
(6)最终分解得到{uk}即为泄漏信号x1(t)的模态分量
(7)同理采用amvmd算法分解泄漏信号x2(t)得到模态分量
2、模态分量自适应提取
对管道泄漏振动信号x1(t)的模态分量
则两泄漏信号的模态分量的时间延迟dk为
为表示分解后的各个模态信号包含的泄漏信息量,定义参数pk为
其中,t为模态信号的时域的长度,γ为指数因子;
对公式(14)中pk进行归一化可得信号清晰度
然后利用最大微分法提取x1(t)和x2(t)的模态分量,最大微分法的提取步骤如下:
(1)对信号清晰度
(2)计算mk相邻的差值erri=mi+1-mi1≤i≤k-1;
(3)寻找err的最大值errmax,提取errmax前面对应x1(t)和x2(t)的模态分量得到
最后,结合模态分量互相时延估计dk,根据两传感器之间的距离和泄漏信号传播的速度计算泄漏信号的时延的范围,选取
3、泄漏定位
对重构有效泄漏信号x′1(t)和x′2(t)进行互相关时延估计如下:
则两泄漏信号的时间延迟d为
泄漏信号模态信号声速模型分别为:
流体主导模态声速为
管壁主导模态声速为
其中,vf为声波在流体介质中的传播速度,vl为声波在管壁材料中的传播速度,二者均可根据材料特性查表获得。b为流体介质的体积弹性模量,a为管壁半径,h为管壁厚度,δ=ω2ρha2,ω为角频率,为管壁材料的密度,ζ为管壁材料的泊松比,e为管壁材料的弹性模量,以上材料参数均可通过查阅材料数据手册获得。
然后,将重构后的有效泄漏信号相干函数的峰值频率作为中心频率,代入速度公式(18)(19)即可确定管壁主导波和流体主导波的传播速度。
流体管道泄漏流体主导波和管壁主导波在管壁径向引起的位移为w1和w2,其中二者之间的比例关系为
其中,vfd、vwd、bf分别表示流体主导模态声速、管壁主导模态声速、管内流体体积弹性模型,β=(2bfa/eh)(1-ν2),将确定的流体主导波和管壁主导波传播速度和材料参数代入公式(18)可确定二者引起的管壁径向位移的比例,从而判断泄漏主导模态类型,从将主导模态的传播速度作为管道泄漏信号的传播速度v。
最后,将管道泄漏信号的传播速度v和两传感器之间的距离l代入泄漏定位公式(21)即可确定泄漏点位置,
其中d为传感器1距离泄漏点的长度,l为两传感器之间的距离,该距离可通过现场测量直接获得。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。