裂隙岩体渗流-传热耦合的复合单元模型

基于复合单元法建立了裂隙岩体渗流-传热耦合的复合单元模型。该模型前处理简便快捷,网格剖分不受限制,可依据裂隙的真实信息自动将其离散在单元内。其次,采用交叉迭代算法,对裂隙岩体的渗流场和温度场进行耦合分析,耦合算法不仅考虑了温度对流体运动黏度的影响,而且可计算裂隙中流体与相邻岩块间渗流-传热过程以及两者间的渗流量和热量交换。通过与已有近似解析解相比较,验证了复合单元耦合算法的可靠性。算例分析表明,渗流-传热耦合作用对裂隙岩体的渗流场和温度场均有一定的影响。分析了不同岩块热传导系数和裂隙开度对热能提取效率的影响,结果显示,岩块热传导系数越大、裂隙开度越大,低温流体从高温岩块中吸取的热能会较多,出口处流体温度下降得较快。     

岩体水力学基础(五)——岩体渗流场与应力场耦合的裂隙网络模型

岩体是由裂隙和其间的岩块组成。当岩块处于近乎隔水状态时，岩体中的渗流可看作为裂隙网络流。本文通过分析岩体裂隙网络渗流场、应力场以及它们之间的相互力学关系，提出了岩体渗流场与应力场耦合的裂隙网络模型及数值计算方法。     

分布热源作用下单裂隙岩体三维水流-传热的半解析计算方法

针对高放射性核废物地下处置库近场饱和裂隙岩体环境,提出一种由分布热源、饱和单裂隙和两侧无限大岩石构成的三维水流-传热简化模型,建立了控制微分方程和基于拉氏变换域格林函数的积分方程;采用矩形单元把裂隙面域离散化,利用极坐标下的解析方法计算包含奇点的单元积分,利用数值方法计算分布热源和不包含奇点的单元积分,建立拉氏变换域的线性代数方程组,求解后,利用拉氏数值逆变换,计算任意时刻裂隙水和岩石的温度分布。对两个无内热源、流场确定的计算模型进行了计算,与仅考虑岩石沿裂隙面法向一维热传导的解析解进行了对比。计算分析了分布热源作用下饱和单裂隙岩体的三维水流-传热特征及其对裂隙水流速、岩石热传导系数和热源热流集度的敏感度。计算结果表明：与直接采用高斯数值积分相比,提出的解析法奇异积分精度较高;就裂隙水温度而言,单裂隙岩体三维水流-传热半解析计算方法与解析法得到的结果基本一致,但由于半解析计算方法考虑了岩石的三维热传导,使得裂隙水的上游温度较低,而下游温度较高;无分布热源作用时,岩石热传导系数越大,裂隙水温度越低;裂隙水流速越大,裂隙进水温度对裂隙水和岩石温度分布的影响越明显;由于受到裂隙水流动传热的作用,分布热源对裂隙水温度和岩石温度的影响在裂隙水流的下游区域比较显著。     

岩体水力学基础（五）

岩体是由裂隙和其间的岩块组成。当岩块御于近乎隔水状态时，岩体中的渗流可看作裂隙网络流。本文通过分析岩体裂隙网络渗流场、应力场以及它们之间的相互力学关系，提出了岩体渗流与应力场耦合的裂隙网络模型及数值计算方法。     

米尺度裂隙岩体模型水流-传热试验的数值模拟分析

