西南山区云南松根土复合体力学特性及其对浅层坡体稳定性的影响

2022-02-19李怡帆张国涛雷鸣宇方海燕李朝月

水土保持通报 2022年6期

李怡帆张国涛雷鸣宇方海燕李朝月

(1.中国科学院地理科学与资源研究所陆地表层格局与模拟重点实验室,北京100101;2.中国科学院大学资源与环境学院,北京100049;3.中国科学院水利部成都山地灾害与环境研究所山地灾害与地表过程重点实验室,四川成都610041)

浅层滑坡是在地壳运动、地形地貌、土壤、气候和植被覆盖等多因素作用下产生的,其中地震活动是滑坡发生的主要原因之一。2008年汶川地震引发了中国西南四川省数千起的山体滑坡[1-3]。开挖建设活动、森林砍伐等人类干扰也将降低浅层边坡土体稳定性[4]。在未来气候变化背景下,极端降水、融雪入渗等可增加边坡失稳的风险和山体滑坡等灾害发生的次数[3,5-6],严重威胁人类的生命安全。

植被作为一种可以进行人工配置和干预的生物工程,在提高土壤抗剪强度、稳定边坡方面的作用已经得到广泛证实[7]。Li等[8]研究得出,植被覆盖是滑坡发生的控制因素,相对于植被覆盖茂密边坡,在植被覆盖稀疏边坡发生滑坡的概率更高。植被对边坡稳定性主要从机械和水文两方面起作用,根系的力学加筋和锚固作用可稳定边坡[9],同时也可通过增加基质吸力、降低土壤含水量、拦截降水等[7,10-11]减小浅层滑坡发生的概率。

由于根系与土壤的变形模量存在巨大差异,因此受到剪切应力时,它们在变形过程中存在相互错动。在这个过程中,根系本身产生的拉应力、根土界面的摩擦阻力以及根系之间的土体侧向约束力发挥作用,从而提高根土复合体的抗剪强度[12]。自然界中,植被类型多样、根系分布复杂,不同类型植被的生理特征不同,根系力学性能存在差异,其对边坡稳定性的影响也不同[13]。根系稳定边坡的作用,与其形态、埋深、结构、生长状态等密切相关[14-16]。不同土层深度的根系土,根系含量、根系分布方式具有明显差别。因此本试验设计不同的根系密度与根系分布方式,试图探究不同深度根系土的抗剪强度差异。同时,滑坡的发生与降雨密不可分,这是由于降雨过程导致土壤含水量变化从而引起边坡失稳。关于土壤含水量与土体抗剪强度之间关系的研究已取得较多成果,然而,综合考虑根系密度、根系分布方式与土壤含水量对根土复合体抗剪强度的影响还需要进行深入的探究。

抗剪强度反映了土体抵抗外力剪切的能力,是衡量土壤力学特性的重要指标之一。在浅层滑坡发生过程中,普遍存在土壤力学特性改变的过程,以抗剪强度为指标探究根土复合体力学特性是研究根系固土、浅层滑坡发生机理的重要手段,对于滑坡预警具有重要指示作用。对于根系固土力学特性的研究,国内外已取得很多成果,但多集中于草本或灌木细小根系的研究,限于试验条件,针对乔木根系固土力学机制、稳定浅层边坡的研究较少,理论滞后于实践[14,17-18]。四川省凉山州全州森林覆盖面积约为2.60×106hm2,主要优势树种云南松、冷杉、丝栗、栎类与华山松分布广泛,其中云南松1.18×106hm2,占该州森林总面积的45.3%。相关研究表明[19],松属通过其侧根牵引对提高根际土层的强度、控制浅层滑坡具有重要意义。四川省地质灾害频发,而针对该地区乔木根系固土的理论探究较为不足,在中国西南四川省山区地质构造活跃、地形急变、土层较薄、降雨分布不均等自然条件下,开展乔木根系对于稳固边坡的研究具有重要意义。

因此,本文选取四川省凉山州西南山区云南松根土复合体开展研究,通过室内试验与数值模拟相结合的方式,探究不同根系密度、根系分布方式与土壤含水量下乔木根土复合体力学特性及浅层边坡稳定性,为西南山区浅层滑坡机理认识、识别预警与防治措施等提供理论基础与数据支持,对揭示大梯度环境植被—水—土耦合下山地灾害防治具有积极意义。

