岩土力学 ›› 2023, Vol. 44 ›› Issue (9): 2757-2766.doi: 10.16285/j.rsm.2022.1524

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

循环剪切作用下三维粗糙裂隙非线性渗流特性数值模拟研究

蔚立元,杨瀚清,王晓琳,刘日成,王蓥森   

  1. 中国矿业大学 深地工程智能建造与健康运维全国重点实验室,江苏 徐州 221116
  • 收稿日期:2022-10-02 接受日期:2023-02-17 出版日期:2023-09-11 发布日期:2023-09-02
  • 通讯作者: 王晓琳,女,1997年生,博士研究生,主要从事单裂隙岩体剪切渗流与传热方面的研究工作。E-mail: wangxiaolin@cumt.edu.cn E-mail:yuliyuan@cumt.edu.cn
  • 作者简介:蔚立元,男,1982年生,博士,研究员,主要从事岩石力学与地下工程方面的教学与研究工作。
  • 基金资助:
    国家自然科学基金项目(No.51979272,No.52179118);江苏省自然科学基金项目(No.BK20211584);江苏省科技厅科技项目?碳达峰碳中和科技创新专项资金项目(No.BK20220025)

Numerical study on nonlinear hydraulic properties of three-dimensional rough joints under cyclic shear conditions

YU Li-yuan, YANG Han-qing, WANG Xiao-lin, LIU Ri-cheng, WANG Ying-sen   

  1. State Key Laboratory of Intelligent Construction and Healthy Operation and Maintenance of Deep Underground Engineering, China University of Mining and Technology, Xuzhou, Jiangsu 221116, China
  • Received:2022-10-02 Accepted:2023-02-17 Online:2023-09-11 Published:2023-09-02
  • Supported by:
    This work was supported by the National Natural Science Foundation of China (51979272, 52179118), the Natural Science Foundation of Jiangsu Province, China (BK20211584) and the Science and Technology Project of Jiangsu Provincial Department of Science and Technology-Carbon Emissions Peak and Carbon Neutrality Science and Technology Innovation Special Fund Project (BK20220025).

摘要: 研究了循环剪切过程中剪切位移和循环次数对三维粗糙裂隙非线性水力特性的影响。首先,基于天然粗糙裂隙面的高程数据生成三维裂隙面,并通过前人的研究结果确定循环剪切过程中的法向位移变化值,建立了16个开度空间模型,从而得到了不同剪切位移和循环次数下接触面积和开度分布情况,发现当循环次数一定时,随着剪切位移的增加,原本吻合的上下表面发生错动,平均开度和各向异性增大,接触面积减小;剪切位移一定时,随着循环次数的增加,上下表面起伏在剪应力作用下被磨平,接触面积增大,平均开度减小。将模型导入数值模拟软件进行计算得到流量数据,结果表明,水力梯度随流量变化规律符合Forchheimer定律。循环1次时,随着剪切位移从4 mm增加到10 mm,线性Fochheimer系数减小了88.3%,非线性Fochheimer系数减小了95.2%;循环2次时,随着剪切位移从4 mm增加到10 mm,线性Fochheimer系数减小了95.4%,非线性Fochheimer系数减小了99.7%,两系数减小的幅度随循环次数增大而增大。进一步分析发现,剪切位移的增加对裂隙渗透性有促进作用,循环次数增加对渗透性起抑制作用。此外,归一化渗透系数随雷诺数变化规律与Zimmerman等提出的预测公式相吻合,且本研究数值模拟计算结果与其公式的拟合系数趋近1,公式中的Fochheimer系数β 随剪切位移增大而减小,随循环次数增大而增大。

关键词: 循环剪切, 剪切渗流, 雷诺数, 归一化渗透系数

Abstract: In this study, the effects of shear displacement and the number of cycles on the nonlinear hydraulic characteristics of three-dimensional rough fracture during cyclic shear are numerically studied. First, a three-dimensional fracture surface is generated based on the elevation data of the natural rough fracture surface, and the normal displacement during the process of cyclic shear is determined through previous studies. Sixteen fracture aperture fields are established so as to obtain the contact area and aperture distribution under different shear displacements and cycles times. It is found that when the number of cycles is constant, the upper and lower surfaces of the original mesh will stagger with the increase of shear displacement, the average mechanical aperture and anisotropy increase, and the contact area decreases; when the shear displacement is constant, with the increase of the number of cycles, the upper and lower surface undulations are smoothed under the action of shear stress, the contact area increases, and the average mechanical aperture decreases. Finally, the model is imported into the numerical simulation software to calculate the flow rate data. The results show that the variation of hydraulic gradient with flow rate conforms to Forchheimer’s law. The linear Fochheimer coefficient decreases by 88.3% and the nonlinear Fochheimer coefficient decreases by 95.2% as the shear displacement increases from 4 mm to 10 mm under 1 cycle; the linear Fochheimer coefficient decreases by 95.4%, the nonlinear Fochheimer coefficient decreases by 99.7% as the shear displacement increases from 4 mm to 10 mm under 2 cycles, and the reduction range of the two coefficients increases with the increase of the number of cycles. Further analysis shows that the increase of shear displacement boosts the fracture permeability, while the increase in number of cycles inhibits the permeability. In addition, the variation of normalized transmissivity coefficient with Reynolds number conforms to the research results of Zimmerman et al., and the numerical simulation results of this study are fitted with the normalized transmissivity coefficient prediction formula proposed by Zimmerman et al. The fitting coefficient approaches to 1, where Fochheimer coefficient β decreases with the increase of shear displacement and increases with the increase of cycle times.

