P-SV反射转换波的倾角时差校正（DMO）方法研究

倾角时差校正（ＤＭＯ）方法对于复杂地质构造条件下的纵波反射波（Ｐ－Ｐ）的多次覆盖地震数据的叠加具有显著的提高信噪比的效果．对于转换反射波（Ｐ－ＳＶ）来说，进行倾角时差校正（ＤＭＯ）处理就更需要．本文首先对Ｐ－ＳＶ波的倾角时差校正（ＤＭＯ）的几何原理进行了分析，然后重点研究用波动方程进行ＤＭＯ的方法，即Ｐ－ＳＶ转换反射波的叠前部分偏移的波动方程法．这种方法也适用于Ｐ－Ｐ反射波的叠前部分偏移处理．研究结果表明，波动方程ＤＭＯ方法在振幅保真方面比几何原理的ＤＭＯ方法有显著改善．     

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

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

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

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

地震叠前深度偏移方法流程及应用

针对复杂介质的成像问题，提出了一套地震资料叠前深度偏移方法流程，主要包括三部分：（１）地震资料精细预处理；（２）速度－深度模型建立；（３）叠前深度偏移成像．以Ｋｉｒｃｈｈｏｆｆ偏移理论为基础，强调地质与地球物理的综合以及地震处理与解释的一体化．在ＺＸ地区成功地实现了二维地震资料叠前深度偏移，所获得的ＮＥ２０６叠前深度偏移剖面揭示了复杂的ＺＸ古潜山及其内部构造，并清楚地展示了逆掩断层的存在．     

磁共振脑血管成像方法的价值探讨

本文作者通过对１２３例病人同时进行３ＤＭＴＴＯＦ（三维磁化传递对时间飞逝法）和３ＤＰＣ（三维相位对比法）、磁共振脑血管成像（ＭＲＡ），对２１例患者同时进行了３ＤＭＴＴＯＦ及３ＤＴＯＦ（三维时间飞逝法）、脑血管成像，将它们对未稍小血管，颈外动脉，静脉显示程度相比较，并进行统计学分析，从而得出了３ＤＭＴＴＯＦ法较普通３ＤＴＯＦ法及３ＤＰＣ法在大脑中动脉岛叶分支、大脑后动脉的顶枕支，距状裂支的显示程度上优越，而且对静脉窦的抑制，颈外动脉及皮层静脉的抑制方面均优于后两者，并对在动脉瘤及动静脉畸形的诊断价值方面进行了探讨。提示３ＤＭＴＴＯＦ对小的动静脉畸形、海绵状血管瘤、小动脉瘤的诊断方面较３ＤＰＣ优越。３ＤＰＣ法对较大的动静脉畸形显示较好，对近年用于临床的ＭＴＴＯＦ法血管成像术的机理进行了探讨。指出这是一项新的磁化递抑制技术，应用于ＭＲＡ取得了较好的效应。     

基于波动方程的广义屏叠前深度偏移

地震波传播算子的计算效率和精度是制约三维叠前深度偏移的关键因素. 广义屏传播算子(GSP, Generalized Screen Propagator)是一种在双域中实现的广角单程波传播算子. 这一方法略去了在非均匀体之间发生的交混回响,但它可以正确处理包括聚焦、衍射、折射和干涉在内的各种多次前向散射现象. 通过背景速度下的相移和扰动速度下的陡倾角校正,广义屏算子能够适应地层速度的强烈横向变化. 这种算子可以直接应用于炮集叠前偏移,通过将广义屏算子作用于双平方根方程,还可以获得一种高效率、高精度的炮检距域叠前深度偏移方法,用于二维共炮检距道集和三维共方位角道集的深度域成像. 本文首先简述了炮检距域广义屏传播算子的理论,进而讨论了共照射角成像(CAI, Common Angle Imaging)条件,由此给出各个不同照射角(炮检距射线参数)下的成像结果,进而得到共照射角像集. 由于照射角和炮检距的对应关系,共照射角像集又为偏移速度分析和AVO(振幅随炮检距变化)分析等提供了有力工具.     

适于窄线三维地震资料的面炮方法

