三维TTI介质弹性波相、群速度的近似配方表征法

TTI介质弹性波相速度可以通过求解Christoffel方程得到,其精确式往往比较复杂,难以运用.为此,从Thomsen各向异性参数表征的TTI介质弹性波精确三维相速度表达式出发,利用近似配方法推导了 TTI介质qP波和qSV波三维相速度的近似公式;将其代入Berryman的群速度表征中,得到了qP波和qSV波群速度的...     

倾斜椭球各向异性介质弹性波传播特征

弹性波的传播特征包括速度和偏振等,它们是研究地震波传播规律的重要参数,在理论研究和实际应用中具有重要作用。椭球各向异性介质是横向各向同性介质的一种特例。从观测坐标系下倾斜椭球各向异性(TEA)介质刚度矩阵出发,根据Christoffel方程推导了倾斜椭球各向异性介质弹性波速度和偏振方向的表达式。理论分析和数值示例表明:①在TEA介质中,qP波和SH波相速度呈各向异性特征,相速度等值面不是椭球面,而两者的群速度等值面均为对称轴倾斜的椭球面;qSV波相速度和群速度完全相同,并且呈各向同性特征。②qP波、qSV波和SH波的偏振方向两两垂直,qP波和qSV波的偏振方向与各向异性参数有关;SH波偏振方向与各向异性参数无关,并且垂直于传播方向和对称轴方向确定的平面。     

三维TTI介质中的弹性波场分解

根据三维TTI介质倾斜对称轴与观测坐标系之间的几何关系,经坐标旋转,可以将其近似看成三维VTI介质,故适用于VTI介质的波场分解方法就可以应用于TTI介质。为此,从三维各向异性弹性波动方程出发,通过求解Christoffel方程的特征向量,构建了三维VTI介质的波场分解算子;再利用坐标旋转推导了三维TTI介质的波场分解算子;最后通过泊松方程实现三维TTI介质的弹性波场分解。引入改进的快速算法,避免了直接求解泊松方程。推导的分解算子考虑了随空间变化的弹性参数和倾斜对称轴,可适用于复杂的三维TTI介质模型。当各向异性参数以及对称轴倾角都为零时,该算子可以退化为各向同性的形式。数值计算结果表明,该TTI介质波场分解算子比各向同性或VTI介质算子能得到更精确的分离结果,可高效实现三维各向异性矢量弹性波场分解。     

双相各向异性介质中的拟慢纵波研究

通过双相各向异性介质理论,利用弹性波传播方程,推导出三维双相各向异性介质中弹性波的Christoffel方程。在Christoffel方程基础上,分别计算并分析了频率、固流相耦合密度、渗透率对双相各向异性介质中拟慢纵波的相速度、衰减的影响。     

三维TTI介质波动方程分解

 波动方程分解是指从弹性波方程中分解出描述各种波场独立传播的波动方程。在各向异性介质中,由于纵波、横波是耦合在一起传播的,所以通常不具备解耦性质。本文从Thomsen弱各向异性近似和声学假设近似两个途径对三维TTI介质弹性波波动方程进行了分解。利用本征值方法求解三维TTI介质弹性波的Christoffel方程,得到描述SH波、qP波和qSV波的精确频散关系方程,通过Thomsen弱各向异性表征理论,推导出了弱各向异性条件下描述qP波和qSV波传播的时空域波动方程;由TTI介质qP波和qSV波耦合的频散关系方程出发,根据声学假设原理,推导出三维TTI介质描述qP波传播的波动方程。数值试算表明,通过波动方程分解获得的三维TTI介质qP波和qSV波波动方程具有较高的精度,为研究三维TTI介质qP波和qSV正演模拟和深度偏移算法奠定了理论基础     

TTI介质qP波伪谱法正演模拟

基于Tsvankin的VTI介质qP波精确相速度公式,利用坐标变换的方法,导出了三维空间TTI介质qP波的精确相速度公式;结合Thomsen的弱各向异性假设以及泰勒展开式,简化了三维TTI介质qP波的精确相速度公式,建立了近似相速度公式;TTI模型的计算结果表明,近似相速度公式能较好地逼近精确相速度公式。从近似相速度公式出发,推导了三维qP波的频散关系以及时间-波数域波动方程,构建了纯qP波时间-波数域波场递推格式,分析了递推公式的稳定性条件。最后,基于伪谱法对三维均匀（VTI、TTI）模型、Hess VTI模型以及BP 2007 TTI模型进行qP波正演模拟。结果表明,利用所提方法得到的波场快照没有伪横波干扰,波场计算结果稳定。     

