双平方根单程波动方程叠前τ偏移方法

本文将常规双平方根（DSR）单程波动方程从深度域变换到双程垂直走时（τ）域，由此推导出可从数学上实现“沉降观测”的单程波DSR传播算子. 其递归波场延拓算法包含波数域针对常速背景的相移处理和空间域针对横向速度扰动的相位校正，可以应对上覆地层速度横向变化对构造成像的影响. 结合零炮检距、零时间成像条件，提出了在τ域进行波场延拓与成像的DSR方程叠前偏移新方法. 为了克服其全三维偏移算法在实际应用中可能面临的困难，本文采用稳相近似，在crossline常炮检距偏移理论基础上推导了实用的共方位角叠前τ偏移方法. 数值试验表明，DSR方程叠前τ偏移在强横向非均匀介质中的成像精度与分辨率优于传统的时间域成像技术.     

共炮检距道集波动方程保幅叠前深度偏移方法

本文提出了一种基于双平方根算子的共炮检距道集波动方程保幅叠前深度偏移方法,将振幅误差补偿作为偏移的一部分与“运动学偏移”一起在偏移过程中实现.其基本内容包括：(1)从保幅的单平方根算子方程出发,推导出由双平方根算子定义的保幅单程波方程;(2)根据地震波摄动理论把速度场分裂为层内常速背景和变速扰动,分别在频率－波数域和频率－空间域求得波场深度延拓的偏移时移量及振幅校正系数,从而得到最终的DSR保幅波场延拓算子;(3)在高频假设条件下,把DSR保幅波场延拓公式中的积分运算进行稳相近似,得到保幅波场延拓的相移公式.理论分析和模型数值试验表明,该方法不但可以使散射能量聚焦、归位,提高成像精度;而且可以输出正确反映地下反射系数的振幅信息,为后续的地震属性分析（如AVO/AVA）提供更真实的地震信息.     

振幅保真的单程波方程偏移理论

本报告综述了近年来发展起来的真振幅单程波方程偏移理论，对各种基于单程波方程的偏移振幅作出了系统的分析介绍.为了得到正确的偏移振幅，叠后偏移之前必须进行球面扩散校正.但是，作为叠后相位移偏移算法的推广，Dubrulle氏的共炮检距偏移除了零炮检距和平层两种特例外，它不能给出正确的偏移振幅.通过分析，我们发现传统的单程波单炮偏移不是保振幅算法.利用真振幅单程波方程分解和校正地表的初值条件，我们可以将常规偏移方法改造为真振幅算法并且证明在高频渐近的意义下，新方法的偏移振幅等价于Kirchhoff反演的结果.进一步研究发现可以利用单平方根算子和双平方根算子输出真振幅共反射角剖面.我们的分析指出，这时正确的成像条件为乘积pUp*D.与真振幅共炮偏移所需的商型成像条件pU/pD相比，共反射角偏移的计算稳定性大为改善并且在实际地震资料处理中有重要的应用.     

利用矩阵分解理论的双程波方程叠前深度偏移方法

叠前深度偏移理论及方法一直是地震数据成像中研究的热点问题.业界对单程波叠前深度偏移方法和逆时深度偏移开展了深入的研究,但对双程波方程波场深度延拓理论及成像方法的研究还鲜有报道.本文以地表记录的波场值为基础,利用单程波传播算子估计波场对深度的偏导数,为在深度域求解双程波方程提供充分的边界条件,并提出利用矩阵分解理论实现双程波方程的波场深度外推.通过对强速度变化介质中传播波场的计算,与传统的单程波偏移方法相比,本文提出的偏移方法计算的波场与常规有限差分技术计算的波场相一致,证明了本方法计算的准确性.通过对SEAM模型的成像,在相同的成像参数下,与传统的单程波偏移算法和逆时深度偏移算法方法相比,本文提出的偏移方法能够提供更少的虚假成像和更清晰的成像结果.本文所提偏移算法具有深度偏移和双程波偏移的双重特色,推动和发展了双程波叠前深度偏移的理论和实践.     

窄方位地震数据双平方根方程偏移方法探讨