1 材料与方法

1.1 研究区样地概况

取样位置位于四川省西南部、金沙江下游的凉山彝族自治州热水河流域内一植被覆盖的斜坡,地处北纬26°03′—29°18′,东经100°03′—103°52′,位于青藏高原与四川盆地中间过渡地带[20]。该州地势西北高、东南低,地质地貌复杂,多山地、高原,州内自然灾害种类多、发生频繁、受灾范围广[21]。凉山州属亚热带季风气候区,年均温10.1～19.3℃,年降雨量776～1 170 mm[22]。冬季干冷,降雨稀少,夏季暖湿,集中了全年90%以上的降雨量,形成干湿季节分明、降雨集中的特点,从而极易发生山洪、滑坡、泥石流等自然灾害[23-24]。州内土壤类型多样,土层较薄,多数厚度不足100 cm[20]。取样地位于大凉山云南松(Pinusyunnanensis)纯林内,树高12.7±1.2 m,林地密度508±77株/hm2,林龄约20 a,林下为自然生长的灌草植物,坡面坡度约35°,土壤类型为壤土。

1.2 采样及制样方法

2019年9月16日,在采样地内斜坡区域选定灌草植被覆盖较少的1 m×1 m取样区域,清除表面植被和表层腐殖质,采集10—50 cm深度范围内土壤放入密封袋内带回实验室。实验室条件下,挑出根系、石块,按照《土工试验方法标准(GB/T50123-1999)》[25]采用室内试验采集土壤常规参数(表1)。将土样烘干至恒定重量,过2 mm筛后平铺在干燥的盘子中备用。

表1 土体常规物理与力学性质指标

采用人工整株挖掘的方法采集云南松竖向主直根系,取样时确保含根系土壤中根系的完整,避免人为破坏根系。随机选取生长正常、茎杆通直均匀、径级2 mm左右的根系,小心去除细根、须根及附着土壤后,洗净,自然风干,剪成长2 cm垂直根系,用密封袋密封,放置于实验室冰箱4℃环境中冷藏保存并尽快制样。

本试验方法严格按照国家规范《土工试验方法标准(GB/T50123-1999)》[25]制备根—土复合体土样,制样过程主要包括土壤含水量控制,压样,削土等。具体操作流程为:

(1)利用取样点土体的干密度、设计含水量(10%,15%,20%)及土样质量,计算所需加水量,将土样平铺于容器中,用喷水工具均匀加水,搅拌均匀使风干土与水充分接触,用保鲜膜封闭样品,静置24 h,使土样含水量均匀。加水量公式为[10]:

式中:mw为所需加水量(g);m为土壤样品质量(g);w0为风干含水量(%);w′为设计含水量(%)。

(2)将根系与土体分层放入直径61.8 mm,高20 mm的圆柱形试样盒中,分3层击实。根据试验设计,制备3种根系分布形式的复合材料(图1),根系垂直分布于材料中心(垂直分布),根系在材料中心水平分布(水平分布),根系交叉分布在样品中心(相交分布)。在每种根系分布方式下,制备3种根系密度的试验样品,即根系含量0.006,0.012和0.018 g/cm3。同时为了考虑土壤含水量对根土复合材料抗剪强度的影响,以5%为变化梯度,制备不同含水量10%,15%,20%3种条件下的试验样品。每种根系分布方式、每级根系密度、每种含水量各制备4个土样,共计108个根土复合土样。此外,利用原状土制备素土试验样品4个,以便比较根系存在产生的抗剪强度差异。

图1 根系排列方式(mm)

1.3 直剪试验

试验于2019年9月在中国科学院成都山地所土工实验室进行。采用应变控制直接剪切装置(图2),在4个竖向压力(50,100,150,200 k Pa)条件下对根土复合体进行直接剪切试验。根据试验结果绘制剪应力与剪切位移散点图,并根据库仑定律计算相应的黏聚力与内摩擦角:

图2 应变控制式直剪仪

式中:τf为土的抗剪强度(kPa);σ为法向应力(kPa);c为土的黏聚力(kPa);φ为土的内摩擦角(°)。

