CN107356677A - 基于旅行时层析成像和逆时偏移成像的超声无损检测方法 - Google Patents
基于旅行时层析成像和逆时偏移成像的超声无损检测方法 Download PDFInfo
- Publication number
- CN107356677A CN107356677A CN201710565844.4A CN201710565844A CN107356677A CN 107356677 A CN107356677 A CN 107356677A CN 201710565844 A CN201710565844 A CN 201710565844A CN 107356677 A CN107356677 A CN 107356677A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- msup
- mtr
- 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/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0654—Imaging
- G01N29/0672—Imaging by acoustic tomography
-
- 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
- G01N2291/0232—Glass, ceramics, concrete or stone
-
- 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
- G01N2291/0234—Metals, e.g. steel
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开的基于旅行时层析成像和逆时偏移成像的超声无损检测方法,包括采用旅行时层析成像方法获取待测结构断面内的超声波速度分布,将其作为逆时偏移成像的初始速度模型,再采用逆时偏移成像方法对待测结构内部缺陷进行成像。本发明可对混凝土内部空洞和钢管混凝土剥离等缺陷进行高分辨率成像,从而获取它们的位置、尺寸和形状等信息,实际应用前景广阔。
Description
技术领域
本发明涉及超声无损检测技术领域,尤其涉及一种基于旅行时层析成像和逆时偏移成像的超声无损检测方法。
背景技术
在无损检测领域当中,超声无损检测技术因其检测对象范围广泛、对人体无害且成本低廉等优点而发展迅猛。它是一种在不破坏和损伤物体的前提下,利用超声波对物体的性能、质量以及物体内是否有缺陷进行检测的技术,现如今在很多的工业部门被广泛使用,如钢铁工业、机器制造工业有关的工业部门,高速发展中的新技术产业如集成电路工业等工业部门。目前大量应用于材料和构件的检测如混凝土检测、木材检测等。但是传统超声无损检测技术诸如基于时频分析技术的超声无损检测技术只能定性地分析是否存在缺陷以及缺陷可能存在的大概区域,并不能定量地去分析缺陷形状及大小并把缺陷准确成像出来。
旅行时层析成像(Travel time Tomography),是一种通过提取的初至时来反演计算区域速度分布的成像方法,被认为是重建近地表模型的有效方法而广泛使用,一般应用在井间地震勘探中,对混凝土以及木材等的缺陷检测却鲜有涉及。同时,基于射线追踪理论的旅行时层析成像对慢度沿炮点到检波器之间的线性积分是一种数学抽象,不能反映真实的物理运动过程,因此会发生低速体偏离、高速区聚集的情况,从而对层析成像的分辨率产生极大的影响。
逆时偏移成像(Reverse Time Migration),简称RTM,是一种基于全波动方程,主要应用在非常规油气勘探和井间地震成像,是目前偏移方法中最准确的成像方法。但是,逆时偏移成像方法对初始模型的精度要求比较高,只有初始模型精度够高,成像效果才能足够显著。
发明内容
本发明的目的在于提供一种基于旅行时层析成像和逆时偏移成像的超声无损检测方法,。
为达成上述技术目的,本发明的技术方案如下:
基于旅行时层析成像和逆时偏移成像的超声无损检测方法,包括采用旅行时层析成像方法获取待测结构断面内的超声波速度分布,将其作为逆时偏移成像的初始速度模型,再采用逆时偏移成像方法对待测结构内部缺陷进行成像。
进一步,所述采用旅行时层析成像方法得到待测结构断面的超声波速度分布作为逆时偏移成像的初始速度模型,具体包括以下步骤:
步骤1、在待测结构四周设置多组超声波激发源和检波器,并在激发源位置处激发超声波信号,在检波器位置处采集超声波传输信号;
步骤2、提取每一道超声波传输信号的初至波到达时间tmea;
步骤3、对结构横断面区域进行离散网格剖分,并按一定的顺序进行编号,设定每个网格内的初始超声波速,得到一个初始慢度模型s0,利用走时和慢度的关系、检波器和激发源的位置以及初至走时tmea进行初次反演,走时计算值和慢度的关系可以表示为:
其中,tn cal表示第n条射线的走时计算值,lnm表示第n条射线在元胞m中的射线长度,sm表示元胞m中的慢度值;
步骤4、对每一个激发源的位置,基于每个网格节点的走时使用梯度下降法进行弯曲射线追踪,得到相对应的系数矩阵L,利用快速推进算法计算已知速度场各网格节点的波前传播时间,二维介质中,传播时间和慢度满足下列程函方程:
其中,T和S分别表示点(x,y)处的波前传播时间和慢度,用有限差分的方式近似上式可以表示为:
其中,Δx和Δy分别为x和y方向上的离散网格间距,并且有
T1=min(Ti-1,j,Ti+1,j)
T2=min(Ti,j-1,Ti,j+1)
其中,Ti-1,j表示第i-1,j个离散网格点处的波前传播时间;
步骤5、依据步骤4,得到初始慢度模型下的走时矩阵t0,计算新的走时矩阵t0和走时t之间的误差Erms,当误差Erms大于允许误差时,继续迭代计算,由走时t0可以得到新的慢度模型s1,以此类推,不断计算走时矩阵ti和走时t之间的误差Erms,迭代反演,不断得到新的慢度矩阵;采用的迭代公式可以用如下的矩阵方程表示为:
其中,Cd为数据协方差矩阵,Cm为模型协方差矩阵,α为正则化因子,Lk为理论计算待测的第k次迭代后的系数矩阵,tmea为测量提取的走时矩阵,sk表示第k次迭代后得到的慢度矩阵,理论计算的时间和测量的走时之间的误差用均方根误差表示:
其中,n表示总的收发道数(收发组合),表示第k+1次迭代后的第i道射线的走时计算值;
步骤6、当误差Erms在允许的范围内或者迭代次数达到要求值时,则停止迭代,得到需要的慢度矩阵si,根据慢度矩阵si得到速度模型,把得到的速度模型以及根据此速度模型得到的密度模型一起作为逆时偏移成像的输入模型;
步骤7、由速度模型获得获得此初始模型的密度分布,一起作为逆时偏移成像的输入模型,在指定激发源位置处放置源信号,求解速度-应力一阶方程使源波场进行正向传播,到达最大时间步tmax时停止正向传播,保存每一个时间步的源波场。
步骤8、进行检波器波场逆时延拓,对单一激励源所对应的所有检波器所采集的超声波信号进行逆时处理后,在这些检波器所在位置激励逆时后的检波器信号,求解速度-应力一阶方程使检波器场反向延拓;
步骤9、应用归一化互相关成像条件对步骤7中的激发源波场和步骤8中的检波器波场应用归一化互相关成像条件,其表达式如下:
其中,Sn(t,z,x)表示由第n个激发器激励的正向延拓源波场,Rn(t,z,x)表示由第n个激发器对应的所有检波器信号逆时后同时激励得到的检波器反向延拓波场,而I(x,z)表示重构的图像,x、z分别表示横纵轴;
步骤10、根据逆时偏移成像结果确定结构内部缺陷的位置、尺寸和形状。
步骤7的速度-应力一阶方程的具体计算公式如下:
其中:ρ为介质的密度,vx、vy、vz分别是x、y、z方向的速度,P是压力,fx、fy、fz分别为x、y、z方向的点力源密度,gp是压力源密度,λ是常数。
采用上述方案后,本发明的积极效果在于:主要对实际工业应用中的缺陷进行无损检测并精确成像。基于旅行时层析成像获得超声波在计算区域中的最佳路径,得到速度分布,并依据此速度分布作为逆时偏移成像的初始模型进行成像,从而定量分析区域内部的缺陷大小和位置。本发明创新性地结合了传统的旅行时层析成像和逆时偏移成像技术,能做到比传统的旅行时层析成像方法精确度更高,比传统的逆时偏移成像方法实用性更强,是一种非常实用的超声无损检测方法,可对混凝土内部空洞和钢管混凝土剥离等缺陷进行高分辨率成像,从而获取它们的位置、尺寸和形状等信息,实际应用前景广阔。
附图说明
图1为本发明超声无损检测方法的流程图;
图2射线在离散元胞中的走时示意图;
图3旅行时层析成像流程图;
图4逆时偏移成像流程图。
具体实施方式
如图1所示,本发明基于旅行时层析成像和逆时偏移成像的超声无损检测方法,具体包括以下步骤:
步骤1、在待测结构四周设置多组超声波激发源和检波器,并在激发源位置处激发超声波信号,在检波器位置处采集超声波传输信号;
步骤2、提取每一道超声波传输信号的初至波到达时间tmea;
步骤3、如图2所示,对结构横断面区域进行离散网格剖分,并按一定的顺序进行编号,设定每个网格内的初始超声波速,得到一个初始慢度模型s0,利用走时和慢度的关系、检波器和激发源的位置以及初至走时tmea进行初次反演,走时计算值和慢度的关系可以表示为:
其中,tn cal表示第n条射线的走时计算值,lnm表示第n条射线在元胞m中的射线长度,sm表示元胞m中的慢度值;
步骤4、对每一个激发源的位置,基于每个网格节点的走时使用梯度下降法进行弯曲射线追踪,得到相对应的系数矩阵L,利用快速推进算法计算已知速度场各网格节点的波前传播时间,二维介质中,传播时间和慢度满足下列程函方程:
其中,T和S分别表示点(x,y)处的波前传播时间和慢度,用有限差分的方式近似上式可以表示为:
其中,Δx和Δy分别为x和y方向上的离散网格间距,并且有
T1=min(Ti-1,j,Ti+1,j)
T2=min(Ti,j-1,Ti,j+1)
其中,Ti-1,j表示第i-1,j个离散网格点处的波前传播时间;
步骤5、如图3所示,依据步骤4,得到初始慢度模型下的走时矩阵t0,计算新的走时矩阵t0和走时t之间的误差Erms,当误差Erms大于允许误差时,继续迭代计算,由走时t0可以得到新的慢度模型s1,以此类推,不断计算走时矩阵ti和走时t之间的误差Erms,迭代反演,不断得到新的慢度矩阵;采用的迭代公式可以用如下的矩阵方程表示为:
其中,Cd为数据协方差矩阵,Cm为模型协方差矩阵,α为正则化因子,Lk为理论计算待测的第k次迭代后的系数矩阵,tmea为测量提取的走时矩阵,sk表示第k次迭代后得到的慢度矩阵,理论计算的时间和测量的走时之间的误差用均方根误差表示:
其中,n表示总的收发道数(收发组合),表示第k+1次迭代后的第i道射线的走时计算值;
步骤6、当误差Erms在允许的范围内或者迭代次数达到要求值时,则停止迭代,得到需要的慢度矩阵si,根据慢度矩阵si得到速度模型,把得到的速度模型以及根据此速度模型得到的密度模型一起作为逆时偏移成像的输入模型;
上述采用逆时偏移成像方法并基于初始速度模型对缺陷位置、形状进行成像,结合图4,具体包括以下步骤:
步骤7、由速度模型获得获得此初始模型的密度分布,一起作为逆时偏移成像的输入模型,在指定激发源位置处放置源信号,求解速度-应力一阶方程使源波场进行正向传播,到达最大时间步tmax时停止正向传播,保存每一个时间步的源波场。采用有限差分时域(FDTD)法求解速度应力的一阶微分方程,这是求解偏微分方程的最常用的数值方法。这种方法在时域上实现了超声波场的正向延拓;其中步骤7的速度-应力一阶方程的具体计算公式如下:
其中:ρ为介质的密度,vx、vy、vz分别是x、y、z方向的速度,P是压力,fx、fy、fz分别为x、y、z方向的点力源密度,gp是压力源密度,λ是常数;
步骤8、进行检波器波场逆时延拓,对单一激励源所对应的所有检波器所采集的超声波信号进行逆时处理后,在这些检波器所在位置激励逆时后的检波器信号,求解速度-应力一阶方程使检波器场反向延拓;
步骤9、应用归一化互相关成像条件对步骤7中的激发源波场和步骤8中的检波器波场应用归一化互相关成像条件,其表达式如下:
其中,Sn(t,z,x)表示由第n个激发器激励的正向延拓源波场,Rn(t,z,x)表示由第n个激发器对应的所有检波器信号逆时后同时激励得到的检波器反向延拓波场,而I(x,z)表示重构的图像,x、z分别表示横纵轴;
步骤10、根据逆时偏移成像结果确定结构内部缺陷的位置、尺寸和形状。
Claims (3)
1.基于旅行时层析成像和逆时偏移成像的超声无损检测方法,其特征在于,包括采用旅行时层析成像方法获取待测结构断面内的超声波速度分布,将其作为逆时偏移成像的初始速度模型,再采用逆时偏移成像方法对待测结构内部缺陷进行成像。
2.如权利要求1所述的基于旅行时层析成像和逆时偏移成像的超声无损检测方法,其特征在于,所述采用旅行时层析成像方法得到待测结构断面的超声波速度分布作为逆时偏移成像的初始速度模型,具体包括以下步骤:
步骤1、在待测结构四周设置多组超声波激发源和检波器,并在激发源位置处激发超声波信号,在检波器位置处采集超声波传输信号;
步骤2、提取每一道超声波传输信号的初至波到达时间tmea;
步骤3、对结构横断面区域进行离散网格剖分,并按一定的顺序进行编号,设定每个网格内的初始超声波速,得到一个初始慢度模型s0,利用走时和慢度的关系、检波器和激发源的位置以及初至走时tmea进行初次反演,走时计算值和慢度的关系可以表示为:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>l</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>l</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mn>1</mn>
<mi>m</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>l</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>l</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mn>2</mn>
<mi>m</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mi>n</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mi>n</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<msub>
<mi>l</mi>
<mrow>
<mi>n</mi>
<mi>m</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>s</mi>
<mi>m</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msup>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mrow>
<mi>c</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
<mrow>
<mi>c</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>...</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<msub>
<mi>t</mi>
<mi>n</mi>
</msub>
<mrow>
<mi>c</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,tn cal表示第n条射线的走时计算值,lnm表示第n条射线在元胞m中的射线长度,sm表示元胞m中的慢度值;
步骤4、对每一个激发源的位置,基于每个网格节点的走时使用梯度下降法进行弯曲射线追踪,得到相对应的系数矩阵L,利用快速推进算法计算已知速度场各网格节点的波前传播时间,二维介质中,传播时间和慢度满足下列程函方程:
<mrow>
<mo>|</mo>
<mo>&dtri;</mo>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
<mo>=</mo>
<msup>
<mi>S</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,T和S分别表示点(x,y)处的波前传播时间和慢度,用有限差分的方式近似▽T,上式可以表示为:
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>x</mi>
</mrow>
</mfrac>
<mo>,</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>y</mi>
</mrow>
</mfrac>
<mo>,</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>=</mo>
<msub>
<msup>
<mi>S</mi>
<mn>2</mn>
</msup>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</mrow>
其中,Δx和Δy分别为x和y方向上的离散网格间距,并且有
T1=min(Ti-1,j,Ti+1,j)
T2=min(Ti,j-1,Ti,j+1)
其中,Ti-1,j表示第i-1,j个离散网格点处的波前传播时间;
步骤5、依据步骤4,得到初始慢度模型下的走时矩阵t0,计算新的走时矩阵t0和走时t之间的误差Erms,当误差Erms大于允许误差时,继续迭代计算,由走时t0可以得到新的慢度模型s1,以此类推,不断计算走时矩阵ti和走时t之间的误差Erms,迭代反演,不断得到新的慢度矩阵;采用的迭代公式可以用如下的矩阵方程表示为:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>C</mi>
<mi>d</mi>
</msub>
<msup>
<mi>L</mi>
<mi>k</mi>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&alpha;C</mi>
<mi>m</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mo>(</mo>
<msup>
<mi>s</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mi>s</mi>
<mi>k</mi>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>C</mi>
<mi>d</mi>
</msub>
<mo>(</mo>
<msup>
<mi>t</mi>
<mrow>
<mi>m</mi>
<mi>e</mi>
<mi>a</mi>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mi>L</mi>
<mi>k</mi>
</msup>
<msup>
<mi>s</mi>
<mi>k</mi>
</msup>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
<mo>...</mo>
<mo>)</mo>
</mrow>
</mrow>
其中,Cd为数据协方差矩阵,Cm为模型协方差矩阵,α为正则化因子,Lk为理论计算待测的第k次迭代后的系数矩阵,tmea为测量提取的走时矩阵,sk表示第k次迭代后得到的慢度矩阵,理论计算的时间和测量的走时之间的误差用均方根误差表示:
<mrow>
<msup>
<msub>
<mi>E</mi>
<mrow>
<mi>r</mi>
<mi>m</mi>
<mi>s</mi>
</mrow>
</msub>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>=</mo>
<msqrt>
<mrow>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mo>|</mo>
<msup>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mrow>
<mi>c</mi>
<mi>a</mi>
<mi>l</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>t</mi>
<mi>i</mi>
</msub>
<mrow>
<mi>m</mi>
<mi>e</mi>
<mi>a</mi>
</mrow>
</msup>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
<mo>...</mo>
<mo>)</mo>
</mrow>
</mrow>
其中,n表示总的收发道数(收发组合),ti cal,k+1表示第k+1次迭代后的第i道射线的走时计算值;
步骤6、当误差Erms在允许的范围内或者迭代次数达到要求值时,则停止迭代,得到需要的慢度矩阵si,根据慢度矩阵si得到速度模型,把得到的速度模型以及根据此速度模型得到的密度模型一起作为逆时偏移成像的输入模型;
步骤7、由步骤6的速度模型获得获得此初始模型的密度分布,一起作为逆时偏移成像的输入模型,在指定激发源位置处放置源信号,求解速度-应力一阶方程使源波场进行正向传播,到达最大时间步tmax时停止正向传播,保存每一个时间步的源波场;
步骤8、进行检波器波场逆时延拓,对单一激励源所对应的所有检波器所采集的超声波信号进行逆时处理后,在这些检波器所在位置激励逆时后的检波器信号,求解速度-应力一阶方程使检波器场反向延拓;
步骤9、应用归一化互相关成像条件对步骤7中的激发源波场和步骤8中的检波器波场应用归一化互相关成像条件,其表达式如下:
<mrow>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mi>n</mi>
</munder>
<mfrac>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mi>t</mi>
</munder>
<msub>
<mi>S</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>z</mi>
<mo>,</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>z</mi>
<mo>,</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mi>t</mi>
</munder>
<msup>
<msub>
<mi>R</mi>
<mi>n</mi>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>,</mo>
<mi>z</mi>
<mo>,</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中,Sn(t,z,x)表示由第n个激发器激励的正向延拓源波场,Rn(t,z,x)表示由第n个激发器对应的所有检波器信号逆时后同时激励得到的检波器反向延拓波场,而I(x,z)表示重构的图像,x、z分别表示横纵轴;
步骤10、根据逆时偏移成像结果确定结构内部缺陷的位置、尺寸和形状。
3.如权利要求2所述的基于旅行时层析成像和逆时偏移成像的超声无损检测方法,其特征在于,步骤7的速度-应力一阶方程的具体计算公式如下:
<mrow>
<mi>&rho;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>v</mi>
<mi>x</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>P</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>x</mi>
</msub>
</mrow>
<mrow>
<mi>&rho;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>v</mi>
<mi>y</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>P</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>y</mi>
</msub>
</mrow>
<mrow>
<mi>&rho;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>v</mi>
<mi>z</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>P</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>z</mi>
</msub>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>P</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>v</mi>
<mi>x</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>v</mi>
<mi>y</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>v</mi>
<mi>z</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>g</mi>
<mi>p</mi>
</msub>
</mrow>
其中:ρ为介质的密度,vx、vy、vz分别是x、y、z方向的速度,P是压力,fx、fy、fz分别为x、y、z方向的点力源密度,gp是压力源密度,λ是常数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710565844.4A CN107356677B (zh) | 2017-07-12 | 2017-07-12 | 基于旅行时层析成像和逆时偏移成像的超声无损检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710565844.4A CN107356677B (zh) | 2017-07-12 | 2017-07-12 | 基于旅行时层析成像和逆时偏移成像的超声无损检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107356677A true CN107356677A (zh) | 2017-11-17 |
CN107356677B CN107356677B (zh) | 2020-02-07 |
Family
ID=60293353
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710565844.4A Active CN107356677B (zh) | 2017-07-12 | 2017-07-12 | 基于旅行时层析成像和逆时偏移成像的超声无损检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107356677B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108320284A (zh) * | 2018-02-06 | 2018-07-24 | 重庆交通大学 | 对钢管混凝土结构界面状态的检测方法及装置 |
CN110470742A (zh) * | 2019-04-25 | 2019-11-19 | 西安交通大学 | 一种管道弯头缺陷的精确检测方法 |
CN110907629A (zh) * | 2019-10-21 | 2020-03-24 | 武汉大学 | 标定钢管混凝土柱约束力与超声波的试验方法 |
CN112014018A (zh) * | 2020-09-01 | 2020-12-01 | 西南交通大学 | 一种基于超声层析成像的应力场测量方法 |
CN113702504A (zh) * | 2020-07-28 | 2021-11-26 | 广州建筑产业研究院有限公司 | 一种套筒灌浆饱满度检测方法及计算机设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1764828A (zh) * | 2003-01-30 | 2006-04-26 | 应用计算机科学促进协会 | 摄录物体声学图像的方法和装置以及相应的计算机程序和相应的计算机可读存储介质 |
CN101171194A (zh) * | 2005-05-02 | 2008-04-30 | 法比奥·泼尼股份公司 | 无芯网材卷筒及其生产机械和方法 |
CN104374828A (zh) * | 2014-11-18 | 2015-02-25 | 上海岩土工程勘察设计研究院有限公司 | 一种隐患探测的超声波层析成像方法 |
US20160376004A1 (en) * | 2015-03-16 | 2016-12-29 | XCraft Enterprises, LLC | Unmanned aerial vehicle with detachable computing device |
-
2017
- 2017-07-12 CN CN201710565844.4A patent/CN107356677B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1764828A (zh) * | 2003-01-30 | 2006-04-26 | 应用计算机科学促进协会 | 摄录物体声学图像的方法和装置以及相应的计算机程序和相应的计算机可读存储介质 |
CN101171194A (zh) * | 2005-05-02 | 2008-04-30 | 法比奥·泼尼股份公司 | 无芯网材卷筒及其生产机械和方法 |
CN104374828A (zh) * | 2014-11-18 | 2015-02-25 | 上海岩土工程勘察设计研究院有限公司 | 一种隐患探测的超声波层析成像方法 |
US20160376004A1 (en) * | 2015-03-16 | 2016-12-29 | XCraft Enterprises, LLC | Unmanned aerial vehicle with detachable computing device |
Non-Patent Citations (1)
Title |
---|
陈可洋: "两种不同偏移起始面逆时成像处理及效果分析――以起伏地表理论模型为例", 《油气藏评价与开发》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108320284A (zh) * | 2018-02-06 | 2018-07-24 | 重庆交通大学 | 对钢管混凝土结构界面状态的检测方法及装置 |
CN110470742A (zh) * | 2019-04-25 | 2019-11-19 | 西安交通大学 | 一种管道弯头缺陷的精确检测方法 |
CN110470742B (zh) * | 2019-04-25 | 2020-06-19 | 西安交通大学 | 一种管道弯头缺陷的精确检测方法 |
CN110907629A (zh) * | 2019-10-21 | 2020-03-24 | 武汉大学 | 标定钢管混凝土柱约束力与超声波的试验方法 |
CN110907629B (zh) * | 2019-10-21 | 2020-10-13 | 武汉大学 | 标定钢管混凝土柱约束力与超声波的试验方法 |
CN113702504A (zh) * | 2020-07-28 | 2021-11-26 | 广州建筑产业研究院有限公司 | 一种套筒灌浆饱满度检测方法及计算机设备 |
CN113702504B (zh) * | 2020-07-28 | 2023-06-20 | 广州建筑产业研究院有限公司 | 一种套筒灌浆饱满度检测方法及计算机设备 |
CN112014018A (zh) * | 2020-09-01 | 2020-12-01 | 西南交通大学 | 一种基于超声层析成像的应力场测量方法 |
CN112014018B (zh) * | 2020-09-01 | 2021-08-27 | 西南交通大学 | 一种基于超声层析成像的应力场测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107356677B (zh) | 2020-02-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107356677A (zh) | 基于旅行时层析成像和逆时偏移成像的超声无损检测方法 | |
Liu et al. | Reverse time migration of acoustic waves for imaging based defects detection for concrete and CFST structures | |
Huthwaite | Evaluation of inversion approaches for guided wave thickness mapping | |
CN102636568B (zh) | 一种检测混凝土内部缺陷的有限元超声成像方法 | |
Chai et al. | Single-side access tomography for evaluating interior defect of concrete | |
CN110346453B (zh) | 一种混凝土结构内缺陷小排列反射回波快速检测方法 | |
Chekroun et al. | Analysis of coherent surface wave dispersion and attenuation for non-destructive testing of concrete | |
US11906472B2 (en) | Non-destructive concrete stress evaluation | |
Luchin et al. | Sonic tomography for masonry walls characterization | |
Jeon et al. | 2D-wavelet wavenumber filtering for structural damage detection using full steady-state wavefield laser scanning | |
CN107741488A (zh) | 隧道施工期岩体完整性系数Kv的计算方法 | |
CN104931584B (zh) | 基于压缩采样理论的超声ct检测方法 | |
CN104181234B (zh) | 一种基于多重信号处理技术的无损检测方法 | |
Wu et al. | Two-step detection of concrete internal condition using array ultrasound and deep learning | |
CN100464185C (zh) | 混凝土超声层析成像算法 | |
Lin et al. | Estimation of wave velocity for ultrasonic imaging of concrete structures based on dispersion analysis | |
Zhu | Non-contact NDT of concrete structures using air coupled sensors | |
Chen et al. | Evaluation of ultrasonic SH-waveform tomography for determining cover thickness and rebar size in concrete structures | |
Wang et al. | High resolution tomography of pipeline using multi-helical Lamb wave based on compressed sensing | |
Jaganathan | Multichannel surface wave analysis of reinforced concrete pipe segments using longitudinal and circumferential waves induced by a point impact | |
Capizzi et al. | Some tests of 3D ultrasonic traveltime tomography on the Eleonora d'Aragona statue (F. Laurana, 1468) | |
CN105403627A (zh) | 一种用于增强超声检测图像横向分辨率的方法 | |
CN100495022C (zh) | 混凝土超声层析成像算法 | |
Cheng et al. | Defects visualization using acoustic emission tomography technique | |
Suryono et al. | Ultrasonic Computed Tomography System for Concrete Inspection |
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 |