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

作者简介:

通讯作者:

中图分类号:

P631

基金项目:

国家自然科学基金项目(面上项目,重点项目,重大项目)


Complex frequency shifted perfectly matched layer for Symplectic Partitioned Runge-Kutta method forward modeling of Ground Penetrating Radar
Author:
Affiliation:

Fund Project:

The National Natural Science Foundation of China (General Program, Key Program, Major Research Plan)

  • 摘要
  • |
  • 图/表
  • |
  • 访问统计
  • |
  • 参考文献
  • |
  • 相似文献
  • |
  • 引证文献
  • |
  • 资源附件
  • |
  • 文章评论
    摘要:

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

    Abstract:

    To eliminate the non-physical reflection at the truncated boundary of GPR forward modeling, it is necessary to deal with the simulated region boundary. At present, the well applied method is to load the perfectly matched layer, and the complex frequency shift perfectly matched layer (CFS-PML), which can enhance the absorption effect of litter wave and prosperous lost wave in low frequency region and suppress the false reflection of signals, is widely used in GPR forward modeling. In this paper, CFS-PML is applied to the numerical simulation of solving the two-dimensional GPR wave equation with the symplectic partitioned Runge-Kutta (SPRK) algorithm, of which the symplectic algorithm can be utilized to realize the symplectic calculation of the time and space coordinates, so as to accomplish the high-precision forward modeling of GPR. Firstly, taking the TM wave as an example, we deduce the SPRK of GPR wave equation solution formula based on CFS-PML. Then, a narrow and long model is utilized to conduct the experiment of the key parameters’ selections of CFS-PML absorbing boundary condition, and the optimal parameters were determined by comparing the reflection error. Moreover, the comparison between CFS-PML and transmission boundary shows that the absorption effect of CFS-PML boundary condition at large Angle grazing is better than that of transmission boundary. Finally, to further verify the effectiveness of our proposed algorithm, the SPRK program loaded with the CFS-PML boundary condition is conduct on a geoelectric model with undulating topography. Then the forward GPR cross-sectional view, which is attained by using both cross-sectional method and wide-angle method, proved that the algorithm can fit the uneven underground interface well.

    参考文献
    相似文献
    引证文献
引用本文

冯德山,谭阳,杨军,等. 基于复频移完全匹配层的辛龙格库塔算法探地雷达正演[J]. 科学技术与工程, 2022, 22(9): 3473-3484.
Feng Deshan, Tan Yang, Yang Jun, et al. Complex frequency shifted perfectly matched layer for Symplectic Partitioned Runge-Kutta method forward modeling of Ground Penetrating Radar[J]. Science Technology and Engineering,2022,22(9):3473-3484.

复制
分享
文章指标
  • 点击次数:
  • 下载次数:
  • HTML阅读次数:
  • 引用次数:
历史
  • 收稿日期:2021-07-26
  • 最后修改日期:2022-03-10
  • 录用日期:2021-11-10
  • 在线发布日期: 2022-03-25
  • 出版日期:
×
律回春渐,新元肇启|《科学技术与工程》编辑部恭祝新岁!
亟待确认版面费归属稿件,敬请作者关注