横向各向同性介质中地震波速度分析及其意义

在地震学有关领域中，地下介质的薄互层特性与裂隙的定向人布在长波长假设下会产生等效的地震各向异性特性，地震波速度随传播方向和深度发生的变化引起了广大地球物理工作者的广泛兴趣，因为速度拾取的正确与否，会直接影响到油气勘探中动校、偏移、ＡＶＯ和时深转换的结果。本文着重讨论了传统速度分析的局限性，分析了薄互层引起的横向各向同性介质中速度分析的现状和结果，实例讨论了各向异性效应对偏移处理的影响，表明各向异性     

基于BISQ机制三维双相正交裂隙各向异性介质衰减及频散方位特性研究

对于复杂裂隙介质,双相各向异性介质模型更能真实地反映实际地层特征.本文基于BISQ机制推导了三维双相各向异性介质的频散方程,给出了用于确定相速度和逆品质因子的波数方程及其表达式,研究了固体骨架各向异性、固流耦合效应各向异性、渗透率各向异性对衰减及频散的影响,并对频散和衰减的方位特性进行了分析,为预测储层流体的存在、分布以及孔隙结构提供有力的理论依据,使得多方位储层地球物理参数联合反演具有可行性.     

川西龙门山及邻区地壳上地幔各向异性地震P波走时层析成像研究

经典的地震层析成像假定介质为各向同性,通过走时反演确定地球介质的非均匀结构,得到的仅是近似地球模型。如果各向异性较强时,应用基于各向同性假设的层析成像算法,可能导致错误的结果。非均匀与各向异性效应是耦合的,如果仅考虑非均匀效应则各向异性结构可能映射到非均匀结构中造成非均匀速度结构的误差,反之亦然。因此,高分辨率的走时层析成像必须同时考虑非均匀和各向异性两种效应。同时反演非均匀和各向异性结构(即各向异性层析成像),不仅可以获得可靠的各向同性速度结构,同时可以确定各向异性结构。相对于S波偏振分析,体波走时方法基于不同的数据集,能提供独立的地震各向异性观测证据。利用地震走时层析成像方法同时研究介质的非均匀性和各向异性,对于认识地球的结构及动力学过程都具有非常重要的意义。     

粘弹各向异性介质中地震波场模拟与特征

通过引入记忆变量,可以避免粘弹性应力-应变关系中的褶积运算,使波场数值模拟易于实现.通过伪谱法对粘弹各向异性介质中的qP波、qS波数值模拟,结合理论分析,研究了粘弹各向异性介质中速度各向异性和衰减各向异性.衰减各向异性要比速度各向异性更为显著,并且qS波比qP波的衰减各向异性明显.粘弹各向异性介质中,粘弹性对波的影响主要在于波的衰减,各向异性主要影响波前面形状.     

偶极各向异性测量中钻孔应力集中的效应

由实验和理论给出了应力集中是如何影响钻孔周围速度场,以及速度场是如何影响偶极各向异性测量的。在低频率,偶极波主要易受远场应力的作用,即标准声波测井解释结果恰能得了最在应力方向。在高频率,偶极波主要受近场应力集中的影响,致命快速极化方向与最小构造应力的方向一致。这睦效应使得在快速极化方向和慢速极化方向测得的偶极波散曲线产季了交叉。根据宽频带偶极资料,偶极波散曲线交叉可被用作为应力诱发各向异性超过弱的     

不同尺度裂缝对弹性波速度和各向异性影响的实验研究

裂缝广泛分布于地球介质中并且具有多尺度的特点,裂缝尺度对于油气勘探和开发有着重要的意义.本文制作了一组含不同长度裂缝的人工岩样,其中三块含裂缝岩样中的裂缝直径分别为2 mm、3 mm和4 mm,裂缝的厚度都约为0.06 mm,裂缝密度大致相同(分别为4.8%、4.86%和4.86%).在岩样含水的条件下测试不同方向上的纵横波速度,实验结果表明,虽然三块裂缝岩样中的裂缝密度大致相同,但是含不同直径裂缝岩样的纵横波速度存在差异.在各个方向上,含数量众多的小尺度裂缝的岩样中纵横波速度都明显低于含少量的大尺度裂缝的岩样中纵横波速度.尤其是对纵波速度和SV波速度,在不同尺度裂缝岩样中的差异更明显.在含数量多的小尺度裂缝的岩样中纵波各向异性和横波各向异性最高,而含少量的大尺度的裂缝的岩样中的纵波各向异性和横波各向异性较低.实验测量结果与Hudson理论模型预测结果进行了对比分析,结果发现Hudson理论考虑到了裂缝尺度对纵波速度和纵波各向异性的影响,但是忽略了其对横波速度和横波各向异性的影响.     

