[发明专利]获得稳定的近最佳匹配层吸收边界条件的方法

权利要求书

1.一种获得稳定的近最佳匹配层吸收边界条件的方法，其特征在于：包括以下步骤：

一、将时间域的弹性波一阶速度-应力方程变换到频率域；

二、对频率域弹性波一阶速度-应力方程进行坐标变换；

三、设置复频移-近最佳匹配层中的复频移拉伸函数参数；

四、将复频移动拉伸函数移入到被求导波场函数上；

五、对复频移拉伸函数进行变量代换重新表示为：并带入到拉伸后的频率域弹性波一阶速度-应力方程中；

六、引入辅助波场变量对拉伸后的频率域弹性波一阶速度-应力方程进行变量代换；

七、将拉伸后的频率域弹性波一阶速度-应力方程及辅助方程反变换到时间域，获得时间域最佳匹配层吸收边界条件的弹性波一阶速度-应力方程及其辅助常微分方程。

2.根据权利要求1所述的一种获得稳定的近最佳匹配层吸收边界条件的方法，其特征在于：步骤一中时间域的弹性波一阶速度-应力方程为二维的。

3.根据权利要求1所述的一种获得稳定的近最佳匹配层吸收边界条件的方法，其特征在于：步骤二中，坐标变换是利用复频移拉伸函数(θ∈{x，z})实现。

4.根据权利要求1所述的一种获得稳定的近最佳匹配层吸收边界条件的方法，其特征在于：步骤三是通过以下步骤完成：

(1)设定参数kθ的值；

(2)按照多轴最佳匹配层同时在多个方向上引入了衰减因子dθ＝dθ(x，z)，表示为：

(3)求得复频移拉伸函数参数αθ＝παmax(1-δθ/Lpml)。

5.根据权利要求4所述的一种获得稳定的近最佳匹配层吸收边界条件的方法，其特征在于：参数kθ为1。

6.根据权利要求1所述的一种获得稳定的近最佳匹配层吸收边界条件的方法，其特征在于：步骤五中，其中变量

7.根据权利要求1所述的一种获得稳定的近最佳匹配层吸收边界条件的方法，其特征在于：步骤四中，复频移拉伸函数为常数。