流体饱和多孔介质黏弹性动力人工边界

基于Biot流体饱和多孔介质本构方程，采用平面波和远场散射波经验叠加来反映外行波传播，以经验参数反映人工边界外行波动的衰减和多角度透射特性。在人工边界处分别施加反映固相和液相介质传播效应的弹簧及阻尼来模拟人工边界以外的无限域介质对来自有限域的外行波的能量的吸收作用。从而形成一种流体饱和多孔介质的黏弹性动力人工边界。数值算例表明：边界的精度和稳定性高于现有的黏性边界、黏弹性人工边界及一阶透射边界。     

波在饱和多孔介质与弹性固体介质交界面上的界面效应

本利用Ｂｉｏｔ饱和多孔介质动力学理论及饱和多孔介质和固体介质交界面的连续条件,推导了平面ＰＩ波以任意角饱和多孔介质向弹性固体介质入射时,反射快Ｐ１波,慢ＰＨ波和ＳＶ波和透射Ｐ波、ＳＶ波系数的计算公式,并给出一算例,本不考虑饱和多孔介质中液体渗流的能量耗散。     

弹性多孔介质中应力—孔隙压耦联理论在新丰江水库地震研究中的应用

本文用二维有限元法,将弹性多孔介质中应力-孔隙压耦联理论应用于新丰江水库诱发地震机制的研究。计算中,不仅考虑了均匀、非均匀模型的耦联,还考虑了断层的性质、产状对耦联作用的影响及库水位随时间的变化。 结果表明:附加应力与孔隙压耦联作用是诱发水库地震的重要因素。新丰江水库地震不仅与库区存在的断裂、高密度基底有关,还与库水的加、卸载方式有关。     

油水两相渗流多孔介质弹性波传播动力学模型研究

弹性波在储层渗流场中的传播与衰减规律是研究波场强化采油动力学机理的重要基础.基于等效流体理论和饱和静态流体弹性波传播Biot理论,建立油水两非混相流体渗流条件下储层多孔介质中弹性波传播的动力学模型,通过算例求解与分析,发现含油水两相渗流储层多孔介质中同时存在着3种纵波P1、P2、P3和1种横波S;受频率和含水饱和度的影响,各波波速和品质因子呈现出不同变化规律,4种体波波速与频率、饱和度正相关,P1、P2波品质因子与饱和度正相关,P3和S波品质因子与饱和度负相关;最后,通过与传统静态弹性波模型结果对比,进一步分析了宏观渗流场对弹性波传播特征的影响规律,为揭示低频人工地震波辅助强化采油技术的动力学机理和工艺参数优化提供了重要理论依据.

     

流体饱和多孔介质波动方程的有限元解法

讨论了流体饱和多孔介质中波传播问题的有限元解法,首先在Ｂｉｏｔ理论的基础 上,概述了数学问题的提法,然后提出了一种新型简便的人工边界上的无反射边界条件,同时 给出了有人工边界时流体饱和多孔介质波动方程的有限元计算公式．数值试验的结果表明, 本文提出的无反射边界条件和数值方法均很有效．     

流体饱和多孔介质波动方程的有限元解法

讨论了流体饱和多孔介质中波传播问题的有限元解法,首先在Ｂｉｏｔ理论的基础 上,概述了数学问题的提法,然后提出了一种新型简便的人工边界上的无反射边界条件,同时 给出了有人工边界时流体饱和多孔介质波动方程的有限元计算公式．数值试验的结果表明, 本文提出的无反射边界条件和数值方法均很有效．     

饱和多孔介质近场波动分析的一种黏弹性人工边界

采用有限模型数值方法求解能量开放系统中的波动问题时,虚拟人工边界的处理方式对计算结果的准确性和精度具有重要的影响.本文针对无限域饱和多孔介质中波传播问题的人工边界处理方式进行了研究,提出了饱和多孔介质近场波动分析的一种黏弹性人工边界处理方法.在考虑多孔介质中固相和液相的相互作用的情况下,通过在人工边界处分别施加反映固相和液相介质波传播效应的弹簧及阻尼来模拟饱和多孔介质中波的能量辐射效应影响.算例表明,本文建议的黏弹性人工边界具有较好的模拟效果.     

含气介质临界点、流体和骨架弹性参数的数值计算方法

