F-X域粘弹性波动方程保幅偏移

提出了一种利用粘弹性声波波动方程进行偏移的新方法。其基本思路是,修改成像条件,使修改后的成像方程中考虑振幅补偿,然后利用粘弹性、单程波动方程在F-X域对震源和接收点波场进行延拓,进而完成散射、透射等补偿。理论数据处理和实际地震资料处理表明,该方法理论基础可靠,处理效果明显,能较好地解决实际地质问题。     

三维叠前波动方程共偏移距拟屏深度偏移

在中点偏移距坐标系中，横向双平方根(DSR)方程为使用屏传播算子研究三维叠前波动方程深度偏移提供了一种方便的框架。共偏移距拟屏深度偏移是偏移共偏移距、共方位角地震资     

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

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

基于波动方程的聚焦点控制照明叠前深度偏移

 基于波动方程的聚焦点控制照明叠前深度偏移技术借助于差分计算，把速度、密度等介质参数的影响体现在差分计算的矩阵方程中，能够自动适应速度场的任意变化，快速傅里叶变换的使用也加速了波场延拓的计算速度。因此此法兼具有限差分偏移方法和傅里叶偏移方法的优点，既可适应速度场的剧烈变化，又可保证对陡倾地层的成像效果，是目前针对复杂构造最有效的成像方法之一。对于单个聚焦点及其周围的成像步骤为：①采用矩形网格情况下绕射走时的有限差分计算方法生成合成算子；②应用合成算子来合成面炮震源和面波记录；③对合成的面炮震源和面波记录做傅里叶有限差分法波动方程叠前深度偏移，得到该聚焦点及其附近区域的成像结果。按照上述成像步骤，将震源波场和炮集记录依据相应的外推公式进行延拓，最终应用成像条件求取成像值。在地质目标处选取多个聚焦点，可以得到面向目标的控制照明偏移成像，在多个层位上选取多个聚焦点进行控制照明叠前深度偏移，可以得到整块区域的成像。通过对Marmousi模型的试算，取得了较好的效果。     

波动方程保幅地震偏移成像方法

地震数据中饱含有丰富的旅行时信息和振幅信息。传统的波动方程地震偏移成像方法忽略了介质参数变化对地震波振幅的改造,不能准确描述地震渡传播的动力学特征。为此,基于严格的解耦理论,从全声波方程出发,进行单程波保幅分解,得到直观、高效率的的单程波动方程,直接面对地震波传播波场的压力分量进行延拓,从而将振幅误差补偿作为偏移的一部分在偏移过程中进行,实现地震偏移成像在获取正确位置的同时也获取真实振幅的信息。数值试算表明,该方法可使散射能量聚焦、归位,提高成像精度;可输出正确反映地下反射系数的振幅信息。     

波动方程保幅地震偏移成像方法

地震数据中饱含有丰富的旅行时信息和振幅信息。传统的波动方程地震偏移成像方法忽略了介质参数变化对地震波振幅的改造,不能准确描述地震渡传播的动力学特征。为此,基于严格的解耦理论,从全声波方程出发,进行单程波保幅分解,得到直观、高效率的的单程波动方程,直接面对地震波传播波场的压力分量进行延拓,从而将振幅误差补偿作为偏移的一部分在偏移过程中进行,实现地震偏移成像在获取正确位置的同时也获取真实振幅的信息。数值试算表明,该方法可使散射能量聚焦、归位,提高成像精度;可输出正确反映地下反射系数的振幅信息。     

保幅单程波叠前深度偏移中成像条件的对比分析

在波动方程炮域偏移中,成像条件的应用对于提供照明补偿和振幅恢复十分重要。尤其是对于波动方程保幅偏移,稳定的成像条件对于有效地恢复介质反射系数必不可少。通过对3界面水平层状模型和Marmousi模型试算,对比分析了各种成像条件,并依据成像结果对各种成像条件做出了评价。测试出最稳定的成像条件是基于平滑窗函数的,使用这种成像条件提高了成像的稳定性,明显改善了保幅效果,而且能够很好地压制偏移噪声。     

