在用有限差分法或有限元法模拟无界区域中的波动时,需要对计算区域的边界做特殊处理,以消除由于把地震波的传播设定在有限区域而产生的边界反射。为了这一目的,人们研究出了多种人工边界处理方法,完全匹配层(PML)吸收边界条件就是理想的方法之一,现已被广泛应用。本文将PML吸收边界条件应用于全波场地震波的数值模拟,数值计算实验表明,对qP波,匹配层的厚度为5个网格间距即可达到要求,而对qSV波与qSH波,为达到理想的吸收效果,匹配层的厚度应当增大,当厚度为13个网格间距时达到了理想的吸收效果。
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.