A 2D coupled hydro-thermal model for the combined finite-discrete element method

Based on the combined finite-discrete element method (FDEM), a two-dimensional coupled hydro-thermal model is proposed. This model can simulate fluid flow and heat transfer in rock masses with arbitrary complex fracture networks. The model consists of three parts: a heat conduction model of the rock matrix, a heat-transfer model of the fluid in the fracture (including the heat conduction and convection of fluid), and a heat exchange model between the fluid and rock at the fracture surface. Three examples with analytical solutions are given to verify the correctness of the coupled model. Finally, the coupled model is applied to hydro-thermal coupling simulations of a rock mass with a fracture network. The temperature field evolution, the effect of thermal conductivity of the rock matrix thermal conductivity and the fracture aperture on the outlet temperature are studied. The coupled model presented in this paper will enable the application of FDEM to study rock rupture driven by the effect of hydro-thermo-mechanical coupling in geomaterials such as in geothermal systems, petroleum engineering, environmental engineering and nuclear waste geological storage.

     

A 2D mixed fracture–pore seepage model and hydromechanical coupling for fractured porous media

A novel two-dimensional mixed fracture–pore seepage model for fluid flow in fractured porous media is presented based on the computational framework of finite-discrete element method (FDEM). The model consists of a porous seepage model in triangular elements bonded by unbroken joint elements, as well as a fracture seepage model in broken joint elements. The principle for determining the fluid exchange coefficient of the unbroken joint element is provided to ensure numerical accuracy and efficiency. The mixed fracture–pore seepage model provides a simple but effective tool for solving fluid flow in fractured porous media. In this paper, examples of 1D and 2D seepage flow in porous media and porous media with a single fracture or multiple fractures are studied. The simulation results of the model match well with theoretical solutions or results obtained by commercial software, which verifies the correctness of the mixed fracture–pore seepage model. Furthermore, combining FDEM mechanical calculation and the mixed fracture–pore seepage model, a coupled hydromechanical model is built to simulate fluid-driven dynamic propagation of cracks in the porous media, as well as its influence on pore seepage and fracture seepage.

     

Seepage-heat transfer coupling process of low temperature return water injected into geothermal reservoir in carbonate rocks in Xian County,China

Fracture seepage and heat transfer in the geothermal reservoir of carbonate rocks after the reinjection of low temperature geothermal return water is a complex coupling process,which is also the frontier of geothermal production and reinjection research. Based on the research of cascade comprehensive development of geothermal resources in Beijing-Tianjin-Hebei (Xian County),the carbonate geothermal reservoir of Wumishan formation in the geothermal field in Xian County is investigated. With the development of the discrete fracture network model and the coupling model of seepage and heat transfer,the numerical solution of seepage field and temperature field with known fracture network is reached using the finite element software COMSOL,and the coupling process of seepage flow and heat in carbonate rocks is revealed. The results show that the distribution of temperature field of fractured rocks in geothermal reservoir of carbonate rocks has strong non-uniformity and anisotropy. The fracture network is interpenetrated,which constitutes the dominant channel of water conduction,and along which the fissure water moves rapidly. Under the influence of convective heat transfer and conductive heat transfer,one of the main factors to be considered in the study of thermal breakthrough is to make the cold front move forward rapidly. When the reinjection and production process continues for a long time and the temperature of the geothermal reservoir on the pumping side drops to a low level,the temperature of bedrocks is still relatively high and continues to supply heat to the fissure water,so that the temperature of the thermal reservoir on the pumping side will not decrease rapidly to the water temperature at the inlet of reinjection,but will gradually decrease after a long period of time,showing an obvious long tail effect. The distribution of fractures will affect the process of seepage and heat transfer in carbonate reservoirs,which should be considered in the study of fluid thermal coupling in carbonate reservoirs.     

考虑土体非均质传热的河岸带水热耦合模型

为弥补当前河岸带水热交换模拟对土体非均质传热考虑不足的缺陷,在饱和-非饱和渗流及多孔介质传热理论基础上,引入土体有效导热系数模型,构建考虑土体非均质传热的河岸带水热耦合模型。结合COMSOL软件的特点,给出河岸带水热耦合模型的实现方法及求解流程,并通过河岸带温度和水位原型观测资料验证和对比分析不同有效导热系数模型下河岸带水热耦合模型的模拟效果。结果表明：与传统不考虑土体非均质传热方法相比,该模型能够较好地反映河岸带水热交换过程。此外,基于Johansen有效导热系数模型的河岸带水热耦合模型模拟效果表现最佳,模拟结果与前人试验结果一致。研究成果可为河岸带水热交换及污染物迁移过程的深入研究提供技术支撑。     

