粘弹垂向各向异性介质地震波场的交错网格数值模拟方法

粘弹垂向各向异性介质地震波场的交错网格数值模拟方法

论文摘要

地震波场数值模拟方法不仅是地球探测的重要基础,也是地球物理勘探的重要组成部分。它具有描述模型简便,参数选择灵活的优点,对研究地震波传播的动力学和运动学规律起着重要作用。当已知地下介质的参数和结构时,可以通过数值计算的方法研究分析地震波在地下介质中的传播规律,从而得到不同时刻的波场快照以及合成地震记录,对全面深入地了解地震波传播规律及传播机理都具有重要意义。波场正演数值模拟方法主要分为三类,其中波动方程求解法是最为重要的一种方法。采用该方法得到的模拟结果可以反映地震波传播的动力学和运动学信息,为研究复杂地质构造的波场提供了基础。该方法在对地震资料进行采集、处理及解释的过程中也起到了重要作用,因此得到了广泛应用。有限差分法是其中一种较为经典的方法,具有计算速度快、内存占用小、精度高等优点。其原理就是将地下地质体简化为一个模型,将该模型区域用有限个网格进行划分;然后使用离散空间和时间来替换连续物理量,用差分和微分的近似关系,并将微分形式的波动方程改变为适合于计算的差分方程;然后给定一个初始条件,利用在空间中前时刻的波场值的分布和变化来推导现在时刻的波场值,从而模拟地震波的传播。在进行数值模拟的过程中,我们一般会假设地下介质是完全弹性的各向同性介质,然而实际的地球介质并不是理想的弹性体,它具有一定的粘性;地震波在介质中传播时能量会被吸收,同时,地震波的高频成分会衰减,振幅也会减小,以至于我们不能直接地得到详细的地下信息。另外,在研究地震传播问题时,通常以均匀各向同性介质来近似实际介质,这种简化会引起相当大的误差。因此我们模拟时采用了粘弹各向异性介质模型来模拟,这样才能更加客观准确。本文以波动理论为基础,采用的是交错网格高阶有限差分数值模拟法,得到适合于各波动方程的高阶交错网格有限差分格式。模拟时采用完全匹配层(PML)吸收边界条件,最终实现了对弹性和粘弹性各向同性介质、VTI介质、HTI介质及TTI介质的正演模拟。同时,通过对比不同品质因子下模拟结果的异同,分析Q值与介质粘弹性的关系。另外,本文对层状介质模型和断层介质模型进行了正演模拟,得到的模拟结果包含丰富的波动信息。模拟结果证明了该模拟方法的正确性和有效性,精度较高,可靠性强。

