A new absorbing boundary method for diffusive-viscous wave equation (DVW) is proposed to deal with the artificial reflections from boundaries caused by a truncated computational domain in seismic wave numerical modeling. The main idea is to make the seismic waves exponentially attenuate with propagation distance by adding a proper absorbing layer around the boundary of computational domain. The solution of the wave equation in infinite homogenous space is obtained by using Fourier transform. Then, auxiliary equations are constructed by introducing decay functions such that the solution of auxiliary equations is the solution of the DVW equation in computational domain and decays exponentially at boundaries. The wave equation is finally solved by using finite-difference method in both homogeneous media and homogeneous layered media. Numerical results are given through comparison of the results from the proposed method and those of non-absorbing boundary condition, and show that boundary reflections are absorbed with the proposed method. The proposed method is also applicable to acoustic equation and Stokes equation.%针对地震波数值模拟中因计算区域有界而产生的人工边界反射问题,提出了一种弥散黏滞性波动方程的吸收边界算法.在所研究的有界计算区域周围加入合适的吸收层,使得吸收层内的地震波随传播距离指数衰减而达到吸收边界反射的目的.利用傅里叶变换求得无界空间该波动方程的解；引入衰减函数构建有界空间相应的辅助方程,合理地选择衰减函数使得计算区域内辅助方程的解为原方程的解,而边界吸收层内的解呈指数衰减；运用有限差分法在均匀介质和均匀层状介质中求解弥散黏滞性波动方程,采用该方法对边界进行处理并与未加边界处理的波场快照和地震记录进行对比.实验结果表明,所提方法处理后的波场快照几乎无边界反射波存在,在空间位置为x=147 m、z=747 m的地震记录中,0.5s处边界反射波的振幅趋于0.所提方法也适用于声波方程和Stokes方程.
展开▼
