煤储层双相EDA介质的地震波场模拟

基于二维-二分量弹性波方程交错网格高阶有限差分法数值模拟,对各向同性介质模型、高裂隙密度含干裂隙和饱和水裂隙煤层各向异性介质和双相介质的二维层状模型进行了弹性波场数值模拟.结果表明:双相EDA介质中存在第二类纵波(慢纵波).各向异性理论和双相介质理论比单相理论要好,单相理论存在衰减问题.各向异性和双相理论在含气情况下比...     

双相介质中地震波的频率-空间域数值模拟

本文利用优化的25点频率-空间域有限差分算法对基于BISQ模型双相各向同性介质中的地震波进行了数值模拟.通过与经典的Biot模型理论模拟结果进行对比,分析了Biot流动（宏观流体流动）和Squirt流动（微观流体流动）耦合作用对地震波在孔隙介质中传播特性的影响.数值模拟在地震频段进行,结果显示：在理想相界和黏滞相界情况下,Squirt流动机制都比Biot流动机制产生了更大的速度频散和能量衰减.其中,在Biot流动和Squirt流动耦合作用下的快P波的速度和振幅小于仅考虑Biot流动影响下快P波速度和振幅,而且慢P波的衰减也更加强烈.本文还研究了地震波在双层双相各向同性介质分界面处的反射和透射特征,双相介质中波的反射与透射现象类似于单相介质的情况.模拟结果表明,利用优化25点频率-空间域有限差分法模拟双相孔隙介质中的地震波场是可行的,这为开展双相孔隙介质全波形反演问题的研究提供了可能.     

基于BISQ模型双相各向同性介质弹性波传播的频率-空间域有限差分模拟

BISQ模型同时包含了合流体孔隙介质中Biot流动和喷射流动两种重要的力学机制.基于BISQ模型的双相各向同性介质弹性波波动方程是一个复变系数偏微分方程组.本文率先建立了该方程的25点频率-空间域有限差分数值解法,在理想相界和黏滞相界情况下,对Biot流动和喷射流动共同作用下的双相各向同性介质中的波传播进行了数值模拟,通过与仅受Biot流动作用下的波场模拟结果的对比,分析了两种机制耦合作用对弹性波传播的影响.同时,本文也研究了波在双层双相各向同性介质分界面处的反射和透射特征.数值模拟结果表明:在Biot流动和喷射流动耦合作用下,双相介质中传播的快P波的速度和振幅都小于仅考虑Biot流动的双相介质中的快P波,且慢P波的衰减也更为强烈,而S波的波速和振幅则无明显差异.这表明局部喷射流动对P波的衰减和频散具有重要影响,而对S波的影响较小;慢P波的强烈衰减使得其在波场快照中无法被观测到,双层双相介质中的波传播现象类似于单相介质的情况.同时本文的研究结果也表明,频率-空间域有限差分法在基于BISQ模型的双相介质中波传播数值模拟中的正确性和有效性,为开展孔隙弹性介质全波形反演问题的研究提供了研究基础.     

基于分数阶时间导数常Q黏弹本构关系的含黏滞流体双相VTI介质中波场数值模拟

分数阶微分算子具有描述历史依赖性和全域相关性的特质,本文利用这种特质描述双相介质固体骨架的黏弹性特征.基于Kjartansson常Q理论将含有分数阶时间导数的黏弹固体骨架各向异性本构关系与双相介质理论有机地结合起来,并引入流变学本构关系描述孔隙流体的黏滞性力学行为,提出一种新的基于分数阶时间导数常Q黏弹本构关系的含黏滞流体双相VTI模型.推导了相应的时间域波传播方程,然后对该方程进行了数值模拟.对整数阶导数采用高阶交错网格有限差分算法,对分数阶时间导数采用短时记忆中心差分算法,进行了不同相界、不同品质因子组及双层地质结构情况下该类介质中波场的数值模拟与特征分析.模拟结果表明：将含有分数阶时间导数的常Q黏弹固体骨架各向异性本构关系及孔隙流体的黏滞性本构关系引入双相介质理论是可行的,二者的结合能更好地反映地下介质的黏弹性特征,对于进一步认识波在黏弹各向异性孔隙介质中的传播机理具有重要意义,为反演和重构地下油气储层和结构奠定正演理论基础.

     

双相各向同性介质弹性波高精度波场分离数值模拟方法

为了便于研究双相介质固流相混合弹性波场中纵横波波场的传播规律,提出了基于交错网格的Biot双相各向同性介质弹性波动方程高精度波场分离正演数值模拟方法.采用高阶交错网格有限差分法来构建一阶双曲型双相各向同性介质弹性波动方程正演算子实现波场正演,并在每一步递推过程中,分别计算出同相和流相分量相应的散度场(纯纵波场)和旋度场...     

常Q衰减介质分数阶波动方程优化有限差分模拟