1.4 单根拉伸试验

选取根长约5～10 cm表皮完好的根系,用游标卡尺测定根样上、中、下3段处的直径,取其平均值作为该根样直径。使用电子万能试验机测量单根抗拉强度。试验时,只有当断裂面处于根段的中间位置时,被认为是由拉力引起的破坏而非其他测试过程中的损伤引起的。单根抗拉强度(F)的计算方法为:

式中:P是极限抗拉力(N);d是平均根径(mm)。

1.5 数值模拟

由于本试验主要分析浅表层土体的应力应变特性,研究根系对表层边坡土体稳定性的影响,因此把边坡底层基岩与表层根土材料区分开。根据浅表层的含义,将其深度设置为2 m,根据取样地坡面情况,设计边坡几何模型,坡高10 m,坡度30°,土体下部为基岩(图3)。在建模过程中,假设素土与含根系土的根土复合体为均质材料,且坡体的初始应力场由坡体自重产生。土体采用常规理想弹塑性体模型,即摩尔库伦单元进行模拟,根据文献[26]与土体物理性质试验,确定土体的材料属性(表2)。

表2 试验土体材料属性

图3 浅层边坡几何模型

根据根土材料直剪试验得到的黏聚力(c)与内摩擦(φ)采用FLAC 3D软件的强度折减法计算边坡安全系数。即通过折减系数F将土的黏聚力(c)与内摩擦(φ)逐渐降低,使某单元应力超过屈服强度,当多余的应力转移到周围单元而产生连续滑动面时,土体发生破坏,此时的F即为边坡的稳定安全系数,其计算公式为:

式中:cm,φm是土强度折减后的黏聚力(k Pa)和内摩擦角(°);Fr是强度折减系数。

1.6 数据处理与分析

数据统计分析采用Excel 2019。显著性检验采用无重复方差分析,借助SPSS 26.0操作完成。

2 结果与讨论

2.1 根系密度对复合体力学特性的影响

不同根系密度下复合体力学特性测试表明,土壤黏聚力与内摩擦角均具有显著的差异性(Pc=0.04＜0.05,Pφ=0.044＜0.05,详见表3—4)。根系分布方式与含水量一定的情况下,在各法向应力下根土复合体的抗剪强度随根系密度的增加呈现增加;相对于含水量10%与15%条件下,含水量20%时不同根系密度之间的抗剪强度差异明显减小。与素土相比,各条件下,黏聚力与内摩擦角均随根系密度的增加而增加,根土复合体黏聚力增量变化范围为1.05～6.05 kPa,内摩擦角增量变化范围为-1.0°～2.6°。根系垂直分布、土壤含水量10%时,根系密度由0.006 g/cm3增至0.018 g/cm3时,黏聚力分别增加了22.93%,28.95%和34.59%,内摩擦角分别增加2.95%,4.80%和5.90%,即在固定的根系分布方式与含水量条件下,随根系密度的增加,黏聚力与内摩擦角均呈增加趋势。在根系相交分布、土壤含水量20%时,黏聚力与内摩擦角随根系密度变化的规律表现不明显,可能是由于根系分布复杂、土壤含水量较高多种因素综合作用下的表现。相比于素土,在根系相交分布、含水量20%条件下黏聚力增量最高,达45.49%;在根系相交分布,含水量10%条件下内摩擦角增量最高,达到9.59%。

表3 不同条件下试样对应法向应力下土体抗剪强度值

分析表明,根系的存在可以提高根土复合体的抗剪强度,并随根系密度的增加而增加,根系增加抗剪强度的效应可能受土壤含水量的影响。这是由于根系的加入使得土体周围的容重增大,增强了土体的摩擦力。同时根系的存在改变了土壤颗粒之间的接触条件,在含有大量根系的复合体中,根系与土体之间产生更多的接触面,进而产生更强的根—土界面间的相互作用力,在进行剪切时,这种作用力转化为根系的抗拉力,因而提高了复合体的抗剪强度[27]。不仅乔木根土复合体具有这种性质,由草本[28]和灌木等[29]植被根系组成的根土复合体同样如此。

