双复杂条件下非倾斜叠加精确束偏移方法及应用I——声波方程

近年来,油气勘探的重心正转向具有复杂地表和复杂地质体的双复杂区域.本文发展了一种精确的双复杂条件下基于地表倾角信息的非倾斜叠加束偏移方法,相对于传统束成像方法无需进行三方面处理:(1)高程静校正;(2)相位校正;(3)束中心与接收点之间关于速度和束出射角的近似替换,因而具有更高的成像精度.通过加拿大逆掩断层模型、中原油田断层模型及实际资料的偏移试算,并与传统束偏移及波动方程偏移成像结果对比可知:本文非近似束偏移方法在近地表、高陡倾等构造处的成像精度、反射界面成像振幅等方面优于传统的偏移方法,以此验证了本文非倾斜叠加精确束偏移方法的正确性、优越性及适应性.     

基于有效邻域波场近似的起伏地表保幅高斯束偏移

随着我国陆上地震勘探向复杂地表探区的转移,高精度、适应性强的地震成像方法在地震资料的处理、解释及后续属性分析、储层预测中具有重要意义.本文基于有效邻域波场近似理论发展了一种成像精度更高且适用于复杂起伏地表条件的叠前保幅高斯束偏移方法.在传统水平地表高斯束偏移的基础上,本文根据中心射线附近有效邻域内高斯束表征的近似波场,导出了起伏地表条件下具有相对振幅保持的高斯束偏移公式,并给出了一种精度更高的旁轴射线传播角度计算方法.同现有的高斯束偏移方法相比,本文方法不仅考虑了起伏地表对高斯束走时的线性影响,而且首次引入了由地表高程差异和近地表速度变化引起的二次时差校正项和振幅校正项,使得成像结果更加准确可靠.两个典型模型算例验证了本文方法的正确性和有效性.     

复杂地表条件下叠前菲涅尔束偏移方法

高斯束偏移虽然克服了Kirchhoff偏移不能处理多波至和单程波动方程偏移不能对陡倾构造准确成像的问题,但在复杂地表条件下其偏移精度取决于所选择的初始束宽度,即当初始宽度较小时,近地表成像精度较高,但此时中深层成像质量较差;反之当初始宽度较大时,中深层成像质量提高,但近地表成像精度降低.针对高斯束偏移中深层和浅层成像精度的矛盾,本文发展了一种适用于陆地复杂地表条件的叠前菲涅尔束偏移方法.基于惠更斯-菲涅尔原理,本文首先给出了菲涅尔束的概念及其表征的格林函数,并采用有效邻域波场近似理论和反褶积成像条件,导出了复杂地表条件下叠前保幅深度偏移公式.最后,针对常规旁轴射线追踪中的数值噪音,给出了一种压制策略.同高斯束偏移相比,本文方法不仅解决了中深层和浅层成像精度的矛盾,而且提高了复杂地表条件下平面波的分解精度,使得偏移结果更加准确可靠.典型的模型算例验证了本文方法的有效性和稳健性.     

高斯束成像方法研究进展综述

兼顾计算效率和成像精度的高斯束偏移成像方法近年来得到广泛应用.一方面高斯束偏移使用动力学射线追踪,解决了传统射线方法的焦散问题,改善了成像效果;另一方面,高斯束方法选择对成像点有贡献的高斯束叠加来计算波场值,无需费时的两点射线追踪,保证了计算效率.本文主要从高斯束叠后偏移、叠前偏移、起伏地表偏移、高斯束逆时偏移、弹性波偏移、各向异性偏移以及高斯束偏移三类优化策略等方面系统地评述了高斯束成像方法的研究进展,并对高斯束偏移方法下一步发展方向进行了初步展望.     

复杂地表条件下保幅高斯束偏移