本文基于Kjartansson常Q模型理论,推导了常Q衰减介质中黏声波和黏弹性波的速度-应力方程,并采用基于二项式窗函数的优化交错网格有限差分方法进行了数值模拟,同时引入不分裂的复频移卷积完全匹配层(CPML)吸收边界条件,以消除边界反射.使用基于自适应时间步长记忆方法的中心差分近似时间分数阶导数,与常用的短时记忆方法相比,提高了波动方程的离散化精度和计算效率.通过对比均匀模型下声波的数值解与解析解,验证了算法的精确性,并进一步分析了不同品质因子下地震波的频散及衰减特征.对BP盐丘模型的数值模拟结果可以较好地反映本文数值方法对复杂介质的适应性及频散压制效果.      

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

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

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

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

各向同性介质中射线理论与有限差分地震波场模拟方法比较研究

地震波场数值模拟不仅是研究复杂地区地震资料采集、处理和解释的有效辅助手段,而且是研究地球深部精细构造和地球深部探测的有效工具.射线理论和波动方程理论是地震波场数值模拟的理论基础.射线理论主要刻画地震波在介质中传播的走时场、地震射线等运动学属性;波动方程理论通过求解波动方程来描述地震波在介质中传播的弹性动力学响应(能量衰减、相位特征、偏振属性、以及全波形等).基于波动方程理论的波场数值模拟由于能够引入丰富的波场信息,使得人们对不同介质中地震波的传播过程有了较全面的了解.本文以二维层状均匀介质模型为例,通过射线追踪法和交错网格有限差分法模拟得到的波场快照图、单炮地震记录剖面、合成理论地震图的分析比较,不但对地震波在各向同性层状均匀介质中的传播规律和特点有了深刻的理解和认识,同时又可以相互验证两种不同方法的正确性和有效性.     

基于旋转交错网格伪谱法和L-S理论的热弹性介质波场模拟与特征分析

基于L-S热弹性理论,采用旋转交错网格伪谱法,实现了均匀各向同性介质的一阶速度-应力-温度微分方程组的数值求解和波场模拟.其中,用时间分裂法解决方程组的刚性问题,用旋转交错伪谱算子计算空间一阶导数,用中心伪谱算子计算空间二阶导数;对于热导率变化比较大的双层介质模型和参考温度随深度按梯度分布的模型,用Crank-Nicolson显式方法取代旋转交错网格伪谱法进行计算;讨论了热耦合波场性质和传播规律,对比了常规伪谱法、交错网格伪谱法和旋转交错网格伪谱法热耦合波场模拟效果.数值模拟结果表明:基于L-S热弹性理论,用时间分裂法结合旋转交错网格伪谱法对均匀各向同性介质的波场计算,能够得到稳定的、高精度的模拟结果,但是在热导率变化剧烈的情况下不能用大时间步长进行求解,而且在参考温度不均匀分布的情况下算法不稳定.本文在网格剖分方式与数值算法的优化组合应用方面进行了探索,为将这些方法推广到孔隙热弹性、热黏弹性和各向异性研究奠定了基础.     

求解双相和黏弹性介质波传播方程的间断有限元方法及其波场模拟

间断有限元（Discontinuous Galerkin：DG）方法具有低数值频散、网格剖分灵活、能模拟地震波在复杂介质中传播等优点.因此,本文将一种新的DG方法推广到双相和黏弹性等复杂介质的地震波场模拟,发展了求解Biot弹性波方程和D'Alembert介质波动方程的DG方法.首先通过引入辅助变量将Biot双相介质弹性波方程和D'Alembert介质波动方程转化为关于时间-空间的一阶偏微分方程组,然后对该方程组进行DG空间离散,得到半离散化的常微分方程组.最后,对此常微分方程组,应用加权的Runge-Kutta格式进行时间推进计算.数值结果表明,DG方法可以有效地求解Biot双相介质弹性波方程和D'Alembert介质波动方程,并能很好地压制因离散求解波动方程而产生的数值频散,获得清晰的各种地震波震相.

     

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

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

黏弹性介质VSP记录模拟及在估算Q值研究中应用

采用标准线性固体模型,本文建立了黏弹性介质完全匹配层吸收边界的高阶速度-应力交错网格有限差分算法,并对黏弹性介质中的地震波传播进行了数值模拟.基于黏弹性波动方程正演模拟提供的零偏VSP全波场数据,本文进行了质心频移法计算Q值的反演分析.结果表明,反射波、转换波及短程多次波对频谱的影响较大,对Q值反演造成一定误差.本文的结论为实际零偏VSP资料估算地层Q值提供了有益的指导.     

起伏地表下基于改进BISQ模型双相介质中曲线网格有限差分法波场模拟

本文以基于改进BISQ模型的二维双相各向同性介质一阶速度-应力方程为基础,推导出了曲线坐标系下对应的方程,然后采用低频散、低耗散的同位网格MacCormack有限差分法来离散方程,并采用紧致的单边MacCormack差分格式结合牵引力镜像法来施加自由地表边界条件,实现了地震波场数值模拟.曲线网格有限差分法采用贴体网格来描述自由表面,地表的网格线紧贴地形,避免了台阶近似造成的数值散射.数值模拟结果表明,在双相介质起伏自由地表和分界面处,各类波型复杂的反射透射规律可以清晰展现,曲线网格有限差分法可以精确地解决地震波在含起伏地表的双相各向同性介质中的传播问题.

     

