双相介质中储层波场特征数值模拟研究

根据双相介质弹性波方程,并对弹性波方程进行分离得到了膨胀波方程。利用高阶差分对膨胀波方程进行数值模拟,边界处理采用完全匹配层的吸收边界条件,算法实现了任意高阶差分和自动加载边界条件的纵波方程正演模拟。计算结果表明,在双相各向同性介质中存在快纵波、慢纵波,两种纵波有明显的区别,快纵波的速度远大于慢纵波的速度。在分界面上,快纵波要产生透射快纵波、透射转换慢纵波,反射快纵波、反射转换慢纵波。慢纵波具有很强的衰减性,耗散系数越小,慢纵波越明显,反之亦然。该算法的高阶差分形式可以显著地降低数值频散,有效提高地震波正演计算的精度,适应性好,操作简便灵活。     

波数域等价波动方程三维弹性波波场分离

将基于传统弹性波波动方程正演模拟或者反向外推得到的波场作为纵横波耦合的混合波场,利用亥姆霍兹分解方法对混合波场进行分离,可得到纯纵波波场和纯横波波场,但是用这种方法分离后的波场和原波场相比产生了振幅畸变。为了消除这一现象,笔者从亥姆霍兹分解方程入手,利用纵波波场旋度为零、横波波场散度为零的特性,在波数域将传统波动方程分解为无旋部分和无散部分,得到波动方程的一种等价表示。数值实例表明,该方法既能获得混合波场、纯纵波波场和纯横波波场,又能使分离前后的波场振幅无畸变。     

饱和多孔介质土动力学理论与数值解法

总结和评价了饱和多孔介质土动力学理论与数值解法在国内外的研究成果 ,着重讨论了饱和多孔介质土动力学的基本方程、饱和土中弹性波的传播特性以及饱和土体动力分析的时域数值解及地震反应等方面的研究现状 ,并指出今后的研究方向 .     

半空间各向同性弹性饱和多孔介质非轴对称动力响应分析

应用弹性饱和多孔介质非轴对称Biot方程的解,研究了半空间饱和多孔介质非轴对称动力响应问题。讨论了半空间饱和多孔介质在边界竖直方向和水平方向,存在非轴对称简谐荷载作用下的边界位移分布,给出边界竖向位移的数值计算结果,并对饱和多孔介质非轴对称动力响应特性进行了分析。     

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

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

多孔介质声波传播

基于二维有限差分给出了Ｂｉｏｔ声波方程的一种数值解方法，模拟了声波在多孔介质中的传播。研究发现，骨架和粘滞性孔隙流体中，同时存在两种纵波，即快纵波和慢纵波，两种纵波骨架振动势函数与流体相对骨架振动势函数相不同，幅度也不同；但具有相同的传播速度，粘滞性和渗透率主要影响能量传播；孔隙度和骨架弹性参数主要影响波速；慢纵波和界面波的存在是地震波能量损失的重要原因。     

多孔介质声波传播

基于二维有限差分给出了Ｂｉｏｔ声波方程的一种数值解方法，模拟了声波在多孔介质中的传播，研究发现，骨架和粘滞性孔隙流体中，同时存在两种纵波，即快纵波和慢纵波，两种纵波骨架振动势流体相对骨架振动势函数相不同，幅度也不同，但具有相同的传播速度，粘滞性和渗透率主要影响能量传播，孔隙度和骨架弹性参数主要影响波速，慢纵波和界面波的存在是地震波能量损失的重要原因。     

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

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

求解弹性饱和多孔介质半空间非轴对称动力问题的传递矩阵法

应用Fourier展开和Hankel变换 ,求解了弹性饱和多孔介质非轴对称Biot波动方程 ,得到半空间问题的一般解 .引入了位移分量和应力分量组合量 ,建立了这些组合量Hankel变换之间的关系 .用传递矩阵方法处理了弹性饱和多孔介质半空间非轴对称动力问题 ,把边界上的应力和位移与介质内的应力和位移用传递矩阵联系起来 ,这种处理方法特别适用于数值分析 .     

储层模型正演模拟及AVO特征研究

为了在多波地震资料解释过程中解决多波信息识别难的问题,基于弹性波正演技术,针对某油田勘探区块的实际储层地质模型,采用高阶交错网格差分技术正演模拟非零偏炮集多波地震记录,并经多波资料处理,得到纵波和转换波CMP和CCP叠前含AVA的动校正道集及叠加剖面。在已知模型的AVA计算、模型层位参数与正演纵、横波波场关系的基础上,对比研究了正演记录中各储层的纵波和转换波对应关系。对比分析结果表明,数值模拟结果与AVA分析结果吻合,同一储层的纵波和转换波波场响应有明显差异。