植被对台风暴雨型滑坡发育的促进作用

孙 强, 张泰丽, 伍剑波, 韩 帅, 赵 阳

(中国地质调查局南京地质调查中心，江苏 南京 210016)

0 引言

我国是全球遭受滑坡灾害影响最为严重的国家之一。据全国地质灾害通报统计，2015—2021年间,我国发生地质灾害34 602起，因灾死亡1 337人,直接经济损失134.9亿元，其中滑坡是最主要的地质灾害类型，约占地质灾害总数的70%[1]。

在台风及期间强降雨的耦合作用下，我国的东南沿海地区常发生大规模的山体滑坡，直接或间接地给群众生命和经济财产造成重大损失。例如： 1996年7月31日至8月1日，台风“贺伯”在台湾省触发滑坡1 315 处、泥石流 20 多处，致使 73 人死亡、463 人受伤，财产损失达 10 亿美元[2]； 2016年9月27日至28日，台风“鲇鱼”在浙江省和福建省诱发大量的山体滑坡、崩塌和泥石流，造成32人死亡、3人失踪[3]。近年来台风暴雨诱发滑坡的数量呈增加趋势[4-5]。

台风暴雨型滑坡的影响因素除了台风风向、风速、风强变频、暴雨强度、降雨时长、总雨量等，还有山体坡度、坡向、岩组结构、构造、植被、高程等特征[6-7]。植被通常被认为对防止滑坡的发生具有积极作用，植被的覆盖率越低，滑坡越发育[8]。前人研究普遍认为植被对边坡稳定性的影响主要与根系有关，根系加固增加了浅层土壤的抗剪强度[9]。同时，植被冠层有利于降雨截留，减少雨水的入渗，增加蒸发，抑制土壤含水率的增加，降低土壤抗剪强度的衰减程度，有利于斜坡的稳定性[10-11]。然而，植被对滑坡发育的影响，目前存在较大的争议[12]。部分研究成果显示，许多滑坡上的植被根部没有生长至滑动面，因此植被反而增加了坡体重量[13]，乔木根-土间隙和根系腐烂后形成的大孔隙也会促进降雨入渗, 对滑坡产生不利影响[14-16]，同时台风的风荷载还会通过植被对滑坡的浅部结构产生不利影响[17]。但由于在台风暴雨作用下植被与土壤之间复杂的耦合关系，前人研究尚未明确台风暴雨条件下植被促进滑坡发育的具体模式及其影响程度。

为了深入分析植被与台风暴雨型滑坡的关系，本文研究了东南沿海飞云江流域由台风暴雨诱发的255处滑坡，分析了滑坡的物质结构、形态特征、植被覆盖类型等。在此基础上开展了风洞物理模拟实验，初步揭示了台风暴雨诱发滑坡过程中植被的促进效应。

1 研究区概况

飞云江流域位于我国浙江省东南部，发源于浙江省景宁县洞宫山白云尖，在瑞安市流入东海，流域面积约3 252 km2，域内最高海拔1 050 m，是该省八大水系之一(图1)。流域多年平均年降雨量1 850 mm，降雨集中分布在4—10月。流域大地构造单元属华南褶皱系，浙东南褶皱带泰顺—祖州断拗的中部，广泛分布着中生代火山岩，燕山期地质构造、火山喷发、岩浆侵入活动强烈，酸性、中酸性侵入岩较发育。区内断裂构造十分发育，影响流域地貌格局的大断裂有温州—镇海大断裂、泰顺—黄岩NE向大断裂及松阳—平阳NW向大断裂。流域属于亚热带常绿阔叶林植被带，主要植被群落有常绿阔叶林、针阔混交林、常绿落叶阔叶混交林、毛竹林、毛竹针叶混交林、毛竹常绿混交林和针叶林等[17]。

图1 飞云江流域地理位置

