VTI介质解耦合声波方程数值模拟

Alkhalifah提出的声学近似方程对于各向异性VTI介质地震资料分析处理是一种较为高效的方法,但由于算法本身所产生的伪横波干扰,以及所带来的稳定性问题制约了该技术在实际生产中的应用.本文介绍了纵横波分离的VTI介质解耦合声波近似方程,以解决伪横波干扰问题.并针对数值求解问题,研究了伪解析法在VTI介质声波近似解耦方程数值模拟中的应用.对比数值模拟结果,该方法准确描述了VTI介质qP波波场特征的同时,消除了伪SV波干扰,使得波场传播更稳定,数值频散也更小.     

VTI介质中准P波方程叠前逆时深度偏移

根据具有垂直对称轴的横向各向同性（VTI）介质中的一阶准P波方程,导出了该方程在交错网格中逆时延拓的高阶有限差分格式,给出了其稳定性条件,采用完全匹配层吸收边界条件解决边界反射问题,分别应用下行波最大能量法和归一化互相关成像条件, 实现了VTI介质中准P波方程的叠前逆时深度偏移．各向异性Marmousi模型的试算结果表明,VTI介质准P波方程叠前逆时深度偏移算法不受地下构造倾角和介质横向速度变化的限制,对复杂模型具有良好的成像能力;应用归一化互相关成像条件能得到更好的成像效果．对比该模型的各向异性和各向同性逆时偏移剖面表明,在各向异性地区采集的纵波数据用各向异性偏移算法理论上能得到更好的成像结果．      

基于GSLS模型TI介质衰减拟声波方程

随着计算机硬件技术的发展以及高分辨率勘探需求的增加,我们希望能够更准确地模拟地下介质,得到更丰富的地层信息.然而,传统的声学假设并不能描述实际地层所存在各向异性和黏滞性,使得成像分辨率较低.为了实现深部储层的高精度成像,本文同时考虑了介质的各向异性和黏滞性,从TI介质弹性波的基本理论出发,结合各向异性GSLS理论,并通过声学近似方法导出基于GSLS模型的各向异性衰减拟声波方程.数值模拟表明该方程既能准确地描述各向异性介质下的准P波运动学规律,又能体现地层的吸收衰减效应;模型逆时偏移结果表明,在实现成像过程中考虑各向异性和黏滞性的影响,能对高陡构造清晰成像,且剖面振幅相对均衡,分辨率较高.     

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

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

VTI介质中的AVO特性

理论和实践已证明,地层中广泛存在着各向异性.研究地层的各向异性,并利用它来提取用各向同性理论无法获得或更为精确地获得某些岩性信息,是一项非常有意义的课题,其中各向异性对AVO的影响是其中的一个重要的内容.本文主要从Kelvin-Chrisroffel方程出发推导出VTI介质中反射与透射系数的计算式,并通过对反射系数与各...     

粘弹性VTI介质地震波模拟特征分析

本文首先利用有限差分法分别对弹性和粘弹性VTI介质进行地震波传播数值模拟,并针对波场快照和波场记录特征,分析不同品质因子组合对波场能量衰减和频率吸收作用的影响.结果表明：对应于膨胀滞弹性形变的品质因子变化主要影响qP波的能量衰减;对应于剪切滞弹性形变的品质因子变化主要影响qSV波的能量衰减;对于qSH波,两个品质因子分别对应于垂直和水平方向的能量衰减;品质因子较小时,qSV波和qSH波的频率向低频方向移动,qP波频率变化不明显.     

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

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

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

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

基于二维VTI拟声波程函方程的椭圆各向异性立体层析反演

基于二维VTI介质拟声波程函方程,应用射线扰动理论建立了该方程控制下的数据空间各参数对于模型空间各个参数之间的线性关系,从而获得二维VTI介质拟声波程函方程的各向异性立体层析核函数.考虑到拟声波近似程函方程中η参数与ε和δ存在强烈耦合,本文首先探讨椭圆各向异性情形,为二维拟声波程函方程椭圆各向异性立体层析算法奠定了理论基础.同时也为日后推广到非椭圆各向异性情况提供了一种获得高质量初始模型的可靠途径.理论数据算例证实了Fréchet核函数求取的正确性以及在此基础上设计的工作流程实现两参数反演的可行性.

     

VTI介质中弹性波的波动特征研究

