地震声波数值模拟中的吸收边界条件

在地震波传播的数值模拟中,在有限区域建立吸收边界条件是一个很重要的问题。在不同类型的吸收边界条件中,比较经典的是Clayton和Engquist提出的Clayton-Engquist吸收边界条件;PML吸收边界条件是Berenger在1994年提出的一种新的人工边界条件。这种边界条件几乎达到零反射。在地震波传播数值模拟中,PML已经被用于速度应力形式的一阶声波方程,而尚未被应用于位移形式的二阶声波方程。推导了二阶位移形式的PML公式,并应用于声波方程的数值模拟中;通过数值例子,与被传统应用的CE吸收边界条件作了比较。