岩土力学 ›› 2025, Vol. 46 ›› Issue (S1): 437-453.doi: 10.16285/j.rsm.2024.0151CSTR: 32223.14.j.rsm.2024.0151

• 数值分析 • 上一篇    下一篇

不同围压条件下花岗岩剪切过程渗透率动态演化规律数值模拟研究

李槟1, 2,申海萌1, 2,李琦1, 2,李霞颖1, 2   

  1. 1. 中国科学院武汉岩土力学研究所 岩土力学与工程安全全国重点试验室,湖北 武汉 430071;2. 中国科学院大学,北京 100049
  • 收稿日期:2024-01-26 接受日期:2024-06-03 出版日期:2025-08-08 发布日期:2025-08-28
  • 通讯作者: 李琦,男,1972年生,博士,研究员,主要从事二氧化碳地质利用与封存方面的研究。E-mail: qli@whrsm.ac.cn
  • 作者简介:李槟,男,2000年生,博士研究生,主要从事岩石力学数值模拟方法研究。E-mail: libin223@mails.ucas.ac.cn
  • 基金资助:
    国家自然科学基金联合基金重点项目(No.U23A20671)

A numerical simulation of dynamic evolution of permeability during granite shear process under different confining pressures

LI Bin1, 2, SHEN Hai-meng1, 2, LI Qi1, 2, LI Xia-ying1, 2   

  1. 1. State Key Laboratory of Geomechanics and Geotechnical Engineering Safety, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan, Hubei 430071, China; 2. University of Chinese Academy of Sciences, Beijing 100049, China
  • Received:2024-01-26 Accepted:2024-06-03 Online:2025-08-08 Published:2025-08-28
  • Supported by:
    This work was supported by the Joint Funds of the National Natural Science Foundation of China (U23A20671).

摘要: 裂隙渗透率的动态演化是工程岩体水力学特性中的一个关键科学问题,了解岩体渗透率动态演化机制,明确影响渗透率动态演化的机制,对工程设计和运行有重要意义。基于离散元和有限元方法,建立粗糙裂隙花岗岩的剪切渗流数值模拟,探究裂隙开度与渗透率在变围压条件下及剪切过程中的动态演化机制,揭示粗糙花岗岩裂隙渗透率在剪切过程中的动态演化规律。开展在1.9~20.0 MPa围压下粗糙裂隙样品的剪切和渗流过程数值模拟,监测粗糙裂隙样品的裂隙开度与渗透率演化行为。数值模拟与试验表现出相同的演化规律,剪切和围压对渗透率的演化有显著影响,围压的大小决定着渗透率的变化趋势。在1.9~10.0 MPa的围压条件下,岩石的渗透率先明显增大,然后在剪切作用后下降;当围压大于10.0 MPa时,裂隙渗透率仅呈下降趋势。剪切数值模拟结果显示,高围压限制了裂隙在剪切过程中的剪胀扩张现象,同时促进了岩石碎屑的形成,降低了裂隙的过水能力。系统提出了基于Barton公式的等效开度修正系数,并给出最终具有物理意义的渗透率计算方法;提出了裂隙等效开度和围压间的演化标准方程,在描述围压与不同裂隙等效开度之间演化关系时,能够较好地预测不同围压下初始、峰值和最小等效开度变化,进而用以计算渗透率。研究采用的剪切渗流数值计算方法为揭示裂隙渗透率动态演化机制提供了一种有效途径。

关键词: 裂隙, 围压, 渗透率, 开度, 数值计算

Abstract: The dynamic evolution of fracture permeability is a critical issue in understanding the hydraulic characteristics of engineering rock masses. Investigating the dynamic evolution mechanisms of rock mass permeability is crucial for engineering design and operation, as it clarifies the factors influencing this process. Using the discrete element method and finite element method, this study establishes a numerical simulation framework for shear seepage in rough granite fractures to investigate the dynamic evolution of fracture aperture and permeability under varying confining pressures during shearing. The study reveals the dynamic evolution law of fracture permeability in rough granite during the shear process. Numerical simulations were conducted to analyze the shear and seepage processes of rough fracture samples under confining pressures ranging from 1.9 to 20.0 MPa, monitoring the evolution of fracture aperture and permeability. The numerical simulation results align with experimental observations, indicating that shear processes and confining pressure conditions significantly influence permeability evolution. Additionally, the magnitude of confining pressure determines the trend of permeability changes. Under confining pressures of 1.9–10.0 MPa, permeability initially increases significantly but decreases after shearing. When confining pressure exceeds 10.0 MPa, fracture permeability exhibits a downward trend throughout. Shear numerical simulations reveal that high confining pressures restrict fracture shear expansion, promote rock debris formation, and reduce the fracture’s water-carrying capacity. This study systematically proposes an equivalent aperture correction coefficient based on the Barton formula, provides a physically meaningful permeability calculation method, and establishes a standard evolution equation linking fracture equivalent aperture and effective confining pressure. These formulations enable accurate predictions of initial, peak, and minimum equivalent apertures under varying confining pressures and facilitate permeability calculations. The shear seepage numerical simulation method offers an effective approach to elucidate the dynamic evolution mechanisms of fracture permeability.

