基于时域积分的泄漏声波信号增强方法-林伟国.pdf
《基于时域积分的泄漏声波信号增强方法-林伟国.pdf》由会员分享,可在线阅读,更多相关《基于时域积分的泄漏声波信号增强方法-林伟国.pdf(8页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、2016 年 3 月 CIESC Journal 797 March 2016第 67卷 第 3期 化 工 学 报 Vol.67 No.3基于时域积分的泄漏声波信号增强方法 林伟国1,王晓东1,吴海燕1,牟长丽2,陈磊1(1北京化工大学信息科学与技术学院,北京 100029;2北京市射线应用研究中心,北京 100015) 摘要 :在管道泄漏检测中,泄漏信号的信噪比是影响误报和漏报的主要因素。传统的小波去噪和 EMD 分解方法受小波基选择、分解尺度选择和重构分量选择的影响难以保证稳定、有效的信号增强效果。通过分析动态压力变送器的传递函数,发现动态压力变送器的输出信号中缺失了反映声波信号低频响应特
2、性的积分项,提出了基于时域积分的、增益可调的泄漏声波信号增强方法,现场实测数据的测试结果表明:基于时域积分的泄漏信号增强方法具有稳定显著的信号增强效果和定位精度,无需复杂的参数寻优,为微小泄漏的检测以及减少误报和漏报提供了有效的技术支持。 关键词 :泄漏声波信号;时域积分;增益可调;信号增强 DOI: 10.11949/j.issn.0438-1157.20151304 中图分类号 : TE 973.6 文献标志码 : A 文章编号 : 0438 1157( 2016) 03 0797 08 Acoustic leak signal enhancement based on time-doma
3、in integration LIN Weiguo1, WANG Xiaodong1, WU Haiyan1, MU Changli2, CHEN Lei1(1College of Information Science and Technology, Beijing University of Chemical Technology, Beijing 100029, China; 2Beijing Radiation Application and Research Center, Beijing 100015, China) Abstract: Signal-to-Ratio (SNR)
4、of leak signal is the main factor influencing on missing alarm and false alarm in pipeline leak detection. Traditional de-noising methods, such as wavelet transform and EMD decomposition, are hard to promise stable and effective signal enhancement because of the selection of wavelet base, decomposit
5、ion scale or reconstruction components. Through analysis of the transfer function of dynamic pressure transducer, it is found that the integral part in the output signal of dynamic pressure transducer, which can reflect the low frequency response characteristics of acoustic signal, is missing. There
6、fore, a gain-adjustable acoustic signal enhancement method based on time-integration is proposed. The test results with field data show that the proposed method has better performance for the signal enhancement with precise location and no complex parameter optimization, and it also provides an effe
7、ctive technical support for small leak detection and decreasing of the missing and false alarm rate. Key words: acoustic leak signal; time-integration; gain-adjustable; signal enhancement 引 言 管道运输是化石能源和危化品输送的主要方式,由于管道老化、腐蚀穿孔及第三方施工破坏等原因,管道泄漏时有发生,造成资源损失、环境污染甚至人员伤亡。因此,管道泄漏检测技术发展至2015-08-13 收到初稿 , 2015-
8、10-13 收到修改稿。 联系人及第一作者 :林伟国( 1968) ,男,教授。 基金项目 :中央高校基本科研业务费( YS0104) ;国民核生化灾害防护国家重点实验室基金项目( SKLNBC2014-10) ;国家自然科学基金项目( 61403017) 。 Received date: 2015-08-13. Corresponding author: Prof. LIN Weiguo, Foundation item: supported by the Fundamental Research Funds for the Central Universities(YS0104), Fo
9、undation of State Key Laboratory of NBC Protection for Civilian(SKLNBC2014-10) and the National Natural Science Foundation of China (61403017). 化 工 学 报 第 67 卷 798 今,依然是国内外学者研究的热点问题之一。 在现有的负压波法、 声波法、 质量 /体积平衡法、分布式光纤法等多种管道泄漏检测方法中,由于声波传感器频率响应宽、灵敏度高、相关函数峰尖等特点1,在泄漏检测灵敏度和定位精度方面声波法是目前最受认可的泄漏检测方法2。 压电式动态压力传
10、感器作为声波传感器的一种3,在管道泄漏检测中得到了广泛应用4-5。不同于负压波法检测的是管道内部绝对压力 (静压) ,声波法检测的是声压的变化而非声压本身2,它反映的是瞬态压力的变化。对于快速发生的泄漏(如爆管) , 由于泄漏点产生的声压变化非常剧烈, 其泄漏检测灵敏度明显高于负压波法;但是对于微小泄漏和缓慢发生的泄漏 (如打孔盗油) , 由于泄漏本身引起的声压变化较小且泄漏声波在传播过程中存在的传播衰减6,使得首末站的声波传感器检测到的泄漏声波信号的信噪比较低,加剧了降低误报和杜绝漏报之间的矛盾。 针对微小泄漏的检测,文献 7-8提出了基于谐波小波分析技术的小泄漏信号识别和负压波拐点准确提取
11、方法,利用谐波小波时频图、时频等高线图及时频剖面图挖据管道泄漏敏感特征,准确提取负压波拐点。文献 9针对微弱负压波信号,应用局部投影将背景噪声和负压波特征信号分解到不同的子空间,通过子空间重构和小波包分析技术,提取弱特征信号,保留负压波拐点信息。上述方法都是针对负压波技术的, 并且方法的工程应用效果没有涉及。 针对泄漏声波信号,文献 10-11对实测信号作EMD 分解, 针对部分本征分量的叠加信号作希尔伯特黄变换,通过实时监测特征频带内信号的 Hilbert谱与正常信号相比是否发生明显变化,实现泄漏的诊断。文献 12利用泄漏声波信号的自相关函数具有的“不可重复”的特征,通过计算自相关函数序列的
12、近似熵,结合 Elman 神经网络,实现噪声淹没的泄漏声波信号的辨识。 微弱泄漏信号检测更直接的方法是信号去噪,从而提高泄漏信号信噪比。文献 13基于噪声白化准则控制检测信号的线性预测,提取出与干扰噪声相关、与泄漏信号无关的信号作为自适应噪声抵消器的参考输入,有效抵消噪声信号。文献 14借助EMD 分解,迭代计算不同 IMF 分量叠加信号的信噪比,实现重构信号信噪比的最大化。文献 15通过小波基的比较选择,实现信号的最优分解。上述这些方法都是通过抑制噪声提高泄漏信号的信噪比,是一种“被动型”的信号增强方法。众所周知,小波去噪是一种非自适应的信号增强方法,选择不同的小波基或滤波尺度可能会取得完全
13、相反的信号增强效果;而 EMD 方法不仅存在边界效应难以有效克服的问题, 而且还要解决选择哪些 IMF 分量重构信号的问题。事实上,由于金属管道及其输送介质具有的低频滤波特性1、泄漏信号传播到首末站时其频域分布与管道内部干扰信号的频域分布已经基本重叠,在抑制噪声的同时泄漏信号的有效成分势必也会受到损失,对于微小泄漏更是如此。 为解决上述问题,本文提出一种泄漏信号“主动”增强方法。首先从动态压力变送器的组成结构出发,通过分析动态压力变送器的传递函数,发现变送器输出信号中缺失了反映声波信号低频响应特性的积分项。基于此,本文提出基于数值积分的信号增强方法,通过在变送器输出端增加可调增益系数的积分项实
14、现对低频信号的放大。基于该方法得到的泄漏信号增强效果明显,信号增强效果不影响泄漏诊断和定位的准确性。 本文所提方法性能稳定,参数选择简单,可以避免传统去噪方法(如小波去噪)中难以进行参数选择的问题。 1 基于时域积分的泄漏信号增强原理 1.1 动态压力变送器传递函数分析 压电式动态压力变送器通常由前置放大、信号调理和信号变送 3 个部分组成16。压电式动态压力传感器利用压电材料(如石英、压电陶瓷等)的压电效应将被测压力转换为电信号,因此压电传感器可以等效为一电荷源。压电传感器输出阻抗极高且信号微弱,管道泄漏信号分布在较低频段,前置放大电路通常选用灵敏度高、低频响应好的电荷放大器。信号调理电路主
15、要完成信号的放大与滤波,电荷信号为微弱信号需要较高的增益,采用二级放大可以有效防止信号的畸变;由于压电式动态压力传感器有较好的高频响应,通常在信号中混有较多高频噪声,因此低通滤波器在变送器中必不可少。信号变送电路主要完成将电压转化为标准的 4 20 mA 电流信号,实现信号远传。动态压力变送器的组成结构如图 1 所示。 设电荷放大器传递函数为 H1(s), 低通滤波电路的传递函数为 H2(s), V/I 转换电路只实现信号类型的转换并不对信号时频域特性产生影响,因此可以与两级放大电路结合,其综合放大倍数设为 k。因 第 3 期 林伟国等:基于时域积分的泄漏声波信号增强方法 799图 1 动态压
16、力变送器组成结构图 Fig.1 Structure diagram of dynamic pressure transducer 此动态压力变送器的传递函数 H(s)可以表示为 12()()()YsHs kH sH sXs=() () (1) 式中, X(s)为变送器的动态压力输入信号 x(t)的拉普拉斯变换, Y(s)为变送器输出信号 y(t)的拉普拉斯变换。 由于泄漏声波信号主要是低频 (准静态)信号,高频成分为干扰信号,通常都是以通过低通滤波器的信号作为处理对象。因此式 (1)可以变换为 12Ys YsHs =kHsXs HsXs=() ()() ()() ()()(2) 式中, X(s
17、)为通过低通滤波器后的动态压力信号的拉普拉斯变换。此时系统的传递函数框图如图2 所示。 图 2 调整后的变送器传递函数框图 Fig.2 Transfer function block diagram after adjustment 电荷放大器的反馈电容 Cf与电阻 Rf决定了电荷放大器的时间常数 Tc,, Tc可以表示为 Tc=RfCf(3) 电荷放大器时间常数的影响近似于一个单端输入的高通滤波器,因此电荷放大器的传递函数可以近似为一阶 RC 高通滤波器16。电荷放大器的传递函数 H1(s)为 1()=1+ccsTHssT(4) 将式( 4)代入到式( 2)中,可以得到经过低通滤波后的动态压
18、力变送器的传递函数为 ()()= =() 1+ccskTYsHsXs sT(5) 经过拉普拉斯反变换可得 111( ) = ( )d + ( ) (0) + (0)t0cx tytytyxkT k k- (6) 由式( 6)可知,经过低通滤波后的动态压力信号 x(t)是由动态压力变送器输出信号 y(t)、 y(t)的积分与初始值 y(0)、 x(0)构成的。 x(t)中包含了动态压力信号的完整信息 (高频成分除外) , 但经过了动态压力变送器后(图 1) ,仅得到了输出信号 y(t)而损失了 y(t)的积分项,因此由 y(t)来表示 x(t)是不完整的。 由文献 17可知,泄漏声波信号的频率成
19、分主要集中在 1 Hz 以内,即使电荷放大器的时间常数Tc较大使得下限截止频率较低,仍会对泄漏信号产生衰减。对式( 6)中的积分项作傅里叶变换可得 1()d ()tFytt Fytjw=(7) 由式( 7)可见:由于低频段(尤其是低于 0.16 Hz 的低频分量)的 w 值较小,式( 7)中除法运算的效果相当于对低频分量进行了放大。因此,在动态压力变送器的输出信号 y(t)的基础上加上 y(t)的积分,有助于克服电荷放大器低频衰减的缺陷,提高泄漏信号的信噪比,实现泄漏信号的主动增强。 1.2 时域积分信号增强实现方法 基于时域积分的信号增强应在信号去噪(低通滤波)的基础上进行,实现步骤如下。
20、( 1)信号去噪,消除高频干扰; ( 2)低阶次的高通滤波,有效消除趋势18; ( 3)时域积分。 在实际处理过程中为进行时域积分,需将式( 6)离散化。本文选择梯形积分方法19并对式( 6)进行离散化有 =2(1) = (1)1()+(1)()=2()= ()+ () , =2,3, ,kiscxyyi yiSy kfTx k Syk yk k N_-(8) 式中, Sy(k)为式( 6)中积分项的梯形积分离散化形式, fs为采样频率, N 为信号样本的数据点数。 与小波去噪等其他信号增强方法不同的是,时域积分方法并不是消除干扰信号,而是对处于低频段的泄漏信号进行补偿, 以达到提高信噪比的效
21、果。在泄漏信号较弱的情况下,可以对输出信号的积分项 Sy(k)进行适当的放大,以实现泄漏信号的增强。为实现对信号增强效果的控制,本文引入低频增益系数 m(m 1),通过调节 m 的值实现对泄漏信号期望的增强效果,此时式( 8)表示为 x(k)=m Sy(k)+y(k) (9) 式中, k=1,2, ,N, Sy(1)=0。 2 基于时域积分的泄漏信号增强方法的现场应用 2.1 信号增强的有效性检验 信号增强的有效性检验主要比较基于小波去 化 工 学 报 第 67 卷 800 噪和基于时域积分的信号增强效果。信号增强效果通常采用信噪比作为评价指标。但由于管道泄漏声波信号为非周期信号,且没有固定的
22、频率点,噪声信号也不符合真正的高斯分布。为了评价信号增强效果,首先根据信号的过零点对信号作正负区间划分20。假设区间划分共得 M 个区间,泄漏信号所在区间的信号峰值 Peaks, 非泄漏信号的区间信号峰值为 PeakN(i),参照信噪比的定义公式,泄漏信号的信噪比可以表示为 22=1PeakSNR = 10lgPeak ( ) / ( )sMLNiiML-(10) 式中, L 为泄漏信号(区间)个数, i 为非泄漏信号的序号, i ML。以式( 10)定义的信噪比对时域积分方法的信号增强效果进行评判,结果如下文所示。 图 3 所示为采集自中石化某油田原油输送管道的模拟泄漏信号经过时域积分信号增
23、强前后的对比图(信号经过 db9 小波去噪,小波分解尺度 6) 。管道全长 26.68 km, 首站压力 1.2 MPa,末站压力 0.2 MPa,管道内径 500 mm, 模拟泄漏孔径 10 mm,模拟泄漏点距离管道首站 3.0 km, 信号采样精度为 12位,采样频率为 100 Hz。由图 3( b)可见,下游泄漏信号并不突出,其几乎被噪声淹没。根据式( 9)经过时域梯形积分信号增强后(时域积分的增益系数 m=1) ,泄漏信号得到了明显的增强 图 3( d) 。根据式( 10)分别计算得到信号增强前后上、下游信号的信号区间个数 M、异常信号个数 L 及泄漏信号的信噪比 SNR 对比如表 1
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 时域 积分 泄漏 声波 信号 增强 方法 林伟国
限制150内