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

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

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

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

基于双线性插值的三维横向各向同性介质初至波射线追踪

给出了基于旅行时双线性插值(LTI)的三维横向各向同性(VTI)介质中的地震初至波旅行时计算及射线追踪方法,通过三维均匀及水平层状VTI介质模型的试算,证实了该方法的准确性及适应性.该方法可用于三维VTI介质中的深度偏移及层析成像.     

非稳态相移法叠前深度偏移

介绍一种能够适应介质速度横向变化的非稳态相移算子及其叠前深度偏移方法．为了克服常规相移偏移算法中要求速度横向不变的缺点,出现了基于非稳态滤波器理论的非稳态相移算子,即PSPI算子、NSPS算子和SNPS算子,其中SNPS算子是将前二者结合起来的一种对称的非稳态相移算子,它比前二者具有更高的精度和稳定性．为了提高运算速度,基于非稳态相移算子的叠前深度偏移算法采取了分片均匀近似的策略,Marmousi模型的叠前深度偏移结果证明了该算法的可行性和有效性。     

基于Walkaway VSP的Thomsen各向异性参数估计与应用(英文)

Thomsen各向异性参数的求取对于正确的时深转换和深度域偏移成像处理至关重要。相比其它各向异性参数估算方法,从VSP资料中更容易获取准确的各向异性参数用于地面地震偏移成像。本文分析研究了利用Walkaway VSP资料估算VTI介质Thomsen各向异性参数的方法,该方法基于VTI介质近偏移距动校正公式利用Walkaway VSP近偏移距初至信息求取各向异性参数δ;基于各向异性介质纵波速度Thomsen近似公式采用射线追踪时差扫描方法求取各向异性参数ε。数值模型正演表明利用该方法估算的各向异性参数误差较小。利用塔里木盆地8个方位的Walkaway VSP实际资料求取了该区深度域Thomsen各向异性参数ε和δ值,同时结合地面三维地震资料建立了较为准确的各向异性深度一速度模型用于叠前深度偏移成像,进一步提高了碳酸盐岩储层的成像精度,减小了目标地质体的深度误差。     

一个近似的VTI介质声波方程

声学近似是提高各向异性介质准纵波数值模拟和偏移处理计算速度的有效方法之一.为了获得VTI(具有垂直对称轴的横向各向同性)介质声波方程,根据VTI介质的精确相速度表达式,用待定系数法简化其中的根号项,得到一个近似相速度公式,由此公式通过反Fourier变换推导出一个时间二阶、空间四阶偏微分方程.定量相速度计算表明,当Th...     

三维各向异性TI介质中的P/S波快速解耦技术及在弹性波逆时偏移中的应用

随着多分量采集技术的发展,弹性波逆时偏移技术在三维各向异性介质复杂地质构造成像中得到了广泛的应用.然而耦合的P波场和S波场,会在传播过程中产生串扰噪声,降低弹性波逆时偏移的成像精度.为了解决这一问题,本研究针对具有倾斜各向异性对称轴的三维横向各向同性(Transverse Isotropy, TI)介质,提出了一种矢量弹性波场快速解耦方法,可以有效提高偏移剖面的成像质量.该方法首先通过坐标转换,将观测系统坐标系的垂直轴旋转到TI介质的对称轴方向,在新坐标系下,根据具有垂直对称轴的三维横向各向同性(Vertical Transverse Isotropy, VTI)介质中的分解算子,推导出三维TI介质解耦算子表达式.接着引入一种在空间域快速计算分解波场的方法,来实现空间域矢量P波场和S波场分离,极大地提高了计算效率.最后,通过点积成像条件,将提出的P/S波分解方法引入到三维TI介质弹性波逆时偏移中,得到高精度的PP和PS成像.与以往的波场分解方法相比,本文方法具有数值稳定和计算效率高的特点.数值算例表明,应用上述三维TI分解算子得到的偏移剖面有效压制了噪声,提高了成像质量.     

VTI介质角度域转换波高斯束偏移成像方法研究

转换波偏移可以利用纵横波波场信息,得到高分辨率的成像结果,从而为油藏描述提供高质量的地震资料.目前的研究主要是利用纵波波场信息进行偏移成像,然而,传统的纵波方法在复杂探区成像时具有一定的局限性.为此,本文在各向异性介质声波射线追踪算法的基础上,推导出各向异性介质转换波射线追踪方程,发展了一种转换波射线追踪算法;并将研究的追踪算法应用到偏移成像中,提出了一种各向异性VTI介质角度域转换波高斯束偏移成像方法.通过各向异性VTI介质断块模型和复杂构造模型试算,说明了本文方法的正确性和有效性.模型试算的结果表明,在考虑地下各向异性时,本文研究的方法具有更好的成像效果,提取的角道集结果可以为偏移速度分析提供依据.     

三维VTI介质qP波方程频率空间域有限差分高阶加权算子

