TM法电磁相位移偏移研究

在地震勘探中，偏移是地震资料的主要处理方法之一。ＭＴ法中的平面电磁波与地震弹性波有本质的区别，但也有很多相似性。导电介质中的平面电磁波在介质的分界面上也发生反射与透射。利用平面电磁波的上行波场Ｕ和下行波场Ｄ在反射面上的相位关系特征，可把地震偏移技术中Ｕ／Ｄ成像理论引入到ＭＴ法数据处理中。相位移法原理简单、速度快，是一种理想的一维偏移方法，但在二维模型数据的偏移处理中处理中还需要作适当的改进。ＰＳＰ     

相移加校正叠前深度偏移及偏移速度分析

从理论上说，叠前深度偏移是众多的地震资料偏移成像方法中最精确的偏移方法。本文针对纵横向变速的复杂地质模型，采用相移加校正方法在ＣＳＰ道集上邮二维叠前深度偏移。该 以相移法煤中上局部速度变化因子校正，并用自动统计子波方法确定震源函数形式，因此该法既不受地层倾角的限制，又可解决复杂构造的横向变速问题。     

MT法电磁相位移偏移研究

在地震勘探中，偏移是地震资料的主要处理方法之一。MT法中的平面电磁波与地震弹性波有本质的区别，但也有很多相似性。导电介质中的平面电磁波在介质的分界面上也发生反射与透射。利用平面电磁波的上行波场U和下行波场D在反射面上的相位关系特征，可把地震偏移技术中U／D成像理论引入到MT法数据处理中。相位移法原理简单、速度快，是一种理想的一维偏移方法，但在二维模型数据的偏移处理中还需要作适当的改进。PSPI（相位移加内插）法可以得到二维U／D成像的较好结果。     

三维波动方程深度偏移和模型正演

构造复杂区的勘探和油田开发对三维地震工作提出了更高的精度要求。三维波动方程深度偏移是解决复杂构造、适度变化(纵、横向)剧烈地区地震波正确归位的唯一手段,它能够提供油气藏圈闭在三维空间的准确分布。利用三维波动方程模拟正演可以对地下地质模型进行数值模拟,为油藏描述提供有效的信息。本文运用分裂算法,将频率(ω)-空间(x)域的精度高、频散低、边界无反射的二维深度波场延拓法推广到三维波场延拓。在每一延拓深度步长内,分别用二维算子对CMP方向(x方向)和测线方向(y方向)依次进行延拓。由于采用了一步法波场延拓,所以这是一种适于速度纵、横向变化的三维叠后波场延拓的准确算法。它既可以用于三维模型的正演,也可用于构造的反演。     

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

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

共炮检距剖面相位移时间偏移

本文提出了一种新的二维频率一波数域共炮检距剖面叠前时间偏移算法。本文方法在常速条件下将零炮检距剖面的叠后相位移公式推广到非零共炮检距剖面，该公式是用数值拟合方法得到的，运用平均速度或叠加速度可使公式的使用范围适用纵向速度变化，理论上本文的公式只在常速下是精确的，但合成记录和野外资料试算表明，即使速度变化较大，该方法也能有效成像，而且更重要的是还能减少计算量。     

三维真振幅有限偏移距偏移

非零偏移距一次反射纵波可描绘成三维时间或深度偏移反射,以致偏移的波场振幅是与角度有关的反射系数的一种度量。最近提出的各种涉及加权绕射叠加的偏移/反演算法都是以Born或Kirchhoff近似法为基础的。本文提出一种三维Kirchhoff型叠前偏移法,该方法中待成象波场的一次反射是由零阶射线近似事先描绘的。因此,恢复与角度有关的反射系数的主要问题就变成了消除一次反射几何扩散因子的问题。为达到该目的而取的权函数不依赖于未知反射层,它能正确地恢复偏移图象中的震源脉冲,而不在乎所用的震源—检波器排列如何以及波场中是否出现焦散线。我们的权函数是用旁轴射线理论计算的,它比得上基于Beylkin行列式的反演积分函数。它们仅相差一个很容易解释的因子。     

相位移的分层模型算法

相位移的分层模型算法是以常途相位移算法为基本组成部分，运用非齐次单程波声波方程，把每一层的层速度作为常速，分层计算波场，即把下一层波场作为上一层的震源函数计算上一层的波场，如此由深至浅向上延拓，直到地面为止。这种方法既克服了常速相移法因速度突变而产生的误差，又保持了常速相移法精度高的优点，非常适合于速度纵横向突变的情况。文中给出的合成记录说明了该方法的可行性。     

