基于弹性波全波形反演的煤层异常体成像研究

发布时间：2020-11-11 14:46

　　 煤炭工业是我国国民经济的重要支柱,煤炭消费在一次性能源结构中比重占60%以上。煤炭生产的地质条件关系到煤矿安全生产,而煤层中地质异常是诱发煤矿安全事故的直接因素。随着多年的煤炭生产,煤炭开采条件逐渐变得复杂和困难,常规地震勘探方法难以满足高效安全开采的地质需求。因此,开展高精度、高分辨率的勘探方法的研究,对实现煤矿井下精细化探测具有重要意义。槽波地震勘探方法是一种较为成熟方法,目前被广泛应用于煤层工作面的探测中。虽然槽波勘探方法具有其它勘探方法无法比拟的优势和精度,但对于煤层小尺度异常以及其它异常体仍然表现不佳或无能为力。地震全波形反演作为一种高精度、高分辨率的地下介质物性参数反演方法,是近年来地震成像领域研究的热点之一。全波形反演利用了地震记录中全部的运动学和动力学信息,因而具有精细刻画复杂地质构造的能力。为解决煤矿井下精细化探测的问题,本文引进了全波形反演方法,探讨在煤层槽波条件下反演的可行性与适应性。通过对波动方程有限差分正演、自由边界条件、弹性波全波形反演基本理论、梯度计算、步长求取、GPU并行加速、煤层常见异常体反演等方面进行详细探讨,取得以下主要成果:(1)通过对传统二维起伏地表网格分类方案进行改进,实现了三维模型隧道或巷道全方位网格点分类策略,解决了含煤层与巷道模型的巷道壁面波的产生问题。(2)设计了基于单GPU的二维弹性波全波形反演实现方案,采用了波场重建方法解决了GPU显存不足的问题;对不同线程块维度配置的计算效率进行对比分析,找到了影响并行算法计算效率的关键因素;利用基于CPU并行和GPU并行的算法进行对比研究,分别得到单GPU卡相对于单CPU与八核CPU正演计算和梯度计算的并行加速比。(3)设计了基于多节点多GPU的三维弹性波全波形反演实现方案。实现了基于模型域分解的并行策略,每一个子模型单独利用一个GPU卡计算,子模型之间通过MPI进行数据通信,从而实现了整个反演过程的双重并行;实现了利用二维线程网格计算三维模型的并行计算策略,以解决线程块共享内存不足的问题;通过对子模型边界区域和内部区域先后计算,采用多个CUDA流共同作业,将子模型边界区域节点间数据交流时间隐藏在计算中,从而实现了GPU粗粒度并行;通过对多GPU并行计算与单GPU并行计算进行对比测试,得到了多GPU相对于单GPU并行加速效率。(4)对常见煤层异常体模型进行理论反演发现,全波形反演算法对煤层槽波的存在具有较好的适应性和稳定性,反演能够较好地恢复陷落柱、断层块、夹矸层以及煤层弯曲构造,反演结果达到较高的分辨率;对比反演结果也能够认识到,由于煤层的存在会对地震波的传播产生限制,对于煤层中的异常体能够较精确地反演,而对于围岩中异常体,距离煤层越近反演分辨率越高,距离越远分辨率越低。(5)观测系统和观测方式对全波形反演结果有一定影响。通过对煤层中不同观测系统的布设对反演结果的影响的研究表明,要想获得更好的反演效果和保证较好的纵向(煤厚方向)分辨率,需要在煤层纵向上布设多条炮线或测线,简单地说,炮点和检波点越多,空间分布范围越大,反演效果越好。本论文的创新之处主要在于:(1)提出将三维弹性波全波形反演方法引进到煤层槽波条件下地质异常的探测中。煤层中断层、夹矸、陷落柱、顶底板起伏等各种异常体理论模型的反演结果表明,全波形反演能够适应煤层槽波条件,反演结果较精确地刻画三维小尺度构造。全波形反演方法能够为煤层的精细化探测提供新的解决方案。(2)发展了传统地表自由表面网格点分类方案,使之能够适用于三维隧道或巷道全方位自由表面网格点模型。(3)设计了基于多节点多GPU的三维弹性波全波形反演并行实现方案。方案中实现了线程级并行和任务级并行组合,线程块优化配置,共享内存、寄存器等资源充分利用,最大程度发挥了GPU计算性能,极大地提高了三维反演问题的计算效率。
【学位单位】：中国地质大学
【学位级别】：博士
【学位年份】：2018
【中图分类】：P631.4
【文章目录】：
作者简历
摘要
abstract
第一章 绪论
    1.1 选题的目的和意义
    1.2 国内外研究现状及存在问题
        1.2.1 槽波勘探技术研究现状
        1.2.2 全波形反演研究现状
        1.2.3 有限差分中GPU并行计算方法研究现状
        1.2.4 存在的问题
    1.3 主要研究内容
    1.4 论文的主要成果及创新点
    1.5 论文的章节安排