波动方程有限差分法是波场模拟的一个重要方法,为解决常规有限差分法存在着数值频散的问题,本文从具有垂直对称轴的三维横向各向同性(VTI)介质频率-空间域qP波动方程出发,在常规差分算子的基础上构造了适合三维VTI介质的频率空间域有限差分优化算子,然后利用最优化理论中的Gauss-Newton法求解了优化算子的系数,使差分方程的相速度与波动方程的相速度尽量吻合,从而在理论上使网格数值频散达到极小,精度对比分析及数值测试表明,有限差分优化算子具有较高的波场数值模拟精度,有效地压制了数值频散现象,为三维VTI介质频率一空间域qP波正演模拟研究提供了理论基础.     

TTI介质隐式有限差分平面波偏移

文章研究了TTI介质隐式有限差分(IFD)波场外推算子和TTI介质平面波偏移.与各向同性和VTI介质相比,TTI介质频散关系要复杂的多,很难得到频散关系的显性表达,也不能再只用对称的偶函数表示,需要加入一个奇函数.文中设计TTI介质IFD波场外推算子并用非线性优化方法求解算子系数.将各向同性介质下平面波偏移理论与TTl...     

Kirchhoff prestack depth migration in 3-D simple models: comparison of triclinic anisotropy with simpler anisotropies

We use Kirchhoff prestack depth migration to calculate migrated sections in 3-D simple anisotropic homogeneous velocity models in order to demonstrate the impact of anisotropy on migrated images. The recorded wave field is generated in models composed of two homogeneous layers separated by one either non-inclined or inclined curved interface. The anisotropy in the upper layer is triclinic. We apply Kirchhoff prestack depth migration to velocity models with different types of anisotropy: a triclinic anisotropic medium, an isotropic medium, transversely isotropic media with a horizontal (HTI) and vertical (VTI) symmetry axis. We observe asymmetry in migration caused by triclinic anisotropy and we show the errors of the migrated interface caused by inaccurate velocity models used for migration. The study is limited to P-waves.     

Modeling and RTM for arbitrary-order pure acoustic wave equation in VTI media using normalized pseudo-analytical method

The conventional pseudo-acoustic wave equations(PWEs) in vertical transversely isotropic(VTI)media may generate SV-wave artifacts and propagation instabilities when anisotropy parameters cannot satisfy the pseudo-acoustic assumption. One solution to these issues is to use pure acoustic anisotropic wave equations, which can produce stable and pure P-wave responses without any SVwave pollutions. The commonly used pure acoustic wave equations(PAWEs) in VTI media are mainly derived from the decoupled P-SV dispersion relation based on first-order Taylor-series expansion(TE), thus they will suffer from accuracy loss in strongly anisotropic media. In this paper, we adopt arbitrary-order TE to expand the square root term in Alkhalifah's accurate acoustic VTI dispersion relation and solve the corresponding PAWE using the normalized pseudoanalytical method(NPAM) based on optimized pseudodifferential operator. Our analysis of phase velocity errors indicates that the accuracy of our new expression is perfectly acceptable for majority anisotropy parameters. The effectiveness of our proposed scheme also can be demonstrated by several numerical examples and reverse-time migration(RTM) result.     

基于平面波加速的VTI介质最小二乘逆时偏移

地震各向异性集中表现为速度各向异性,势必影响地震波运动学特征.传统声波逆时偏移(RTM)和最小二乘逆时偏移(LSRTM)没有考虑介质各向异性特征,导致反射波不能正确归位、同相轴出现扭曲及寻优速度慢或不收敛等,VTI介质逆时偏移(VTI-RTM)矫正了声波成像的不足,但仍存在低频干扰严重、中深部成像不佳、振幅保持差等缺陷.为此,本文首先实现了VTI介质最小二乘逆时偏移(VTI-LSRTM)方法,为了节省I/O及内存需求并提高效率,进一步引入平面波编码技术,提出了一种基于平面波加速的VTI介质最小二乘逆时偏移(VTI-PLSRTM)策略.在此基础上开展了简单模型及复杂Marmousi模型成像试验,并与标准逆时偏移剖面对比表明:本方法能够校正各向异性造成的相位畸变,且在迭代中自动压制串扰及低频噪声、补偿中深部能量,是一种兼具质量与效率的保幅成像策略;对速度误差的敏感性测试说明该方法需要相对正确的偏移速度及Thomsen参数模型.     

Anisotropic migration velocity analysis: Application to a data set from West Africa

