三维VSP数据高效偏移成像的超道集方法

当前的三维VSP地震数据偏移成像实现都是在共炮点道集或共检波点道集中逐个道集循环进行的,计算效率相对较低.根据三维VSP观测系统中炮点和检波点布置的特殊性和地震波场满足线性叠加的特性,本文提出了一种三维VSP数据的高效偏移成像方法,即首先通过对三维VSP共接收点道集进行地震数据的广义合成得到一种超道集,然后在共接收点道集的波场深度外推过程中逐步应用多震源波场对超道集进行偏移成像,即利用一次波场深度外推循环完成对所有共检波点道集数据的偏移成像.通过三维VSP模型数据与实际地震数据的偏移成像试验验证了这种高效的超道集偏移成像方法可取得与常规共检波点道集相当的偏移成像效果,还具有极高的计算效率,其计算量与单个共检波点道集的偏移成像计算量相当.     

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

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

共反射角叠前偏移成像研究及应用

共偏移距道集已被广泛地应用于地震速度建模及振幅随偏移距变化(AVO)的研究中,但复杂构造及射线多路径产生的共偏移距道集不保幅性等一系列缺陷给AVO研究带来很大的困难.共反射角道集包含有能反映地下速度和岩性变化的信息,更有利于速度模型优化、地震振幅属性分析及地下岩性和断裂的研究.本文通过研究共反射角深度偏移方法和理论,完善了基于目标的共反射角深度偏移技术,提出了获得相对保幅共反射角道集方法.该方法克服了共偏移距域道集在复杂介质中遇到的困难,更能有效地反映波场和地质结构方面的信息.通过理论模型数据进行了试算,并采用实际地震数据对此方法进行了验证,在陡倾角成像方面取得较好效果.     

广义地震数据合成及其偏移成像

根据地震波场的线性叠加原理,提出了对地震共炮道集及其震源进行线性叠加的一般方案——广义地震数据合成的方法.利用这个方法,可以根据不同的地质情况和要求得到各种不同的人工合成地震数据道集和震源,如平面波数据道集和震源、局部平面波（束）数据道集和震源以及面向目标的人工合成地震数据道集和震源.对于人工合成地震数据道集的偏移成像可应用单平方根方程实现.不同的合成地震数据道集具有不同偏移成像特性：平面波数据道集具有很高的计算效率,局部平面波数据道集具有很好的方向性,面向目标的合成地震数据道集具有很好的面向目标特性.     

面炮成像、控制照明与AVA道集

基于波场延拓的叠前深度偏移是实现复杂构造地质体成像的最可靠方法,但存在着计算量大、对观测系统适应性差等缺点.面炮偏移是波动方程实现精确叠前成像的另一类方法,具有较高的计算效率,不存在偏移孔径问题,而且可以通过控制照明方法,解决平面波在目标区域的能量补偿问题.本文采用面炮成像技术进行叠前深度偏移,通过对面炮震源下行波场的质量控制和优选射线个数和范围,以达到最佳的成像效果.采用控制照明技术,较大地提高了目标地层的成像精度.与此同时,得到振幅随入射角变化(AVA)道集,有利于叠前振幅解释和储层岩性预测.数据实验表明面炮成像技术是一种快速有效的方法,其成像精度与单平方根算子的共炮点道集偏移和双平方根算子的共中心点道集偏移相当,但在计算速度上要快得多,而且易于并行计算.     

多波高斯束叠前深度偏移速度分析(英文)

多分量地震资料叠前深度偏移技术可以对地下复杂地质构造进行更准确的成像,精确成像的前提是获取准确的纵横波偏移速度。本文采用高斯束偏移方法对多波地震数据进行偏移速度分析,首先分别给出纵波和转换波共偏移距域高斯束叠前深度偏移方法原理,在此基础上抽取纵波和转换波偏移距域共成像点道集;然后根据共成像点道集拉平准则,分别对纵波和横波速度进行更新;当两种波成像深度不一致时,对纵波和转换波成像剖面进行深度匹配,完成高精度的纵横波偏移速度分析。模型数据和实际资料试算表明,该方法是一种有效的多波偏移速度分析方法。     

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

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