Thermo-mechanical performance of layered transversely isotropic media around a cylindrical/tubular heat source

We develop a new numerical model based on a precise integration method to investigate the coupled thermo-mechanical performance of layered transversely isotropic media around a cylindrical/tubular heat source. To obtain the relational matrices of the extended precise integration method, we first convert the governing equations of the problem into ordinary differential matrix equations through the Laplace–Hankel transform. Then, the cylindrical heat source is divided into a series of plane heat sources, and the plane temperature load term is added to the state vector between layer elements. By combining the layer elements, we build a layered transversely isotropic numerical model containing a cylindrical heat source in the transformed domain. Finally, we solve the model in the transformed domain and obtain the solution of the problem in the real domain through the Laplace–Hankel transform inversion. The accuracy of this method is verified by comparing the solutions with the results of the analytical method and the finite element method. Then, we study the influence of the anisotropy of thermal parameters, the embedded depth, the length/radius ratio, the type of heat source and the stratification of the medium on the thermo-mechanical coupled performance.

     

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

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

热-应力-损伤耦合作用下深埋隧洞围岩稳定性分析

以热力学和弹塑性力学理论为基础,分析岩石热-力完全耦合作用及其对力学参数和热特性参数的影响,建立了岩石热-力-损伤耦合模型及其参数演化方程,以ABAQUS软件为平台对其进行二次开发,并通过典型算例验证了岩石热-力完全耦合的重要性。然后以某深埋软岩隧洞为例,研究温度和开挖卸载共同作用下的隧洞围岩力学行为和损伤过程。计算结果表明：温度对岩石的力学性质和损伤演化过程影响显著,开挖损伤和热应力诱发的损伤对围岩热力学参数的影响不可忽略;所提出的力学模型可以有效反映围岩损伤演化、调热圈演化以及热力学参数演化,具有一定的借鉴作用。     

Numerical modelling of heat shock-assisted rock fracture

This paper presents a numerical investigation on the effects of thermal shock as a pretreatment of rock prior to comminution. More specifically, the effect of heat shock-induced cracks on the uniaxial compressive strength of rock is numerically studied. The chosen constitutive model of rock employs a (strong) embedded discontinuity finite element formulation to describe cracks. The thermomechanical problem that governs the heat shock pretreatment of rocks is considered as an uncoupled problem because of a highly dominating role of the external heat influx. Two solution methods of the global problem are presented: an explicit-explicit dynamic scheme and an implicit-implicit quasi-static scheme. The model performance is tested in simulations on heterogeneous numerical rock samples subjected first to a heat shock pretreatment and then to a mechanical compression test. According to the results, the compressive strength of intact granite rock having the axial splitting failure mode can be substantially reduced by heat shock pretreatment.     

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

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

Thermal stresses in borehole heat exchangers

The backfilling materials of borehole heat exchangers (BHE), particularly the grout material, must provide a suitable thermal contact and ensure durability to the induced thermal stresses because of the heat loading. In this paper, the thermal stresses that occurred in BHEs because of heat injection or extraction is investigated with an analytical solution of a hollow cylinder model that is adapted for time‐dependent heat loading, the geometry of a BHE, and the thermo‐mechanical properties of surrounding ground conditions. Firstly, the hollow cylinder model is solved with the considered boundary conditions in 2D plane stress. Secondly, the temperature differences at the inner and outer circles of the cylinder are evaluated with the heat line source models for continuous and discontinuous loading to observe the impact of the heat loading schedule. The developed analytical solution for thermal stress investigation is validated with numerical models. It is demonstrated that the analytical solutions agree well with numerical results for two types of BHE configurations (co‐axial and single U‐shaped pipes). Furthermore, the calculated maximum stresses are compared with the tensile strength of grout materials obtained from Brazilian tests. It is predicted that the thermal contraction of the grout, partially constrained by the surrounding rock, generates tensile stresses that may lead to cracking in the BHE. According to the results, the stiffness of rock has a primary role on the developed tensile stresses, and the relationship between the thermal conductivity of the ground and of the grout induces a proportional impact on the magnitude of thermal stresses. Copyright © 2015 John Wiley & Sons, Ltd.     

