管道导波无损检测回波信号处理方法研究
- 格式:pdf
- 大小:633.81 KB
- 文档页数:5
管道导波无损检测回波信号处理方法研究王悦民,曾 帆,孙丰瑞(海军工程大学船舶与动力学院,武汉430033)摘 要:简述了导波在管道中传播的频散特性,给出了管道导波无损检测实验信号特征。
对管道导波检测的回波信号首先进行Hilbert 变换,然后对变换后的信号分别进行加窗函数滤波和小波变换,处理后的缺陷信号明显易于辨别。
经过分析比较,结果表明了这两种组合方法用于管道导波无损检测信号处理的效果基本一样,在实际检测应用中,选取其中任一组合方法都是可行的。
关键词:导波;Hilbert 变换;Kaiser 窗;小波变换中图分类号:TB553 文献标志码:A 文章编号:1009-3486(2008)05-0018-05Study on a method of echo signal processing for guidedw aves nondestructive testing of pipesWAN G Yue 2min ,ZEN G Fan ,SUN Feng 2rui(College of Naval Architect ure and Power ,Naval Univ.of Engineering ,Wuhan 430033,China )Abstract :The dispersion characteristic of guided wave p ropagation in pipes was described and t he ex 2perimental signals feat ure of guided waves ND T of pipes was given.The Hilbert t ransform met hod was used for p rocessing t he echo signals detected by guided waves in pipe ,and t hen t he transformed signals were p rocessed respectively wit h digital window f unction filter and wavelet t ransform me 2t hods.The defect signals could be easily distinguished after being processed.The analysis and com 2pare result s show t hat t he effect by using two kinds of combined met hods to process t he signals of guided waves ND T of pipes are nearly t he same.So selecting eit her of t he combined signal p rocessing met hods is feasible in t he p ractical testing.K ey w ords :guided wave ;Hilbert t ransform ;Kaiser window ;wavelet t ransform利用超声导波对管道进行无损检测时,检测的信号会受到多个因素影响,如被检测管道表面或端头的反射、空间电磁干扰等,若不采取相应的措施,这些干扰信号甚至会淹没有用信号,检测效果不理想。
因此,除采取电磁屏蔽措施外,利用有效的信号处理方法对增强目标信号、滤除噪声是非常必要的。
超声导波在管道中传播时具有频散特性和多模式特征,虽然通过计算导波在管道中传播的频散曲线、选择合适的激励频率可以减弱频散现象,但对于任何给定的激励频率在管道中至少可产生两种模态的导波,并且随着频厚积的增加,产生的导波模态还会增加。
另外,导波在管道中传播时,遇到不连续处或者随传播距离的增加,波包都会发生改变,信号幅值也将减小等,这些因素导致管端和缺陷处产生的回波信号非常复杂。
目前,导波信号处理的方法主要有傅立叶变换、时频分析、小波变换、希尔波特变换等。
Alleyne D N 等[1]利用二维傅立叶变换将多种模式导波信号的重叠部分分离出来。
Kwun H 等[2] 第20卷 第5期 2008年10月 海军工程大学学报 J OURNAL OF NAVAL UN IV ERSIT Y OF EN GIN EERIN G Vol.20 No.5 Oct.2008 3收稿日期:2008204208;修回日期:2008206230。
基金项目:中国博士后科学基金资助项目(20060390865);湖北省自然科学基金资助项目(2007ABA224)。
作者简介:王悦民(1969-),男,副教授,博士,主要研究方向为舰船动力及热力系统的监测、控制与故障诊断,E 2mail :yuemin_wang @ 。
利用时频分析的方法———短时傅立叶变换对圆管中纵向导波发散特性进行了分析,指出纵波存在“低通”或“高通”截止频率特性。
Siqueira M H S [3]等用小波分析方法对管道导波检测信号进行分析,指出经小波变换后信噪比得到改善,其值可达12dB 。
针对管道导波检测回波信号的特点,王悦民[4]选用Hilbert 变换对回波信号进行了处理,结果表明,Hilbert 变换法对导波信号的分辨比较有效。
吴斌等[5]综述了超声导波无损检测中的信号处理研究进展,指出信号处理在提高信号的可读性、提高缺陷的检测与识别能力、表征被检对象特性等方面具有重要作用。
文中首先简述了管道导波频散特性和实验检测信号特征,然后利用Hilbert 变换对检测信号进行处理,最后分别采用加窗函数滤波和小波变换方法对Hilbert 变换后的信号作进一步处理,分析了处理的效果。
1 管道导波频散特性和实验检测信号特征1.1 管道导波频散特性导波频散(dispersion )是指导波的传播速度随着频率、构件几何边界形状的变化而变化,以及导波会以多个模式同时传播,当遇到构件不连续处还会发生模式转换。
图1给出了长为6.72m ,外径为38mm ,壁厚为2.5mm 钢管的导波频散曲线。
为了标记导波在管道中存在的不同模态,图1中分别用L (0,n ),T (0,n )和F (m ,n )表示对应的纵向、扭转和弯曲模态导波,其中:m 表示导波沿圆管周向的位移和应力谐变量阶数;n 表示模数。
从图中频散曲线可看出,在频率小于100k Hz 的范围内,同时有L (0,1),L (0,2),T (0,1),F (1,1),F (1,2),F (1,3)模式导波出现;对于任一给定频率,至少存在两个以上的导波模态,随着频率的增大导波模态也迅速增多。
从图1(a )群速度频散曲线中,还可以看出L (0,1)和L (0,2)模式导波分别有“低通”和“高通”截止频率出现。
图1 钢管导波频散曲线Fig.1 Dispersion curves of guided wave in steel pipe在理想的情况下,希望在管道中仅产生单一模态的导波,若导波遇到缺陷产生反射,所得到的回波将较容易分析。
多模态的存在大大增加了导波技术应用于管道检测的复杂性。
在实际检测应用中,对于给定的被检测管道,选择适当的导波激励频率非常重要,这就要求选择的激励频率不在截止频率区。
另外,为了提高检测的灵敏度或信噪比,尽可能抑制其他模式导波的产生,从而激励出单一模式的导波,并且使激励的导波频率在非频散区域内。
1.2 实验检测信号特征取长为6.72m 、外径为38mm 、壁厚为2.5mm 的钢管为实验样本,在钢管上钻3个通孔,模拟通孔缺陷。
孔大小分别为<12mm 、<10mm 和<5mm ,其中,孔<12mm 位于钢管两端的中间处,而孔<10mm 和<5mm 分别位于钢管的三等分处,即与钢管两端的距离均为2.24m (见图2)。
该实验装置设置・91・ 第5期 王悦民等:管道导波无损检测回波信号处理方法研究 的具体内容可参考文献[6],选用频率为33kHz 的纵向导波激励。
检测信号波形如图3所示。
从图3可看出,由于受空间磁场干扰以及接收传感器自身的噪声影响,接收的时域信号虽然能够反映出缺陷信号特征,但噪声信号很大,为了使缺陷信号更易于分辨并实现缺陷定量化研究,必须对检测信号作数字处理。
2 检测信号的H ilbert 变换基于磁致伸缩效应的管道无损检测技术得到的磁致伸缩脉冲回波信号是一种窄带信号,包含管道状态的重要信息,必须对这些窄带信号进行分析处理,以从中提取有用的信息。
在信号处理中,Hilbert 变换是获得窄带信号包络的有效方法。
利Hilbert 变换,可以构造出相应的解析信号,使其仅含正频率成分,从而可降低信号的抽样率。
还可以应用解析表达式中实部与虚部之间的关系,定义出任意时刻的瞬时包络、瞬时相位及瞬时频率,从而解决短信号和复杂信号瞬时特征的提取问题。
对于一实函数,f (t )的Hilbert 变换定义为f ∧(t )=H [f (t )]=1π∫+∞-∞f (τ)t -τd τ=f (t )31πt 。
(1)式中:H 为Hilbert 变换;f ∧(t )为f (t )的Hilbert 变换;3为卷积。
Hilbert 变换相当于通过一个冲激响应为1πt的线性网络,该网络也称为Hilbert 滤波器。
f (t )与其Hilbert 变换f ∧(t )正交。
图4 Hilbert 变换后的信号波形 Fig.4 Signals processed by Hilbert transform对于窄带信号x (t )=a (t )co s θ(t ),如果引入其正交分量x ∧(t )=a (t )sin θ(t ),可以将它们组成一个复信号:u (t )=a (t )co s θ(t )+ja (t )sin θ(t )=a (t )e j θ(t )。
(2)这样信号x (t )的瞬时包络a (t )、相位φ(t )、频率ω(t )可以分别表示如下:a (t )=x 2(t )+x ∧2(t ),(3)φ(t )=arctan x ∧(t )x (t ),(4)ω(t )=12π・d φ(t )d t。
(5) 对图3中有3个孔缺陷的钢管检测信号进行Hilbert 变换,变换后的信号如图4所示。
从信号变换结果可以看出,经过Hilbert 变换后的钢管缺陷信号分辨更清楚,而信号的幅值、相位和频率保持相对・02・海 军 工 程 大 学 学 报 第20卷 不变。
因此,在时域中对管道磁致伸缩脉冲回波信号的分析和识别,Hilbert 变换是一种非常简便和实用的方法。