三维VSP数据的波动方程偏移成像

常规的三维VSP地震数据偏移成像大都使用基于射线理论的Kirchhoff积分方法在共炮点道集中进行，因此对地下速度横向变化大的复杂构造区适应性差、效率很低。基于三维VSP地面激发井中接收观测系统的特殊性，利用地震波场的互易性原理，提出了一种在共接收点道集中进行三维VSP地震数据波动方程偏移成像的方法，不仅极大地提高了三维VSP偏移成像的计算效率，而且还适用于速度横向变化大的复杂构造区。为了利用时间域偏移成像方法对速度模型的相对不敏感特性，在三维VSP偏移成像中，借助于波动方程叠前深度偏移成像算法的概念，还提出了一种可以满足时间域速度模型横向变化要求的波动方程叠前时间偏移成像方法，拓展了波动方程叠前时间偏移成像的应用范围。     

相移插值法偏移中一种新的变延拓步长方法

在叙述相移插值法波动方程偏移原理、对该偏移过程中的波场外推公式及最终成像方程进行适当的改进变形后,重新标定了向下外推算子,使得成像过程可借助于 FFT 来完成,并给出了一种不同于时间域内插的新的变延拓步长的方法。本方法与传统的相移插值法偏移中时间域内插值法变延拓步长相比,向下外推步长可远远大干地震波视波长的四分之一,从而满意地解决了相位插值偏移算法效率低的问题。文中给出的理论模型计算与实际煤田及陡倾角地震数据偏移处理的结果表明:本方法是一种快速经济的,适合于陡倾角、复杂地质构造归位的 f-k 域偏移成像技术。     

三维一步法全倾角偏移

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

三维混合延拓一步法波动方程深度偏移

对于二维陡倾斜地层，相位移延拓方法成像精度高、稳定性好、速度快，但难以适应速度的横向变化；频率、空间域的有限差分算法简单、稳定性好、能适应速度的纵、横向变化，但成像精度低、运算速度慢；而用快速的45°有限差分法与相位移延拓法相结合，则能使整个延拓过程既适应速度的纵、横向变化，又能够使得陡倾角地层精确归位。本文把这种混合延拓思想拓展成三维相位移加有限差分混合延拓法，并进一步提出了能适应速度纵、横向变化，构造任意复杂介质的三维一步法深度偏移。数值模拟试验结果表明，该方法具有移速度快、精度高的特点。     

叠后三维一步法偏移处理

在实现地震资料的全三维处理中，叠后三维一步法偏移是至关重要的一个步骤。本文通过对一块浅海资料的一步法和两步偏移处理结果的对比，说明一步法偏移在改善陡倾角成像、确定断层位置和构造方面有良好的效果。     

三维VSP资料波动方程叠前深度偏移研究

已有的一些三维VSP资料处理方法，存在保幅性较差、对复杂构造不能正确归位、成像分辨率低、噪声强和计算效率低等问题，为此，提出了三维VSP资料波动方程叠前深度偏移成像方法，其实质是把地表数据的三维单程波动方程叠前深度偏移成像方法推广应用到三维VSP数据。对方法进行了阐述，其基本原理是将三维VSP数据按井中检波器的位置抽成共检波点道集，对共检波点道集进行三维叠前深度偏移。实现过程为将VSP波场分解为背景波场和扰动波场，背景波场采用相移法偏移，扰动波场采用频率空间域有限差分法偏移，然后对偏移后的上下波场进行相关处理，获得成像结果。数值模拟证明方法正确可行，实际应用表明方法实用。     

三维地震资料的全三维处理方法及其在江苏地区的应用

以往江苏油田对三维地震资料的处理，没有采用全三维处理方法，致使三地震资料处理的效果不理想。为此，我们研究了一套三维地震资料的全三处理方法软件，主要包括三地震预处理，三维速度分析和三波动方程空间域一步法偏移、将这套方法用于该区复杂地质构造地区的三地震资料处理，明显地提高了三维地震资料的处理精度和剖面的信噪比。     

横向各向同性介质中的相移法偏移

本文在非垂直对称轴横向各向同性（ＴＩ）介质情况下，引用描述平面声波传播的克里斯托弗尔方程，推导出ＴＩ介质中ＳＨ波，ｑＰ波和ｑＳＶ波的频散关系式，然后根据单程标量波动方程相移法偏移的原理，讨论了这三种波偏移向下延拓算子的解析表达式以及在ＴＩ介质中相移法偏移的实现方法，理论脉冲响应的计算证明了该方法的正确性。