Key words: cyclic shear, shear fluid flow, Reynolds number, normalized transmissivity coefficient

中图分类号: 

  • TU411
[1] 秦悠, 杜歆宇, 马维嘉, 吴琪, 陈国兴, . 不同循环应力路径下饱和珊瑚砂超静孔压增长的改进应力模型[J]. 岩土力学, 2023, 44(6): 1729-1738.
[2] 尹敬涵, 崔臻, 盛谦, 陈健, 张茂础, . 循环荷载作用下岩石劈裂结构面剪切力学特性演化与影响因素研究[J]. 岩土力学, 2023, 44(1): 109-118.
[3] 孙杰豪, 郭保华, 田世轩, 程坦, . 峰前循环剪切作用下岩石节理剪切力学特性[J]. 岩土力学, 2022, 43(S2): 52-62.
[4] 侯娟, 邢行, 徐东, 陆向前, . 土工膜-土工布界面动力剪切特性试验研究[J]. 岩土力学, 2022, 43(2): 365-376.
[5] 王柳江, 刘斯宏, 赵志杰, 沈超敏, 鲁洋. 土工袋界面动力特性的循环直剪试验研究[J]. 岩土力学, 2021, 42(6): 1625-1634.
[6] 刘新荣, 许彬, 周小涵, 谢应坤, 何春梅, 黄俊辉, . 软弱层峰前循环剪切宏细观累积损伤机制研究[J]. 岩土力学, 2021, 42(5): 1291-1303.
[7] 芮圣洁, 国振, 王立忠, 周文杰, 李雨杰, . 钙质砂与钢界面循环剪切刚度与阻尼比的试验研究[J]. 岩土力学, 2020, 41(1): 78-86.
[8] 张强, 李小春, 周英博, 石露, 白冰, . 高压孔隙CO2/水作用下完整四川三叠系 砂岩剪切特性的试验研究[J]. 岩土力学, 2019, 40(8): 3028-3036.
[9] 许 江, 瞿佳美, 刘义鑫, 彭守建, 王 威, 吴善康, . 循环剪切荷载作用下充填物对结构面 剪切特性影响试验研究[J]. 岩土力学, 2019, 40(5): 1627-1637.
[10] 李书兆, 王忠畅, 贾 旭, 贺林林, . 软黏土中张紧式吸力锚循环承载力简化计算方法[J]. 岩土力学, 2019, 40(5): 1704-1712.
[11] 吴 野,王 胤,杨 庆. 考虑钙质砂细观颗粒形状影响的液体拖曳力系数试验[J]. , 2018, 39(9): 3203-3212.
[12] 熊 峰,孙 昊,姜清辉,叶祖洋,薛道锐,刘乳燕,. 粗糙岩石裂隙低速非线性渗流模型及试验验证[J]. , 2018, 39(9): 3294-3302.
[13] 刘飞禹,施 静,王 军,蔡袁强,. 三明治形加筋土筋-土界面动力剪切特性[J]. , 2018, 39(6): 1991-1998.
[14] 刘飞禹,王 攀,王 军,蔡袁强,. 颗粒粒径对格栅-土界面静、动力直剪特性的影响[J]. , 2017, 38(1): 150-156.
[15] 刘飞禹,王 攀,王 军,胡秀青,蔡袁强,. 筋-土界面循环剪切刚度和阻尼比的试验研究[J]. , 2016, 37(S1): 159-165.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 姚仰平,侯 伟. 土的基本力学特性及其弹塑性描述[J]. , 2009, 30(10): 2881 -2902 .
[2] 徐金明,羌培,张鹏飞. 粉质黏土图像的纹理特征分析[J]. , 2009, 30(10): 2903 -2907 .
[3] 向天兵,冯夏庭,陈炳瑞,江 权,张传庆. 三向应力状态下单结构面岩石试样破坏机制与真三轴试验研究[J]. , 2009, 30(10): 2908 -2916 .
[4] 石玉玲,门玉明,彭建兵,黄强兵,刘洪佳. 地裂缝对不同结构形式桥梁桥面的破坏试验研究[J]. , 2009, 30(10): 2917 -2922 .
[5] 夏栋舟,何益斌,刘建华. 土-结构动力相互作用体系阻尼及地震反应分析[J]. , 2009, 30(10): 2923 -2928 .
[6] 徐速超,冯夏庭,陈炳瑞. 矽卡岩单轴循环加卸载试验及声发射特性研究[J]. , 2009, 30(10): 2929 -2934 .
[7] 张力霆,齐清兰,魏静,霍倩,周国斌. 淤填黏土固结过程中孔隙比的变化规律[J]. , 2009, 30(10): 2935 -2939 .
[8] 张其一. 复合加载模式下地基失效机制研究[J]. , 2009, 30(10): 2940 -2944 .
[9] 易 俊,姜永东,鲜学福,罗 云,张 瑜. 声场促进煤层气渗流的应力-温度-渗流压力场的流固动态耦合模型[J]. , 2009, 30(10): 2945 -2949 .
[10] 陶干强,杨仕教,任凤玉. 崩落矿岩散粒体流动性能试验研究[J]. , 2009, 30(10): 2950 -2954 .