基于三阶Adams格式求解声波方程多步算法.docVIP

创新项目论文 一种基于三阶Adams格式的求解声波方程的多步算法 摘 要 一个准确、高效、低数值频散的正演算法能够提高反演精度、加快反演收敛速度，因此研究地震波场正演模拟技术具有重要意义。区别于传统的空间离散方法，利用空间插值, 用网格点处的函数值及其梯度共同逼近空间高阶偏导数的方法称为近似解析离散化方法。声波方程通过变换，并采用近似解析离散化方法进行空间离散，从而转变成为一个半离散化的常微分方程组再利用三阶显式Adams格式进行时间推进求解半离散化的常微分方程组，从而得到了一个新的求解声波方程的有限差分方法（AD-STEM）。对AD-STEM进行了理论误差和数值误差分析、计算效率比较和数值波场模拟。研究表明，与传统法AD-LWC比较， AD-STEM方法数值精度更高，数值频散更低，更高效，且与解析解匹配更好。AD-STEM法能够通过压制数值频散而提高计算效率在无可见数值频散的条件下，AD-STEM的计算速度是AD-LWC的1.88倍，而存储量只有其72%，更适合在粗网格下进行大规模地震波场数值模拟。 关键词：近似解析离散化；三阶Adams；数值频散；有限差分  目 录 1 绪论 1 1.2粘弹性介质国内外研究现状 1.3有限差分国内外研究现状 1.4本文主要研究内容 2 粘弹性介质的基本模型 3 3 方法介绍 3 3.1 Stereo