黏弹TTI介质中旋转交错网格高阶有限差分数值模拟

以Carcione黏弹各向异性理论为基础,给出了适用于黏弹性具有任意倾斜对称轴横向各向同性介质(黏弹TTI介质)的二维三分量一阶速度-应力方程,采用旋转交错网格任意偶数阶精度有限差分格式求解该方程,并推导出了二维黏弹TTI介质完全匹配层(PML)吸收边界条件公式和相应的旋转交错网格任意偶数阶精度有限差分格式,实现了该类介质的地震波场数值模拟.数值模拟结果表明:该方法模拟精度高,边界吸收效果好,可以得到高精度的波场快照和合成记录;并且波场快照和合成记录能较好地反映地下介质的各向异性特征和黏弹性特征.     

包含Biot机制和分数阶黏弹性机制的波传播模型

含流体孔隙介质中的波能量耗散通常由多种力学机制造成.传统Biot理论中的能量耗散仅仅考虑了固流两相相对运动引起的摩擦耗散,无法准确预测波在孔隙介质中低频段出现的高频散与强衰减现象.为了建立一个能准确预测地震波频段高频散与强衰减现象的动力学模型,我们在Biot理论的基础上引入黏弹性机制,并利用分数阶导数刻画黏弹性本构关系,最终获得了一种新的孔隙介质波传播模型.与传统的Biot模型相比,新模型考虑了含流体孔隙介质中固体骨架的内耗散,对波能量耗散的刻画更为精准.通过数值算例,我们研究了分数阶导数的阶数参数对快P波和S波频散和衰减的影响,并通过来自不同地区且具有不同物理性质的几组流体饱和岩芯实验数据,对比研究了新模型的有效性.结果表明,文章提出的新模型能更准确地预测快P波和S波在低频段出现的高频散和强衰减现象.     

基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移

时间域常Q黏声波方程,由于含分数阶时间导数项,数值求解需要大量内存,计算效率低,不利于地震偏移的实施.通过一系列近似,可将该方程简化为介质频散效应和衰减效应解耦的分数阶拉普拉斯算子黏声波方程,数值求解内存需求少,计算效率高.本文采用交错网格有限差分逼近时间导数,改进的伪谱法计算空间导数,PML吸收边界去除边界反射,对该方程进行数值离散和地震正演模拟,开展地震数据的黏声介质逆时偏移,实现波场逆时延拓过程中同时完成频散校正和衰减补偿.改善深层构造的成像精度,数值结果表明,基于分数阶拉普拉斯算子解耦的黏声介质地震正演模拟与逆时偏移可大幅度提高地震模拟计算效率,偏移剖面明显优于常规声波偏移剖面,极大改善深层构造的成像品质.     

黏弹介质波动方程有限差分解的稳定性研究

稳定性问题是地震波数值模拟的一个重要问题.基于地震波传播理论,从黏弹介质本构方程出发,对矩形网格下不同黏弹模型波动方程有限差分解的稳定性进行了理论分析,导出了Kelvin-Voigt黏弹模型和Maxwell黏弹模型在任意空间差分精度下稳定性条件的表达式;给出了品质因子Q≥5时的简化式,并通过数值算例验证了理论研究结论的正确性;总结了地震波速度、频率、空间网格大小、差分系数以及品质因子与稳定性条件的关系;通过误差分析给出了近似公式的使用条件.     

黏弹双相介质中的松弛骨架模型

本文基于Biot理论,考虑了多孔介质中固体骨架的松弛特征,引入了纵波品质因子、横波品质因子与耗散品质因子三参数来描述黏弹双相介质波动方程.采用虚谱法在地震频段进行了波场模拟,模拟结果表明：松弛骨架机制不仅适用于高频段,也可用于解释地震频段下的弹性波衰减现象,以描述固体微细颗粒的中观松弛特征.结合小生境遗传算法对三层模型24个介质参数进行了反演,反演结果表明：无噪波场的反演结果具有较高精度,对于含噪波场,取值在奇异点附近的介质参数反演精度降低.最后,对中国东部某地区的实际资料进行了浅层参数的反演,得到了该地区的表层固体体积模量、固体密度以及品质因子.     

