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

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