用小孔径波动方程叠前深度偏移完成大面积数据的成像，不仅总的计算量较大，而且成像畸变较严重。为了进行大孔径波动方程叠前深度偏移，必须将实际的小孔径资料合成大孔径资料，并解决大孔径资料的成像方法。AMO等合成孔径方法，数据合成方式简单，但成像方法是近似的且较复杂。面炮方法从理论上讲数据合成方法较简单，成像方法是精确的且较简单，但其应用适用于远道振幅比近道弱很多的情况。在窄线三维地震资料情况，其联络测线方向不易满足面炮应用条件。针对这种情况，本文将最小二乘原理应用于波场成像，提出适用窄线三维地震资料的面炮方法。为了解三维波动方程叠前深度偏移的面炮方法解决复杂油气藏的可行性、所需计算资源和成像效果，我们利用大庆SW地区三维地震资料进行了实验研究。通过研究，我们获得了面炮方法的实际效果，获得了所需计算资源的测试数据。对实际效果的分析对比表明，面炮方法是一种三维波动方程叠前深度偏移的快速方法，有较好应用前景。     

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

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

拉东投影法三维叠前深度偏移

对地下地质构造进行正确成像是地震勘探的最终目的,由于三维地震资料采集不可能都沿垂直构造走向的方向进行,为地震资料的三维处理带来了许多困难. 本文将三维叠后拉东投影偏移思想应用于三维叠前处理中,提出了三维叠前投影偏移算法. 利用拉东投影变换的原理,将整个三维叠前数据体投影到一系列各方向的径向线上,各方位角的构造都包含在其中某条或多条径向剖面上. 投影完成后,形成一系列的独立的二维叠前测线,可采用二维叠前深度偏移成像方法来实现各径向线的叠前偏移,当各径向剖面偏移完成后,在时间切片上进行反投影,从而最终形成三维叠前深度偏移结果. 实际应用表明,用本方法进行三维叠前深度偏移时,深度偏移剖面对横向分辨率有所提高,对陡地层和小断层的成像效果有所改善.     

特征高斯波包叠前深度偏移方法

高斯波包（Gaussian packet）传播算子可在局部时空域高效地计算局部波包的传播.高斯波包叠前深度偏移的基础是在Gabor变换域描述观测数据,再利用高斯波包传播算子计算炮点波场和检波点波场,两者相关即可得到偏移结果.利用炮道集的局部τ-p特征在Gabor变换域表达观测数据,可以仅关注部分高斯波包框架函数上的数据投影,这样既实现了波场的压缩存储,同时可利用高斯波包传播算子反传框架函数以实现整个炮道集的快速反传.这些综合了观测数据局部τ-p特征的高斯波包函数称为特征高斯波包（characteristic Gaussian packet,CGP）,相应的波场反传称为特征高斯波包反传.理论及数值分析证明了上述特征高斯波包反传方法是有效且快速的.炮点正传波场也利用高斯波包传播算子模拟.利用互相关成像条件可实现特征高斯波包叠前深度偏移（characteristic Gaussian packet pre-stack depth migration,CGPM）.由于高斯波包传播算子描述了局部方向及局部空间的波场,所以CGPM可以自然地提取角度域成像道集（ADCIG）,并易于实现面向目标叠前深度偏移,从而作为偏移引擎为偏移速度分析（MVA）服务.数值实验证明了CGPM和面向目标CGPM的有效性和实用性.     

双平方根方程三维叠前深度偏移

从双平方根（DSR）形式的波动方程出发，基于沉降观测概念和地震波扰动理论，介绍了深度域的DSR全偏移算子及共成像道集的生成方法. 根据三维地震数据的方位角特征，通过对全偏移算子的稳相近似，依次导出了适应于零方位角道集、Cross line共偏移距道集以及共偏移距矢量道集的偏移算子. 理论分析与合成数据的数值试验表明，DSR全偏移算子、共方位角偏移算子对介质速度变化的适应性很强，而其余两种偏移算子仅适用于缓变速情况.     

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

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

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

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

三维偏移距平面波有限差分叠前时间偏移

本文提出了中点-半偏移距域内的三维偏移距平面波(offset plane-wave)方程,并给出了其有限差分解法.偏移距平面波可通过对CMP道集进行平面波分解(倾斜叠加或线性Radon变换)生成,然而这样做会产生严重的噪音干扰.本文提出了局部倾斜叠加方法(local slant-stacking)来消除离散线性Rado...     

基于叠前方位振幅的大港—埕海地区奥陶系风化壳裂缝储层的叠前预测

P波通过裂缝地层传播时,其响应会随着地面观测方位的变化而变化,表现出明显的方位各向异性.这种随方位变化的裂缝响应信息(如振幅等)广泛存在于叠前方位一偏移距二维空间域中,检测此域中的裂缝响应,通过椭圆拟合的方法识别各向异性的存在,进而计算裂缝的方向——椭圆长轴的方向,和裂缝密度——椭圆长轴与短轴的比.本文研究出了一种基于常规叠前地震反射振幅来检测裂缝的三维叠前裂缝储层预测方法.它采用面元组合的方法扩大方位角和偏移距范围;采用多级数据控制技术获取所需的信息;采用稳妥的技术步骤来预测.通过大港油田实际工区的应用,结果表明本文研究出的方法技术应用效果好、精度高、经济实用.     

