叠前逆时深度偏移低频噪声压制策略研究

低波数噪音严重影响了逆时偏移的成像质量,选择合适的低波数噪音压制方法是进行逆时偏移时必然涉及到的一个重要方面.本文首先详细分析了逆时偏移中低波数噪音产生的本质原因,然后根据噪音的产生机制对目前各类低波数噪音压制方法进行了分析和归类:第一类低波数噪音压制策略的核心是在构建波场中避免发生背向反射;第二类低波数噪音压制策略的核心是进行选择性成像.最后通过对压噪方法进行合理的选择和搭配,构建了在不同实现条件和成像需求情况下逆时偏移去噪方法的最佳抉择策略.结合数值算例对几种压噪方案进行了测试和比较,实例表明了对噪音压制方法进行合理的搭配可以获得高信噪比逆时偏移成像结果.     

面炮叠前深度偏移与控制照明技术

基于波场延拓的叠前深度偏移是实现复杂构造地质体成像的最可靠方法,但存在着计算量大、对观测系统适应性差等缺点。面炮偏移是波动方程实现精确叠前成像的另一类方法,具有较高的计算效率,不存在偏移孔径问题,而且可以通过控制照明方法,解决平面波在目标区域的能量补偿问题。本文采用面炮成像技术进行叠前深度偏移成像,通过对面炮震源下行波场的质量控制以及射线参数的个数和范围的选取,以达到最佳的成像效果。采用不同深度点上的控制照明技术,较大地提高了目标地层的成像精度。数据实验表明面炮成像技术是一种快速有效的方法,其成像精度与单平方根算子的共炮点道集偏移和双平方根算子的共中心点道集偏移相当,但在计算速度上要快得多,而且易于并行计算。     

高精度混合法叠前深度偏移及其并行实现

叠前深度偏移是复杂构造成像的重要手段.本文基于波场分裂理论,首先给出了波场延拓的一般方程,即一个上下行波的耦合方程组.通常所用的波场延拓方程就是该耦合方程组的特例.根据平方根算子的近似,推导了一种新的高精度混合法偏移方法,用分裂法即可进行计算.通过对一个横向强烈变速模型的叠后偏移及Marmousi复杂构造模型的叠前偏移计算,说明了方法的有效性,精度较高.采用MPI并行编程实现并行计算,提高了计算效率.该方法可用于对横向强烈变速的复杂构造的精确成像.     

高精度混合法叠前深度偏移及其并行实现

叠前深度偏移是复杂构造成像的重要手段.本文基于波场分裂理论,首先给出了波场延拓的一般方程,即一个上下行波的耦合方程组.通常所用的波场延拓方程就是该耦合方程组的特例.根据平方根算子的近似,推导了一种新的高精度混合法偏移方法,用分裂法即可进行计算.通过对一个横向强烈变速模型的叠后偏移及Marmousi复杂构造模型的叠前偏移计算,说明了方法的有效性,精度较高.采用MPI并行编程实现并行计算,提高了计算效率.该方法可用于对横向强烈变速的复杂构造的精确成像.     

双复杂条件下的波动方程叠前深度偏移

复杂地表条件下的地震勘探越来越被人们所关注.双复杂条件下的叠前深度偏移方法是解决复杂地表条件和复杂地质构造成像的有效手段."波场上延"法能实现由非水平观测界面开始的偏移过程,解决复杂地表对地下构造成像的影响.复杂理论模型的试算以及实际资料处理表明,"波场上延"方法较好地克服了起伏地形对地下构造成像的影响,取得了令人满意的效果,实现了波动方程基准面校正和深度成像的有机结合.     

层状各向异性介质转换波克希霍夫叠前时间偏移

