Influence of hydraulic-mechanical coupling and dry-wet cycle effect on surficial layer stability of residual soil slopes
-
摘要: 降雨入渗过程中，土体吸力降低，体积明显改变。天然浅层土体长期受到季节性气候变化的影响，因此，开展水-力耦合及干湿交替对浅层残积土坡稳定性影响的数值分析，分析浅层土坡孔隙水压力、湿润锋及安全系数的时空演变规律，并对水-力耦合及干湿交替条件下的浅层土坡失稳破坏机制进行探讨显得尤为必要。研究结果表明：随着干湿循环次数的增加，水-力耦合分析下孔隙水压力以及湿润锋的迁移速度增加更快，边坡也更易失稳破坏；干湿交替初期，雨水入渗易引起地下水位上升，边坡可因正孔隙水压力的增加而失稳；干湿交替后期，湿润锋的快速推进加剧基质吸力迅速丧失及土体强度下降，边坡安全系数显著降低，发生失稳破坏的时间缩短。因此，可将湿润锋处的安全系数(局部最小值)作为控制边坡长期稳定性的临界值。Abstract: During rainfall infiltration, soil suction decreased and volume changed significantly. Since the natural shallow soil was under the influence of seasonal climate changes for a long time, numerical analysis of the influence of hydraulic-mechanical coupling and dry-wet alternation on the stability of shallow residual soil slopes was carried out. The temporal and spatial evolution law of pore water pressure, wetting front and safety factor of shallow soil slopes were analyzed deeply, then, the failure mechanisms of soil slopes under hydraulic-mechanical coupling and alternating dry-wet conditions were further discussed. Results show that with the increase of dry-wet cycles, the migration velocity of wetting front and pore water pressure increased more quickly, and the slope was more unstable under hydro-mechanical coupling analysis. At the early stage of the dry-wet cycle, the infiltration of rainwater would easily cause the groundwater level to rise, and the slope might lose its stability due to the increase of positive pore water pressure. In the later stage of dry-wet cycle, the rapid advance of wetting front accelerated the rapid loss of matrix suction and the decrease of soil strength, the safety factor of slope was significantly reduced and the time of failure was shorter. Therefore, the safety factor (local minimum) at the wetting front could be used as the critical value to control the long-term stability of the slope.
-
0. 引言
福建山区公路建设中存有大量的高陡边坡，路基的变形、稳定及病害处治等备受关注[1]。天然浅层土体长期受季节性气候变化的影响，土体原生结构发生改变，导致其抵抗外荷载的能力显著下降[2]。残积土滑坡作为福建省最主要的地质灾害，易受降雨的影响。雨水入渗过程中，因边坡内部水/气两相流发生运动，导致土体的微观结构、力学特性和有效应力状态发生改变，弱化土体的持水性能、渗流特性及边坡稳定性，引起坡体不可逆变形的发生以及滑面的形成与贯通。因此，开展水-力耦合及干湿交替对浅层非饱和残积土坡稳定性影响的研究，对此类路基边坡灾害防治与减灾研究具有重要的科学意义和现实意义。
非饱和土中的基质吸力对土的抗剪强度及路基边坡的稳定性有贡献。雨季浅层土坡失稳与湿润峰迁移引起的吸力损失有关，降雨条件下边坡暂态饱和区分布特征与演化规律以及地下水位变动所引起的坡体内基质吸力减小、抗剪强度降低和孔隙水压力提高是影响边坡稳定性的主要因素[3-8]。Song等[9]基于一维降雨入渗模拟试验，研究不同降雨强度作用下边坡内部岩土体基质吸力变化规律；曾铃等[10]研究了降雨入渗条件下坡体暂态饱和区的形成条件及多因素影响下暂态饱和区的分布；邱祥等[11]基于饱和-非饱和渗流理论，分析了不同降雨强度、饱和渗透系数、表面吸力、坡度条件下边坡降雨入渗影响区与暂态饱和区变化规律。
部分学者采用水-力耦合方法来考虑渗流过程中的土体变形。郑开欢等[8]建立温度场-渗流场-应力场耦合有限元模型，结合排土场边坡渗流场、应力场以及稳定性响应规律，提出基于边坡危险滑动带来讨论非均质边坡失稳特征。张良以等[12]基于膨胀应变与基质吸力增量的线性关系，在非饱和流-固耦合模型中引入膨胀性指标，建立一种适用于膨胀土工程的非饱和渗流场-应力场-膨胀应变场的多场耦合数值计算方法；熊勇林等[13]建立了一个基于有限元-有限差分理论的水-土-气三相渗流-变形耦合有限元数值分析程序，该程序可反映降雨入渗引起的边坡变形，边坡破坏先由坡趾面开始，逐渐向坡顶延伸最终形成连贯塑性剪切带。
以上研究都是基于非饱和土力学理论开展相关的渗流分析和水-力耦合分析，但对干湿循环效应下土体的变形和耦合行为对残积土边坡稳定性的影响仍需进一步解释。因此，本文利用Geo-studio软件的SIGMAR/W模块和SEEP/W模块模拟长时小雨下典型残积土边坡的水-力耦合及渗流效应，重点分析水-力耦合及干湿循环条件下坡体内部孔隙水压力、湿润锋及边坡安全系数的时空演变规律，以期为残积土边坡的长期稳定性评估提供参考。
1. 浅层滑坡稳定性计算方法
降雨诱发浅层土坡失稳的滑面常平行于斜坡表面，考虑滑面具有较大的长深比，可进一步假定破坏的土块沿滑动方向无限延伸，按无限边坡进行分析。因此，可将斜坡的安全系数定义为土体的抗剪强度与滑面上的下滑应力的比值。
非饱和土抗剪强度中考虑用体积含水率函数来描述基质吸力对强度的贡献，形式如下[14]：
$$ \tau _{\text{f}} = {c\text{'}} + {\text{(}}{\sigma _{\text{n}}} - {u_{\text{a}}}{\text{)tan}}{\varphi \text{'}} + {\text{(}}{u_{\text{a}}}{ - }{u_{\text{w}}})\left[\left(\frac{{{\theta _{\text{w}}} - {\theta _{\text{r}}}}}{{{\theta _{\text{s}}} - {\theta _{\text{r}}}}}\right)\tan{\varphi \text{'}}\right] $$ (1) 式中：
$\tau _{\text{f}}$ ——土体抗剪强度/kPa；σn——法向应力/kPa；
ua、uw——孔隙气压力/kPa和孔隙水压力/kPa；
c'——有效黏聚力/kPa；
φ'——土体的有效内摩擦角/（°）；
θw、θr、θs——土体的体积含水量、残余体积含水率 和饱和体积含水率。
根据式（1）及图1的计算简图[15]，无限边坡安全系数可表示为：
$$ {F_s} = \frac{{{c\text{'}}}}{{\gamma {\textit{z}}\cos \beta \sin \beta }} + \frac{{{\text{tan}}{\varphi\text{'}}}}{{\tan \beta }} + \frac{{({u_{\text{a}}}{{ - }}{u_{\text{w}}})\left(\dfrac{{{\theta _{\text{w}}} - {\theta _{\text{r}}}}}{{{\theta _{\text{s}}} - {\theta _{\text{r}}}}}\right)\tan{\varphi\text{'}}}}{{\gamma {\textit{z}}\cos \beta \sin \beta }} $$ (2) 式中：γ——土体重度/（kN·m−3）；
β——坡角/（°）；
z——坡体深度/m。
公式（2）表明边坡安全系数由与土体有效黏聚力、有效内摩擦角和基质吸力相关的变量组成。饱和状态下，
$ u_w$ ≥0，且$({\theta _{\text{w}}} - {\theta _{\text{r}}})/({\theta _{\text{s}}} - \theta _{\text{r}})$ =1，此时，土体非饱和状态下的抗剪强度向饱和状态时的抗剪强度转变，有效应力准则仍然适用。因此，对于处于水位以下的滑面仍然可以采用公式（2）对其进行边坡稳定性分析。水-力耦合分析采用Geo-Studio中的Seep/W及Sigma/W模块完成。Seep/W基于加权余量的伽辽金有限元法，可对土体稳态或瞬态渗流进行分析，降雨条件下坡体内渗流场时刻改变，涉及饱和区和非饱和区的土体，需考虑瞬态渗流分析。Sigma/W可以使用一个完全耦合来求解土的固结问题，完全耦合分析需要对每个网格中的节点求解应力-变形和渗流连续性方程，以此给出位移和孔隙水压力的变化。Seep/W和Sigma/W数值分析所需的相关方程可参考相应手册[16-17]。
2. 工程概况及计算参数
研究区位于福州福清福政大道至福人大道公路工程，该地段属剥蚀低丘，沿路线呈波状起伏展布，间夹宽缓的沟谷。该填方路堤段在原有基础上最大填筑高度约7 m，最终可形成一高12 m左右的高边坡，坡率为1∶1.5，钻孔资料显示该地段以残积黏性土为主，地下水位线位于2 m高程，该填方边坡典型剖面如图2所示。现场及室内模型试验表明[1, 18]：土坡失稳破坏的范围集中在1.5～3.0 m。因此，选择0～3 m的浅层土坡作为分析对象，并将边坡简化为单一土层。
2.1 计算参数
X射线衍射表明残积黏性土体以石英、高岭石、钾长石等为主，残积土的颗粒级配曲线及孔径分布曲线见图3，土体液、塑限分别为43%、24%。干湿循环条件下残积土土-水特征曲线（Soil WaterCharacteristic Curve，SWCC）可由应力相关的SWCC压力板仪系统测得（图4），天然浅层土体含水率变化大，现场勘查显示残积土质量含水率在17%至饱和含水率之间。残积土SWCC测定时，先对抽真空饱和后的原状土加压至1000 kPa（此时对应的含水率约为17%），待试样排水稳定后，逐渐降低气压，此时土体缓慢吸水，重复上述步骤，直至完成目标干湿循环次数（1、2、3、5次循环）。整个干湿循环试验过程中，记录数据频次为每隔12 h一次，当24 h内试样吸水或排水的体积不超过土样体积的0.05%时，可认为吸力达到平衡[19]。渗流分析所需的渗透系数函数（HCF）由室内试验成果及参考类似工程[20-21]综合确定（图4），数值分析所需的参数见表1。土体SWCC曲线受孔径分布的影响，浅层土体长期处于湿热环境，土体裂隙发育，土体的孔径分布曲线（图3）表明其为双孔隙结构，土体内部所包含的大孔隙（裂隙）和小孔隙（土颗粒胶结孔隙）导致SWCC曲线呈典型的双峰特性（图4）。既有研究表明[22-26]：随着干湿交替的进行，土体饱和渗透系数（ksat）逐渐增加、有效黏聚力（c'）波动下降，而有效内摩擦角（φ'）则基本保持不变。由文献[19，22，23，25]及室内试验可得考虑水-力耦合及干湿循环的数值计算参数（表2），不同干湿循环次数数值计算所需的HCF函数与图4形式一致，这里仅考虑因干湿交替引起土体开裂而导致饱和渗透系数的变化。
![]()
图 3 残积土颗粒级配曲线与孔径分布曲Figure 3. Particle gradation curve and pore size distribution curve of residual soil表 1 数值计算所需参数取值Table 1. Parameter value required for numerical calculation性质 单位 数值 重度γ kN·m−3 18.5 杨氏模量 Esat kPa 2300 泊松比 μ − 0.4 有效强度指标 c' kPa 16 φ' ° 15 饱和渗透系数 ksat m·d−1 1.8×10−3 饱和体积含水率 θs − 0.51 残余体积含水率 θr − 0.15 表 2 考虑水-力耦合及干湿循环的数值计算方案Table 2. Numerical calculation scheme considering hydraulic-coupling and dry-wet cycle干湿循环次数 分析类型 SWCC
（图4）ksat
/（m·d−1）c'
/kPaφ'
/kPa1 耦合 曲线1 1.8×10−3 16 15 2 耦合 曲线2 2.2×10−3 14 3 耦合 曲线3 2.4×10−3 12 5 耦合 曲线4 2.8×10−3 10 数值分析所需的非饱和条件下杨氏模量需通过构建孔隙比与双应力状态变量之间的关系曲面确定，但对于文中所研究的低渗透性的黏性土，建立上述曲面需耗费大量的时间。因此，文中基于SWCC曲线，采用文献[27]的半经验模型来预测杨氏模量随吸力的非线性变化，具体如下：
$$ {E_{{\text{unsat}}}} = {E_{{\text{sat}}}}\left[ {1 + \alpha \frac{{({u_{\text{a}}} - {u_{\text{w}}})}}{{{P_{\text{a}}}/101.3}}S_{\text{r}}^{\beta}} \right] $$ (3) 式中：Eunsat——非饱和条件下的杨氏模量/kPa；
Esat——饱和条件下的杨氏模量/kPa；
ua−uw——基质吸力/kPa；
Pa——用于维持方程左右两端单位一致的大气压力/101.3 kPa；
Sr——土体饱和度；
α、β——拟合参数。
根据文献[28]给出的不同土体特性及应力条件下拟合参数α、β的取值，低应力条件下拟合参数α介于0.05～0.15。因此，综合考虑土体塑性指数（文中为17）及文献[28]的经验取值建议，当拟合参数α和β定为0.1和1.2时，可较好的估算残积黏性土的杨氏模量（图5）。
![]()
图 5 半经验模型下土体Eunsat随吸力的变化曲线Figure 5. Variation curve of soil Eunsat with suction under semi-empirical model由公式（3）可知：耦合分析中所用的杨氏模量具有时空效应，该模量随深度及降雨入渗时间的不同而不同。目前的SIGMA/W商业软件无法直接定义Eunsat随吸力的变化，Eunsat随吸力变化的精确解可通过数值积分技术（如：Euler积分求解）来获取，这里采用类似于有限元数值分析中非线性材料的切线刚度方法来实现Eunsat随基质吸力的变化。首先，将浅层土体分割成多个细土层（这里划分每个土层厚度为0.15 m，见图2）；其次，利用公式（3）计算上述各划分土层孔隙水压力时空分布曲线所对应的杨氏模量Eunsat的大小（图6）；最后，将图6的不同降雨时刻下Eunsat随基质吸力的变化结果更新于耦合分析中，以上耦合非线性方程的近似增量求解方法可较为准确的体现Eunsat与基质吸力的关系。为便于后续分析，本文选取的计算点有（图2）：上部、中部和下部土层，坐标范围依次为：（16，24）～（16，21）、（25，18）～（25，15）、（34，12）～（34，9）。
![]()
图 6 中部土体Eunsat随时间和深度的变化曲线Figure 6. Variation curve of soil Eunsat with time and depth2.2 计算工况及边界条件
依据气象部门对研究区的月降水量统计（图7），该地区年总降雨量为1530 mm，降水主要集中于梅雨季和台风季，其中每年的5—8月（降雨量约占年总降雨量的55 %）为滑坡、崩塌、泥石流易发时段，尤其是经历长时梅雨季节后（前期有效降雨），后续更易出现台风暴雨诱发型滑坡[19，29]。因此，将年平均降雨量的55%（雨强大小约为2.4×10−3 m/d）和3个月持时（梅雨季节为3—5月）作为耦合（考虑干湿循环效应）分析的流量边界（简化雨强略大于土体饱和渗透系数，按流量边界考虑，流量为简化雨强大小）。模型两侧饱和区（地下水位线以下范围）作为水头边界条件，模型两侧非饱和区为零流量边界，模型底部为不透水边界。长时小雨过程中，可能在坡表处形成暂态饱和区，即：地下水渗流可能沿着坡表溢出，因此坡表可定为渗流面。
3. 数值计算结果分析
3.1 干湿循环及耦合条件对孔隙水压力的影响
图8为考虑耦合及干湿循环次数下浅层土体孔隙水压力随时间和深度的变化趋势图。长时小雨条件下，边坡上部土体吸力显著下降，局部易出现暂态饱和区。降雨60 d后（图8），不同干湿循环次数下的边坡表面负孔隙水压力均增加至零，而坡体下部土层仍保留部分静水条件下的吸力。同一深度下（如：深度1.95 m处），干湿循环次数N=1，2，3，5对应的土体吸力分别为（图8）：81.5 kPa、56 kPa、26.5 kPa、10.2 kPa，这表明高干湿循环次数下湿润表层土体所需的时间更少。
![]()
图 8 干湿循环对浅层土体孔隙水压力空间分布的影响Figure 8. Influence of dry-wet cycles on the spatial distribution of pore water pressure in shallow soils60 d降雨后，中部土体在1.95 m处的体积含水率变化更快，此时干湿循环次数N=1，2，3，5对应的土体体积含水率依次为：0.3182，0.3118，0.3253，0.3786（图9），由于干湿交替引起土体内部结构损伤，导致不同体积含水率对应的土体基质吸力不一样，如N=1，2时，虽然第二次的体积含水率小于第一次，但第一次干湿循环下土体的内部吸力却大于第二次干湿循环；随着干湿循环的推进，边坡土体逐渐软化，坡体内部体积进一步改变，从而引起土体吸力大幅度下降，导致湿润锋推进得更快（图9）。因此，在同一时间段内，随着干湿交替的进行，湿润锋的快速推进可引起越深位置的坡体孔隙水压力发生改变。
![]()
图 9 干湿循环对浅层土体体积含水率空间分布的影响Figure 9. Influence of dry-wet cycles on the spatial distribution of shallow soil volumetric water content3.2 耦合条件下干湿循环次数对湿润锋的影响
图10为不同降雨历时及不同干湿循环次数下土层湿润锋推进变化曲线。随着降雨的持续进行，雨水可逐渐入渗至土层深处，如干湿循环次数N=1，5时，整个降雨周期内湿润锋的位置依次为：0.45 m→0.75 m→1.2 m→1.35 m→1.65 m→1.95 m、0.75 m→1.2 m→1.65 m→1.95 m→2.4 m→2.55 m。干湿循环次数下越大，整个降雨周期内湿润锋可推进至更深位置，如：降雨历时为75 d时，干湿循环次数N=1，5时湿润锋的值依次为：1.65 m、2.4 m，增加幅度为45.5%。不同降雨持时后（15 d、30 d、60 d、75 d、90 d）干湿循环，湿润锋推进范围内的土体基本达到饱和状态，如整个降雨和干湿循环周期后，坡体的暂态饱和区域由0～0.45 m发展至0～2.55 m；而湿润锋前沿土体的吸力损失随着降雨和干湿循环的进行而快速丧失，土体含水率逐渐处于较高水平。
![]()
图 10 干湿循环对湿润锋推进的影响Figure 10. Influence of the dry-wet cycle on the wet front propulsion图11为不同降雨时刻下中部土层湿润锋处基质吸力随深度的变化曲线，Line1和Line2分别为降雨历时60 d和30 d时湿润锋的推进曲线，Line3为同一深度处（即：降雨历时60 d湿润锋所处的位置）土体的基质吸力值。同一时间下，干湿循环次数下土体润湿前深度所具有的基质吸力较低（Line1和Line2）；同一湿润锋深度下，随着降雨时间的推进以及干湿循环次数的增加，基质吸力下降的时空分布更为显著（Line2和Line3），如：湿润锋为1.75 m、1.95 m时，坡前吸力变化分别为：91.3 kPa→6.2 kPa、75.1 kPa→10.2 kPa。综上可知，雨水入渗速度与干湿循环次数有关，随着干湿循环次数的增加，湿润锋前沿更深，因湿润锋深度附近的基质吸力下降可导致安全系数进一步降低。
![]()
图 11 不同降雨时刻土层湿润锋处吸力随深度的变化曲线Figure 11. Change curve of suction at the wet front of the soil layer with depth at different rainfall moments3.3 耦合条件下干湿循环对边坡安全系数的影响
图12为不同降雨时刻下干湿循环次数对边坡安全系数的影响。整体上看，同一干湿循环次数下，湿润锋处的边坡安全系数随着降雨的进行而持续降低，但降幅有所减小，如：第1次干湿循环，整个降雨周期内边坡湿润锋处的安全系数分别为：4.92→3.11→2.24→1.96→1.66→1.53，降幅依次为：36.8%、27.9%、12.3%、15.5%、7.8%。同一降雨时刻下（降雨初期），湿润锋处的边坡安全系数随干湿循环次数的增加呈持续下降趋势，但在降雨后期，边坡安全系数降幅减小，如：降雨15 d、90 d，干湿循环次数N=1，2，3，5对应湿润锋处的安全系数分别为：4.92→4.26→2.94→2.16、1.53→1.28→1.12→1.02，下降幅度依次为：13.4%→30.9%→26.5%、16.3%→12.7%→9.0%，尤其是第5次干湿循环后，边坡安全系数为1.02（湿润锋推进的位置为2.55 m），处于临界失稳状态。
降雨60 d时，不同深度的边坡安全系数随干湿循环次数的变化曲线可见图13。由图13可知：上部土层因应力水平低而具有相对较高的安全系数，随着降雨的持续进行，湿润锋逐渐向坡体内部推进，安全系数随之下降；受干湿交替的影响，雨水入渗引起湿润锋的推进速度和深度更快，直至可能出现局部最低的安全系数。值得注意的是：第1次干湿循环时，土体湿润锋处的安全系数为1.964，但整个浅层边坡的最低安全系数仍然处于最深处，由于土层上覆应力的影响，深部土体的剪应力更高，当雨水入渗不足以引起加剧土体强度的迅速下降时，此时边坡的安全系数由土体所具有的初期应力水平决定；但随着干湿交替的加剧，此时因雨水入渗引起深部土体基质吸力的大幅度下降，导致边坡的安全系数可在湿润锋的推进处达到最低值。
4. 讨论
Li等[30]指出降雨入渗下，只有当土体的有效内摩擦角小于边坡坡角时，湿润锋的迁移才可能导致浅层边坡的失稳破坏。当边坡坡角增加时，边坡浅层土体的稳定性会显著下降，即：陡坡可加大浅层滑坡的发生。由于湿润锋推进的过程中孔隙水压力为负值，因此，浅层边坡失稳破坏的条件为有效黏聚力和有效内摩擦角（式（2）的前面两项）对安全系数的贡献小于1。图14为式（2）前两项计算出的安全系数与深度关系曲线以及边坡破坏（Fs=1）所需的孔隙水压力大小，由图14可知：因湿润锋迁移导致边坡失稳所处的位置为2.4 m以下（区域3）。在区域3中，深度越大，边坡破坏所需的基质吸力也越大。在区域2内（深度范围1.65～2.4 m），正的孔隙水压力（地下水位上升）可导致边坡失稳破坏。在区域1内（深度范围0～1.65 m），由于边坡破坏所需的正孔隙水压力大于静水条件下的孔隙水压力，区域1范围内的边坡处于稳定状态。
![]()
图 14 降雨诱发土坡失稳的破坏机制（第3次干湿循环）Figure 14. Failure mechanism of soil slope instability induced by rainfall (third dry-wet cycle)图15为干湿循环次数对降雨诱发土坡失稳机制的影响。第1次干湿循环下，地下水位逐渐上升，此时非饱和区域逐渐向饱和区域过渡，土坡可因边坡内正孔隙水压力的增加而失稳（区域2，即：坡体位置为2.4～2.7 m）；随着干湿循环次数的增加，边坡失稳破坏不仅可因正孔隙水压力增加引起，还可因湿润锋快速推进、土体基质吸力损失而产生（即：区域3），并呈现区域3范围逐渐扩大的特点。第2、3、5次干湿循环下，区域3的范围依次为2.7～3.0 m、2.4～3.0 m、1.95～3.0 m，区域2的范围依次为：1.8～2.7 m、1.65～2.4 m、1.35～1.95 m，区域1的范围分别为：0～1.8 m、0～1.65 m、0～1.35 m，区域1和区域2的范围随干湿循环次数的增加而减小，表明边坡的长期稳定性应由湿润锋的迁移所引起的基质吸力和土体抗剪强度损失所控制，此时湿润锋前缘的安全系数为过渡区域内的局部最小值。因此，可以将其作为控制边坡长期稳定性的临界值。
![]()
图 15 干湿循环次数下降雨诱发土坡失稳机制Figure 15. The mechanism of soil slope instability induced by rainfall under the effects of dry-wet cycles综上所述，当湿润锋传播至区域3时，文中所分析的边坡才可能发生破坏；此外，破坏滑动面越深，所需的负孔隙水压力越小。因此，湿润锋迁移的速率对于浅层土坡的稳定性起着关键性的作用，湿润锋前沿越早入渗至区域3，土坡失稳所需的时间就越短。以上分析进一步解释了耦合效应及干湿循环次数对边坡长期稳定性的不利影响，单一考虑水力耦合效应而不考虑干湿交替的影响会过低的估计湿润锋的传播速度，导致具有更高的安全系数。同时，也应注意到土体所处实际湿热（或风干）环境，易引起其表层产生干缩裂缝，有利于雨水入渗，因此，采用压力板仪试验（即：封闭环境）获取的残积土干湿循环水力特性数据与实际有出入，由此计算得到的边坡浅层稳定性会较高。
5. 结语
（1）第1次干湿循环下，坡体深处因具有高水平的剪应力而具有更低的安全系数。降雨过程中，湿润锋推进深度内边坡安全系数随着深度的增加而减小，与数值分析得到的孔隙水压力变化一致。
（2）随着干湿循环次数的增加，湿润锋的迁移速度进一步加快，孔隙水压力显著增加，水-力耦合分析下浅层残积土边坡更易失稳。干湿交替初期边坡可因正孔隙水压力(雨水入渗引起的地下水位上升)的增加而失稳，而干湿交替后期边坡可因湿润锋的快速推进(吸力损失)而破坏。
（3）干湿循环后期，耦合分析条件下的湿润锋迁移速率更快，土坡失稳所需的时间越短。边坡的长期稳定性受湿润锋迁移所引起的吸力和强度下降控制，可将湿润锋前缘的安全系数（局部最小值）作为评估边坡长期稳定性的临界值。
（4）文中虽然考虑了干湿循环次数的增加引起土体的开裂，进而导致饱和渗透系数变化，但未考虑体积变化对土体水力传导特性的影响。降雨入渗过程，土体孔隙比的变化可增加土体的透水性。因此，在耦合分析中，湿润峰前沿深度的增长速度可能更快。此外，可进一步开展降雨-蒸发循环条件（实际风干环境下）残积土边坡破坏机理试验研究，结合相似理论及数值模拟结果，对比分析湿润锋迁移速率及其引起基质吸力和强度下降规律，明确水力耦合及干湿循环条件下浅层残积土边坡的失稳机制。
-
![]()
图 3 残积土颗粒级配曲线与孔径分布曲
Figure 3. Particle gradation curve and pore size distribution curve of residual soil
![]()
图 5 半经验模型下土体Eunsat随吸力的变化曲线
Figure 5. Variation curve of soil Eunsat with suction under semi-empirical model
![]()
图 6 中部土体Eunsat随时间和深度的变化曲线
Figure 6. Variation curve of soil Eunsat with time and depth
![]()
图 8 干湿循环对浅层土体孔隙水压力空间分布的影响
Figure 8. Influence of dry-wet cycles on the spatial distribution of pore water pressure in shallow soils
![]()
图 9 干湿循环对浅层土体体积含水率空间分布的影响
Figure 9. Influence of dry-wet cycles on the spatial distribution of shallow soil volumetric water content
![]()
图 10 干湿循环对湿润锋推进的影响
Figure 10. Influence of the dry-wet cycle on the wet front propulsion
![]()
图 11 不同降雨时刻土层湿润锋处吸力随深度的变化曲线
Figure 11. Change curve of suction at the wet front of the soil layer with depth at different rainfall moments
![]()
图 14 降雨诱发土坡失稳的破坏机制（第3次干湿循环）
Figure 14. Failure mechanism of soil slope instability induced by rainfall (third dry-wet cycle)
![]()
图 15 干湿循环次数下降雨诱发土坡失稳机制
Figure 15. The mechanism of soil slope instability induced by rainfall under the effects of dry-wet cycles
表 1 数值计算所需参数取值
Table 1 Parameter value required for numerical calculation
性质 单位 数值 重度γ kN·m−3 18.5 杨氏模量 Esat kPa 2300 泊松比 μ − 0.4 有效强度指标 c' kPa 16 φ' ° 15 饱和渗透系数 ksat m·d−1 1.8×10−3 饱和体积含水率 θs − 0.51 残余体积含水率 θr − 0.15 
下载: 导出CSV
-
[1] 蔡沛辰,阙云,李显. 非饱和花岗岩残积土水-气两相驱替过程数值模拟[J]. 水文地质工程地质,2021,48(6):54 − 63. [CAI Peichen,QUE Yun,LI Xian. Numerical simulation of water-gas two-phase displacement process in unsaturated granite residual soil[J]. Hydrogeology & Engineering Geology,2021,48(6):54 − 63. (in Chinese with English abstract) DOI: 10.16030/j.cnki.issn.1000-3665.202010017 CAI Peichen, QUE Yun, LI Xian. Numerical simulation of water-gas two-phase displacement process in unsaturated granite residual soil[J]. Hydrogeology & Engineering Geology, 2021, 48（6）: 54-63. （in Chinese with English abstract）]https: //kns. cnki. net/KCMS/detail/detail. aspx?filename=SWDG202106006&dbname=CJFD&dbcode=CJFQhttps: //oversea.cnki.net/KCMS/detail/detail.aspx?filename=SWDG202106006&dbname=CJFD&dbcode=CJFQ DOI: 10.16030/j.cnki.issn.1000-3665.202010017
[2] ZENG H,TANG C S,CHENG Q,et al. Coupling effects of interfacial friction and layer thickness on soil desiccation cracking behavior[J]. Engineering Geology,2019,260:105220. DOI: 10.1016/j.enggeo.2019.105220
[3] 蔡荣坤,戴自航,徐根连,等. 降雨对花岗岩风化层路堑边坡滑动模式影响—以福建云平高速云霄段为例[J]. 中国地质灾害与防治学报,2021,32(2):27 − 35. [CAI Rongkun,DAI Zihang,XU Genlian,et al. Influence of rainfall on sliding modes of cutting slope of weathered granite stratum:taking Yunxiao section in the Yunping freeway in Fujian for example[J]. The Chinese Journal of Geological Hazard and Control,2021,32(2):27 − 35. (in Chinese with English abstract) CAI Rongkun, DAI Zihang, XU Genlian, et al. Influence of rainfall on sliding modes of cutting slope of weathered granite stratum: taking Yunxiao section in the Yunping freeway in Fujian for example[J]. The Chinese Journal of Geological Hazard and Control, 2021, 32（2）: 27-35. （in Chinese with English abstract）]https: //kns. cnki. net/KCMS/detail/detail. aspx?filename=ZGDH202102004&dbname=CJFD&dbcode=CJFQhttps: //oversea.cnki.net/KCMS/detail/detail.aspx?filename=ZGDH202102004&dbname=CJFD&dbcode=CJFQ
[4] TANG G P,HUANG J S,SHENG D C,et al. Stability analysis of unsaturated soil slopes under random rainfall patterns[J]. Engineering Geology,2018,245:322 − 332. DOI: 10.1016/j.enggeo.2018.09.013
[5] WANG Y X,CHAI J R,CAO J,et al. Effects of seepage on a three-layered slope and its stability analysis under rainfall conditions[J]. Natural Hazards,2020,102(3):1269 − 1278. DOI: 10.1007/s11069-020-03966-1
[6] 饶鸿,王金淑,赵志明,等. 基于有限元软件自定义本构模型的膨胀土边坡降雨入渗分析[J]. 水文地质工程地质,2021,48(1):154 − 162. [RAO Hong,WANG Jinshu,ZHAO Zhiming,et al. An analysis of rainfall infiltration of expansive soil slope based on the finite element software custom constitutive model[J]. Hydrogeology & Engineering Geology,2021,48(1):154 − 162. (in Chinese with English abstract) RAO Hong, WANG Jinshu, ZHAO Zhiming, et al. An analysis of rainfall infiltration of expansive soil slope based on the finite element software custom constitutive model[J]. Hydrogeology & Engineering Geology, 2021, 48（1）: 154-162. （in Chinese with English abstract）]https: //kns. cnki. net/KCMS/detail/detail. aspx?filename=SWDG202101019&dbname=CJFD&dbcode=CJFQhttps: //oversea.cnki.net/KCMS/detail/detail.aspx?filename=SWDG202101019&dbname=CJFD&dbcode=CJFQ
[7] ZENG L,BIAN H B,SHI Z N,et al. Forming condition of transient saturated zone and its distribution in residual slope under rainfall conditions[J]. Journal of Central South University,2017,24(8):1866 − 1880. DOI: 10.1007/s11771-017-3594-6
[8] 郑开欢,罗周全,江宏. 天气因素对排土场生态边坡稳定性的影响[J]. 中国地质灾害与防治学报,2018,29(6):97 − 102. [ZHENG Kaihuan,LUO Zhouquan,JIANG Hong. Weather factors' influence on the stability of ecological slopes of waste dump[J]. The Chinese Journal of Geological Hazard and Control,2018,29(6):97 − 102. (in Chinese with English abstract) ZHENG Kaihuan, LUO Zhouquan, JIANG Hong. Weather factors' influence on the stability of ecological slopes of waste dump[J]. The Chinese Journal of Geological Hazard and Control, 2018, 29（6）: 97-102. （in Chinese with English abstract）]https: //kns. cnki. net/KCMS/detail/detail. aspx?filename=ZGDH201806014&dbname=CJFD&dbcode=CJFQhttps: //oversea.cnki.net/KCMS/detail/detail.aspx?filename=ZGDH201806014&dbname=CJFD&dbcode=CJFQ
[9] SONG S,BROCCA L,WANG W,et al. Testing the potential of soil moisture observations to estimate rainfall in a soil tank experiment[J]. Journal of Hydrology,2020,581:124368. DOI: 10.1016/j.jhydrol.2019.124368
[10] 曾铃,史振宁,付宏渊,等. 降雨入渗对边坡暂态饱和区分布特征的影响[J]. 中国公路学报,2017,30(1):25 − 34. [ZENG Ling,SHI Zhenning,FU Hongyuan,et al. Influence of rainfall infiltration on distribution characteristics of slope transient saturated zone[J]. China Journal of Highway and Transport,2017,30(1):25 − 34. (in Chinese with English abstract) DOI: 10.3969/j.issn.1001-7372.2017.01.004 ZENG Ling, SHI Zhenning, FU Hongyuan, et al. Influence of rainfall infiltration on distribution characteristics of slope transient saturated zone[J]. China Journal of Highway and Transport, 2017, 30（1）: 25-34. （in Chinese with English abstract）]https: //kns. cnki. net/KCMS/detail/detail. aspx?filename=ZGGL201701004&dbname=CJFD&dbcode=CJFQhttps: //oversea.cnki.net/KCMS/detail/detail.aspx?filename=ZGGL201701004&dbname=CJFD&dbcode=CJFQ DOI: 10.3969/j.issn.1001-7372.2017.01.004
[11] 邱祥,蒋煌斌,欧健,等. 降雨条件下边坡暂态饱和区形成条件与演化特征数值分析[J]. 水利学报,2020,51(12):1525 − 1535. [QIU Xiang,JIANG Huangbin,OU Jian,et al. Numerical analysis of formation conditions and evolution characteristics of transient saturation zone of a slope under rainfall conditions[J]. Journal of Hydraulic Engineering,2020,51(12):1525 − 1535. (in Chinese with English abstract) DOI: 10.13243/j.cnki.slxb.20200254 QIU Xiang, JIANG Huangbin, OU Jian, et al. Numerical analysis of formation conditions and evolution characteristics of transient saturation zone of a slope under rainfall conditions[J]. Journal of Hydraulic Engineering, 2020, 51（12）: 1525-1535. （in Chinese with English abstract）]https: //kns. cnki. net/KCMS/detail/detail. aspx?filename=SLXB202012010&dbname=CJFD&dbcode=CJFQhttps: //oversea.cnki.net/KCMS/detail/detail.aspx?filename=SLXB202012010&dbname=CJFD&dbcode=CJFQ DOI: 10.13243/j.cnki.slxb.20200254
[12] 张良以,陈铁林,张顶立. 降雨诱发膨胀土边坡渐进破坏研究[J]. 岩土工程学报,2019,41(1):70 − 77. [ZHANG Liangyi,CHEN Tielin,ZHANG Dingli. Progressive failure of expansive soil slopes under rainfall[J]. Chinese Journal of Geotechnical Engineering,2019,41(1):70 − 77. (in Chinese with English abstract) ZHANG Liangyi, CHEN Tielin, ZHANG Dingli. Progressive failure of expansive soil slopes under rainfall[J]. Chinese Journal of Geotechnical Engineering, 2019, 41（1）: 70-77. （in Chinese with English abstract）]https: //kns. cnki. net/KCMS/detail/detail. aspx?filename=YTGC201901010&dbname=CJFD&dbcode=CJFQhttps: //oversea.cnki.net/KCMS/detail/detail.aspx?filename=YTGC201901010&dbname=CJFD&dbcode=CJFQ
[13] 熊勇林,朱合华,叶冠林,等. 降雨入渗引起非饱和土边坡破坏的水-土-气三相渗流-变形耦合有限元分析[J]. 岩土力学,2017,38(1):284 − 290. [XIONG Yonglin,ZHU Hehua,YE Guanlin,et al. Analysis of failure of unsaturated soil slope due to rainfall based on soil-water-air seepage-deformation coupling FEM[J]. Rock and Soil Mechanics,2017,38(1):284 − 290. (in Chinese with English abstract) DOI: 10.16285/j.rsm.2017.01.036 XIONG Yonglin, ZHU Hehua, YE Guanlin, et al. Analysis of failure of unsaturated soil slope due to rainfall based on soil-water-air seepage-deformation coupling FEM[J]. Rock and Soil Mechanics, 2017, 38（1）: 284-290. （in Chinese with English abstract）]https: //kns. cnki. net/KCMS/detail/detail. aspx?filename=YTLX201701038&dbname=CJFD&dbcode=CJFQhttps: //oversea.cnki.net/KCMS/detail/detail.aspx?filename=YTLX201701038&dbname=CJFD&dbcode=CJFQ DOI: 10.16285/j.rsm.2017.01.036
[14] ZHAI Q,RAHARDJO H,SATYANAGA A,et al. Estimation of unsaturated shear strength from soil-water characteristic curve[J]. Acta Geotechnica,2019,14(6):1977 − 1990. DOI: 10.1007/s11440-019-00785-y
[15] 叶万军,张宇鹏. 长期降雨作用下黄土边坡失稳模型试验[J]. 中国科技论文,2021,16(6):603 − 609. [YE Wanjun,ZHANG Yupeng. Model test study on instability of loess slopes under long-term rainfall[J]. China Sciencepaper,2021,16(6):603 − 609. (in Chinese with English abstract) DOI: 10.3969/j.issn.2095-2783.2021.06.006 YE Wanjun, ZHANG Yupeng. Model test study on instability of loess slopes under long-term rainfall[J]. China Sciencepaper, 2021, 16（6）: 603-609. （in Chinese with English abstract）]https: //kns. cnki. net/KCMS/detail/detail. aspx?filename=ZKZX202106006&dbname=CJFD&dbcode=CJFQhttps: //oversea.cnki.net/KCMS/detail/detail.aspx?filename=ZKZX202106006&dbname=CJFD&dbcode=CJFQ DOI: 10.3969/j.issn.2095-2783.2021.06.006
[16] INTERNATIONAL G S. SEEP/W for finite element seepage analysis, version 3 and 4[M/OL]. GEO-SLOPE International Ltd;. https://www.researchgate.net/publication/37407388_seepw_for_finite_element_seepage_analysis_version_3_and_4
[17] GEOSLOPE international Ltd. Sigma/W user’s guide for stress-deformation analysis. Calgary, Alta: GEO-SLOPE International Ltd; 2007.
[18] 许旭堂,简文彬,吴能森,等. 降雨诱发残积土坡失稳的模型试验[J]. 中国公路学报,2018,31(2):270 − 279. [XU Xutang,JIAN Wenbin,WU Nengsen,et al. Model test of rainfall-induced residual soil slope failure[J]. China Journal of Highway and Transport,2018,31(2):270 − 279. (in Chinese with English abstract) DOI: 10.3969/j.issn.1001-7372.2018.02.029 XU Xutang, JIAN Wenbin, WU Nengsen, et al. Model test of rainfall-induced residual soil slope failure[J]. China Journal of Highway and Transport, 2018, 31（2）: 270-279. （in Chinese with English abstract）]https: //kns. cnki. net/KCMS/detail/detail. aspx?filename=ZGGL201802030&dbname=CJFD&dbcode=CJFQhttps: //oversea.cnki.net/KCMS/detail/detail.aspx?filename=ZGGL201802030&dbname=CJFD&dbcode=CJFQ DOI: 10.3969/j.issn.1001-7372.2018.02.029
[19] 许旭堂, 简文彬. 土坡对降雨入渗的响应及其失稳演变[M]. 北京: 人民交通出版社, 2020 XU Xutang, JIAN Wenbin. Research on the response and failing process of unsaturated soil slope under rainfall infiltration[M]. Beijing: China Communications Press, 2020. （in Chinese）
[20] YANG K H,UZUOKA R,THUO J N,et al. Coupled hydro-mechanical analysis of two unstable unsaturated slopes subject to rainfall infiltration[J]. Engineering Geology,2017,216:13 − 30. DOI: 10.1016/j.enggeo.2016.11.006
[21] HAN Z,VANAPALLI S,ZOU W L. Integrated approaches for predicting soil-water characteristic curve and resilient modulus of compacted fine-grained subgrade soils[J]. Canadian Geotechnical Journal,2016,54(5):646 − 663. DOI: 10.1139/cgj-2016-0349
[22] XU X T,SHAO L J,HUANG J B,et al. Effect of wet-dry cycles on shear strength of residual soil[J]. Soils and Foundations,2021,61(3):782 − 797. DOI: 10.1016/j.sandf.2021.03.001
[23] XU X T,JIAN W B,WU N S,et al. Void ratio-dependent water retention model for a deformable residual clay[J]. International Journal of Geomechanics,2020,20(8):04020131. DOI: 10.1061/(ASCE)GM.1943-5622.0001773
[24] 安然,孔令伟,黎澄生,等. 炎热多雨气候下花岗岩残积土的强度衰减与微结构损伤规律[J]. 岩石力学与工程学报,2020,39(9):1902 − 1911. [AN Ran,KONG Lingwei,LI Chengsheng,et al. Strength attenuation and microstructure damage of granite residual soils under hot and rainy weather[J]. Chinese Journal of Rock Mechanics and Engineering,2020,39(9):1902 − 1911. (in Chinese with English abstract) DOI: 10.13722/j.cnki.jrme.2020.0073 AN Ran, KONG Lingwei, LI Chengsheng, et al. Strength attenuation and microstructure damage of granite residual soils under hot and rainy weather[J]. Chinese Journal of Rock Mechanics and Engineering, 2020, 39（9）: 1902-1911. （in Chinese with English abstract）]https: //kns. cnki. net/KCMS/detail/detail. aspx?filename=YSLX202009017&dbname=CJFD&dbcode=CJFQhttps: //oversea.cnki.net/KCMS/detail/detail.aspx?filename=YSLX202009017&dbname=CJFD&dbcode=CJFQ DOI: 10.13722/j.cnki.jrme.2020.0073
[25] 简文彬,胡海瑞,罗阳华,等. 干湿循环下花岗岩残积土强度衰减试验研究[J]. 工程地质学报,2017,25(3):592 − 597. [JIAN Wenbin,HU Hairui,LUO Yanghua,et al. Experimental study on deterioration of granitic residual soil strength in wetting-drying cycles[J]. Journal of Engineering Geology,2017,25(3):592 − 597. (in Chinese with English abstract) JIAN Wenbin, HU Hairui, LUO Yanghua, et al. Experimental study on deterioration of granitic residual soil strength in wetting-drying cycles[J]. Journal of Engineering Geology, 2017, 25（3）: 592-597. （in Chinese with English abstract）]https: //kns. cnki. net/KCMS/detail/detail. aspx?filename=GCDZ201703003&dbname=CJFD&dbcode=CJFQhttps: //oversea.cnki.net/KCMS/detail/detail.aspx?filename=GCDZ201703003&dbname=CJFD&dbcode=CJFQ
[26] 刘越,陈东霞,王晖,等. 干湿循环下考虑裂隙发育的残积土边坡响应分析[J]. 岩土力学,2021,42(7):1933 − 1943. [LIU Yue,CHEN Dongxia,WANG Hui,et al. Response analysis of residual soil slope considering crack development under drying-wetting cycles[J]. Rock and Soil Mechanics,2021,42(7):1933 − 1943. (in Chinese with English abstract) LIU Yue, CHEN Dongxia, WANG Hui, et al. Response analysis of residual soil slope considering crack development under drying-wetting cycles[J]. Rock and Soil Mechanics, 2021, 42（7）: 1933-1943. （in Chinese with English abstract）]https: //kns. cnki. net/KCMS/detail/detail. aspx?filename=YTLX202107018&dbname=CJFD&dbcode=CJFQhttps: //oversea.cnki.net/KCMS/detail/detail.aspx?filename=YTLX202107018&dbname=CJFD&dbcode=CJFQ
[27] OH W T,VANAPALLI S K,PUPPALA A J. Semi-empirical model for the prediction of modulus of elasticity for unsaturated soils[J]. Canadian Geotechnical Journal,2009,46(8):903 − 914. DOI: 10.1139/t09-030
[28] ADEM H H, VANAPALLI S K. A simple method for prediction of the modulus of elasticity of unsaturated expansive soils. In: Khalili N, Russell A, Khoshghalb A, editors. Proceedings of the 6th international conference on unsaturated soils, Sydney, Australia. Unsaturated soils: research & applications. CRC Press; 2014.
[29] 唐栋,李典庆,周创兵,等. 考虑前期降雨过程的边坡稳定性分析[J]. 岩土力学,2013,34(11):3239 − 3248. [TANG Dong,LI Dianqing,ZHOU Chuangbing,et al. Slope stability analysis considering antecedent rainfall process[J]. Rock and Soil Mechanics,2013,34(11):3239 − 3248. (in Chinese with English abstract) TANG Dong, LI Dianqing, ZHOU Chuangbing, et al. Slope stability analysis considering antecedent rainfall process[J]. Rock and Soil Mechanics, 2013, 34（11）: 3239-3248. （in Chinese with English abstract）]https: //kns. cnki. net/KCMS/detail/detail. aspx?filename=YTLX201311032&dbname=CJFD&dbcode=CJFQhttps: //oversea.cnki.net/KCMS/detail/detail.aspx?filename=YTLX201311032&dbname=CJFD&dbcode=CJFQ
[30] LI J L,ZHANG L. Soil-water characteristic curve and permeability function for unsaturated cracked soil[J]. Canadian Geotechnical Journal,2011,48(7):1010 − 1031. DOI: 10.1139/t11-027
-
期刊类型引用(5)
1. 魏平新，郑志文，周志华，李秀娟，廖忠浈，刘任鸿. 广东省暴雨型浅层滑坡灾害动力预警模型与气象风险预警研究. 中国地质灾害与防治学报. 2024(02): 30-39 . 
本站查看
2. 林诗哲，胡新丽，张海燕，李宁杰，刘欣宇. 白鹤滩水电站库区红层砂岩干湿循环下结构劣化及渗透性演化规律. 中国地质灾害与防治学报. 2024(05): 67-77 . 本站查看
3. 陆烨，黄越川. 基于DIC分析的含石量对碎石土边坡稳定性影响. 中国地质灾害与防治学报. 2023(04): 49-57 . 本站查看
4. 柴少波，宋浪，周炜，付晓东，周永强. 水岩作用对充填节理岩石劣化的影响. 吉林大学学报(地球科学版). 2023(05): 1510-1520 . 百度学术
5. 王康，畅俊斌，李晓科，朱文峰，卢啸，刘慧. 基于水文过程和应力应变耦合的陕北黄土滑坡复活机理分析——以延安二庄科滑坡为例. 中国地质灾害与防治学报. 2023(06): 47-56 . 本站查看
其他类型引用(2)
计量
- 文章访问数: 344
- HTML全文浏览量: 263
- PDF下载量: 233
- 被引次数: 7
邮件订阅
RSS