基于复频移完全匹配层的辛龙格库塔算法探地雷达正演

为了消除探地雷达(ground penetrating radar, GPR)正演模拟的截断边界处的非物理反射,需要对模拟区域边界进行特殊处理。目前应用效果较好的边界处理方法为加载完全匹配层吸收边界条件。而复频移完全匹配层(complex frequency shift perfectly matched layer, CFS-PML)可以加强对低频区凋落波、倏失波的吸收效果,还可以压制信号的虚假反射,被引入到GPR正演边界条件处理中。将CFS-PML应用于辛龙格库塔求解二维GPR波动方程的数值模拟中,利用辛算法可以实现对时间坐标和空间坐标的保辛计算,从而实现GPR高精度正演模拟。首先,以横磁(transverse magnetic, TM)波为例,推导了基于CFS-PML的辛龙格库塔GPR波动方程求解公式。其次,利用狭长模型开展了CFS-PML吸收边界条件中的关键参数的选取试验,通过对比反射误差的大小确定了最优参数。此外,通过CFS-PML吸收边界条件与透射边界条件的对比实验,说明大角度掠射下CFS-PML边界条件较透射边界具有更好的吸收效果。最后,为了进一步验证该算法的效果,应用加...     

高阶复合完全匹配层的应用研究

研究了应用于时域算法的一种高阶完全匹配层(Perfectly Matched Layer,PML)方法,并提出了一种复合PML研究方法。在高阶PML原理的基础上,给出了复合PML的研究思路,及其应用于显式时域有限差分方法的过程。采用高阶PML和复合PML计算了波导的反射系数,结果表明有限元计算的结果是准确而稳定的。对比计算结果,可以看出与普通PML相比,高阶PML和复合PML对隐失波和各个频段的传输波都能起到良好的吸收效果,为有限元算法的广泛应用奠定了良好的基础。     

基于时域有限差分法的随钻声波测井仪隔声体隔声效果的数值模拟

随钻测井环境下,钻铤纵模式波会严重干扰地层纵波的测量,因此需要使用隔声体来减弱钻铤波对地层纵波测量的干扰.采用圆柱坐标系下的交错网格技术,自行开发了一种可用于设计随钻声波测井仪器隔声体的三维时域有限差分程序.利用该程序数值模拟采用不同结构的隔声体时,硬地层井眼中的随钻单极子声波测井响应,并对这几种不同结构的隔声体性能进行评价.结果表明:采用不同结构的隔声体,对钻铤纵模式波的时延量和隔声量存在很大差异;与其他结构隔声体相比,周期性交错排列的凹槽结构的隔声体在测量地层纵波速度时具有较好的效果.