各向同性地层中声波测井的有限差分模拟-毕业(学术)论文设计.docVIP

MACROBUTTON MTEditEquationSection2 SEQ MTEqn \r \h \* MERGEFORMAT SEQ MTSec \r 1 \h \* MERGEFORMAT SEQ MTChap \r 2 \h \* MERGEFORMAT 第1章?绪 论 1.1 课题背景和意义 弹性波是应力波的一种，扰动或外力作用引起的应力和应变在弹性介质中传递的形式。弹性介质中质点间存在着相互作用的弹性力。某一质点因受到扰动或外力的作用而离开平衡位置后，弹性恢复力使该质点发生振动，从而引起周围质点的位移和振动，于是振动就在弹性介质中传播，并伴随有能量的传递。在振动所到之处应力和应变就会发生变化。弹性波理论被广泛应用于地球物理和资源勘探、材料的无损探伤、工程结构的抗震抗爆、岩土动力学等领域。比如，利用不同介质中弹性波性质的变化，可以反演地层参数，获得地层信息；利用弹性波遇界面后的反射和折射等现象，可以探测材料中是否有损伤情况。为了研究各种工程问题所涉及的弹性波问题，需要对这些问题中的弹性波传播进行数值模拟。由于弹性波在复杂介质结构中传播时的解析解，往往很难获得，这就需要利用有限差分、有限元、边界元等数值算法来实现。 本文研究既有流体又有固体的复杂结构中弹性波传播的有限差分模拟问题。有限差分(Finite Difference, FD)是出现最早、发展最成熟的一种基于微分方程的数值算法，因其具有方法简单、概念清晰等特点被广泛应用于各种波场问题的数值分析。该方法的基本原理是[1]：将求解域划分为差分网格，用有限个网格节点代替连续的求解域，进而采用离散的代数形式的有限差分方程近似代替连续性的微分方程，在代数方程中将空间各点待求量的值与其临近点的值联系起来。 与有限元相比，在对弹性波传播问题的模拟方面，有限差分法更容易实现，而且计算效率明显更高。 本文以声波测井为背景，用时域有限差分法模拟无限大固体空间中井孔内弹性波的传播。由于计算空间不能无限大，需要添加吸收边界，本文使用目前公认吸收效果最好的PML吸收边界条件来处理截断边界引起的波反射，并采用不分裂场量的PML设置方式。PML层是包围有限计算空间的、各向异性的薄层，它可以毫无反射地吸收入射波，并具有很高的耗散性。 需要说明的是，尽管本文针对工程上的声波测井进行有限差分的数值模拟，所实现的算法和程序可以很方便地修改从而用于其它有关固体中弹性波传播的工程问题。 1.2 国内外研究现况 传统的有限差分法主要用于求解标量问题，多用于求解电位场、温度场等分布问题，当然也能用于模拟声波、弹性波的传播问题[2-5]。为了更方便、直接地求解矢量波场问题，1966年Yee[6]设计了新的空间网格，同时求解电磁场所有分量，这种方法被称为时域有限差分(Finite Difference Time Domain, FDTD)方法。该方法对电磁场的电场和磁场分量在空间和时间上采取交替抽样的离散方式，每一个电场（或磁场）分量周围有四个磁场（或电场）分量环绕，应用这种离散方式将含时间变量的电磁场麦克斯韦旋度方程转化为差分格式，并在时间轴上逐步推进地求解空间电磁场。Yee提出的这种抽样方式后来被称为Yee元胞。此后，这种抽样方式被移植到声波、弹性波模拟问题中，发展出了交错网格的速度-应力有限差分法(Velocity-Stress Finite Difference, VSFD)[ 7-11]。 在FDTD中，由于计算机存储量的限制，包含有限差分的网格必须是有限的，比如FDTD网格可以位于一个盒子内。为了使所求的解唯一，则必须在这一盒子的表面加上一定的边界条件。但在很多应用中，介质是无限延续的，在此情况下如果在盒子边界上加以边界条件就会引起反射，而在实际的物理过程中是不应该有这些反射的。因此，为了模拟无限大的介质，应在盒子的表面加以吸收边界条件(Absorbing boundary condition, ABC)以消除反射。早期的ABC包括：单向波ABC，比如Clayton和Engquist[12]，Reynolds[13]等；海绵层ABC，比如Cerjan等[14]。尽管这些ABC成功应用在许多数值模拟问题中，但它们只在特定的入射角和频率范围内有效，且吸收效率不高。 1994年，Berenger[15] 提出了在计算区域以外设置完全匹配层(Perfectly matched layer, PML)的概念,开辟了ABC研究的新方向，并获得巨大的成功。与上述ABC完全不同， PML利用非物理媒质构成计算区域的边界，并具有如下特性：平面波可以无反射地射入其中，并在内部按指数规律衰减。PML的吸收效率比上述ABC高出几个数量级。 上述PML的应用中，对PML区域中的场进行了裂化处理，以避免进行时