目前临界孔隙度及孔隙介质的研究主要为测试分析的实验方法,如何运用数值计算方法求取孔隙流体介质临界点和组分弹性参数一直是备受关注的课题.本文提出了求取孔隙介质临界点、孔隙内流体和骨架弹性参数的数值计算公式及方法,并结合含气样品测试数据实现了这种数值计算.文章首先给出了3个包含孔隙度的线性方程,通过有机组合每个线性方程中的两个系数得到求取有关弹性参数的数值计算公式.然后,详细阐述了计算的方法和步骤以及需要注意的问题,把室内含气砂岩样品实测的综合介质密度、纵波速度和横波速度作为数值计算的输入数据,求得了临界点和流体及骨架等相关弹性参数的具体数值.通过比较分析本文数学方法求得的计算数据与实验方法测得的实测数据,表明了该数值计算公式的正确性和数值实现方法的有效性.     

孔隙介质弹性波频散—衰减理论模型

储层地球物理学中,孔隙介质的各类弹性波模型常用于了解地层岩石物理性质.本文介绍了含油、气、水等物质的多相孔隙介质弹性波频散和衰减研究进展,给出了流体饱和与部分饱和孔隙介质中波传播的物理模型综述.根据孔隙介质中的固、流体分布情况,从相关基础理论和实验研究工作等方面出发,在宏观、微观和介观尺度上对流体替换、Biot孔隙力学、喷射流、Biot喷射流(BISQ)、等效球体癍块饱和、双重孔隙介质局部流动等现有主要合流体孔隙介质速度频散和衰减理论进行了回顾.研究表明,应力松弛过程是弹性波频散和衰减的基本机理,该过程由平衡特征时间刻画.该特征时间与孔隙介质的渗透率、流体粘性和体积模量紧密相关.当波频较低时,特征时间小于波周期,压力平衡得以发生,可以用等效流体模型描述波速;反之,当波频较高时,局部压差始终保持较高水平。整个骨架体积模量升高,等效模型面临困难,发展出斑块饱和模型.在分析了各类模型理论框架适用性以及所面临困难后,我们对未来研究方向给出了一些有意义的探讨.     

饱和多孔介质三维时域黏弹性人工边界与动力反应分析的显式有限元法

本文基于Biot的饱和多孔介质本构方程,考察具有辐射阻尼的外行球面波,推导了饱和多孔介质三维黏弹性人工边界的法向和切向边界方程;在已有的饱和多孔介质二维显式有限元数值计算方法基础上,提出该理论的三维方法,并开发了实现该三维方法的有限元程序.算例表明饱和多孔介质三维时域黏弹性人工边界与动力反应分析的显式有限元法具有较好的精度和稳定性.     

流体饱和多孔隙介质弹性波方程边界元解法研究

基于流体饱和多孔隙各向同性介质模型,本文首先推导了流体饱和多孔隙介质中弹性波传播的频率域系统动力方程及边界积分方程,然后给出了流体饱和多孔隙介质弹性波方程的基本解,最后,利用本文给出的边界元方法对流体饱和多孔隙各向同性介质中的弹性波传播进行了数值模拟.结果表明：不论是从固相位移,还是液相位移的地震合成记录都能看到明显的慢速P波,本文提出的流体饱和多孔隙介质弹性波边界元法是有效可行的.     

层状各向同性多孔介质中井间地震的数值计算

依据Ｂｉｏｔ理论，提出计算理论地震图的相应传播矩阵算法，并针对井间地震震源和检波器的不同配置进行数值计算，结果表明，本文算法高频可达１０００Ｈ２以上．文中算例说明利用震源和检波器相对于储油层的不同位置，可从反射震相记录波形上的差异，判定储油层存在与否．     

各向异性介质中的弹性阻抗及其反演

地震反演已成为油藏描述中的重要组成部分.绝大多数的常规地震反演是叠后地震数据体上进行,很少考虑各向异性存在的情况.随着勘探开发的发展,地震各向异性和叠前地震波阻抗反演引起了人们极大关注.本文在各向同性介质中弹性阻抗研究基础上,推导出了各向异性介质中的弹性阻抗方程,提出了地震各向异性介质中用弹性阻抗进行储层参数描述的技术路线和框架,并对反演过程中存在的问题进行了有益探讨.     

含气介质临界点、流体和骨架弹性参数的线性数值计算方法(英文)

