ISSN 1003-8035 CN 11-2852/P
  • 中国科技核心期刊
  • CSCD收录期刊
  • Caj-cd规范获奖期刊
  • Scopus 收录期刊
  • DOAJ 收录期刊
  • GeoRef收录期刊
欢迎扫码关注“i环境微平台”

不同降雨条件下非饱和裂隙土优势流入渗特征研究

黄晓虎, 张芷宁, 张鹏, 魏兆亨, 陈建伟

黄晓虎,张芷宁,张鹏,等. 不同降雨条件下非饱和裂隙土优势流入渗特征研究[J]. 中国地质灾害与防治学报,2024,35(0): 1-12. DOI: 10.16031/j.cnki.issn.1003-8035.202403005
引用本文: 黄晓虎,张芷宁,张鹏,等. 不同降雨条件下非饱和裂隙土优势流入渗特征研究[J]. 中国地质灾害与防治学报,2024,35(0): 1-12. DOI: 10.16031/j.cnki.issn.1003-8035.202403005
HUANG Xiaohu,ZHANG Zhining,ZHANG Peng,et al. Study on the dominant infiltration characteristics of unsaturated fissure soil under different rainfall conditions[J]. The Chinese Journal of Geological Hazard and Control,2024,35(0): 1-12. DOI: 10.16031/j.cnki.issn.1003-8035.202403005
Citation: HUANG Xiaohu,ZHANG Zhining,ZHANG Peng,et al. Study on the dominant infiltration characteristics of unsaturated fissure soil under different rainfall conditions[J]. The Chinese Journal of Geological Hazard and Control,2024,35(0): 1-12. DOI: 10.16031/j.cnki.issn.1003-8035.202403005

不同降雨条件下非饱和裂隙土优势流入渗特征研究

基金项目: 国家自然科学基金青年基金项目(42007237)资助.
详细信息
    作者简介:

    黄晓虎(1987—),男,博士,副教授,硕士生导师,主要从事地质灾害预测与防治方面的研究。E-mail:88569096@qq.com

    通讯作者:

    张 鹏(1982—),男,博士,副教授,硕士生导师,主要从事地质灾害成灾机理与风险分析方面的研究。E-mail:pzhang@ctgu.edu.cn