论文目录

  • 摘要
  • abstract
  • 第1章 绪论
  •   1.1 选题背景与研究意义
  •   1.2 粘弹介质地震波场模拟国内外研究现状
  •   1.3 地震波有限差分算法国内外研究现状
  •   1.4 论文的主要内容
  •   1.5 本章小结
  • 第2章 各向异性介质波场模拟理论
  •   2.1 各向同性介质弹性波动方程
  •     2.1.1 位移波动方程
  •     2.1.2 一阶速度-应力波动方程
  •   2.2 各向异性介质基本理论及分类
  •   2.3 各向异性介质弹性波动方程
  •   2.4 本章小结
  • 第3章 各向异性粘弹性波动理论
  •   3.1 粘弹性介质基本性质
  •   3.2 粘弹性介质的两大基本模型
  •     3.2.1 Maxwell模型
  •     3.2.2 kelvin模型
  •   3.3 地层吸收衰减相关参数
  •     3.3.1 吸收系数
  •     3.3.2 对数衰减量
  •     3.3.3 衰减系数
  •     3.3.4 品质因子
  •   3.4 粘弹性介质波动方程
  •     3.4.1 粘弹性各向同性介质波动方程
  •     3.4.2 粘弹性TTI介质波动方程
  •   3.5 本章小结
  • 第4章 交错网格有限差分法
  •   4.1 空间上的2M阶差分近似
  •   4.2 时间上的2N阶差分近似
  •   4.3 交错网格差分格式
  •   4.4 本章小结
  • 第5章 地震波场模拟的初边值条件
  •   5.1 初始条件
  •   5.2 震源函数的选择
  •   5.3 稳定性的分析
  •   5.4 完全匹配层吸收边界条件
  •   5.5 本章小结
  • 第6章 典型介质模型下的数值结果及分析
  •   6.1 弹性及粘弹性各向同性介质波场数值模拟
  •     6.1.1 粘弹性各向同性介质有限差分格式
  •     6.1.2 弹性及粘弹性各向同性介质模型试算及对比
  •   6.2 弹性及粘弹性VTI介质波场数值模拟
  •     6.2.1 弹性VTI介质有限差分格式及模型试算
  •     6.2.2 粘弹性VTI介质模型试算
  •   6.3 粘弹性HTI介质波场数值模拟
  •   6.4 弹性TTI介质波场数值模拟
  •     6.4.1 均匀模型
  •     6.4.2 层状模型
  •   6.5 粘弹性TTI介质波场数值模拟
  •     6.5.1 均匀模型
  •     6.5.2 断层模型
  •   6.6 本章小结
  • 第7章 结论与展望
  • 参考文献
  • 作者简介及在学期间所取得的科研成果
  • 致谢
  • 文章来源

    类型: 硕士论文

    作者: 隋竞函

    导师: 刘财

    关键词: 完全匹配层,交错网格,有限差分,粘弹性介质,各向异性介质

    来源: 吉林大学

    年度: 2019

    分类: 基础科学,工程科技Ⅰ辑

    专业: 地质学,地球物理学,矿业工程

    单位: 吉林大学

    分类号: P631.4

    总页数: 74

    文件大小: 2957K

    下载量: 159

    相关论文文献

    • [1].各向异性介质中的电磁辐射场[J]. 咸阳师范学院学报 2019(02)
    • [2].电各向异性介质中带电圆环的电势和场强[J]. 广西物理 2014(04)
    • [3].电各向异性介质中泊松方程的具体形式[J]. 泉州师范学院学报 2015(02)
    • [4].格林函数法求电各向异性介质中的静电势[J]. 内江师范学院学报 2015(12)
    • [5].各向异性介质中大地电磁正演研究综述[J]. 地球物理学进展 2011(06)
    • [6].线电荷在分层电各向异性介质中激发的电势[J]. 郑州师范教育 2015(06)
    • [7].各向异性介质转换波速度逐层单参数分析方法[J]. 物探与化探 2011(02)
    • [8].各向异性介质圆柱的散射特性研究[J]. 装备环境工程 2011(06)
    • [9].各向异性介质中S波接收函数成像方法[J]. 地球物理学报 2020(05)
    • [10].电偶极子的各向异性辐射特性[J]. 装备环境工程 2017(07)
    • [11].任意各向异性介质相(群)速度的计算[J]. 地球物理学报 2012(10)
    • [12].各向异性介质椭球的瑞利散射特性[J]. 重庆理工大学学报(自然科学) 2011(06)
    • [13].磁各向异性介质中圆电流外部的磁场[J]. 海南大学学报(自然科学版) 2011(02)
    • [14].各向异性介质椭球内电场的研究[J]. 装备环境工程 2010(04)
    • [15].非椭圆各向异性介质高阶优化走时计算方法及应用[J]. 地球物理学进展 2012(03)
    • [16].基于多尺度理论的电各向异性介质椭球内电场[J]. 强激光与粒子束 2011(04)
    • [17].层状各向异性介质转换波克希霍夫叠前时间偏移[J]. 地球物理学报 2009(12)
    • [18].水平层状电各向异性介质大地电磁正演研究[J]. 物探化探计算技术 2015(03)
    • [19].各向异性介质岩石物理模型方法研究[J]. 化工管理 2016(29)
    • [20].应力诱导各向异性介质中地震波的传播[J]. 华北地震科学 2012(03)
    • [21].弹性各向异性介质岩石体波相速度计算[J]. 地震地磁观测与研究 2011(06)
    • [22].旋转轴对称弱各向异性介质模型的射线追踪方法[J]. 地球物理学进展 2013(06)
    • [23].三维各向异性介质中的波动方程叠前深度偏移方法[J]. 地球物理学报 2011(11)
    • [24].涂覆单轴双各向异性介质圆柱的电磁散射特性[J]. 电脑知识与技术 2017(31)
    • [25].三维方位各向异性介质数值模拟及波场分析[J]. 世界地质 2010(01)
    • [26].各向异性介质多波速度分析方法及软件实现[J]. 物探化探计算技术 2010(04)
    • [27].复杂各向异性介质初至波射线追踪[J]. 石油物探 2016(01)
    • [28].各向异性介质断层震源等效体力[J]. 地壳构造与地壳应力文集 2011(00)
    • [29].强各向异性介质的微地震数据反演的陷阱(英文)[J]. Applied Geophysics 2016(02)
    • [30].大地电磁各向异性模型解析解及其讨论[J]. 中国科学:地球科学 2013(09)

    标签:;  ;  ;  ;  ;  

    粘弹垂向各向异性介质地震波场的交错网格数值模拟方法
    下载Doc文档

    猜你喜欢