在工作站上实现的3D保幅叠前深度偏移方法

本文提出了一种基于射线+波恩近似的3D保幅叠前深度偏移算法(PAPsDM)。射线+波恩近似反演能够定量地恢复模型的扰动。格林函数则采用波前构造法(WFC),通过3D非均匀平滑速度场的动态射线追踪来估算。这种算法是在SunSparc 20工作站上实现的,工作站的CPU效率和内存需求应引起特别注意。我们将这种方法应用于一个3D海上数据体(13GB),获得一个7×1×1km的偏移成像结果大约需一周的CPU时间。     

分步傅立叶波动方程三维叠前深度偏移算法及应用

把分步傅立叶波场外推算子推广到三维叠前深度偏移成象,得出频率波数域与空间频率域混合域中运算的三维波场延拓算子,提出一种基于波动方程波场延拓的炮域三维叠前深度偏移算法.用SEG/EAGA盐丘模型标准数据试算,获得了良好的成象效果.对胜利探区埕北地区的海上实际资料进行了试处理,取得了良好的效果.     

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

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

波动方程叠前深度偏移角度道集求取方法

针对常规直接转换法和直射线叠前时间偏移(PSTM)法提取的角度道集无法满足地质构造复杂地区的AVA分析与反演精度的问题,研究基于共炮检距道集波动方程保幅叠前深度偏移方法,提取角度域共成像点道集。首先,对地震数据进行叠前深度偏移得到炮检距域共成像点道集;其次,在频率-波数域对炮检距域共成像点道集利用快速插值映射法抽取角度域共成像点道集。通过与直接转换法和直射线PSTM法提取的角度道集比较,分析不同方法求得的角度道集信息的准确性和AVA特征。数值试算和实际资料处理结果表明:共炮检距道集波动方程保幅叠前深度偏移方法求取的角度道集的振幅相对保真,角度范围更广,更适用于AVA分析。     

山前带逆掩构造保幅FFD叠前深度偏移

山前带逆掩构造区上覆高速地层对地震波传播能量具有比较严重的屏蔽作用,因而导致逆掩下伏构造成像比较困难。如今广泛采用的Kirchhoff积分法叠前深度偏移由于其自身算法的限制难以对其精确成像。基于此,本文引入了保幅FFD(傅里叶有限差分)叠前深度偏移方法,通过在传统延拓算子的基础上增加振幅恢复项求得保幅延拓算子,并给出保幅偏移的成像条件,从而实现了保幅FFD叠前深度偏移。数值试验及实际资料处理结果均表明:该方法能够有效补偿地震波场传播中的能量损失,提高逆掩下伏构造的成像质量。     

3D VSP波动方程深度偏移成像方法研究

介绍了一种基于波场外推的波动方程深度偏移成像方法。该方法不仅考虑了波场在地下传播过程中的绕射和折射效应,而且没有Kirchhoff方法中的高频近似假设和在复杂构造区的焦散及多路径现象。该方法在胜利垦71井区3DVSP实际资料处理中取得了比较好的成像效果。     

波动方程可变参数偏移方法

在波动方程偏移方法中，采用保种计算方法是提高地震资料处理效果的关键，而在计算方法中使用何种差分方程取代偏微分方程则是问题的核心。笔者在分析现有波动方程偏移方法“ＷＥＭＩＧ”程序过程中发现，用差分方程近似取代偏微分方程时，不考虑地质构造的特点，一律采用固定不变的权重比例系     

三维波动方程P-R分裂偏移