基于分数阶时间导数的黏弹性衰减VTI介质中平面波传播

目前在地震勘探频带范围内通常假设品质因子Q与频率无关,且呈衰减各向同性.事实上,相比较速度各向异性,介质的衰减各向异性同样不可忽视.本文将衰减各向异性和速度各向异性二者与常Q模型相结合,建立了黏弹性衰减VTI介质模型,并基于分数阶时间导数理论,给出了对应的本构关系和波动方程.利用均匀平面波分析和Poynting定理,推导出准压缩波qP、准剪切波qSV和纯剪切波SH的复速度、相速度、能量速度以及品质因子的解析表达式.对模型的正确性进行了数值验证,并分析了qP,qSV和SH波在介质中的传播特性.数值试验结果表明:本模型能够实现理想的恒定Q行为,表现了品质因子和速度的各向异性特征,显示出黏弹性增强将导致能量速度和相速度的频散曲线变化剧烈;速度和衰减各向异性参数与传播角度之间的耦合效应对qP,qSV和SH波的速度和能量影响明显;qP,qSV和SH波的频散曲线和波前面随着衰减各向异性强度的改变发生显著变化,其中耦合在一起的qP和qSV波变化趋势相同,而SH波与它们呈现相反的变化规律.本研究为从常Q模型角度分析地震波在衰减各向异性黏弹性介质中的传播特征奠定了理论基础.     

基于Zener线性体的黏弹孔隙介质衰减频散特征分析

地震波传播经过储层时,表现出明显的衰减和频散特征.该现象的主要成因应该包含两种机制:一是固体骨架的黏弹机制,二是孔隙流体的流动机制.这两种机制在单独研究时具有各自的局限性.Zener线性体模型描述的黏弹机制将储层等效为一个固体,忽略了孔隙中流体流动造成的地震波衰减;而BISQ模型考虑了储层的双相结构以及宏观及微观的孔隙流体流动机制,但其模拟衰减频带处于高频,不能解释在地震频带明显的衰减现象.本文分析了Zener模型和BISQ模型的作用机制,并将两种机制结合,推导了Zener-BISQ模型的本构方程.通过求解新模型的平面波解,分析了其相速度和逆品质因子的衰减和频散特征.最后,通过3个模型的数值实验研究了Zener模型其中的作用和控制方式.     

常Q滞弹性介质地震波动数值模拟——时域本构优化逼近

在地震波动模拟中计入常Q滞弹性阻尼,可有效降低模拟波形的误差.就时域有限差分和有限元模拟而言,常基于广义标准线性体建立阻尼介质的时域本构逼近.广义标准线性体由若干标准线性体并联得到,增加标准线性体个数能有效提高模拟精度,但计算量及计算存储将成倍增长.目前尚未有普适的标准线性体个数优化取值方案.本文基于广义标准线性体参数的非线性优化拟合方法,详细分析了时域本构逼近误差的影响因素,清楚揭示了逼近误差仅取决于频带宽度,与频带上下限取值无关这一特性,阐明了构建具有普适性标准线性体个数优化取值方案的可行性.论证了波形模拟精度主要取决于波传播距离与模拟波长的比值以及标准线性体的个数取值.综合考虑上述两个控制因素,结合在波动正反演问题中广为采纳的波形时频误差衡量准则,对不同Q值介质给出了标准线性体个数优化取值表.进一步,本文提出采用不同个数标准线性体以近似不同Q值的阻尼介质时域本构,解决了以往波动数值模拟中统一采用相同个数标准线性体而导致的计算量及计算存储浪费或模拟精度低下等问题,并基于数值实验验证了这一方法的精度.本文工作对推进滞弹性介质波动数值模拟及其在全波形反演问题中的应用具有理论价值和实践意义.     

黏弹各向异性介质中波的反射与透射问题分析