基于L1-L1范数稀疏表示的共偏移距道集地震数据重建方法

传统地震数据稀疏重建方法面临着:(1)叠前共炮点道集或CMP道集反射波为双曲线型同相轴,地震数据重建会损害有效波;(2)地震信号存在噪声和畸变,要求重建方法具有较好的噪声鲁棒性.针对这两个问题,提出一种基于L_1-L_1范数稀疏表示的共偏移距道集地震数据重建方法.该方法利用了共偏移距道集中地震波为水平同相轴,无道间时差,满足空间重建要求,和L_1-L_1范数稀疏表示具有较好的噪声鲁棒性.首先抽取共偏移距道集地震数据,并根据地震采集信息构造复合采样矩阵,然后采用L_1-L_1范数稀疏表示对数据稀疏重建后,再将数据反变换回共炮点道集或CMP道集,能够同时实现地震信号稀疏重建和随机噪声压制.理论模型和实际数据试算结果验证所提方法具有较好重建精度和噪声鲁棒性.     

Kirchhoff积分叠前时间偏移全方位角度道集生成方法与实现

方位角度域共成像点道集能够客观反映地下介质的速度、各向异性参数异常以及振幅随角度变化（AVA）和裂缝信息。传统Kirchhoff PSTM通常输出偏移距域共成像点道集,对于速度分析、各向异性分析、AVA分析、裂缝识别等均存在诸多不便。本文提出了基于走时梯度的Kirchhoff叠前时间偏移全方位角度集输出方法并提出工业上切实可行的实现方案。通过走时场梯度计算波场传播方向矢量,形成能够反映观测系统参数和波场传播情况的全方位角度域共成像点道集。为了在大规模地震数据Kirchhoff积分叠前时间偏移中输出全方位角度道集,本文给出基于输入道方式的偏移实现方法,采用逐条inline线进行线偏移成像,从而大大降低了全方位角度道集输出对计算机内存的压力,显著提高了Kirchhoff积分时间偏移输出全方位角度道集的可行性。三维盐丘模型测试和海上某区块三维实际资料试验证明了本文方法的正确性。      

角度域共成像点道集剩余曲率法偏移速度分析

叠前深度偏移速度分析是地震数据处理方法研究的重点,是影响地震成像效果的关键技术之一。本文的基于角度域共成像点道集剩余曲率的偏移速度分析方法,旨在通过角度域共成像点道集(ADCIGs)来建立局部入射角与偏移深度之间的联系,通过角度道集剩余速度分析来修正、更新偏移速度。不同于偏移距域共成像点道集(ODCIGs),ADCIGs较少地受到由于多路径问题而产生的波场干涉假象的影响,更适合于偏移速度分析。本文在推导出基于ADCIGs的水平地层深度剩余量公式的基础上,提出了基于ADCIGs的沿层扫描法速度更新流程,凹陷模型和实际资料的试处理表明了本方法的正确性和有效性。     

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.     

用于AVO分析的振幅保真平面波叠前时间偏移

To support amplitude variation with offset （AVO） analysis in complex structure areas, we introduce an amplitude-preserving plane-wave prestack time migration approach based on the double-square-root wave equation in media with little lateral velocity variation. In its implementation, a data mapping algorithm is used to obtain offset-plane-wave data sets from the common-midpoint gathers followed by a non-recursive phase-shift solution with amplitude correction to generate common-image gathers in offset-ray-parameter domain and a structural image. Theoretical model tests and a real data example show that our prestack time migration approach is helpful for AVO analysis in complex geological environments.     

基于Gabor-Daubechies小波束叠前深度偏移的角度域共成像道集

