渗流场与温度场耦合分析的岩体裂隙网络数学模型.doc

PAGE 6 岩体裂隙网络渗流-温度耦合机理分析与模型建立 许增光1，柴军瑞1,2，程汉鼎1 ( 1. 西安理工大学 水利水电学院,陕西 西安 710048； 2. 三峡大学 土木水电学院,湖北 宜昌 443002) 1 引 言 近年来，由于新能源、新技术的发展，国内外在裂隙岩体基础上开辟了许多新的研究领域。核废料处理，我国目前正致力于核电建设，然而高放射性核废料地下深埋处理已成为举世关注的岩石力学难题。地热开发领域，在地表以下2000米~4000米的部位，埋藏有能量丰富的热干岩。然而我国对这方面的基础研究和现场试验才刚刚开始。水电开发领域，我国的水电资源开发率仅为16%左右，发展潜力巨大。但今后开发重点是西南、西北地区，这些地区的特点是地质环境复杂，山高谷深、岩溶发育、深厚覆盖层等 [1~3]。随着这些新领域的发展，越来越多的岩土工程问题随之而出现，尤其是岩体裂隙网络渗流场与温度场的耦合分析，现在还处于机理研究阶段[4~5]。 本文基于渗流温度基本理论和相互影响关系之上，分析岩体裂隙网络渗流场、水流温度场和岩块温度场之间的耦合机理，建立三者的耦合模型。 2 基本假设 (1) 假设岩块不透水，渗透水流只沿岩体裂隙运动； (2) 不考虑由于温度不均引起水的密度变化而造成的自然热对流作用，忽略热机械弥散作用和热辐射作用[6]； (3) 温度的变化不引起水的相变（保持液态不变）。 3 耦合水流温度场的渗流场分析 3.1 水流温度对渗流场的影响机理分析 由于假设岩块不透水，因此渗流场分析区域为由裂隙组成的岩体裂隙网络系统。水流温度对渗流场的影响作用包括以下两个方面。 （1）水温对渗透系数的影响 由于水流运动粘滞系数是温度的函数，一般采用以下经验公式[7~8]： （1） 式中 为水的运动粘滞系数，单位为cm2/s；为水温，单位为℃。 将式（1）代入裂隙水流渗透系数中可得： （2） 式中 为渗透系数；为水的动力粘滞系数；为水的容重；为隙宽。 （2）温度梯度引起水的流动 由温度梯度而引起的水流流量为 [9]： （3） 式中 为温度梯度引起的水流流量，为温差作用下的水流扩散率[10]。 3.2 耦合水流温度场的渗流场数学模型 由3.1分析可得耦合水流温度场的岩体裂隙网络渗流场数学模型为： （4） 式中 为贮存系数；为渗流源汇项。 4 耦合渗流场与岩块温度场的水流温度场分析 图1 热量运移单元体示意图Fig.1 Element sketch of heat transport对于水流温度场而言，其研究区域为水流所在区域，即由裂隙组成的裂隙网络系统。从研究热动力学方程出发，选取直角坐标系，在裂隙水流中选取空间一无限小六面体，如图1。假设坐标轴方向与主渗透方向一致，六面体边长分别为，，，且与坐标轴平行，作为平衡单元体，则在时间内引起单元体内温度变化的作用主要有：热量随水质点一起运移的对流作用；水流的热传导作用 图1 热量运移单元体示意图 Fig.1 Element sketch of heat transport （1）对流作用 由于对流作用而引起的流入流出单元体的热量差为： （5） 式中 为水流温度；为水的比热；为水的密度；，，分别为水的渗透流速在三个方向的分量。 （2）水流的热传导作用 由热传导作用引起的流入和流出单元体的热量差为： （6） 式中 为水的热传导系数。 （3）岩块与水流之间的热量运移 假设岩块温度为，容易得出由于岩块与水流之间的热量交换而进入单元体内水流中的热量为： （7） 式中 为热交换系数。时热量进入水中，时水中的热量进入岩块。 另一方面，在时间内在单元体内温度变化所增加的热量为： （8） 根据能量守恒定律，可以得到耦合渗流场、岩块温度场的水流温度场数学模型： （9） 式中 为水温源汇项。 5 耦合水流温度场的岩块温度场分析 当岩块表面与水流接触时，通过岩块表面的热流密度，与岩块表面温度与水流温度之差成正比， （10） 式中 ，，为边界表面向外法线的方向余弦。 对于岩块温度场，其研究区域为由裂隙所分割开的岩块。将式（10）代入岩块温度场数学模型中[11]，可得耦合水流