基于正交贴体网格的VTI介质地震波数值模拟  

Numerical Simulation of Seismic Wave in VTI Media Based on Orthogonal Body-fitted Grid

作　　者：刘志强 李钢柱 黄磊[1] 许磊 李文宝[1] LIU Zhi-qiang;LI Gang-zhu;HUANG Lei;XU Lei;LI Wen-bao(Water Conservancy and Civil Engineering College,Inner Mongolia Agricultural University,Hohhot 010018,Inner Mongolia,China)

机构地区：[1]内蒙古农业大学水利与土木建筑工程学院,内蒙古呼和浩特010018

出　　处：《地球科学与环境学报》2023年第2期168-178,共11页Journal of Earth Sciences and Environment

基　　金：内蒙古自治区高等学校科学技术研究项目(NJZZ19042);内蒙古自治区自然科学基金项目(2021BS04007,2021MS04024);内蒙古自治区科技计划项目(2020GG0076)。

摘　　要：为了解决规则网格有限差分法在处理起伏地表自由边界条件时需进行复杂坐标转换和插值运算的问题,以及对起伏地表进行阶梯状离散近似所产生的虚假散射波问题,将计算流体力学中的正交贴体网格生成技术引入到起伏地表下VTI介质的网格剖分中,采用基于同位网格的DRP/opt MacCormack有限差分法模拟起伏地表下VTI介质中的地震波场。网格的正交性使实施自由边界条件时无需做复杂的坐标转换和插值运算。数值算例表明,本文提出的网格剖分法能够有效消除阶梯状网格产生的虚假散射波,从而提高起伏地表下VTI介质地震波数值模拟精度。In order to solve the complex coordinate transformation and interpolation operations when the regular grid finite difference method deals with the free boundary conditions of irregular surface,as well as the false scattered waves generated by the stepped discrete approximation of irregular surface,the orthogonal body-fitted gird generation technique in computational fluid mechanics was introduced into the mesh generation of VTI media under irregular surface.The DRP/opt MacCormack finite difference method based on the collocated grid was used to simulate seismic waves in VTI media under the irregular surface.The orthogonality of the grid eliminates the need for complex coordinate transformation and interpolation operations in the implementation of free boundary conditions.Numerical examples show that the proposed grid generation method can effectively eliminate the false scattered waves generated by the stepped grid,thus improve the accuracy of seismic numerical simulation in VTI media under the irregular surface.

关 键 词：正交贴体网格 VTI介质 地震波 数值模拟 起伏地表 自由边界条件 有限差分法 阶梯状网格 

分 类 号：P315.31[天文地球—地震学]