简单回顾双平方根（DSR）方程偏移一般理论之后,介绍了适应偏移速度分析的角度域成像快速算法.针对地震勘探中窄方位三维地震数据成像问题,探讨了限定输入数据体或限定波场传播方向的DSR方程偏移方法.数值试验表明,共方位角偏移等DSR方程叠前深度偏移方法在盐体或礁体、古潜山等强横向变速地区具有很好的应用前景.     

角度域叠前时间偏移振幅补偿方法

本文从单程波方程出发,推导出角度域叠前时间偏移成像幅值计算方法.文章首先在时间域实现了速度梯度引起的幅值透射损失补偿算法;借鉴波动方程偏移的反褶积成像条件,对地震波能量随传播距离的几何扩散进行部分补偿;针对观测系统中炮点覆盖引起的角度域非均匀照明提出了相应的振幅补偿方法.三维水平层状理论模型成像结果表明获得的角道集具有较好的AVA特性.     

波动方程叠前深度偏移直接产生角道集

基于单程波深度延拓算法,得到震源波场在成像点的入射角度,结合对地层倾角的估计,获得入射地震波与界面法线的夹角.通过运用"保幅"的反褶积成像条件和考虑累加的炮数,解决了炮点覆盖不均匀导致的成像幅值误差问题,进而建议了炮域波动方程叠前深度偏移直接产生角道集的方法和流程.与基于空间移动或时移成像条件的波动方程叠前深度偏移提取角道集的方法相比,本文建议的方法只需少量额外的存储空间,又可补偿观测系统非均匀覆盖对成像幅值的影响;其增加的计算量与炮域偏移算法相比几乎可以忽略.文中算例表明,本文方法提取的角道集可为叠前反演提供较精确的AVO振幅特性.此外,就改善地震成像效果本身而言,提取角道集使得可在波动方程叠前深度偏移中应用剩余动校和拉伸切除技术,从而可更好地保持高频成分并提高成像的信噪比.     

频率-空间域有限差分法叠前深度偏移

为了处理横向强变速介质中的深度成像问题,本文提出一种基于共炮道集的优化系数的傍轴近似方程叠前深度偏移算子,并在基于反射系数估算的成像条件下,可实现叠前深度偏移成像.该算子具有方程阶数低且能对陡倾角成像的特征,并采用有限差分法波场延拓,能适应速度场的任意变化.当在频率-空间域进行计算时,相对于纯粹的时间-空间域有限差分算法有计算效率高、成像方便的优点.脉冲响应测试和对Marmousi模型进行的叠前深度偏移结果表明,该偏移方法在强横向变速情况下具有非常好的成像效果.     

基于Gabor-Daubechies小波束叠前深度偏移的角度域共成像道集

传统炮检距域共像集（CIG）在复杂介质中因波传播的多路径而存在反射体位置不确定的问题. 角度域CIG由于克服了这一缺陷而逐步成为速度分析、AVA以及振幅保真偏移成像等研究的主要手段. 以波动理论为基础的地震偏移成像方法的发展为获得高质量的角度域CIG提供了可靠的实现途径. 其中,基于波场局域化分解和传播的小波束域波场延拓和偏移成像方法,因其波场分解基本函数和传播算子在空间和方向上的双重局域特性,而成为角度相关分析研究的有效工具. 本文在采用Gabor Daubechies框架分解的小波束叠前角度域偏移成像基础上,利用不同的叠加方法由局部角度域像矩阵得到了反射角域CIG（CRAIG）和倾角域CIG（CDAIG）. 以SEG EAGE二维盐体模型为例,通过对CRAIG和CDAIG的对比,探讨了这两种角度域CIG的特点及其在地震偏移成像中的潜在应用.     

起伏地表条件下分块法插值叠前深度偏移

为实现复杂起伏地表条件下基于单程波动方程的高效叠前深度偏移,本文采用"逐步一累加"思想,应用傅里叶有限差分法(FFD)算法实现复杂起伏地表务件下的波动方程叠前深度偏移.为提高效率,本文提出在整个成像实现过程中,将整个成像区域分两块区别对待,在地表高程最低点所在的水平面到基准面的区域内,应用常规精细延拓步长进行波场延拓及...     

基于波动方程的广义屏叠前深度偏移

