共炮检距道集波动方程保幅叠前深度偏移方法

本文提出了一种基于双平方根算子的共炮检距道集波动方程保幅叠前深度偏移方法,将振幅误差补偿作为偏移的一部分与“运动学偏移”一起在偏移过程中实现.其基本内容包括：(1)从保幅的单平方根算子方程出发,推导出由双平方根算子定义的保幅单程波方程;(2)根据地震波摄动理论把速度场分裂为层内常速背景和变速扰动,分别在频率－波数域和频率－空间域求得波场深度延拓的偏移时移量及振幅校正系数,从而得到最终的DSR保幅波场延拓算子;(3)在高频假设条件下,把DSR保幅波场延拓公式中的积分运算进行稳相近似,得到保幅波场延拓的相移公式.理论分析和模型数值试验表明,该方法不但可以使散射能量聚焦、归位,提高成像精度;而且可以输出正确反映地下反射系数的振幅信息,为后续的地震属性分析（如AVO/AVA）提供更真实的地震信息.     

大步长波场深度延拓的理论

波场延拓是地震偏移成像的基础. 快速进行目标区波场延拓对石油勘探中急需发展的深部地震勘探和无组合海量地震数据的成像有重要意义. 在目标区成像中，目前已有的波场延拓方法，包括基于走时计算的Dix方法和射线追踪方法，以及基于小步长波场递推的方法，在适应复杂介质、计算精度和计算效率的某一方面还不能完全满足实际需要. 本文提出一种基于“算子相位”李代数积分的快速计算延拓算子的方法，称为大步长波场延拓方法. 在该方法中，指向目标区的波场延拓算子象征的复相位被表示成波数的线性组合. 线性组合的系数是层速度函数及其导数的深度积分，计算和存储较为方便. 波场延拓算子通过相移算子加校正的方法，利用快速Fourier变换在空间域和波数域予以实现. 利用动力学等价关系导出了便于计算的表达式. 本文比较了算子主象征函数用一步法展开和用两步法展开的精度，从而说明大步长方法的精度要高于递推方法. 在横向和纵向线性变化介质中，将大步长方法的脉冲响应与递推法做了比较，说明大步长延拓算子的走时精度主要取决于相移因子中的横向变速校正项；且在各种近似下，大步长算子发生的频散都非常小.     

VTI介质qP波方程高精度有限差分算子

波动方程有限差分法是一种使用广泛的地震波数值模拟方法.但是有限差分法本身固有存在着数值频散问题，会降低地震波场模拟的精度与分辨率.为了克服常规有限差分算子的数值频散，本文针对VTI介质地震波数值模拟问题，构造了频率-空间域qP波波动方程高精度有限差分优化算子，根据最优化理论中高斯-牛顿法确定了高精度有限差分算子的优化系数.利用常规差分算子和高精度优化差分算子对归一化相速度的频散关系精度进行了对比分析，并对均匀各向同性介质和均匀VTI介质中的qP波地震波场进行了有限差分数值模拟，通过频散关系精度分析和波场数值模拟结果表明：有限差分优化算子具有较高的波场数值模拟精度，有效压制了传统有限差分算子数值模拟中的数值频散现象，提高了有限差分算子精度，为VTI介质频率-空间域qP波正演模拟奠定了基础.     

大步长单程波算子相位积分实现方式初探

大步长单程波算子是穿过厚层的单程波算子的积分,在理论上它可以通过频率波数域表达式和波数域到空间域变换来实现,目前,通过李代数积分和指数映射的研究,已经导出频率波数域表达式,需要研究波数域到空间域变换的实现方式,波数域到空间域的变换可以通过鞍点法和相位积分来实现,本文研究了相位积分的三种实现途径:分裂方法、相移加插值方法和佳格点方法,在相移加插值方法中提出一种波场"角度域插值"的插值方法.     

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

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

最优可分表示法三维波场延拓

