›› 2011, Vol. 32 ›› Issue (2): 333-340.

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

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

项彦勇,郭家奇   

  1. 北京交通大学 土木建筑学院,北京 100044
  • 收稿日期:2009-09-25 出版日期:2011-02-10 发布日期:2011-02-16
  • 作者简介:项彦勇,男,1959年生,博士,教授,博士生导师,从事岩土工程方面的教学与研究。
  • 基金资助:

    国家自然科学基金资助(No. 50778014)。

A Laplace transform and Green function method for calculation of water flow and heat transfer in fractured rocks

XIANG Yan-yong, GUO Jia-qi   

  1. School of Civil Engineering, Beijing Jiaotong University, Beijing 100044, China
  • Received:2009-09-25 Online:2011-02-10 Published:2011-02-16

摘要:

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

关键词: 核废物地下处置库, 裂隙岩体, 渗流, 传热, 格林函数, 离散拉氏变换

Abstract:

For performance assessment of nuclear waste repositories in fractured rocks, a simplified conceptual model and mathematical formulation is proposed; and a Laplace transform and Green function semi-analytical method is developed for calculation of flow and heat transfer in single-fracture rock with distributed heat source. Taking into account of distributed heat source and 2D thermal conduction in single-fracture rock of infinite extent, the mathematical model is formulated and solved by using a Green function approach, in which a fundamental solution of the governing differential equations after Laplace transform is employed. The singularities in the integral equation are handled through analytical integration; and a numerical procedure is developed to solve the transient temperature distributions in fracture water and rock matrix. Numerical examples are provided for illustration of the proposed method with comparison of an analytical solution based on 1D rock thermal conduction, as well as for features of flow and heat transfer in single-fracture rock and the sensitivities to fracture aperture, rock thermal conductivity and heat source intensity. The calculations indicate: (1) The temperature of water in the fracture calculated from using the semi-analytical method is lower than that calculated the analytical solution, due to the fact that the former method takes account 2D thermal conduction in the rock matrix, whereas the latter assumes 1D conduction. (2) Temperatures in fracture water and the rock matrix are more sensitive, in a relative sense, to the fracture aperture and the heat source intensity than to the thermal conductivity of the rock matrix. The proposed model and solution method may serve, among other possible applications, as a foundation for semi-analytical calculation of flow and heat transfer in multiple fracture rocks with distributed heat sources.

Key words: nuclear waste repository, fractured rock, groundwater flow, heat transfer, Green function, discrete Laplace transform

中图分类号: 

  • TU 452
[1] 苏永华, 李诚诚. 强降雨下基于Green-Ampt模型的边坡稳定性分析[J]. 岩土力学, 2020, 41(2): 389-398.
[2] 刘成禹, 陈博文, 罗洪林, 阮家椿, . 满流条件下管道破损诱发渗流侵蚀的试验研究[J]. 岩土力学, 2020, 41(1): 1-10.
[3] 刘忠玉, 夏洋洋, 张家超, 朱新牧. 考虑Hansbo渗流的饱和黏土 一维弹黏塑性固结分析[J]. 岩土力学, 2020, 41(1): 11-22.
[4] 夏才初, 喻强锋, 钱 鑫, 桂 洋, 庄小清. 常法向刚度条件下岩石节理剪切−渗 流特性试验研究[J]. 岩土力学, 2020, 41(1): 57-66.
[5] 肖瑶, 邓华锋, 李建林, 支永艳, 熊雨. 长期浸泡作用下灌浆加固裂隙岩体劣化效应[J]. 岩土力学, 2019, 40(S1): 143-151.
[6] 支永艳, 邓华锋, 肖瑶, 段玲玲, 蔡佳, 李建林. 微生物灌浆加固裂隙岩体的渗流特性分析[J]. 岩土力学, 2019, 40(S1): 237-244.
[7] 许江, 邬君宇, 刘义鑫, 雷娇, . 不同充填度下岩体剪切−渗流耦合试验研究[J]. 岩土力学, 2019, 40(9): 3416-3424.
[8] 刘伟俊, 张晋勋, 单仁亮, 杨昊, 梁辰, . 渗流作用下北京砂卵石地层多排管局部 水平冻结体温度场试验[J]. 岩土力学, 2019, 40(9): 3425-3434.
[9] 李博, 黄嘉伦, 钟振, 邹良超, . 三维交叉裂隙渗流传质特性数值模拟[J]. 岩土力学, 2019, 40(9): 3670-3768.
[10] 张强, 李小春, 周英博, 石露, 白冰, . 高压孔隙CO2/水作用下完整四川三叠系 砂岩剪切特性的试验研究[J]. 岩土力学, 2019, 40(8): 3028-3036.
[11] 王鹏飞, 谭文辉, 马学文, 李子建, 刘景军, 武洋帆, . 不同粗糙度和隙宽贯通充填裂隙 渗流特性试验研究[J]. 岩土力学, 2019, 40(8): 3062-3070.
[12] 张天军, 庞明坤, 蒋兴科, 彭文清, 纪翔, . 负压对抽采钻孔孔周煤体瓦斯渗流特性的影响[J]. 岩土力学, 2019, 40(7): 2517-2524.
[13] 曹洪, 胡瑶, 骆冠勇. 滤管两端均不在含水层层面的承压不 完整井近似计算方法研究[J]. 岩土力学, 2019, 40(7): 2774-2780.
[14] 刘忠玉, 崔鹏陆, 郑占垒, 夏洋洋, 张家超. 基于非牛顿指数渗流和分数阶Merchant模型的 一维流变固结分析[J]. 岩土力学, 2019, 40(6): 2029-2038.
[15] 刘 健, 陈 亮, 王春萍, 马利科, 王 驹. 一种非稳态气体渗流条件下岩石渗透特性 参数计算方法及应用[J]. 岩土力学, 2019, 40(5): 1721-1730.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 黄平路,陈从新,肖国峰,林 健. 复杂地质条件下矿山地下开采地表变形规律的研究[J]. , 2009, 30(10): 3020 -3024 .
[2] 孙 钧,戚玉亮. 隧道围岩稳定性正算反演分析研究 ——以厦门海底隧道穿越风化深槽施工安全监控为例介绍[J]. , 2010, 31(8): 2353 -2360 .
[3] 侯兴民,孔令召,陈建立. 地基阻尼比测试的一种计算方法[J]. , 2010, 31(9): 2995 -2999 .
[4] 张建刚,何 川,杨 征. 大断面宽幅盾构管片三维内力分布分析[J]. , 2009, 30(7): 2058 -2062 .
[5] 陈 芳,秦 昊. 细观尺度下岩石沿晶断裂应力强度因子计算研究[J]. , 2011, 32(3): 941 -945 .
[6] 陈国庆 ,苏国韶 ,江 权 ,李天斌. 大型地下洞室高边墙的拉应变准则研究[J]. , 2011, 32(S1): 603 -0608 .
[7] 文海家,张永兴,陈 云. 基于三维地质信息模型的边坡风险分析[J]. , 2009, 30(S2): 367 -370 .
[8] 彭林军 ,赵晓东 ,李术才. 深部开采地表沉陷规律模拟研究[J]. , 2011, 32(6): 1910 -1914 .
[9] 杨志全 ,侯克鹏 ,郭婷婷 ,马 秋. 黏度时变性宾汉体浆液的柱-半球形渗透注浆机制研究[J]. , 2011, 32(9): 2697 -2703 .
[10] 张忠苗,谢志专,赵玉勃,李慧明. 高临界循环荷载水平下钻孔桩的抗拔性状分析[J]. , 2012, 33(2): 343 -348 .