地震波传播算子的计算效率和精度是制约三维叠前深度偏移的关键因素. 广义屏传播算子(GSP, Generalized Screen Propagator)是一种在双域中实现的广角单程波传播算子. 这一方法略去了在非均匀体之间发生的交混回响,但它可以正确处理包括聚焦、衍射、折射和干涉在内的各种多次前向散射现象. 通过背景速度下的相移和扰动速度下的陡倾角校正,广义屏算子能够适应地层速度的强烈横向变化. 这种算子可以直接应用于炮集叠前偏移,通过将广义屏算子作用于双平方根方程,还可以获得一种高效率、高精度的炮检距域叠前深度偏移方法,用于二维共炮检距道集和三维共方位角道集的深度域成像. 本文首先简述了炮检距域广义屏传播算子的理论,进而讨论了共照射角成像(CAI, Common Angle Imaging)条件,由此给出各个不同照射角(炮检距射线参数)下的成像结果,进而得到共照射角像集. 由于照射角和炮检距的对应关系,共照射角像集又为偏移速度分析和AVO(振幅随炮检距变化)分析等提供了有力工具.     

基于共聚焦点道集的叠前深度偏移

共聚焦点(CFP)偏移技术是一种基于等时原理,将Kirchhoff积分法的一步偏移分两步聚焦（即激发聚焦和检波聚焦）来完成的叠前地震成像方法.该方法借助于逆时聚集算子和共聚焦点道集来实现叠前偏移成像.基于共聚焦点道集的叠前深度偏移是把基于共炮集的深度偏移的算法引入到CFP技术上来,基于波场延拓的理论来实现偏移成像,该方法首先生成共聚焦点道集,然后基于面炮合成的理论合成聚焦震源,最终通过相关成像来实现叠前偏移成像.该方法选取较少的聚焦点就可以实现对于地下构造的偏移成像,和炮域波动方程偏移相比,其计算效率得到了提高.通过模型试算和实际资料的试处理,验证了该方法在实现叠前深度偏移成像上的正确性和有效性.     

三维叠前偏移新技术

首次采用空间（ｘ，ｙ）坐标旋转的方法推导出了适应于任何炮检方向并与炮检方向相关的三维波动方程，理论上定量地说明了炮检方向及炮检距与偏移响应的关系，进而证明３Ｄ－ＤＭＯ（叠前部分偏移）只需在炮检方向上进行二维ＤＭＯ即可得以实现。理论模型试算结果表明，该方法正确描述了波的实际传播过程，实现了三维叠前全偏移及３Ｄ－ＤＭＯ处理。     

两种类型共炮点域相关型叠前深度偏移方法分析

研究不同偏移方法的成像机理是实现复杂构造高精度成像的前提,研究了两类共炮点域的相关型叠前深度偏移成像方法：基于单程波动方程的叠前深度成像方法和基于双程波动方程的叠前深度成像方法,同时对比了它们在计算效率、数据存储量、成像精度、成像机理、速度敏感性等方面的差异及其共性。以复杂构造模型为例,采用了傅里叶有限差分法（FFD）和逆时成像法（RTM）,这两种方法实现了121个共炮点道集的叠前深度偏移成像处理。计算结果表明,当速度准确时,两种深度域波动方程成像方法均可以恢复出各个地质反射界面,其中逆时偏移对陡倾角成像效果显著,当速度存在百分比误差和随机扰动情况时,逆时成像结果要差于单程波方法,因此,逆时偏移方法对速度的敏感性较大,且低频噪声较为严重。     

基于局部余弦基小波束的观测系统沉降法叠前深度偏移

将局部余弦基小波束波场分解、传播与观测系统沉降法叠前深度偏移相结合,推导了源-检波器观测系统沉降法传播算子.本算法中,先对频率域的共点源和共点检波器道集做局部余弦小波束分解,然后分别沿共小波束源和共小波束检波器在深度方向延拓得到下一层波场.每个深度的波场,都等效于把源和检波器放在该层后所能接收到的地震记录,每点的像值由炮点和检波点重合时的零时刻波场值给出.通过二维SEG/EAGE盐丘模型和Marmousi模型的偏移成像结果验证该方法理论推导的正确性.另外,结果显示该方法继承了小波束域波场延拓在速度扰动较大情况下波传播及成像精度高的优点.     