黏性土失水开裂多场耦合离散元数值模拟

黏性土在失水过程中逐渐变形开裂,裂隙相互交错形成网络,这一过程涉及水热力等多场耦合的作用。本文基于南京大学自主研发的三维离散元软件MatDEM,采用紧密堆积离散元模型对土体开裂进行模拟。在此模型中,每个离散元单元代表一定体积的土颗粒、孔隙和孔隙水的集合体。单元具有含水量属性,并采用有限差分思想计算水分运移量,实现了水分场模拟。同时,考虑水分场对土体抗拉强度等力学性质的影响,建立水分场和应力场的耦合。在数值模拟中,假定土体表面水分以一定速率蒸发,根据试验数据由含水量计算单元直径和力学参数,从而实现黏土蒸发失水、收缩和开裂变形过程模拟。数值模拟与前人室内试验结果基本一致,能较好地再现开裂过程中的各个阶段。本文为研究多场作用下土体变形破坏模拟提供了一个新的思路。     

跨季节蓄热型地源热泵热传递规律研究

基于有限元分析法,建立了准三维非稳态传热模型。在试验验证的基础上,分析了跨季节蓄热型地源热泵蓄热过程中土壤温度、单位井深换热量、热作用半径随热泵运行时间的变化规律,讨论了土壤结构、入口水温、入口质量流量、热泵运行模式等对土壤传热规律的影响,并研究了土壤热平衡问题。结果表明:同一半径不同深度处,土壤温度增长幅度随土壤热扩散率的增加而增大;土壤热作用半径随热泵运行时间及入口水温的升高而增加,并逐渐趋于稳定;间歇运行模式下,地埋管附近土壤温度及换热量均呈波动式变化,且当径向距离大于0.3 m时,与连续模式一致;在满足换热量的情况下,流体质量流量不宜过大;系统运行一个周期（360 d）后土壤温度基本可以恢复,且流体入口温度不宜低于40 ℃。     

饱和稀疏裂隙岩体三维水流-传热过程中位移和应力的一种半解析计算方法

针对高放射核废深地质处置库近场环境,建立分布热源作用下饱和裂隙岩体三维水流-传热过程中位移和应力的一种半解析计算方法：采用Goodier热弹性位移势和Laplace变换计算由温度梯度产生的温梯位移和应力;考虑单一裂隙的情况,利用经典弹性力学的Boussinesq解和Cerruti解计算为满足边界条件的约束位移和应力,与温梯位移和应力叠加,可得总体热位移和应力;把裂隙面离散为矩形单元集合,采用极坐标系下的解析法计算包含奇点的单元积分,采用数值法计算与分布热源有关和不含奇点的单元积分。与基于裂隙面法向一维热传导假设的一种解析解对比,结果表明,半解析法与解析法的计算结果基本一致,但由于半解析法考虑岩石的三维热传导,因温度时空分布和演变的不同而导致不同的温梯应力。针对一个假想单裂隙岩体三维水流-传热过程,计算温梯位移和应力、约束位移和应力、总体位移和应力;结果表明,裂隙水流-传热可能对位移和应力的分布和演变有显著影响,距离分布热源较近的岩石因升温膨胀受到约束而出现压应力,而距离分布热源较远的岩石则可能因协调收缩受到约束而出现拉应力。     

套管式地埋管换热器热短路及换热性能

套管式地埋管换热器是深层地源热泵系统常用的换热装置。基于流体流动换热方程,建立套管式地埋管换热器与周围岩体之间的传热模型。以第一个供暖季为例,分析内管导热系数和循环水流量对换热性能的影响,并引入换热器效能对热短路现象进行评估,研究结果显示:内管导热系数越大,热短路现象越显著;热短路使内外管中循环水温差降低,管内出现热堆积,导致换热器换热功率降低;套管式地埋管换热器的换热功率随循环水流量的增大逐渐增大;内外管之间存在热短路时,出口水温随循环水流量的增大先升高后降低,随着流量增大,换热器效能增大。研究成果可为深层地源热泵系统中地埋管换热器的设计提供借鉴。      

Numerical analysis of thermal fracturing in subsurface cold water injection by finite element methods