本文提出一种新的有限差分三维波动方程一步偏移方法,称为P-R分裂法。该方法是将三维偏移方程在τ(或ζ)方向分两步分别以二分之一步长进行延拓。在前个步长内将关于x的二阶偏导数取成隐格式而将关于y的二阶偏导数取成显格式;为了保持对称性,在后半个步长内,将关于y的二阶偏导数取成隐格式而将关于x的二阶偏导数取成显格式,经一定的简化后,可以得到一个在形式上非常简单的三维差分方程。该差分方程具有两个重要特点：①比Clearbout分裂方法差分精度高,频散现象小;②该差分方程可以局部地分解成一系列二维求解方程,它比马在田四点平均算法求解方便,计算量小。实际处理结果表明,该方法能使倾斜反射波正确归位,断层和断点显示清晰可靠,从而在很大程度上改善了三维偏移剖面的质量。     

矩形网格三维显式差分偏移

由Hale提出的基于McClellan变换的三维显式差分偏移算法能够对速度横向变化剧烈的陡倾角复杂介质偏移成像，该算法计算效率高，但只适用于方形网格上采集的资料。对于矩形网格的地震记录，通常需将其先插值到方形网格上，然后再应用基于McClellan变换的三维显式差分偏移算法进行偏移成像，这样既降低了计算效率，又增加了对内存的需求。为此，提出了一种修正McClellan变换的三维显式差分偏移算法，即将插值嵌于外推算子中，无需对数据重采样，可以直接对矩形网格资料进行偏移成像。理论分析与数值实验结果表明，应用修正方法偏移成像的结果与基于数据插值的显式偏移算法结果相吻合，但计算量显著降低。     

The application of amplitude-preserved processing and migration for carbonate reservoir prediction in the Tarim Basin, China

Conventional seismic exploration method based on post-stack data usually fails to identify the distribution of fractured and caved carbonate reservoirs in the Tarim Basin,so the rich pre-stack information should be applied to the prediction of carbonate reservoirs.Amplitude-preserved seismic data processing is the foundation.In this paper,according to the feature of desert seismic data (including weak reflection,fast attenuation of high frequency components,strong coherent noises,low S/N and resolution),a set of amplitude-preserved processing techniques is applied and a reasonable processing flow is formed to obtain the high quality data.After implementing a set of pre-stack amplitude-preserved processing,we test and define the kernel parameters of amplitude-preserved Kirchhoff PSTM (pre-stack time migration) and subsequent gathers processing,in order to obtain the amplitude-preserved gathers used to the isotropic pre-stack inversion for the identification of caved reservoirs.The AVO characteristics of obtained gathers fit well with the synthetic gathers from logging data,and it proves that the processing above is amplitudepreserved.The azimuthal processing techniques,including azimuth division and binning enlargement,are implemented for amplitude-preserved azimuthal gathers with the uniform fold.They can be used in the anisotropic inversion to detect effective fractures.The processing techniques and flows are applied to the field seismic data,and are proved available for providing the amplitude-preserved gathers for carbonate reservoir prediction in the Tarim Basin.     

三维一步法全倾角偏移

本文采用时间层分裂法实现三难一步法全倾角滤动方程有限差分偏移。该方法对和构造走向及倾角都无特殊限制。文中阐述了差分议程理论及实现方法，并对山区三维、煤田三维以及海陆交互区三维地震数据进行了处理。结果表明，本方法的偏移成像剖面比两步法偏移成像剖面分辨率高、界面归位合理、断点清楚，取得了明显的效果。     

频率域电磁场偏移成像

电磁场偏移是研究观测电磁场的上行波场，是以反向扩散传播方式使各个场源逐渐归位的。若向下传播的偏移场的相位与原来上行波场的相位在介质的分界面上相同，则通过电磁场偏移成像能够地球内部反射场的相位，本文讨论的问题是建立在激发源为平面电磁波基础之上的，当然此地方也可以应用于一般地面激发的可控制问题，文中采用有限差分方法进行偏移成像，不仅可以确定地下电阻率异常的位置，而且还能反映电阻率异常的相对大小，理论试