【摘要】 在对弹性波方程进行数值模拟时 ,低阶差分格式往往产生严重的数值频散 ,高阶显示差分格式需要用较多的网格点 ,不利于边界的处理。而紧致差分格式吸收了它们的优点 ,弥补了它们的不足。为此该文应用紧致差分格式的思想 ,发展了二维情况下弹性波方程初值问题的紧致差分方法 ,研究了它的稳定性 ,并用 Fourier方法分析了显示差分格式和紧致差分格式的相速度误差 ,最后利用紧致差分方法在粗网格条件下对地震波传播进行了数值模拟 ,并同五点四阶中心差分方法的计算结果进行了对比。结果表明 ,求解弹性波方程的紧致差分方法有效 ,且具有比同网格点差分格式更高的计算精度和较小的数值频散。更多还原

【Abstract】 Simulation of wave motion in elastic media using finite difference methods generate too many phase velocity errors which make the results unusable or require a wide stencil which make simulation of the boundary condition very difficult. The paper presents a compact difference scheme (CDS) for two-dimensional elastic wave equations and analysis its stability and phase velocity error. Comparision of seismic propagation simulations using the CDS with a large grid site, with results of an explicit centered finite difference method shows that the CDS for the elastic wave equations has higher accuracy and less numerical dispersion.更多还原