Handling dip discrimination phenomenon in common‐reflection‐surface stack via combination of output‐imaging‐scheme and migration/demigration

In the application of a conventional common‐reflection‐surface (CRS) stack, it is well‐known that only one optimum stacking operator is determined for each zero‐offset sample to be simulated. As a result, the conflicting dip situations are not taken into account and only the most prominent event contributes to any a particular stack sample. In this paper, we name this phenomenon caused by conflicting dip problems as ‘dip discrimination phenomenon’. This phenomenon is not welcome because it not only leads to the loss of weak reflections and tips of diffractions in the final zero‐offset‐CRS stacked section but also to a deteriorated quality in subsequent migration. The common‐reflection‐surface stack with the output imaging scheme (CRS‐OIS) is a novel technique to implement a CRS stack based on a unified Kirchhoff imaging approach. As far as dealing with conflicting dip problems is concerned, the CRS‐OIS is a better option than a conventional CRS stack. However, we think the CRS‐OIS can do more in this aspect. In this paper, we propose a workflow to handle the dip discrimination phenomenon based on a cascaded implementation of prestack time migration, CRS‐OIS and prestack time demigration. Firstly, a common offset prestack time migration is implemented. Then, a CRS‐OIS is applied to the time‐migrated common offset gather. Afterwards, a prestack time demigration is performed to reconstruct each unmigrated common offset gather with its reflections being greatly enhanced and diffractions being well preserved. Compared with existing techniques dealing with conflicting dip problems, the technique presented in this paper preserves most of the diffractions and accounts for reflections from all possible dips properly. More importantly, both the post‐stacked data set and prestacked data set can be of much better quality after the implementation of the presented scheme. It serves as a promising alternative to other techniques except that it cannot provide the typical CRS wavefield attributes. The numerical tests on a synthetic Marmousi data set and a real 2D marine data set demonstrated its effectiveness and robustness.     

Double-square-root one-way wave equation prestack tau migration in heterogeneous media

In this paper, source‐receiver migration based on the double‐square‐root one‐way wave equation is modified to operate in the two‐way vertical traveltime (τ) domain. This tau migration method includes reasonable treatment for media with lateral inhomogeneity. It is implemented by recursive wavefield extrapolation with a frequency‐wavenumber domain phase shift in a constant background medium, followed by a phase correction in the frequency‐space domain, which accommodates moderate lateral velocity variations. More advanced τ‐domain double‐square‐root wave propagators have been conceptually discussed in this paper for migration in media with stronger lateral velocity variations. To address the problems that the full 3D double‐square‐root equation prestack tau migration could meet in practical applications, we present a method for downward continuing common‐azimuth data, which is based on a stationary‐phase approximation of the full 3D migration operator in the theoretical frame of prestack tau migration of cross‐line constant offset data. Migrations of synthetic data sets show that our tau migration approach has good performance in strong contrast media. The real data example demonstrates that common‐azimuth prestack tau migration has improved the delineation of the geological structures and stratigraphic configurations in a complex fault area. Prestack tau migration has some inherent robust characteristics usually associated with prestack time migration. It follows a velocity‐independent anti‐aliasing criterion that generally leads to reduction of the computation cost for typical vertical velocity variations. Moreover, this τ‐domain source‐receiver migration method has features that could be of help to speed up the convergence of the velocity estimation.     

3D post‐stack interval velocity analysis with effective use of datuming

Interval velocity analysis using post‐stack data has always been a desire, mainly for 3D data sets. In this study we present a method that uses the unique characteristics of migrated diffractions to enable interval velocity analysis from three‐dimensional zero‐offset time data. The idea is to perform a standard three‐dimensional prestack depth migration on stack cubes and generate three‐dimensional common image gathers that show great sensitivity to velocity errors. An efficient ‘top‐down’ scheme for updating the velocity is used to build the model. The effectiveness of the method is related to the incorporation of wave equation based post‐stack datuming in the model building process. The proposed method relies on the ability to identify diffractions along redatumed zero‐offset data and to analyse their flatness in the migrated local angle domain. The method can be considered as an additional tool for a complete, prestack depth migration based interval velocity analysis.     

用共方位角偏移实现叠前二维资料三维化成像

波动方程偏移保持了波场动力学特征,依此本文应用共方位角叠前深度偏移技术来实现二维资料三维化处理.本方法能使共反射点偏移成像与叠前数据插值一步完成,较好地解决二维资料处理中的成像点不能准确归位、侧面反射波难处理、主测线与联络线不闭合等问题.