第二章 弹性波方程正演模拟
    2.1 弹性波方程及数值离散
    2.2 数值频散及稳定性条件
    2.3 CPML边界条件
    2.4 自由边界条件
    2.5 本章小结
第三章 全波形反演理论方法
    3.1 非线性反演理论
    3.2 梯度计算
    3.3 步长估计
    3.4 多尺度反演
    3.5 波场重建
    3.6 本章小结
第四章 基于GPU的全波形反演算法并行实现
    4.1 CUDA编程模型
    4.2 单GPU2D全波形反演实现及性能分析
        4.2.1 并行实现
        4.2.2 性能分析
    4.3 多GPU3D全波形反演实现及性能分析
        4.3.1 并行实现
        4.3.2 性能分析
    4.4 2D弹性波全波形反演实例
    4.5 3D弹性波全波形反演实例
    4.6 本章小结
第五章 煤层地质异常体成像
    5.1 煤层中主要地质异常体
    5.2 槽波的形成
    5.3 槽波波场特征
    5.4 煤层地质异常体反演
        5.4.1 模型及观测系统设置
        5.4.2 陷落柱模型反演
        5.4.3 断层块模型反演
        5.4.4 夹矸模型反演
        5.4.5 煤层起伏模型反演
    5.5 本章小结
第六章 结论与建议
    6.1 结论
    6.2 建议
致谢
参考文献

【相似文献】

相关期刊论文 前10条

1 本刊编辑部;;写在“全波形反演”专题前面的话[J];石油物探;2017年01期

2 高福春;Paul WILLIAMSON;R.Gerhard PRATT;;全波形反演的一个新目标函数:数据域中的微分相似优化[J];石油物探;2017年01期

3 张广智;姜岚杰;孙昌路;黄义双;;基于照明预处理的分步多参数时间域声波全波形反演方法研究[J];石油物探;2017年01期

4 魏哲枫;朱成宏;陈业全;;三维全波形反演高效异构并行计算[J];石油物探;2017年01期

5 周斯琛;李振春;张敏;张凯;;基于截断牛顿法的频率域全波形反演方法[J];物探与化探;2017年01期

6 柯睿;杨艳;;多频点同时反演策略在频域全波形反演中的应用和实现[J];半导体光电;2017年01期

7 张文生;罗嘉;;时间域两网格全波形反演[J];数值计算与计算机应用;2016年01期

8 夏冬明;宋鹏;谭军;李金山;张晓波;钟梦轩;;基于归一化局部互相关算法的相位编码全波形反演[J];中国海洋大学学报(自然科学版);2017年07期

9 王庆;张建中;黄忠来;;时间域地震全波形反演方法进展[J];地球物理学进展;2015年06期

10 孟旭;刘四新;吴俊军;傅磊;邓莉;;时间域跨孔雷达全波形反演及实际应用[J];世界地质;2016年01期

相关博士学位论文 前10条

1 蒋锦朋;基于弹性波全波形反演的煤层异常体成像研究[D];中国地质大学;2018年

2 雷涛;时间域全波形反演方法及最优化理论研究[D];中国地质大学;2017年

3 陈汉明;波动方程数值模拟与粘滞波形反演方法研究[D];中国石油大学（北京）;2017年

4 张庆臣;时间域弹性波动方程全波形反演方法研究[D];中国石油大学（北京）;2017年

5 刘玉金;复杂介质条件下基于反演的地震成像方法研究[D];中国石油大学(华东);2014年

6 任乾慈;基于全波形反演的探地雷达数据和地震数据交叉梯度联合反演[D];吉林大学;2017年

7 陈国新;地震数据波形反演方法理论研究[D];浙江大学;2018年

8 段超然;基于初始模型优化的稀疏高效全波形反演方法[D];吉林大学;2018年

9 王志强;非光滑约束的高精度全波形反演方法研究[D];吉林大学;2018年

10 韩淼;深水区地震全波形反演策略与应用[D];吉林大学;2014年

相关硕士学位论文 前10条

1 张天泽;基于小波变换包络的时间域全波形反演研究[D];吉林大学;2018年

2 肖增佳;VTI介质弹性波全波形反演[D];东北石油大学;2018年

3 李松龄;时域全波形反演及其敏感性分析研究[D];东北石油大学;2018年

4 向世明;基于边界引导的全波形反演算法[D];中国科学技术大学;2018年

5 胥洋;基于快速推进法（FMM）的层析反演研究[D];中国石油大学（北京）;2017年

6 尹畅;基于包络数据的全波形反演实现技术[D];吉林大学;2018年

7 李铭;基于频率域全波形反演的隧道前方岩体波速预测方法研究[D];山东大学;2018年

8 韩雨桐;起伏地表声波全波形反演方法研究[D];中国石油大学(华东);2016年

9 孙昌路;面向储层的全波形反演方法研究[D];中国石油大学(华东);2016年

10 尹力;低频缺失下全波形反演应用研究[D];中国石油大学(华东);2015年

本文编号：2879326