在克希霍夫叠前时间偏移处理中,地震波走时的计算方法是决定大偏移距地震资料成像品质的重要因素.在常规的三维转换波各向异性叠前时间偏移公式中,走时的计算是基于等效单层各向异性介质的非双曲线方法.用这种方法处理的成像道集,在偏移/深度比超过一定阈值后,成像道集中的反射同相轴将出现过偏现象,这种偏移不平的同相轴将影响偏移叠加的最佳响应,使得偏移成像波组呈低频化特征,最终降低三维转换波偏移成像质量.我们采用层状介质的走时计算方法代替常规算法,并且利用了常规方法的转换波各向异性偏移速度模型.基于层状介质的算法能够提高大偏移距转换波走时计算精度,克服中浅地层大偏移距远道成像道集中反射同相轴逐渐上翘的问题.两个地区的三维转换波资料处理结果证实,基于层状各向异性介质的转换波克希霍夫叠前时间偏移方法,明显改善了反射成像剖面的连续性和分辨率,提高成像剖面构造的可解释性.     

角度域弹性波Kirchhoff叠前深度偏移速度分析方法

为提高地震成像结果的准确性并真实反映实际地震波场在介质中的传播特性,应该充分利用多分量地震数据的矢量特征进行弹性波成像,其中,最为棘手的问题是纵横波偏移速度场的确定,为此,本文提出了直接利用多分量地震数据进行弹性波角度域偏移速度分析的方法.基于空移成像条件的弹性波Kirchhoff偏移方程提取了弹性波局部偏移距域共成像...     

基于拟线性Born近似的叠前深度偏移方法

在频率波数域和频率空间域实现了一种基于拟线性Born近似的叠前深度偏移方法,并在2-D空间进行了Marmousi模型炮集数据的处理.通过与Split-StepFourier、Phase-Screen和稳定的Born近似叠前深度偏移等方法比较,认为基于拟线性Born近似的叠前深度偏移方法不仅在效果上要优于前三者,而且还能更好地处理速度横向变化.在散射波场计算中,使用了一个更稳定的散射波场计算公式,扩大了拟线性Born近似的应用范围,使基于拟线性Born近似的叠前深度偏移方法能够适应更强的横向速度变化.     

地震照明叠前深度偏移方法综述

对地震照明叠前深度偏移的基本概念、实现方法进行了分类和阐述.地震照明叠前深度偏移是通过对震源和炮记录进行合理的选择和合成,从而进行地震照明成像的一种有效方法.其可以分为平面波偏移和小束波偏移.理论模型的处理效果证明地震照明叠前深度偏移成像技术有很高的计算效率,并且还可以提高地下特定目标的成像质量.     

地震叠前深度偏移方法流程及应用

针对复杂介质的成像问题,提出了一套地震资料叠前深度偏移方法流程,主要包括三部分：（１）地震资料精细预处理;（２）速度－深度模型建立;（３）叠前深度偏移成像．以Ｋｉｒｃｈｈｏｆｆ偏移理论为基础,强调地质与地球物理的综合以及地震处理与解释的一体化．在ＺＸ地区成功地实现了二维地震资料叠前深度偏移,所获得的ＮＥ２０６叠前深度偏移剖面揭示了复杂的ＺＸ古潜山及其内部构造,并清楚地展示了逆掩断层的存在．     

复杂介质下保真振幅Kirchhoff深度偏移

本文详细讨论地下复杂介质地震成像问题.当速度模型含有强横向变化时,射线场会折叠形成多次走时,在这种情况下,文中证明在保真振幅叠前深度偏移时,应当考虑所有的到时才能得到定量的像.通过实例证明,在多次走时情况下,常规的共炮集、共检波器集都存在强的假像;由此提出共角度偏移的概念.最后,在多次走时下,讨论多次偏移算子(走时、相移、振幅等)的高效重构,给出了三维叠前深度偏移的快速算法.     

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

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

Prestack reverse-time depth migration of arbitrarily wide-angle wave equations

Based on arbitrarily wide-angle wave equations,a reverse-time propagation scheme is developed by substituting the partial derivatives of depth and time with central differences. The partial derivative of horizontal direction is replaced with high order difference. The imaging condition is computed by solving the eikonal equations. On the basis of above techniques,a prestack reverse-time depth migration algorithm is developed. The processing exam-ples of synthetic data show that the method can remove unwanted internal reflections and decrease the migration noise. The method also has the advantage of fidelity and is applicable of dip angle reflector imaging.     

