非均匀介质中交错网格高阶有限差分数值模拟

地震波场的数值模拟一直是地球物理学的一个重要的研究领域,而在数值正演模拟方法的研究中,计算精度和计算效率是评价该方法有效性及优越性的二个关键问题。这里从一阶速度—应力弹性波动方程出发,着重介绍如何构造离散化模型的网格,如何求解空间导数,如何选取边界条件等内容,从而更有效地提高数值计算的精度与计算效率。文中构造了不同类型的介质模型,并在交错网格中,利用高阶有限差分模拟非均匀介质的波场传播。模拟结果表明,该方法实现简单,具有很好地稳定性和较高的精度,能够直观、高效地反映出介质中波场的传播规律。     

天然地震数值模拟与实例对比

为了解决长时间正演模拟的不稳定性问题, 实现天然地震波场长时间数值模拟.在高阶有限差分数值模拟的基础上, 给出了多轴向完全匹配吸收边界(M-PML)二维划分方式, 讨论了M-PML吸收边界转换系数P的取值对其吸收能力与截断误差对数值模拟有效信号的影响.实现了二维弹性波高阶有限差分的长时间数值模拟.采用2013年12月16日湖北省巴东县地震及余震信息, 对该地震数据进行了数值模拟, 并与武汉、秭归两台站数据进行了分析和讨论, 验证该正演模拟方法的正确性和有效性.为天然地震波场传播规律、成像和震相识别及震源定位等研究提供了基础.      

大地电磁三维正演并行算法研究

大地电磁法三维正演算法计算量大,用传统的串行程序计算相当耗时。而三维正演是逐个频率按顺序计算的,并行性好,适合并行运算。结合MPI自身的优越性,在深入分析大地电磁三维正演串行程序实现流程的基础上,确定了并行计算的思路,实现了三维正演的并行计算。通过三个理论模型对实现的三维正演并行程序进行了试算,分析对比了在多种情况下程序的执行效率。测试结果表明,所实现的三维正演并行程序运行结果正确,效率提高明显。此思路可为解决其它地球物理超大计算量问题所借鉴。     

黏弹TTI介质中旋转交错网格高阶有限差分数值模拟

以Carcione黏弹各向异性理论为基础,给出了适用于黏弹性具有任意倾斜对称轴横向各向同性介质(黏弹TTI介质)的二维三分量一阶速度-应力方程,采用旋转交错网格任意偶数阶精度有限差分格式求解该方程,并推导出了二维黏弹TTI介质完全匹配层(PML)吸收边界条件公式和相应的旋转交错网格任意偶数阶精度有限差分格式,实现了该类介质的地震波场数值模拟.数值模拟结果表明:该方法模拟精度高,边界吸收效果好,可以得到高精度的波场快照和合成记录;并且波场快照和合成记录能较好地反映地下介质的各向异性特征和黏弹性特征.     

含煤地层各向异性介质有限差分数值模拟

煤储层是一种双重孔隙介质,其孔隙由基质孔隙和裂隙组成,并具有自身独特的割理系统。传统的煤田地震勘探将煤层假设为单一的各向同性介质,没有考虑其中裂隙对煤层的影响。利用交错网格高阶有限差分算法,分别对不同厚度VTI和HTI型煤层的弹性波传播特征进行数值模拟,获得波场快照和含煤层状介质的地震模拟记录。从波场快照中可以看出,弹性波在VTI和HTI型煤层中传播存在明显差异;煤层厚度变化对波场的敏感性要强于煤层类型的变化。      

地震波数值模拟中的最优Shannon 奇异核褶积微分算子

如何有效协调地震波数值模拟过程中的计算精度和计算效率,一直是许多地球物理工作者密切关注且努力寻求解决的问题.这一问题求解的好坏,将直接影响到人们正确认识地下构造的地球物理响应.本文从离散Shannon奇异核理论出发,推导了基于Shannon奇异核的交错网格褶积微分算子,分析并求出了影响算子精度窗函数的最优化参数.通过将...     

交错网格高阶差分法三维弹性波数值模拟

本文应用交错网格高阶有限差分方法模拟弹性波在三维各向同性介质中的传播。采用时间上二阶、空间上高阶近似的交错网格高阶差分公式求解三维弹性波位移-应力方程,并在计算边界处应用基于傍轴近似法得到的三维弹性波方程吸收边界条件。在此基础上进行了三维盐丘地质模型的地震波传播数值模拟试算。试算结果表明该方法模拟精度高,在很大程度上减小了数值频散,绕射波更加丰富,而且适用于介质速度具有纵向变化和横向变化的情况。     

用交错网格有限差分法计算三维频率域电磁响应

用交错网格有限差分法(SFD)，实现了三维频率域电磁场响应 的数值模拟. 该方法适用于任何方向的磁偶极子源. 经与解析方法、积分方程等 其他方法的计算结果对比表明，交错网格有限差分法结合散度校正和不完全乔累斯基分解预 处理的双共轭梯度迭代方法进行正演计算，速度快、精度高、结果稳定，能适应三维复杂介 质的数值模拟，为三维电磁反演奠定了基础.     

基于双相介质理论的AVO属性交汇图识别天然气水合物

天然气水合物(gas hydrates)作为一种新的潜在能源 是目前世界范围内研究的热点和难点之一。目前 普遍认为拟海底反射(简称BSR)是天然气水合物存 在的重要标志之一,但有关水合物的识别仍有很 多问题需要解决。例如世界各地已发现水合物的 例证表明在没有BSR的情况下仍然存在水合物,在 没有BSR的情况下如何识别水合物仍是没有解决 的研究课题。本文以天然气水合物的识别为研究 目的,以AVO属性交汇图为研究手段,通过双相介 质正演理论模型研究了当地层含有天然气水合物 或游离气时的AVO属性特征,并与由双相介质的 交错网格有限差分法得到的合成地震记录所反演 得出的AVO交汇图进行了对比,得到高度一致的 对比结果,表明利用AVO属性交汇图是识别天然 气水合物和游离气并能估测其含量的有效技术。     

基于改进型TSP观测装置的地震波场数值模拟

随着隧道技术的快速发展,工程上对地质超前预报的可行性和准确性也提出了更高的要求。首先分析了现有以及改进后的定向探测系统的优缺点;然后利用一阶速度—应力弹性波动方程和高阶交错网格法,导出用于改进后的TSP地震波场数值模拟的差分计算方程,并构建倾斜软弱夹层和空洞地质模型,对上述模型进行正演模拟,获得相应的波场快照和单炮记录,进而对地震波场进行分析。结果表明:改进后的TSP观测装置易于实现;有限差分法可有效地模拟改进后的TSP地震波场,其模拟结果清晰可靠地揭示了隧道地质中地震波的传播规律。