任意空间取向TI介质中长排列非双曲动校正公式优化(英文)

常规长排列非双曲动校正公式是在VTI介质中得到的,它不能满足任意空间取向TI(ATI)条件下的扩展.本文以VTI介质中非双曲动校正公式为基础,基于我们推导得出的ATI介质中精确四次时差系数解析解和NMO速度解析解,给出ATI介质中长排列优化的非双曲动校正公式.通过与各向异性射线追踪方法计算所得出的"精确走时"结果对比,研究表明优化后的非双曲动校正公式能精确地描述任意强弱、ATI介质中随测线方位变化的走时曲线,可以用来替代耗时、多偏移距、多方位的射线追踪方法正演拟合ATI介质中长偏移距反射走时,为利用非双曲时距的各向异性参数反演提供理论基础性认识。     

几种反射波时距方程的比较

在地震资料处理中，速度分析和成像技术极为重要.常规地震资料处理方法是利用双曲线方程来描述反射波时距曲线规律，此方程随着地层非均质性、各向异性和排列长度的增加，其误差变大.目前发展的反射波非双曲时距方程，主要有基于层状各向同性模型的非双曲时距方程、基于均匀弱各向异性模型的时距方程、基于速度随炮检距变化模型的时距方程、基于线性连续速度模型的时距方程.本文针对三个典型的模型：层状均匀各向同性模型、层状弱各向异性模型和层状非均匀模型，对这几种时距方程进行了精度比较与分析，得出了一些有益的结论.最后，从不同角度说明了应用这几种方程的合理性.     

各向异性Radon变换及其在多次波压制中的应用

即使采用分辨率很高的双曲Radon变换,对速度各向异性发育介质及长偏移距情况下的地震数据,其Radon域内能量仍不收敛.为了克服此难题,我们在Radon变换的积分路径中考虑了非双曲走时的影响,通过引入非双曲时差公式中的各向异性非椭圆率η参数,可以准确描述出长偏移距条件下来自同一层位的时距曲线,并推导了由偏移距、慢度、非椭圆率三参数控制的积分曲线正反变换公式,我们称之为各向异性Radon变换.离散化求解时,各向异性Radon变换是时变的,频率域快速算法已不适用,本文采用了最优相似系数加权Gauss-Seidel迭代算法,保持其计算精度的同时也有较高的计算效率.将此方法应用在模型数据以及实际长偏移距海上地震数据的多次波压制处理中,收到了较好的处理效果.     

VTI介质长偏移距非双曲动校正公式优化

常规Alkhalifah动校正公式精度低,不能精确描述各向异性介质长偏移距地震反射同相轴的时距关系.本文以提高VTI介质长偏移距地震资料动校正公式的精度为目标,在分析VTI介质常规动校正方程的基础上,根据误差最小原理建立优化校正系数图版,实现对常规动校正公式大偏移距误差的修正,建立最优化校正Alkhalifah动校正方程,实现了对VTI介质长偏移距地震资料常规动校正方程的改进.之后由Fomel群速度公式导出高精度VTI模型长偏移距时距函数,提出了高精度VTI介质长偏移距地震资料动校正方程.将以上的动校正方程用于各向异性参数反演,模型计算表明最优化校正Alkhalifah动校正方程的反演精度是常规长偏移距动校正方程反演精度的2～4倍,高精度动校正方程的反演精度是常规动校正方程反演精度的2～8倍.     

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

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

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

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

三维VTI介质中波动方程深度偏移的最优分裂Fourier方法

从含Thomsen各向异性参数的qP波相速度表示式出发,建立并求解三维VTI介质中的频散方程,得到三维VTI介质中的相移算子,进而将以相移算子为基础的最优分裂Fourier方法推广到三维VTI介质,发展了一个三维VTI介质的深度偏移方法.文中使用的各向异性介质的速度模型与现行的各向异性构造的速度估计方法一致,将各向同性、弱各向异性及强各向异性统一在一个模型中.文中提出的偏移算法对相移法引入了高阶校正项来补偿介质横向变化的影响,使该方法可应用于横向非均匀VTI介质的陡角度成像,文中给出的偏移脉冲响应很好地证明了这一点.     

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

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