Key words: fracture, confining pressure, permeability, aperture, numerical simulation

中图分类号: TU452
[1] 王宁博, 姚仰平, 刘林, 李翔宇, 毛安琪, 李宁, . 考虑围压效应的砂土统一硬化模型[J]. 岩土力学, 2025, 46(S1): 297-308.
[2] 苗日成, 唐贝, 祁飞, 江志安, 崔溦, . 随机裂隙岩体滚刀破岩过程离散元模拟研究[J]. 岩土力学, 2025, 46(S1): 541-552.
[3] 刘一鸣, 李振, 冯国瑞, 杨鹏, 白锦文, 黄炳雄, 李东, . 循环加卸载下裂隙砂岩声−热响应特征及前兆规律[J]. 岩土力学, 2025, 46(9): 2773-2791.
[4] 王昕琪, 冯子军, 陈正男, 高祺, 阴伟涛, 靳佩桦, 李玉彬, . 超临界水作用下花岗岩裂隙渗流特性演化规律研究[J]. 岩土力学, 2025, 46(9): 2847-2858.
[5] 王邺, 王述红, 张泽, 韩博文, 杨润生, . 裂隙及水幕影响下的水封洞库油品运移范围研究[J]. 岩土力学, 2025, 46(9): 2907-2928.
[6] 孙闯, 蒲云波, 敖云鹤, 陶琦, . 冻融饱水裂隙砂岩力学特性及细观破裂演化特征研究[J]. 岩土力学, 2025, 46(8): 2339-2352.
[7] 王炳文, 刘臣毅, 康明超, 李乾龙, 杨雷, 周森林, 钱垒. 考虑应力的非闭合裂隙岩体损伤机制及裂纹扩展研究[J]. 岩土力学, 2025, 46(8): 2409-2420.
[8] 高帅, 任富强, 常远. 预制单裂隙砂岩注浆过程中的浆液流动特性研究[J]. 岩土力学, 2025, 46(8): 2505-2515.
[9] 熊飞, 刘新荣, 刘文武, 钟祖良, 杨忠平, 王南云, 王浩, 薛熠, . 含深大裂隙陡崖式岩溶山体采动破坏与失稳机制[J]. 岩土力学, 2025, 46(8): 2516-2531.
[10] 雷瑞德, 顾清恒, 胡超, 贺培, 周林森, . 裂隙砂岩声发射信号特征及破裂前兆识别研究[J]. 岩土力学, 2025, 46(7): 2023-2038.
[11] 刘杰, 钟伦伟, 钟振, 张宪尚, 隆清明, 应鹏, . 恒定法向刚度边界条件下受剪岩石裂隙非线性渗流特性的数值模拟研究[J]. 岩土力学, 2025, 46(6): 1919-1933.
[12] 姜海波, 卢燕, 李琳, 张军, . 干湿-冻融作用输水渠道膨胀土强度特性及损伤演化规律研究[J]. 岩土力学, 2025, 46(5): 1356-1367.
[13] 欧阳淼, 张红日, 王桂尧, 邓人睿, 郭鸥, 汪磊, 高游, . 基于响应面法的生物基质改良膨胀土配比优化研究[J]. 岩土力学, 2025, 46(5): 1368-1378.
[14] 包卫星, 田磊, 吴谦, 黄志明, 张智勇. 伊犁黄土水分入渗特征及其对湿陷变形的影响[J]. 岩土力学, 2025, 46(5): 1379-1391.
[15] 张天军, 田嘉伟, 张磊, 庞明坤, 潘红宇, 孟伟, 贺绥男, . 循环载荷下破碎煤体渗透率及迂曲度演化研究[J]. 岩土力学, 2025, 46(5): 1409-1418.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!