波动方程地震偏移成像方法的现状与进展.pdfVIP

维普资讯 第 17卷 第 4期 地 球 物 理 学 进 展 V01．17 No．4 2002年 l2月 (582～591) PR0GRESS IN GEOPHYSICS Dec．2002 波动方程地震偏移成像方法的现状与进展 刘喜武 刘 洪 (1．中国科学院地质与地球物理研究所 ，北京 100101； 2．中国科学 院研究生 院，北京 100039) [摘 要 】 综述 了波动方程地震偏移成像方法研 究的现状 ，指 出了各种算法 的 问题和 解决方案 ； 给 出了相空间小波分析步进算法偏移成像 以及保幅保结构李群算法偏移成像 的最新进展 ；对波动 方程偏移成像发展趋势进行展 望 ． 【关键词 】 波动方程地震偏移成像 ；算法 ；相空间；小波分析 ；李群 ；研 究进展 [中图分类号 】 P315 [文献标识码 】 A [文章编号 】 1004-2903(2002)04-0582-10 0 引 言 地震偏移是一种将地震信息进行重排的反演运算 ，以便使地震波能量归位到其空间的 真实位置 ，获取地下真实构造 图像 ¨J．除了深度域构造成像外 ，地震偏移还为其它特殊处理 提供振幅、相位等信息 ，用于速度估计和属性分析 ，建立在波动方程基础上的地震偏移成像 技术代表 了地震处理的极致 J． 地震偏移最初是在水平迭加基础上进行 的，目的是使倾斜界面共深度 “映像 ”聚焦 ，使绕 射波归位 ，即将能量还原到它们正确位置上 ．早期人工偏移是按照偏移空间的时距关系作 图；若将共深度点剖面看作一系列绕射点组成 的源反射 ，可用计算机实现对这些绕射点的偏 移 ，即建立在射线理论基础上的绕射扫描迭加方法 以及后来 的Kirehhoff偏移 ．20世纪 7O年 代初美 国斯坦福大学 以J．F．Claerbout为首的 SEP研究小组第一个对标量波动方程提出了 有限差分近似解法 ．4J，实现 了地震偏移 ．此后建立在波动方程基础上 的地震偏移成像方法 如有限差分法 、Kirehhoff积分法 、F-K方法及其各种变形等方法广泛应用 J．爆炸反射面模 型 J(b~ewenthaleta1．，1976)为波动方程偏移成像条件 的建立奠定 了理论基础 ．由于波动方 程描述地震波地下传播规律，因而波动方程偏移一方面可以解决复杂介质条件下成像 问题 ， 另一方面保持 了波场的动力学特征 ． 地震偏移各种方法最初是作为时间偏移方法 出现 的，目的是满足二维时间叠加剖面成 图需要 ，后来为满足横 向变速情况下成像精度需要 ，发展 了深度偏移方法 ’8(Hubral1977， Earner1981)．近 20年来，偏移方法又发展到了三维和叠前偏移[，三维叠前深度偏移代表地 震偏移 的发展水平 ．当今各种各样 的偏移技术方法极为丰富，如时间偏移 、深度偏移 、二维 、 三维、迭前、迭后；如使用共炮集、共方位角道集 、面炮方法实现等等 ．偏移算法也多种多样 ． 实际应用 中根据具体情况和要求选取相适应 的方法 ．考虑如此复杂而庞大的内容 ，本文只对 波动方程地震偏移中核心算法——波场延拓和成像 的现状与进展进行阐述 ，并对其发展趋 势进行展 望 ． [收稿 日期] 2002-01—18； [修 回日期] 2002-04—04． [基金项 目] 中国科学院知识创新工程项 目(KZCXL-SW-18)资助 ． [作者筒介] 刘喜武，1971年生 ，黑龙江讷河人 ，中国科学院地质与地球物理研究所博士研究生，主要从事地震信号处理 和成像研 究 ． 维普资讯 4期 刘喜武 ，等 ：波动方程地震偏移成像方法 的现状与进展 ·583 。 1 波动方程地震偏移成像方法 的研究现状 波动方程地震偏移成像的各种方法都是建立在波场反 向外推基础上 ，按照算法实现 的 原理可 以分为两大类 ：基于射线理论 的偏移方法和基于波场延拓 的偏移方法 ¨ ．基于射线 的Kirchhoff积分类方法