临界孔隙度及其孔隙介质的研究均为测试分析的实验方法,如何运用数值计算方法求取孔隙流体介质临界点、流体和骨架的弹性参数是备受关注的课题。本文提出了求取临界点、流体和骨架弹性参数的数值计算公式及方法,并结合含气样品测试数据实现了这种数值计算。首先,基于孔隙度≯为白变量,而密度ρ,密度与横波速度平方的乘积ρVs^2和密度与纵波速度平方的乘积ρVp^2均为函数的三个线性方程,从每个线性方程中的两个系数的有机组合得到求取有关弹性参的线性数值计算公式。然后,详细阐述了数值计算的步骤和方法以及需要注意的问题,把室内含气砂岩数的线性数值计算公式。然后,详细阐述了数值计算的步骤和方法以及需要注意的问题,把室内含气砂岩样品测试的整体介质密度、纵波速度和横波速度作为数值计算的输入数据,求取了临界点和流体及骨架共计11个弹性参数的具体数值。通过数学方法的计算数据与实验方法的测试数据的比较分析,表明了本文数值计算公式的正确性和实现方法的有效性。本文提出的求取含气介质临界点、流体和骨架弹性参数的“数值计算公式及方法”含义清晰且形式简洁,为孔隙流体介质的数值计算分析和流体属性研究提供了可能的新方法。     

饱和多孔介质多相流体系中孔隙压传播数值分析

孔隙压波传播与衰减对解释深部流体在地震孕育、触发和断裂活动过程中的作用机制,以及石油与天然成藏机理等方面具有重要的学术和应用价值.目前孔隙压波的理论与计算以孔隙弹性介质理论、BISQ 理论为基础.Biot孔隙弹性理论仅考虑流体沿着波的传播方向发生流动,严重低估了Biot损失.Squirt机理考虑了微裂隙中粘性流体的局部流动,能够合理地解释超声波的衰减特征.     

地壳介质弹性、电性各向异性理论及对地震过程的联合解释

本文共分 1 0章 :第 1章为引论 ;第 2～ 6章为基础理论 ;第 7～ 1 0章为应用研究。主要研究内容可分为 5个方面。1 )第 1章引论部分综述了地壳介质各向异性现象和分类 ,从发展沿革、理论方法及成果、存在的问题等几个方面分别论述了地壳介质弹性各向异性和电性各向异性。在此基础上提出了本文的研究思路和主题 :将地震学中的S波分裂法与电法中的大地电磁测深法联合起来解释地震过程中的各向异性动态特征与发震应力场的关系。围绕这个主题设定了 2个核心 :一是建立地壳介质电性各向异性微观模型 ;二是基于 EDA(广泛扩容各向异性 )裂隙模型…     

声波在两种多孔介质界面上的反射和透射

本文导出了声波在两种多孔介质界面上反射、透射的一般计算公式.作为例子,数值计算了P₁波入射于界面时,P₁、P₂和S波的反射、透射系数与声波频率、入射角等量之间的关系.结果表明,各种模式波的反射、透射系数与入射角、多孔介质性质有关,在Biot特征频率附近与频率有关,并用界面两侧的法向能流相等验证了结果的正确性.若把多孔介质当作均匀固体处理,将会得到显著不同的结果.     

江苏地区介质非弹性衰减和场地响应研究

选取江苏数字地震台网在2001年1月-2010年7月记录到43次高质量数字地震波形,采用Atkinson(1992)和Moya(2000)的方法,反演得到江苏地区的介质非弹性衰减特性和场地响应.     

三维层状孔隙介质中弹性波的一种积分表达式II：正确性检验和数值模拟实验

基于前一篇文章中得到的关于三维层状孔隙介质中弹性波场的积分形式半解析解,本文通过离散波数法开展了数值模拟.将全空间均匀孔隙介质中单力点源和爆炸点源作用下弹性波场的解析解和我们的数值模拟结果进行对比,发现两者是完全一致的.而在一个两层半空间模型下的数值模拟,验证了固相位移Green函数的9组空间互易性情况.通过以上两种对比检验,验证了半解析解理论公式、数值模拟方法以及相应程序代码的正确性和可靠性.随后利用敏感度分析研究了不同的介质参数变化对爆炸点源在界面上会产生的反射波场的影响.通过垂直地震剖面模型的数值模拟,发现弹性波场能很好地反映孔隙介质物理性质的变化,同时也讨论了动力协调这一孔隙介质中的特殊现象.我们发展的基于半解析解的数值模拟方法可以为三维层状孔隙介质中弹性波传播特征的研究提供一种可供选择的有效工具和手段.