适于复杂介质的高精度波场延拓算子是叠前深度偏移研究的重要内容。本文采用最优可分表示方法,运用正反傅立叶变换构造了三维单程波场延拓算子,算子实现了波数域变量与空间(速度)域变量分离。波数域内进行相移计算,在空间域对因介质横向变速引起的时移作修正。脉冲响应显示在区域内各速度的脉冲计算值与理论值基本一致,说明最优可分表示法叠前深度偏移可适用于强变速条件下复杂介质的成像需求。SEG／EAGE模型和实测数据的成像结果验证了本文方法对复杂构造的成像能力。     

单程波算子积分解的象征表示

单程波波场延拓算子在地震偏移成像中有重要应用.单程波波场延拓算子按其实现方式可分为Kirchhoff积分、空间隐式有限差分和Fourier变换方法，他们代表了算子的不同表示方法，当截断使用这些方法时会得到不同的精度.象征表示对这些方法的导出和精度分析有重要作用.算子作用于正弦波函数所得函数称为算子的象征.算子的象征是褶积算子Fourier变换的推广.Fourier变换方法则直接用象征函数的可分表示求出.空间隐式有限差分则可以用象征函数的Padè近似或部分分式导出.单程波算子在深度域的积分称为单程波算子积分解.本文推导了单程波算子积分解的象征表达式，给出了算子象征的代数运算的头几阶表达式，这些表达式还未在前人文献中发现.Kirchhoff积分所需格林函数可以通过象征函数和鞍点法导出.基于积分解的象征表达式给出了非对称走时公式，对改善Kirchhoff积分的聚焦性能有重要意义.     

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

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

轴对称各向异性介质波动方程深度偏移的对称非平稳相移方法

由所建立的三维qP波相速度表示式出发，导出并解析求解各向异性介质中的频散方程，得到三维各向异性介质中的相移算子，进而将以相移算子为基础的对称非平稳相移方法推广到各向异性介质，发展了一个三维各向异性介质的深度偏移方法. 文中使用的各向异性介质的速度模型与现行的各向异性构造的速度估计方法一致，将各向同性、弱各向异性及强各向异性统一在一个模型中. 所建立的各向异性介质对称非平稳相移波场延拓算子可以同时适应速度及各向异性参数横向变化；文中给出的算例虽然是针对二维VTI介质的，但所提出的算法同样适用于三维TI介质.     

三维复杂构造中地震波模拟的单程波方法

复杂构造中单程波与双程波方法模拟结果的比较表明,就地震勘探中主要关心的一次反射波而言,单程波算法已具有足够的精度. 使用单程波方程将极大地减少数值计算的计算量,同时对介质的几何和物理参数建模也降低了要求. 单程波算法可视为深度偏移的“逆运算”,这样可以很好地借用已知的深度偏移方法及其程序系统. 基于计算效率和计算精度的双重考虑,本文在介质速度结构较复杂时采用显式短算子波场延拓方法,而在介质速度结构相对简单时采用分裂步相移法. 反射系数的计算中考虑了其随入射角的变化.     

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

基于三维VTI各向异性介质的频散关系,构建波数项和空间项分离的单程波算子表达式,以优化算法,确定算子的待定系数,实现广角逼近三维VTI介质的广义相移算子,发展了可灵活处理强或弱各向异性介质的波动方程叠前深度偏移方法.文中同时也针对其工业应用建议了三维VTI各向异性介质中可提高计算效率的频率相关变步长波场深度延拓算法及稀...     

3D TABLE-DRIVEN MIGRATION1

