CN107356677B - 基于旅行时层析成像和逆时偏移成像的超声无损检测方法 - Google Patents

Abstract

本发明公开的基于旅行时层析成像和逆时偏移成像的超声无损检测方法，包括采用旅行时层析成像方法获取待测结构断面内的超声波速度分布，将其作为逆时偏移成像的初始速度模型，再采用逆时偏移成像方法对待测结构内部缺陷进行成像。本发明可对混凝土内部空洞和钢管混凝土剥离等缺陷进行高分辨率成像，从而获取它们的位置、尺寸和形状等信息，实际应用前景广阔。

Description

基于旅行时层析成像和逆时偏移成像的超声无损检测方法

技术领域

本发明涉及超声无损检测技术领域，尤其涉及一种基于旅行时层析成像和逆时偏移成像的超声无损检测方法。

背景技术

在无损检测领域当中，超声无损检测技术因其检测对象范围广泛、对人体无害且成本低廉等优点而发展迅猛。它是一种在不破坏和损伤物体的前提下，利用超声波对物体的性能、质量以及物体内是否有缺陷进行检测的技术，现如今在很多的工业部门被广泛使用，如钢铁工业、机器制造工业有关的工业部门，高速发展中的新技术产业如集成电路工业等工业部门。目前大量应用于材料和构件的检测如混凝土检测、木材检测等。但是传统超声无损检测技术诸如基于时频分析技术的超声无损检测技术只能定性地分析是否存在缺陷以及缺陷可能存在的大概区域，并不能定量地去分析缺陷形状及大小并把缺陷准确成像出来。

旅行时层析成像(Travel time Tomography)，是一种通过提取的初至时来反演计算区域速度分布的成像方法，被认为是重建近地表模型的有效方法而广泛使用，一般应用在井间地震勘探中，对混凝土以及木材等的缺陷检测却鲜有涉及。同时，基于射线追踪理论的旅行时层析成像对慢度沿炮点到检波器之间的线性积分是一种数学抽象，不能反映真实的物理运动过程，因此会发生低速体偏离、高速区聚集的情况，从而对层析成像的分辨率产生极大的影响。

逆时偏移成像(Reverse Time Migration)，简称RTM，是一种基于全波动方程，主要应用在非常规油气勘探和井间地震成像，是目前偏移方法中最准确的成像方法。但是，逆时偏移成像方法对初始模型的精度要求比较高，只有初始模型精度够高，成像效果才能足够显著。

发明内容

本发明的目的在于提供一种基于旅行时层析成像和逆时偏移成像的超声无损检测方法，。

为达成上述技术目的，本发明的技术方案如下：

基于旅行时层析成像和逆时偏移成像的超声无损检测方法，包括采用旅行时层析成像方法获取待测结构断面内的超声波速度分布，将其作为逆时偏移成像的初始速度模型，再采用逆时偏移成像方法对待测结构内部缺陷进行成像。

进一步，所述采用旅行时层析成像方法得到待测结构断面的超声波速度分布作为逆时偏移成像的初始速度模型，具体包括以下步骤：

步骤1、在待测结构四周设置多组超声波激发源和检波器，并在激发源位置处激发超声波信号，在检波器位置处采集超声波传输信号；

步骤2、提取每一道超声波传输信号的初至波到达时间t_mea；

步骤3、对结构横断面区域进行离散网格剖分，并按一定的顺序进行编号，设定每个网格内的初始超声波速，得到一个初始慢度模型s⁰，利用走时和慢度的关系、检波器和激发源的位置以及初至走时t_mea进行初次反演，走时计算值和慢度的关系可以表示为：

其中，t_n ^cal表示第n条射线的走时计算值，l_nm表示第n条射线在元胞m中的射线长度，s_m表示元胞m中的慢度值；

步骤4、对每一个激发源的位置，基于每个网格节点的走时使用梯度下降法进行弯曲射线追踪，得到相对应的系数矩阵L，利用快速推进算法计算已知速度场各网格节点的波前传播时间，二维介质中，传播时间和慢度满足下列程函方程：

其中，T和S分别表示点(x,y)处的波前传播时间和慢度，用有限差分的方式近似

上式可以表示为：

其中，Δx和Δy分别为x和y方向上的离散网格间距，并且有

