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

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

Kirchhoff型真振幅反偏移加权函数分析

Kirchhoff型真振幅反偏移的基本目的是通过等时线叠加得到偏移的输入场,是将深度偏移的结果再反映射到时间域去,而真振幅加权函数是Kirchhoff型真振幅反偏移的核心内容,它没有唯一的数学形式,目前发表在文献当中的加权函数公式是在物理直观下的最优表述.本文在前人研究的基础上推导了真振幅反偏移加权函数公式,对影响加权...     

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

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

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

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

几类叠前深度偏移方法的研究现状

叠前深度偏移是解决复杂地质构造成像的有力手段,关于叠前深度偏移方法的研究一直是勘探地球物理研究的热点和难点。本文在收集、整理国内外有关研究资料的基础上,有选择地对若干类型方法研究背景及现状、基本原理及方法特点等进行分析,目的在于作为研究工作的起点,把握叠前深度偏移方法的发展趋势,以期生成适合我国复杂油气储集特点的地震成像方法及技术手段。     

复杂观测系统下的三维波动方程叠前深度偏移

波动方程三维叠前深度偏移是近年应复杂地质构造与地震岩性成像需求而发展的一项关键技术。此项技术与集群式并行机的结合，更是将其价格性能比较低至工业生产规模应用水平，然而，波动方程三维叠前深度偏移的实用化还需要诸如地震资料网络化、并行计算负载平衡等一系列配套技术。本文针对油田实际资料，试验了其应用波动方程三维叠前深度偏移的规则化技术，并结合地震炮集数据的特点，在自组装的集群式并行机上，解决了其节点间的负载平衡问题。本项工作有助于推近复杂观测系统下的地震数据实现三维波动方程叠前深度偏移。     

三维非均匀介质中真振幅地震偏移算子研究

利用三维非均匀介质中的波动方程,进行振幅保真波场偏移算子分解,得到用于真振幅偏移的单程波方程. 经过数学推理,导出裂步Fourier法真振幅偏移和Fourier有限差分法真振幅偏移的算子方程,并给出其具体的实现过程.     

共聚焦点叠前深度偏移速度估计方法研究

基于共聚焦点技术可以将叠前深度偏移分为激发聚焦与接收聚焦两个过程,通过逆时聚焦算子的作用实现.由速度模型驱动逆时聚焦算子的更新,合成DTS面板,根据等时原理,通过DTS面板中时移曲线的曲率变化来交互的更新速度场.但是当地下介质稍为复杂的情况下,从面板中拾取分析点对应的时移曲线变得非常困难.针对这个问题,本文通过对分析点在限定深度和速度范围内合成的DTS面板零时刻振幅进行叠加,由其值的变化来确定分析点的合理位置及上层介质速度.     

波动方程叠前深度偏移及并行计算在南黄海地区的应用

近年来发展起来的微分解型叠前深度偏移成像技术，在复杂油气勘探的实践中不断完善，已经成为复杂地质目标偏移成像必不可缺的工具．同时，随着计算机技术的飞速发展，使得并行算法能够极大地提高数值计算的效率．本文以南黄海地区为例，描述了运用二维波动方程叠前深度偏移辛几何算法及并行运算的过程并给出结果．     

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

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

Dreamlet source‐receiver survey sinking prestack depth migration

Survey sinking migration downward continues the entire surface observed multi‐shot data to the subsurface step by step recursively. Reflected energy from reflectors at current depth appear at zero time and zero offset in the extrapolated wavefield. The data (seismic records) of t > 0 at this depth are equivalent to the data acquired by a survey system deployed at this depth. This is the reason to name the process ‘survey sinking’. The records of negative time need not to be further propagated since they carry no information to image structures beneath the new survey system. In this paper, we combine survey sinking with dreamlet migration. The dreamlet migration method decomposes the seismic wavefield and one‐way wave propagator by complete time‐space localized bases. The localization on time gives flexibility on time‐varying operations during depth extrapolation. In dreamlet survey sinking migration, it only keeps the data for imaging the structures beneath the sunk survey system and gets rid of the data already used to image structures above it. The deeper the depth is, the shorter is the valid time records of the remaining data and less computation is needed for one depth step continuation. For data decomposition, in addition to time axis, dreamlet survey sinking also decomposes the data for source and receiver gathers, which is a fully localized decomposition of prestack seismic data. A three‐scatter model is first used to demonstrate the computational feature and principle of this method. Tests on the two‐dimensional SEG/EAGE salt model show that with reduced data sets the proposed method can still obtain good imaging quality on complex geology structures and a strong velocity contrast environment.     

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

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

