声波和弹性波叠后逆时深度偏移

对比声波和弹性波叠后逆时深度偏移的成像效果,从计算波动理论出发,采用凹陷模型制作了基于爆炸反射界面原理的零偏移距剖面。并根据零时间逆时成像条件,制作了较高精度的声波、弹性波叠后逆时成像剖面。逆时延拓数值试算表明,采用两种方法均能恢复出准确的地质界面,其中声波在倾斜界面情况下与实际地质模型有所偏移,没有准确实现地震波场归位;而弹性波垂直分量得到了准确成像,倾斜界面清晰,而且垂直分量直接显示出断点和倾斜界面位置,弥补了声波深度成像的不足。图2参8     

基于高阶有限差分的波动方程叠前逆时偏移方法

高阶有限差分波动方程叠前逆时偏移方法的基本原理是,从二维1阶双曲型地震波波动方程出发,通过时间上的2阶差分和空间上的任意高阶交错网格差分对该方程进行离散,得到高阶有限差分逆时偏移算子。采用单频双程波动方程计算初至走时,并以此作为叠前逆时深度偏移的成像条件。设计了凹陷模型和断陷模型,对方法进行了验证。模拟结果表明,基于高阶有限差分的波动方程叠前逆时深度偏移方法可以对复杂地层进行准确成像,并能消除由逆时偏移引起的低频噪声。     

各向异性弹性波叠后逆时深度偏移

大量数据表明地下岩层的地震各向异性是普遍存在的。然而在常规偏移中，通常视地下介质为各向同性介质，这必然导致一定的误差。研究了各向异性弹性波的有限元逆时延拓算法及有关人工吸收边界条件等问题，实现了各向异性弹性波有限元逆时深度偏移系统。2个理论模型的自激自收资料实算表明，该系统能够清晰准确地进行偏移成像。另外，对这两个模型资料的各向同性移证实，常规各向同性偏移剖面上存在较大误差。     

波方程傅氏有限差分法叠前深度偏移

通过对几种叠前深度偏移算法进行对比，分析了各自的优缺点。着重阐述了近年来发展的波动方程傅氏有限差分算法的原理及实现过程，该方法将波场延拓算子分解成频率波数域和频率空间域的三个算子进行运算。结合了相移法和有限差分算法的优点，克服了两种算法的不足。与目前广泛应用的波动方程Kirchhoff积分法叠前深度偏移相比，具有成像精度高，保持地震波动力学特征等优点。应用FFD偏移成像原理，研制了波动方程傅氏有限差分法叠前深度偏移软件，在Marmousi模型上成功地进行了FFD叠前深度偏移处理，取得了理想的成像效果。     

时限时移相关法叠前逆时成像条件及其应用

地震波逆时偏移法是目前成像精度最高的成像技术。分析了相关法叠前逆时成像的基本原理,提出了时限时移相关法叠前逆时成像条件。运用高阶交错网格有限差分法,按照一定的观测方式合成了基于Marmousi模型的123个炮集记录,并应用时限时移相关法叠前逆时成像条件进行了叠前逆时深度偏移,将获得的剖面与反射系数深度剖面进行了对比。数值计算结果表明,应用时限时移相关法叠前逆时成像条件获得的偏移剖面准确反映了Marmousi模型中的各个地质反射界面及其反射系数的变化特征。与时空移相关法叠前逆时成像条件相比,时限时移相关法叠前逆时成像条件能够大大减少计算量和存储量,提高了计算效率。     

两种声波方程的叠前逆时深度偏移及比较

 本文分别从双程声波方程和任意广角声波方程出发，推导了地震波逆时传播的高阶有限差分格式，给出了逆时偏移的稳定性条件、边界条件与成像条件的计算方法，实现了这两种声波方程在二维各向同性介质中的逆时偏移。上述两种方程对理论数据的偏移结果表明：双程声波方程的逆时偏移无倾角限制，但层间反射会严重影响偏移剖面质量，而任意广角声波方程则可消除层间反射的影响，得到较好的成像效果。     