TI介质各向异性速度多参数分析

本文采用横向各向同性（TI）介质弹性参数的Anderson表征方式,利用小偏移距同类反射波或转换波信息重建纵横波速,以中长排列同类反射波或转换反射波信息重建各向异性因子图像,分步进行深度域TI介质中P、SV波多参数各向异性速度分析,以理论模型验证了方法的可行性.最后给出宽角反射PP与SS波折合剖面的速度与各向异性因子解释结果.     

基于高密度速度分析的天然气水合物地震资料各向异性参数分析

在天然气水合物地震资料处理中,经过基于双曲线的高密度速度分析后,远炮检距数据仍然存在较大的剩余动校正量,分析认为:火成岩各向异性和垂向速度变化引起非双曲线时差,利用近炮检距地震数据得到均方根动校正速度,同时采用扫描法和旅行时法分别得到可以相互验证的各向异性参数,用提取的各向异性参数进行各向异性速度分析,剩余动校正得到消除,成像效果有显著改善。     

甘肃及邻近地区的各向异性介质MT反演及该区地壳深部应力场和形变带

分析了某些条件下地球介质的电导率各向异性行为，并应用一维层状对称各向异性介质大地电磁资料的反演方法，对甘肃及其邻近地区实测ＭＴ资料的各向异性模型进行反演解释，由此确定地壳深部应力场和形变带的方向。     

第二届国际地震各向异性专题讨论会

人们普通认为,地壳和上地幔中存在P波各向异性,即P波速度随方位而变化。现已观测到分层介质、韧性岩石内排列成行的晶体和脆性岩石内排列成行的裂隙的各向异性效应,并由此进行推测,使人们获得地球构造和地质力学过程方面的资料。地壳和上地幔P波地震各向异性是1982年在苏联苏兹达尔召开的第一     

TI介质偏移速度建模研究

经过多年的研究发展,各向异性叠前深度偏移算法已经趋于完善.然而,在地震资料处理过程中导致成像结果不理想的主要原因还是由于建立的地层参数场不够精确.当地层参数接近其真实值时,基于波动方程的剩余曲率建模方法由于不受构造的影响,能够在各向异性和横向变速介质中进行速度分析,所以得到了广泛的研究.本文从偏移结果中抽取共成像道集,然后通过交互运用叠前深度偏移和参数更新实现各向异性偏移速度建模.对理论模型和实际资料进行的试算表明,该方法具有较强的适应性,能极大改善VTI介质反射界面成像效果和分辨率.     

弹性各向异性介质岩石体波相速度计算

综述弹性波相速度计算的若干严格数值方法和近似解析方法,对数值算例及对不同方法作了比较.结果显示,不同数值方法得到结果一致,其适用于计算任意复杂各向异性介质岩石的相速度,但不便从物理意义的角度去分析物理参数对速度的影响;而近似解析解虽然精确程度各有不同,但是便于从物理意义的角度去进行速度分析,并在此基础上解析求取群速度这...     

应力诱导各向异性介质中地震波的传播

主要讨论了应力变化如何影响各向异性介质中波速度的问题。推导了一般各向异性介质在初始应力下的Christoffel方程,得到介质中3种波的相速度和初始应力的关系表达式;通过实验数据验证了单轴应力能够诱导各向异性,当施加单轴应力时,速度在沿应力的方向增加最大,在垂直应力的方向增加最小,实验结果与理论推导一致;用Christoffel方程的数值解模拟在3种对称情况下的弹性各向异性介质中初始应力对波速度的影响。数值结果表明:初始应力对各向异性介质中波传播速度的影响,随着各向异性强度的增加而增大,而且速度越慢,影响越大。     

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

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