本文收集分析了飞云江流域台风暴雨诱发滑坡的调查资料。滑坡的规模整体较小，大多在1 000 m3以内，属于典型的小型浅层土质滑坡。发生在台风期间的滑坡结构具有明显的二元特征，即下部为基岩岩层，上部覆盖了一定厚度的松散层，厚度大部分小于5 m，处于1～2 m之间。滑体总体上是由基岩原地风化而成的含黏土、砂质黏土、粉土、亚砂土、碎石砂土等的第四系残坡积层或全风化层组成。滑坡多沿松散层与下伏基岩的接触界面或沿土层内部软弱面发生。滑坡面坡度集中在20°～50°之间，约占70%，其中在30°附近的分布最为集中(图2)。滑坡面形态主要为平直状。从时间上看，台风暴雨诱发的滑坡发生时间与台风的大风、强降雨过程大致同步，无明显的滞后现象。

图2 滑坡面规模与坡度特征

滑坡发生地坡面生长的植被类型多样，多为毛竹、灌木、乔木、草以及果树、农作物等的集合体，本研究将其划分为乔木、灌木、竹类及草地4种主要类型。据不完全统计，台风暴雨诱发的滑坡发生地坡面上的植被类型以乔木、竹类居多，草地、灌木次之(图3)，且坡面植被整体覆盖较好。

(a) 滑坡发生地植被分布野外照片 (b) 滑坡发生地坡面植被占比

2 技术方法

为分析验证台风条件下植被对滑坡发育的影响，本文进行了风荷载下植被对斜坡的结构损伤实验。实验平台采用中南大学风洞实验室。风洞洞体为全钢结构，实验段宽12 m、高3.5 m、长18 m，风速范围0～20 m/s，湍流度小于1%，并配备有高频电子压力扫描阀、微型压力传感器和动态数据采集分析系统等仪器设备。

2.1 相似条件

本研究根据台风暴雨型滑坡的规模、物质结构、气象条件和实验设备条件等因素，综合确定相似比和模型参数。模型的设计比例为1∶16，即相似系数n的取值为16。实验设计的三维物理模型中，主要参数的相似系数满足以下条件：

2.2 材料制作

坡体材料采用河砂、膨胀土、黏土、重晶石粉进行配比调制，力学参数参考区域内的经验值[19]。在实验室内进行若干组土工实验，直到材料各项参数接近经验值，最终确定河砂、膨胀土、黏土、重晶石粉的质量比例为39∶10∶39∶12。坡体模型材料各项物理力学参数见表1。

表1 坡体模型材料物理力学参数

2.3 模型植被

乔木高度高，风荷载作用更加明显，飞云江流域的乔木类型多为中乔，树高多在10 m左右，主根系深度一般在3 m以内，因此本次实验选择高度0.6 m，根系深度约0.15 m的植被模型。根据乔木的外形特征，制作了满足要求的模型树。模型树的树干和树叶为塑料材质，树根采用细钢丝制作，树叶尺寸约2 cm×3 cm，树冠的直径约30 cm； 树干的直径约2～4 cm，模型树根部两侧各有4条根系(共8条根系)，根系的展布范围约15 cm。

2.4 模型搭建

模型箱为钢制结构，长×宽×高为1.38 m×1.00 m×0.80 m。侧壁采用有机玻璃板，底部为建筑木板，设置模型斜坡坡度为35°。斜坡模型上布置了12棵模型树，从坡顶到坡脚共4排，每排3棵，编号见图4。模型树之间左右间隔0.4 m，前后间隔0.3 m。将植被模型固定在模型箱后，从下往上逐渐填入坡体材料，并逐层击实(图4)。完成后的植被模型覆盖度(植被的叶、茎、枝在地面的垂直投影面积占坡体总面积的百分比)约77%，与飞云江流域的植被覆盖度基本相符。

风洞实验室内不能产生任何颗粒物，因此，实验前用塑料薄膜覆盖滑坡模型表面，防止因风速过大及植被破坏导致颗粒物的产生。

图4 建成后的坡体模型