黏弹各向异性介质中传播不均匀波，其反射、透射模式不仅与介质分界面两侧速度对比有关，还与品质因子Q的对比有关. 用伪谱技术模拟黏弹各向异性介质分界面上波的反射、透射，并与弹性各向异性介质、黏弹各向同性介质和弹性各向同性介质的模拟结果做比较. 计算平面波的反射、透射系数，分析介质的黏弹性和各向异性对反射、透射系数的影响. 数值模拟了一个三层介质模型中的波场，分析两个分界面上产生的反射波的特征. 黏弹各向异性介质中，qS波比qP波衰减程度大.     

双相各向异性介质中弹性波方程的有限元解法及波场模拟

基于双相各向异性介质模型,首先推导了双相各向异性介质中弹性波传播的动力学方程及其Galerkin变分方程和有限元运动方程,然后给出了孔隙弹性波方程的有限元数值解法以及二维双相PTL介质中波场模拟的人为吸收边界条件. 最后,利用本文给出的有限元方法对双相PTL介质和双相各向同性介质中的弹性波传播进行了数值模拟. 结果表明:有限元方法和吸收边界条件有效、可行,在理想相界条件下,不论是从固体位移,还是从流体位移的波场快照都能看到明显的慢速拟P波;在黏滞相界情况下,能否观察到慢速拟P波,与含流体地层介质的耗散性质有关.对实际含流体介质,从流体位移分量的波场快照比从固体位移波场快照更容易观察到慢速拟P波.     

非均匀黏弹性介质中声波模拟的非规则网格方法

本文研究了满足衡Q模型的黏弹性介质中声波模拟的时间域数值解法.通过采用有理式基函数描述频率相关的体积模量,使模型满足地震波黏性吸收的衡Q要求.结合弹性波模拟的非规则、非结构网格方法——格子法,本文提出了一种非均匀黏弹性介质中声波模拟的非规则网格方法.非规则、非结构网格的使用,可以精细地刻画地下介质中复杂速度、Q值的变化,及界面的形状.通过和解析解的比较及复杂模型算例分析,验证了该方法的精度及有效性.     

裂缝诱导双相HTI介质地震波场错格伪谱法模拟与波场特征分析

裂缝诱导的双相具有水平对称轴的横向各向同性(HTI)介质模型是由一组平行排列的垂直裂缝嵌入到统计各向同性的流体饱和多孔隙岩石中而组成的,它综合考虑了裂缝型储层岩石的各向异性和孔隙性.高精度的地震波场数值模拟技术是研究该介质中地震波传播规律的主要方法.本文结合错格伪谱法和时间分裂法,求解描述该介质中地震波传播的一阶速度-应力方程.模拟了单层和双层模型中的地震波场,并对其进行了特征分析.研究结果表明:错格伪谱法能有效消除标准网格伪谱法波场模拟结果中出现的数值伪影现象,与时间分裂法结合能够获得稳定的、高精度的模拟结果;裂缝诱导双相HTI介质中的地震波场兼具裂缝各向异性介质和双相介质中传播的地震波的波场特征.     

方位各向异性粘弹性介质波场数值模拟

当地震信号通过复杂地球介质时，地层除了表现为各向异性，还表现为内在的粘弹性特征.因此，为准确描述地震波在地球介质中的传播特征，理想的地球介质模型应该能够模拟岩石的各向异性特征和衰减特征.本文给出了各向异性粘弹性介质模型的波动方程及其差分格式，并利用有限差分法实现了地震波波场数值模拟.结果表明了该介质模型中地震波场特征与各向异性主轴方位和介质的粘滞性参数之间的关系.     

基于改进BISQ机制的双相HTI介质波传播伪谱法模拟与特征分析

