双相组分孔隙介质波场数值模拟

孔隙介质理论的研究以往常用试验分析方法,本文详细介绍了双相组分孔隙介质理论,基于宏观上等效的一系列体积分数加权平均组分分量方程,得到双相组分孔隙介质整体弹性矩阵与所有组分弹性矩阵之间的体积分数及耦合系数加权关系,得出各弹性参数的表达式;结合弹性波动力学中的Cauchy、Navier和本构三个方程,得到弹性波在等效后的组分孔隙介质中传播满足的波动方程,应用交错网格有限差分法求解该波动方程;采用不同孔隙度双相组分孔隙介质模型波场数值模拟,精确得到了混合波场;总结了双相组分型弹性孔隙流体介质中地震波传播的特点和规律。     

双相介质波动方程正演模拟研究

本文在详细地介绍了双相介质和小波分析的有关理论的基础之上,对Biot方程进行了差分数值模拟正反演研究。利用有限差分法对双相介质Biot理论的波动方程进行了正演模拟,该方法对地震勘探中的双相介质方程的正演数值模拟取得较好的结果,为今后继续从事双相介质的工作打下了一定的基础。     

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

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

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

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

随机介质模型波场的有效统计特征

通过弹性波方程的交错网格有限差分正演,模拟了平面地震波在单层矩形介质模型中的传播及其自激自收时间记录.将理论记录剖面(垂直和水平分量)分割成5个不同的时间区段.在5个不同的时间区段上,分别计算剖面的统计特征:横向中心波数、纵向中心频率和波场能量相对值.这样,对应每一个弹性随机介质模型,均可计算得到15个不同的波场特征量.进一步引入了相对标准差的概念,并由此得到了若干有效波场特征量.通过正演模拟,依次研究了当介质的自相关长度变化、介质的自相关长度和粗糙度因子变化及介质的自相关长度和扰动标准差变化时的波场特征.讨论了介质模型发生变化时,对应的波场特征量的变化特点.     

PSTM四层介质理论探测场的数值模拟

使用最新光子扫描隧道显微镜探测场的理论计算公式，采用比较符合实际探测的理论模型，对ＰＳＴＭ探测场进行了计算机数据模拟，对探测场分布随扫描控制参量的变化规律进行了系统的研究，     

改进BISQ模型的双相介质地震波场数值模拟及频散校正

从改进BlSQ模型的双相介质所对应的速度-应力运动方程出发,构建2×2N阶交错网格有限差分模拟算法:同时,为压制模拟过程中的数值频散现象,采用通量校正传输(FCT)技术获得带FCT修正的交错网格有限差分模拟算法,对各向同性双相介质进行地震波场数值模拟.研究结果表明:(1)波场存在快纵波、慢纵波和横波等波场特征,并与模型的理论响应相符,说明交错网格模拟算法具有正确性和可行性;(2)采用FCT修正的交错网模拟算法能够有效压制数值频散,并保留真实的波场特征.     

TI介质地震波场数值模拟边界条件处理

地震波场数值模拟研究中,边界条件处理非常重要,它直接影响模拟地震数据的精度.针对TI(Transversely Isotropic)介质一阶速度—应力弹性波动方程组,利用时间域分裂法给出了PML(Perfectly matched layer,完全匹配层)吸收边界条件,其基本思想是在所研究的区域外引入吸收层,地震波从有效计算区域传播到吸收层时不发生任何反射,且在吸收层内按传播距离呈指数衰减,不发生反射.利用交错网格高阶差分技术对PML吸收边界条件求解,对比不同边界条件数值模拟结果表明,PML吸收边界条件较应用广泛的Clayton吸收边界有更好的吸收效果,能更好地满足地震数值模拟的需要.     

Biot方程差分数值模拟反演研究

本文对Biot方程进行了差分数值模拟反演研究。利用小波多尺度方法结合时卷正则化方法对于四层地层模型的孔隙率进行了反演计算,计算结果表明小波多尺度反演方法是一种快速的非线性的全局反演方法。这种方法结果稳定、有效、简洁、实用。通过增加15%噪声证实,上述方法的抗噪性也良好。     

黄土塬斜坡带地震波激发数值模拟与波场特征

黄土塬斜坡带的地震资料品质最差,研究程度极低。采用先进的交错网格、谱元法和高阶有限差分波场模拟算法,对斜坡带地震波激发进行了数值模拟,准确地反映了地震波场特征,从理论上优选出地震波激发参数。经试验提出一套超多井组合地震波激发技术,取得了高精度的三维地震资料。     

热采井筒应力的数值模拟分析

将注汽热采井井筒温度场传输到井筒应力模型中,用有限元软件ANSYS进行了计算。计算结果表明,使用常温套管弹性模量对套管极限温度载荷进行计算,所得极限温度比实际极限温度要低,因此能提高设计结果的安全性。增大水泥高温热膨胀系数有利于保护套管和水泥环。     

热采井筒应力的数值模拟分析