Investigating a time‐shift extended imaging condition in a Kirchhoff pre‐stack depth migration algorithm

Extracting true amplitude versus angle common image gathers is one of the key objectives in seismic processing and imaging. This is achievable to different degrees using different migration techniques (e.g., Kirchhoff, wavefield extrapolation, and reverse time migration techniques) and is a common tool in exploration, but the costs can vary depending on the selected migration algorithm and the desired accuracy. Here, we investigate the possibility of combining the local‐shift imaging condition, specifically the time‐shift extended imaging condition, for angle gathers with a Kirchhoff migration. The aims are not to replace the more accurate full‐wavefield migration but to offer a cheaper alternative where ray‐based methods are applicable and to use Kirchhoff time‐lag common image gathers to help bridge the gap between the traditional offset common image gathers and reverse time migration angle gathers; finally, given the higher level of summation inside the extended imaging migration, we wish to understand the impact on the amplitude versus angle response. The implementation of the time‐shift imaging condition along with the computational cost is discussed, and results of four different datasets are presented. The four example datasets, two synthetic, one land acquisition, and a marine dataset, have been migrated using a Kirchhoff offset method, a Kirchhoff time‐shift method, and, for comparison, a reverse time migration algorithm. The results show that the time‐shift imaging condition at zero time lag is equivalent to the full offset stack as expected. The output gathers are cleaner and more consistent in the time‐lag‐derived angle gathers, but the conversion from time lag to angle can be considered a post‐processing step. The main difference arises in the amplitude versus offset/angle distribution where the responses are different and dramatically so for the land data. The results from the synthetics and real data show that a Kirchhoff migration with an extended imaging condition is capable of generating subsurface angle gathers. The same disadvantages with a ray‐based approach will apply using the extended imaging condition relative to a wave equation angle gather solution. Nevertheless, using this approach allows one to explore the relationship between the velocity model and focusing of the reflected energy, to use the Radon transformation to remove noise and multiples, and to generate consistent products from a ray‐based migration and a full‐wave equation migration, which can then be interchanged depending on the process under study.     

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

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

非零偏VSP弹性波叠前逆时深度偏移技术探讨

非零偏VSP地震资料是一种多分量资料,处理非零偏VSP资料,弹性波叠前逆时深度偏移技术无疑是最适合的处理技术.本文从二维各向同性介质的弹性波波动方程出发,研究了对非零偏VSP资料进行叠前逆时深度偏移的偏移算法,讨论了逆时传播过程中的边值问题和数值频散问题及其相应的解决方案;采用求解程函方程计算得到地下各点的地震波初至时间作为成像时间,实现了非零偏VSP资料的叠前逆时深度偏移.最后进行了模型试算和非零偏VSP地震资料的试处理,结果表明该方法不受地层倾角限制,较适用于高陡构造地区或介质横向速度变化较大地区的非零偏VSP地震资料处理.     

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

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

TI介质局部角度域射线追踪与叠前深度偏移成像

研究与实践表明,对于长偏移距、宽方位地震数据,忽略各向异性会明显降低成像质量,影响储层预测与描述的精度.针对典型的横向各向同性(TI)介质,本文面向深度域构造成像与偏移速度分析的需要,研究基于射线理论的局部角度域叠前深度偏移成像方法.它除了像传统Kirchhoff叠前深度偏移那样输出成像剖面和炮检距域的共成像点道集,还遵循地震波在成像点处的局部方向特征、基于扩展的脉冲响应叠加原理获得入射角度域和照明角度域的成像结果.为了方便快捷地实现TI介质射线走时与局部角度信息的计算,文中讨论和对比了两种改进的射线追踪方法:一种采用从经典各向异性介质射线方程演变而来的由相速度表征的简便形式;另一种采用由对称轴垂直的TI(即VTI)介质声学近似qP波波动方程推导出来的射线方程.文中通过坐标旋转将其扩展到了对称轴倾斜的TI(即TTI)介质.国际上通用的理论模型合成数据偏移试验表明,本文方法既适用于复杂构造成像,又可为TI介质深度域偏移速度分析与模型建立提供高效的偏移引擎.     

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

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

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

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

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

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