三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟

三维数值模拟是研究各向异性介质中复杂弹性波波场特征和传播规律的重要手段。根据Taylor展开式,推导出了交错网格一阶空间导数的任意偶数阶精度有限差分近似式和相应的差分系数计算式,给出了三维各向异性介质中弹性波一阶双曲型应力速度方程交错网格任意偶数阶精度差分格式和稳定性条件,并推导出了三维各向异性介质PML法吸收层边界条件公式和相应的交错网格差分格式。对方位各向异性介质模型和正交各向异性介质模型中弹性波的传播进行的三维数值模拟结果表明,弹性波在三维各向异性介质中传播时存在拟P波、拟SV波和拟SH波,并出现横波分裂、横波分裂盲区、波面三分叉等特殊现象,另外,弹性波场在空间是变化的,其拟P波、拟SV波和拟SH波的耦合关系比较复杂。     

三维双相各向异性介质弹性波方程交错网格高阶有限差分法模拟

基于Biot理论，给出了三维双相各向异性介质应力一速度弹性波方程交错网格任意偶阶精度有限差分解法，并对三维双相各向异性（TIH）介质中弹性波场进行了模拟。结果表明，弹性波在三维双相各向异性介质中传播时存在快纵波qP1、快横波qS1、慢横波qS2和慢纵波qP2，并清楚地观测到。T横波分裂、横波分裂盲点、波面三分叉等特殊现象。快纵波在固相和流相中的相位相同；而慢纵波在固相和流相中的相位相反，即慢纵波在流相中振幅大，而在固相中的振幅较小。另外，三维双相各向异性介质中弹性波场的快纵波和快横波的耦合关系、波的类型、能量分布和相位等都是在三维空间中变化的。     

三维双相各向异性介质弹性波方程交错网格高阶有限差分法模拟

基于Biot理论,给出了三维双相各向异性介质应力-速度弹性波方程交错网格任意偶阶精度有限差分解法,并对三维双相各向异性(TIH)介质中弹性波场进行了模拟.结果表明,弹性波在三维双相各向异性介质中传播时存在快纵波qP1、快横波qS1、慢横波qS2和慢纵波qP2,并清楚地观测到了横波分裂、横波分裂盲点、波面三分叉等特殊现象.快纵波在固相和流相中的相位相同;而慢纵波在固相和流相中的相位相反,即慢纵波在流相中振幅大,而在固相中的振幅较小.另外,三维双相各向异性介质中弹性波场的快纵波和快横波的耦合关系、波的类型、能量分布和相位等都是在三维空间中变化的.     

非均匀介质地震波传播交错网格高阶有限差分法模拟

采用常规的二阶声波方程有限差分方法对于非均匀介质进行了数值模拟时，其数值模拟精度较低。而采用一阶双曲型标量波动方程，则无须对介质的弹性常数进行空间求导。根据Taylor级数展开式，推导出了交错网格一阶空间导数的任意偶数阶精度展开式和相应差分系数计算式以及一阶双曲型标量波动方程交错网格任意偶数阶精度差分格式，并给出了该差分算法的稳定性条件。用该差分算法对均匀介质模型、非均匀介质模型和Mannousi模型进行了数值模拟试验，并与伪谱法进行了对比。结果表明，一阶双曲型标量波动方程交错网格高阶差分法的模拟精度与伪谱法的精度非常接近，计算效率高，且适合于模拟非均匀介质、复杂构造和复杂地质体的地震波场。     

非均匀介质地震波传播交错网格高阶有限差分法模拟

采用常规的二阶声波方程有限差分方法对于非均匀介质进行了数值模拟时 ,其数值模拟精度较低。而采用一阶双曲型标量波动方程 ,则无须对介质的弹性常数进行空间求导。根据Taylor级数展开式 ,推导出了交错网格一阶空间导数的任意偶数阶精度展开式和相应差分系数计算式以及一阶双曲型标量波动方程交错网格任意偶数阶精度差分格式 ,并给出了该差分算法的稳定性条件。用该差分算法对均匀介质模型、非均匀介质模型和Marmousi模型进行了数值模拟试验 ,并与伪谱法进行了对比。结果表明 ,一阶双曲型标量波动方程交错网格高阶差分法的模拟精度与伪谱法的精度非常接近 ,计算效率高 ,且适合于模拟非均匀介质、复杂构造和复杂地质体的地震波场     

频率域弹性波场模拟中高阶有限差分法的精度对比与改进策略

推导基于二维交错网格的二阶和四阶频率域弹性波有限差分算子,并结合最优化差分系数和质量加权平均方法压制数值各向异性。基于均匀各向同性介质比较不同差分算子的模拟精度和频散关系,并详细比较二阶算子和四阶算子的精度和误差量级。结果表明,当采用交错网格最优化四阶差分算子,单位横波波长内含3个及以上网格点时,可将模拟误差控制在2%以内。同时,相对于常规网格,基于交错网格的差分算子可用于对含流体介质模型的模拟。     

任意偶数阶精度交错网格声波方程数值模拟

以波动理论为理论基础,利用泰勒级数展开推导出了一阶应力-速度声波方程组的空间任意偶数阶精度交错网格差分格式。选用衰减边界条件,进行了边界效果对比及差分精度选取测试。结果表明差分精度越高,频散越弱,数值模拟的效果越好。应用空间八阶、时间二阶精度差分格式,实现了各向同性介质模型的声波方程数值模拟。正演模拟结果可以清晰地反映出声波波场中波传播的运动学和动力学特征,可推广应用于三维声波高精度波场正演模拟中。     

基于波场分离的双相各向同性介质正演模拟

基于双相各向同性介质一阶速度-应力方程,利用旋度和散度理论的波场分离方法进行纵横波场解耦,借助高阶交错网格有限差分技术,推导出了其任意偶数阶精度交错网格有限差分格式,最后运用空间八阶时间二阶差分精度及PML边界条件,完成了双相各向同性介质中地震波场正演模拟。结果表明,该方法可准确获得双相各向同性介质固、流相中的高精度纯纵波波场和纯横波波场,有效描述双相各向同性介质中的波场传播规律。同时得出:双相各向同性介质中快、慢纵波相互耦合,无法单独分离;也指出了分界面处波的传播和转换关系,是油气勘探的重要理论依据。     

任意裂隙方位双相HTI介质的地震波方程及正演模拟

基于BISQ方程的高频极限,导出了适用于任意方位角的双相HTI介质波动方程,推导了该方程在交错网格空间中求解的高阶有限差分格式和对应的完全匹配层(PML)吸收边界条件,实现了该类介质的地震波场正演模拟。模拟结果表明,该方法能准确模拟地震波在双相HTI介质中的传播过程,得到高精度的波场快照和合成记录。分析得出:在双相各向异性介质中,弹性波的传播仍表现出各向异性特征;孔隙度、粘滞系数和固液耦合密度主要影响慢纵波的振幅和形态。     

方位各向异性介质的多尺度有限差分法波场模拟

采用以裂缝走向为方位角的方位各向异性介质模型可以较好地描述裂缝性油气藏的实际情况,对弹性波在此介质中的传播过程进行准确的数值模拟有助于提高油气开采的准确程度。该文采用紧支集正交小波基对空间域进行多尺度离散,采用二阶精度有限差分算子对时域离散,推导得到了多尺度有限差分方法正演模拟的递推公式,并实现了相应的波传过程数值模拟。数值结果准确地反映了方位各向异性介质中波场的变化过程,可以清晰地观察到横波分裂和方位特征差异等现象。