波动特征研究对各向异性介质的弹性参数反演,速度反演、NMO,DMO、时间域的偏移以及地震资料的精确解释有着重要意义.本文由Christoffel方程出发推导出二维VTI介质中XOZ对称平面内弹性波相速度与各向异性参数的关系式以及慢度面的计算式并通过理论模型分析得出以下结论：1)ε值决定了P波的水平相速度的大小,反映了P波水平传播的相速度与垂直传播的相速度的差异;δ值决定了P波垂直方向附近相速度的各向异性程度的大小.2)VTI各向异性介质中纵波偏振方向不再与慢度面是平行而是有一定的夹角,且SV波的慢度面不再完全与偏振方向垂直的.3)准SH 波的波慢度面呈椭圆,准P波与准SV波慢度面为不规则椭圆.     

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

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

     

VTI介质qP波方程高精度有限差分算子

波动方程有限差分法是一种使用广泛的地震波数值模拟方法.但是有限差分法本身固有存在着数值频散问题,会降低地震波场模拟的精度与分辨率.为了克服常规有限差分算子的数值频散,本文针对VTI介质地震波数值模拟问题,构造了频率-空间域qP波波动方程高精度有限差分优化算子,根据最优化理论中高斯-牛顿法确定了高精度有限差分算子的优化系数.利用常规差分算子和高精度优化差分算子对归一化相速度的频散关系精度进行了对比分析,并对均匀各向同性介质和均匀VTI介质中的qP波地震波场进行了有限差分数值模拟,通过频散关系精度分析和波场数值模拟结果表明：有限差分优化算子具有较高的波场数值模拟精度,有效压制了传统有限差分算子数值模拟中的数值频散现象,提高了有限差分算子精度,为VTI介质频率-空间域qP波正演模拟奠定了基础.     

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

基于三维VTI各向异性介质的频散关系,构建波数项和空间项分离的单程波算子表达式,以优化算法,确定算子的待定系数,实现广角逼近三维VTI介质的广义相移算子,发展了可灵活处理强或弱各向异性介质的波动方程叠前深度偏移方法.文中同时也针对其工业应用建议了三维VTI各向异性介质中可提高计算效率的频率相关变步长波场深度延拓算法及稀疏采样情况下可实现陡倾角构造正确成像的反假频波场延拓算法.SEG二维Hess VTI各向异性介质理论模型数据及野外数据集观测系统下的三维脉冲响应计算表明,本文提出了一种具有工业应用潜力的三维各向异性波动方程叠前深度方法.     

VTI介质中P-SV波转换点与各向异性参数关系

对于转换波地震勘探中的转换点位置这个重要问题,提出转换点位置不仅与纵横波速度比,偏移距深度比以及源检距有关,还与地下介质的各向异性的性质有关,计算了忽略地下介质的各向异性影响对转换点的确定带来的严重误差,从而影响地下地质体的精确成像.通过对层状VTI介质中的转换点近似方程的推导过程,引入该方程不同于传统方程的导出是对层状各向同性介质而言,该方程通过引入各向异性参数,使我们对转换波可以有进一步的认识,拓展了转换波处理中各向异性的应用.该方程对于偏移距深度比小于3.0的情况是比较准确的,这对于大偏移距转换波勘探具有实际意义.     

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.     

复杂介质中平面波模拟的单程波方法

地震勘探中平面波分解技术在地震资料处理与反演领域均得以广泛应用,通过正演算法模拟平面波记录有助于理解和检验平面波分解的效果和精度.本文给出一种用以模拟复杂介质中平面波的数值算法,考虑了振幅随入射角度的变化.与前人研究的平面波模拟算法相比,该算法基于波动方程的单程波解法,不但适用于横向非均匀介质,在地层倾角水平条件下还能够准确模拟振幅随入射角关系.通过MArmousi模型的计算实例,证明了本文方法的可行性.     

黏弹VTI介质井间地震高斯束波场正演数值模拟

为了研究黏弹VTI介质井间地震波的波场特征,理论分析地震波在黏弹VTI介质中传播特性的基础上,编程实现高斯射线束方法对黏弹VTI介质井间地震波场的正演模拟.为将黏弹VTI介质和完全弹性各向同性介质正演记录进行对比,将黏弹VTI介质的黏弹性参数和各向异性参数设定为零,得到完全弹性各向同性介质下的正演结果.两种介质正演结果对比显示的差异与地震波理论完全一致,证明了所研究的黏弹VTI介质井间地震高斯束正演数值模拟的正确性,为井间地震复杂介质地震波场研究提供了借鉴.