煤层地震波场特征分析及波形反演方法研究

煤层地震波场特征分析及波形反演方法研究

论文摘要

煤炭作为我国能源结构的重要支撑之一,在未来社会和经济发展中仍将发挥重要作用。近年来,煤田勘探已逐渐向深部煤层以及复杂构造区域转移,这对精细地震处理解释技术提出了越来越高的要求,而精细地震处理解释技术的实现依赖于地震全波场信息的高效利用,但当前煤田地震勘探依旧集中于利用一次反射波反演方法获取煤层相关介质参数,因此迫切需要发展面向煤田的地震全波场反演方法。全波形反演是一种高精度的定量化地震反演方法,基于波动方程充分利用地震全波场信息,获取到深度域的地层介质参数,其将成为煤田勘探的重要发展趋势,这是本文研究的出发点。全波形反演的基本理念于上世纪八十年代被提出,但受制于计算量大、稳定性差等因素制约。直至近年来计算机软硬件性能不断提升,全波形反演方法才开始在石油勘探领域迅速发展,目前局限于模型测试和海上资料应用。相比于油气储层来说,其一,煤层埋深浅,煤田地震资料易受噪音干扰从而导致地震数据品质较差;其二,煤层是薄低速层,且往往以多套煤互层形式存在,如此导致煤层的地震波形响应规律比油气储层更为复杂;其三,由于岩溶坍塌,煤系地层遭到破坏形成陷落柱,这种小尺度煤层特殊地质构造的反演刻画比常规油气田断裂系统的反演刻画更为困难。因此,将全波形反演技术理念应用到煤田地震勘探时,已然不能直接套用当前针对油气田勘探发展出的全波形反演方法和技术策略,为此面向煤层研究地震波传播规律并研究配套的全波形反演方法技术就成了关键,这也就是本文拟解决的关键问题。在全波形反演中,反演初始模型至关重要,往往由地震构造解释和测井资料插值获取,其正确性决定了全波形反演迭代收敛的速度和稳定性,对煤田弹性波全波形反演来说,其初始模型建立除需要常规纵波测井资料外,还需要横波测井资料,但实际生产中因成本问题而导致横波测井资料较少,往往通过经验关系拟合横波速度,存在较大误差。为提高建模精度,本文从岩石物理测试出发获取到代表性煤层的弹性参数和各向异性参数,并基于人工智能算法建立了基于常规测井参数的横波速度测井曲线重构模型,并验证了其正确性,创新性地提出了基于人工智能的横波速度估计方法,这为煤田地震波场分析和弹性波全波形反演奠定了良好的基础。煤田地震资料全波形反演方法策略的建立依赖于煤田地震波场的正演模拟分析,为此,本文基于粘弹性各向异性波动方程交错网格有限差分正演模拟方法,综合分析了煤层层厚、粘弹性、各向异性、互层以及陷落柱地质构造的地震波场响应特征,研究表明,煤层强反射波形是由顶底界面反射调谐而成,在多套煤互层情况下,调谐影响规律更为复杂,尽管通过时间域波形已经不能有效的分辨煤层,但在不同的频率尺度上,煤层特征会得到不同程度的体现;在一定煤层厚度条件下,粘弹性和各向异性对煤层地震反射波场的影响可以忽略不计;陷落柱等典型煤田地质构造会因断层或者断点绕射使地震波场更为复杂。这些规律和认识为后续煤田地震全波形反演方法的研究提供了方向。考虑到薄互煤系地层反射波形调谐规律的频率响应差异,本文研究了面向薄互层地质结构的多尺度包络-波形反演方法;针对煤田陷落柱断层绕射波等产生的地震波场复杂性,本文提出了基于非精确牛顿优化算法的稳定收敛时间域全波形反演方法;考虑到煤层埋深浅且厚度薄不易估计,本文发展了井约束的正则化全波形反演方法。综合考虑煤田煤系地层特征,将以上三种技术思路互补融合,综合形成了面向煤田勘探的组合优化地震全波形反演方法,模型测试表明该方法能有效的反演刻画出陷落柱构造和煤层构造。该方法为煤层地震反演解释提供了技术支撑,也为今后煤田实际地震资料全波形反演研究奠定了良好的基础。

