三维复杂介质的单程波算法

摘要

随着勘探目标区变的更复杂,地震勘探方法必须考虑地球介质的三维非均匀影响;而为了更好发展三维地震资料处理和成像方法,必须正确认识地震波在复杂的三维地下构造的传播路径及运动学和动力学特征.只有这样准确地分析油气储集体所产生的反射地震波场特性,才能更有效地依据地震记录对储层进行识别和圈划.有限差分法可容易处理介质纵向和横向非均匀,是一种较有效的复杂介质构造的地震波模拟方法,但就现行计算机的计算速度和内存而言,尚不能应用三维有限差分法解决勘探规模的地震波模拟问题.因此必须发展新的高效的地震波模拟方法.现行的地震资料处理主要是基于分析和利用一次反射波,而多次波和表面波等其它波相均被视为噪音.因此就研究地下油气储集体的地震响应及验证地质建模而言,正确地模拟复杂介质构造的一次反射波已满足需要.该文应用波动方程的单程波算法发展了一个三维复杂构造的地震波模拟方法.这一方法可正确地模拟地震勘探中最关心的一次反射波.它最重要的应用是生成零偏移距叠加剖面;在这一计算中该方法所需的计算量比双程波方法计算一个单炮计录还要少许多倍.二维复杂构造下双程波有限差分法的单炮记录与该文单程波模拟结果的比较表明,该文建议的单程波算法有很高的精度.这一方法在程序实现上可大量地借用现行的深度偏移方法和程序系统,易于实现并行计算.除了所需计算量很小外,这一方法对速度模型不象双程波算法那样敏感,而这一点对地质模型的地震波模拟尤为重要.因为实际上我们仅能得到近似的地下速度构造,由较小的速度差异而引出的较多差异的地震响应将给判断地质建模优劣带来困难.基于计算效率和计算精度的双重考虑,该文在介质速度结构较复杂射采用精度较高的显式短算子波场延拓方法,而在介质速度结构相对简单时采用计算效率较高的分裂步相移方法,建议了一个将结合显式短算子法与分裂步相移方法的三维单程波计算策略.文中进一步发展了显示短算子的单程波算法,提出了一个加权二次规划的显式短算子设计方法,这一算子设计方法更好地保证了显示短算子算法的稳定性,使得可用更短的算子取得较高的精度,从而使三维显式短算子单程波算法有更高的计算效率和精度.该文算法已应用于实际地震勘探工程,完成的东部某油田地质模型的三维地震模拟为解释人员正确认识该地区的地震特征很好的帮助.

目录

文摘

英文文摘

0前言

1弹性波理论在地震勘探中的应用

1.1地震勘探方法简述

1.2地震波数值模拟及其应用

1.3地震波成像的偏移方法

1.4本文研究的重点和章节的安排

2地震波传播的单程波和双程波理论

2.1声波和弹性波方程

2.2声波方程的单程波近似

3地震波数值模拟

3.1地震波数值模拟方法简述

3.2基于单程波算法的地震波数值模拟方法

3.2.1叠加剖面计算

3.2.2炮记录计算

3.3界面反射系数计算

3.3.1入射角度的确定

3.3.2反射系数计算

4混合的显式短算子和分裂步相移单程波算法

4.1二维显式短算子单程波算法

4.1.1显式短算子单程波算法

4.1.2一维显式短算子的设计

4.1.3显式短算子方法的程序实现及算子表的建立和应用

4.2分裂步相移单程波算法

4.3三维显式短算子单程波算法

4.3.1轴对称算子及McClellan变换

4.3.2直接利用二维显示算子进行三维深度延拓

4.4混合的显式短算子和分裂步相移方法及其实现策略

5数值算例及讨论

5.1单程波与双程波方法结果比较

5.2三维地质模型的模拟结果

5.2.1三维地质模型的零偏移距叠加剖面计算

5.2.2特定储层的零偏移距叠加剖面计算

6结论与展望

参考文献

致谢