T₁＝min(T_i-1,j,T_i+1,j)

T₂＝min(T_i,j-1,T_i,j+1)

其中，T_i-1,j表示第i-1,j个离散网格点处的波前传播时间；

步骤5、依据步骤4，得到初始慢度模型下的走时矩阵t⁰，计算新的走时矩阵t⁰和走时t之间的误差E_rms，当误差E_rms大于允许误差时，继续迭代计算，由走时t⁰可以得到新的慢度模型s¹，以此类推，不断计算走时矩阵tⁱ和走时t之间的误差E_rms，迭代反演，不断得到新的慢度矩阵；采用的迭代公式可以用如下的矩阵方程表示为：

其中，C_d为数据协方差矩阵，C_m为模型协方差矩阵，α为正则化因子，L^k为理论计算待测的第k次迭代后的系数矩阵，t^mea为测量提取的走时矩阵，s^k表示第k次迭代后得到的慢度矩阵，理论计算的时间和测量的走时之间的误差用均方根误差表示：

其中，n表示总的收发道数(收发组合)，

表示第k+1次迭代后的第i道射线的走时计算值；

步骤6、当误差E_rms在允许的范围内或者迭代次数达到要求值时，则停止迭代，得到需要的慢度矩阵sⁱ，根据慢度矩阵sⁱ得到速度模型，把得到的速度模型以及根据此速度模型得到的密度模型一起作为逆时偏移成像的输入模型；

步骤7、由速度模型获得获得此初始模型的密度分布，一起作为逆时偏移成像的输入模型，在指定激发源位置处放置源信号，求解速度-应力一阶方程使源波场进行正向传播，到达最大时间步t_max时停止正向传播，保存每一个时间步的源波场。

步骤8、进行检波器波场逆时延拓，对单一激励源所对应的所有检波器所采集的超声波信号进行逆时处理后，在这些检波器所在位置激励逆时后的检波器信号，求解速度-应力一阶方程使检波器场反向延拓；

步骤9、应用归一化互相关成像条件对步骤7中的激发源波场和步骤8中的检波器波场应用归一化互相关成像条件，其表达式如下：

其中，S_n(t,z,x)表示由第n个激发器激励的正向延拓源波场，R_n(t,z,x)表示由第n个激发器对应的所有检波器信号逆时后同时激励得到的检波器反向延拓波场，而I(x,z)表示重构的图像，x、z分别表示横纵轴；

步骤10、根据逆时偏移成像结果确定结构内部缺陷的位置、尺寸和形状。

步骤7的速度-应力一阶方程的具体计算公式如下：

其中：ρ为介质的密度，v_x、v_y、v_z分别是x、y、z方向的速度，P是压力，f_x、f_y、f_z分别为x、y、z方向的点力源密度，g_p是压力源密度，λ是常数。

采用上述方案后，本发明的积极效果在于：主要对实际工业应用中的缺陷进行无损检测并精确成像。基于旅行时层析成像获得超声波在计算区域中的最佳路径，得到速度分布，并依据此速度分布作为逆时偏移成像的初始模型进行成像，从而定量分析区域内部的缺陷大小和位置。本发明创新性地结合了传统的旅行时层析成像和逆时偏移成像技术，能做到比传统的旅行时层析成像方法精确度更高，比传统的逆时偏移成像方法实用性更强，是一种非常实用的超声无损检测方法，可对混凝土内部空洞和钢管混凝土剥离等缺陷进行高分辨率成像，从而获取它们的位置、尺寸和形状等信息，实际应用前景广阔。

附图说明

图1为本发明超声无损检测方法的流程图；

图2射线在离散元胞中的走时示意图；

图3旅行时层析成像流程图；

图4逆时偏移成像流程图。

具体实施方式

如图1所示，本发明基于旅行时层析成像和逆时偏移成像的超声无损检测方法，具体包括以下步骤：

步骤1、在待测结构四周设置多组超声波激发源和检波器，并在激发源位置处激发超声波信号，在检波器位置处采集超声波传输信号；

步骤2、提取每一道超声波传输信号的初至波到达时间t_mea；

步骤3、如图2所示，对结构横断面区域进行离散网格剖分，并按一定的顺序进行编号，设定每个网格内的初始超声波速，得到一个初始慢度模型s⁰，利用走时和慢度的关系、检波器和激发源的位置以及初至走时t_mea进行初次反演，走时计算值和慢度的关系可以表示为：