Integration of poromechanics and fracture mechanics plays an important role in understanding a series of thermal fracturing phenomena in subsurface porous media such as cold water flooding for enhanced oil recovery, produced‐water reinjection for waste disposal, cold water injection for geothermal energy extraction, and CO₂ injection for geosequestration. Thermal fracturing modeling is important to prevent the potential risks when fractures propagate into undesired zones, and it involves the coupling of heat transfer, mass transport, and stress change as well as the fracture propagation. Analytical method, finite element method, and finite difference method as well as boundary element method have been used to perform the thermal fracturing modeling considering different degrees and combinations of coupling. In this paper, extended finite element method is employed for the thermal fracturing modeling in a fully coupled fashion with remeshing avoided, and the stabilized finite element method is employed to account for the convection‐dominated heat transfer in the fracturing process with numerical oscillation circumvented. With the thermal fracturing model, a hypothetical numerical experiment on cold water injection into a deep warm aquifer is conducted. Results show that parameters such as injection rate, injection temperature, aquifer stiffness, and permeability can affect the fracture development in different ways and extended finite element method and stabilized finite element method provide effective tools for thermal fracturing simulation. Copyright © 2012 John Wiley & Sons, Ltd.     

小尺度含水层热量运移试验研究及热弥散效应评估

将小尺度含水层热机械弥散系数模型应用于对流弥散传热过程中,推导了该条件下对流弥散热量运移的解析解,结合试验进行验证,对含水层热弥散效应进行评估,结果表明,热机械弥散系数等于1×10-2 W/(m?°C)可以作为热弥散对温度场影响的临界点,从而将热机械弥散系数的分布划分为不可忽略的三角区域和可以忽略的多边形区域;明确了自然含水层结构条件下纵向热弥散度范围,从小尺度的热弥散研究结果来看,纵向热弥散度最大值为厘米数量级,它与野外大尺度条件下热弥散度的研究成果有着显著差异,表明热弥散尺度效应的存在,这将是进一步开展研究工作的方向。     

随机固-热耦合模型与岩石热破裂数值试验研究

1. 大连轻工业学院 信息科学与工程学院,大连 116034;2. 辽宁工程技术大学 力学与工程科学系,辽宁 阜新 123000     

裂隙岩体不稳定温度场的复合单元算法研究

基于复合单元法,结合三维热传导-对流方程和“充填模型”,提出了裂隙岩体不稳定温度场的复合单元模型。该模型前处理简便快捷,计算网格生成时无需考虑裂隙的存在,网格剖分不受限制,随后利用复合单元前处理程序,依据裂隙的位置和方位将其自动离散在单元内。对常规热传导-对流方程进行自伴随性调整,应用变分原理,推导出裂隙岩体不稳定温度场的复合单元算法,该算法可分别计算出岩块子单元和裂隙的温度值,且可真实反映裂隙中水流与相邻岩块间的热能量交换规律。将复合单元数值模型计算的不稳定温度场结果与相应的实测数据进行对比分析可知,数值计算结果与实测数据基本一致,验证了裂隙岩体不稳定温度场复合单元算法的可靠性与有效性。算例分析表明,裂隙中水流与相邻岩块间有明显的热传导和热对流作用。     

Theoretical and numerical analyses of pore-fluid flow focused heat transfer around geological faults and large cracks

In this paper, theoretical and numerical methods are used to investigate pore-fluid flow focused temperature distribution patterns around geological faults and cracks of any length-scales in hydrothermal systems. If the far field inflow is uniform and the long axis of an elliptical fault of any length-scale is parallel to the far field inflow direction, a complete set of analytical solutions has been presented for the pore-fluid velocity, stream function and excess pore-fluid pressure around the elliptical fault embedded in fluid-saturated porous media. Because the analytical solutions are explicitly expressed in the conventional Cartesian coordinate system, not only can they be used to gain a theoretical insight into the pore-fluid flow patterns around geological faults and large cracks, but also they can be used as valuable benchmark solutions for validating any numerical methods. After a finite element computational model is validated by comparing the numerical solutions with the present analytical solutions, it is used to investigate pore-fluid flow focused heat transfer around geological faults in hydrothermal systems. Some interesting conclusions in relation to the effects of geological faults on pore-fluid flow focused heat transfer have been made through both the theoretical and the numerical analyses.