中国安全科学学报 ›› 2025, Vol. 35 ›› Issue (S2): 50-57.doi: 10.16265/j.cnki.issn1003-3033.2025.S2.0018

• 安全工程技术 • 上一篇    下一篇

基于双重渗透介质模型的土质边坡降雨入渗动态特征

孙长坤1,2(), 吴海飞3,4,**(), 张亦海3,4, 杨江林1, 元宇1,5, 袁云星3,4   

  1. 1 云南黄金矿业集团股份有限公司, 云南 昆明 650299
    2 北京科技大学 金属矿山开采安全与灾害防治全国重点实验室, 北京 100083
    3 中国安全生产科学研究院, 北京 100012
    4 中安国泰(北京)科技发展有限公司, 北京 102209
    5 鹤庆北衙矿业有限公司, 云南 鹤庆 671500
  • 收稿日期:2025-08-14 出版日期:2026-02-04
  • 通信作者:
    **吴海飞(1999—),男,安徽铜陵人,硕士,主要从事露天矿山安全监测、边坡稳定性分析等方面的工作。E-mail:
  • 作者简介:

    孙长坤 (1987—),男,吉林长春人,博士研究生,高级工程师,主要从事采矿工艺、岩石力学、露天矿山安全监测、边坡稳定性分析等方面的工作。E-mail:

    张亦海 高级工程师

    袁云星 工程师

  • 基金资助:
    国家重点研发计划重点专项(2022YFC3004601)

Study on dynamic rainfall infiltration characteristics of soil slopes based on a dual-permeability model

SUN Changkun1,2(), WU Haifei3,4,**(), ZHANG Yihai3,4, YANG Jianglin1, YUAN Yu1,5, YUAN Yunxing3,4   

  1. 1 Yunnan Gold Mining Group Co., Ltd., Kunming Yunnan 650299, China
    2 State Key Laboratory of Metal Mine Mining Safety and Disaster Prevention and Control, University of Science and Technology Beijing, Beijing 100083, China
    3 China Academy of Safety Science and Technology, Beijing 100012, China
    4 Cathay Safety Technology Co., Ltd., Beijing 102209, China
    5 Heqing Beiya Mining Co., Ltd., Heqing Yunnan 671500, China
  • Received:2025-08-14 Published:2026-02-04

摘要:

为提高滑坡预测预警能力、优化工程防灾减灾设计,基于双重渗透介质理论,构建考虑基质-裂隙耦合效应的非饱和渗流模型。模型充分考虑降雨强度、裂隙渗透系数、基质渗透系数以及土体初始含水率等因素的影响,以真实反映水分在基质与裂隙域之间的竞争、迁移与耦合作用。并基于改进的降雨边界,结合土体非饱和渗流特性,利用COMSOL Multiphysics有限元数值计算软件为平台,综合分析水分在基质与裂隙域之间的竞争与耦合作用,揭示双重渗透介质作用下均质土体的降雨入渗特征。结果表明:裂隙的存在显著加快雨水的初期入渗速率,并在一定条件下形成“快速通道”,诱发局部饱和现象;同时,基质吸力对水分再分布起主导作用,基质域与裂隙域在水分获取与传输上的竞争与耦合作用,对后期降雨入渗时空分布及孔压演化具有决定性影响。

关键词: 双重渗透介质, 降雨入渗, 非饱和渗流, 基质域, 裂隙域

Abstract:

To enhance landslide prediction and early warning capabilities and to optimize engineering disaster prevention and mitigation design, a non-saturated seepage model incorporating the coupling effect between matrix and fractures was developed based on dual-permeability theory. The model systematically accounted for rainfall intensity, rainfall duration, fracture permeability, spatial distribution and connectivity of fractures, matrix permeability, and the initial water content of the soil, enabling a realistic representation of water competition, migration, and coupling between the matrix and fracture domains. By using the improved rainfall boundary condition and considering the characteristics of unsaturated soil flow, COMSOL Multiphysics finite element software was employed as the computational platform, and the interactions between water in the matrix and fracture domains were comprehensively analyzed to reveal the rainfall infiltration characteristics of homogeneous soil. The results show that the presence of fractures significantly accelerates the initial infiltration rate and, under specific conditions, forms preferential flow paths, inducing local saturation; meanwhile, matrix suction plays a dominant role in water redistribution, and the competition and coupling between the matrix and fracture domains have a decisive effect on the spatiotemporal distribution of rainfall infiltration and the evolution of pore pressure.

Key words: dual-permeability medium, rainfall infiltration, unsaturated seepage, matrix domain, fracture domain

中图分类号: