基于SPML和海绵边界的伪谱法弹性波模拟复合吸收边界条件

2018年 57卷 第No. 1期
阅读:105
查看详情
A composite absorbing boundary based on the SPML and sponge absorbing boundary for pseudo-spectral elastic wave modeling
(1.中国海洋大学海洋地球科学学院海底科学与探测技术教育部重点实验室,山东青岛266100;2.青岛海洋科学与技术国家实验室海洋矿产资源评价与探测技术功能实验室,山东青岛266061;3.中国石油化工股份有限公司胜利油田分公司物探研究院,山东东营257022;4.中国石油集团东方地球物理勘探有限责任公司辽河物探分公司,辽宁盘锦124010)
(1.Key Laboratory of Seabed Science and Detection Technology MOE,College of Marine Geosciences,Ocean University of China,Qingdao 266100,China;2.Evaluation and Detection Technology Laboratory of Marine Mineral Resources,Qingdao National Laboratory for Marine Science and Technology,Qingdao 266061,China;3.Geophysical Research Institute of Shengli Oilfield,SINOPEC,Dongying 257022,China;4.Liaohe Geophysical Research Institute,Bureau of Geophysical Prospecting INC.,CNPC,Panjin 124010,China)

弹性波数值模拟中分裂格式的完全匹配层吸收边界(SPML)难以吸收与边界接近平行传播的地震波,并会产生数值噪声,降低了地震波模拟的精度。针对此问题,提出了复合吸收边界条件,将SPML边界条件与海绵吸收边界条件组合成边界的内外层,实现了对与边界接近平行传播地震波的有效吸收。针对数值噪声沿边界传播的特点,对复合边界中的海绵边界吸收系数进行了改进,提高了海绵边界对平行于边界方向传播的地震波的衰减能力,进一步提高了对入射至边界处地震波的吸收效果。数值测试结果显示,不论对简单速度模型还是复杂速度模型,复合吸收边界对地震波的吸收效果优于常规SPML边界的吸收效果,且不会产生明显的数值噪声。

In pseudo-spectral elastic wave modeling,the conventional split perfectly matched layer (SPML) is less effective in attenuating seismic waves that propagate almost parallel to the boundary and generates numerical noise,reducing the numerical simulation accuracy.Here,we propose a composite absorbing boundary (CAB) by adding a sponge-absorbing boundary layer outside the SPML absorbing boundary to attenuate the unabsorbed seismic waves through the SPML boundary.We modified the sponge boundary absorption coefficient in the CAB to avoid numerical noise generation.The numerical experiments show that the CAB has a higher attenuation capability than the conventional SPML for seismic waves traveling parallel to the boundary,and no significant noise is generated around the boundary.Since the proposed CAB adds an absorption layer outside the original PML and does not introduce additional complex mathematical transformations,it can be easily applied to other types of numerical simulations of seismic waves,such as finite difference and finite element methods.

完全匹配层; 海绵吸收边界; 复合吸收边界; 伪谱法; 弹性波正演模拟; 弹性波方程; 数值噪声;
perfectly matched layer,; sponge absorbing boundary,; composite absorbing boundary,; pseudo-spectral,; elastic wave modeling,; elastic wave equation,; numerical noise;

国家科技重大专项(2016ZX05006002-001)、国家自然科学基金重点项目(41230318)、山东省自然科学基金(ZR2014DM006)和教育部留学回国人员科研启动基金(教外司留[2015]1098号)共同资助。

10.3969/j.issn.1000-1441.2018.01.013