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

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

基于时空局域化dreamlet单程波算子的观测系统沉降法偏移

Dreamlet偏移成像目的是探索一类能够对地震波场和单程波传播算子同时分解和压缩的理论和方法,也即实现在压缩域的传播与成像.地震数据在时间和空间的非平稳性质,决定了要实现地震数据的有效稀疏表示,分解方法必须在时间和空间上同时具有局域化性质.Dreamlet由时间和空间局部分解原子的张量积构成,可以看作一种脉冲-小波束形式的波场分解原子.时空局域化的dreamlet单程波传播算子在对波场沿深度方向延拓时,地震数据在时间轴上总是向同一方向流动.随着深度的增加,部分用于成像浅层结构的数据归位至其空间位置后被dreamlet算子丢弃,波场的有效记录时间变短,每一步用于波场延拓的计算量也相应下降.为了充分发挥这一优势,本文介绍dreamlet传播算子与观测系统沉降法偏移相结合的理论与方法.观测系统沉降法偏移每一步都将记录到的所有数据向下延拓,沉降后的波场等效于把源和检波器都放置在目标深度所能接收的反射数据.Dreamlet观测系统沉降过程只保留用于成像观测系统下部地质结构的有效数据,自动丢弃已经用于成像观测系统上部而对下部成像没有贡献的信号.本文通过二维SEG/EAGE叠后和Marmousi叠前数据算例展示了dreamlet传播算子应用于观测系统沉降法偏移的这一特点.数值算例结果显示,在不影响成像质量的前提下,该偏移方法能够有效减少传播数据量,为发展一种快速高效的偏移方法提供了新的思路.

     

波动方程三维叠前深度偏移并行计算流程探索

波动方程三维叠前深度偏移研究,已成为近年众多油气科研生产单位的研究焦点,其中关键问题之一是借助于软硬件及其匹配条件的进步,以求降低庞大的计算量障碍.本文借助于地震波场成像基本算法的内在并行性,运用集群式并行机结构的软硬件特点,提出了两种算法处理流程对策,以甚低的资金投入实现了在工业生产规模上的应用与效果.经油田实际资料计算,结果表明文中提出的分解算法具有甚高的加速比,良好的扩展性及适应性.     

三维叠前深度偏移的准三维算法研究

介绍了目前叠前偏移方法的研究状况，引入准三维算法的概念并对其进行了讨论，同时对共方位角偏移算子作了详细的介绍，最后，给出了运用共方位角偏移算子在实际中的应用结果，认为共方位角方法具有运算速度快，适应范围大的优点，有较好的研究开发前景。     

常用叠前深度偏移方法特点分析与实例对比

目前常用叠前深度偏移主要有四种:基于射线理论的Kirchhoff偏移和高斯束偏移以及基于波动方程理论的单程波偏移和逆时偏移.要想在实际应用中取得好的成像效果,必须选取合理、有效的叠前深度偏移算法.本文首先简单阐述了这四种常用叠前深度偏移方法的基本原理,对比分析了四种方法在陡倾角、速度依赖性、起伏地表、非规则数据、低信噪比数据等方面的优劣性以及计算效率和计算精度的差异,然后利用模型试算和实际资料处理分析说明这四种叠前深度偏移方法的成像效果,最后给出不同情况下叠前深度偏移方法的选取建议.     

补偿吸收衰减的高角度叠前深度偏移方法(英文)

实际介质的吸收衰减使地震波能量减弱,主频降低,影响分辨率。在偏移过程中,可以实现吸收衰减的补偿,提高地震资料的分辨率。本文基于高角度时间-空间域单程波方程,推导出能补偿介质吸收衰减的叠前深度偏移波场外推公式。设计了一个复杂构造的地质模型,用声波方程模拟了无吸收衰减和有吸收衰减介质存在时的两种地震炮集记录。用没有补偿吸收衰减和本文讨论的补偿吸收衰减的偏移方法分别处理了两种地震资料,得到无吸收衰减、吸收衰减未得到补偿和得到补偿的三个偏移叠加剖面。补偿后的偏移剖面与没有吸收衰减的地震记录的偏移剖面几乎相同,而有吸收衰减的地震记录用无吸收衰减补偿的偏移方法处理,由于介质的吸收衰减,使吸收衰减层下方的反射信号衰减,主频降低,能量减弱。实际资料的处理结果表明,补偿吸收的偏移剖面的分辨率得到了明显改善。     

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

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

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

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

三维各向异性介质中的波动方程叠前深度偏移方法

基于三维VTI各向异性介质的频散关系,构建波数项和空间项分离的单程波算子表达式,以优化算法,确定算子的待定系数,实现广角逼近三维VTI介质的广义相移算子,发展了可灵活处理强或弱各向异性介质的波动方程叠前深度偏移方法.文中同时也针对其工业应用建议了三维VTI各向异性介质中可提高计算效率的频率相关变步长波场深度延拓算法及稀疏采样情况下可实现陡倾角构造正确成像的反假频波场延拓算法.SEG二维Hess VTI各向异性介质理论模型数据及野外数据集观测系统下的三维脉冲响应计算表明,本文提出了一种具有工业应用潜力的三维各向异性波动方程叠前深度方法.     

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

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

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深度偏移

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

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

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

傅里叶有限差分法保幅叠前深度偏移方法

地震数据中饱含有丰富的走时信息和振幅信息. 为解决传统偏移方法中几何扩散和入射角变化引起的振幅误差问题,本文提出了一种实用的波动方程保幅地震偏移方法. 该方法从全声波方程出发进行单程波保幅分解,得到直观、高效率的直接面对地震波传播波场的压力分量进行延拓的保幅偏移单程波方程,进而推导出一个含有6项的傅里叶有限差分法保幅偏移的算子方程;修改边界条件和成像条件,使修改后的边界条件和成像方程中考虑振幅补偿,从而从三方面补偿几何扩散损失和入射角变化对振幅的影响. 脉冲响应测试、单炮记录的数值试验以及SEG/EAGE盐丘模型的叠前偏移结果表明,该方法不但可以使散射能量聚焦、归位,提高成像精度;而且可以输出正确反映地下反射系数的振幅信息,为后续的地震属性分析(如AVO/AVA)提供更真实的地震信息.     

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

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

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

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