波动方程方法是解决地震正反演问题的基本工具之一.无单元法作为一种新兴的偏微分方程数值计算方法,已经在材料力学、热传导等领域取得了显著的成功.由于抛弃了单元的概念及采用滑动最小二秉的拟合方法,使得无单元法具有前处理简单、精度高、独立变量解高次连续等优点.本文首先介绍无单元法求解波动方程的原理,指出影响其精度的主要因素.在算例的基础上详细讨论了无单元法用于实际波动问题的效果,并进一步尝试利用无单元法进行地震波数值模拟和反演成像的研究.模型计算的结果表明,无单元法能够较好的处理地震模拟和成像问题,精度和稳定性是令人满意的.
Wave equation method is one of the fundamental techniques for seismic modeling and imaging. In this paper we discuss a rather new numerical method-the Element-Free Method ( EFM ), which has demonstrated quite successful in elasticity, heat conduction and fatigue crack growth problems. The key point of this method is the absence of elements, which means potential cost savings. Besides, the Moving Least Squares (MLS) approximation in the EFM design leads to high accuracy and smooth derivatives. The basic principle and a simple example of vibrant film are discussed in details in this paper. Finally, some more examples will be given to show the good performance of EFM in seismic modeling and imaging problems.