CN110243945A - 基于合成孔径聚焦与模式转换波的超声tofd盲区抑制方法 - Google Patents
基于合成孔径聚焦与模式转换波的超声tofd盲区抑制方法 Download PDFInfo
- Publication number
- CN110243945A CN110243945A CN201910669559.6A CN201910669559A CN110243945A CN 110243945 A CN110243945 A CN 110243945A CN 201910669559 A CN201910669559 A CN 201910669559A CN 110243945 A CN110243945 A CN 110243945A
- Authority
- CN
- China
- Prior art keywords
- tofd
- wave
- point
- scanning
- blind area
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0654—Imaging
- G01N29/069—Defect imaging, localisation and sizing using, e.g. time of flight diffraction [TOFD], synthetic aperture focusing technique [SAFT], Amplituden-Laufzeit-Ortskurven [ALOK] technique
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/023—Solids
Landscapes
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Signal Processing (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
一种基于合成孔径聚焦与模式转换波的超声TOFD盲区抑制方法,属于无损检测技术领域。该方法采用由超声探伤仪、TOFD探头、倾斜有机玻璃楔块和扫查装置组成的TOFD检测系统,沿待测工件表面实施扫查与图像采集,获得不同扫查位置处的A扫描信号集合。根据费马定理求解各位置处的模式转换波最短传播声时,并对所有A扫描信号进行时间延迟和幅值叠加处理,通过逐点聚焦实现TOFD扫查图像重建,削弱横向冗余信号,改善成像质量。利用重建图像中模式转换波显示深度与缺陷端点深度之间的相关性,计算盲区范围内缺陷的端点深度。与现有的模式转换波方法相比,该方法能够在抑制TOFD盲区的同时,提高缺陷定量精度,具有较高的工程应用和推广价值。
Description
技术领域
本发明涉及一种基于合成孔径聚焦与模式转换波的超声TOFD盲区抑制方法,其属于无损检测技术领域。
背景技术
超声衍射时差法(Time of Flight Diffraction,TOFD)是一种利用缺陷端点衍射纵波时间差实施缺陷定量的无损检测方法,具有定量精度高、缺陷检出率高等优点,被广泛应用于石油化工、电力及天然气等领域的各向同性材料检测。然而,对近表面缺陷进行超声TOFD扫查时,由于超声波脉冲具有一定宽度,使得缺陷端点衍射波信号容易湮没在直通波中,形成近表面检测盲区,难以对盲区内缺陷进行深度定量。
为解决上述问题,研究人员根据横波传播声时大于纵波的特性,提出利用缺陷端点衍射横波,即模式转换波实施TOFD盲区内缺陷定量检测。然而,声束扩散导致不同扫查位置均能接收到模式转换波,使得扫查图像中存在横向冗余信号,干扰模式转换波识别,降低缺陷定量精度。本发明提出了一种基于合成孔径聚焦与模式转换波的超声TOFD盲区抑制方法,可削弱扫查图像中模式转换波横向冗余信号,改善成像质量,在减小TOFD盲区范围的同时,提高缺陷定量精度。
发明内容
本发明提供一种基于合成孔径聚焦与模式转换波的超声TOFD盲区抑制方法。其目的是针对TOFD盲区内缺陷扫查图像中模式转换波存在横向冗余信号的问题,结合合成孔径聚焦技术对模式转换波实施逐点聚焦成像,并利用重建图像中模式转换波显示深度与缺陷端点深度的相关性,实现TOFD盲区内缺陷端点深度的精确定量。
本发明采用的技术方案是:一种基于合成孔径聚焦与模式转换波的超声TOFD盲区抑制方法,针对TOFD盲区内缺陷扫查图像中模式转换波存在横向冗余信号的问题,结合合成孔径聚焦技术对模式转换波实施逐点聚焦成像,并利用重建图像中模式转换波显示深度与缺陷端点深度的相关性,实现TOFD盲区内缺陷的深度定量;所述方法采用如下步骤:
(a)TOFD检测参数确定
根据待测工件的材料、几何尺寸及检测范围选取TOFD检测参数,包括TOFD探头频率、楔块斜楔角、探头中心间距、检测增益、采样频率、电子扫查步进;
(b)扫查数据采集
采用步骤(a)中确定的TOFD检测参数,控制TOFD探头沿工件表面进行扫查,扫查步长为ΔS,步进次数为N-1,共获得N个A扫描信号构成的数据集;
(c)坐标系建立及被检测区域网格划分
以待测工件表面为X轴,扫查方向为X轴正向,工件深度方向为Y轴正向建立坐标系,并设扫查开始时发射探头声束出射点位置为坐标原点;将被检区域划分为矩形网格,各网格节点作为图像重建点;任意图像重建点的横、纵坐标分别定义为x和y;
(d)模式转换波延迟时间计算
设TOFD探头中心间距为2S,发射探头的声束出射点到任意图像重建点(x,y)的距离为L1,该图像重建点到接收探头的距离为L2;根据TOFD探头中心间距与图像重建点(x,y)的几何关系可知,当进行到第k次步进时,0≤k≤N-1,L1和L2分别为
根据费马定理知,声波沿传播时间最短的路径传播;设cL和cS分别为纵波和横波声速,则第k次步进时,TOFD探头入射纵波在图像重建点产生的散射横波,即模式转换波的传播时间tk为:
进一步地,定义tk中的最小值为tmin,0≤k≤N-1,则模式转换波传播时间相对于tmin的时间延迟dtk为:
dtk=tk-tmin (4)
(e)合成孔径聚焦图像重建
针对图像重建点x,y,对不同扫查位置处采集的A扫描信号施加时间延迟dtk并进行信号叠加处理,得到该点的信号幅值I(x,y)为:
式中,fk()为第k次步进时的探头接收信号,t为时间;
对被检区域内每个网格点重复上述过程,得到合成孔径聚焦后的重建图像,图中模式转换波的横向冗余信号得到削弱;
(f)TOFD盲区内缺陷端点深度定量
设待测工件内部缺陷位于TOFD盲区内,其端点深度为d;从重建图像中读取两侧模式转换波幅值最高点间距L,以及幅值最高点显示深度d′,则利用式(6)反演得到缺陷端点深度d;
本发明的有益效果是:这种基于合成孔径聚焦与模式转换波的超声TOFD盲区抑制方法,结合合成孔径聚焦技术削弱模式转换波的横向冗余信号,改善成像质量;利用重建图像中模式转换波显示深度与缺陷端点深度之间的相关性,计算TOFD盲区范围内缺陷的端点深度。与现有的模式转换波方法相比,该方法能够在抑制TOFD盲区的同时,提高缺陷定量精度,具有较高的工程应用和推广价值。
附图说明
下面结合附图和实例对本发明做进一步说明。
图1是采用的TOFD检测系统示意图。
图2是加工底面开口槽的碳钢平板试块图。
图3是原始的TOFD非平行扫查图像。
图4是合成孔径聚焦处理后的重建图像。
具体实施方式
基于合成孔径聚焦与模式转换波的超声TOFD盲区抑制方法,采用的超声检测系统如图1所示,其中包括超声探伤仪、一对TOFD探头、一对纵波角度楔块和扫查装置。具体检测及处理步骤如下:
(a)试验对象如图2所示碳钢平板试块,其壁厚为30.0mm,材料纵波声速为5890m/s,横波声速为3230m/s,试块内部加工了距表面3.0mm的底面开口槽。
(b)采用标称频率5MHz的TOFD探头实施检测。主要检测参数包括:纵波楔块角度70°、探头中心间距2S=30.0mm、采样频率100MHz、检测增益60dB、扫查步进ΔS=0.40mm、A扫描时间窗口起始位置设置为直通波到达接收探头之前,则计算可得TOFD近表面盲区理论深度约为6.1mm。
(c)采用步骤(b)中确定的TOFD检测参数,利用TOFD探头对碳钢平板试块中的底面开口槽进行非平行扫查,获得164个A扫描信号构成的数据集,并通过探伤仪A/D转换器以.txt文件形式导出。图3给出对应的非平行扫查图像,显然底面开口槽端点衍射波与直通波混叠,且模式转换波存在横向冗余信号。此时,直接从原始图像中读取两侧模式转换波幅值最高点间距L为28.2mm,幅值最高点显示深度d′为9.76mm,代入式(6)反演计算得到底面开口槽端点深度为2.64mm,绝对误差为0.36mm。
(d)根据费马定理,计算任意图像重建点对应不同扫查位置处的模式转换波最短传播声时。基于合成孔径聚焦成像原理,逐点实施时间延迟和幅值叠加处理,得到各点的合成信号,重建TOFD扫查图像,如图4所示。显然,经过合成孔径聚焦处理后,模式转换波横向冗余信号得到削弱。
(e)如图4所示,从重建图像中读取两侧模式转换波幅值最高点间距L为27.3mm,幅值最高点显示深度d′为9.95mm,代入式(6)反演计算得到底面开口槽端点深度为2.79mm,绝对误差为0.21mm。显然,利用该方法可将TOFD盲区从6.1mm抑制到3.0mm以内,且缺陷深度定量相对误差从处理前的12.0%下降至7.0%。
Claims (1)
1.一种基于合成孔径聚焦与模式转换波的超声TOFD盲区抑制方法,针对TOFD盲区内缺陷扫查图像中模式转换波存在横向冗余信号的问题,结合合成孔径聚焦技术对模式转换波实施逐点聚焦成像,并利用重建图像中模式转换波显示深度与缺陷端点深度的相关性,实现TOFD盲区内缺陷的深度定量;其特征在于,所述方法采用如下步骤:
(a)TOFD检测参数确定
根据待测工件的材料、几何尺寸及检测范围选取TOFD检测参数,包括TOFD探头频率、楔块斜楔角、探头中心间距、检测增益、采样频率、电子扫查步进;
(b)扫查数据采集
采用步骤(a)中确定的TOFD检测参数,控制TOFD探头沿工件表面进行扫查,扫查步长为ΔS,步进次数为N-1,共获得N个A扫描信号构成的数据集;
(c)坐标系建立及被检测区域网格划分
以待测工件表面为X轴,扫查方向为X轴正向,工件深度方向为Y轴正向建立坐标系,并设扫查开始时发射探头声束出射点位置为坐标原点;将被检区域划分为矩形网格,各网格节点作为图像重建点;任意图像重建点的横、纵坐标分别定义为x和y;
(d)模式转换波延迟时间计算
设TOFD探头中心间距为2S,发射探头的声束出射点到任意图像重建点(x,y)的距离为L1,该图像重建点到接收探头的距离为L2;根据TOFD探头中心间距与图像重建点(x,y)的几何关系可知,当进行到第k次步进时,0≤k≤N-1,L1和L2分别为
根据费马定理知,声波沿传播时间最短的路径传播;设cL和cS分别为纵波和横波声速,则第k次步进时,TOFD探头入射纵波在图像重建点产生的散射横波,即模式转换波的传播时间tk为:
进一步地,定义tk中的最小值为tmin,0≤k≤N-1,则模式转换波传播时间相对于tmin的时间延迟dtk为:
dtk=tk-tmin (4)
(e)合成孔径聚焦图像重建
针对图像重建点x,y,对不同扫查位置处采集的A扫描信号施加时间延迟dtk并进行信号叠加处理,得到该点的信号幅值I(x,y)为:
式中,fk()为第k次步进时的探头接收信号,t为时间;
对被检区域内每个网格点重复上述过程,得到合成孔径聚焦后的重建图像,图中模式转换波的横向冗余信号得到削弱;
(f)TOFD盲区内缺陷端点深度定量
设待测工件内部缺陷位于TOFD盲区内,其端点深度为d;从重建图像中读取两侧模式转换波幅值最高点间距L,以及幅值最高点显示深度d′,则利用式(6)反演得到缺陷端点深度d;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910669559.6A CN110243945B (zh) | 2019-07-23 | 2019-07-23 | 基于合成孔径聚焦与模式转换波的超声tofd盲区抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910669559.6A CN110243945B (zh) | 2019-07-23 | 2019-07-23 | 基于合成孔径聚焦与模式转换波的超声tofd盲区抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110243945A true CN110243945A (zh) | 2019-09-17 |
CN110243945B CN110243945B (zh) | 2021-04-30 |
Family
ID=67893348
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910669559.6A Active CN110243945B (zh) | 2019-07-23 | 2019-07-23 | 基于合成孔径聚焦与模式转换波的超声tofd盲区抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110243945B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113358751A (zh) * | 2021-06-01 | 2021-09-07 | 中车青岛四方机车车辆股份有限公司 | 工件缺陷的检测方法、装置和系统 |
CN113484421A (zh) * | 2021-07-01 | 2021-10-08 | 中国工程物理研究院机械制造工艺研究所 | 基于合成孔径的激光超声内部缺陷多模成像方法及系统 |
CN113758617A (zh) * | 2021-09-10 | 2021-12-07 | 哈尔滨工业大学 | 一种基于宽带扫频信号频域计算的应力梯度高效无损检测系统及其检测方法 |
CN114674928A (zh) * | 2022-03-21 | 2022-06-28 | 大连理工大学 | 基于二阶模式波的薄板结构超声tofd检测盲区抑制方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102207488A (zh) * | 2011-03-29 | 2011-10-05 | 北京理工大学 | 横波tofd缺陷定位方法 |
CN102393422A (zh) * | 2011-08-22 | 2012-03-28 | 江苏省产品质量监督检验研究院 | 基于超声tofd的缺陷离线判别方法 |
CN103149274A (zh) * | 2013-01-28 | 2013-06-12 | 中国科学院声学研究所 | 一种混凝土缺陷检测方法 |
CN104458912A (zh) * | 2014-12-17 | 2015-03-25 | 国家电网公司 | 瓷绝缘子用小角度纵波斜探头入射点入射角及延时测量法 |
CN105388218A (zh) * | 2015-12-10 | 2016-03-09 | 大连理工大学 | 一种用于粗晶奥氏体不锈钢焊缝超声检测的图像降噪方法 |
CN107271370A (zh) * | 2017-07-04 | 2017-10-20 | 九江学院 | 一种基于材料内部缺陷检测的激光超声检测系统及其方法 |
US20180106765A1 (en) * | 2016-10-18 | 2018-04-19 | The Aerospace Corporation | Ultrasonic inspection technique to assess bond quality in composite structures |
-
2019
- 2019-07-23 CN CN201910669559.6A patent/CN110243945B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102207488A (zh) * | 2011-03-29 | 2011-10-05 | 北京理工大学 | 横波tofd缺陷定位方法 |
CN102393422A (zh) * | 2011-08-22 | 2012-03-28 | 江苏省产品质量监督检验研究院 | 基于超声tofd的缺陷离线判别方法 |
CN103149274A (zh) * | 2013-01-28 | 2013-06-12 | 中国科学院声学研究所 | 一种混凝土缺陷检测方法 |
CN104458912A (zh) * | 2014-12-17 | 2015-03-25 | 国家电网公司 | 瓷绝缘子用小角度纵波斜探头入射点入射角及延时测量法 |
CN105388218A (zh) * | 2015-12-10 | 2016-03-09 | 大连理工大学 | 一种用于粗晶奥氏体不锈钢焊缝超声检测的图像降噪方法 |
US20180106765A1 (en) * | 2016-10-18 | 2018-04-19 | The Aerospace Corporation | Ultrasonic inspection technique to assess bond quality in composite structures |
CN107271370A (zh) * | 2017-07-04 | 2017-10-20 | 九江学院 | 一种基于材料内部缺陷检测的激光超声检测系统及其方法 |
Non-Patent Citations (2)
Title |
---|
S. CONG ET AL: "Analysis on Ultrasonic TOFD Imaging Testing for Ultra-thick-walled EBW Joint of Aluminum Alloy", 《PROCEDIA ENGINEERING》 * |
丁宁 等: "基于波型转换的TOFD近表面盲区抑制研究", 《机械工程学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113358751A (zh) * | 2021-06-01 | 2021-09-07 | 中车青岛四方机车车辆股份有限公司 | 工件缺陷的检测方法、装置和系统 |
CN113358751B (zh) * | 2021-06-01 | 2022-09-06 | 中车青岛四方机车车辆股份有限公司 | 工件缺陷的检测方法、装置和系统 |
CN113484421A (zh) * | 2021-07-01 | 2021-10-08 | 中国工程物理研究院机械制造工艺研究所 | 基于合成孔径的激光超声内部缺陷多模成像方法及系统 |
CN113484421B (zh) * | 2021-07-01 | 2023-04-28 | 中国工程物理研究院机械制造工艺研究所 | 基于合成孔径的激光超声内部缺陷多模成像方法及系统 |
CN113758617A (zh) * | 2021-09-10 | 2021-12-07 | 哈尔滨工业大学 | 一种基于宽带扫频信号频域计算的应力梯度高效无损检测系统及其检测方法 |
CN114674928A (zh) * | 2022-03-21 | 2022-06-28 | 大连理工大学 | 基于二阶模式波的薄板结构超声tofd检测盲区抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110243945B (zh) | 2021-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110243945A (zh) | 基于合成孔径聚焦与模式转换波的超声tofd盲区抑制方法 | |
CN104034287B (zh) | 一种弹性各向异性金属基体热障涂层厚度超声测量方法 | |
WO2022151530A1 (zh) | 基于单角度楔块的临界折射纵波多材料检测系统及其声速测量方法 | |
CN111337171B (zh) | 一种应用于临界折射纵波应力检测的声时差测量方法 | |
CA2616900C (en) | Method for error-free checking of tubes for surface faults | |
Espinosa et al. | Accuracy on the time-of-flight estimation for ultrasonic waves applied to non-destructive evaluation of standing trees: a comparative experimental study | |
CN106198739A (zh) | 一种基于波型转换的tofd近表面盲区缺陷定位检测方法 | |
CN103543208B (zh) | 基于频谱分析原理减小tofd检测近表面盲区的方法 | |
CN104698089A (zh) | 一种适用于倾斜裂纹定量和成像的超声相对时间传播技术 | |
CN104897777A (zh) | 基于Burg算法的自回归谱外推技术提高TOFD检测纵向分辨率的方法 | |
CN111521136B (zh) | 一种基于水平剪切波的钢筋混凝土结构裂缝深度检测方法及检测装置 | |
CN111174894B (zh) | 一种激光超声横波声速测量方法 | |
Freeseman et al. | Characterization of concrete at various freeze-thaw damage conditions using SH-waves | |
CN114674928A (zh) | 基于二阶模式波的薄板结构超声tofd检测盲区抑制方法 | |
CN105044215A (zh) | 一种非破坏性的材料声速现场测量方法 | |
CN110440896B (zh) | 一种超声波测量系统及测量方法 | |
JP2011047763A (ja) | 超音波診断装置 | |
Sun et al. | Damage identification in thick steel beam based on guided ultrasonic waves | |
JP4405821B2 (ja) | 超音波信号検出方法及び装置 | |
JP3036387B2 (ja) | 超音波探傷方法及び装置 | |
US4596142A (en) | Ultrasonic resonance for detecting changes in elastic properties | |
Li et al. | Research on the imaging of concrete defect based on the pulse compression technique | |
CN114397365A (zh) | 一种超声波检测钢砼结构缺陷方法 | |
Bottiglieri et al. | Corrective Techniques for the Ultrasonic Nondestructive Evaluation of Ceramic Materials | |
EP0605720B1 (en) | Inspection apparatus |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |