一阶速度-应力Biot双相各向同性介质弹性波波场分离数值模拟

提出一种等价的一阶双曲型速度-应力Biot双相各向同性介质弹性波波动方程,以实现双相介质混合波场中纯快慢纵波和纯横波波场分离的问题.应用散度和旋度理论证明双相介质等价方程波场分离的可行性,采用高阶交错网格有限差分法构建高精度正演算子,推导其PML吸收边界条件和稳定性条件,并对均匀双相介质和层状非均匀双相介质模型进行数值模拟试验,准确得到固流相混合弹性波场、被完全分离的纯纵横波波场,同时边界吸收效果良好,数值模拟精度较高.计算结果还表明,固流相中的快慢纵波相互伴生因而无法实现分离,且归属于纯纵波波场,流相慢纵波能量比固相慢纵波能量强,这对认识双相介质弹性波的传播规律以及完善双相介质理论具有重要意义.

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

为了便于研究双相介质固流相混合弹性波场中纵横波波场的传播规律,提出了基于交错网格的Biot双相各向同性介质弹性波动方程高精度波场分离正演数值模拟方法。采用高阶交错网格有限差分法来构建一阶双曲型双相各向同性介质弹性波动方程正演算子实现波场正演,并在每一步递推过程中,分别计算出固相和流相分量相应的散度场(纯纵波场)和旋度场(纯横波场)。以均匀介质和层状介质模型为例进行波场分离数值模拟试验,计算结果表明,该方法能够从固流相的任一分量混合弹性波场中准确地分离出纯纵波波场和纯横波波场,具有计算精度和效率高的特点,同时指出快慢纵波相互伴生因而无法实现分离。

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

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

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

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

各向异性双相介质中多极源声波测井理论研究

将各向异性地层模拟为横向各向同性双相介质 ,使其对称主轴与井孔中心轴线相垂直 .本文针对井孔模型下的多极源声波测井问题提出了一种解析的摄动方法 .将各向异性双相介质看作是由对应的各向同性双相介质加上摄动部分而组成 ,并引入 5个关于弹性模量偏离各向同性时的摄动量 ,对对称主轴与井轴垂直的各向异性双相介质的多极源声波测井问题进行了系统的理论研究和分析 ,严格求解了井孔内外的零级和一级近似声波场 .发现在n级多极源激发下 ,井孔中不仅出现n级多极声场 。

双相TI介质中弹性波交错网格高阶有限差分法数值模拟

应用双相介质波动方程,推导了双相横向各向同性介质(TI)中波动方程的有限差分格式,对双相TI介质中弹性波有限差分数值进行了模拟。结果表明,弹性波在双相TI介质中传播时,除了存在常规的快纵波(qP1)和横波以外,还存在慢纵波(qP2)。并且慢纵波的速度明显小于快纵波,而且受耗散系数的影响衰减地很快,所以在实际中很难观测到慢纵波。快纵波在固相和流相中相位相同,而慢纵波在固相和流相中的相位相反。慢纵波在流相中振幅大,而在固相中的振幅较小。

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

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