高斯束偏移是一种准确、灵活、高效的深度域成像方法,其不但具有接近于波动方程偏移的成像精度,还保留了Kirchhoff偏移灵活、高效的特点以及对复杂地表条件良好的适应性.本文提出了一种适用于复杂地表条件的且具有相对振幅保持特点的高斯束偏移方法.通过考虑地表高程、倾角以及实际的道间距等信息,推导了基于高斯束表示的波场反向延拓公式,并结合反褶积成像条件,得到了复杂地表条件下的共炮域保幅高斯束偏移公式.同原有方法相比,本文方法不但可以直接在起伏的地表面进行局部平面波的分解,具有更高的成像精度,而且可以得到反映地下随角度变化反射系数的成像结果.数值模型的试算验证了上述结论.     

复杂地表条件下高斯波束叠前深度偏移(英文)

在复杂地表条件的区域,地震数据的采集和处理是一项极大的挑战。虽然可以通过静校正来消除起伏地表的影响,然而当地表高程以及近地表速度剧烈变化时,简单的垂直时移对地震波场造成的畸变会严重降低偏移成像的质量。基于射线的偏移方法可以直接在起伏地表面进行波场的延拓成像,是解决上述问题的有效手段。本文针对复杂地表条件下的高斯波束叠前深度偏移进行研究,对倾斜叠加公式进行修改,使之包含地表高程以及速度的信息,通过直接在复杂地表面进行平面波的合成,得到了一种具有更高成像精度的改进方法。首先简单介绍常规高斯波束偏移的基本原理和计算流程,并以此为基础,给出复杂地表条件下高斯波束偏移原有的实现方法以及本文的改进方法,最后通过模型和实际资料的试算验证本文方法的有效性。     

起伏地表QR径向基函数有限差分及其在弹性波逆时偏移中的应用

弹性波逆时偏移不受倾角和偏移孔径的限制,能够实现任意复杂构造的高精度多波成像,是目前最精确的多分量资料偏移成像方法之一.逆时偏移算法的核心是波场延拓,传统波场延拓以水平基准面为边界条件,基于固定采样步长进行规则网格剖分,采用阶梯近似法处理起伏地表和复杂构造界面时会产生台阶散射,严重影响起伏地表复杂构造的成像精度.基于无网格节点模型,定量分析了弹性波模拟中径向基函数有限差分法的频散关系和稳定性条件.基于此,提出一种基于QR径向基函数的高精度有限差分方法,并提出一种优化的起伏地表自适应节点剖分方法,推导了精确的无网格自由边界条件和弹性波无网格混合吸收边界条件,形成了新的基于无网格的起伏地表弹性波数值模拟方法.此外,本文将此无网格径向基函数有限差分方法应用于精确的纵横波场矢量分解公式,实现了起伏地表弹性波逆时偏移成像.通过对高斯山丘模型,起伏凹陷模型和起伏地表Marmousi-2模型进行数值试算,验证了本文方法的有效性和可行性.     

基于汉宁窗函数滤波时间域高斯束成像方法

实际地震数据中通常存在背景噪声及人为干扰信号,传统射线类方法对低信噪比数据成像适应性较差,基于此本文发展了一种适应于低信噪比数据的时间域高斯束成像方法.相比于频率域类束方法,新方法的主要优点包括:(1)频率域高斯束偏移方法需计算每个成像点处各个频率的格林函数,基于时间域高斯束偏移方法无需上述过程,大大提高了计算效率;(2)时间域高斯束公式推导中,直接对地震记录进行局部倾斜叠加并加入汉宁窗滤波处理,对低信噪比数据具有更好的适应性.在实现成像方法的基础上,本文通过对不同信噪比Marmousi炮记录进行模型试算及含噪声实际资料成像试处理表明:时间域高斯束偏移方法能够较为明显地提高低信噪比数据的偏移效果,提高成像方法对实际数据的适应性.     

接收函数的克希霍夫２D偏移方法

本文将地震反射资料成像的Kirchhoff偏移方法引入到接收函数成像研究，提出了一种接收函数的Kirchhoff偏移方法，以适应介质速度的横向变化，提高接收函数成像的精度和分辨率. 模型检验表明，与传统的水平分层的共转换点偏移叠加方法相比，Kirchhoff偏移方法能够实现转换波的有效归位，消除水平分层共转换点偏移叠加造成的假象. 一个很重要的现象是，倾斜或弯曲界面的成像明显受控于远震的入射方向，其成像精度的提高有赖于不同方向接收函数的联合成像.     