波动方程叠前深度偏移方法综述

:对波动方程叠前深度偏移的基本概念、实现方法进行了分类和阐述。波动方程叠前深度偏移方法主要分为2类：一类是单平方根方程偏移,在偏移过程中,上、下行波分别向下延拓,并通过互相关成像条件来提取成像值；另一类是基于“沉降观测”概念的双平方根偏移,在偏移过程中,炮点和检波点同时向下外推,当两者重合时（零偏移距）,零时间的波场值就作为该空间点的成像值。对共炮真振幅偏移进行了阐述,并指出也可以在角度域道集实现该算法。理论模型的处理效果证明，波动方程叠前深度偏移成像技术是解决强横向变速情况下复杂构造成像的一种有效手段。     

双相各向异性介质弹性波传播交错网格高阶有限差分法模拟

含流体孔隙各向异性介质中的弹性波传播问题是目前石油勘探和地震学研究的难点之一。基于Biot理论,本文提出了二维双相任意倾斜各向异性介质三分量弹性波方程交错网格任意偶阶精度有限差分解法,并对均匀及两层双相VTI介质和TTI介质中的弹性波场进行了模拟。波场快照和合成记录表明:当横波源的偏振方向与TI介质的对称轴存在一定夹角、且耗散系数较小时,在固相与流相中存在快横波、慢横波、快纵波和慢纵波;双相TI介质中横波分裂现象及波前面尖角和三分叉现象与TI介质中的基本相同,只是黏滞相界中弹性波速度有所减小,且具有频散和衰减特性;从直达波转换为反射波后,该波型具有反射波的性质。     

标量声波波动方程高阶交错网格有限差分法

为进一步提高标量声波波动方程正演数值模拟精度,减小数值频散和计算量,改善边界吸收效果,提出将高阶交错网格有限差分法直接用于求解标量声波波动方程的方法,并推导出了时间和空间导数的高阶交错网格有限差分格式,给出了数值模拟计算所需的稳定性条件和吸收边界条件.模型试算结果表明,与传统的中心有限差分法相比,本文方法不仅可以快速提高时空导数的差分精度,而且具有频散小、计算效率高等特点,可推广应用于三维标量声波高精度正演数值模拟和叠后逆时深度偏移.     

双程声波方程叠前逆时深度偏移的成像条件

 成像条件是决定波动方程叠前逆时深度偏移效果的重要因素之一。本文从构造成像角度研究了不同的成像条件计算方法在双程声波方程叠前逆时深度偏移中的应用效果。模型试验表明,不论采用何种成像条件,基于双程波方程的叠前逆时深度偏移均会产生较强的层间反射,该反射会恶化浅层成像结果。应用波场归一化互相关成像条件时,层间反射的影响相对较弱;基于波场互相关的逆时偏移成像条件对深部构造的成像能力总体上优于求解程函方程得到的成像条件;而归一化互相关成像条件能够在对地震波进行偏移成像的同时补偿深层能量,在相同条件下,它对深部地层的成像能力更强。     

弹性波Kirchhoff叠前深度偏移速度分析

本文提出的弹性波偏移速度分析方法是基于Kirchhoff叠前深度偏移形成纵波和转换波炮检距域共成像点道集拉平准则,分别对纵波和横波偏移速度进行更新。当两分量地震数据成像深度不-致时,通过调整横波偏移速度进行深度匹配,完成高精度的弹性波场偏移速度分析。文中分别给出速度更新及深度匹配方法。模型数据和实际资料试算结果表明了该方法的可行性和有效性。     

任意起伏地表弹性波方程交错网格高阶有限差分法数值模拟