论文目录

  • 摘要
  • Abstract
  • 1 引言
  •   1.1 研究背景及意义
  •   1.2 国内外研究现状
  •     1.2.1 波动方程正演模拟研究现状
  •     1.2.2 全波形反演方法研究现状
  •     1.2.3 煤田地震勘探现状
  •   1.3 研究内容及技术路线
  •     1.3.1 研究内容
  •     1.3.2 技术路线
  •     1.3.3 主要创新点
  • 2 地震波场正演和全波形反演理论
  •   2.1 波动方程正演理论
  •     2.1.1 弹性波波动方程
  •     2.1.2 波动方程的有限差分求解
  •     2.1.3 震源函数
  •     2.1.4 边界条件
  •     2.1.5 稳定性条件
  •     2.1.6 数值频散
  •   2.2 全波形反演理论
  •     2.2.1 全波形反演基本流程
  •     2.2.2 非线性反演理论
  •     2.2.3 伴随状态法求取梯度
  •     2.2.4 迭代步长的求取
  •     2.2.5 Hessian矩阵预处理方法
  •   2.3 时间域全波形反演中的关键问题
  •     2.3.1 低频波场
  •     2.3.2 初始模型
  •   2.4 MPI并行计算
  •     2.4.1 MPI并行模式
  •     2.4.2 区域分解下的并行框架
  •     2.4.3 并行试算分析
  •   2.5 本章小结
  • 3 基于岩石物理和人工智能的煤层弹性参数获取
  •   3.1 岩石物理测试
  •     3.1.1 煤岩岩石物理实验
  •     3.1.2 煤岩的弹性参数及各向异性参数提取
  •   3.2 基于测井参数的横波速度预测
  •     3.2.1 人工智能理论与方法
  •     3.2.2 横波速度测井曲线构建
  •     3.2.3 模型试算
  •   3.3 本章小结
  • 4 煤层地震波场数值模拟
  •   4.1 粘弹各向异性弹性波动方程
  •   4.2 薄互煤层调谐波形响应特征分析
  •     4.2.1 煤层层厚影响分析
  •     4.2.2 薄互煤层影响分析
  •   4.3 粘弹各向异性煤层地震波场特征分析
  •     4.3.1 各向异性煤介质地震波场特征分析
  •     4.3.2 粘弹煤介质地震波场特征分析
  •   4.4 煤层隐蔽致灾模型地震波场特征分析
  •   4.5 本章小结
  • 5 井约束时间域多尺度包络-波形反演
  •   5.1 基于非精确牛顿优化算法的全波形反演
  •   5.2 基于包络-波形的全波形反演
  •     5.2.1 包络-波形反演方法
  •     5.2.2 包络数据提取
  •     5.2.3 Marmousi2模型测试
  •   5.3 基于低通滤波技术的多尺度全波形反演
  •     5.3.1 Butterworth低通滤波器
  •     5.3.2 基于低通滤波的多尺度反演方法
  •     5.3.3 Marmousi2模型测试
  •   5.4 基于正则化的井约束全波形反演
  •     5.4.1 井数据约束基本方法
  •     5.4.2 井约束全波形反演策略
  •     5.4.3 Marmousi2模型测试
  •   5.5 组合优化方法的全波形反演
  •   5.6 基于陷落柱模型的全波形反演
  •   5.7 本章小结
  • 6 结论与展望
  •   6.1 结论
  •   6.2 展望
  • 参考文献
  • 致谢
  • 作者简介
  • 文章来源

    类型: 博士论文

    作者: 王攀

    导师: 彭苏萍,杜文凤

    关键词: 煤田地震勘探,人工智能,时间域弹性波全波形反演,非精确牛顿,多尺度

    来源: 中国矿业大学(北京)

    年度: 2019

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

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

    单位: 中国矿业大学(北京)

    分类号: P631.4;P618.11

    DOI: 10.27624/d.cnki.gzkbu.2019.000192

    总页数: 210

    文件大小: 21665K

    下载量: 76

    相关论文文献

    标签:;  ;  ;  ;  ;  

    煤层地震波场特征分析及波形反演方法研究
    下载Doc文档

    猜你喜欢