起伏地表地震资料成像方法分析

西部地区地表起伏剧烈,往往伴随着地下构造复杂,给地震偏移处理工作带来了极大的困难,导致后续的解释结果与实际构造存在一定的差异.因此,复杂地区的地震资料成像问题显得尤为重要.文章通过对复杂地表模型的正演模拟以及山区实际地震资料的应用两方面的分析,获得了适合复杂地表的成像方法.结果表明:针对复杂地区,叠加技术成像比偏移技术成像更能真实地反映地下构造形态,是一种提高复杂地区构造成像精度的有效方法;同时,将叠加剖面和偏移剖面两者有机结合,对复杂地表资料进行解释分析,更能体现地下构造的真实形态,为煤矿工作的安全开采提供了有利依据.     

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

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

基于局部倾斜叠加的快速射线束叠前时间偏移

宽方位高密度地震勘探可以有效的提高地震资料的空间分辨率和裂缝预测的精度,但地震道数的增加也大幅度提升了地震数据的处理成本.为提高海量地震数据偏移处理的计算效率,本文发展了一种快速射线束叠前时间偏移方法.该方法首先根据给定的射线束中心间隔将炮记录划分为一系列数据子集,然后利用倾斜叠加将数据子集分解为不同方向的平面波,最后根据射线束中心到地下成像点的双程走时和射线参数拾取相应的平面波振幅累加到成像点上.同Kirchhoff叠前时间偏移相比,本文方法不但保持了一致的成像精度,且由于仅需在稀疏的射线束中心位置进行成像累加运算,计算效率得到了大幅度的提升.文中所给出的模型和实际资料的测试结果验证了本文方法的正确性和有效性.     

三维共偏移距叠前地震数据反假频射线束形成偏移成像

对稀疏/非规则采样或者低信噪比数据,射线束提取困难并伴随有假频产生,对叠加剖面和道集造成严重干扰.为了提升射线束偏移在稀疏和低信噪比地震数据采集中的成像效果,本文提出基于三角滤波的局部倾斜叠加波束形成偏移假频压制方法.射线束偏移首先将地震数据划分为超道集,经过部分NMO后转化为以射线束中心定义的共偏移距数据,倾斜叠加和反假频操作均在局部共中心点坐标上实现.时间域倾斜叠加是对地震数据的时移累加操作,三角低通滤波同样可以在时间域完成,在对地震数据进行因果和反因果积分后,亦为地震数据的时移累加.因此,三角低通滤波与倾斜叠加可在时间域结合同时完成,避免了频域滤波的正反傅里叶变换.本文在反假频公式中加入权重系数,用以对反假频的程度进行控制,达到分辨率和噪声压制的最佳折衷.以某海上三维实际数据为例,文中展示了反假频射线束形成对偏移叠加剖面和共成像点偏移距道集中的噪声进行了有效压制.     

基于起伏地表的层析速度反演方法研究(英文)

叠前深度偏移速度分析是地震数据处理方法研究的重点,是影响地震成像效果的关键技术之一。对于地表起伏、地下构造复杂的双复杂地区,常规的叠前偏移速度分析方法是将起伏面校正到固定基准面上或进行表层建模再进行偏移与速度分析。本文提出的基于起伏地表的层析速度反演方法从起伏地表直接进行速度场的更新,可以提高层析速度反演的精度与效率。首先介绍基于起伏地表的角度域共成像点道集的提取方法,以此为基础阐述了起伏地表层析速度反演方法。起伏地表模型和实际资料试算验证了本文方法的有效性。     

三维最小二乘弹性高斯束叠前深度偏移

与声波高斯束成像相比,弹性波高斯束偏移更适用于复杂油气藏多波多分量地震勘探.但是由于观测系统的局限性和深部构造的复杂性,该方法同样存在成像分辨率低、照明不均衡等问题.本文结合最小二乘偏移和弹性高斯束偏移的优势,提出了一种通过弹性高斯束叠加构建Born正演(反偏移)算子和偏移算子的三维最小二乘叠前深度偏移方法.依据最小二乘反演理论,建立基于反偏移数据与实际观测数据残差的目标函数,采用共轭梯度算法迭代更新来建立地下真实的反射率.与传统弹性高斯束偏移方法相比,该方法不仅提高了成像分辨率,而且使复杂构造特别是陡倾角地层的成像照明也得到了补偿.理论模型测试结果证明了本文方法的可行性和有效性.     

