岩体结构与岩体水力耦合计算模型

水力耦合作用是岩体中力学过程与渗流过程相互作用的物理过程。水力耦合机理的理解是水力耦合分析的关键问题,其耦合机理是由岩体结构特性决定的。在分析岩体水力耦合过程基础上,根据岩体的基本结构及代表性单元体（REV）是否存在提出了建立水力耦合模型的方法。当裂隙岩体中不存在代表性单元体（REV）时,提出了裂隙岩体多重介质流固耦合分析的全耦合数学模型,给出流固耦合模型数值方法求解的数学模型及有限元计算表达式。程序编制和验证工作正在进行中。     

基于粒径分维值检验的松散介质渗透系数空间估值——以河北平原为例

为了克服松散含水介质非均质性各向异性对介质渗透系数的求取所带来的困难,以水动力条件为线索,利用互相关法寻找同一时代地层及同一时代地表水的主、次流方向.在主流方向上,建立了石家庄-束鹿-衡水段介质颗粒粒径随距离递变的数学模型,并以此为基础推求无钻孔地段的粒径.同时,运用分形理论,计算了中、上更新统地层主流方向平均粒径空间分布的分维值,并将其作为检验粒径估值正确与否的判据.通过不同岩性的渗透系数-深度模型,实现了用平均粒径预测渗透系数的全过程.     

Isogeometric analysis of THM coupled processes in ground freezing

An isogeometric analysis (IGA) based numerical model is presented for simulation of thermo-hydro-mechanically (THM) coupled processes in ground freezing. The momentum, mass and energy conservation equations are derived based on porous media theory. The governing equations are supplemented by a saturation curve, a hydraulic conductivity model and constitutive equations. Variational and Galerkin formulation results in a highly nonlinear system of equations, which are solved using Newton-Raphson iteration. Numerical examples on isothermal consolidation in plane strain, one-dimensional freezing and heave due to a chilled pipeline are presented. Reasonably good agreements were observed between the IGA based heave simulations and experimental results.     

Application of stress-pore pressure coupling theory for porous media to the Xinfengjiang reservoir earthquakes

Theory of the coupling of stress-pore pressure in the saturated, elastic porous media is used in the study of the formation mechanism of the Xinfengjiang reservoir-induced earthquakes. Based on the results, it is believed that compared with the mechanism of additional stress in the vicinity of the reservoir, the mechanism of the coupling of additional stress and pore pressure may be more well-founded for the occurrence of reservoir-induced earthquakes.     

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

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

Velocities of elastic waves andv P/v S ratios in dry and water-saturated rock samples at high pressure

Elastic wave velocities and the ratiov _P/v _S were studied for dry and initially saturated samples of carbonate and crystalline rocks at pressures to 2 kbar. In initially saturated samplesv _P increases in crystalline rock, whereas in sedimentary rock it may either increase or decrease with increasing pressure. Under the same conditionsv _S remains approximately constant in crystalline rocks and decreases in sedimentary samples. The ratiov _P/v _S as a function of pressure increases in dry rocks and decreases in saturated ones. Saturated samples always have higherv _P/v _S ratios than dry ones. It appears that the geometry of the pore space influences the acoustic properties of saturated crystalline rocks.     

辽宁地区波速比分布特征及构造意义

利用辽宁数字台网2001年以来的观测报告,采用单台和达法计算辽宁测震台网34个台站的平均波速比,重点分析辽宁地区波速比的空间分布特征.结果显示,辽宁地区的波速比空间分布呈横向不均匀性,下辽河盆地的波速比平均值低于整个辽宁地区波速比的均值,而辽西隆起和辽东隆起区的波速比却明显高于下辽河盆地的波速比和辽宁地区的均值.研究也...     

周期性层状含孔隙、裂隙介质模型纵波衰减特征

地震波在含孔隙、裂隙斑块饱和介质传播过程中会诱发多个尺度孔隙流体流动而产生衰减和速度频散.在含有宏观尺度“Biot流”和介观尺度“局域流”衰减诱导机制的周期性层状孔隙介质模型基础上,引入了微观尺度硬币型和尖灭型裂隙“喷射流”的影响,构建了周期性层状含孔隙、裂隙介质模型.利用双解耦弹性波动方程的方法数值计算了该模型地震频带的纵波衰减和速度频散并与周期性层状孔隙介质模型做了对比研究.分析了该模型在不同裂隙参数（裂隙密度、裂隙纵横比）及裂隙体积含量下的纵波衰减和频散特征,裂隙密度越高对于纵波衰减和频散的影响越大,裂隙纵横比越小,由裂隙引起的纵波衰减部分向高频段移动,裂隙体积含量越少,纵波衰减先降低后小幅增加再降低,频散速度增加,并逐渐接近于周期性层状孔隙介质模型的纵波衰减和频散速度曲线.最后研究了周期性层状含孔隙、裂隙介质模型有效平面波模量的高低频极限以及流固相对位移在该模型中的分布特征.     

含泥质低孔渗各向异性黏弹性介质中的波频散和衰减研究

本文定义了各向异性黏弹性参数修正因子,并将其引入到黏弹性模型中以体现泥质含量对黏弹性机制的影响,同时将波传播过程中孔隙介质骨架黏弹性力学机制与两种孔隙流体流动力学机制（Biot流动和喷射流动机制）有机地统一起来处理,从而给出了描述含泥质低孔渗孔隙各向异性介质中波传播规律的黏弹性Biot/squirt (BISQ)模型.数值计算结果表明,入射波的方位角、各向异性渗透率以及泥质含量等对含流体复杂孔隙介质中波频散和衰减的影响具有显著的方位各向异性特征,在低频范围内（地震波勘探频率）黏弹性力学机制对波传播能量的衰减起主导作用.     

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

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