基于最佳平方逼近的TTI介质解耦qP波与qSV波逆时偏移

在实际地震资料处理中,若忽略各向异性影响、TTI介质对称轴倾角变化,会造成目标体成像位置偏差,甚至由焦散问题造成振幅不均衡。传统TI介质耦合拟声波方程正演模拟与逆时偏移存在伪横波干扰,且在各向异性参数ε<δ时数值模拟不稳定。为此,借鉴前人经验,首先使用最佳平方逼近（OQA）近似VTI介质精确相速度公式,分别得到较高精度的近似qP波、qSV波相速度公式。然后从近似相速度公式出发,推导了二维TTI介质解耦qP波、qSV波波动方程。为实现高效、稳定的数值模拟,使用有限差分-伪谱混合法求解波动方程。最终,建立了高精度、高效率、稳定的TTI介质正演模拟与逆时偏移成像算法。数值模拟实验表明,所提方法可实现复杂TTI介质偏移成像。     

TTI介质qP波方程频率—空间域加权平均有限差分算子

 波动方程有限差分方法能够较精确地模拟任意非均匀介质中的地震波场，但它本身存在着数值频散问题。在具有倾斜对称轴的横向各向同性介质(TTI介质)地震波正演模拟中，为了解决常规有限差分算子的数值频散问题，本文构造了频率—空间域qP波方程加权平均有限差分算子，求取了归一化相速度，并根据最优化理论中的高斯—牛顿法确定了加权平均差分算子的最优加权系数。利用常规差分算子和加权平均差分算子对归一化相速度进行了频散分析，并对均匀TTI介质(包括各向同性介质和椭圆各向异性介质)中的qP波地震波场进行了有限差分数值模拟。结果表明：加权平均有限差分算子具有较高的数值精度，能有效地压制常规有限差分算子的数值频散，为TTI介质频率—空间域qP波正演模拟奠定了基础。     

TTI介质一阶qP波方程交错网格伪谱法正演模拟

受TTI声介质强各向异性和对称轴极化倾角急剧变化的影响,传统有限差分法正演模拟易出现不稳定现象。本文基于TTI介质二阶耦合准纵(qP)波方程,推导了TTI介质一阶速度-应力方程,构建了交错网格中基于完全匹配层(PML)吸收边界条件的高阶有限差分格式;通过引入基于PML原理的一阶qP波方程伪谱法递推格式,实现了TTI介质qP波波场的稳定延拓。对二维TTI介质Wedge和BP模型的正演模拟试算结果表明:基于伪谱法递推格式实现的TTI介质正演模拟法,能较准确且稳定地模拟强各向异性及倾角变化剧烈区域的地震波场,验证了本文方法的有效性和较强适应性。     

TTI介质相速度对弹性模量参数的敏感性研究

速度敏感函数是相(群)速度对各向异性介质参数的偏微分表达式,它是各向异性介质走时扰动方程的重要组成成分,所以具有重要的研究价值。由于群速度是相速度的函数,所以利用相速度偏微分表达速度敏感函数更加简便。参考前人基于VTI介质模型的速度敏感函数计算方法,推导出了基于TTI介质模型以相速度偏微分表达的速度敏感函数计算方法,并对VTI介质和TTI介质的qP波、qSV波和qSH波分别模拟计算了速度敏感函数。对速度敏感函数的特性分析结果表明,不同弹性参数的速度敏感性不同,同一弹性参数不同观测角度的速度敏感性也不同,VTI介质的速度敏感函数存在对称性,而TTI介质的速度敏感函数相对复杂得多。最后利用TI介质反射波非线性走时反演算法对层间块状异常体模型进行了走时反演,反演结果验证了速度敏感函数计算方法的正确性,也证实了速度敏感函数对反演结果具有预测和指示作用。     

TTI介质弹性波相速度的一般性近似