基于一阶速度-应力方程的VTI介质最小二乘逆时偏移

地下地层普遍存在各向异性,忽略介质各向异性会导致速度估计不准确,成像精度下降.基于二阶声波方程的最小二乘逆时偏移忽略了介质各向异性及密度变化的影响,致使模拟地震数据与实际观测数据不匹配,影响收敛速度和反演成像质量.VTI介质一阶速度-应力方程能较好适应各向异性变密度情况,为此,本文首先从VTI介质一阶速度-应力方程出发,进行波动方程线性化;其次推导了相应的扰动方程和伴随方程,并通过伴随状态法得到梯度更新公式;最终形成基于一阶方程的LSRTM算法理论及实现流程.在实现算法的基础上,通过数值试算及成像结果对比,验证了本文算法在处理变密度和VTI介质时的有效性和优越性.偏移速度以及各向异性Thomsen参数误差的敏感性测试及误差收敛曲线对比结果进一步表明:速度及Thomsen参数对成像结果存在明显影响,其中速度敏感性最强,参数epsilon次之,参数delta的敏感性最弱.     

VTI介质SH-SH波地震反演方法研究

近年来,随着横波可控震源技术的发展,国内外已经实现了纯横波地震勘探.相较于地震纵波,地震横波对横波阻抗、横波速度尤其是各向异性参数变化更为敏感,因此地震横波能够用来更好地估算上述地层参数.VTI(具有垂直对称轴的横向各向同性)是地层介质中广泛存在的一种各向异性形式,对振幅随偏移距变化(AVO)影响显著.本文提出了一种改进的VTI介质SH-SH波反射系数近似公式,新公式具有较高的精度且仅包含两项待求参数：横波阻抗和水平横波(SH波)速度.基于新方程建立了VTI介质SH-SH波反演方法,该方法相比VTI介质的PP波反演不确定性明显下降,同时降低了常规PP波各向异性反演对大角度数据的要求.为了获得独立的横波各向异性参数,进一步地提出了一种基于岩石物理关系的横波各向异性参数估算方法.合成地震数据测试和柴达木盆地九分量地震勘探实际地震数据应用结果表明,新方法能够准确地预测地层的横波阻抗、水平横波(SH波)速度、各向异性参数,为各向异性地层的岩性解释和油气储层预测提供了可靠的解决方案,从而深化了横波地震勘探的应用潜力.     

VTI介质中地震波反射波合成记录的方法研究

在向向异性介质中，由于地震波的相速度和群速度有较大差异，相角和群角不同，相速度、群速度及群角相与角之间的关系比较复杂，因此与各同性相比较，计算地震波相速度，群速度及群角就更加困难，其合成地震记录的难度也就之增大，本文根据VTI（具有垂直对称轴的横向各向同性）介质中地震波的运动学特征，应用射线追踪方法，计算了VTI介质中的地震波的旅行时，并合成了反 地震记录，这为反射地震波的旅行时及速度分析提供了基础。     

Pure P- and S-wave equations in transversely isotropic media

Pure-mode wave propagation is important for applications ranging from imaging to avoiding parameter tradeoff in waveform inversion. Although seismic anisotropy is an elastic phenomenon, pseudo-acoustic approximations are routinely used to avoid the high computational cost and difficulty in decoupling wave modes to obtain interpretable seismic images. However, such approximations may result in inaccuracies in characterizing anisotropic wave propagation. We propose new pure-mode equations for P- and S-waves resulting in an artefact-free solution in transversely isotropic medium with a vertical symmetry axis. Our approximations are more accurate than other known approximations as they are not based on weak anisotropy assumptions. Therefore, the S-wave approximation can reproduce the group velocity triplications in strongly anisotropic media. The proposed approximations can be used for accurate modelling and imaging of pure P- and S-waves in transversely isotropic media.     