复杂介质多尺度变网格地震波数值模拟

在应用有限差分法地震波场数值模拟过程中,模拟精度和计算效率是关键问题之一,特别是对于速度变化剧烈的含小尺度地质体的数值模拟尤为重要.为了既能精细刻画介质的局部结构,又能保证模拟的时效性,本文在传统交错网格有限差分算法的基础上,将变网格有限差分方法引入到交错网格高阶差分数值模拟中,对交错网格空间算法进行了改进,并避免了因插值因素影响模拟精度和计算效率.本文采用完全匹配层吸收边界条件实现对可变网格数值模拟的边界处理,根据完全匹配层吸收边界原理和参量的分裂思想,在不同的方向分别加相应的阻尼衰减因子,解决边界的反射问题.此外还分析了数值频散对于模拟结果的影响.文中分别设计了地堑和生物礁两个介质模型对可变交错网格算法的有效性和稳定性进行了验证.数值模拟结果表明,可变交错网格有限差分算法数值模拟能够对地下介质物性的空间变化进行准确刻画,进一步增强数值模拟对复杂介质的适应性,同时为地震数据的波场成像、纵横波联合解释等提供可靠依据.     

基于横向各向同性BISQ方程的弹性波传播数值模拟

Biot流动和喷射流动是含流体多孔隙介质中流体流动的两种重要力学机制. 近年来，利用同时处理这两种力学机制的BISQ(Biot-Squirt)模型，弹性波衰减和频散的问题已被广泛研究；然而基于BISQ方程的波场数值模拟尚未见到公开的报道．本文从BISQ方程出发，利用交错网格方法对横向各向同性孔隙介质中不同频率和相界情况，以及双层介质中的弹性波传播进行数值模拟，研究了在同时考虑两种流动机制作用情况下地震波和声波的传播特性及传播过程中出现的各种波动现象．      

基于高阶广义标准线性体模型的三维黏弹性介质弹性波正演模拟

吸收衰减是地震波在实际地球介质中传播的固有特征.在实际应用中,通常假设表征吸收衰减特征的品质因子Q在地震频带范围内不随频率变化.高阶广义流变模型能够在时间域内精确的表征品质因子Q不随频率变化的特征,为黏弹性介质波动方程精细模拟奠定了基础.基于广义标准线性体模型理论,采用最小二乘拟合方法对Q值不随频率变化特征进行拟合,分析了不同阶次广义标准线性体模型对黏弹性介质Q值特征的拟合程度,在权衡计算精度和三维计算量的基础上,确定了五阶广义标准线性体模型并建立了相应的三维黏弹性波的速度-应力方程,结合CFS-PML边界条件开展了高精度三维黏弹性波正演模拟.通过均匀介质正演模拟,验证算法的正确性,明确了地震波的传播时的吸收衰减特征,对三维盐丘模型进行数值模拟,表明了五阶广义标准线性体可以精确的模拟黏弹性介质地震波吸收衰减特征.     

基于拟微分算子的横向各向同性介质黏声波方程

由于地层的各向异性与黏滞性,地震波在传播过程中表现出方向上的各向异性和振幅值的吸收衰减特征.若简单地将地下介质视为各向同性,实际地质结构的地震波响应可能会被曲解,甚至忽略了有效信息.为了得到深部储层的高精度成像,本文将标准线性固体模型(SLS)扩展到各向异性介质中,并结合各向异性拟声波方程,推导出时间空间域的各向异性介质黏声波方程组,该方程中用伪微分算子来表述地震波的衰减特征.在逆时偏移的过程中,为了解决高频不稳定问题,文中通过引入规则化算子来构建稳定的逆时传播算子.模型的正演模拟表明该方程既能准确的描述各向异性介质下的标量波的传播规律,又能体现地层的吸收衰减效应.     

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

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

基于Biot-Squirt方程的波场模拟

Biot流动和喷射流动是含流体多孔隙介质中流体流动的两种重要力学机制,对地震波和声波的传播均产生重要影响. Dvorkin和Nur提出了同时包含Biot流动和喷射流动力学机制的统一的BISQ(Biot-Squirt)模型,基于这一模型,尽管有关弹性波在多孔隙介质中的衰减和频散问题已被广泛研究,然而,基于BISQ波传播方程的波场数值模拟至今仍未见报道. 本文从同时包含两种力学机制的孔隙弹性波方程出发,利用FCT有限差分法对含流体孔隙各向同性介质中的地震波和声波进行了数值模拟,并与基于Biot流动的Biot理论之模拟结果进行比较. 数值模拟结果表明:同时包含Biot流动和喷射流动影响的地震波和声波速度比仅包含Biot流动作用的地震波和声波速度慢,慢P波的衰减比根据Biot理论模拟的慢P波衰减更强.