其中，t_n ^cal表示第n条射线的走时计算值，l_nm表示第n条射线在元胞m中的射线长度，s_m表示元胞m中的慢度值；

步骤4、对每一个激发源的位置，基于每个网格节点的走时使用梯度下降法进行弯曲射线追踪，得到相对应的系数矩阵L，利用快速推进算法计算已知速度场各网格节点的波前传播时间，二维介质中，传播时间和慢度满足下列程函方程：

其中，T和S分别表示点(x,y)处的波前传播时间和慢度，用有限差分的方式近似

上式可以表示为：

其中，Δx和Δy分别为x和y方向上的离散网格间距，并且有

T₁＝min(T_i-1,j,T_i+1,j)

T₂＝min(T_i,j-1,T_i,j+1)

其中，T_i-1,j表示第i-1,j个离散网格点处的波前传播时间；

步骤5、如图3所示，依据步骤4，得到初始慢度模型下的走时矩阵t⁰，计算新的走时矩阵t⁰和走时t之间的误差E_rms，当误差E_rms大于允许误差时，继续迭代计算，由走时t⁰可以得到新的慢度模型s¹，以此类推，不断计算走时矩阵tⁱ和走时t之间的误差E_rms，迭代反演，不断得到新的慢度矩阵；采用的迭代公式可以用如下的矩阵方程表示为：

其中，C_d为数据协方差矩阵，C_m为模型协方差矩阵，α为正则化因子，L^k为理论计算待测的第k次迭代后的系数矩阵，t^mea为测量提取的走时矩阵，s^k表示第k次迭代后得到的慢度矩阵，理论计算的时间和测量的走时之间的误差用均方根误差表示：

其中，n表示总的收发道数(收发组合)，表示第k+1次迭代后的第i道射线的走时计算值；

步骤6、当误差E_rms在允许的范围内或者迭代次数达到要求值时，则停止迭代，得到需要的慢度矩阵sⁱ，根据慢度矩阵sⁱ得到速度模型，把得到的速度模型以及根据此速度模型得到的密度模型一起作为逆时偏移成像的输入模型；

上述采用逆时偏移成像方法并基于初始速度模型对缺陷位置、形状进行成像，结合图4，具体包括以下步骤：

步骤7、由速度模型获得获得此初始模型的密度分布，一起作为逆时偏移成像的输入模型，在指定激发源位置处放置源信号，求解速度-应力一阶方程使源波场进行正向传播，到达最大时间步t_max时停止正向传播，保存每一个时间步的源波场。采用有限差分时域(FDTD)法求解速度应力的一阶微分方程，这是求解偏微分方程的最常用的数值方法。这种方法在时域上实现了超声波场的正向延拓；其中步骤7的速度-应力一阶方程的具体计算公式如下：

其中：ρ为介质的密度，v_x、v_y、v_z分别是x、y、z方向的速度，P是压力，f_x、f_y、f_z分别为x、y、z方向的点力源密度，g_p是压力源密度，λ是常数；

步骤8、进行检波器波场逆时延拓，对单一激励源所对应的所有检波器所采集的超声波信号进行逆时处理后，在这些检波器所在位置激励逆时后的检波器信号，求解速度-应力一阶方程使检波器场反向延拓；

步骤9、应用归一化互相关成像条件对步骤7中的激发源波场和步骤8中的检波器波场应用归一化互相关成像条件，其表达式如下：

其中，S_n(t,z,x)表示由第n个激发器激励的正向延拓源波场，R_n(t,z,x)表示由第n个激发器对应的所有检波器信号逆时后同时激励得到的检波器反向延拓波场，而I(x,z)表示重构的图像，x、z分别表示横纵轴；

步骤10、根据逆时偏移成像结果确定结构内部缺陷的位置、尺寸和形状。

Claims (2)

1.基于旅行时层析成像和逆时偏移成像的超声无损检测方法，其特征在于，包括采用旅行时层析成像方法获取待测结构断面内的超声波速度分布，将其作为逆时偏移成像的初始速度模型，再采用逆时偏移成像方法对待测结构内部缺陷进行成像；