横向各向同性介质拟声波方程及其在逆时偏移中的应用

地下岩石的速度各向异性影响地震波的传播与成像.横向各向同性(TI)介质为最普遍的等效各向异性模型.引入TI介质拟声波方程可以避免复杂的弹性波方程求解以及各向异性介质波场分离,以满足对纵波成像的实际需要.本文从垂直横向各向同性(VTI)介质弹性波方程出发,推导出正应力表达的拟声波方程以及相应的纵波分量的表达式,进而分析从频散关系得到的拟声波方程的物理意义,而后将拟声波方程扩展到更一般的倾斜横向各向同性(TTI)介质中.波前快照与群速度平面的对比验证了拟声波方程可以很好地近似描述qP波的运动学特征.在此基础上,将拟声波方程应用在逆时偏移中并与其特例声波近似方程进行对比,讨论了计算效率、稳定性等实际问题.数值试验表明VTI介质情况下采用声波近似方程可以提高计算效率,而TTI介质qP-qSV波方程则在效率相当的情况下可以保证稳定性.SEG/HESS模型和逆冲模型逆时偏移试验验证了本文TI介质拟声波方程的实用性.     

Post‐stack diffraction imaging in vertical transverse isotropy media using non‐hyperbolic moveout approximations

Diffractions carry valuable information about local discontinuities and small‐scale objects in the subsurface. They are still not commonly used in the process of geological interpretation. Many diffraction imaging techniques have been developed and applied for isotropic media, whereas relatively few techniques have been developed for anisotropic media. Ignoring anisotropy can result in low‐resolution images with wrongly positioned or spurious diffractors. In this article, we suggest taking anisotropy into account in two‐dimensional post‐stack domain by considering P‐wave non‐hyperbolic diffraction traveltime approximations for vertical transverse isotropy media, previously developed for reflection seismology. The accuracy of the final images is directly connected to the accuracy of the diffraction traveltime approximations. We quantified the accuracy of six different approximations, including hyperbolic moveout approximation, by the application of a post‐stack diffraction imaging technique on two‐dimensional synthetic data examples.     

Moveout approximation for vertical seismic profile geometry in a 2D model with anisotropic layers

I introduce a new explicit form of vertical seismic profile (VSP) traveltime approximation for a 2D model with non‐horizontal boundaries and anisotropic layers. The goal of the new approximation is to dramatically decrease the cost of time calculations by reducing the number of calculated rays in a complex multi‐layered anisotropic model for VSP walkaway data with many sources. This traveltime approximation extends the generalized moveout approximation proposed by Fomel and Stovas. The new equation is designed for borehole seismic geometry where the receivers are placed in a well while the sources are on the surface. For this, the time‐offset function is presented as a sum of odd and even functions. Coefficients in this approximation are determined by calculating the traveltime and its first‐ and second‐order derivatives at five specific rays. Once these coefficients are determined, the traveltimes at other rays are calculated by this approximation. Testing this new approximation on a 2D anisotropic model with dipping boundaries shows its very high accuracy for offsets three times the reflector depths. The new approximation can be used for 2D anisotropic models with tilted symmetry axes for practical VSP geometry calculations. The new explicit approximation eliminates the need of massive ray tracing in a complicated velocity model for multi‐source VSP surveys. This method is designed not for NMO correction but for replacing conventional ray tracing for time calculations.     

Traveltime approximation in vertical transversely isotropic layered media