任意起伏地表弹性波数值模拟主要涉及两个问题，其一是如何求解弹性波方程；其二是如何处理自由边界条件。本文首先从泰勒级数展开式出发，推导出交错网格一阶空间导数的任意偶数阶精度展开式和相应差分系数计算式以及一阶双曲型应力一速度弹性波方程交错网格任意偶数阶精度差分格式求解方程；然后采用将零速度法和广义虚像法相结合的方法来处理自由边界，并在自由边界上采用四阶精度差分格式；运用上述方法对光滑起伏模型和任意起伏模型进行了数值模拟试验。结果表明，本文所述方法稳定性好、模拟精度高，且适合任意起伏地表弹性波波场模拟。     

二维Offset平面波有限差分法叠前时间偏移

在岩性油气藏勘探中,保振幅或保波形成像是进行储层特征分析和参数估计的基础.由于Kirchhoff积分叠前时间偏移、F-K域波动方程叠前时间偏移不适应弱横向变速介质情形,因此,基于偏移距(Offset)域平面波偏移思想,提出了Offset平面波有限差分法叠前时间偏移方法.讨论了CMP道集平面波分解原理,推导了Offset平面波方程及其有限差分解法,给出了角度道集生成方法.利用大庆油田的一个横向缓变的速度模型,对Offset平面波有限差分法叠前时间偏移方法的有效性进行了验证.结果表明,该方法适应缓横向变速介质情形,波场外推算子保幅,计算效率高,可以为AVP/AVA分析提供较好的基础数据.     

地震波旅行时计算方法及其模型试验分析

地震波旅行时的计算是地震波偏移成像技术和地震层析成像技术的核心内容,同时也是实现保幅叠前偏移成像的重要因素。推导了地震波波动方程高阶有限差分法计算格式,提出了最大绝对振幅能量节点初至旅行时拾取的方法。均匀介质模型试验表明,该方法与用解析解、程函方程法计算的旅行时等值线基本重合。基于层状介质模型、凹陷模型和Marmousi模型对方法的适用性进行了试验,将由该方法计算得到的初至旅行时作为叠前逆时成像条件,对凹陷模型和Marmousi模型的数值模拟记录进行了叠前逆时偏移成像,偏移叠加剖面很好地反映了模型的地质特征。     

碳酸盐岩区单程波方程保幅叠前深度偏移

本文首先细致分析了保幅偏移方程各项的物理含义;然后介绍了理论反射系数的计算以及近似表达式;通过简单模型的计算,对比了成像结果与理论反射系数间的关系;最后将单程波保幅偏移技术应用于中国南方A探区碳酸盐岩地震资料的处理,实际处理结果表明本文方法对碳酸盐岩区地震资料的深层能量进行了有效补偿,中、深层信息更加丰富,可为AVO异常测定提供高质量的CIP道集。     

三维各向同性介质弹性波方程交错网格高阶有限差分法模拟

给出了三维各向同性介质中一阶应力一速度弹性波方程交错网格任意偶数阶精度有限差分格式，推导出了三维各向同性弹性介质完全匹配层吸收边界条件公式和相应的交错网格高阶有限差分格式。对三维French模型进行了弹性波模拟，结果表明，该方法模拟精度高，边界吸收效果好。在三维French模型的XOZ平面和yoz平面的弹性波场快照中可以见到断面反射波、侧面反射波和散射波等波场特征。多次侧面反射波和多次散射波说明不同三维构造所形成的波场是相互影响的。     

三维混合延拓一步法波动方程深度偏移

对于二维陡倾斜地层，相位移延拓方法成像精度高、稳定性好、速度快，但难以适应速度的横向变化；频率、空间域的有限差分算法简单、稳定性好、能适应速度的纵、横向变化，但成像精度低、运算速度慢；而用快速的45°有限差分法与相位移延拓法相结合，则能使整个延拓过程既适应速度的纵、横向变化，又能够使得陡倾角地层精确归位。本文把这种混合延拓思想拓展成三维相位移加有限差分混合延拓法，并进一步提出了能适应速度纵、横向变化，构造任意复杂介质的三维一步法深度偏移。数值模拟试验结果表明，该方法具有移速度快、精度高的特点。