根系样本的根径范围为0.5～9.4 mm,抗拉强度范围为6.50～169.85 MPa,根系抗拉强度与根径之间具有明显的幂函数负相关关系(图4),其抗拉强度与根系直径拟合方程为:

图4 根系抗拉强度和根径的关系图及拟合曲线

随根径增大,抗拉强度减小,当根径小于2 mm时,抗拉强度递减速率最快,在2～10 mm根径范围内,根系抗拉强度随根径变化平缓。这说明细根相比粗根具有更强的抗拉强度,在土体中的加筋作用更为明显,其将土体的剪应力转化为自身的拉应力从而增强土体抗剪强度的作用更为显著。众多研究也得出了类似的结论,李牧阳[30]对油松、华山松、日本落叶松和锐齿槲栎的研究证实了根系抗拉强度随根径增加呈负指数减少;李佳[31]对陇南地区4种灌木研究也得出了类似的结论。

2.2 根系分布方式对复合体力学特性的影响

在固定的根系密度与土壤含水量条件下,不同根系分布方式下根土复合体的力学特性有明显的规律。土壤含水量10%时,随根系密度的增加,根系垂直分布方式下,黏聚力分别为16.35,17.15和17.9 kPa,内摩擦角分别为27.9°,28.4°和28.7°,根系水平分布方式下,黏聚力分别为15.65,16.8和18.7 k Pa,内摩擦角分别为27.5°,28.1°和28.6°(表4)。根系相交分布方式下,黏聚力分别为17.7,17.95和18.9,内摩擦角分别为28.8°,28.6°和29.7°。由此可知,根系垂直分布方式与水平分布方式对增强力学特性的效果相差不大,根系相交分布方式下,黏聚力与内摩擦角均高于同条件下垂直与水平分布的数值。在抗剪强度上,同样具有这样的规律;且相比于含水量10%与15%条件下,含水量20%时,根系相交分布下的抗剪强度与另外两种分布方式下抗剪强度的差异要明显。

表4 不同条件下试样黏聚力与内摩擦角

研究表明,根系相交分布方式下,抗剪强度、黏聚力与内摩擦角均明显高于同条件下根系水平分布与垂直分布,相较于根系单一分布方式,复杂根系分布方式对提高根土复合体抗剪强度的作用更显著。这主要是相交分布方式下垂直方向的根系起到锚固作用,水平方向的根系起到牵拉作用[32],在剪切过程中,土粒之间的摩擦力转化为根系的剪切与拉伸,根系的存在多角度上增强了复合体的抗剪强度,使复合体的理化性质发生显著的变化,即根系纵横交错状态下固土护坡效果最为显著。已有研究[14]表明,在水平根、垂直根和复合根3种形态的根型中,复合根最适合于增强土体强度、保持土体稳定。

2.3 土壤含水量对复合体力学特性的影响

在根系垂直分布、根系密度0.018 g/cm3时,根土复合体的黏聚力与内摩擦角随含水量的增加而减小,如土壤含水量从10%增至20%,黏聚力减少10.91%,内摩擦角减少4.79%。根系水平分布、根系密度0.018 g/cm3条件下,根土复合体的黏聚力与内摩擦角也存在相似的变化规律。然而,根系相交分布方式下,黏聚力和内摩擦角不随含水量的增加而减小,而是当含水量由10%增至15%时,黏聚力和内摩擦角减小,含水量由15%增至20%时,黏聚力和内摩擦角增加。由此可知,根系水平或垂直分布条件下,黏聚力与内摩擦角随含水量的增加而减小,根系相交分布条件下,黏聚力与内摩擦角随含水量的增加,均表现为先减小后增加的趋势。同样地,根系水平与垂直分布方式下,根土复合体的抗剪强度随含水量的增加而减小,在根系相交分布方式下,这种规律表现不明显。