将注汽热采井井筒温度场传输到井筒应力模型中 ,用有限元软件ANSYS进行了计算。计算结果表明 ,使用常温套管弹性模量对套管极限温度载荷进行计算 ,所得极限温度比实际极限温度要低 ,因此能提高设计结果的安全性。增大水泥高温热膨胀系数有利于保护套管和水泥环。     

双重介质油藏数值模拟并行算法研究

裂缝型碳酸盐岩油藏非均质性极强,数值模拟运算时间长,迫切需要一种切实可行的解决方法。介绍了区域分解法在数值模拟中的应用方法,提供了油藏数值模拟串行算法并行化的一般策略与方法。基于已开发的双重介质油藏数值模拟串行算法及建立的Windows平台PC集群网络并行计算环境,利用消息传递界面（MPI库）实现了并行过程中所有消息的传递功能。通过对国内某实际碳酸盐岩油田的模拟表明,并行后的算法是成功的,2个PC加速比高达1.86,4个PC达到3.46,8个PC达到6.89。通过模拟还发现,并行数值模拟时,存在最佳的子区域规模和数目。该算法解决了大规模碳酸盐岩油藏数值模拟时间长的问题,具有实际的应用价值。     

论多孔介质中流体流动问题的数值模拟方法

简要介绍了多孔介质中流动问题的数值模拟方法及其发展现状 ,对各种方法的优劣及适用范围进行了评价。对多孔介质纯流体耦合流动区域的数值模拟方法及发展现状进行了综合对比。多孔介质模型方法可以用于多孔介质中流体流动问题的数值模拟 ,并且非常适于进行多孔介质纯流体耦合流动区域的数值模拟。     

论多孔介质中流体流动问题的数值模拟方法

简要介绍了多孔介质中流动问题的数值模拟方法及其发展现状,对各种方法的优劣及适用范围进行了评价。对多孔介质－纯液体耦合流动区域的数值模拟方法及发展现状进行了综合对比。多孔介质模型方法可以用于多孔介质中液体流动问题的数值模拟,并且非常适于进行多孔介质－纯液体耦合流动区域的数值模拟。     

多孔介质两相流的守恒迎风有限元法

本文讨论油藏模拟中描述多孔介质两相不可压非混溶流动的偏微分方程组。压力方程用Galerkin方法,而对流占优的饱和度方程用一类部分迎风有限元法,数值解满足离散的质量守恒原理,并以L^∞（0,T;L^2(Ω））范收敛于原解。     

凝析气田自流注气的数值模拟研究

凝析气藏通常采用衰竭式开发,或采用循环注气以保持压力的方式进行开发。利用气藏气对凝析气藏进行自流注气无疑是一种经济有效的开采方法。将三维三相拟四组分模型改进为气藏井筒气藏模型,定量描述了自流注气过程中自流注气量及气藏压力的变化。研究表明,自流注气是两个气藏压力平衡过程。数值模拟的结果表明,压力平衡时间与气藏系统的产气速度、两个气藏的储量、地层压差,以及边底水能量大小有关。自流注气过程中的自流注气量是动态变化的,其变化符合产量的调和递减规律。该研究成果可为实际气田自流注气开采方案设计提供理论依据     

锥形旋转射流井底流场的数值模拟

针对旋转水射流代替机械钻头破岩钻径向水平井的实际要求,运用数值模拟方法研究了旋转射流在井底的流动规律。结果表明,旋转射流的流动不同于普通圆射流,它的每一个质点都具有三维速度,射流能量不是集中在射流中心部位,而是集中在一定半径的圆环上,同时具有明显的扩散性,并在中心形成回流区。随着旋流强度的增大,射流扩散能力增强。数值模拟得到的流动规律与实验结果完全一致。     

含薄互储集层半空间中点源激发的地震波场的理论求解

应用 Biot 双相介质的波动理论,严格推导了含有双相介质层(油储集层)时点源激发的表面位移的 Fouri-er-Bessel 积分解.自然地将储集层岩性参致(孔隙度、粘滞系数、渗透率等)纳入积分解中.     

三维波动方程的数值频散关系及其叠前和叠后数值模拟

采用适于三维各向同性介质正演的有限差分格式 ,对声波方程进行了差分离散 ,导出了差分格式的数值频散关系。研究表明 ,三维正演的数值频散主要受差分精度及网格空间步长的影响 ,尤以后者影响更甚。在相同的条件下 ,剖分步长越大及差分精度越低 ,数值频散现象越明显。研究还表明 ,剖分网格有低通滤波作用 ,因而在正演时应使激发源的主要频率远离网格的截止频率 ,以保证精度。在此基础上对某地质模型进行了叠前和叠后数值模拟。计算结果表明 ,叠前模拟能够比较真实地反映地震波在地下介质中的传播规律 ,但计算量过大 ;叠后模拟时 ,爆炸界面法能较好地反映地下反射界面形状 ,但不宜计算较大的模型 ;平面波照射法可用于计算较大的模型 ,且能较好地反映波的动力学特征 ,但不能精确地反映界面形态