温度梯度水分扩散对双重孔隙岩体中THM耦合影响的有限元分析

建立了一种饱和-非饱和遍有节理岩体的双重孔隙-裂隙介质热-水-应力耦合模型,并研制出相应的二维有限元程序.通过一个假定的位于非饱和双重孔隙-裂隙岩体中的高放废物地质处置库算例,就温度梯度水分扩散系数不同的三种工况,考察了岩体中的温度、孔隙水压力、饱和度、地下水流速和主应力的变化、分布情况.结果显示:各工况计算域中温度场及应力场基本相同,当岩体温度梯度水分扩散系数较大时,近场的负孔隙水压力上升到很高的数值,负裂隙水压力有所下降,饱和度亦有相应的变化,当温度梯度水分扩散系数小到一定程度后,其影响也将逐渐消失.     

温度梯度水分扩散系数对非饱和双重孔隙岩体中THM耦合影响的有限元分析

建立了一种饱和–非饱和遍有节理岩体的双重孔隙—裂隙介质热—水—应力耦合模型,并研制出相应的二维有限元程序。通过一个假定的位于非饱和双重孔隙—裂隙岩体中的高放废物地质处置库算例,就温度梯度水分扩散系数不同的三种工况,考察了岩体中的温度、孔隙水压力、饱和度、地下水流速和主应力的变化、分布情况。结果显示：各工况计算域中温度场及应力场基本相同, 当岩体温度梯度水分扩散系数较大时,近场的负孔隙水压力上升到很高的数值,负裂隙水压力有所下降,饱和度亦有相应的变化,当温度梯度水分扩散系数小到一定程度后,其影响也将逐渐消失。     

露天矿边坡饱和渗透系数随机场数值模型研究

为分析强降雨入渗及渗透系数空间变异性对闭坑露天矿边坡渗流场的影响程度,基于非饱和渗流理论和随机场理论,采用非侵入式随机方法,通过FISH语言编写非饱和区单元饱和度、渗透系数与基质吸力的修正函数,建立饱和-非饱和渗流随机场模型,开展强降雨作用的高大陡深岩质边坡渗流特征研究。研究结果表明,修正饱和-非饱和渗流随机场模型能够准确地描述露天矿边坡的降雨入渗过程。且降雨入渗主要影响到露天矿边坡浅层渗流场。随降雨持时变化,在坡面与地下水位线之间形成包围的且逐渐缩小的非饱和区。坡面最早出现暂态饱和区且厚度逐渐增加,但增幅逐渐放缓。同时坡面点孔隙水压力最早达到稳定值零,离坡面越远的点孔隙水压力达到稳定值零的用时会越长。该结论可为闭坑露天矿边坡的地质灾害风险预测提供参考意义。     

非饱和渗流下土石坝坝坡可靠性分析

将饱和-非饱和渗流场统一起来,建立了饱和-非饱和渗流的微分方程.根据极限平衡原理,考虑基质吸力对抗剪强度的影响,采用延伸的摩尔-库仑破坏准则,建立边坡在非饱和土条件下的安全系数公式.运用饱和-非饱和渗流有限元程序,对土石坝工程进行了计算和分析,探讨了稳态和瞬态下库水位变化时,坝体内浸润面的变化、瞬态孔隙水压力分布情况.进一步采用蒙特卡罗法把渗流计算所得的瞬态孔隙水压力分布用于坝坡可靠性分析,计算坝坡的瞬态可靠指标,对土石坝在非饱和、非稳定渗流作用下,不同水位下降速度对上游坝坡可靠性的影响进行了研究.结果表明:库水位下降对坝坡产生饱和-非饱和渗流作用,渗流场受降速的影响大,库水位下降速度增大,坝坡稳定性降低;坝坡可靠指标与库水位的关系呈抛物线变化趋势,渗流浸润线开始比较陡,最后趋于平缓,而坝坡可靠指标表现为开始减小,经过一段时间达到最小值,然后逐渐增大.     

岩体系统渗流场与应力场耦合的广义双重介质模型的应用研究

本文通过实验资料分析, 提出了裂隙岩体应力与渗流之间的分形几何关系式。以研究区岩体系统结构为基础, 提出了岩体系统应力场与渗流场耦合的广义双重介质模型, 并将该模型应用于某水利水电工程的渗流与控制中, 对坝区岩体进行了应力场和渗流场的耦合分析, 取得了比较满意的结果。     