弹性波速度是弹性波动力学研究的核心内容之一。倾斜横向各向同性(TTI)介质弹性波相速度的理论精确公式比较复杂。由于地球介质弱各向异性的普遍性,在理论研究和实际应用中常常采用弱各向异性假设进行简化,但是面对页岩等中强各向异性介质时,弱各向异性近似不能满足精度要求,甚至不再适用。基于前人的研究成果,从TTI介质弹性波相速度理论精确公式出发,利用近似的配方法将精确公式内的D项近似表示为两项的平方和,从而实现多项式的开方运算。在此基础上,推导了TTI介质qP波和qSV波相速度的一般性近似表达式。理论分析和数值计算结果表明:TTI介质弹性波相速度一般性近似表达式由于在推导过程中没有直接利用弱各向异性的假设,因此不仅适用于弱各向异性,也适用于中强各向异性的情况,即在较大范围内,一般性近似的精度高于弱各向异性近似的精度。并且对于同一种TTI介质,qP波的一般性近似式精度高于qSV波的精度。     

双相各向异性介质弹性波传播交错网格高阶有限差分法模拟

含流体孔隙各向异性介质中的弹性波传播问题是目前石油勘探和地震学研究的难点之一。基于Biot理论,本文提出了二维双相任意倾斜各向异性介质三分量弹性波方程交错网格任意偶阶精度有限差分解法,并对均匀及两层双相VTI介质和TTI介质中的弹性波场进行了模拟。波场快照和合成记录表明:当横波源的偏振方向与TI介质的对称轴存在一定夹角、且耗散系数较小时,在固相与流相中存在快横波、慢横波、快纵波和慢纵波;双相TI介质中横波分裂现象及波前面尖角和三分叉现象与TI介质中的基本相同,只是黏滞相界中弹性波速度有所减小,且具有频散和衰减特性;从直达波转换为反射波后,该波型具有反射波的性质。     

TTI介质一阶qP波稳定方程波场数值模拟及逆时偏移

相比于各向异性介质弹性波波动方程,利用声学近似的各向异性介质qP波方程进行波场数值模拟及逆时偏移更具优势。常规的声学近似方法往往会造成非均匀TTI介质中倾角剧变区域出现数值不稳定。为此,基于精确的TTI介质qP-qSV波耦合频散关系,首先引入一个各向异性控制参数σ,推导了新的TTI介质二阶qP波稳定方程,并通过引入波场的伪速度分量,将其转换为等价的一阶应力-速度形式波动方程。然后,利用优化的最小二乘交错网格高阶有限差分(LS-SGFD)方法数值求解TTI介质一阶qP波稳定方程,构建波场延拓算子,实现了精确的各向异性介质波场模拟及逆时偏移成像。模型试算结果表明,TTI介质一阶qP波方程能够稳定地模拟qP波的波场传播特征,利用优化的LS-SGFD方法能够有效地提高波场模拟的精度,进一步可以改善偏移成像质量。     

正交介质弹性本构关系及角散的研究

本文首先讨论了双层EDA裂隙介质的各向异性弹性本构关系，并导出了基于Postma方法的正交EDA介质的等效弹性常数求解公式和角散公式；然后简要介绍了正演伪谱法算法；最后分别利用角散公式和伪谱法对正交介质的波速与走时进行了数值模拟计算，获得了一致的结果，即正交介质的S波分裂现象比其他介质类型明显。     

TTI介质各向异性伪谱法逆时偏移

在Thomsen精确相速度方程和弱各向异性假设前提下,将VTI相速度方程拓展为TTI介质相速度方程,进而得到TTI介质的P波方程。利用伪谱算法对P波方程进行了正演和逆时偏移,并和各向同性F-X域深度偏移效果进行了对比,结果表明:①对于各向异性介质,利用各向同性的偏移算法得到的结果会有误差,成像结果有明显的畸变;②在考虑效率和精度的前提下,伪谱法由于其利用快速二维傅里叶变换解决空间导数的计算问题,这对于实际应用中计算成本颇高的逆时偏移是一个较理想的选择;③推导的近似方程在弱各向异性条件下可以较好地应用于TTI介质中。     

倾斜横向各向同性介质中弹性波波场伪谱法模拟

各向异性可以提供有关地层构造和岩石性质的大量信息,通过对各向异性地层的波场特性研究,就可以根据地震波的不同特征来分析各种类型的各向异性的成因,所以这是认识地震各向异性的重要途径之一.根据Taylor展开式,给出了倾斜各向异性介质中的位移波动方程,并应用伪谱法进行数值模拟.模拟结果表明(1)TTI(倾斜横向各向同性)介质中的弹性波包括qP波、qSV波、qSH波、qPSV波、qPSH波,透射波有P波、SH波、SV波和转换波,同时波前的对称轴方向发生偏转.(2)伪谱法能较好地胜任正演模拟,能完全推广到其他各向异性介质波场模拟中.     