双复杂介质条件下频率空间域有限差分法保幅偏移

油气勘探的重点正转向复杂地表条件和复杂地质条件的区域.双复杂条件下的叠前深度偏移是解决复杂地表条件和复杂地质构造成像的有效手段.基于“逐步累加”的“直接下延”法是解决复杂地表成像的有效手段,能够较好地消除地形起伏的影响.波动方程频率空间域有限差分(xwfd)叠前深度偏移对介质速度横向变化有较强的适应性,适宜于复杂构造的偏移成像,同其他常规波动方程深度偏移一样,常规的xwfd偏移方法,主要也是针对相位进行波场延拓,没有对振幅做任何处理.我们基于保幅单程波方程,推导出了基于xwfd的保幅波场延拓算子,针对xwfd求解时引入误差的影响,我们在xwfd保幅波场延拓过程中加入了误差补偿,实现了带误差补偿的xwfd保幅偏移.基于带误差补偿的xwfd保幅算子,应用适合起伏地表的直接下延法,对双复杂介质模型和实际资料进行了试算,改善了双复杂介质的成像效果.其中,误差补偿可以在若干个外推步长上进行,所以相对于保幅傅里叶有限差分(ffd)法偏移来说,该方法在改善成像质量的同时,也具有较高的运算效率.     

起伏地表条件下基于复Pade逼近的叠前深度偏移

叠前深度偏移是解决复杂地表和复杂构造地震成像的有效技术,而波场“直接下延”法实现了复杂地表条件下的地震成像.基于上述成果,结合高精度的波场延拓算子,本文提出了一种新的叠前深度偏移方法,这种方法是在波场延拓时,对声波方程中的平方根项进行复Pade逼近,通过推导得到基于复Pade逼近的傅里叶有限差分算子,结合波场“直接下延”法,实现了起伏地表条件下的叠前深度偏移,该算法减少了偏移噪音,从而得到准确、稳定的偏移成像结果.通过理论模型试算和实际资料试处理,验证了该方法的有效性.     

基于单程波方程的角度域保幅偏移(英文)

传统叠前深度偏移只能够提供地下的构造信息,但工业界在需要构造信息的同时还要与地下界面反射系数成比例的振幅信息。最近几年,基于单程波方程的保幅叠前深度偏移算法有了一定的发展,但是,基于炮域、单程波的保幅型叠前深度偏移必须应用反褶积型的成像条件,这种成像条件在构造复杂、速度变化剧烈的地区会出现不稳定现象。基于角度域的保幅深度偏移克服了这一不稳定性缺点的同时,还域的保幅深度偏移,模型和实际资料的试算分析验证该思路方法的正确性和有效性。     

粘弹性介质单程波法非零偏移距地震数值模拟与偏移

地震资料分辨率降低,得不到深层介质的精确信息实际上是由于大地吸收效应的影响.同时与双程波动方程相比单程波动方程避免了多次波的干扰并且计算效率高、占用内存少.本文首先基于开尔芬粘弹性介质模型将品质因子与单程波分步傅立叶法波场延拓算子相结合,实现了粘弹性介质波场延拓,从而将单程波弹性介质波场延拓推广到了粘弹性介质.然后在定位原理,数学检波器原理以及等时叠加原理的基础之上实现了粘弹性介质非零偏移距叠前正演模拟.最后将数值模拟得到的正演记录进行弹性偏移和粘弹性偏移并进行对比分析.通过数值算例可以看出,粘弹性介质叠前正演深层的反射波能量减弱,同相轴变粗,频带变窄,主频减小,分辨率降低;粘弹性偏移不但实现了振幅的恢复,而且同时偏移剖面的垂向空间分辨率也得到了提高.     

角度域叠前时间偏移:非均匀介质中的单程波方法

从单程波方程出发,推导出角度域叠前时间偏移的走时、入射波与反射波夹角、成像幅值计算方法;构建了可直接生成角度域成像道集的叠前时间偏移方法与偏移流程.文中定量分析了速度梯度对走时、角度、幅值的影响,给出了可更好考虑介质非均匀性的角度、幅值计算方法.理论模型以及实际数据验证了本文方法的有效性.