无单元法用于地震波波动方程模拟与成像

波动方程方法是解决地震正反演问题的基本工具之一.无单元法作为一种新兴的偏微分方程数值计算方法,已经在材料力学、热传导等领域取得了显著的成功.由于抛弃了单元的概念及采用滑动最小二乘的拟合方法,使得无单元法具有前处理简单、精度高、独立变量解高次连续等优点.本文首先介绍无单元法求解波动方程的原理,指出影响其精度的主要因素.在算例的基础上详细讨论了无单元法用于实际波动问题的效果,并进一步尝试利用无单元法进行地震波数值模拟和反演成像的研究.模型计算的结果表明,无单元法能够较好的处理地震模拟和成像问题,精度和稳定性是令人满意的.