双参数展开CRP叠加和速度分析方法研究

椭圆展开共反射点(CRP)方法可以获得比常规倾角时差校正(DMO)方法更近似的零偏移距时间剖面和相应CRP速度场.大量研究和实践证实,在非均质性较弱的地区,该方法取得的成果显著.但由于该方法没有考虑速度的横向变化和转换波等情况,当地下介质存在较强非均质性时,该方法不再准确,需要引进反映速度横向变化的双参数(上行波与下行波的平均速度和速度比)进行改进.本文详细推导了引入双参数后的叠加和速度分析算法,并通过数值模型和地震资料处理证实,修正后的算法可以更好地解决地质复杂地区速度建模和叠加成像问题.     

非零井源距VSP逆高斯束叠加成像

VSP-CDP叠加成像是非零井源距VSP数据成像的重要技术手段,但是当地下构造比较复杂时,基于水平层状介质假设的叠加成像方法误差较大.本文使用基于高斯射线束的叠加成像方法对复杂构造VSP数据进行成像,根据高斯束有效邻域波场近似理论提出逆高斯叠加权函数,并使用该权函数对VSP数据进行VSP-CDP转换.同常规的VSP-CDP叠加成像方法相比,本文方法不仅解决了复杂构造非零井源距VSP数据成像问题,而且可有效加密反射点覆盖次数从而提高成像的横向分辨率.模型算例及实际资料验证了本文方法的有效性和稳定性.     

波动方程基准面延拓研究进展

复杂的地表采集条件是地震偏移成像必须面对的问题．地形高程起伏变化较大,以及近地表低速层速度的横向变化,是导致反射地震数据非双曲时差特征的主要原因,严重影响了后续速度估计和偏移成像的精度．对于复杂地表条件,基于地表一致性的时移法静校正很难满足高精度资料处理的要求．波动方程基准面延拓方法不受地表一致性限制,很好地保持了波传播的动力学特征,而且在改善成像质量方面效果显著,应用广泛．本文探讨了消除地表因素影响的基准面延拓方法,对一些意义比较突出的方法进行了详细介绍．     

TI介质局部角度域高斯束叠前深度偏移成像

各向异性射线理论基础上的局部角度域叠前深度偏移方法能够为深度域构造成像与基于角道集的层析反演提供有力支撑,但是对于复杂地质构造而言,高斯度叠前深度偏移在不失高效、灵活等特点的情况下,具有明显的精度优势.为此,本文研究局部角度域理论框架下的高斯束叠前深度偏移方法.为提高算法效率与实用性,文中讨论了一种从经典弹性参数表征的各向异性介质运动学和动力学射线方程演变而来的由相速度表征的简便形式,并提出了一种比较经济的各向异性高斯束近似合成方案.结合地震波局部角度域成像原理,讨论一种适合高斯束偏移的角度参数计算方法.国际上通用的理论模型合成数据试验表明：相比局部角度域Kirchhoff叠前深度偏移成像方法,本文方法具有更高的成像精度与抗噪能力,既适用于复杂构造成像,也可为TI介质深度域偏移速度分析与模型建立提供高效的偏移引擎.     

输出道成像方式的共反射面元叠加方法I——理论

共反射面元（Common Reflection Surface）叠加是一种独立于宏观速度模型的零偏移距成像方法，该方法属于典型的克希霍夫型成像方法. 根据成像方式的不同，克希霍夫型成像方法可以分为两大类：输出道成像方式和输入道成像方式. 考察共反射面元叠加方法，它属于输入道成像方式. 本文基于理论模型数据，实现了输出道成缘方式的CRS叠加方法. 相比传统的输入道成像方式，它具有能够保证大偏移距反射信息的成像精度和计算效率较高的优点，而且更加容易推广到三维情形.