Study on the dominant infiltration characteristics of unsaturated fissure soil under different rainfall conditions

  • 摘要:

    本文将膨胀土滑坡中的降雨入渗简化为垂直的斜坡土柱模型,利用Richards方程基本理论,基于离散孔隙-裂隙介质模型,结合短时强降雨和长时弱降雨两种典型降雨工况,以不同裂隙深度为主要变量,从入渗过程、饱和度剖面、入渗率、表面压力及入渗量来探讨裂隙在发育演化过程中的优势流入渗特征。结果表明,(1)裂隙处,降雨在裂隙末端形成暂饱和区,随后向裂隙两侧水平入渗,饱和方式为由下向上饱和;无裂隙处降雨主要沿土柱表面向下垂直入渗,为由上向下饱和。(2)短时强降雨作用下,随着裂隙深度增加,优势流控制着降雨入渗的过程,入渗量和入渗深度与裂隙深度呈正比,入渗深度与裂隙深度基本一致。(3)长时弱降雨作用下,优势流入渗效果减弱,优势流现象产生较慢,降雨沿裂隙入渗速度与在基质孔隙入渗速度基本一致,入渗主要表现为沿土柱表面向下的垂直入渗。

    Abstract:

    This paper simplifies the rainfall infiltration in expanding soil landslide into a vertical slope soil column model. Utilizing the basic theory of the Richards equation and based on a discrete pore-fracture medium model, combined with two typical rainfall conditions of short-term heavy rainfall and long-term light rainfall, with different fracture depths as the main variable, the study explores the preferential flow infiltration characteristics during the development and evolution of fractures from the dominant permeability, saturation profile, permeability rate, surface pressure, and permeability in the development of permeability characteristics. The results show that, (1) at the fissure, the rainfall forms a temporary saturation area at the end of the fissure, then horizontally infiltrates to both sides of the fissure, with the saturation mode being bottom-up. at non-fracture sites, rainfall mainly vertically infiltrates along the soil column surface, with the saturation mode being top-down.(2) Under short-term heavy rainfall, with the increase of fissure depth, the dominant flow controls the process of rainfall infiltration. The infiltration amount and infiltration depth are directly proportional to the fissure depth, and the infiltration depth is basically consistent with the fissure depth.(3) Under long-term weak rainfall, the infiltration effect of dominant inflow is weakened, and the phenomenon of dominant flow is slow. The infiltration speed of rainfall along the fissure is basically the same as that of the matrix pore, and the infiltration is mainly manifested as vertical infiltration downward along the surface of the soil column.

  • 由于西藏的特殊地理位置,为支撑西部社会经济发展、保障国防安全,近年来国家加大了对西藏的建设力度,其中标志性的重大建设工程即为川藏铁路。川藏铁路规划线路穿越青藏高原东缘高寒带岛状山地冻土区,除山峰为多年冻土外,大多段为季节性冻土区,冻融环境相对复杂[1],反复冻融循环势必对滑坡治理工程产生损伤,进而影响滑坡治理工程的有效性。

    目前滑坡治理工程冻融效应主要集中在挡墙及锚索的冻融效应。挡墙的冻融效应方面,国内的研究工作比较多,梁波等[2-3]对青藏铁路L型挡土墙的土压力进行了一个冻融循环的实测;汪恩良等[4] 通过模拟野外实测降温过程,进行加筋挡土墙在开敞条件下的冻融试验研究;赵坚等[5]根据我国近年来挡土墙抗冻技术研究成果,阐述了季节冻土区挡土墙抗冻结构设计的要求、荷载组合、冻胀力的设计取值及墙体强度和稳定的验算方法。季节性冻土冻融对锚索的影响方面,国内主要集中于对北方地区越冬基坑支护锚杆的监测及分析(ZHANG等[6]、赵延凯等[7]),认为冻季锚杆内力增加,暖季土体融化松散,容易出现裂缝坍塌,郭红仙等[8]在文献分析的基础上论述了土钉内力的冻融效应。国外研究边坡锚索工程冻融效应相对较早,如德国[9]、美国[10]、法国[11]均对边坡锚索工程进行越冬内力、位移监测,结果表明冻期锚杆内力可达到初始预应力的2~4倍,解冻后,锚杆内力将减少。

    上述滑坡治理工程冻融效应研究主要集中在青藏线及北方地区,极少涉及青藏高原东南部地区,针对川藏铁路规划线路季节性冻土地区的滑坡防治工程冻融效应分析及其对策的研究也未见报道。文章在对川西地区、昌都地区G317、G318沿线滑坡防治工程进行调查的基础上,分析总结季节性冻土区滑坡防治工程可能存在的问题及设计时所需注意的要点,为川藏铁路沿线滑坡治理工程及拟修建施工便道的边坡支护提供参考。

    川藏铁路规划线路穿越青藏高原东缘高寒带岛状山地冻土区,依次经过四川盆地、川西高山峡谷区、川西山地区、藏东南横断山高山峡谷区、藏南谷地区等5个地貌单元。从东向西依次跨越岷江、大渡河、雅砻江、金沙江、澜沧江、怒江、雅鲁藏布江等7 条大江大河,沿途地形跌宕起伏,岭谷高差可达5000 m。除山峰为多年冻土外,大多段为季节性冻土区,该区季节性冻土区下限约为海拔3000~3300 m,属于稳定型、南方型和亚热带型季节性冻土[12]

    徐斅祖、郭东信(1982年)编制了1∶400万中国冻土分布图;李树德、程国栋(1996年)编制了1∶300万青藏高原冻土图,为文中的研究提供了基础。川藏铁路沿线地区冻土分布、标准冻深等值线及本次案例点位置见图1,地形变化及重要站点季节性冻土厚度见图2。与滑坡防治工程直接相关的冻融环境因素除海拔、季节性冻土深度以外,还有冻结指数、冻结时间及混凝土年有害循环次数。以昌都地区为例,根据气象站点多年资料,冻期从10月开始,11月冻土深度为17 cm,12月冻土深度为62 cm,1—2月冻结深度达到最大,3—4月开始解冻,5月消失。昌都冻结指数(年内日平均温度中的负温累计绝对值)为232,为轻冻区;其混凝土年有害循环次数为131次[13],年有害循环次数相对较多,因此混凝土的冻害环境也相对复杂。

    图  1  研究区冻土分布图
    Figure  1.  Distribution of frozen soil in the study area
    图  2  川藏铁路沿线季节性冻土分布图
    Figure  2.  Seasonal frozen soil distribution along the Sichuan-Tibet railway

    本次选取川西、昌都季节性冻土区为研究对象,重点沿国道G317、G318 线进行调查。案例调查共涉及滑坡防治工程25处,共发现37处分项工程产生不同程度的冻融破坏。根据现场调查情况初步将冻融效应划分为四类情形:A、微弱:防治工程结构无明显冻伤或冻胀沉融变形;B、一般:防治工程结构完整,仅浅表层冻伤或产生微小裂隙及形变,不影响工程结构强度;C、中等:由于冻胀作用使得防治工程结构出现不可逆的较小结构性变形,损伤面较为局限,结构受力满足强度设计安全标准,且无恶化趋势,但进一步扩展将影响工程结构;D、严重:由于冻胀融沉循环作用使得防治工程结构出现不可逆的整体冻胀形变或是大面积严重直接冻害,并有加速劣化的趋势。

    对37处分项工程分析,滑坡防治工程冻融破坏方式及影响因素见图3,具体表现为:

    图  3  滑坡防治工程冻融破坏方式及影响因素
    Figure  3.  Freeze damage and influencing factors of landslide control projects

    ①挡墙冻融效应及破坏主要表现为:挡墙直接冻害,挡墙中下部出水口处最为严重,主要分布于海拔4000 ~4400 m;由于冻胀力远大于土压力使得挡墙产生贯通裂缝而失效,以竖向张拉裂缝为主。

    ②排水工程冻融效应及破坏主要表现为:高海拔(4200 ~4500 m)高寒、地下水埋深浅、细粒较多(黏土)的滑坡排水工程冻融效应明显。由于不均匀冻胀融沉,使得排水沟多沿横断面产生贯通裂缝,在地形转折处变形更为严重;由于冻土产生的侧胀力,使得排水沟侧壁垮塌。

    ③锚索冻融效应及破坏主要表现为:锚索内力因冻胀变形将增大,从而破坏锚头及格构;由于原始坡面或格构与挡墙相交处填土的自身沉降及融沉效应,或是融雪时对格构基础的破坏,使得锚索或锚杆的预应力损失;在海拔较高地区(3500 ~4400 m),存在格构直接冻害及冻土冻胀开裂的情形。

    ④锚喷工程冻融效应及破坏较为普遍,主要原因是锚喷层太薄、未加钢筋,且阻碍了坡体水的排泄,在坡面含水率较高、细粒物质较多的局部,容易冻裂,垮塌。

    ⑤桩板冻融效应及破坏冻融效应微弱,而桩间板的冻害与挡墙类似,且外挂板的受力模式更容易被破坏。损坏形式为受水平冻胀作用产生拉裂纹;混凝土致密性差,冻融作用破坏表层。主要原因:桩间土回填不密实,没有反滤层或土体饱水等,见图3(a)。

    ⑥标准冻深的影响:在冻胀率能相近的情况下,显然季节性冻土标准冻深越深,其冻胀融沉变形量越大。反复作用下,其对防治工程的劣化效应越明显。从图3(b)可见,产生冻融损伤的滑坡防治工程多分布于冻深大于60 cm的地区,冻深越深则防治工程冻融破坏程度越严重。

    ⑦岩土体性质的影响:川西地区、昌都地区G317、G318沿线除高原面上发育大面积黏土外,其余地区以粗粒碎石土为主,碎石土虽为弱冻胀土,但其含水量高、汇水面积大,地下水埋深较浅时,其对滑坡防治工程的冻融效应是不可忽视的,见图3(c)。

    ⑧高程的影响:从图3(d)可见,冻融破坏程度为“中等”及以上的案例多分布于海拔3700 ~4500 m。总体而言,海拔越高,气温越低,冻结指数越大,直接冻害越严重。

    调查中发现大量的混凝土结构产生严重的直接冻害,冻融循环作用对混凝土造成的不利影响主要体现在两个方面:表面剥落以及冻胀开裂,会出现混凝土质量损失和对相对动弹性模量产生影响,而对混凝土强度的影响则更大。

    混凝土冻害机理十分复杂,目前并没有公认的能解释全部冻害现象的理论。混凝土是由水泥、水和粗、细骨料组成的毛细孔多孔体,包括气孔、毛细孔、凝胶孔等。当环境温度降低到−1~−1.9℃时,混凝土孔隙中的水由大孔开始结冰,逐渐扩展到较小的孔,当温度在−l2℃时,毛细孔中的水开始结冰,而凝胶孔中的水一般不结冰。根据静水压假说[14-17],水转化为冰时体积膨胀9% ,迫使未结冰的孔溶液从结冰区向外迁移,因而产生静水压力,静水压力随孔隙水流程长度增加而增加,当孔隙水的流程长度大于某极限长度时,静水压力将超过混凝土的抗拉强度,从而造成破坏。

    季节性冻土地区滑坡防治工程应考虑采用抗冻水泥混凝土和抗冻水泥砂浆,其抗冻等级可参考《JTG/T D31-06—2017季节性冻土地区公路设计与施工技术规范》[13],根据其所处地区的有害冻融次数、饱水状态、有盐/无盐条件及设计年限来确定,抗冻等级可分为F100、F150、F200、F250、F300、F350、F400及F450。对于C60以下混凝土可考虑适量的引气措施来提高混凝土的抗冻性能。

    封闭系统冻胀是指冻胀过程中,没有外部水补给,仅由土体中原驻水引起的冻胀。挡墙墙后的土体冻胀属于双向冻结,当挡墙基础埋深足够时,垂直于地面的冻胀力只能抬高墙后土体,对墙无害,但垂直墙身的水平向冻胀力是直接造成挡土墙断裂、位移过大等主要因素。特别是翼墙拐角处,由于应力集中的原故,破坏更为严重。

    关于水平冻胀力的研究,前人已取得大量的成果,认为水平冻胀力远大于主动土压力,在特定条件下可达到10倍以上,《GBT 50662—2011 水工建筑物抗冰冻设计规范》[18]中对挡墙水平冻胀力大小及分布计算也作了相关规定。

    文中拟进行一个简单算例,以直观比较挡墙水平冻胀力与主动土压力大小。某挡土墙墙高5 m,墙背竖直,墙后填土为黏土,内摩擦角20°,黏聚力10 kPa,重度18 kN/m3,墙后填土冻胀级别为Ⅲ级,冻胀量12 cm。

    按照相关规范条例,其主动土压力与水平冻胀力分布见图4,其中Ht为填土高度,β 为与冻胀量相关的系数,Zd为挡墙基础埋深。多边形DEFG为水平冻胀力分布,其合力为169.3 kN/m;三角形ABC为主动土压力分布,其合力为51.4 kN/m。可见,该挡墙水平冻胀力大于主动土压力,数值上约为3.3倍。如果重力式挡墙进行抗滑、抗倾覆设计及墙身强度校核时未充分考虑水平冻胀力的影响,则挡墙势必造成整体稳定性问题或者墙身强度问题,挡墙可能产生倾斜、墙顶位移过大,或是墙身负弯矩过大,产生贯通的竖向裂纹等,影响挡土墙支挡功能的正常发挥。

    图  4  挡墙土压力及水平冻胀力分布图
    Figure  4.  Distribution of horizontal frost-heaving force

    开放系统冻胀是指冻胀过程中,除土体中有原驻水外,还有外部水源源不断地迁移、聚集至冻结面处,进而发生的更严重的冻胀。公路修建时,往往会破坏原有含水层,形成透水层、含水层及相对不透水层的三层体结构,挡墙排水孔多设置于排水层中,由于负温的影响,冬季泄水孔常被堵塞,透水层也开始冻结,在含水层顶部形成一个冰盖,而含水层中的水分不断向泄水孔周围迁移,使得冻胀相比较封闭系统而言显得更为剧烈,从而破坏挡墙工程(图5)。

    图  5  挡墙开敞式冻胀示意图
    Figure  5.  Freezing damage to retaining walls

    调查中发现,高海拔高寒、地下水埋深较浅、细粒较多的滑坡排水工程冻融效应较为明显,应当引起岩土工程师的重视。季节性冻土区排水沟冻胀破坏主要表现为侧壁垮塌及沿排水沟横剖面产生贯通裂缝,往往在地形转折处变形更为严重。一般需设置地表排水工程的滑坡其地面径流条件较好,土体常处于潮湿状态,含水量较高,因此其冻胀变形较大,故排水沟侧壁直接受水平冻胀力作用也较大,当冻胀荷载大于排水沟侧壁强度时,容易引起垮塌(图6)。

    图  6  排水沟冻胀破坏示意图
    Figure  6.  Freezing damage of drainage ditch

    排水沟基底主要受法向冻胀力,其方向垂直于等温面的,也近似垂直于排水沟基底。由于滑坡物质空间分布、含水率分布等不均,因此可能产生差异较大的法向冻胀力,使得排水沟整体上抬,造成排水沟内部产生内力集中,其基础底板可能产生较大的负弯矩,当超过其极限值时便开始开裂。在春季冻土融沉时,由于排水沟设置部位的地形突变、地基土性质及含水率的差异等,使得排水沟基础底板靠自重不能与融化的地基土同步下沉、复位,进一步扩展裂缝,造成下渗流量大,冲刷地基,在反复冻胀融沉作用下,排水沟可能产生破坏,进而影响滑坡的整体稳定性。

    热融滑塌是季节性冻土区边坡失稳的常见形式。季节性冻土区水文地质条件不良的粉质或黏性土滑坡,秋末多雨,结冻前土体原始含水量大, 土层冻结具备开放系统条件,地下水位接近或位于冻层,汇水条件好,大量地表水向边坡渗透聚集,使得冻融层湿度较高。春季时冻融层土体开始解冻,融化后土体的抗剪强度参数会有大幅度的折减。根据《JTG/T D31-06-2017季节性冻土地区公路设计与施工技术规范》[13],抗融滑稳定性计算公式为:

    Fs=Guw+qL)cosαtanφ+cL(G+qL)sinα (1)

    式中:Fs——抗融滑稳定性系数;

    G——土体自重;

    uw——水的扬压力;

    q——坡面砌体荷载;

    L——滑面长;

    α——滑面倾角;

    cφ——折减后的土体抗剪强度系数。

    黏性土的黏聚力折减系数为0.3~0.6,内摩擦角折减系数为0.15~0.16,同时受未及时排出的溶解水的影响,易在冻融层附近产生平直滑面型的滑移破坏,融滑深度为冻结线深度的1/3~2/3,多为0.1~0.5 m。

    由于排水沟的设计并未考虑支挡功能,热融滑塌的产生会对排水沟产生严重的影响。由于热融滑塌,会使得设置于相对稳定坡体上部的排水沟基础脱空,逐渐产生破坏,设置于坡脚处的排水沟将受到浅层融滑的剩余下滑力,因而产生挤压变形(图7图8) 。

    图  7  热融滑塌对排水沟的破坏照片
    Figure  7.  Damage caused by hot melt slump
    图  8  热融滑塌对排水沟的破坏示意图
    Figure  8.  Destruction of drainage ditch by hot melt slump

    以左贡县田妥镇斜马通村G318 公路K3663+750处滑坡为例[19],该滑坡位于西藏昌都地区左贡县田妥镇,海拔高度4100 m。冬春季气候干燥寒冷。年平均气温为6.2℃, 1月份平均气温为−6.7℃,年最大冻土深度 1 m。滑坡及周围坡体地层岩性主要为第四系全新统崩坡积物、第四系全新统冲洪积砂、卵石,以及下伏的三叠系上统巴贡群砂岩夹页岩及泥灰岩。该滑坡防治工程对锚索进行了为期2年的长期监测,数据见表1。可见由于滑坡变形或冻胀影响,最大锚固力相对于锁定锚固力分别增长了25.52%、41.38%、46.03 %及37.89%,最大者增长了300 kN,则锚索(杆)很可能被拉断或者锚头破坏,滑坡有失稳的可能。而当温度变暖时,由于土体融沉作用,锚索预应力相对于最大锚固力均有较大程度的降低,最低者降低320 kN,小于初始锚固力。

    表  1  锚索监测数据[19]
    Table  1.  Monitoring data of anchor cables[19]
    编号锁定锚固力
    /kN
    最大锚固力
    /kN
    最大增长率
    /%
    稳定值
    /kN
    339#48260525.52600
    345#63890241.38620
    357#63092046.03600
    358#805111037.89950
    下载: 导出CSV 
    | 显示表格

    根据典型案例调查,季节性冻土区锚索工程的冻融效应主要体现在以下方面: ①土体冻胀时,会产生指向坡外的变形,而锚索框格约束这种变形,势必在框格基底产生较大的法向应力,关于法向应力大小,规范中认为在基础不变形的情况下,最大可达到210 kPa。法向应力增强使得框格梁负弯矩过大,产生贯通裂缝。②锚索约束并不能完全消除冻胀变形,锚索也会产生指向坡外的变形,因而使得内力增大,由于锚索的弹模较大,即使产生很小的位移,也会引起内力大幅度增加,主要表现为破坏锚索外锚头。③由于土体融缩变形,锚索也随之回缩,同时地基土融化时,其地基承载力也相应降低,因此锚索预应力将产生较大的预应力损失,甚至存在锚具脱空预应力完全丧失的情况。

    为减缓季节性冻土区滑坡防治工程冻融效应,可采取以下措施:

    (1)框格锚索工程可采用耗能结构及应力保持装置,对框格梁适当提高标号,并采取保温措施。

    (2)排水沟减少地表裸露的方式,改用渗水盲沟、保温沟、冷暖沟的形式;排水沟基础采取保温措施,地形转折处应加强截面强度;泉水出露处采用集水井。

    (3)抗滑挡墙及桩间板适当提高混凝土标号应做好疏干、排水工作。桩间板应采取内挂的方式。

    (4)锚喷工程应与主动网联合使用。

    针对季节性冻土针对冻融条件下岩土体冻胀、热融将造成锚固工程预应力陡增、损失的问题,除控制填土细粒含量及对锚索(杆)及框格梁采取一定的安全系数外,项目组提出一种可以抵抗冻土冻胀、融缩变形的锚索结构(图9),最小化其预应力损失,使其满足冻融环境下的深层滑坡的防治要求:

    (1)利用屈服锚具调节较大的由冻胀变形或滑坡滑动变形预应力增加,可兼顾抗震设计。

    (2)利用碟簧蓄能系统,抵消部分预应力损失,也可抵抗较小的由变形引起的预应力增加。

    (3)行程指标:50 mm。吨位要求:50 t,4根钢绞线。

    图  9  双向调节预应力锚索结构图
    1—碟簧耗能系统;2—注浆管;3—保护套筒;4—退位卡套;5—安全卡套;6—封底板;7—钢绞线; 8—锚具;9—架线环;10—导向帽;11—承力转换板;12—自由段; 13—锚固段
    Figure  9.  Structure of bi-directional adjusting prestressed anchor cable

    50 t级双向自动调节式预应力锚索额定工作锚固力为400 kN,在冻结期冻胀作用下,当锚固力小于设计屈服力500 kN时,应力保持装置发挥作用,大大降低由于冻胀产生的预应力增量,当锚固力大于设计屈服力500 kN时,屈服装置发挥作用,可维持500 kN锚固力;在融化期融沉作用下,应力保持装置发挥作用,能补偿大部分由于融沉作用产生的预应力损失;在地震或者大变形作用下,屈服装置发挥作用。

    (1)川藏铁路规划线路穿越青藏高原东缘高寒带岛状山地冻土区,除山峰为多年冻土外,大多段为季节性冻土区,冻融环境相对复杂,滑坡治理中应关注冻融对工程的影响。

    (2)通过对已建工程G317、G318滑坡及边坡治理工程运行现状的实地调研,总结了冻融条件下各常用工程的破坏形式及受力机理,认为开放式冻胀较封闭式冻胀对防治工程的影响更大;并指出季节性冻土区滑坡防治应注意的水平冻胀力、材料抗冻等级等问题。

    (3)根据冻融条件下滑坡防治工程破坏形式,针对性提出对锚索采用耗能结构,对锚喷工程提出锚喷与主动网联合使用,对排水工程应当采取保温措施且在地形转折处应加强等。

    (4)文中以案例调查为基础来探讨滑坡防治工程的冻融效应,以期引起岩土工程师对这一问题的关注,但冻融效应涉及到冻胀理论、土体冻胀性能、与结构协调变形以及防治工程劣化机理等,尚需更多实测数据来支撑。

  • 图  1   优势流入渗的无限边坡模型

    Figure  1.   Infinite slope model of dominant infiltration and seepage

    图  2   数值模拟模型

    Figure  2.   Numerical simulation model

    图  3   网格划分

    Figure  3.   Illustration of Mesh division

    图  4   研究区降雨情况

    Figure  4.   Rainfall situation in the study area

    图  5   短时强降雨下压力水头分布随降雨时间的变化过程

    Figure  5.   Variation process of pressure head distribution with rainfall time under short-term heavy rainfall

    图  6   长时弱降雨下压力水头分布随降雨时间的变化过程

    Figure  6.   Variation process of pressure head distribution with rainfall time under long-time weak rainfall

    图  7   短时强降雨下饱和度剖面随降雨时间的变化过程

    Figure  7.   Variation process of saturation profile with rainfall time under short-term heavy rainfall

    图  8   长时弱降雨下饱和度剖面随降雨时间的变化过程

    Figure  8.   Variation process of saturation profile with rainfall time under long-time weak rainfall

    图  9   短时强降雨下入渗率随降雨时间的变化曲线

    Figure  9.   Variation curve of infiltration rate with rainfall time under short-term heavy rainfall

    图  10   长时弱降雨下入渗率随降雨时间的变化曲线

    Figure  10.   Variation curve of infiltration rate with rainfall time under long-term weak rainfall

    图  11   短时强降雨下入渗量及表面压力变化曲线

    Figure  11.   Variation curve of infiltration volume and surface pressure under short-term heavy rainfall

    图  12   长时弱降雨下入渗量及表面压力变化曲线

    Figure  12.   Variation curve of infiltration volume and surface pressure under long-term weak rainfall

    图  13   不同工况下降雨入渗深度变化曲线

    Figure  13.   Variation curve of rainfall infiltration depth under different working conditions

    图  14   不同工况下总入渗量变化曲线

    Figure  14.   Variation curve of total infiltration volume under different working conditions

    图  15   短时强降雨作用下优势流入渗过程示意图

    Figure  15.   Schematic diagram of dominant inflow and infiltration process under short-term heavy rainfall

    图  16   长时弱降雨作用下优势流入渗过程示意图

    Figure  16.   Schematic diagram of dominant inflow and infiltration process under long-term weak rainfall

    表  1   模拟工况

    Table  1   Simulated working conditions

    工况 编号 裂隙深度(m) 裂隙宽度(mm) kf(m/s) 降雨
    短时
    强降雨
    1 无裂隙 - - 120 mm/d,
    持续48 h
    2 0.5 5 0.032
    3 1 10 0.131
    4 2 20 0.527
    长时
    弱降雨
    5 无裂隙 - - 12 mm/d,
    持续480 h
    6 0.5 5 0.032
    7 1 10 0.131
    8 2 20 0.527
    下载: 导出CSV

    表  2   模拟参数

    Table  2   Simulation parameters

    符号 名称 数值 来源
    ks 饱和渗透系数 5.787×10−8 m/s 勘查资料
    α V-G模型本构参数 0.272 经验取值
    n V-G模型本构参数 2.3 经验取值
    θs 饱和含水率 32.6% 室内试验
    θr 残余含水率 6.4% V-G模型拟合
    α B-C模型本构参数 6.8 经验取值
    n B-C模型本构参数 1 经验取值
    下载: 导出CSV
  • [1]

    FOX G A,WILSON G V. The role of subsurface flow in hillslope and stream bank erosion:a review[J]. Soil Science Society of America Journal,2010,74(3):717 − 733. DOI: 10.2136/sssaj2009.0319

    [2]

    CHEN H X,ZHANG L M. A physically-based distributed cell model for predicting regional rainfall-induced shallow slope failures[J]. Engineering Geology,2014,176:79 − 92. DOI: 10.1016/j.enggeo.2014.04.011

    [3] 邓威,肖世国. 含裂隙近水平红层软岩边坡渗透稳定性模型试验[J]. 水文地质工程地质,2024,51(1):57 − 68. [DENG Wei,XIAO Shiguo. Model test on stability of soft rock slopes composed of nearly horizontal redbeds with cracks[J]. Hydrogeology & Engineering Geology,2024,51(1):57 − 68. (in Chinese with English abstract)]

    DENG Wei, XIAO Shiguo. Model test on stability of soft rock slopes composed of nearly horizontal redbeds with cracks[J]. Hydrogeology & Engineering Geology, 2024, 51(1): 57 − 68. (in Chinese with English abstract)

    [4] 吴玮江,宋丙辉,刘迪,等. 黄土塬区包气带水分运移特征研究[J]. 水文地质工程地质,2023,50(3):12 − 22. [WU Weijiang,SONG Binghui,LIU Di,et al. Research on the characteristics of water transport in the aeration zone of loess tableland[J]. Hydrogeology & Engineering Geology,2023,50(3):12 − 22. (in Chinese with English abstract)]

    WU Weijiang, SONG Binghui, LIU Di, et al. Research on the characteristics of water transport in the aeration zone of loess tableland[J]. Hydrogeology & Engineering Geology, 2023, 50(3): 12 − 22. (in Chinese with English abstract)

    [5]

    VAN ASCH T W J,HENDRIKS M R,HESSEL R,et al. Hydrological triggering conditions of landslides in varved clays in the French Alps[J]. Engineering Geology,1996,42(4):239 − 251. DOI: 10.1016/0013-7952(95)00082-8

    [6]

    HENCHER S R. Preferential flow paths through soil and rock and their association with landslides[J]. Hydrological Processes,2010,24(12):1610 − 1630. DOI: 10.1002/hyp.7721

    [7] 杨宗佶,蔡焕,雷小芹,等. 非饱和地震滑坡堆积体降雨破坏水-力耦合行为试验[J]. 岩土力学,2019,40(5):1869 − 1880. [YANG Zongji,CAI Huan,LEI Xiaoqin,et al. Experiment on hydro-mechanical behavior of unsaturated gravelly soil slope[J]. Rock and Soil Mechanics,2019,40(5):1869 − 1880. (in Chinese with English abstract)]

    YANG Zongji, CAI Huan, LEI Xiaoqin, et al. Experiment on hydro-mechanical behavior of unsaturated gravelly soil slope[J]. Rock and Soil Mechanics, 2019, 40(5): 1869 − 1880. (in Chinese with English abstract)

    [8]

    KRZEMINSKA D M,BOGAARD T A,MALET J P,et al. A model of hydrological and mechanical feedbacks of preferential fissure flow in a slow-moving landslide[J]. Hydrology and Earth System Sciences,2013,17(3):947 − 959. DOI: 10.5194/hess-17-947-2013

    [9] 刘向峰,于冰,郝国亮,等. 降雨作用下多软弱夹层边坡滑移机理分析[J]. 中国地质灾害与防治学报,2023,34(0):1 − 12. [LIU Xiangfeng,YU Bing,HAO Guoliang,et al. Analysis of sliding mechanism of the cut slopes with multi-weak interlayers under rainfall[J]. The Chinese Journal of Geological Hazard and Control,2023,34(0):1 − 12. (in Chinese with English abstract)]

    LIU Xiangfeng, YU Bing, HAO Guoliang, et al. Analysis of sliding mechanism of the cut slopes with multi-weak interlayers under rainfall[J]. The Chinese Journal of Geological Hazard and Control, 2023, 34(0): 1 − 12. (in Chinese with English abstract)

    [10] 李斌,郝继锋,鞠远江,等. 渠坡非饱和膨胀土含水率与强度关系试验研究[J]. 水文地质工程地质,2022,49(5):129 − 136. [LI Bin,HAO Jifeng,JU Yuanjiang,et al. An experimental study of the relationship between water content and strength of unsaturated expansive soil on canal slope[J]. Hydrogeology & Engineering Geology,2022,49(5):129 − 136. (in Chinese with English abstract)]

    LI Bin, HAO Jifeng, JU Yuanjiang, et al. An experimental study of the relationship between water content and strength of unsaturated expansive soil on canal slope[J]. Hydrogeology & Engineering Geology, 2022, 49(5): 129 − 136. (in Chinese with English abstract)

    [11] 陈亮胜,韦秉旭,廖欢,等. 膨胀土边坡非饱和渗流及渐进性破坏耦合分析[J]. 水文地质工程地质,2020,47(4):132 − 140. [CHEN Liangsheng,WEI Bingxu,LIAO Huan,et al. A coupling analysis of unsaturated seepage and progressive failure of an expansive soil slope[J]. Hydrogeology & Engineering Geology,2020,47(4):132 − 140. (in Chinese with English abstract)]

    CHEN Liangsheng, WEI Bingxu, LIAO Huan, et al. A coupling analysis of unsaturated seepage and progressive failure of an expansive soil slope[J]. Hydrogeology & Engineering Geology, 2020, 47(4): 132 − 140. (in Chinese with English abstract)

    [12] 曹玲. 膨胀土的裂隙演化机制及其边坡稳定性分析[D]. 上海:上海交通大学,2016. [CAO Ling. Evolution law of Cracks&Slope stability analysis of expansive soil[D]. Shanghai:Shanghai Jiao Tong University,2016. (in Chinese with English abstract)]

    CAO Ling. Evolution law of Cracks&Slope stability analysis of expansive soil[D]. Shanghai: Shanghai Jiao Tong University, 2016. (in Chinese with English abstract)

    [13]

    PEI Pei,ZHAO Yanlin,NI Pengpeng,et al. A protective measure for expansive soil slopes based on moisture content control[J]. Engineering Geology,2020,269:105527. DOI: 10.1016/j.enggeo.2020.105527

    [14]

    INDRAWAN I G B,RAHARDJO H,LEONG E C. Drying and wetting characteristics of a two-layer soil column[J]. Canadian Geotechnical Journal,2007,44(1):20 − 32. DOI: 10.1139/t06-090

    [15] 朱伟,陈学东,钟小春. 降雨入渗规律的实测与分析[J]. 岩土力学,2006,27(11):1873 − 1879. [ZHU Wei,CHEN Xuedong,ZHONG Xiaochun. Observation and analysis of rainfall infiltration[J]. Rock and Soil Mechanics,2006,27(11):1873 − 1879. (in Chinese with English abstract)] DOI: 10.3969/j.issn.1000-7598.2006.11.003

    ZHU Wei, CHEN Xuedong, ZHONG Xiaochun. Observation and analysis of rainfall infiltration[J]. Rock and Soil Mechanics, 2006, 27(11): 1873 − 1879. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-7598.2006.11.003

    [16] 马佳. 裂土优势流与边坡稳定性分析方法[D]. 武汉:中国科学院研究生院(武汉岩土力学研究所),2007. [MA Jia. Preferential flow and stability analysis method for fissure clay slopes[D]. Wuhan:Institute of Rock and Soil Mechanics,Chinese Academy of Sciences,2007. (in Chinese with English abstract)]

    MA Jia. Preferential flow and stability analysis method for fissure clay slopes[D]. Wuhan: Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, 2007. (in Chinese with English abstract)

    [17]

    BROOKS R H,COREY A T. Properties of porous media affecting fluid flow[J]. Journal of the Irrigation and Drainage Division,1966,92(2):61 − 88. DOI: 10.1061/JRCEA4.0000425

    [18] 周创兵,熊文林. 岩石节理的渗流广义立方定理[J]. 岩土力学,1996,17(4):1 − 7. [ZHOU Chuangbing,XIONG Wenlin. Generalized cubic theorem of seepage in rock joints[J]. Rock and Soil Mechanics,1996,17(4):1 − 7. (in Chinese)]

    ZHOU Chuangbing, XIONG Wenlin. Generalized cubic theorem of seepage in rock joints[J]. Rock and Soil Mechanics, 1996, 17(4): 1 − 7. (in Chinese)

    [19]

    XIAO Weimin,XIA Caichu,WEI Wang,et al. Combined effect of tortuosity and surface roughness on estimation of flow rate through a single rough joint[J]. Journal of Geophysics and Engineering,2013,10(4):045015. DOI: 10.1088/1742-2132/10/4/045015

    [20]

    BOUDREAU B P. The diffusive tortuosity of fine-grained unlithified sediments[J]. Geochimica et Cosmochimica Acta,1996,60(16):3139 − 3142. DOI: 10.1016/0016-7037(96)00158-5

    [21]

    BELEM T,HOMAND-ETIENNE F,SOULEY M. Quantitative parameters for rock joint surface roughness[J]. Rock Mechanics and Rock Engineering,2000,33(4):217 − 242. DOI: 10.1007/s006030070001

    [22] 赵宽耀,许强,刘方洲,等. 黄土中优势通道渗流特征研究[J]. 岩土工程学报,2020,42(5):941 − 950. [ZHAO Kuanyao,XU Qiang,LIU Fangzhou,et al. Seepage characteristics of preferential flow in loess[J]. Chinese Journal of Geotechnical Engineering,2020,42(5):941 − 950. (in Chinese with English abstract)] DOI: 10.11779/CJGE202005017

    ZHAO Kuanyao, XU Qiang, LIU Fangzhou, et al. Seepage characteristics of preferential flow in loess[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(5): 941 − 950. (in Chinese with English abstract) DOI: 10.11779/CJGE202005017

    [23] 年庚乾,陈忠辉,周子涵,等. 基于双重介质模型的裂隙岩质边坡渗流及稳定性分析[J]. 煤炭学报,2020,45(增刊2):736 − 746. [NIAN Gengqian,CHEN Zhonghui,ZHOU Zihan,et al. Seepage and stability analysis of fractured rock slope based on dual medium model[J]. Journal of China Coal Society,2020,45(Sup 2):736 − 746. (in Chinese)]

    NIAN Gengqian, CHEN Zhonghui, ZHOU Zihan, et al. Seepage and stability analysis of fractured rock slope based on dual medium model[J]. Journal of China Coal Society, 2020, 45(Sup 2): 736 − 746. (in Chinese)

    [24] 侯晓萍,樊恒辉. 基于COMSOL Multiphysics的非饱和裂隙土降雨入渗特性研究[J]. 岩土力学,2022,43(2):563 − 572. [HOU Xiaoping,FAN Henghui. Study on rainfall infiltration characteristics of unsaturated fractured soil based on COMSOL Multiphysics[J]. Rock and Soil Mechanics,2022,43(2):563 − 572. (in Chinese with English abstract)]

    HOU Xiaoping, FAN Henghui. Study on rainfall infiltration characteristics of unsaturated fractured soil based on COMSOL Multiphysics[J]. Rock and Soil Mechanics, 2022, 43(2): 563 − 572. (in Chinese with English abstract)

    [25] 姚云琦,曾润强,马建花,等. 考虑优势流作用的降雨入渗边坡可靠度分析[J]. 岩土力学,2022,43(8):2305 − 2316. [YAO Yunqi,ZENG Runqiang,MA Jianhua,et al. Reliability analysis of slope under rainfall infiltration considering preferential flow model[J]. Rock and Soil Mechanics,2022,43(8):2305 − 2316. (in Chinese with English abstract)]

    YAO Yunqi, ZENG Runqiang, MA Jianhua, et al. Reliability analysis of slope under rainfall infiltration considering preferential flow model[J]. Rock and Soil Mechanics, 2022, 43(8): 2305 − 2316. (in Chinese with English abstract)

    [26] 胡静,陈胜宏. 渗流分析中排水孔模拟的空气单元法[J]. 岩土力学,2003,24(2):281 − 283. [HU Jing,CHEN Shenghong. Air element methad for modeling drainage holes in seepage analysis[J]. Rock and Soil Mechanics,2003,24(2):281 − 283. (in Chinese with English abstract)] DOI: 10.3969/j.issn.1000-7598.2003.02.029

    HU Jing, CHEN Shenghong. Air element methad for modeling drainage holes in seepage analysis[J]. Rock and Soil Mechanics, 2003, 24(2): 281 − 283. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-7598.2003.02.029

    [27] 侯晓萍,徐青,陈胜宏. 用空气单元法求解渗流场的逸出边界问题[J]. 岩土力学,2015,36(8):2345 − 2351. [HOU Xiaoping,XU Qing,CHEN Shenghong. An air element method for overflow boundary problem in seepage analysis[J]. Rock and Soil Mechanics,2015,36(8):2345 − 2351. (in Chinese with English abstract)]

    HOU Xiaoping, XU Qing, CHEN Shenghong. An air element method for overflow boundary problem in seepage analysis[J]. Rock and Soil Mechanics, 2015, 36(8): 2345 − 2351. (in Chinese with English abstract)

    [28]

    LEE Lee min,GOFAR N,RAHARDJO H. A simple model for preliminary evaluation of rainfall-induced slope instability[J]. Engineering Geology,2009,108(3/4):272 − 285.

    [29]

    SILLERS W S,FREDLUND D G. Statistical assessment of soil-water characteristic curve models for geotechnical engineering[J]. Canadian Geotechnical Journal,2001,38(6):1297 − 1313. DOI: 10.1139/t01-066

    [30] 李林翠. 压实黄土持水和降雨入渗特性试验研究[D]. 西安:长安大学,2021. [LI Lincui. An experimental study on soil water retension characteristics and rainfall infiltration characteristics of the compacted soils[D]. Xi’an:Changan University,2021. (in Chinese with English abstract)]

    LI Lincui. An experimental study on soil water retension characteristics and rainfall infiltration characteristics of the compacted soils[D]. Xi’an: Changan University, 2021. (in Chinese with English abstract)

    [31] 黄晓虎,魏兆亨,易武等. 裂隙优势流入渗诱发堆积层滑坡浅层破坏机理研究[J/OL]. 岩土工程学报. [Huang Xiaohu,Wei Zhaoheng,Yi Wu. et al. Study on the mechanism of fracture preferential flow infiltration-induced shallow damage of colluvial landslides,Chinese Journal of Geotechnical Engineering. (in Chinese with English abstract)]

    Huang Xiaohu, Wei Zhaoheng, Yi Wu. et al. Study on the mechanism of fracture preferential flow infiltration-induced shallow damage of colluvial landslides, Chinese Journal of Geotechnical Engineering. (in Chinese with English abstract)

    [32]

    MEIN R G,LARSON C L. Modeling infiltration during a steady rain[J]. Water Resources Research,1973,9(2):384 − 394. DOI: 10.1029/WR009i002p00384

  • 期刊类型引用(3)

    1. 李聪然,马明星,潘懿,宋军涛,杜晓飞,刘文进,岳祖润. 季冻区露天矿山边坡稳定性分析及加固措施研究. 河北地质大学学报. 2024(03): 73-79 . 百度学术
    2. 王谊,吴浩秋,彭善平,全永琴,蒋峻峰. 高原冻土滑坡治理措施技术分析. 价值工程. 2024(29): 4-7 . 百度学术
    3. 黄德昕,温韬,贾文君,胡明毅. 基于CiteSpace的国内岩石力学研究热点及趋势分析. 资源环境与工程. 2023(06): 792-800 . 百度学术

    其他类型引用(4)

图(16)  /  表(2)
计量
  • 文章访问数:  141
  • HTML全文浏览量:  83
  • PDF下载量:  0
  • 被引次数: 7
出版历程
  • 收稿日期:  2024-03-01
  • 修回日期:  2024-03-12
  • 录用日期:  2024-05-23
  • 网络出版日期:  2024-07-02

目录

/

返回文章
返回