An efficient full 3D wavefield extrapolation technique is presented. The method can be used for any type of subsurface structure and the degree of accuracy and dip-angle performance are user-defined. The extrapolation is performed in the space-frequency domain as a space-dependent spatial convolution with recursive Kirchhoff extrapolation operators. To get a high level of efficiency the operators are optimized such that they have the smallest possible size for a specified accuracy and dip-angle performance. As both accuracy and maximum dip-angle are input parameters for the operator calculation, the method offers the possibility of a trade-off between these quantities and efficiency. The operators are calculated in advance and stored in a table for a range of wavenumbers. Once they have been calculated they can be used many times. At the basis of the operator design is the well-known phase-shift operator. Although this operator is exact for homogeneous media only, it is assumed that it may be applied locally in case of inhomogeneities. Lateral velocity variations can then be handled by choosing the extrapolation operator according to the local value of the velocity. Optionally the operators can be designed such that they act as spatially variant high-cut filters. This means that the evanescent field can be suppressed in one pass with the extrapolation. The extrapolation method can be used both in prestack and post-stack applications. In this paper we use it in zero-offset migration. Tests on 2D and 3D synthetic and 2D real data show the excellent quality of the method. The full 3D result is much better then the result of two-pass migration, which has been applied to the same data. The implementation yields a code that is fully vectorizable, which makes the method very suitable for vector computers.     

Cost-effective 3D one-pass depth migration1

A crucial point in the processing of 3D seismic data is the migration step, both because of its 3D nature and the computational cost involved. The efficiency and accuracy of 3D migration are determined by the wavefield extrapolation technique employed. Wavefield extrapolation based on second-order differential operators of variable-length is very efficient and accurate at the same time. Compared to migration based on the McClellan transform and operator splitting, the use of variable-length second-order differential operators offers significant advantages. The 3D migration operator has an almost perfect circular symmetry. No positioning errors in the 45° azimuth between the in-line and cross-line directions are evident. The method is, in practice, only limited by spatial aliasing and does not require expensive interpolation of data to reduce numerical artifacts. This reduces the computational cost of 3D one-pass depth migration by a large factor.     

A self-adaptive algorithm for choosing reference velocities in the presence of lateral velocity variations

Based on perturbation theory, the wave equation extrapolation operator with mixed domains has the ability to deal with lateral velocity variations. It is the image method that has undergone much research in seismology. All extrapolation operators face the problem of choosing the reference velocity due to continuation in depth. The wavefield extrapolation operator with a single reference velocity is suitable for media with weak lateral variation. The multi-reference velocity extrapolation operator can cope with severe lateral velocity variations and improve image accuracy. However, the calculation cost is large. We present a self-adaptive approach to automatically determine the number of selected reference velocities according to the complexity of structure and the given velocity threshold value. The approach can be used to construct the SSF, FFD, WXFD, and GSP multi-reference velocity wavefield extrapolation image algorithms. The result of a salt-dome model data test demonstrates that the self-adoptive multi-reference wavefield extrapolation algorithm has the ability to deal with severe lateral velocity variations and can also be used for structure edge detection. The method is flexible and computationally cost-effective.     

自适应横向变速参考速度选择方法

基于扰动理论的混合域波动方程波场外推算子,具有一定介质横向速度变化适应能力,是反射地震学中研究较多的成像方法。此类波场外推算子沿深度层进行波场外推,都存在参考速度选择问题。单参考速度波场外推算子,适应地下介质弱横向变速,而多参考速度波场外推算子可以提高横向变速的适应能力和成像精度,但要以大量计算为代价。本文提出的自适应多参考速度选择策略,根据外推层地质构造的复杂度和给定的速度门槛值自动选择参考速度个数,利用该策略构造混合域SSF、FFD、WXFD和GSP等多参考速度波场外推成像算法。盐丘模型理论数据测试结果表明,自适应多参考速度波场外推算法具有强横向变速适应能力和较高成像精度。     

波动方程深度偏移的频率相关变步长延拓方法