为了研究高放射性核废物地下处置库近场的水流-传热耦合问题,采用国内高放废物地下处置库预选场址--甘肃北山地区的花岗岩石块体,加工组合成米尺度的规则裂隙岩体模型,设置边界热源和裂隙水流,试验模拟裂隙水水流与传热之间的相互作用。作为该室内模型试验的前期理论研究,采用等效孔隙介质数值模型,着重分析了裂隙开度、裂隙流量和热源功率对流场和温度场的影响。在设定条件下,计算分析表明：热传导和裂隙水水流由热源作用初期的不耦合很快转化为耦合;不流动的裂隙水主要表现为热存储和热传导,而流动的裂隙水还引起流动传热和水与岩石之间的对流换热,使岩体温度场明显不同于单纯热传导的情况;如果保持裂隙水流量不变,则裂隙开度的变化对水流-传热影响不大;如果保持裂隙水流速不变,则裂隙开度的变化对水流-传热影响显著;热源功率越大,通过裂隙水的热流量越大,裂隙水压强越大,而当温度超过100 ℃时,裂隙水会因汽化而压强显著增大;加热7 d时,热量的输入和输出几乎相等,裂隙水流带走的热量接近热源供给的热量,模型系统基本达到了热平衡。     

裂隙水流运动的入渗试验及数值模拟研究

裂隙对喀斯特水流起着重要控制作用。选取具有垂向与水平向交叉裂隙的岩石剖面,做单环注水入渗试验,基于Navier-Stokes(N-S)方程组建立入渗水流数值模型,运用试错法,以实测单环水位变化为目标推求裂隙水力隙宽,计算入渗水量在不同裂隙中的比例;利用立方定律,估计出裂隙饱和渗透系数,基于地下水流连续性方程构建裂隙水流数值模型;比较两种数值模型计算结果。模拟结果表明:N-S方法推求的裂隙等效水力宽度远小于实测裂隙宽度;裂隙水流大部分沿垂向裂隙下渗,但也有部分水流沿横向裂隙渗流;立方定律与N-S方法数值模拟结果差异不大,前者估计的裂隙下渗率稍小于后者,相对而言,N-S方法能更精确地描述裂隙水流运动过程。     

采用复合单元法模拟裂隙多孔介质变饱和流动

采用复合单元法建立了模拟裂隙多孔介质变饱和流动的数值模型。该模型具有以下特点：裂隙不需要离散成特定单元，而是根据几何位置插入到孔隙基质单元中形成复合单元；在复合单元中，分别建立裂隙流和孔隙基质流的计算方程，二者通过裂隙?基质界面产生联系并整合成复合单元方程；复合单元方程具有和常规有限单元方程相同的格式，因此，可以使用常规有限单元方程的求解技术。采用欠松弛迭代、集中质量矩阵以及自适应时步调节等技术，开发了裂隙多孔介质变饱和流动计算程序。通过模拟一维干土入渗和复杂裂隙含水层内的流动问题，验证了该模型的合理性和适用性。模拟结果为进一步认识非饱和裂隙含水层地下水流动特性提供了理论依据。     

裂隙岩体流-热耦合传热的三维数值模拟分析

通过对潘西煤矿水文地质条件的分析,基于裂隙岩体的流-热耦合数学模型,描述了裂隙岩体渗流场分布和水流及岩体的温度场分布,并结合边界条件及计算参数对裂隙岩体的流-热耦合传热进行了数值模拟和分析。数值模拟结果表明,岩体内裂隙水流所引发的热量迁移,对裂隙岩体的温度场分布有重要影响。断裂带及地下水流的存在改变了岩体的原有温度场分布。在渗流初期,温度梯度矢量沿渗流方向向两侧岩体方向流动,由于两侧岩体的渗透性系数低于断裂带处的渗透性系数,右侧等温线及温度梯度矢量方向逐渐向渗流方向移动,改变了两侧岩体的温度场分布。通过对断裂带内裂隙水流渗透性系数的折减,分析渗透性系数发生变化时对岩体温度场分布的影响,渗透性系数越大,伴随的热量迁移增大,对岩体的温度场分布的影响也越大。     

填砂裂隙岩体渗流传热模型试验与数值模拟