传统炮检距域共像集（CIG）在复杂介质中因波传播的多路径而存在反射体位置不确定的问题. 角度域CIG由于克服了这一缺陷而逐步成为速度分析、AVA以及振幅保真偏移成像等研究的主要手段. 以波动理论为基础的地震偏移成像方法的发展为获得高质量的角度域CIG提供了可靠的实现途径. 其中,基于波场局域化分解和传播的小波束域波场延拓和偏移成像方法,因其波场分解基本函数和传播算子在空间和方向上的双重局域特性,而成为角度相关分析研究的有效工具. 本文在采用Gabor Daubechies框架分解的小波束叠前角度域偏移成像基础上,利用不同的叠加方法由局部角度域像矩阵得到了反射角域CIG（CRAIG）和倾角域CIG（CDAIG）. 以SEG EAGE二维盐体模型为例,通过对CRAIG和CDAIG的对比,探讨了这两种角度域CIG的特点及其在地震偏移成像中的潜在应用.     

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

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

面炮叠前深度偏移与控制照明技术

基于波场延拓的叠前深度偏移是实现复杂构造地质体成像的最可靠方法,但存在着计算量大、对观测系统适应性差等缺点。面炮偏移是波动方程实现精确叠前成像的另一类方法,具有较高的计算效率,不存在偏移孔径问题,而且可以通过控制照明方法,解决平面波在目标区域的能量补偿问题。本文采用面炮成像技术进行叠前深度偏移成像,通过对面炮震源下行波场的质量控制以及射线参数的个数和范围的选取,以达到最佳的成像效果。采用不同深度点上的控制照明技术,较大地提高了目标地层的成像精度。数据实验表明面炮成像技术是一种快速有效的方法,其成像精度与单平方根算子的共炮点道集偏移和双平方根算子的共中心点道集偏移相当,但在计算速度上要快得多,而且易于并行计算。     

频率-空间域有限差分法叠前深度偏移

为了处理横向强变速介质中的深度成像问题,本文提出一种基于共炮道集的优化系数的傍轴近似方程叠前深度偏移算子,并在基于反射系数估算的成像条件下,可实现叠前深度偏移成像.该算子具有方程阶数低且能对陡倾角成像的特征,并采用有限差分法波场延拓,能适应速度场的任意变化.当在频率-空间域进行计算时,相对于纯粹的时间-空间域有限差分算法有计算效率高、成像方便的优点.脉冲响应测试和对Marmousi模型进行的叠前深度偏移结果表明,该偏移方法在强横向变速情况下具有非常好的成像效果.     

Double plane‐wave reverse‐time migration

We develop a new time‐domain reverse‐time migration method called double plane‐wave reverse‐time migration that uses plane‐wave transformed gathers. Original shot gathers with appropriate data acquisition geometry are double slant stacked into the double plane‐wave domain with minimal slant stacking artefacts. The range of plane‐wave components needed for migration can be determined by estimating the maximum time dips present in shot gathers. This reduces the total number of input traces for migration and increases migration efficiency. Unlike the pre‐stack shot‐profile reverse‐time migration where the number of forward propagations is proportional to the number of shots, the number of forward propagations needed for the proposed method remains constant and is relatively small even for large seismic datasets. Therefore, the proposed method can improve the efficiency of the migration and be suitable for migrating large datasets. Double plane‐wave reverse‐time migration can be performed for selected plane‐wave components to obtain subsurface interfaces with different dips, which makes the migration method target oriented. This feature also makes the method a useful tool for migration velocity analysis. For example, we are able to promptly obtain trial images with nearly horizontal interfaces and adjust velocity models according to common image gathers. Seismic signal coming from steeply dipping interfaces can be included into the migration to build images with more detailed structures and higher spatial resolution as better velocity models become available. Illumination compensation imaging conditions for the proposed method are also introduced to obtain images with balanced amplitudes.     

波动方程偏移速度建模:一种直接反演方法

基于炮域波动方程叠前深度偏移所产出的角道集,发展了不需要多次迭代的波动方程偏移速度建模方法.文中分析了炮域波动方程偏移生成角道集的方法,给出了非均匀介质中角道集同相轴曲率与速度误差的定量关系,发展了基于同相轴曲率的速度模型直接更新算法.在初始模型较合理的情况下,应用一次炮域波动方程偏移计算,即可得到较准确的速度模型.这...