起伏地表弹性波数值模拟的不等距差分法:以黄土塬区为例  

Non-equidistant difference method for elastic wave numerical simulation with irregular surface:a case study of loess tableland region

作　　者：汪登科 孙章庆[1] 韩复兴[1] 胡自多[2] 王瑞湖 陆济璞 岑文攀 高正辉[1] WANG Dengke;SUN Zhangqing;HAN Fuxing;HU Ziduo;WANG Ruihu;LU Jipu;CEN Wenpan;GAO Zhenghui(College for Geo-Exploration Science and Technology,Jilin University,Changchun 130026,China;Research Institute of Petroleum Exploration and Development,Northwest PetroChina,Lanzhou 730020,China;Geological Survey Institute of Guangxi Zhuang Autonomous Region,Nanning 530031,China;Bureau of Geology and Mineral Exploration and Development of Guangxi Zhuang Autonomous Region,Nanning 530023,China)

机构地区：[1]吉林大学地球探测科学与技术学院,长春130026 [2]中国石油勘探开发研究院西北分院,兰州730020 [3]广西壮族自治区地质调查院,南宁530031 [4]广西壮族自治区地质矿产勘查开发局,南宁530023

出　　处：《世界地质》2022年第4期842-852,872,共12页World Geology

基　　金：广西重点研发计划项目(2021AB30011);国家自然科学基金项目(42074150,41404085)联合资助。

摘　　要：针对包含复杂近地表及起伏地形的黄土塬区的正演模拟问题,提出一种弹性波不等距差分数值模拟方法。该方法不需要在地表处进行阶梯状离散,且在计算中可以将自由地表边界条件融入不等距差分格式中,能够提高计算效率。采用一个倾斜模型进行模拟,经与特征变量法对比,得到相同的波场,验证了该方法的有效性。对实际复杂黄土塬区模型进行试算,得到其波场快照与地震记录,证明了在复杂黄土塬区的适用性和可能的应用前景。In order to solve the seismic forward modeling of loess tableland region with complex near-surface and irregular surface,the non-equidistant difference method for elastic wave numerical simulation is proposed.This method does not need to perform staircase-shaped discretization at the surface,and the free surface boundary condition can be integrated into the non-equidistant difference scheme to improve the computational efficiency.An inclined model is used for simulation,and we get the same wavefield by comparing it with the characteristic variable method,which shows the effectiveness.The wavefield snapshot and seismograms are obtained by actual model of the complex loess tableland region,which proves the applicability and possible application prospects.

关 键 词：弹性波 不等距差分法 差分格式 边界条件 黄土塬 

分 类 号：P631.4[天文地球—地质矿产勘探]