选取中国高放射核废物地下处置库重要预选场区--甘肃北山地区的花岗岩,加工组合成规则裂隙岩体,将垂直裂隙用粒径为0.5～0.63 mm的砂土填充,进行了裂隙水渗流传热试验;对模型试验进行了数值模拟,进而计算分析了热源温度、裂隙水流速和裂隙开度变化对裂隙岩体模型稳态温度场的影响。模型试验表明,当热源温度维持在120 ℃时,裂隙水仍无相变,裂隙岩体模型稳态温度场分布规律与热源温度为95 ℃时一致;热源温度越高,热源的水平影响距离越大,模型达到稳态需要的时间越长;裂隙填砂加强了裂隙两侧岩石之间的热传导,热源的水平影响距离和模型到达稳态需要的时间均明显大于无填充裂隙岩体模型的情况。模型试验得到的岩体模型温度场与数值计算得到的岩体模型温度场规律一致。试验过程中裂隙岩体模型在边界上存在一些热量散失,无法与数值计算中的绝热边界条件等同,致使试验数据低于数值计算值,并且热源温度越高,两者之间的差异越大。模型试验和数值计算均表明,邻近热源侧的裂隙水渗流对模型的温度场分布起控制作用,而远离热源侧的裂隙水渗流则主要影响该侧的边界温度和模型达到稳态所需要的时间。数值参数敏感性分析表明,裂隙水流速与裂隙开度越大,裂隙水对水平传热的阻滞作用越明显。     

北方岩溶含水层双重介质水流模型及其应用探讨

本文将北方岩溶介质划分为两种基本类型:溶隙网络—溶孔型和溶隙网络—微裂隙型。概化出裂隙岩溶介质的双重介质水流模型。解决了下面几个问题:（1）认识北方岩溶介质所具有的双重介质性;（2）分别对孔隙和微裂隙所构成的岩块域和溶隙网络域建立数学模型,并用岩块和裂隙的水量交换项将两者耦合起来。模型中岩块并未被概化成某一特殊几何形态如球体或平行板,而是按其实际产出形态处理;（3）在岩块中也应用了二维水流模型。以往作者多用一维流方程,目的在于求其解析解,这有许多不足之处;（4）在裂隙和岩块域均采用有限单元法,将两个域的模型联立求数值解,克服了解析法或解析—数值法联合的局限性。 在给定参数情况下,本文进行了典型水文地质条件下的计算,证明了模型是能够较好地刻画溶隙网络—溶孔或微裂隙型介质的特征和水流特点的。同时,本文成功地应用了济南西南郊岩溶含水层的实际资料进行了模拟计算,对双重介质方法的实用进行了听的尝试。     

Algorithm of Coupled Normal Stress and Fluid Flow in Fractured Rock Mass by the Composite Element Method

This paper presents a composite element algorithm of coupled normal stress and fluid flow process for fractured rock mass, developed from the composite element method (CEM). The coupled relation between the fracture flow and normal stress makes use of the “filled model”, which examines the asperities in the fracture as a layer of granular medium having high porosity and being clipped by the two parallel plates. The existence of fractures is not considered in the mesh generation, but it will be considered explicitly in the mapped composite element. The coupled normal stress and fluid flow process has been simulated by applying a cross iterative algorithm between the two fields. The proposed algorithm considers not only the flow through the fractures, but also the flow exchange between fractures and the surrounding rock blocks. In addition, it can be used for both the filled and non-filled fractures. The verification of the proposed algorithm has been conducted through the illustration of three examples by comparison with the conventional finite element method (FEM), from which the advantages and reliability of the proposed algorithm have been shown clearly.     

裂隙岩体变形模量尺寸效应研究Ⅰ:有限元法

通过对裂隙交叉和非贯通裂隙进行的网格剖分研究,建立了随机裂隙岩体网格生成方法,并编制了裂隙岩体有限元网格自动生成程序。根据锦屏电站辅助洞白山组大理岩裂隙分布调查结果,用蒙特卡洛法生成随机裂隙网络,进行了岩石-裂隙二元岩体网格自动剖分,将试验室获得的岩石和裂隙的力学参数分别赋予岩体模型中的岩石和裂隙,通过有限元方法研究了裂隙岩体等效变形模量的尺寸效应和各向异性。尺寸效应研究结果表明,所研究的裂隙岩体的等效变形模量表征单元体REV为8 m;各向异性研究结果表明,在REV尺寸下,所研究的裂隙岩体各个方向变形模量可以通过柔度张量进行拟合,且张量拟合误差不超过8%。     