水力压裂中裂纹扩展的数值模拟

水力压裂是在高压粘滞流体或清水作用下地层内裂缝起裂与扩展的过程。由于包含岩石断裂和流-固耦合等复杂问题,对该过程的数值模拟具有相当大的挑战性。本文建立基于有限元与离散元混合方法的裂纹模型,模拟岩石裂纹扩展,实现了连续向非连续的转化;建立双重介质流动模型,裂隙流作为孔隙渗流的压力边界,孔隙渗流反作用裂隙的压力求解,处理了流体在基岩与人工裂缝中的协调流动;将裂纹模型与流体流动模式进行结合,建立断裂-应力-渗流耦合形式的力学模型,进一步分析了水力压裂的基本过程,综合多种数值计算方法,编写程序,在验证岩体裂纹模型与双重介质流动模型有效性的基础上,对压裂过程进行复现,将模拟结果与文献结果进行了对比,并讨论了所构建模型的优缺点。     

非饱和土半空间Lamb问题及能量传输特性

地球表面绝大多数土层处于非饱和状态, 故采用传统饱和两相介质理论进行动力学分析时, 结果往往与实际情况不符. 针对这一问题, 本文以非饱和半空间作为研究对象, 基于连续介质力学和多孔介质理论, 考虑非饱和多孔介质中各相的质量守恒方程、动量守恒方程、本构方程以及有效应力原理等基本方程, 建立了以骨架位移、孔隙水压力和孔隙气压力为基本未知量的动力学控制方程. 针对非饱和半空间表面在竖向集中简谐荷载作用下的动力学响应及能量传输问题, 建立了频域内经典Lamb问题的轴对称计算模型, 采用Helmholtz分解法, 通过引入势函数Φ和Ψ表示骨架的位移分量, 结合本构方程获得了不同边界条件下半空间表面位移场和能量场等物理量的解析解答, 并通过数值算例对荷载参数(激振频率)、材料参数(饱和度、渗透系数)等影响因素进行了分析与讨论. 结果表明: (1)饱和度的升高或者激振频率下降, 都会提高非饱和半空间的表面位移幅值; (2)当渗透系数下降至一临界值时, 地表位移幅值会趋于一极限值, 并且透水(气)边界与不透水(气)边界条件下渗透系数的影响表现出明显的差异性.      

基于孔与裂隙网络模型的平行微裂隙对驱油的影响规律研究

认识双重多孔介质中油水两相微观渗流机制是回答形成什么类型的裂隙网络可提高油藏采收率的关键. 微裂隙的分布可以提高多孔介质的绝对渗透率,但对于基质孔隙中的流体介质,微裂隙的存在会引起多孔介质中局部流体压力和流场的变化,导致局部流动以微裂隙流动为主,甚至出现窜流现象,降低驱油效率. 本文基于孔与裂隙双重网络模型,在网络进口设定两条平行等长且具有一定间隔的微裂隙,分析微裂隙的相对间隔(微裂隙之间距离/喉道长度)和微裂隙相对长度(微裂隙长度/喉道长度)对于微观渗流特征的影响. 结果表明:随微裂隙相对长度的增加,出现驱油效率逐渐降低,相对渗透率曲线中的油水共渗区水饱和度和等渗点增加,油水两相的共渗范围减小等现象;随着微裂隙之间相对间隔增大,周围越来越多的基质孔穴间的压力差减小,在毛管压力的限制下,驱替相绕过这些区域,而导致水窜现象.      

裂隙岩体渗透系数确定方法研究

裂隙岩体渗透系数以及渗透主方向的确定对研究岩体渗透性大小及各向异性具有重要意义。高放废物地质处置库介质岩体的渗透性能将直接影响其使用安全性。本文运用离散裂隙网络模拟的方法对我国高放废物处置库甘肃北山预选区3#钻孔附近裂隙岩体进行了渗透性质分析。通过对3#钻孔1715～1780m段压水试验数据的反演,标定了离散裂隙网络渗流模型中的裂隙渗透参数(导水系数T)。利用标定的离散裂隙网络模型对场区裂隙岩体进行了渗流模拟,确定了该区域裂隙岩体的渗流表征单元体(REV)的尺寸大小以及渗透主值和主渗透方向。运用离散裂隙网络模型计算得出的渗透主值的几何均值与现场压水试验计算结果较接近,证明了计算结果的有效性。     