土壤的抗剪强度主要取决于黏聚力和摩擦力,黏聚力与法向应力无关,主要由基质吸力或负孔隙水压力产生[33],受黏粒含量、含水量影响;摩擦力与法向应力有关,来源于剪切面土粒之间的摩擦和互相嵌入所产生的咬合力,受应力状态、土壤颗粒级配、土粒表面粗糙度等影响。本试验采用过2 mm筛的重塑土,故土壤细颗粒含量高,比表面积大,含水量提高使土壤颗粒之间的水膜增厚,土粒之间基质吸力减小,故黏聚力降低;当含水量继续增加时,水分产生的吸力对土壤颗粒的牵引力持续增强,因此表现为黏聚力增大。水分在土粒中表现为润滑作用,从而使摩擦力减小;高含水量可以增强团粒之间的咬合作用,有利于提高摩擦力。在根系相交分布方式下,复杂根系的加入改变了根土复合体物理力学性质,故在剪切过程中,表现为黏聚力与内摩擦角随含水量的增加,表现为先减小后增加的趋势。

灌木和草本根系与土壤之间的摩擦效应同样具有这样的规律。Li等[10]研究结果与本文一致,灌木根际土壤的黏聚力随土壤含水量的增加显著降低,内摩擦角降低但不显著。Wang等[34]根据根系拔出试验得出,牧草籼稻根系与土壤之间的界面摩擦效应随土壤水分含量的增加而迅速降低。但本研究并没有得到王晨沣等[33]研究得出的黏聚力随含水量增加先增加后减少的结论。这可能是由于本试验中各种性质根系的加入,使得不同含水量下根土复合体的抗剪强度变化不同于土壤。

2.4 边坡稳定分析

基于不同根土特性下数值模拟表明(表5),边坡安全系数模拟结果与复合体力学特性表现一致。根系存在下的边坡安全系数均高于素土条件下的边坡安全系数,同时边坡安全系数随着根系密度的增加而提高。因此,在进行护坡植被选择时,可将根系是否发达作为选择指标进行考虑。以往的研究多认为植被可以控制滑坡发生,然而有学者通过黄土高原降雨条件下植被覆盖斜坡上滑坡的失效机制研究,发现植被根系增加浅层土壤含水量而不能控制滑坡侵蚀[35]。同时,边坡植被增加了坡体自重载荷,进而增加了坡体失稳的风险。因而,在气候、地形、土壤等条件存在显著差异的地区,如西南山区、黄土高原和华北土石山区等,植被对浅层滑坡的影响存在差异。因此,如何从水—土—植被耦合的角度多层面定量表达植被对浅层滑坡的作用,还需要进一步的研究。

表5 不同根土特性条件下边坡安全系数对比

边坡安全系数在根系水平分布与垂直分布方式下相差不大,根系相交分布时,边坡安全系数明显高于同条件下其他两种分布方式,说明不同形态根系对土体稳定影响很大。因此,护坡植被配置时,纵向主根系与横向须根系发达的植被搭配护坡效果比单一植被下要好。选择不同根系形态的护坡植物,采用不同种类植物混交,形成植物根系纵横交错状态更有利于减少滑坡等山地灾害的发生[15]。

基于根系相交分布方式、最大和最小根系密度等多工况下的边坡x向位移云图,边坡安全系数随含水量增加的变化规律与土壤抗剪强度一致。在根系水平分布与垂直分布方式下,边坡安全系数随含水量的增加而减小。这也揭示了降雨降低边坡稳定性、诱发浅层滑坡的一个重要原因:降雨入渗会增加土壤载荷和孔隙压力[35-37],减弱土粒聚集性[7],减少浅坡土壤的基质吸力[11],降低根系抗拉强度[37]和土壤抗剪强度。在相交分布方式下,边坡安全系数随含水量增加表现出先减小后增加的趋势。在不同根系排列方式之下,边坡稳定性随含水量变化的不同表现可能是由于根系排列方式造成的根土复合体的物理力学性质不同造成的。

土壤水分是植物生长的重要基础,同时植被会对土壤水分分布和含量产生影响[38-41],因此单纯进行不同含水量下重塑土根土复合体的抗剪强度试验无法确定自然界中不同类型植被地下根系与土壤水分变化如何影响边坡稳定性[10],需要将野外原位观测试验与室内试验结合展开进一步研究。同时在本试验中,根土复合体是通过室内人工添加根系配置的,这样制成重塑土根土复合体与野外原状根土的理化性质存在一定差别,如没有考虑在根系生长过程中根系分泌物、根瘤菌、根菌等对土壤和根土复合体理化性质的影响,因此,需要通过野外采集原状土进行试验以得到更符合真实情况的结果。