时-空局域化地震波传播方法:Dreamlet叠前深度偏移

提出了一种在时间和空间上完全局域化的波场分解和传播算法─dreamlet偏移方法.Dreamlet是一种脉冲-小波束形式的波场分解原子,它利用多维局部分解变换,把时空域波场映射到局部时间-频率-空间-波数相空间,并用局部相空间的传播算子(dreamlet算子)沿深度延拓.本文利用多维局部余弦变换实现dreamlet算法,分解后的波场系数和传播算子不仅有很好的稀疏性,且均为实数,也即波的传播和成像过程完全在实数域实现.文中推导了局部余弦基dreamlet波场分解和传播算子理论公式并将其应用于叠前深度偏移.在dreamlet相空间波的传播过程为稀疏矩阵相乘,而且延拓后的地表数据波场的有效时间长度随深度的增加不断减小,从而可以减少需要传播的波场系数.二维SEG/EAGE盐丘和SIGSBEE模型算例验证了理论推导的正确性,成像结果显示该方法在横向速度变化剧烈情况下有很好的精度.     

时-空局域化地震波传播方法:Dreamlet叠前深度偏移

提出了一种在时间和空间上完全局域化的波场分解和传播算法─dreamlet偏移方法.Dreamlet是一种脉冲-小波束形式的波场分解原子,它利用多维局部分解变换,把时空域波场映射到局部时间-频率-空间-波数相空间,并用局部相空间的传播算子(dreamlet算子)沿深度延拓.本文利用多维局部余弦变换实现dreamlet算法,分解后的波场系数和传播算子不仅有很好的稀疏性,且均为实数,也即波的传播和成像过程完全在实数域实现.文中推导了局部余弦基dreamlet波场分解和传播算子理论公式并将其应用于叠前深度偏移.在dreamlet相空间波的传播过程为稀疏矩阵相乘,而且延拓后的地表数据波场的有效时间长度随深度的增加不断减小,从而可以减少需要传播的波场系数.二维SEG/EAGE盐丘和SIGSBEE模型算例验证了理论推导的正确性,成像结果显示该方法在横向速度变化剧烈情况下有很好的精度.     

VTI介质中准P波方程叠前逆时深度偏移

根据具有垂直对称轴的横向各向同性（VTI）介质中的一阶准P波方程,导出了该方程在交错网格中逆时延拓的高阶有限差分格式,给出了其稳定性条件,采用完全匹配层吸收边界条件解决边界反射问题,分别应用下行波最大能量法和归一化互相关成像条件, 实现了VTI介质中准P波方程的叠前逆时深度偏移．各向异性Marmousi模型的试算结果表明,VTI介质准P波方程叠前逆时深度偏移算法不受地下构造倾角和介质横向速度变化的限制,对复杂模型具有良好的成像能力;应用归一化互相关成像条件能得到更好的成像效果．对比该模型的各向异性和各向同性逆时偏移剖面表明,在各向异性地区采集的纵波数据用各向异性偏移算法理论上能得到更好的成像结果．      

盐丘模型 推覆体模型 波动方程 叠前深度偏移

三维波动方程叠前深度偏移是复杂介质中进行构造成像、弹性参数反演的重要环节.由于其技术实现不仅涉及波场延拓理论的创新,而且需要大规模计算,因而研究难度较大. 本文以实验效果的取得为目的,完整地实现了SEG/EAEG盐丘和推覆体模型的三维波动方程辛几何算法的叠前深度偏移成像计算. 文中详细考察了所研制的波动方程三维叠前深度偏移软件系统及其对复杂地质构造的成像能力. 具体包括：1）对于盐丘模型,文中讨论了成像参数的选择、地震子波对成像精度的影响、完成二维及三维叠前深度偏移的比较;2）对推覆体模型,文中进行了脉冲响应测试;3）由两个模型的成像结果可见本文的波动方程三维叠前深度偏移软件系统已具有适应强速度横向变化、复杂构造的成像能力.     

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

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