基于孔与裂隙网络模型的平行微裂隙对驱油的影响规律研究

认识双重多孔介质中油水两相微观渗流机制是回答形成什么类型的裂隙网络可提高油藏采收率的关键.微裂隙的分布可以提高多孔介质的绝对渗透率,但对于基质孔隙中的流体介质,微裂隙的存在会引起多孔介质中局部流体压力和流场的变化,导致局部流动以微裂隙流动为主,甚至出现窜流现象,降低驱油效率.本文基于孔与裂隙双重网络模型,在网络进口设定两条平行等长且具有一定间隔的微裂隙,分析微裂隙的相对间隔(微裂隙之间距离/喉道长度)和微裂隙相对长度(微裂隙长度/喉道长度)对于微观渗流特征的影响.结果表明:随微裂隙相对长度的增加,出现驱油效率逐渐降低,相对渗透率曲线中的油水共渗区水饱和度和等渗点增加,油水两相的共渗范围减小等现象;随着微裂隙之间相对间隔增大,周围越来越多的基质孔穴间的压力差减小,在毛管压力的限制下,驱替相绕过这些区域,而导致水窜现象.     

FEM ANALYSIS FOR INFLUENCES OF A FAULT ON COUPLED T-H-M-M PROCESS IN DUAL-POROSITY ROCK MASS

For the case in which a large geological structure like fault existing within the surrounding rock mass in the near field of a repository for high-level radioactive nuclear waste, one kind of coupled thermo-hydro-mechanical-migratory model of dual-porosity medium for saturatedunsaturated ubiquitous-joint rock mass was established. In the present model, the seepage field and the concentration field are double, but the stress field and the temperature field are single, and the influences of sets, spaces, angles, continuity ratios, stiffness of fractures on the constitutive relationship of the medium can be considered. At the same time, a two-dimensional program of finite element method was developed. Taking a hypothetical nuclear waste repository located at a rock mass being unsaturated dual-porosity medium as a calculation example, the FEM analysis for thermo-hydro-mechanical-migratory coupling were carried out under the condition of radioactive nuclide leaking for the cases with and without a fault, and the temperatures, pore pressures, flow velocities, nuclide concentrations and principal stresses in the rock mass were investigated. The results show that the fracture water in the fault flows is basically along the fault direction, and its flow velocity is almost three orders of magnitude higher than that of fracture water in rock mass; the nuclide concentration in the fault is also much higher than that without fault, and the nuclides move along the fault faster; moreover, the fault has obvious influences on the pore pressures and the principal stresses in the rock mass.     

Dual-Porosity Coarse-Scale Modeling and Simulation of Highly Heterogeneous Geomodels

Accurate upscaling of highly heterogeneous subsurface reservoirs remains a challenge in the context of modeling of flow and transport. In this work, we address this challenge with emphasis on the representation of the displacement efficiency in coarse-scale modeling. We propose a dual-porosity upscaling approach to handle displacement calculations in high resolution and highly heterogeneous formations. In this approach, the pore space is arranged into two levels of porosity based on flow contribution, and a dual-porosity dual-permeability flow model is adapted for coarse-scale flow simulation. The approach uses fine-scale streamline information to transform a heterogeneous geomodel into a coarse dual-continuum model that preserves the global flow pathways adequately. The performance of the proposed technique is demonstrated for two heterogeneous reservoirs using both black oil (waterflooding) and compositional (gas injection) modeling approaches. We demonstrate that the coarse dual-porosity models predict the breakthrough times accurately and reproduce the post-breakthrough responses adequately. This is in contrast to conventional single-porosity upscaling techniques that overestimate breakthrough times and displacement efficiencies (sweep). By preserving large-scale heterogeneities, coarse dual-porosity models are demonstrated to be significantly less sensitive to the level of upscaling, when compared to conventional single-porosity upscaling. Accordingly, the proposed upscaling approach is a relevant and suitable technique for upscaling of highly heterogeneous geomodels.     

Large-scale averaging analysis of multiphase flow in fractured reservoirs

The method of large-scale averaging is applied to derive and analyze a dual-porosity model of multiphase flow in naturally fractured reservoirs. The dual-porosity model contains the usual equations based on Darcy's law, and the coupling terms representing the fluid transfer between the matrix and the fractures. Both quasisteady and transient closure schemes are considered to obtain and analyze the fracture and matrix permeability tensors and the fluid transfer terms. The techniques developed here are not restricted to regular geometric fractures. Computational work aimed at showing the implications of the theory behind the derivation of the present dual-porosity model is also described. In particular, comparisons among the dual-porosity model, the single porosity model, and other dual-porosity models are presented through numerical experiments.     

