资源描述：
《波动方程叠前深度偏移方法综述》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、第30卷第3期勘探地球物理进展Vol.30,No.32007年6月ProgressinExplorationGeophysicsJun.,2007文章编号:16718585(2007)03015309波动方程叠前深度偏移方法综述马淑芳,李振春(中国石油大学(华东)地球资源与信息学院,山东东营257061)摘要:对波动方程叠前深度偏移的基本概念、实现方法进行了分类和阐述。波动方程叠前深度偏移方法主要分为2类:一类是单平方根方程偏移,在偏移过程中,上、下行波分别向下延拓,并通过互相关成像条件来提取成像值;另一类是基于沉降观测概念的双平方根偏移
2、,在偏移过程中,炮点和检波点同时向下外推,当两者重合时(零偏移距),零时间的波场值就作为该空间点的成像值。对共炮真振幅偏移进行了阐述,并指出也可以在角度域道集实现该算法。理论模型的处理效果证明,波动方程叠前深度偏移成像技术是解决强横向变速情况下复杂构造成像的一种有效手段。关键词:波动方程;叠前深度偏移;真振幅偏移中图分类号:P631.443文献标识码:A!!自20世纪70年代Claerbout提出著名的15∀测概念的双平方根(DSR)偏移。在其偏移过程偏移方程以来,偏移成像研究的重点转移到以波动中,只对上行波场进行延拓,相当于向地下延拓(沉方程为
3、理论基础的研究。波动方程叠前深度偏移降)震源与接收点,当二者重合时(零偏移距),零时成像技术是解决强横向变速情况下复杂地质体成间的波场值就作为该空间点的成像值。像的有效方法。基于计算方法的不同,可分为有限同时真振幅波动方程偏移能实现对球面扩散差分偏移法和Fourier偏移法。前者借助于差分损失的补偿,恢复期望的反射系数,还能在进行构计算,把速度、密度等介质参数的影响体现在差分造成像的同时给出震源子波的信息,所以真振幅偏计算的矩阵方程中,这类方法能自动适应速度场的移对于AVO、AVA和偏移后的振幅解释具有特[12,13]任意变化。后者一般借助快速F
4、ourier变换来进别重要的意义。21世纪初,Zhang等在真振行波场延拓计算。最早期的Fourier偏移方法应幅波动方程偏移方面的研究取得了较大的进展,该当是相移偏移,该算法精度高,无频散现象,但要求方法从改进型单程波方程入手,通过反假频和成像速度不能有横向变化,因而在实际生产中的应用受条件的改善实现真正意义上的保幅偏移成像。到了限制。后来提出了适应速度横向缓慢变化的[1]各种修正方法,如相移加插值法(PSPI)和扩展1!单平方根方程波场外推理论方法[2]Stolt方法等。近年来,结合有限差分法和Fourier法各自的优点,发展了各种各样的混
5、合法,如单平方根算子炮域叠前深度偏移是将炮点和[3][4]分步Fourier法(SSF)、伪屏法(GSP)及Fou检波点分别向下外推。检波点外推输入数据为地[5]rier有限差分法(FFD)等。为了提高对横向局表接收的共炮点道集,采用上行波方程外推地下波[6]部变速的适应能力,Hill和BrandsbergDahl场。对于炮点外推,需要人为给定一个地震子波作[7]等提出了高频方法,同时也提高了计算效率。为震源函数,采用下行波方程外推地下波场,得到[8][9]Chen等提出的小束波偏移和Margrave等提单炮传播波场。在同一物理点上,根据成像
6、条件二出的非稳相相移偏移等新方法也在研究中。李录者互相关叠加求和,实现炮集单平方根算子的叠前[10]明等进一步研究了转换波叠前深度偏移方法。深度偏移成像。[11]左建军等利用Gazdag相移法实现了井间地震波动方程叠前深度偏移的核心问题是波场外成像,获取了高分辨率的成像剖面。推算子。一般要求算子能适应陡倾角反射的成像典型的波动方程偏移方法有2类。一类是单及较剧烈的横向速度变化,同时还要具有较高的计平方根方程偏移,如炮集偏移和合成面炮偏移。在算效率。波场外推算子的研究伴随着偏移成像理其偏移过程中,上、下行波基于各自的单程波方程收稿日期:20060
7、710;改回日期:20070310。分别进行延拓,并通过2个延拓波场的互相关(零第一作者简介:马淑芳(1984),女,在读硕士,主要从事地震波传时间条件)来提取成像值。另一类是基于沉降观播与地震成像理论与方法的研究工作。154勘探地球物理进展第30卷论发展的全过程。以下介绍几种常见的单平方根1.1.3!Fourier有限差分偏移算子波动方程叠前深度偏移算子。基于对相移处理和差分处理方法的优、缺点的[5]1.1!波动方程波场外推算子认识和理解,Ristow等提出的Fourier有限差分1.1.1!频率空间域有限差分偏移算子法偏移融合了上述
8、2种方法的优点。该方法是在早期的有限差分偏移方法主要在时间空间域分步Fourier算子的基础上,增加了关于二阶扰动[14,15]上实现