Numerical estimation of REV and permeability tensor for fractured rock masses by composite element method

The Monte Carlo method is used to generate parent stochastic discrete fracture network, from which a series of fractured rock samples of different sizes and orientations are extracted. The fracture network combined with a regular grid forms composite element mesh of the fractured rock sample, in which each composite element is composed of sub‐elements incised by fracture segments. The composite element method (CEM) for the seepage is implemented to obtain the nodal hydraulic potential as well as the seepage flow rates through the fractured rock samples. The application of CEM enables a large quantity of stochastic tests for the fractured rock samples because the pre‐process is facilitated greatly. By changing the sizes and orientations of the samples, the analysis of the seepage characteristics is realized to evaluate the variation of the permeability components, the existence of the permeability tensor and the representative element volume. The feasibility and effectiveness are illustrated in a numerical example. Copyright © 2008 John Wiley & Sons, Ltd.     

A finite element method for modeling coupled flow and deformation in porous fractured media

Modeling the flow in highly fractured porous media by finite element method (FEM) has met two difficulties: mesh generation for fractured domains and a rigorous formulation of the flow problem accounting for fracture/matrix, fracture/fracture, and fracture/boundary fluid mass exchanges. Based on the recent theoretical progress for mass balance conditions in multifractured porous bodies, the governing equations for coupled flow and deformation in these bodies are first established in this paper. A weak formulation for this problem is then established allowing to build a FEM. Taking benefit from recent development of mesh‐generating tools for fractured media, this weak formulation has been implemented in a numerical code and applied to some typical problems of hydromechanical coupling in fractured porous media. It is shown that in this way, the FEM that has proved its efficiency to model hydromechanical phenomena in porous media is extended with all its performances (calculation time, couplings, and nonlinearities) to fractured porous media. Copyright © 2015 John Wiley & Sons, Ltd.     

FDEM-TH3D: A three-dimensional coupled hydrothermal model for fractured rock

This paper proposes a three-dimensional coupled hydrothermal model for fractured rock based on the finite-discrete element method to simulate fluid flow and heat transport. The 3D coupled hydrothermal model is composed of three main parts: a heat conduction model for the rock matrix, a heat transfer model for the fluid in the fractures (including heat conduction and heat convection), and a heat exchange model between the rock matrix and the fluid in the fractures. Four examples with analytical solutions are provided to verify the model. A heat exchange experiment of circulating water in a cylindrical granite sample with one fracture is simulated. The simulation results agree well with the experimental results. The effects of the fracture aperture, fluid viscosity, and pressure difference on the heat exchange between the fluid and rock are studied. Finally, an application concerned with heat transport and fluid flow in fractured rock is presented. The simulation results indicate that the 3D fully coupled hydrothermal model can capture the fluid flow and temperature evolution of rocks and fluids.     

分布热源作用下裂隙岩体渗流-传热的拉氏变换-格林函数半解析计算方法

以裂隙岩体高放射性核废物地下处置库性能评估为目标,提出了分布热源作用下单裂隙岩体渗流-传热的简化概念模型、控制微分方程和拉氏变换-格林函数半解析法,为进一步采用半解析法计算分布热源作用下多裂隙岩体的渗流-传热问题奠定了基础。针对单裂隙岩体的渗流-传热问题,建立考虑岩石内热源和二维热传导的控制微分方程,利用拉氏变换域微分方程的基本解建立格林函数积分方程,采用解析法处理其中的奇点,通过数值积分和拉氏数值逆变换求解,计算任意时刻裂隙水和岩石的温度分布。通过算例,与基于岩石一维热传导假定的解析解进行了对比,并计算分析了分布热源作用下单裂隙岩体的渗流-传热特征及其对裂隙开度、岩石热传导系数和热流集度的敏感度。算例表明,（1）就裂隙水温度而言,由于考虑了岩石的二维热传导,拉氏变换-格林函数半解析解小于基于岩石一维热传导假定的解析解;（2）裂隙水温度和岩石温度对裂隙开度和热流集度的敏感度较大,对岩石热传导系数的敏感度较小。     

