›› 2014, Vol. 35 ›› Issue (10): 2845-2854.

• 基础理论与实验研究 • 上一篇    下一篇

壁面局部接触裂隙岩体水流-传热的一种理论模型及计算分析

项彦勇,任 鹏   

  1. 北京交通大学 土木工程学院,北京 100044
  • 收稿日期:2013-06-25 出版日期:2014-10-11 发布日期:2014-10-27
  • 作者简介:项彦勇,男,1959年生,博士,教授,从事岩土工程与地下工程方面的研究和教学工作。
  • 基金资助:
    高等学校博士学科点专项科研基金(No. 20120009110022);国家自然科学基金(No. 50778014,No. 51378055)。

A theoretical model and its calculation analysis of water flow and heat transfer in a fractured rock with local fracture wall asperity contacts

XIANG Yan-yong,REN Peng   

  1. School of Civil Engineering, Beijing Jiaotong University, Beijing 100044, China
  • Received:2013-06-25 Online:2014-10-11 Published:2014-10-27

摘要: 核废物地质处置、地热开发、石油开采等工程领域都可能涉及稀疏裂隙岩体中的水流-传热过程。现有的裂隙岩体水流-传热理论模型和计算方法基本上都是以平行光滑壁面裂隙模型为基础的,没有考虑裂隙的壁面局部接触对水流、水-岩热交换以及岩体传热的影响。针对粗糙壁面裂隙水流过程,阐述了基于Stokes方程的Reynolds润滑方程及Hele-Shaw裂隙模型,采用MATLAB软件中的PDE工具求解,并与Walsh的等效水力开度公式进行对比;分析壁面局部接触裂隙水流-传热与填充裂隙水流-传热的相似性,提出了瞬时局部热平衡假设的适用条件,并在裂隙局部接触体传热满足Biot数条件的前提下,计算分析裂隙局部接触体与水流之间的局部热平衡时间及其影响因素;在裂隙局部接触体与水流之间满足瞬时热平衡假设的前提下,利用填充裂隙水流-传热的解析解,计算了壁面局部接触裂隙水及两侧岩石的温度分布,并分析了裂隙局部接触面积率、裂隙开度、裂隙水平均流速对岩石温度和裂隙水温度的影响特征,结果表明:(1)在设定条件下,由于裂隙局部接触体与裂隙水流之间的热交换,裂隙水流对其两侧岩石温度的影响范围随接触面积率的增大而减小,裂隙两侧岩石对裂隙水流温度的影响程度随接触面积率的增大而增大;(2)裂隙开度和裂隙水流速对岩石温度和裂隙水温度的影响方式的影响是一致的,即由于裂隙水流量随裂隙开度和裂隙水流速的增大而增大,裂隙水流对其两侧岩石温度的影响范围随裂隙开度和裂隙水流速的增大而增大,裂隙两侧岩石对裂隙水流温度的影响程度随裂隙开度和裂隙水流速的增大而减小。

关键词: 壁面局部接触裂隙, Hele-Shaw模型, 水流-传热, 瞬时局部热平衡, 参数分析

Abstract: Water flow and heat transfer processes in fractured rocks may be encountered in several engineering disciplines, such as geological disposal of radioactive waste, geothermal development, petroleum exploration, etc. Currently existent theoretical models and calculation methods for water flow and heat transfer in fractured rocks are almost exclusively based on the assumption of parallel smooth fracture walls, without considering the effects of local fracture wall asperity contacts on water flow, water-rock heat exchange and heat transfer. The following studies are presented in this paper. Firstly, dealing with water flow in a fracture with rough walls, the Reynolds lubrication equation and the Hele-Shaw model of a fracture with local wall asperity contacts based on the Stokes equations are expounded; solutions are sought by using the PDE tool of the MATLAB software, and compared with the Walsh formula of the effective hydraulic aperture of a rough walled fracture. Secondly, the similarity between the water flow and heat transfer in a fractured rock with local fracture wall asperity contacts and the water flow and heat transfer in a fractured rock with fills is analyzed; the condition for applicability of the instantaneous local thermal equilibrium is proposed; and the time and its influencing parameters for local thermal equilibrium between the local fracture wall asperity contacts and the flowing water in the fracture are computationally analyzed under the premise that the Biot number condition can be satisfied by the local fracture wall asperity contacts. Lastly, assuming instantaneous local thermal equilibrium between the local fracture wall asperity contacts and the flowing water in the fracture, the distributions of temperatures in the flowing water and the rock matrix are calculated using an existing analytical solution for water flow and heat transfer in a fractured rock with fills, and the effects of the ratio of local asperity contact area, the fracture aperture and the averaged water velocity on the temperature distributions are analyzed, which indicate, under the prescribed conditions, (1) that due to the heat exchange between the local fracture wall asperity contacts and the flowing water in the fracture, the range of influence of the flowing water on the temperature of the rock matrix decreases with the increase of the ratio of local asperity contact area, whereas the degree of influence of the rock matrix on the temperature of the flowing water in the fracture increases as the ratio of local asperity contact area increases; (2) the influences of the fracture aperture and the water velocity on the temperatures of the rock matrix and the flowing water in the fracture are consistent in that the range of influence of the flowing water on the temperature of the rock matrix increases with the increases of the fracture aperture and the water velocity, whereas the degree of influence of the rock matrix on the temperature of the flowing water decreases as the fracture aperture and the water velocity increase.

