全波场地震波数值模拟中的边界处理实践研究.pdfVIP

学兔兔 全波场地震波数值模拟中的边界处理实践研究 刘保童 全波场地震波数值模拟中的边界处理实践研究* 刘保童 (天水师范学院物理与信息科学学院 甘肃天水，741001) 摘 要：在用有限差分法或有限元法模拟无界区域中的波动时，需要对计算区域的边界做特殊处理，以消除由于把地 震波的传播设定在有限区域而产生的边界反射。为了这一目的，人们研究出了多种人工边界处理方法，完全匹配层(PML) 吸收边界条件就是理想的方法之一，现已被广泛应用。本文将PML吸收边界条件应用于全波场地震波的数值模拟，数值计 算实验表明，对 qP 波，匹配层的厚度为 5 个网格间距即可达到要求，而对 qSV 波与 qSH 波，为达到理想的吸收效果，匹配 层的厚度应当增大，当厚度为13个网格间距时达到了理想的吸收效果。 关键词：数值模拟；吸收边界；完全匹配层(PML)；各向异性 Abstract：The simulation of waves by finite-difference or finite-element methods in unbounded domains requires a special treat- ment for the boundaries of the necessarily truncated computational domain. For this purpose, many processing approaches of bound- ary are worked out, one of ideal methods is PML, which has been used widely. This paper uses PML absorbing boundary condition in the simulation of full seismic wavefield. Numerical experiments showed that the thickness of PML for qSV and qSH waves is larg- er than P wave, 13 grids is appropriate. Key words: Numerical simulation; Absorbing boundary; PML; Anisotropy 中图分类号：P631.4 文献标识码：A 文章编号：1001-9227 （2014） 01-0032-03 0 引 言 地震波场数值模拟是研究地震波传播规律、识别各种波的 重要方法，对采集设计、资料处理起着非常重要的作用[1]。在 用有限差分法或有限元法模拟无界区域中的波动??，需要对计 算区域的边界做特殊处理[2]。为了这一目的，人们研究出了多 种人工边界处理方法，代表性的有：Clayton 和 Engquist 的吸收 边界条件[3]，Reynolds 的透明边界条件[4]，Cerjan 等人提出的衰 减 边 界 [5]， Berenger 提 出 的 完 全 匹 配 层 (Perfec tly Matc hed 图1 PML吸收边界示意图 Layer，PML)吸收边界[6]，等等。PML吸收边界是Bere