一种基于LFM信号互相关相位特性的超声波测厚方法技术

本发明专利技术公开了一种基于LFM信号互相关相位特性的超声波测厚方法，该方法是采用LFM(线性调频)信号作为超声波测厚的发射信号，在接收到回波信号后对回波信号进行匹配滤波，基于LFM信号互相关相位特性估计LFM声信号在油气输送管道管壁内外表面之间的传播时间，最后结合预置的声速求得管壁厚度。本发明专利技术对比现有技术，有测量精度高、采样率要求低、发射功率要求低、抗噪声干扰能力强、多径容忍度高等特点。多径容忍度高等特点。多径容忍度高等特点。

【技术实现步骤摘要】
一种基于LFM信号互相关相位特性的超声波测厚方法


[0001]本专利技术涉及超声波测厚
，具体涉及一种基于LFM信号互相关相位特性的超声波测厚方法。

技术介绍

[0002]在石油化工行业，需要对油气输送管道进行定期的检测，以防止管道因外界电化学等作用而遭受腐蚀。管道因腐蚀而产生破损，具有一定的危险性，轻则产生油气泄漏而污染环境，重则产生爆炸而危害人们的生命财产安全。在针对管道腐蚀问题通常采用超声波测厚技术对管壁厚度进行定期检测，根据管壁厚度的变化程度，判定油气输送管道的腐蚀程度。目前，超声波测厚技术已经在长输管线和海上管道得到广泛应用。超声波测厚技术是一种主动式无损检测技术，具有定位精度较高、检测成本较低、无需拆除防腐层，可实现管壁100％检测等特点。
[0003]超声测厚技术按测量原理不同可分为共振法、兰姆波法和脉冲发射法三种。其中的共振法以及兰姆波法都对待测物体表面光洁度要求较高，不适用于对油气输送管道这类带漆面且表面可能存在腐蚀现象的材料进行检测。脉冲发射法的原理是通过估计短脉冲声信号在油气输送管道管壁内外表面之间往返一次的传播时间来求得工件的厚度。相对于前两种方法而言，脉冲反射法原理简单，实现方便，对管道的表面光洁度要求不高，可测粗糙表面、凹面以及带漆面材料，适用范围广，是现在应用最为广泛的一种方法。
[0004]但是，这类采用短脉冲信号作为发射信号的超声测厚技术通常具有发射信号持续时间较短、占空比较低等特点。根据数字信号处理理论，基于这类采用短脉冲信号作为发射信号的超声测厚技术有以下缺陷：
[0005](1)测量精度低，根据香农信息论和参数估计理论，单一频率的短脉冲信号的时宽—带宽积小，回波能量小，回波检测的精度低。目前常用的超声波测厚装置的测量误差一般为5％，测量精度有待提高。
[0006](2)探测深度短，脉冲持续时间越短，所携带的能量就越低，易受到信号幅度衰减的影响，从而导致较短的探测深度，不利于测量深埋管道的管壁厚度。
[0007](3)对噪声敏感，对测量环境有一定的要求，因为采用短脉冲信号作为发射信号的超声测厚技术对声信号在管壁内的传播时延估计是基于声压(即幅度)特征的，在信噪比较低的环境下，容易受噪声干扰。
[0008](4)多径容忍度低，单一频率的短脉冲信号受不同回波路径而产生的多径干扰影响较大，抗干扰能力弱，易造成回波信号的误检或漏检现象。

技术实现思路