Thermo-hydro-mechanical modeling of fracturing porous media with two-phase fluid flow using X-FEM technique

In this paper, a fully coupled thermo-hydro-mechanical model is presented for two-phase fluid flow and heat transfer in fractured/fracturing porous media using the extended finite element method. In the fractured porous medium, the traction, heat, and mass transfer between the fracture space and the surrounding media are coupled. The wetting and nonwetting fluid phases are water and gas, which are assumed to be immiscible, and no phase-change is considered. The system of coupled equations consists of the linear momentum balance of solid phase, wetting and nonwetting fluid continuities, and thermal energy conservation. The main variables used to solve the system of equations are solid phase displacement, wetting fluid pressure, capillary pressure, and temperature. The fracture is assumed to impose the strong discontinuity in the displacement field and weak discontinuities in the fluid pressure, capillary pressure, and temperature fields. The mode I fracture propagation is employed using a cohesive fracture model. Finally, several numerical examples are solved to illustrate the capability of the proposed computational algorithm. It is shown that the effect of thermal expansion on the effective stress can influence the rate of fracture propagation and the injection pressure in hydraulic fracturing process. Moreover, the effect of thermal loading is investigated properly on fracture opening and fluids flow in unsaturated porous media, and the convective heat transfer within the fracture is captured successfully. It is shown how the proposed computational model is capable of modeling the fully coupled thermal fracture propagation in unsaturated porous media.     

裂隙岩体渗流耦合传热分析

以地下裂隙岩体在裂隙水—孔隙水和温度场之间耦合作用为研究对象,对热和流体流动控制方程采用有限容积数值方法进行离散求解,设置了六种裂隙水—孔隙水流速方案,给出了部分无量纲温度场,并分析了传热与流动原因。分析结果表明：岩体内裂隙水—孔隙水引发的热质迁移对裂隙岩体的温度场分布有重要影响;当裂隙岩体内发生地下裂隙水—孔隙水渗流、及热量的转移时,会产生渗流场、温度场之间的耦合作用;裂隙内水流渗透速度是影响岩体温度的主要因素,孔隙内水流渗透速度是影响岩体温度的次要因素,温差主要发生在裂隙水边界层处。     

单轴压缩条件下岩溶化裂隙岩体损伤破坏特征研究

岩溶化裂隙岩体是普遍发育于自然界中具有初始损伤的岩体。为了研究岩溶化裂隙岩体损伤破坏特征,本文以贵州某地赋存的溶蚀岩体为研究对象,运用损伤力学理论构建岩溶化裂隙岩体在单轴压缩条件下的损伤演化模型,并建立岩溶化裂隙岩体损伤演化方程。采用颗粒流数值软件进行单轴压缩数值试验,进一步研究岩溶化裂隙岩体试件在单轴压缩条件下的损伤演化特征,分析岩溶化裂隙岩体的微观损伤特征。结果表明：岩溶化裂隙岩体的初始损伤主要包括溶蚀损伤和裂隙损伤。岩溶化裂隙岩体的初始损伤随着溶蚀率的增加而增加,最终增加速率趋于平缓;岩溶化裂隙岩体的损伤演化曲线均呈“S”型变化,先缓慢增加,再迅速增加,最后缓慢增加至损伤值1;岩溶化裂隙岩体存在异构特征,导致破坏裂隙起源于具有初始损伤的溶蚀孔洞和裂隙处,随后裂隙经历萌发、扩张和剪切作用、数量和长度增加以及裂隙贯通4个阶段后发生宏观破坏。