岩石裂隙压剪变形破坏与非线性渗流特性

作     者：李博 汪佳飞 刘日成 伍法权 LI Bo;WANG Jiafei;LIU Richeng;WU Faquan

作者机构：绍兴文理学院土木工程学院浙江绍兴312000 中国矿业大学深部岩土力学与地下工程国家重点实验室江苏徐州221116 

出 版 物：《工程科学与技术》 (Advanced Engineering Sciences)

年 卷 期：2021年第53卷第6期

页      面：103-112页

核心收录：

学科分类：08[工学] 080104[工学-工程力学] 0801[工学-力学（可授工学、理学学位）] 

基　　金：浙江省自然科学基金项目(LR19E090001) 国家自然科学基金项目(42077252 51979272 41831290) 中央高校基本科研业务专项(FRF–BD–20–01A)。 

主　　题：单裂隙 法向应力 剪切位移 裂隙变形 非线性渗流 

摘      要：探索岩体中流体的流动机理对解决隧道开挖和地下工程施工过程中的突涌水问题具有重要意义,定量描述岩石单裂隙受压剪共同作用下的渗流特性是理解裂隙渗流机理的基础,但已有研究尚未建立应力、位移、开度、惯性系数等关键控制参数之间的定量关系。本文制作了含单裂隙的花岗岩试样,采用固定剪切位移后逐渐增加法向应力的方式开展变水头渗流试验,得到裂隙渗透率的演化规律;通过高精度3维轮廓仪获取裂隙表面的形貌数据,在自主开发的基于接触力学变分原理框架的裂隙受压变形计算程序的基础上,表征不同剪切错位和不同应力条件下裂隙面的变形特性和内部空腔结构演化规律;提取裂隙开度分布数据,在COMSOL软件中求解Navier–Stocks方程,实施不同剪切位移和受压状态下的裂隙非线性渗流数值模拟;定量分析剪切位移、法向应力、空腔几何特征与非线性渗流控制参数之间的关系。结果表明:试验得到的裂隙表面破坏区域与数值模拟结果基本吻合,验证了裂隙变形计算程序的可靠性。法向应力、剪切位移分别与力学开度呈幂函数递减和指数函数递增的关系,且剪切位移的增加会导致裂隙内的接触面更集中。Forchheimer方程中的惯性系数B和临界水力梯度J_(c)可以定量刻画裂隙的非线性渗流特性,B和J_(c)与剪切位移之间呈幂函数递减关系,且随着剪切位移的增大,J_(c)和B的增加速率和波动范围逐渐降低:当剪切位移从2 mm增加到8 mm时,J_(c)的波动范围从6.10×10^(–3)减少到1.20×10^(–3),降低了80.32%;B的波动范围从2.97×10^(14) Pa·s^(2)·m^(–7)减少到2.43×10^(13) Pa·s^(2)·m^(–7),降低了91.28%。裂隙开度的相对标准偏差RSD与惯性系数B和临界水力梯度J_(c)之间均存在相似的幂函数关系,并据此建立了计算渗流从线性转换为非线性临界点的预测公式。