基于裂隙岩体温度场分布的等效隙宽反演.pdfVIP

基于裂隙岩体温度场分布的等效隙宽反演 白兰兰’ 陈建生 董海洲 陈亮 (河海大学岩土工程研究所。南京210098) 摘要：基于裂隙岩体的渗流场和温度场的相互作用机理，分析了平板单裂隙渗流作用下的温度场分布特征。从 流体力学、传热学的基本原理出发，根据现场温度场分布，推导出裂隙内的渗流量和等效隙宽，提出了通过温 度场求解等效水力隙宽的新方法，并利用FEMLAB模拟不同等效隙宽影响下裂隙岩体的温度场分布，为渗流 计算提供了新思路。最后，结合工程实例，验证了模型的正确性。 关键词：渗流场温度场对流换热等效渗透系数等效水力隙宽 1 引言 21，对于裂隙岩体，由于其天然的宏 目前研究渗流场与温度场相互作用时，通常假设介质均匀、连续[L 观和微观结构的非均质性、各向异性以及不连续性，使用这种等效模型显然不大合适。裂隙岩体地区的大裂 隙或裂隙密集区域一般存在强渗流，因此可以将小孔隙均匀分布的区域等效为连续的多孔介质；将大的裂隙 区域等效为裂隙结构介质。①当岩体孑L隙率较小，而裂隙较发达时，认为渗流主要发生在裂隙通道内，只考 虑集中裂隙通道对温度场的影响。②对于孔隙和裂隙都比较发达的岩体，认为孔隙均匀渗流与裂隙集中渗流 同时存在，研究渗流、管流耦合作用下的温度场分布。 单裂隙作为裂隙岩体的基本组成单位，研究其渗流特性以及与之相平衡的温度场分布，是研究复杂裂隙 岩体耦合作用的基础。在处理裂隙渗流时，一般根据实际的粗糙程度，引进平均隙宽，来计算其等效渗透系 数，现有的描述裂隙面粗糙性的方法有凸起高度表征法、节理粗糙度系数JRC表征法【3】和分数维表征法¨3。 但以上方法对于表面粗糙性的现场测试要求较高，参数获取较困难，而利用温度场来研究等效隙宽的方法至 今还未见到，本文从温度场与渗流场的相互作用出发，提出了利用温度场来推求等效隙宽的方法。 2裂隙岩体地区渗流场和温度场分布 有渗流存在的裂隙岩体地区，流动的地下水与其所接触的岩体表面之间将发生热量传递(即表面对流换 热)，一方面地下水不断流入，水与岩土体之间的温差(地下水体温度一般较岩土体温度低)致使岩土体与水 之间发生热量交换(以热传导为主)；另一方面，地下水渗流过程中存在热对流效应。上述两者共同作用， 最终使得岩土体温度场分布总体衰减，温度降低。存在较强集中渗漏的裂隙岩体，流体与固体壁面间的温度 重新分配是渗流影响下裂隙岩体温度场形成的主要因素。 温度场也会对渗流场产生一定的影响，裂隙岩体中温度场分布的差异性改变，新温度场内部各部分之间 存在温差，由温差引起的热应力作用下岩体结构特征的改变(裂隙同岩体中不连续结构面的张开度、粗糙度及 相互延伸交切状况)导致裂隙岩体渗透性能变化，进而影响围岩体中地下水渗流场的分布。 3单裂隙渗流下的温度场分析 首先针对单裂隙进行研究，假设单裂隙中流体为层流，对于管流和平板裂隙模型，根据现场物探资料， 确定裂隙的中心位置，按照平均隙宽建立流体流动传热模型。 3．1 裂隙剖面速度分布 3．1．1管流速度剖面 从目前研究较多的充分发展的管流出发，对于充分发展管流可得到圆管中的流速【6】： 材～。[1-(}2】 (1) ·180· 这是抛物线速度分布定律，也称Poiseuille Flow。其中“：一生生是中心线速度∥为动力粘滞系数，，． dX ‘4／a 为到圆管中心的距离，R为圆管半径。 ‘ 仍假设管流为层流，由达西定律矿=肜r知，，：：一生竺相当于上式中的一三萼．因此可考虑应用等效 。d工 yo工 渗透系数来研究裂隙通道中的管道流。由上述推导知：平均等效渗透系数