各向异性介质弹性波高斯束偏移方法

弹性波高斯束偏移是一种兼顾成像精度与计算效率的多分量地震成像方法,但当前的研究多集中在各向同性介质声波、弹性波成像和各向异性介质声波成像,针对各向异性介质弹性波偏移的文献相对较少。为此,在前人的研究基础上,发展了各向异性介质弹性波射线追踪方程,以二维各向异性弹性波Kirchhoff-Helmholtz积分为基础,利用弹性动力学表征的格林张量推导了各向异性介质中弹性波波场正、反向延拓公式,提出了一种针对各向异性介质多分量地震数据的成像方法。通过在成像公式中引入权函数,很好地压制各向异性介质成像中不同波型引起的串扰。通过引入符号函数,解决了转换波成像过程中的极性反转问题。TTI介质洼陷模型和TTI介质多层构造模型试算结果表明：与各向同性介质弹性波高斯束偏移方法相比,各向异性介质弹性波高斯束偏移成像结果的信噪比较高、同相轴连续性较好、构造位置更准确;使用多分量地震记录,利用纵、横波波场信息成像,有效提高了成像分辨率;与传统的基于标量波场的成像方法相比,所提方法的适应性更好。     

简化的三维TTI介质黏滞纯声波方程及其数值模拟

地球介质普遍具有非弹性和各向异性特征,因此在研究地震波在地下空间传播时应该同时考虑各向异性和黏滞性。在各向异性介质波场模拟、偏移成像和波形反演中,目前主要采用的是伪声波方程,该类方程是在直接将横波速度设置为0的基础上发展的,当介质参数不满足假设条件时容易产生伪横波数值干扰及模拟不稳定。考虑到伪声波方程存在的问题,文中应用泊松算子和有限差分相结合的策略求解高精度的三维TTI介质纯声波方程。同时,考虑到衰减介质对地震波振幅和相位的影响,在各向同性黏滞声波方程的基础上,推导了一种简化的三维TTI介质黏滞纯声波方程,该方程能够模拟纯声波的相位畸变和振幅衰减。应用三维层状模型、TTI楔状体模型和改进的Marmousi模型验证了方法的有效性和适用性。     

基于Lebedev网格的TTI介质二维三分量正演模拟

传统交错网格有限差分法是研究地震波传播规律的一种较为常用的数值模拟方法,但是在交错网格中每个变量的不同分量都是交错定义的,对于没有定义的点需要进行变量插值,从而降低了模拟精度。为此,在前人的研究基础上,推导了TTI介质二维三分量的应力—速度方程,采用Lebedev网格对其进行了高精度的差分离散处理,避免了传统交错网格在处理各向异性介质时波场插值引起的误差,提高了模拟精度,并将多轴完全匹配层吸收边界（M-PML）引入了Lebedev网格。分别对单层TTI介质和含透镜体的复杂TTI介质模型进行了正演模拟,结果表明：①由于Lebedev网格对各向异性介质的弹性波方程做离散时无需进行波场插值,与传统交错网格有限差分方法相比,模拟精度更高;对TTI介质进行模拟时可以清晰地观察到纵波、快横波和慢横波,并且快、慢横波的偏振方向相反,在单炮记录中观测到的三种波的速度特征也符合波场传播规律。②引入多轴完全匹配层（M-PML）吸收边界条件后,在不影响模拟效果的情况下边界反射现象被有效地压制。     

双相各向异性介质地震波传播特征

实际裂隙性油气储层表现出双相各向异性特征。本文基于岩石物理理论求取裂隙介质各向异性等效弹性参数,结合各向异性BISQ理论,建立了裂隙诱导双相各向异性介质理论模型。从双相各向异性介质的弹性波方程出发,推导出平面波的Christoffel方程,求解该方程获得了三维双相完全各向异性介质中qP1、qP2、qSV和qSH波的波数方程,据此波数方程计算四类波的相速度和衰减因子。最后针对碎屑岩水平裂隙诱导双相各向异性介质,分析了泥质含量、裂隙密度和形状、孔隙度、含水饱和度、流体黏滞系数、渗透率等参数对弹性波沿不同方向传播时相速度和能量衰减特征的影响,可为储层的孔隙、裂隙结构及含油气性等分析提供理论依据。