为监测不同风速作用下坡面和植被上的风压值，模型中共安装42处风压监测管。在滑坡模型的坡脚、坡面中部以及坡顶位置分别布置3处风压监测点(图5)； 模型中共有12棵模型树，每棵树的树干底部、中部和树冠处分别布置3处风压监测管，离坡面高度分别为0.1 m、0.2 m、0.35 m，以此来监测实验过程中植被不同高度受到的风压(图5)。

图5 模型风压监测点设置

2.5 实验工况

根据飞云江流域历次台风的风速统计资料，本次实验设置4种等级的台风工况，分别为8～9级(热带风暴)、12～13级(台风)、16级(超强台风)以及超过16级台风的极端工况。根据风速相似比确定实验风速分别为8.3 m/s、10.3 m/s、13.4 m/s和17 m/s。由于风洞实验室的限制，每种工况的实验时间为1 h。台风与坡体模型坡向的夹角设为0°。每种工况结束后，重新夯实土体，并将坡体植被重新恢复至原样，4种台风工况见表2。

表2 实验台风工况

3 实验结果分析

3.1 风荷载对斜坡推力的影响

受限于风洞实验室的安全要求，模型设置为迎风坡，即风力是作用于下滑力的反方向。但是通过分析风荷载的大小，同样能反映实际条件下顺风坡风荷载对斜坡下滑力的贡献。

按照传统的剩余推力法(不考虑降雨、地下水位及其他外部荷载)仅考虑滑体自重导致的下滑力、摩阻力以及滑面黏聚力的作用，在不同风荷载作用下，模型斜坡稳定性系数的计算公式为

。

(1)

式中：Fs为稳定性系数；WH为滑体自重，N；α为斜坡坡度，(°)；c为滑带黏聚力，Pa；φ为滑带内摩擦角，(°)；L为斜坡长度,m；fp为植被的风荷载(顺坡向为正，逆坡向为负)，N。

无风条件下，斜坡原始的抗滑力为1 394 N，下滑力为942 N，fp为0 N，斜坡稳定性系数为1.48。

当风速为8.3 m/s时，树冠的平均风压为50 Pa，树干中部平均风压为36 Pa，根部底部平均风压为32 Pa(图6)。将每棵树上不同高度处受到的风压值乘以其对应的宽度，可得到随高度变化的风荷载曲线。从植被的风荷载曲线看，在台风作用下，树冠处受到的风荷载远大于树干部位(图7)。

植被受风荷载作用于斜坡的推力为41 N，占抗滑力的2.9%。假设风向为顺风坡(不考虑风力角度的影响)，下滑力将增大至983 N，稳定系数下降至1.38。

图6 风速8.3m/s时11号树的风压变化

图7 风速8.3m/s时植被模型的风荷载变化

同理，当实验风速分别为10.3 m/s、13.4 m/s和17 m/s时，植被受风荷载作用于斜坡的推力分别增大至75 N、110 N和186 N，达到斜坡抗滑力的5.4%、7.9%和13.3%，稳定性系数降至1.34、1.24和1.17。

3.2 斜坡结构的破坏情况

在工况一、工况二和工况三条件下，随着风速的增大，虽然植被的摇摆程度逐渐增强，但实验结束后，坡体和植被未见明显的损伤特征。在工况四条件下，当风速达到17 m/s时，坡体结构发生了变形和破坏，主要表现为植被倾倒、土体隆起、裂缝产生。实验结束后对植被的倾倒角度进行了测量，斜坡的12棵模型树中倾倒的树木共7棵，模型的前两排树和最后一排的3号树均发生了明显的倾倒，其中第一排树的倾倒角度更大(约10°)，第二排树的倾倒角度稍小(约5°)，最后一排的3号树树冠处受到的风压较大，向后倾倒的角度为8°(图8、表3)。

图8 第4组实验后植被模型倾倒状况

表3 植被模型倾倒角度统计

部分植被模型下部的坡面发生隆起，隆起位置在第1排10号和11号树的根系附近，这两棵树的根土结构被破坏，隆起高度约2 cm，根系即将被拔起(图9)。这种现象与实际情况相符，野外调查过程中，在强大的风力作用下，高大的乔灌木往往被连根拔起[20]。