发展了波动方程深度延拓的频率相关变步长深度延拓方法和表驱动的单点波场插值技术.前者通过减少深度延拓的次数减少了波动方程深度偏移的计算量,而后者用很少的计算量实现了等间距、理想采样的深度成像.就同一偏移方法,采用频率相关变步长深度延拓加单点插值,其计算量大约是常规的等间距采样延拓方法的三分之一,但两者的成像效果基本相同.文中以最优分裂Fourier方法为例,用二维理论数据（Marmousi模型）和三维实际地震资料验证了这一方法,但这一方法可适用于各类频率域波动方程深度偏移方法.     

Depth imaging in anisotropic media by symmetric non-stationary phase shift

We present a new depth‐imaging method for seismic data in heterogeneous anisotropic media. This recursive explicit method uses a non‐stationary extrapolation operator to allow lateral velocity variation, and it uses the relationship between phase angle and the spectral coordinates of seismic data to allow velocity variation with phase angle. A qualitative comparison of migration impulse responses suggests that, for an equivalent cost, the symmetric non‐stationary phase‐shift (SNPS) operator is superior to the phase‐shift plus interpolation (PSPI) operator, for very large depth intervals. To demonstrate the potential of the new method, seismic data from a physical model acquired over a transversely isotropic medium are imaged using a shot‐record migration based on the SNPS operator.     

TOWARDS OPTIMUM ONE-WAY WAVE PROPAGATION1

Numerical wavefield extrapolation represents the backbone of any algorithm for depth migration pre- or post-stack. For such depth imaging techniques to yield reliable and interpretable results, the underlying wavefield extrapolation algorithm must propagate the waves through inhomogeneous media with a minimum of numerically induced distortion, over a range of frequencies and angles of propagation. A review of finite-difference (FD) approximations to the acoustic one-way wave equation in the space-frequency domain is presented. A straightforward generalization of the conventional FD formulation leads to an algorithm where the wavefield is continued downwards with space-variant symmetric convolutional operators. The operators can be precomputed and made accessible in tables such that the ratio between the temporal frequency and the local velocity is used to determine the correct operator at each grid point during the downward continuation. Convolutional operators are designed to fit the desired dispersion relation over a range of frequencies and angles of propagation such that the resulting numerical distortion is minimized. The optimization is constrained to ensure that evanescent energy and waves propagating at angles higher than the maximum design angle are attenuated in each extrapolation step. The resulting operators may be viewed as optimally truncated and bandlimited spatial versions of the familiar phase shift operator. They are unconditionally stable and can be applied explicitly. This results in a simple wave propagation algorithm, eminently suited for implementation on pipelined computers and on large parallel computing systems.     

Wavefield extrapolation and prestack depth migration in anelastic inhomogeneous media

Wavefield depth extrapolation and prestack depth migration in complex anelastic media are studied. Kjartansson's frequency‐independent Q law is used to describe the absorption of seismic energy. The macromodel used is analogous to the macromodel used for current migration schemes except that an additional frequency‐independent Q macromodel needs to be provided. Absorption in the forward one‐way propagator is introduced by assuming a complex phase velocity, and the inverse one‐way propagator is obtained using the reciprocity theorem for one‐way wavefields in dissipative media. The stability of the inverse propagator is achieved by limiting the angle of propagation of wavefields. A table‐driven explicit operator scheme for imaging complex 2D anelastic media is presented. High‐accuracy, short convolution operators are designed by the weighted least‐squares method, and two kinds of imaging conditions are proposed. Numerical examples of depth extrapolation in laterally varying media, the migration of a spatial impulse with dispersion as well as shot record depth migration demonstrate the potential of the proposed explicit forward operator, the explicit inverse operator and the prestack depth migration scheme, respectively.     

黏弹性VTI介质频率空间域准P波正演模拟

有限差分方法是波场数值模拟的一个重要方法,时间域有限差分计算方法因按时间片递推计算,每个时间片的舍入误差会累积到下一片中,当时间片较多,最终会导致累积误差太大.而频率域计算是按频率片对空间网格进行整体求解方程组,其计算误差分配到了每个网格点上,并且各个频率片之间是独立计算的,因此不存在累计误差,而且在频率-空间域更易于...