所述采用旅行时层析成像方法得到待测结构断面的超声波速度分布作为逆时偏移成像的初始速度模型，具体包括以下步骤：

步骤1、在待测结构四周设置多组超声波激发源和检波器，并在激发源位置处激发超声波信号，在检波器位置处采集超声波传输信号；

步骤2、提取每一道超声波传输信号的初至波到达时间t_mea；

步骤3、对结构横断面区域进行离散网格剖分，并按一定的顺序进行编号，设定每个网格内的初始超声波速，得到一个初始慢度模型s⁰，利用走时和慢度的关系、检波器和激发源的位置以及初至走时t_mea进行初次反演，走时计算值和慢度的关系表示为：

其中，t_n ^cal表示第n条射线的走时计算值，l_nm表示第n条射线在元胞m中的射线长度，s_m表示元胞m中的慢度值；

步骤4、对每一个激发源的位置，基于每个网格节点的走时使用梯度下降法进行弯曲射线追踪，得到相对应的系数矩阵L，利用快速推进算法计算已知速度场各网格节点的波前传播时间，二维介质中，传播时间和慢度满足下列程函方程：

其中，T和S分别表示点(x,y)处的波前传播时间和慢度，用有限差分的方式近似

上式表示为：

其中，Δx和Δy分别为x和y方向上的离散网格间距，并且有

T₁＝min(T_i-1,j,T_i+1,j)

T₂＝min(T_i,j-1,T_i,j+1)

其中，T_i-1,j表示第i-1,j个离散网格点处的波前传播时间；

步骤5、依据步骤4，得到初始慢度模型下的走时矩阵t⁰，计算新的走时矩阵t⁰和走时t之间的误差E_rms，当误差E_rms大于允许误差时，继续迭代计算，由走时t⁰得到新的慢度模型s¹，以此类推，不断计算走时矩阵tⁱ和走时t之间的误差E_rms，迭代反演，不断得到新的慢度矩阵；采用的迭代公式用如下的矩阵方程表示为：

其中，C_d为数据协方差矩阵，C_m为模型协方差矩阵，α为正则化因子，L^k为理论计算待测的第k次迭代后的系数矩阵，t^mea为测量提取的走时矩阵，s^k表示第k次迭代后得到的慢度矩阵，理论计算的时间和测量的走时之间的误差用均方根误差表示：

其中，n表示总的收发道数，t_i ^cal,k+1表示第k+1次迭代后的第i道射线的走时计算值；

步骤6、当误差E_rms在允许的范围内或者迭代次数达到要求值时，则停止迭代，得到需要的慢度矩阵sⁱ，根据慢度矩阵sⁱ得到速度模型，把得到的速度模型以及根据此速度模型得到的密度模型一起作为逆时偏移成像的输入模型；

步骤7、由步骤6的速度模型获得此初始模型的密度分布，一起作为逆时偏移成像的输入模型，在指定激发源位置处放置源信号，求解速度-应力一阶方程使源波场进行正向传播，到达最大时间步t_max时停止正向传播，保存每一个时间步的源波场；

步骤8、进行检波器波场逆时延拓，对单一激励源所对应的所有检波器所采集的超声波信号进行逆时处理后，在这些检波器所在位置激励逆时后的检波器信号，求解速度-应力一阶方程使检波器场反向延拓；

步骤9、应用归一化互相关成像条件对步骤7中的激发源波场和步骤8中的检波器波场应用归一化互相关成像条件，其表达式如下：

其中，S_n(t,z,x)表示由第n个激发器激励的正向延拓源波场，R_n(t,z,x)表示由第n个激发器对应的所有检波器信号逆时后同时激励得到的检波器反向延拓波场，而I(x,z)表示重构的图像，x、z分别表示横纵轴；

步骤10、根据逆时偏移成像结果确定结构内部缺陷的位置、尺寸和形状。

2.如权利要求1所述的基于旅行时层析成像和逆时偏移成像的超声无损检测方法，其特征在于，步骤7的速度-应力一阶方程的具体计算公式如下：

其中：ρ为介质的密度，v_x、v_y、v_z分别是x、y、z方向的速度，P是压力，f_x、f_y、f_z分别为x、y、z方向的点力源密度，g_p是压力源密度，λ是常数。

Images