图9 10～12号树根部裂缝分布

坡体模型上第一排的10号、11号、12号树和最后一排的3号树的根部土体附近出现了裂缝，裂缝长约5～12 cm，宽0.2～0.5 cm(图10)。与斜坡失稳变形导致的拉裂缝、臌胀裂缝具有平行排列的特征不同，台风作用下植被撕拽导致的裂缝以树的根部为中心，主要分布在树的前侧，向四周呈放射状。

图10 3号树附近的裂缝分布

从上述变形破坏的现象中不难发现，斜坡的前排和最后排坡顶位置的植被根部土体最容易发生破坏。出现以上现象是由于前排植被在台风过程中首当其冲，受到的风荷载也较大，作用于根部土体的力矩更大，容易造成土体的拉裂。实验过程中发现与前3排植被相比，坡顶植被枝叶摆动剧烈，摇晃幅度更大，可能与位置较高，后部无其他物体阻挡而产生复杂的空气扰动有关。

3.3 斜坡材料渗透系数变化特征

在每种工况实验结束后，在坡顶2号树和3号树之间的位置取原状土样进行渗透系数的测试。圆筒的直径15 cm，高20 cm，每次取样深度为15 cm，获取的渗透系数能够较好地代表坡体材料的渗透能力(图11)。

与原始状况相比，在8.3 m/s、10.3 m/s、13.4 m/s的风速条件下，坡体材料的渗透系数无较大变化，在17 m/s以上的风速条件下，由于斜坡植被发生了明显倾斜，根部土体产生了大量裂缝，土体结构发生了明显破坏，渗透性急剧增大至原来的10倍以上。受限于实验条件，无法观察到植被风荷载造成了坡体结构性损伤，进而导致降雨入渗加剧的现象，但研究表明斜坡岩土体裂缝的增加和渗透系数的增大对斜坡的稳定性会产生显著的影响[21-23]。

图11 渗透系数变化曲线

本文通过风洞实验初步揭示了台风、植被耦合作用对斜坡的稳定性影响。台风暴雨对滑坡的诱发可以分为两个方面： 一方面台风、植被摇曳导致土体松动、裂缝产生和渗透性增加，进一步增加了降雨入渗软化土体的作用; 另一方面台风荷载通过植被作用影响斜坡附加的下滑推力，实验证明了在强烈的台风作用下，不能忽视台风-植被-斜坡的力学机制。

5 结论与现存问题

5.1 结论

(1)飞流江流域台风暴雨型滑坡基本上是第四系残坡积层的浅层土质滑坡，坡度多在30°左右，规模较小，发生时间与台风的大风、暴雨过程密切相关，无明显的滞后性。

(2)研究区台风暴雨诱发的滑坡坡面植被类型以乔木、竹类居多，草本、灌木次之，且坡面植被整体覆盖较好。

(3)在台风和植被的耦合作用下，风荷载最大可达滑坡抗滑力的13%，对滑坡稳定性有直接影响。当风向和滑向一致时，强台风条件下风荷载通过植被作用会导致斜坡稳定系数下降0.28，超强台风条件下会导致斜坡稳定性下降0.31。

(4)台风期间，植被摆动撕拽会导致植被根部产生大量裂缝，破坏根部土体的结构，造成土体物理力学指标的变化，渗透系数会增大至原来的10倍以上。

5.2 现存问题

本文对台风暴雨型滑坡的演化机制研究仅是初步的和尝试性的。植被对台风风荷载的传导过程是非常复杂的，如台风过程中存在脉冲风力、风振等情况。本次实验将风荷载简化为简单的静力作用，未对风荷载的力矩做进一步分析，而植被力矩是影响土体拉裂的重要因素。同时，受限于实验条件，实验中没有模拟台风伴随的降雨过程，对风、雨对滑坡的耦合作用机制缺乏讨论。因此，要阐明台风、暴雨、植被对滑坡的耦合作用，仍需进行深入的理论分析和大量的实验工作。