时域全波形反演及其敏感性分析研究
发布时间：2020-10-25 04:21
获得准确的地震波传播速度是地震资料处理中的重要环节之一,速度参数的精准性对高精度地震偏移成像的刻画起到了至关重要的作用。常规确定速度参数的方法有速度谱扫描、射线走时层析技术等。然而,上述常规方法只能得到一个粗略的速度模型,影响了偏移成像的精度。全波形反演是基于双程波动方程,通过匹配接收记录中的全波形数据重建地下参数的方法,该方法能重建出高精度的速度模型。但是,全波形反演是一个计算量巨大、高度病态的非线性反问题。尤其需要高精度高效率的波场模拟和富含低频信息的初始速度模型。针对高精度波场模拟问题,本文首先从波场模拟的精度和效率出发,推导了基于双程波动方程的高阶有限差分的差分系数;分析了波场模拟中的稳定性条件和频散关系;为消除由于人工截断而引起的边界反射问题,对比分析了Clayton吸收边界和完全匹配吸收(PML)边界的吸收效果;研究了基于GPU/CPU并行加速的波场模拟策略,再经存储器和归约算法优化后,大幅度的提升了全波形反演算法的计算效率。全波形反演的梯度计算需要求取所有时刻的瞬时震源波场和伴随反传的检波点波场,存储量巨大。为降低波场存储量,本文研究了仅存储部分震源波场的波场重建策略,与存储所有震源波场相比该策略可有效降低百分之九十的存储需求,缓解了全波形反演算法中波场存储的巨大压力;为了缓解地震资料中低频信息的不足的问题,研究了一种基于时间二阶积分的全波形反演策略。该策略能有效地提高地震资料中的低频信息比重。在一定程度上避免了由于初始模型不准确引起的周波跳跃现象,改善了速度反演的精度,并利用Marmousci2模型证明了该方法的正确性;最后利用了多组数值实验,验证了初始模型、震源子波、观测方式等敏感性因素对全波形反演结果的影响。
【学位单位】：东北石油大学
【学位级别】：硕士
【学位年份】：2018
【中图分类】：P631.4
【部分图文】：
动方程差分稳定性条件，对比分析了 Clayton 吸收边界和完全匹配吸收（PML）边界对边界反射的吸收效果；分析了波场模拟中参数选择对频散的影响；（2）研究了 GPU/CPU 协同并行加速波场模拟策略和提高矩阵求和效率的归约算法；（3）介绍了全波形反演的基本理论，研究了时间域全波形反演算法的流程。推导了全波形反演中梯度、步长等关键参数；利用保存部分波场的波场重建策略重建了全部震源波场，并利用数值实验证明了波场重建的准确性；（4）应用水平层状模型检验了时间域全波形反演方法的实现流程，应用 Marmousi2模型测试了时间域全波形反演方法的有效性和稳定性，并分析了时间域域全波形反演的反演结果与迭代次数的关系；（5）研究了基于时间二重积分的全波形反演方法；（6）对全波形反演中的初始模型、震源子波、观测方式等敏感性因素进行了测试和分析。本文的技术路线如图 0.1 所示。
第一章 波动方程正演模拟2 222 22( , , ) 2 ( , , ) ( , , ) ( , , )( , , )Nii NNii Nv tP x z i t P x z t P x z t C P x i x z txv tC P x z i z tz（为验证有限差分算法波场延拓的正确性，本文采用上述公式对常速度均匀介质波场模拟计算，其正演参数如下：在模型中央使用主频为 20Hz 的 ricker 子波激速度为 2000m/s，水平和垂直维度计算点数目均为 400，水平和垂直空间网格 10m，时间方向采样率为 1ms。采用空间 10 阶差分算法。试算结果如图 1. 1.1a 为地震波传播至 0.8s 时刻的波场快照，图 1.1b 为地震波传播至 1.2s 时刻照）：
削弱由于设置人工边界而产生的实际测量的真实解。 RobertClayton 和 Engquist[51]提paraxialapproximations，PA），只模拟向边界外侧传播的地震波震波的入射角度，当地震波入收边界计算公式为：2 2 22 212p p v px t v t z 可以写成相似于（1-10）的形式.2。图中白色区域为双程波动方单程波场方程计算区域其上下左
【参考文献】
本文编号：2855454
【学位单位】：东北石油大学
【学位级别】：硕士
【学位年份】：2018
【中图分类】：P631.4
【部分图文】：
动方程差分稳定性条件，对比分析了 Clayton 吸收边界和完全匹配吸收（PML）边界对边界反射的吸收效果；分析了波场模拟中参数选择对频散的影响；（2）研究了 GPU/CPU 协同并行加速波场模拟策略和提高矩阵求和效率的归约算法；（3）介绍了全波形反演的基本理论，研究了时间域全波形反演算法的流程。推导了全波形反演中梯度、步长等关键参数；利用保存部分波场的波场重建策略重建了全部震源波场，并利用数值实验证明了波场重建的准确性；（4）应用水平层状模型检验了时间域全波形反演方法的实现流程，应用 Marmousi2模型测试了时间域全波形反演方法的有效性和稳定性，并分析了时间域域全波形反演的反演结果与迭代次数的关系；（5）研究了基于时间二重积分的全波形反演方法；（6）对全波形反演中的初始模型、震源子波、观测方式等敏感性因素进行了测试和分析。本文的技术路线如图 0.1 所示。
第一章 波动方程正演模拟2 222 22( , , ) 2 ( , , ) ( , , ) ( , , )( , , )Nii NNii Nv tP x z i t P x z t P x z t C P x i x z txv tC P x z i z tz（为验证有限差分算法波场延拓的正确性，本文采用上述公式对常速度均匀介质波场模拟计算，其正演参数如下：在模型中央使用主频为 20Hz 的 ricker 子波激速度为 2000m/s，水平和垂直维度计算点数目均为 400，水平和垂直空间网格 10m，时间方向采样率为 1ms。采用空间 10 阶差分算法。试算结果如图 1. 1.1a 为地震波传播至 0.8s 时刻的波场快照，图 1.1b 为地震波传播至 1.2s 时刻照）：
削弱由于设置人工边界而产生的实际测量的真实解。 RobertClayton 和 Engquist[51]提paraxialapproximations，PA），只模拟向边界外侧传播的地震波震波的入射角度，当地震波入收边界计算公式为：2 2 22 212p p v px t v t z 可以写成相似于（1-10）的形式.2。图中白色区域为双程波动方单程波场方程计算区域其上下左
【参考文献】
相关期刊论文 前10条
1 王庆;张建中;黄忠来;;时间域地震全波形反演方法进展[J];地球物理学进展;2015年06期
2 黄超;董良国;迟本鑫;;低频缺失情况下的弹性波包络反演（英文）[J];Applied Geophysics;2015年03期
3 敖瑞德;董良国;迟本鑫;;不依赖子波、基于包络的FWI初始模型建立方法研究[J];地球物理学报;2015年06期
4 张文生;罗嘉;滕吉文;;频率多尺度全波形速度反演[J];地球物理学报;2015年01期
5 刘洋;;波动方程时空域有限差分数值解及吸收边界条件研究进展[J];石油地球物理勘探;2014年01期
6 陈永芮;李振春;秦宁;张凯;;波动方程双向照明优化的全波形反演[J];地球物理学进展;2013年06期
7 高凤霞;刘财;冯晅;鹿琪;王典;;几种优化方法在频率域全波形反演中的应用效果及对比分析研究[J];地球物理学进展;2013年04期
8 胡光辉;贾春梅;夏洪瑞;贺剑波;宋林;沈忠秋;;三维声波全波形反演的实现与验证[J];石油物探;2013年04期
9 王保利;高静怀;陈文超;张唤兰;;地震叠前逆时偏移的有效边界存储策略[J];地球物理学报;2012年07期
10 刘国峰;刘洪;孟小红;闫浩飞;;频率域波形反演中与频率相关的影响因素分析[J];地球物理学报;2012年04期
相关硕士学位论文 前1条
1 宫启明;基于非精确全牛顿算法VTI介质多参数全波形反演[D];吉林大学;2017年
本文编号：2855454
本文链接：https://www.wllwen.com/kejilunwen/diqiudizhi/2855454.html