Transient Pressure Behavior of Reservoirs with Discrete Conductive Faults and Fractures

Fractures and faults are common features of many well-known reservoirs. They create traps, serve as conduits to oil and gas migration, and can behave as barriers or baffles to fluid flow. Naturally fractured reservoirs consist of fractures in igneous, metamorphic, sedimentary rocks (matrix), and formations. In most sedimentary formations both fractures and matrix contribute to flow and storage, but in igneous and metamorphic rocks only fractures contribute to flow and storage, and the matrix has almost zero permeability and porosity. In this study, we present a mesh-free semianalytical solution for pressure transient behavior in a 2D infinite reservoir containing a network of discrete and/or connected finite- and infinite-conductivity fractures. The proposed solution methodology is based on an analytical-element method and thus can be easily extended to incorporate other reservoir features such as sealing or leaky faults, domains with altered petrophysical properties (for example, fluid permeability or reservoir porosity), and complicated reservoir boundaries. It is shown that the pressure behavior of discretely fractured reservoirs is considerably different from the well-known Warren and Root dual-porosity reservoir model behavior. The pressure behavior of discretely fractured reservoirs shows many different flow regimes depending on fracture distribution, its intensity and conductivity. In some cases, they also exhibit a dual-porosity reservoir model behavior.     

Numerical investigation of dual-porosity model with transient transfer function based on discrete-fracture model

Based on the characteristics of fractures in naturally fractured reservoir and a discrete-fracture model, a fracture network numerical well test model is developed.Bottom hole pressure response curves and the pressure field are obtained by solving the model equations with the finite-element method. By analyzing bottom hole pressure curves and the fluid flow in the pressure field, seven flow stages can be recognized on the curves. An upscaling method is developed to compare with the dual-porosity model(DPM). The comparisons results show that the DPM overestimates the inter-porosity coefficient λ and the storage factor ω. The analysis results show that fracture conductivity plays a leading role in the fluid flow. Matrix permeability influences the beginning time of flow from the matrix to fractures. Fractures density is another important parameter controlling the flow. The fracture linear flow is hidden under the large fracture density.The pressure propagation is slower in the direction of larger fracture density.     

页岩气藏的双重介质-离散裂缝模型

考虑页岩气藏开发中渗流的多尺度效应,提出了一个基于裂缝-孔隙双重介质的离散裂缝模型．在该模型中,基质、天然裂缝和人工压裂裂缝采用各自控制方程独立计算,不同介质之间通过流量交换相互关联．为分析模型可靠性,分别和基于渗透率粗化及压裂裂缝导流能力无穷大的模型对比．数值算例显示,伴随着网格细分,该模型与精确渗透率粗化模型具有相同计算精度,两者收敛速度均较快,但该模型易推广到多相流动问题,而等压模型对产量将有所高估．研究了地质参数和工艺参数对气井产量的影响规律．计算结果表明天然裂缝渗透率及基质孔隙扩散系数对产气速率有着重要影响,产气速率伴随着人工压裂裂缝导流能力、长度以及数目的增加而增加,但是增加幅度会逐步趋缓．     

论“岩体结构控制论”

经过长期实践和研究,作者于1984年提出岩体结构控制论是岩体力学的基础理论,并全面、系统地以岩体结构控制论为指导研究了岩体变形、岩体破坏及岩体力学性质的基本规律;提出岩体变形系山岩体材料变彤和岩体结构变形共同贡献的,岩体破坏系受岩体材料破坏和岩体结构破坏控制的;岩体力学性质不仅决定于岩体材料力学的性质,而且受控于岩体结构力学效应及环境因素力学效应。在此基础上,作者提出了岩体可以划分为连续介质、碎裂介质,块裂介质及板裂介质四种岩体力学介质,从而建立了完整的岩体结构力学理论体系。     

EXACT SOLUTION AND ITS BEHAVIOR CHARACTERISTIC OF NONLINEAR DUAL-POROSITY MODEL

IntroductionIn the computation of petroleum reservoir engineering design,the nonlinear quadraticgradient term is neglected by assuming small pressure gradient or small compressibility.Theassumption of small pressure gradient may cause significant errors i…