TI介质局部角度域射线追踪与叠前深度偏移成像

研究与实践表明,对于长偏移距、宽方位地震数据,忽略各向异性会明显降低成像质量,影响储层预测与描述的精度.针对典型的横向各向同性(TI)介质,本文面向深度域构造成像与偏移速度分析的需要,研究基于射线理论的局部角度域叠前深度偏移成像方法.它除了像传统Kirchhoff叠前深度偏移那样输出成像剖面和炮检距域的共成像点道集,还遵循地震波在成像点处的局部方向特征、基于扩展的脉冲响应叠加原理获得入射角度域和照明角度域的成像结果.为了方便快捷地实现TI介质射线走时与局部角度信息的计算,文中讨论和对比了两种改进的射线追踪方法:一种采用从经典各向异性介质射线方程演变而来的由相速度表征的简便形式;另一种采用由对称轴垂直的TI(即VTI)介质声学近似qP波波动方程推导出来的射线方程.文中通过坐标旋转将其扩展到了对称轴倾斜的TI(即TTI)介质.国际上通用的理论模型合成数据偏移试验表明,本文方法既适用于复杂构造成像,又可为TI介质深度域偏移速度分析与模型建立提供高效的偏移引擎.     

由VSP反射波旅行时重建各向异性层速度结构

本文提出-种利用有偏VSP资料反射波旅行时信息重建椭圆各向异性介质中水平向与垂直向速度的方法。其中,地下介质假定为层状椭圆各向异性介质。反射波旅行时间采用射线追踪理论及几何关系计算得到,反演中的线性方程组采用奇异值分解(SVD)技术进行求解。 方法检测时,我们对各向同性介质及椭圆各向异胜介质情况下有限差分法正演模拟的深井有偏移距VSP地震资料分别进行各向同性和各向异性方法反演成像。结果表明,本文所述方法较之各向同性介质模型反演方法对介质类型有很好的适用性,同时也说明了本方法的司行性。最后,我们分别介绍了对实际有偏VSP资料反演得到的地下介质的速度结构图像。     

基于虚拟偏移距方法的各向异性转换波保幅叠前时间偏移

In this paper, we use the method of pseudo-offset migration （POM） to complete converted wave pre-stack time migration with amplitude-preservation in an anisotropic medium. The method maps the original traces into common conversion scatter point （CCSP） gathers directly by POM, which simplifies the conventional processing procedure for converted waves. The POM gather fold and SNR are high, which is favorable for velocity analysis and especially suitable for seismic data with low SNR. We used equivalent anisotropic theory to compute anisotropic parameters. Based on the scattering wave traveltime equation in a VTI medium, the POM pseudo-offset migration in anisotropic media was deduced. By amplitude-preserving POM gather mapping, velocity analysis, stack processing, and so on, the anisotropic migration results were acquired. The forward modeling computation and actual data processing demonstrate the validity of converted wave pre-stack time migration with amplitude-preservation using the anisotropic POM method.     

Anisotropic velocity analysis1

Results from walkaway VSP and shale laboratory experiments show that shale anisotropy can be significantly anelliptic. Heterogeneity and anellipticity both lead to non-hyperbolic moveout curves and the resulting ambiguity in velocity analysis is investigated for the case of a factorizable anisotropic medium with a linear dependence of velocity on depth. More information can be obtained if there are several reflectors. The method of Dellinger et al. for anisotropic velocity analysis in layered transversely isotropic media is examined and is shown to be restricted to media having relatively small anellipticity. A new scheme, based on an expansion of the inverse-squared group velocity in spherical harmonics, is presented. This scheme can be used for larger anellipticity, and is applicable for horizontal layers having monoclinic symmetry with the symmetry plane parallel to the layers. The method is applied to invert the results of anisotropic ray tracing on a model Sand/shale sequence. For transversely isotropic media with small anisotropy, the scheme reduces to the method of Byun et al. and Byun and Corrigan. The expansion in spherical harmonics allows the P-phase slowness surface of each layer to be determined in analytic form from the layer parameters obtained by inversion without the need to assume that the anisotropy is weak.