Although it is widely recognized that anisotropy can have a significant influence on the focusing and positioning of migrated reflection events, conventional depth imaging methods still operate with isotropic velocity fields. Here, we present an application of a 2D migration velocity analysis (MVA) algorithm, designed for factorized v(x, z) VTI (transversely isotropic with a vertical symmetry axis) media, to an offshore data set from West Africa. By approximating the subsurface with factorized VTI blocks, it is possible to decouple the spatial variations in the vertical velocity from the anisotropic parameters with minimal a priori information. Since our method accounts for lateral velocity variation, it produces more accurate estimates of the anisotropic parameters than those previously obtained with time‐domain techniques. The values of the anellipticity parameter η found for the massive shales exceed 0.2, which confirms that ignoring anisotropy in the study area can lead to substantial imaging distortions, such as mis‐stacking and mispositioning of dipping events. While some of these distortions can be removed by using anisotropic time processing, further marked improvement in image quality is achieved by prestack depth migration with the estimated factorized VTI model. In particular, many fault planes, including antithetic faults in the shallow part of the section, are better focused by the anisotropic depth‐migration algorithm and appear more continuous. Anisotropic depth migration facilitates structural interpretation by eliminating false dips at the bottom of the section and improving the images of a number of gently dipping features. One of the main difficulties in anisotropic MVA is the need to use a priori information for constraining the vertical velocity. In this case study, we successfully reconstructed the time–depth curve from reflection data by assuming that the vertical velocity is a continuous function of depth and estimating the vertical and lateral velocity gradients in each factorized block. If the subsurface contains strong boundaries with jumps in velocity, knowledge of the vertical velocity at a single point in a layer is sufficient for our algorithm to determine all relevant layer parameters.     

The stability problem of reverse time migration for viscoacoustic VTI media

In real strata anisotropy and viscosity extensively exists. They degraded waveforms in amplitude, resulting in which reducing of image resolution. To obtain high-precision imaging of deep reservoirs, we extended the separated viscous and anisotropic reverse time migration (RTM) to a stable viscoacoustic anisotropic RTM for vertical transverse isotropic (VTI) media, based on single generalized standard and linear solid (GSLS) media theory.. We used a pseudo-spectral method to develop the numerical simulation. By introducing a regularization operator to eliminate the high-frequency instability problem, we built a stable inverse propagator and achieved viscoacoustic VTI media RTM. High-resolution imaging results were obtained after correcting for the effects of anisotropy and viscosity. Synthetic tests verify the validity and accuracy of algorithm.     

Traveltime computation with the linearized eikonal equation for anisotropic media

A linearized eikonal equation is developed for transversely isotropic (TI) media with a vertical symmetry axis (VTI). It is linear with respect to perturbations in the horizontal velocity or the anisotropy parameter η. An iterative linearization of the eikonal equation is used as the basis for an algorithm of finite-difference traveltime computations. A practical implementation of this iterative technique is to start with a background model that consists of an elliptically anisotropic, inhomogeneous medium, since traveltimes for this type of medium can be calculated efficiently using eikonal solvers, such as the fast marching method. This constrains the perturbation to changes in the anisotropy parameter η (the parameter most responsible for imaging improvements in anisotropic media). The iterative implementation includes repetitive calculation of η from traveltimes, which is then used to evaluate the perturbation needed for the next round of traveltime calculations using the linearized eikonal equation. Unlike isotropic media, interpolation is needed to estimate η in areas where the traveltime field is independent of η, such as areas where the wave propagates vertically.
Typically, two to three iterations can give sufficient accuracy in traveltimes for imaging applications. The cost of each iteration is slightly less than the cost of a typical eikonal solver. However, this method will ultimately provide traveltime solutions for VTI media. The main limitation of the method is that some smoothness of the medium is required for the iterative implementation to work, especially since we evaluate derivatives of the traveltime field as part of the iterative approach. If a single perturbation is sufficient for the traveltime calculation, which may be the case for weak anisotropy, no smoothness of the medium is necessary. Numerical tests demonstrate the robustness and efficiency of this approach.     

Acoustic VTI modeling and pre-stack reverse-time migration based on the time–space domain staggered-grid finite-difference method

Reverse-time migration (RTM) is based on seismic numerical modeling algorithms, and the accuracy and efficiency of RTM strongly depend on the algorithm used for numerical solution of wave equations. Finite-difference (FD) methods have been widely used to solve the wave equation in seismic numerical modeling and RTM. In this paper, we derive a series of time–space domain staggered-grid FD coefficients for acoustic vertical transversely isotropic (VTI) equations, and adopt these difference coefficients to solve the equations, then analyze the numerical dispersion and stability, and compare the time–space domain staggered-grid FD method with the conventional method. The numerical analysis results demonstrate that the time–space domain staggered-grid FD method has greater accuracy and better stability than the conventional method under the same discretizations. Moreover, we implement the pre-stack acoustic VTI RTM by the conventional and time–space domain high-order staggered-grid FD methods, respectively. The migration results reveal that the time–space domain staggered-grid FD method can provide clearer and more accurate image with little influence on computational efficiency, and the new FD method can adopt a larger time step to reduce the computation time and preserve the imaging accuracy as well in RTM. Meanwhile, when considering the anisotropy in RTM for the VTI model, the imaging quality of the acoustic VTI RTM is better than that of the acoustic isotropic RTM.