改进BISQ(Biot-Squirt)机制在不引入特征喷流长度的情况下,将含流体孔隙介质中Biot流动和喷射流动两种重要的力学机制有机地结合起来,且各相关参数具有明确物理意义和可实现性.本文将改进BISQ机制一维孔隙流体压力公式推广到三维具有水平对称轴横向各向同性介质(HTI介质)情况,结合裂缝各向异性理论,给出了基于改进BISQ机制的双相HTI介质模型及其二维三分量波传播方程,采用伪谱法求解该方程,进行了不同相界、不同频率以及双层地质结构情况下该类介质中波场的数值模拟与特征分析.数值模拟结果表明:伪谱法模拟精度高,压制网格频散效果好,可以得到高精度的波场快照和合成记录;基于改进BISQ机制的双相HTI介质模型兼具裂缝各向异性特征和孔隙弹性特征,其为从双相各向异性理论角度深入研究裂缝性储层的地震响应奠定了理论基础.     

黏弹性VTI介质频率空间域准P波正演模拟

有限差分方法是波场数值模拟的一个重要方法,时间域有限差分计算方法因按时间片递推计算,每个时间片的舍入误差会累积到下一片中,当时间片较多,最终会导致累积误差太大.而频率域计算是按频率片对空间网格进行整体求解方程组,其计算误差分配到了每个网格点上,并且各个频率片之间是独立计算的,因此不存在累计误差,而且在频率-空间域更易于...     

Generalized migration in frequency-wavenumber domain (MGF-K) in anisotropic media

In this paper, the background of MGF-K migration in dual domain (wavenumber-frequency K-F and space-time) in anisotropic media is presented. Algorithms for poststack (zero-offset) and prestack migration are based on downward extrapolation of acoustic wavefield by shift-phase with correction filter for lateral variability of medium’s parameters. In anisotropic media, the vertical wavenumber was determined from full elastic wavefield equations for two dimensional (2D) tilted transverse isotropy (TTI) model. The method was tested on a synthetic wavefield for TTI anticlinal model (zero-offset section) and on strongly inhomogeneous vertical transverse isotropy (VTI) Marmousi model. In both cases, the proper imaging of assumed media was obtained.     

Seismic waves in cracked two-phase anisotropic medium with fourfold rotation symmetry

IntroductionBasedontheresearchresultSofmanyintemalandextemalseismologiSts,mediuminseis-mogeniczonemaybethecrackedinordermediumsaturatedwithfluid(airorwater),namely,two-phaseanisotropicmedium.Inseismicprospecting,oilorgasreservoirrockmayaPproachtothecrackedtwo-phaseanisotropicmediumtoo.Therefore,thefiJrtherstudyonpropagaiontheoryofseismicwavesintWo-phaseanisotropicmediumpossessestheoreticalsignificanceandaPPlicationprospects-Inthispaper,onthebasisofthetheoryofseismicwavesinanuid-sbedporousso…     

交错网格与旋转交错网格对VTI介质波场分离的影响分析

利用交错网格有限差分和旋转交错网格有限差分进行各向异性介质弹性波场数值模拟时, 质点振动速度分量与应力张量的网格节点定义方式均不相同, 从而对各向异性波场分离效果产生不同的影响. 针对这一问题, 本文以具有垂直对称轴的横向各向同性(VTI)介质的波场分离为例, 首先分析了两种网格的参数定义方式以及VTI介质波场的分离过程; 其次, 详细研究和分析了这两种网格的参数定义方式对各向异性介质波场分离的影响, 并依据波前面连续性以及波场分离效果等方面, 通过数值模拟实验对该影响进行分析验证. 结果表明, 旋转交错网格的参数定义方式更有利于进行各向异性介质波场数值模拟和波场分离.      

裂缝诱导各向异性双孔隙介质波场传播特征

基于裂缝诱导各向异性和双相介质理论,对裂缝诱导的具有水平对称轴的横向各向同性（HTI）双孔隙介质的本构关系进行了研究,与等效连续介质模型相结合,综合考虑裂缝系统和基质孔隙系统的两种孔隙度和两种渗透率参数,得到裂缝诱导HTI双孔隙介质的等效孔隙度和等效渗透率,进而得到介质的运动平衡方程;并进一步推导出介质的一阶速度-应力方程.采用交错网格高阶有限差分法对模型进行了数值模拟,结果揭示了介质中两套系统的存在对其波场传播特征的影响,为进一步研究实际地球介质的波场特征奠定了基础.