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

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

双相各向异性介质中波场速度与衰减特征分析

 基于BISQ模型,分别推导出二维平面内双相各向异性介质中拟快纵波(qP1)、拟慢纵波(qP2)、拟SV波(qSV)和拟SH波(qSH)的波数方程。通过计算机数值计算得到四种类型波的相速度与逆品质因子,在极坐标系中观察地震波速度和衰减随波传播方向的变化规律,发现四种波的各向异性特征表现不同,且衰减各向异性特征表现得比相速度更为强烈。文中分析了孔隙度、黏滞系数、渗透率和喷射流长度等参数在垂直方向以及全传播方位角范围对四种类型波的相速度和衰减的影响,发现这种影响也存在各向异性特征,还分析了不同参数条件下波场特征随频率变化的特性,得出一系列有意义的结论。     

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

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

三维TTI介质波动方程分解

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

双相PTL介质地震波场模拟

从基于Biot理论建立的双相各向异性弹性波方程出发,推导出二维双相PTL介质弹性波方程。在此基础上,高阶实现了高阶有限差分数值模拟,并对双相PTL介质模型和实际测井数据模型进行模拟试算。基于模拟结果,分析了各类波的传播规律,讨论了弹性参数和耗散系数对地震波传播的影响。结果表明：在双相PTL介质中存在快纵波、慢纵波和SV波,并具有各向异性,两种纵波有明显的区别,快纵波的速度远大于慢纵波的速度;慢纵波具有很强的衰减性,能否观察到慢纵波是与含流体地层介质的耗散性质有关。     

VTI介质qP波广义高阶屏单程传播算子

 本文从VTI介质弹性波波动方程入手，通过本征值方法及声学理论近似建立了VTI介质qP波频散关系方程，再根据介质扰动理论深入研究了VTI介质qP波广义高阶屏单程传播算子。将VTI介质参数场分解为均匀各向异性速度场、速度扰动场以及各向异性参数扰动场，这样qP波广义高阶屏单程传播算子就由背景介质波场自由传播项和各向异性扰动项两部分组成。背景场的计算在波数域完成，而非均匀各向异性扰动项(包括慢度扰动项、各向异性参数扰动项)在空间域计算。当各向异性参数横向变化较为剧烈时，采用高阶屏算子能提高单程qP波高角度传播精度。均匀各向异性介质和复杂各向异性介质模型脉冲响应试验表明：qP波广义高阶屏算子不仅在大角度传播时能够满足精度要求，而且能够适应任意强横向变速、强横向各向异性扰动时的情况，适用于复杂VTI介质qP波地震深度偏移。     

椭球各向异性介质弹性波完全解耦波动方程

波动方程解耦是指从弹性波方程中解耦出描述各种波型独立传播的波动方程，在地震波数值模拟、地震偏移成像和多分量地震学理论研究中具有非常重要的作用。在大多数各向异性介质中，纵波、横波是耦合在一起传播，通常不具备精确解耦性质，但椭球各向异性（EA）介质却是一个例外。首先从均匀EA介质弹性波精确频散关系方程出发，利用因式分解方法将其分解为三个解耦的频散方程；然后利用傅里叶反变换得到了均匀EA介质完全精确解耦的qP波、qSV波和SH波的波动方程。理论公式和数值模拟均表明：在均匀EA介质中，qP波、qSV波和SH波可以完全精确解耦，三种波可通过解耦的波动方程独立传播；qP波和SH波的波前面为椭球状，qSV波的波前面为球状且与各向异性参数无关；三个完全解耦的波动方程不仅适用于弱各向异性EA介质，也同样适合强各向异性EA介质。     

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

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

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

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

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

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

裂缝性多孔介质纵波频变特性研究