Key words: fracture with local wall asperity contacts, Hele-Shaw model, water flow and heat transfer, instantaneous local thermal equilibrium, parametric analysis

中图分类号: 

  • TU 45
[1] 杨鹏, 蒲诃夫, 宋丁豹. 竖井地基的大应变固结分析[J]. 岩土力学, 2019, 40(10): 4049-4056.
[2] 李 磊,蒋明镜,张伏光, . 深部岩石考虑残余强度时三轴试验离散元定量模拟及参数分析[J]. , 2018, 39(3): 1082-1090.
[3] 高俊义,项彦勇. 稀疏不规则裂隙岩体模型三维水流-传热对温度和应力影响的试验研究[J]. , 2016, 37(11): 3145-3155.
[4] 邱清文 ,张文杰 ,程泽海 . 湿润地区垃圾填埋场蒸发蒸腾覆盖层参数分析[J]. , 2012, 33(S1): 283-289.
[5] 靳建明 ,梁仕华 . 成层Gibson地基中单桩沉降的非线性分析[J]. , 2012, 33(6): 1857-1863.
[6] 朱 磊 ,龚晓南 ,邢 伟. 土钉支护基坑抗隆起稳定性计算方法研究[J]. , 2012, 33(1): 167-170.
[7] 周万欢 ,殷建华. 上覆压力和剪胀作用下土钉抗拔的有限元模拟[J]. , 2011, 32(S1): 691-0696.
[8] 林宇亮,杨果林,赵炼恒. 地震条件下挡墙后黏性土主动土压力研究[J]. , 2011, 32(8): 2479-2486.
[9] 卢 正 ,姚海林 ,罗海宁 ,杨 洋 ,杨明亮 . 双参数黏弹性地基上连续配筋混凝土路面振动参数分析[J]. , 2008, 29(8): 2177-2182.
[10] 赵 颖,梁 冰,薛 强,刘 磊,王永波. 土壤水环境中污染物运移双点吸附解吸动力学模型[J]. , 2007, 28(12): 2574-2578.
[11] 李成柱,周志芳. 抽水引致的地面沉降的参数敏感性分析[J]. , 2006, 27(S2): 81-85.
[12] 张 云. 修正剑桥模型参数对计算结果的影响[J]. , 2006, 27(3): 441-444.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 程 涛,晏克勤. 应力路径对地表变形特性影响的数值模拟[J]. , 2010, 31(2): 661 -666 .
[2] 关云飞,高 峰,赵维炳,俞 缙. ANSYS软件中修正剑桥模型的二次开发[J]. , 2010, 31(3): 976 -980 .
[3] 何先龙,赵立珍. 基于多重互相关函数分析剪切波速[J]. , 2010, 31(8): 2541 -2545 .
[4] 张玉成,杨光华,姜 燕,姚 捷,史永胜. 沉管隧道基槽爆破施工对既有堤岸稳定性影响的数值仿真分析[J]. , 2010, 31(S1): 349 -356 .
[5] 肖 琳,杨成奎,胡增辉,李晓昭,李 茉. 地铁隧道围岩内温度分布规律的模型试验及其热导率反算研究[J]. , 2010, 31(S2): 86 -91 .
[6] 陈敬虞,龚晓南,邓亚虹. 软黏土层一维有限应变固结的超静孔压消散研究[J]. , 2009, 30(1): 191 -195 .
[7] 潘 岳,张 勇,王志强. 煤与瓦斯突出中单个煤壳解体突出的突变理论分析[J]. , 2009, 30(3): 595 -602 .
[8] 付成华,周洪波,陈胜宏. 混凝土喷层支护节理岩体等效力学模型及其应用[J]. , 2009, 30(7): 1967 -1973 .
[9] 张春会,赵全胜. 基于ARCGIS的矿山开采沉陷灾害预警系统[J]. , 2009, 30(7): 2197 -2202 .
[10] 崔 溦,宋慧芳,张社荣,闫澍旺. 爆炸荷载作用下土中爆坑形成的数值模拟[J]. , 2011, 32(8): 2523 -2528 .