共查询到16条相似文献，搜索用时 500 毫秒
1.
针对高放射性核废物地下处置库近场饱和裂隙岩体环境,提出一种由分布热源、饱和单裂隙和两侧无限大岩石构成的三维水流-传热简化模型,建立了控制微分方程和基于拉氏变换域格林函数的积分方程;采用矩形单元把裂隙面域离散化,利用极坐标下的解析方法计算包含奇点的单元积分,利用数值方法计算分布热源和不包含奇点的单元积分,建立拉氏变换域的线性代数方程组,求解后,利用拉氏数值逆变换,计算任意时刻裂隙水和岩石的温度分布。对两个无内热源、流场确定的计算模型进行了计算,与仅考虑岩石沿裂隙面法向一维热传导的解析解进行了对比。计算分析了分布热源作用下饱和单裂隙岩体的三维水流-传热特征及其对裂隙水流速、岩石热传导系数和热源热流集度的敏感度。计算结果表明：与直接采用高斯数值积分相比,提出的解析法奇异积分精度较高;就裂隙水温度而言,单裂隙岩体三维水流-传热半解析计算方法与解析法得到的结果基本一致,但由于半解析计算方法考虑了岩石的三维热传导,使得裂隙水的上游温度较低,而下游温度较高;无分布热源作用时,岩石热传导系数越大,裂隙水温度越低;裂隙水流速越大,裂隙进水温度对裂隙水和岩石温度分布的影响越明显;由于受到裂隙水流动传热的作用,分布热源对裂隙水温度和岩石温度的影响在裂隙水流的下游区域比较显著。 相似文献
2.
分布热源作用下裂隙岩体渗流-传热的拉氏变换-格林函数半解析计算方法 总被引：1，自引：0，他引：1
以裂隙岩体高放射性核废物地下处置库性能评估为目标,提出了分布热源作用下单裂隙岩体渗流-传热的简化概念模型、控制微分方程和拉氏变换-格林函数半解析法,为进一步采用半解析法计算分布热源作用下多裂隙岩体的渗流-传热问题奠定了基础。针对单裂隙岩体的渗流-传热问题,建立考虑岩石内热源和二维热传导的控制微分方程,利用拉氏变换域微分方程的基本解建立格林函数积分方程,采用解析法处理其中的奇点,通过数值积分和拉氏数值逆变换求解,计算任意时刻裂隙水和岩石的温度分布。通过算例,与基于岩石一维热传导假定的解析解进行了对比,并计算分析了分布热源作用下单裂隙岩体的渗流-传热特征及其对裂隙开度、岩石热传导系数和热流集度的敏感度。算例表明,（1）就裂隙水温度而言,由于考虑了岩石的二维热传导,拉氏变换-格林函数半解析解小于基于岩石一维热传导假定的解析解;（2）裂隙水温度和岩石温度对裂隙开度和热流集度的敏感度较大,对岩石热传导系数的敏感度较小。 相似文献
3.
米尺度裂隙岩体模型水流-传热试验的数值模拟分析 总被引：1，自引：0，他引：1
为了研究高放射性核废物地下处置库近场的水流-传热耦合问题,采用国内高放废物地下处置库预选场址--甘肃北山地区的花岗岩石块体,加工组合成米尺度的规则裂隙岩体模型,设置边界热源和裂隙水流,试验模拟裂隙水水流与传热之间的相互作用。作为该室内模型试验的前期理论研究,采用等效孔隙介质数值模型,着重分析了裂隙开度、裂隙流量和热源功率对流场和温度场的影响。在设定条件下,计算分析表明：热传导和裂隙水水流由热源作用初期的不耦合很快转化为耦合;不流动的裂隙水主要表现为热存储和热传导,而流动的裂隙水还引起流动传热和水与岩石之间的对流换热,使岩体温度场明显不同于单纯热传导的情况;如果保持裂隙水流量不变,则裂隙开度的变化对水流-传热影响不大;如果保持裂隙水流速不变,则裂隙开度的变化对水流-传热影响显著;热源功率越大,通过裂隙水的热流量越大,裂隙水压强越大,而当温度超过100 ℃时,裂隙水会因汽化而压强显著增大;加热7 d时,热量的输入和输出几乎相等,裂隙水流带走的热量接近热源供给的热量,模型系统基本达到了热平衡。 相似文献
4.
选取高放射核废物处置库重要预选场区甘肃北山地区的花岗岩,制作750 mm（宽）×300 mm（厚）×1 000 mm（高）的稀疏不规则裂隙岩体模型,该模型由18块花岗岩和竖向与斜向各两条裂隙组成,在裂隙及岩石内部埋置温度传感器、水压计、直角应变花,并在模型一侧设置局部热源,研究热源温度和裂隙水流速对岩石温度和应力的影响。结果表明,竖裂隙水主要从顶部进水口流向底部出水口,斜裂隙水主要从侧部进水口流向侧部出水口,竖裂隙与斜裂隙在交汇处存在微小流量交换;由于热源处在两条斜裂隙进水口之间,并且斜裂隙长度小于竖裂隙,岩石热传导与斜裂隙水流对岩石温度分布起控制作用,竖裂隙水流对岩石横向热传导起阻滞作用;由于热传导和水流传热的不规则性,上层岩石形成从左向右为主的传热路径,中层和下层岩石形成从上向下为主的传热路径;由于上、下层岩石温度梯度较小,岩石收缩受热拉应力,而中层岩石温度梯度较大,岩石膨胀受热压应力,大主应力的方向大致垂直于斜裂隙面与竖裂隙面的交线,岩石应力增量随斜平面方向的温度梯度增大而增大;热源温度越高,裂隙水流速越低,岩石温度越高、岩石应力越大,系统达到稳态需要的时间越长。 相似文献
5.
裂隙岩体渗流-传热耦合的复合单元模型 总被引：1，自引：0，他引：1
《岩土力学》2016,(1)
基于复合单元法建立了裂隙岩体渗流-传热耦合的复合单元模型。该模型前处理简便快捷,网格剖分不受限制,可依据裂隙的真实信息自动将其离散在单元内。其次,采用交叉迭代算法,对裂隙岩体的渗流场和温度场进行耦合分析,耦合算法不仅考虑了温度对流体运动黏度的影响,而且可计算裂隙中流体与相邻岩块间渗流-传热过程以及两者间的渗流量和热量交换。通过与已有近似解析解相比较,验证了复合单元耦合算法的可靠性。算例分析表明,渗流-传热耦合作用对裂隙岩体的渗流场和温度场均有一定的影响。分析了不同岩块热传导系数和裂隙开度对热能提取效率的影响,结果显示,岩块热传导系数越大、裂隙开度越大,低温流体从高温岩块中吸取的热能会较多,出口处流体温度下降得较快。 相似文献
6.
基于复合单元法,结合三维热传导-对流方程和“充填模型”,提出了裂隙岩体不稳定温度场的复合单元模型。该模型前处理简便快捷,计算网格生成时无需考虑裂隙的存在,网格剖分不受限制,随后利用复合单元前处理程序,依据裂隙的位置和方位将其自动离散在单元内。对常规热传导-对流方程进行自伴随性调整,应用变分原理,推导出裂隙岩体不稳定温度场的复合单元算法,该算法可分别计算出岩块子单元和裂隙的温度值,且可真实反映裂隙中水流与相邻岩块间的热能量交换规律。将复合单元数值模型计算的不稳定温度场结果与相应的实测数据进行对比分析可知,数值计算结果与实测数据基本一致,验证了裂隙岩体不稳定温度场复合单元算法的可靠性与有效性。算例分析表明,裂隙中水流与相邻岩块间有明显的热传导和热对流作用。 相似文献
7.
《铀矿地质》2020,(4)
为研究正交与非正交稀疏裂隙岩体水流-传热对温度的影响,采用三维离散单元法程序(3 Dimension Distinct Element Code,简称3DEC)分析正交与非正交稀疏裂隙岩体水流-传热对温度场影响规律。研究表明:正交裂隙水流显著改变了岩体温度场;由于交叉裂隙水流动传热,致使裂隙水流动区域等温线出现明显的断续态,从瞬态到稳态,岩体温度升高,温度梯度逐渐减小。对比正交裂隙岩体模型,当横裂隙向下倾斜,纵裂隙靠近热源时,裂隙岩体温度升高;对比非正交裂隙岩体模型,当横裂隙倾斜度不变,纵裂隙远离热源时,裂隙岩体温度降低。非正交裂隙出水口水温高于正交裂隙出水口水温,纵裂隙靠近热源时裂隙出水口水温高于纵裂隙远离热源时裂隙出水口水温;非正交裂隙岩体模型比正交裂隙岩体模型达到稳态所需要的时间更短,纵裂隙靠近热源时的模型比纵裂隙远离热源时的模型达到稳态所需时间更长。 相似文献
8.
填砂裂隙岩体渗流传热模型试验与数值模拟 总被引：1，自引：0，他引：1
选取中国高放射核废物地下处置库重要预选场区--甘肃北山地区的花岗岩,加工组合成规则裂隙岩体,将垂直裂隙用粒径为0.5～0.63 mm的砂土填充,进行了裂隙水渗流传热试验;对模型试验进行了数值模拟,进而计算分析了热源温度、裂隙水流速和裂隙开度变化对裂隙岩体模型稳态温度场的影响。模型试验表明,当热源温度维持在120 ℃时,裂隙水仍无相变,裂隙岩体模型稳态温度场分布规律与热源温度为95 ℃时一致;热源温度越高,热源的水平影响距离越大,模型达到稳态需要的时间越长;裂隙填砂加强了裂隙两侧岩石之间的热传导,热源的水平影响距离和模型到达稳态需要的时间均明显大于无填充裂隙岩体模型的情况。模型试验得到的岩体模型温度场与数值计算得到的岩体模型温度场规律一致。试验过程中裂隙岩体模型在边界上存在一些热量散失,无法与数值计算中的绝热边界条件等同,致使试验数据低于数值计算值,并且热源温度越高,两者之间的差异越大。模型试验和数值计算均表明,邻近热源侧的裂隙水渗流对模型的温度场分布起控制作用,而远离热源侧的裂隙水渗流则主要影响该侧的边界温度和模型达到稳态所需要的时间。数值参数敏感性分析表明,裂隙水流速与裂隙开度越大,裂隙水对水平传热的阻滞作用越明显。 相似文献
9.
核废物地质处置、地热开发、石油开采等工程领域都可能涉及稀疏裂隙岩体中的水流-传热过程。现有的裂隙岩体水流-传热理论模型和计算方法基本上都是以平行光滑壁面裂隙模型为基础的,没有考虑裂隙的壁面局部接触对水流、水-岩热交换以及岩体传热的影响。针对粗糙壁面裂隙水流过程,阐述了基于Stokes方程的Reynolds润滑方程及Hele-Shaw裂隙模型,采用MATLAB软件中的PDE工具求解,并与Walsh的等效水力开度公式进行对比;分析壁面局部接触裂隙水流-传热与填充裂隙水流-传热的相似性,提出了瞬时局部热平衡假设的适用条件,并在裂隙局部接触体传热满足Biot数条件的前提下,计算分析裂隙局部接触体与水流之间的局部热平衡时间及其影响因素;在裂隙局部接触体与水流之间满足瞬时热平衡假设的前提下,利用填充裂隙水流-传热的解析解,计算了壁面局部接触裂隙水及两侧岩石的温度分布,并分析了裂隙局部接触面积率、裂隙开度、裂隙水平均流速对岩石温度和裂隙水温度的影响特征,结果表明：（1）在设定条件下,由于裂隙局部接触体与裂隙水流之间的热交换,裂隙水流对其两侧岩石温度的影响范围随接触面积率的增大而减小,裂隙两侧岩石对裂隙水流温度的影响程度随接触面积率的增大而增大;（2）裂隙开度和裂隙水流速对岩石温度和裂隙水温度的影响方式的影响是一致的,即由于裂隙水流量随裂隙开度和裂隙水流速的增大而增大,裂隙水流对其两侧岩石温度的影响范围随裂隙开度和裂隙水流速的增大而增大,裂隙两侧岩石对裂隙水流温度的影响程度随裂隙开度和裂隙水流速的增大而减小。 相似文献
10.
裂隙岩体因含有发育程度不同的裂隙、节理和断层等不连续面,致其渗透性具有各向异性、不连续性等特点,因此传统的有限元法对分布密集的裂隙岩体渗流场求解有一定的难度。本文提出了采用无单元Glaerkin法求解有自由面裂隙渗流问题,并推导了无单元法求解渗流场的基本方程和积分格式,给出了应用罚函数法处理渗流边界条件和自由面处理方法。采用IDL语言编制了二维无单元法计算软件LIDAREFM。文中以北京怀柔桥梓镇某裂隙岩体边坡渗流场计算为例,研究了复杂裂隙共同作用下渗流场特性和自由面分布,讨论了不同开度、不同连通程度的裂隙对渗流场的影响。研究结果表明:无单元法可以较好地解决有密集裂隙的岩体渗流场的求解问题,实现了裂隙处结点任意加密以及积分网格的独立布置,避免了对有自由面和裂隙穿越的子域的重新处理,简化了渗流问题的求解过程。 相似文献
11.
A semi‐analytical approach is developed for modeling 3D heat transfer in sparsely fractured rocks with prescribed water flow and heat source. The governing differential equations are formulated, and the corresponding integral equations over the fracture faces and the distributed heat source are established in the Laplace transformed domain using the Green function method with local systems of coordinates. The algebraic equations of the Laplace transformed temperatures of water in the fractures are formed by dividing the integrals into elemental ones; in particular, the fracture faces are discretized into rectangular elements, over which the integrations are carried out either analytically for singular integrals when the base point is involved or numerically for regular integrals when otherwise. The solutions of the algebraic equations are inverted numerically to obtain the real‐time temperatures of water in the fractures, which may be employed to calculate the temperatures at prescribed locations of the rock matrix. Three example calculations are presented to illustrate the workability of the developed approach. The calculations found that water flux in the fractures may decrease the rate of temperature rise in regions close to the distributed heat source and increase the rate of temperature rise in regions downstream away from the distributed heat source and that the temperature distribution and evolvement in a sparsely fractured rock mass may be significantly influenced by water flow exchange at intersection of fractures. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
12.
Poro‐mechanical and thermo‐mechanical processes change the fracture aperture and thus affect the water flow pattern in the fracture during the cold water injection into enhanced geothermal systems (EGS). In addition, the stresses generated by these processes contribute to the phenomenon of reservoir seismicity. In this paper, we present a three‐dimensional (3D) partially coupled poro‐thermoelastic model to investigate the poroelastic and thermoelastic effects of cold water injection in EGS. In the model, the lubrication fluid flow and the convective heat transfer in the fracture are modeled by the finite element method, while the pore fluid diffusion and heat conductive transfer in the reservoir matrix are assumed to be 3D and modeled by the boundary integral equation method without the need to discretize the reservoir. The stresses at the fracture surface and in the reservoir matrix are obtained from the numerical model and can be used to assess the variation of in situ stress and induced seismicty with injection/extraction. Application of the model shows that rock cooling induces large tensile stresses and increases fracture conductivity, whereas the rock dilation caused by fluid leakoff decreases fracture aperture and increases compressive total stresses around the injection zone. However, increases in pore pressure reduce the effective stresses and can contribute to rock failure, fracture slip, and microseismic activity. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
13.
Pushpendra Sharma Amit Kumar Verma Anil Negi Manish Kumar Jha Pradeep Gautam 《Arabian Journal of Geosciences》2018,11(15):431
Thermoelastic deformation of rock significantly affects the stability of rock slope because thermoelastic strains may cause fracture propagation under favorable condition of failure. Rock slope stability depends on the balance between shear stress and shear resistance along the plane of weakness. Due to warming of rock slopes by heat transfer phenomena, viz. conduction and convection, considerable change in induced stresses (normal and shear) and resistance takes place which further causes instability in rock slope. In this paper, a two-dimensional finite element model has been used to simulate the stability of jointed rock slope containing crack in its upper surface. Four different cases have been simulated on the basis of infilling material (air, water, ice, water and ice) in the crack. Stability of rock slope is examined in terms of shear displacement and factor of safety for different thermal conditions of slope surface. A comparative study has been done for the four cases of infilling material in the crack. The various affecting parameters, viz. shear displacement, factor of safety, shear strength along the joint, and different surface temperature conditions, are illustrated by means of graphs. It has been found that the values of horizontal and vertical displacements are in the range of millimeters. The maximum values of horizontal and vertical displacements are 2.17 mm. Moreover, the maximum values of vertical compressive and tensile stresses are 15.4 MPa and 4.45 MPa respectively for the said four cases. According to the infilling material in the crack, the stability of the rock slope for the given geometry of slope is found in the following order: crack filled with ice < crack filled with ice and water < crack filled with water < empty crack. Validations of numerical results have been done from previous studies, and it has been found that the trends of normal stress, shear strength, and shear displacement along the joint are well matched. 相似文献
14.
为研究填充裂隙水流速度对岩体温度和应力的影响,选取甘肃北山地区的花岗岩,制作了稀疏非正交裂隙岩体模
型,采用河砂填充裂隙后进行模型试验;并对模型试验进行离散元数值模拟,分析了模型试验所测岩体温度和热应力与数
值模拟结果的差异和原因。结果显示：填砂裂隙强化了裂隙介质的热导能力,无填充时岩体温度和热应力比填砂时低;模
型试验和离散元模拟均表明,岩体温度和岩体应力随裂隙水流速度增大而减小,但是系统达到稳态所需要的时间变短;模
型试验中斜裂隙水流对温度场起主要作用;由于现有3DEC软件不能考虑水的热物性参数随温度的变化,进而产生自然对
流换热,斜裂隙水流和靠近热源侧的竖裂隙水流对温度场起主要作用。 相似文献
15.
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. 相似文献
16.
The Manh Vu Jean Sulem Didier Subrin Nathalie Monin 《Rock Mechanics and Rock Engineering》2013,46(2):213-229
A semi-analytical solution based on the transfer matrix technique is proposed to analyze the stresses and displacements in a two-dimensional circular opening excavated in transversely isotropic formation with non-linear behavior. A non-isotropic far field can be accounted for and the process of excavation is simulated by progressive reduction of the internal radial stress. A hyperbolic stress–strain law is proposed to take into account the non-linear behavior of the rock. The model contains seven independent parameters corresponding to the five elastic constants of an elastic material with transverse isotropy and to the friction coefficient and cohesion along the parallel joints (weakness planes). This approach is based on the discretization of the space into concentric rings. It requires the establishment of elementary solutions corresponding to the stress and displacement fields inside each ring for given conditions at its boundaries. These solutions, based on complex variable theory, are obtained in the form of infinite series. The appropriate number of terms to be kept for acceptable approximation is discussed. This non-linear model is applied to back analyze the convergence measurements of Saint-Martin-la-Porte access gallery. Short-term and long-term ground parameters are evaluated. 相似文献