各向异性介质地震波衰减的研究一般只针对横向各向同性介质,即包含一组垂向裂缝。实际储层中往往存在多组正交或斜交裂缝,因此研究含多组裂缝储层的地震波频散和衰减非常必要。Chapman建立了包含两组裂缝的中观尺度喷射流模型,但是没有讨论两组裂缝对特征频率、频散和衰减幅度的影响。为此,建立了计算含多组裂缝储层纵波频变特性的方法。首先对Chapman模型数值模拟,讨论裂缝参数对特征频率和频散、衰减幅度的作用;然后分析Chapman模型和标准线性固体模型之间的关系;最后建立利用广义标准线性固体模型计算包含多组裂缝储层的频变纵波模量的方法,即利用Chapman模型和各向异性Gassmann方程得到低频模量,利用Chapman模型得到高频模量。所提方法利用低频、高频极限模量和模量损失表征地震波的频散和衰减特性,其中模量损失可表征每组裂缝的贡献,并得出以下认识：1裂缝方位和地震波入射方向只影响频散幅度和衰减幅度,而不影响频散频带和衰减频带;2不同裂缝发育情况的纵波速度V_P(或纵波逆品质因子Q_P^-1)曲线间存在代数关系,因此可以利用单组裂缝的V_P(或Q_P^-1)曲线表示两组甚至多组裂缝的V_P(或Q_P^-1)曲线;3当地震波传播方向与裂缝平行时,V_P随频率的变化很小,P波衰减最小。所提方法对于研究各向异性介质的地震波场特征及地震响应机理具有一定实际意义。     

各向异性含气砂岩模型正演及AVO响应特征分析

为认识和研究各向异性含气砂岩储层中弹性波的多波波场特征,以利于多波地震资料解释。在各向异性弹性波方程基础上,采用高阶交错网格有限差分技术模拟地震波在各向异性含气砂岩模型中的传播,得到该类储层的弹性波场。同时结合AVO分析技术分析了各向异性系数对多波波场特征的影响,比较了地震波在各向同性含气砂岩模型和各向异性含气砂岩模型中的波场响应异同。对比研究结果表明,地层的各向异性对波场特征和AVO响应有明显影响,并且这种影响在纵波和转换波上表现不一致;正演记录中的振幅特征与AVO异常响应是吻合的,说明将正演模拟与AVO分析相结合可为准确识别多波波场提供一种思路。     

双相黏弹VTI介质一阶速度—应力方程正演模拟及双程波照明研究

地层岩石具有各向异性、黏弹性及双相性等特性,因此建立更加精确的地下介质模型,研究地震波传播规律,对认识复杂介质中地震波的传播特征和实际资料的解释有着重要的实际意义.为此,基于黏弹性广义标准线性体(GSLS)模型,首先推导了双相黏弹VTI介质的一阶速度—应力方程并进行正演模拟,与单相弹性各向同性介质相比,双相黏弹VTI介...     

黏弹性介质地震波吸收衰减研究进展

实际地下介质的复杂多样对地震探测技术提出严峻的挑战.黏弹性介质理论的发展及其应用,使得波场衰减信息的利用成为地下油气储层预测的重要技术之一.在概要回顾介质非弹性吸收衰减的一些经典黏弹性模型、黏弹孔隙介质理论和黏弹各向异性理论的基础上,重点介绍了介质中地震波衰减机理的研究进展、衰减各向异性以及强衰减介质理论的提出和发展现...     

双相各向异性介质正演研究

建立了一个具有双相各向异性介质特征同时可考查慢纵波行为的特殊模型,应用交错网格有限差分方法模拟并分析了此模型的合成地震记录;此外,通过改变双相各向异性介质的固流耦合密度和孔隙度,对比并分析了介质参数对慢纵波的影响。     

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

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

各向异性粘弹性孔隙介质伪谱法波场正演模拟

在用地震资料描述油气储层时应综合考虑地球介质所表现出来的各向异性、粘弹性以及孔隙特征。本文建立了各向异性粘弹性孔隙介质模型,导出了各向异性粘弹性孔隙介质的弹性波波动方程,采用伪谱法进行了正演模拟,并分析了其波场特征。该研究有益于加深对地震波在实际地球介质中传播规律的认识。