The well‐known asymptotic fractional four‐parameter traveltime approximation and the five‐parameter generalised traveltime approximation in stratified multi‐layer transversely isotropic elastic media with a vertical axis of symmetry have been widely used for pure‐mode and converted waves. The first three parameters of these traveltime expansions are zero‐offset traveltime, normal moveout velocity, and quartic coefficient, ensuring high accuracy of traveltimes at short offsets. The additional parameter within the four‐parameter approximation is an effective horizontal velocity accounting for large offsets, which is important to avoid traveltime divergence at large offsets. The two additional parameters in the above‐mentioned five‐parameter approximation ensure higher accuracy up to a given large finite offset with an exact match at this offset. In this paper, we propose two alternative five‐parameter traveltime approximations, which can be considered extensions of the four‐parameter approximation and an alternative to the five‐parameter approximation previously mentioned. The first three short‐offset parameters are the same as before, but the two additional long‐offset parameters are different and have specific physical meaning. One of them describes the propagation in the high‐velocity layer of the overburden (nearly horizontal propagation in the case of very large offsets), and the other characterises the intercept time corresponding to the critical slowness that includes contributions of the lower velocity layers only. Unlike the above‐mentioned approximations, both of the proposed traveltime approximations converge to the theoretical (asymptotic) linear traveltime at the limit case of very large (“infinite”) offsets. Their accuracy for moderate to very large offsets, for quasi‐compressional waves, converted waves, and shear waves polarised in the horizontal plane, is extremely high in cases where the overburden model contains at least one layer with a dominant higher velocity compared with the other layers. We consider the implementation of the proposed traveltime approximations in all classes of problems in which the above‐mentioned approximations are used, such as reflection and diffraction analysis and imaging.     

VTI介质纯P波混合法正演模拟及稳定性分析

各向异性介质纯P波方程完全不受横波的干扰,在一定程度上可以减缓由于介质各向异性引起的数值不稳定,本文推导了具有垂直对称轴的横向各向同性(VTI)介质纯P波一阶速度-应力方程.由于纯P波方程存在一个分数形式的伪微分算子,无法直接采用有限差分法求解.针对该问题,本文采用伪谱法和高阶有限差分法联合求解波动方程,重点分析了混合法求解纯P波一阶速度-应力方程的稳定性问题,并给出了混合法求解纯P波方程的稳定性条件.数值模拟结果表明纯P波方程伪谱法和高阶有限差分混合法能够进行复杂介质的正演模拟,在强变速度、变密度的地球介质中仍然具有较好的稳定性.     

多层垂直对称轴横向各向同性介质精确走时计算

给出了计算多层垂直对称轴横向各向同性（VTI）介质精确射线路径和走时的方法,所用的体波相速度公式、群速度公式和Snell定律都是严格的显式解析公式. 任意基本波的射线路径和走时计算问题都可以转化成一个等效的透射问题,再用文中的公式来计算,具体实现方法用一个多次波和一个首波的实例给出. 最后分别用精确公式和Thomsen近似公式计算了相同模型相同基本波的走时曲线. 比较两者计算结果可发现, 近似公式反复使用会使误差积累,同时揭示了近似公式适用范围的局限性,强调了使用近似公式需要注意其适用范围的重要性.      

Velocity anisotropy characteristics and pre-stack time migration imaging of the marine sedimentary strata in the South Yellow Sea Basin

The South Yellow Sea is a superimposed basin overlying Mesozoic-Cenozoic continental sediments, which in turn overlie Paleozoic-Mesozoic marine deposits that are now the target of hydrocarbon exploration. Strongly modified by multiple tectonic events, the marine sediments feature a large tectonic relief, with obvious horizontal anisotropy in seismic velocity, which significantly affects the seismic image quality. In this study, the sedimentary velocity anisotropy and its influence on image quality were analyzed using an analytical theory method, assuming transversely isotropic medium with vertical axis of symmetry (VTI), and using seismic and well-log data. Additionally, an anisotropic prestack time migration was used for the imaging of the field data. The results showed that the anisotropic pre-stack time migration processing could be used to significantly improve the accuracy of the seismic images in areas with distinct faults, offering clear images of accurately located fault planes and fault edges, thereby improving the lateral resolution of the seismic data and its signal-to-noise ratio.