基于辛几何格式的叠前逆时偏移成像的工程结构无损检测研究论文.docx

基于辛几何格式的叠前逆时偏移成像的工程结构无损检测研究论文偏移技术可以根据高频近似构造射线理论的方法成像，也可以从波动方程的角度成像。基于射线理论的方法计算量小，能够正确显示内部缺陷的尺寸，不能够正确显示内部结构的方位；基于波动方程的方法成像精度（包括位置精度）很高，但计算量大。偏移技术在地球物理勘探方面有广泛而成熟的应用，近几年在无损检测方面也取得了一些可喜的成果。根据对波动方程的处理方法的不同，波动方程偏移可归纳为三类：空间时间域直接离散法、频率波数域的二维傅里叶变换法（F-K）和空间时间域格林函数克希霍夫积分法（Kirchhoff）。有限差分法属于第一类算法。从波动方程求解方法的稳定性来看，此类算法不是根据能量泛函推导而来，人为规定离散点与离散点之间的差分关系，在长期跟踪计算能力上表现不佳。本文将研究空间时间域的另一种算法，即以辛几何数值模拟算为基础，结合逆时偏移思想，构建一种新的叠前偏移成像算法,并给出典型的数值模拟实例和混凝土结构成像实验结果。辛几何格式的逆时偏移辛几何格式的叠前逆时偏移理论的数学基础是辛几何格式的数值模拟方法和逆时偏移成像条件，其主要技术手段是对在结构表面采集到超声波信号进行处理，使反射超声波归位并呈现出结构内部反射面的影像。数值模拟方法应力波数值模拟是指在已知介质结构模型的情况下，研究应力波在地下各种介质中传播规律的一种数值模拟方法，其理论基础就是表征应力波在介质中传播的应力波理论。而辛几何数值模拟方法是对于某连续体结构，在空间域采用传统的某种能量泛函方法如有限元、无网格伽辽金方法、边界元法等进行空间离散，把无限自由度问题转化为有限自由度问题，从而得到时间域的连续常微分方程，再利用辛几何精细积分算法进行逐步积分，获得任意时刻所有空间变量的位移场。波动方程偏移成像技术的逆时偏移成像条件对于不同的发收波场，有不同的成像条件。零发收距（自发自收）观测数据的成像条件GAMcmechan等在研究零炮检距观测数据逆时外推偏移算法时，根据爆炸反射界面原理，提出了“零时间”成像条件。共发射点观测数据的成像条件FChangW等针对共炮点数据，将“零时间”成像条件加以推广，提出了“激发时间”成像条件。它一般应用于非零发收距观测数据的叠前逆时偏移计算中。这种共发射点观测数据的成像有两种实现方法：零延迟互相关成像和激发时成像。零延迟互相关成像的基本过程是，将各时刻的正时与逆时声场对应相乘，然后将各时刻的相乘所获得声场的值进行累加，即得到成像剖面。激发时成像的基本过程是，从逆时波场中依照震源到成像点的单程旅行时提取相应的波场，即得到该点的成像值，各点成像的值的总和即是偏移剖面。建立有限元模型按照20mm的尺度划分网格，保存数据,并求出有限元离散方程的系数矩阵。根据激励条件，设置好载荷向量，长度为n.划分网格的细密程度根据计算机的内存和实际计算需要综合确定。在本实验中，由于采用普通PC计算，所以采用了计算机能够计算的最细密网格划分方式。建立一阶哈密顿方程用勒让德变换将拉格朗日体系下的n维方程转化为哈密尔顿体系下的`2n维方程，同时方程的阶次由2降为1,为方程在时域的展开求解提供了便利条件。获取正向波场时间步长为0.1s,在某点设置冲击波载荷，频率为200MHz,通过哈密顿方程求解出正向波场。其中，任意时刻任一点的声波信号幅值在第6个步骤中有应用，表面测线上节点的信号在第5个步骤中应用。获取逆向波场获取逆向波场的具体步骤是：在某点进行激励得到的波场在表面测线上被记录下来后，沿着时间轴逆转，把测线上记录下来的所有信号当作作用在表面上节点上的激励，再通过计算获取新的声波场。根据成像条件获取缺陷图像将每一个激励信号得到的正向波场和衍生出的逆向波场点乘，然后把所有信号得到的点乘结果叠加，便得到成像。窄裂缝的位置能够准确显示，方向显图2含窄裂缝结构的辛偏移成像结果Fig.2Symplecticmigrationimagingresultsofstructurewiththenarrowcrack示不够准确，这是由于缺陷的大小与超声激励主波长具有相同量级，部分超声波会发生绕射现象，反射波减弱，因此只能显示缺陷中部一部分的反射面，而缺陷边缘产生绕射的部分，未能有效显示。另外，底面的位置也能准确显示，但经过底面但不能通过一次反射到达表面的部分，不能够使底面的这一部分正确成像。因此缺陷正下方、正下方外围一部分区域，以及底面的外边缘不能显示，或显示较弱。几何模型在实际工程中，混凝土结构的质量参差不一，又常受到复杂的荷载，因此所发生的缺陷的形状也千姿百态。将四个具有低声阻抗性质的薄膜裹紧木板预埋入混凝土，模拟混凝土构件中的两个水平缺陷和两个倾斜缺陷。实验和数据采集利用数字信号发生器、功率放大器和数字示波器组成的测试系统和1GHz高频超声换能器进行实