波动方程叠前深度成像理论和应用的研究.pdf

摘要 摘要 波动方程叠前深度偏移成像是复杂介质成像的有效手段。波动方程叠前深 度偏移理论方法已逐渐成熟。但在解决横向强变速问题时，仍存在许多技术难 题。尤其是陆上数据的叠前深度成像存在更多的技术瓶颈。 波动方程叠前深度偏移的最核心工作是利用波场外推算子进行波场深度外 推。基于摄动理论的混合域单平方根传播算子都是对单平方根方程的近似，存 在近似误差。各传播算子都涉及到参考速度的选择。近似条件要求局部参考速 度与实际速度相差不能太大。否则在强横向变速区域，会出现计算的不稳定， 计算精度降低。因此，采用多个参考速度是必要的。本文针对混合域传播算子 的特点，提出自适应多参考速度选择策略，与混合域波场外推算子结合，形成 适应强横向变速的混合域波场外推算法。 单平方根有限差分算子的数值计算中存在高角度成像误差、强横向变速情 况下可能出现数值计算的不稳定及计算误差大等固有缺陷。三维问题的方向分 裂有限差分计算会引入数值各向异性。本章在理论研究的基础上，总结前人的 方法优势，提出四程分裂有限差分解法加有限差分补偿的波场外推方法，并完 成了详细的理论推导。它可以解决剧烈横向变速及复杂构造情况下的地震波成 像问题，认为是目前最完善的空间频率域三维有限差分波场外推方法。 理论上，全三维双平方根波动方程波场外推进行全三维叠前深度偏移是可 行的。但由于受到目前计算机系统能力的限制，还无法用于实际。共方位角双 平方根波场外推算子，是全三维双平方根波动方程稳相近似得到的。它拥有较 高的计算效率，但由于稳相近似降维处理，波场外推精度和对横向变速适用能 力有所降低。通过研究提出适应于强横向变速介质的自适应多参考速度选择中 点一炮检距域广义屏波场外推算法和有限差分补偿的广义屏算法，并进行了系数 的优化。 实验表明，自适应多参考速度选择策略的混合域各种波场外推算子、四程 分裂有限差分解法加有限差分补偿的波场外推方法和三维共方位角双平方根广 义屏波场外推方法，均可提高对横向变速的适应能力，有效地提高成像精度。 波动方程叠前深度偏移的应用是有条件的。为促进其实际应用的发展，针 摘要 对性地研究了波动方程叠前深度偏移密切相关的并行计算和数据规则化等配套 技术。 单平方根波动方程波场外推和共角度域叠前深度偏移计算量十分巨大，必 须借助计算速度高的微机集群和并行计算技术，才能促进其实际应用。在研究 波动方程波场外推算法的基础上，提出了适合多任务多处理器的分布式微机集 群处理系统主从模式负载动态平衡并行算法。目前的微机集群并行处理系统， 没有断点恢复功能，作业必须从头开始。系统级断点恢复技术，是世界级的研 究难题。本文通过研究提出软件级断点恢复方法，在实际应用中获得良好的应 用效果。 野外数据采集方式是数据不规则产生的原因，不可避免的波动方程叠前深 度偏移处理造成影响。在数据规则化理论研究的基础上，分类总结已有的数据 规则化方法。重点研究基于模型驱动的数据规则化方法和基于数据驱动的Radon 变换道插值方法。提出反假频积分法正DMO加f-k反DMO方位角校正方法，软件 研发及实际应用。 最后，本文从理论程序的软件化、处理流程建立和实际资料处理等几方面， 进行了应用研究。将单平方根混合域炮道集叠前深度偏移程序软件化，建立单 平方根混合域炮道集叠前深度偏移及相关处理流程，应用自主研发的处理系统 软件--STseis，在国内首次用三维炮域波动方程叠前深度偏移方法，完成胜利 探区140平方公里的实际资料处理，获得较好的处理效果。同时将三维共方位 角波动方程深度偏移算法纳入了STseis软件系统，建立了合理的三维共方位角 波动方程深度偏移处理流程，在国内首次用波动方程深度偏移技术完成了130 平方公里以上的三维地震资料处理，也取得了明显的应用效果。 本文的主要创新点： 1．提出自适应多参考速度选择策略，改进混合域单平方根和双平方根波动 方程波场外推算子。 2．提出主从模式负载平衡动态任务分配并行算法和软件级断点恢复方法。 3．将研究成果软件化，形成具有自主知识产权的叠前成像软件系统的核心 模块。在国内首次用波动方程叠前偏移技术完成300平方公里的三维地震资料 的处理，使自主创新技术实用化。 关键词：波动方程叠前深度偏移，单平方根算子，双平方根算子，并