于基三阶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.1选题背景和研究意义 1.2粘弹性介质国内外研究现状 1.3有限差分国内外研究现状 1.4本文主要研究内容 2 粘弹性介质的基本模型 3 3 方法介绍 3 3.1 Stereo-modeling方法简介 3 3.2 Lax-Wendroff correction方法简介 5 4 粘弹性介质中的波场数值模拟 6 4.1 波场快照 6 4.2 波形图 6 4.3 SEG模型的地表地震记录 10 5 结论 13 6 参考文献 15 1 绪论 1.1 选题背景和研究意义 随着我国石油和天然气工业的迅速发展，油气勘探的精度和难度也在不断地加大。这就要求我们必须深入了解地震波在地下介质中的传播规律，以便我们能更加准确地进行地下油气的勘探。地震波数值模拟技术是目前地震勘探领域人们研究地震传播规律的一个重要手段。传统地震波数值模拟技术，一般假设地层介质是均匀和完全弹性的，但是我们发现：在这种假设情况下，数值模拟的结果和我们在野外观察到的实际记录有很大差别。造成这种差别的主要原因是因为实际的地层介质并不是均匀和完全弹性的，我们应该用粘弹性介质来代替，这对于我们研究地震波波场的传播规律具有了非常重要的理论和实际意义。 1.2 粘弹性介质国内外研究现状 为了对实际地震资料的解释以及地震波的偏移反演等处理提供可靠依据，我们需要建立不同的介质模型来模拟地下介质，发展相应的理论和勘探方法。从以往的模拟结果来看，以经典弹性波理论为基础得到的理论记录实际地震波在大地中传播时所接收到的记录之间存在很大差别，为了消除这种差别，我们通常要对球面扩散、界面的透射损失等进行补偿，但是在补偿过后，理论记录上来自深层的地震波振幅仍然较小、低频成分较丰富，这说明地下介质对地震波有吸收衰减的作用，特别是对高频成分，吸收更为严重。地震波衰减受多种因素影响，事实上现在没有一种机理可以描述所有环境条件下产生的损耗，介质的粘滞性就是其中的一个主要原因，为了弄清楚地震波在地层中的衰减机制，很多研究人员开始对粘弹性介质的理论和应用方法进行研究。 国外研究现状 早在1845年Stokes就开始研究粘弹性介质中的地震波，他提出要考虑具有耗损的固体就必须使这种固体具有类似于粘滞液体的性质，因此他给固体的剪切模量附加上一个剪切粘滞，并建立了由于粘滞型内摩擦引起的能量耗损的Stokes粘弹性波动方程。但是粘弹介质理论的发展却是非常缓慢的，直到上个世纪40年代，美国地球物理学家N.K.Ricker首次在粘弹介质理论中完整详细的描述了地震波在地下介质传播中时的衰减问题。其后，人们纷纷开始研究粘弹介质中地震波的传播理论，提出了一系列的理论和研究方法。Aki和Richards(1980)建立了一种新的粘弹性介质模型——标准线性固体模型，相对于开尔芬—佛格特模型和麦克斯韦尔模型，标准线性固体模型更加适合描述地震波衰减的物理机制。粘弹介质地震波数值模拟技术的研究是从上世纪八十年代开始的，Day和Minster (1984)第一次成功的在粘弹性介质中利用Pade近似法进行2-D时间域数值模拟。EmmeriCh和Korn (1987)发现这种方法存在明显的缺陷即质量低劣和计算效率低，于是他们提出了一种新的近似称之为“广义标准线性固体”(GSLS)，同时对粘弹性模量的有理式进行推导，并发展二维有限差分算法使之适合标量波的传播，提高