[0009]本专利技术目的在于克服采用短脉冲信号作为发射信号的超声测厚技术的不足，提出一种基于LFM信号互相关相位特性的超声波测厚方法。该方法是采用LFM(线性调频)信号作为超声波测厚方法的发射信号，在接收到回波信号后对回波信号进行匹配滤波，基于LFM信
号互相关相位特性估计LFM声信号在油气输送管道管壁内外表面之间的传播时间，最后结合预置的声速求得管壁厚度。
[0010]本专利技术的目的可以通过采取如下技术方案达到：
[0011]一种基于LFM信号互相关相位特性的超声波测厚方法，应用于由数据处理单元、D/A转换模块、功率放大器、发射超声探头、接收超声探头、低噪声放大器、A/D转换模块组成的超声波测厚装置。其中，数据处理单元与D/A转换模块、功率放大器、发射超声探头依次相连组成发射信号通道。接收超声探头、低噪声放大器、A/D转换模块和数据处理单元依次相连组成接收信号通道。数据处理单元完成超声波测厚装置中的多项工作，包括：产生LFM信号作为发射信号；对回波信号进行匹配滤波，得到回波信号和发射信号的互相关函数；基于LFM信号互相关相位特性估计LFM声信号在管壁内外表面之间的传播时间；结合预置的声速求得油气输送管道管壁厚度。D/A转换模块将数据处理单元产生的数字信号转换为模拟电信号。功率放大器将模拟电信号的功率放大。发射超声探头将放大的模拟电信号转化成LFM声信号发射至油气输送管道管壁内部。接收超声探头接收从石油管道管壁内外表面反射的LFM声信号并将其转化为模拟电信号。低噪声放大器对模拟电信号进行放大。A/D转换模块将模拟电信号转化为数字信号并作为回波信号发送到数据处理单元。
[0012]该超声波测厚方法包括以下具体步骤：
[0013]S1、利用线性调频的调制方式产生LFM信号作为发射信号x[n]，x[n]的表达式为：
[0014][0015]其中，A为发射信号的幅度，f0与f1分别是发射信号的起始频率与终止频率，T为信号时长，F
S
为采样率，μ＝(f
1-f0)/(F
s2
T)，μ为调频斜率，N＝F
S
T，N为样本总数，发射信号x[n]的频率随着时间的变化在一定范围内呈线性增长；
[0016]S2、发射信号x[n]经过发射信号通道转化为LFM声信号并发射至待测油气输送管道的管壁；
[0017]S3、从油气输送管道管壁的内、外表面反射回的回波声信号经过接收信号通道转化为回波信号y[n]，回波信号y[n]是由两重具有不同时延的LFM信号叠加而成，设两重具有不同时延的LFM信号的衰减幅度分别为a1、a2，时延估计值分别为t1、t2，对应的分数阶坐标其中T
S
为采样间隔，得到y[n]的表达式为：
[0018][0019]S4、对回波信号y[n]进行匹配滤波，得到发射信号x[n]与回波信号y[n]的互相关函数Q(n)，Q(n)的表达式为：
[0020][0021]S5、搜索互相关函数Q(n)的第一主峰值点，得到其对应的离散坐标并计算得到粗时延估计值表达式为：
[0022][0023]其中，T
S
为采样间隔，
[0024]S6、引入分数阶估计参数α，得到精时延估计值t1与粗时延估计值关系式如下：
[0025][0026]其中，t
α
为残留时延，t
α
＝αT
S
，
[0027]构建LFM信号与其单一时延信号做互相关得到的一般互相关函数R(n,n1)表达式：
[0028][0029]将一般互相关函数R(n,n1)坐标左移使得峰值点位于零点，再取左移后的一般互相关函数R(n,n1)的相位部分，构建相位函数)的相位部分，构建相位函数表达式如下：
[0030][0031]同样对匹配滤波得到的发射信号和回波信号的互相关函数Q(n)进行左移并取相位部分，得到实际相位函数Q2(n)，Q2(n)表达式如下：
[0032][0033]构建代价函数J：
[0034][0035]Δn取值1，2，
…
，10，令得到对应的残留时延t
α
，为简化表达令：
[0036][0037][0038][0039][0040]得到残留时延t
α
的表达式为：
[0041][0042]最终得到精时延估计值
[0043]S7、重复步骤S5和S6，搜索互相关函数Q(n)的第二主峰值点，得到精时延估计值t2；...

【技术保护点】

【技术特征摘要】
1.一种基于LFM信号互相关相位特性的超声波测厚方法，其特征在于，所述的超声波测厚方法包括以下具体步骤：S1、利用线性调频的调制方式产生LFM信号作为发射信号x[n]，x[n]的表达式为：其中，A为发射信号的幅度，f0与f1分别是发射信号的起始频率与终止频率，T为信号时长，F
S
为采样率，μ为调频斜率，N＝F
S
T，N为样本总数，发射信号x[n]的频率随着时间的变化在一定范围内呈线性增长；S2、发射信号x[n]经过发射信号通道转化为LFM声信号并发射至待测油气输送管道的管壁；S3、从油气输送管道管壁的内、外表面反射回的回波声信号经过接收信号通道转化为回波信号y[n]，回波信号y[n]是由两重具有不同时延的LFM信号叠加而成，设两重具有不同时延的LFM信号的衰减幅度分别为a1、a2，时延估计值分别为t1、t2，对应的分数阶坐标其中T
S
为采样间隔，得到y[n]的表达式为：S4、对回波信号y[n]进行匹配滤波，得到发射信号x[n]与回波信号y[n]的互相关函数Q(n)，Q(n)的表达式为：S5、搜索互相关函数Q(n)的第一主峰值点，得到其对应的离散坐标并计算得到粗时延估计值表达式为：其中，T
S
为采样间隔，S6、引入分数阶估计参数α，得到精时延估计值t1与粗时延估计值t
1*
关系式如下：其中，t
α
...

【专利技术属性】
技术研发人员：杨萃，孙家豪，余华，宁更新，陈芳炯，季飞，
申请(专利权)人：华南理工大学，
类型：发明
国别省市：