降雨诱发边坡破坏的数值模拟.pdf

第 29 卷第 11 期2023 年 11 月水利科技与经济Water Conservancy Science and Technology and EconomyVol.29 No.11November,2023收稿日期 2023-05-06作者简介 袁文辉(1968-),男,江西南昌人,工程师,研究方向为堤防、水闸管理.doi:10.3969/j.issn.1006-7175.2023.11.003降雨诱发边坡破坏的数值模拟袁文辉(南昌市城市防洪事务中心,南昌 330006)摘 要土壤水分特征、强度和应力应变特性等岩土工程性质以及气候条件,已成为边坡稳定性数值模拟的重要参数。为了探究降雨强度和降雨时间对边坡稳定性影响,通过考虑水力和强度参数进行数值模拟,采用有限差分法,评估降雨期间非饱和边坡的稳定性,并通过参数研究,分析土壤渗透性、SWCC 和黏聚力对边坡稳定性的影响。结果表明,不同的降雨强度会产生不同的边坡破坏模式;渗透性和 SWCC 特性在非饱和土壤的渗透过程中产生超孔隙水压力;黏聚力对非饱和边坡的安全系数有显著影响。关键词边坡稳定性;降雨强度;数值模拟中图分类号 TU432 文献标识码 A 文章编号 1006-7175(2023)11-0011-041 概 述暴雨引起的边坡破坏是重要的岩土工程问题,负孔隙水压力(基质吸力)在暴雨诱发的边坡失稳中起着重要作用。其中,导致边坡破坏的机制是雨水开始渗入非饱和土壤时,引起负孔隙水压力的损失。为此,可以采取非饱和土壤的边坡稳定性分析。同时,为了更好地理解土壤中水的密度与水相能量状态之间的关系,还需要绘制土壤-水特征曲线(SWCC)。模拟非饱和土壤中的水分入渗,除了 SWCC 参数外,降雨情景(即降雨强度和持续时间)、土壤渗透性和土壤中的初始条件也是非饱和边坡中降雨入渗数值模拟所需的参数。目前,许多学者针对这些参数的影响进行了相关研究。王花蕊1利用灰色系统理论,探究了降雨强度、内摩擦角、内聚力对边坡稳定性的影响。陈乔等2利用降雨型滑坡监测模型,对降雨强度、孔隙水压力、含水率对边坡坡失稳破坏成因进行了分析。刘平军3分析了边坡稳定性系数与降雨强度的关系,并推导了边坡稳定性公式。刘亚龙4通过研究发现,浅层边坡滑动主要出现在强降雨和暴雨中后期,且边坡稳定性与降雨方式、强度、土体入渗能力有关。本文利用模型的几何、水力和岩土性质作为数值模拟的输入参数,通过参数研究土壤渗透性、SWCC 和黏聚力对边坡稳定性的影响。结果表明,不同的降雨强度会产生不同的边坡破坏模式,渗透性和 SWCC 在非饱和土壤的地面入渗行为中起着重要作用,土壤黏聚力也是非饱和边坡安全的一个重要控制因素。2 数值模拟分析2.1 构建模型选取水库上游典型边坡截面,目标斜坡高度30m,坡度 3.3 1。研究边坡的代表性截面见图 1。图 1 研究边坡的代表性截面11第 29 卷第 11 期2023 年 11 月水利科技与经济Water Conservancy Science and Technology and EconomyVol.29 No.11November,2023研究参数为渗透性、降雨强度和持续时间、土壤饱和度和黏聚力。为了更准确地预测降雨期间的浅层边坡破坏,应考虑选择内聚值作为非饱和边坡的初始属性。根据统一的土壤分类系统,该区域的风化花岗岩土壤可归类为 SM,沿坡面分布,深度在210m 之间,假定该边坡为风化花岗岩土层的均质边坡。该区域土壤的物理力学性质和 SWCC参数见表 1。表 1 土壤物理力学性质和土壤水分性质项目物理力学性质参数值土壤干容重/tm-31.48Eff 凝聚力/kPa0Eff 摩擦角/()28.30杨氏模量/Pa8.4E7位置比0.32初始饱和度0.5饱和渗透率/cmh-11.08有空性0.1土壤-水a/m-11.91性质n1.503 本研究包括渗透分析和边坡稳定性分析,可与 FLAC2D-两相流配置同时进行模拟。Van Ge-nuchten 方程式(1980)用于研究不同降雨条件下雨水对边坡土壤性质的影响;非饱和土的抗剪强度折减法用于估算雨水渗入土壤的量对边坡安全系数(FOS)的影响。2.2 水流过非饱和土壤在 FLAC 中,Van Genuchten(1980)基于毛细管模型建立一组封闭式方程,用于表示非饱和土壤的水力特性。为了实现由于不同的坡度角,雨水沿坡面渗入土壤的量的可变性,编辑 FLAC 中提供的调整流量 FISH。非饱和土壤的水力特性公式如下:Se=Sw-Sr1-Sr11+()nm(1)K()=KsKr=KsS0.5e1-(1-S1/me)m2(2)式中:Se为有效饱和度;Sw为饱和度;Sr为剩余饱和度;n、m 为曲线拟合系数,m=1-1/n,n1;Ks、Kr为饱和相对渗透率。2.3 非饱和土的抗剪强度折减方法使用 Bishop 有效应力(1959),稳定性分析中采用的非饱和土的抗剪强度方程如下:max=(-ua)tan+(ua-uw)tan+c(3)或者弗雷德隆德和拉哈德乔(1993 年)所示:max=(-ua)tan+(ua-uw)tanb+c(4)这里 Vanapalli 等(1996 年)将定义为:=(Sw-Sr)/(100-Sr)(5)式中:Sw、Sr、ua、uw、分别为饱和度、残余饱和度、气压、水压和总应力。在利用抗剪强度折减技术计算非饱和边坡的安全系数时,采用折减后的抗剪强度参数 Ctrial和 trial进行一系列模拟,试验定义如下:Ctrial=1FtrialC(6)trial=arctan1Ftrialtan()(7)式中:Ftrial为试验的安全系数,开始剪切强度降低 Ftrial设置为足够小,以确保系统稳定;Ftrial逐渐增加,直到发生破坏,FOS 是影响其破坏的关键因素之一。2.4 可变参数分析中的可变参数包括降雨强度分布、饱和度、渗透系数、SWCC 参数和土壤黏聚力。2.4.1 降雨模式本研究假设不同的降雨强度。为了估算土壤的渗透能力,使用与土壤渗透性相等的降雨强度,并选择与流域当前最高降雨强度相似的降雨强度 11cm/h。沿斜坡降雨 110h,预计 FOS 增量为 10h。2.4.2 饱和渗透系数使用两种不同的渗透率,研究其对渗透过程和正孔隙水压力随时间的发展影响。在使用的渗透率中,分别为 1.08 和 10.8cm/h。2.4.3 SWCC 参数SWCC 提供了土壤中水的质量/体积与水相能量状态之间的概念性理解。根据 SWCC 实验,、n 和 r被确定并用作 FLAC 中非饱和土的输入特性。n 值用于计算 Van Genuchten 方程中的参数,(m-1)代表参考孔隙压力 P0(kPa)。另一个拟合参数是剩余饱和度 r,与孔隙压力计算中使用的有效饱和度有关。将残余含水量定义为土壤水分从主要通过毛细作用保持在土壤中到主要通过吸附力保持在土壤中的含水量。在本研究中,为了估计该参数的影响,使用 FOS 获得的0 和 0.11 的值。其中,零剩余饱和度是一种特殊21袁文辉:降雨诱发边坡破坏的数值模拟第 11 期情况,因为饱和度等于有效饱和度。2.4.4 饱和度为了估计饱和度对边坡失稳阈值随时间的影响,假设不同的饱和度值。这些估计值分别为0.4、0.5 和 0.6,其中 0.5 为实验值。2.4.5 土壤黏聚力土壤黏聚力被理解为将土壤颗粒固定在一起的力。对于非饱和土,黏聚力是联锁、物理、化学作用产生的有效内聚力以及由颗粒间力或吸力和表面张力控制的表观内聚力的总和。为了通过数值模拟了解所选内聚力对降雨条件下边坡破坏的合理预测影响,假设 5 个值。假定的黏聚力值为 0、5、10、20 和 30kPa。3 结果和分析3.1 渗透性的影响 为了探究降雨 20h 后孔隙压力随深度的变化,两种不同的渗透性值 k1(k1=1.08cm/h)和k2(k2=10.8cm/h)导致土壤在降雨条件下的不同行为。在这种情况下,假设土壤的内聚力和饱和度分别为 10kPa 和 0.5kPa。随着渗透性的增加,更多的水渗入边坡。对于渗透值较高的土壤,正孔隙压力在边坡表面附近显著发展。降雨20h 后孔隙压力随深度的变化见图 2;雨水渗透对边坡安全系数的影响见图 3。图 2 降雨 20h 后孔隙压力变化图 3 雨水渗透对边坡安全系数的影响由图 2 和图 3 可知,由于渗透率较低,基质吸力显著降低。由图 2 和图 3 能够估计给定降雨时间后形成的湿润带,对于降雨为 k1 的边坡,可以估计约 6m 的湿润深度。对于低饱和渗透系数,边坡内的负孔隙压力下降得更慢,任何变化只能在地表附近观察到。因此,在渗透率较低的情况下,可以看到 FOS 的稳步下降。另外,在受降雨影响的高渗透性土壤中,FOS 显著下降。3.2 降雨强度的影响为了估计降雨强度对边坡 FOS 的影响,考虑低渗透性土壤 k1,该边坡受到不同降雨强度的影响:r1=1.08cm/h 和 r2=10.8cm/h。不同降雨强度下的安全系数见图 4。图 4 不同降雨强度下的安全系数从图 4 可以看出,降雨越大,FOS 下降越快,导致边坡失稳的速度越快。这种现象还可以用雨水渗透引起的负孔隙压力变化来解释,降雨量越大,渗入边坡的水量越大,进而导致孔隙压力显著增加。3.3 初始饱和度和剩余饱和度的影响土壤饱和度是影响边坡初始 FOS 值的主要因素。如果渗透系数和降雨强度影响降雨期间 FOS的下降,则饱和度直接决定降雨开始前边坡的稳定性阈值。对于不同的饱和度,可得到不同的自由空间。由图 5 和图 6 可知,降雨期间 FOS 的变化在低降雨强度或高降雨强度下都有类似的趋势。当考虑残余饱和度时,也可以看到类似的效果。图 5 降雨为 1.08cm/h 时的安全系数31第 29 卷第 11 期2023 年 11 月水利科技与经济Water Conservancy Science and Technology and EconomyVol.29 No.11November,2023图 6 降雨为 10.8cm/h 时的安全系数由图 7 可以看出,当考虑剩余饱和度时,可得到较小的有效饱和度。而剩余饱和度的存在表明,有效饱和度影响孔隙压力分布,进而影响获得的 FOS。有效饱和度越低,坡度越稳定。图 7 剩余饱和度对安全系数的影响3.4 土壤黏聚力的影响为研究黏性的影响,采用饱和度为 0.5、渗透性为 1.08cm/h、黏性变化较大的土壤,有效黏聚力对安全系数的影响见图 8。为了确认孔隙压力的变化是否在这种情况下起重要作用,分析降雨70h 后的孔隙压力和相应内聚值的边坡破坏面,见图 9。图 8 有效黏聚力对安全系数的影响由图 8 可知,降雨期间,FOS 减少。在给定的降雨强度下,零内聚力导致 5h 后边坡失稳,FOS小于 1.3,而 70h 的降雨导致 5kPa 内聚力的边坡发生破坏。但对于更高的内聚力,FOS 的下降速度平缓。图 9 降雨 70h 后边坡破坏面由图 9 可知,对于所有黏聚力的情况,孔隙压力的分布都是相似的,表明在所有情况下,渗入土壤的水量是相同的。降雨 70h 后相应内聚值的边坡破坏面,在斜坡上降雨 100h(1.08 cm/h)后,20和 30kPa 内聚力下获得的 FOS 仍然高于 1.8。4 结 论 通过数值模拟,本文研究了几个水文参数和抗剪强度参数对锦水河典型土边坡降雨诱发边坡失稳的影响,分析了饱和度、渗透率和内聚力在控制初始 FOS 和渗透过程中的作用。结果表明,饱和度、渗透性和降雨强度对水渗入土壤和负孔隙压力分布有影响。随着渗透性变高,大量雨水渗入土壤,进而增加正孔隙压力。所有参数(即降雨强度、降雨分布和饱和渗透系数)之间存在显著的相互作用,表明在对非饱和边坡进行瞬态分析时,需要对土壤的水力特性进行适当的选择。参考文献1 王花蕊.降雨诱发黄土边坡失稳破坏机制研究J.陕西水利,2019(2):27-28,31.2 陈乔,徐烽淋,朱洪林,等.基于降雨滑坡模型的两种边坡破坏模式研究J.水利水电技术,2018,49(4):138-144.3 刘平军.降雨条件下顺层岩质边坡破坏模式及稳定性研究J.交通科技,2015(5):80-83.4 刘亚龙.关于边坡浅层滑动的影响因素及研究J.西部探矿工程,2022,34(2):20-22,26.5 蒋庆华.锦水河年径流量及设计洪水计算探讨J.内蒙古水利,2021(10):25-27.41