0 引言

玄武岩残积土是玄武岩经物理和化学风化后而残留在原地的残积物，西部交通工程建设过程中常遇到此类土。受西南地区湿热、多雨的季风气候影响，雨水的入渗和蒸发使其常处于含水率变动的非饱和状态，加之玄武岩自身节理、裂隙发育^[1]，为雨水入渗提供了有利条件，长期的降雨会使坡面出现坍塌、滑坡等灾害。

非饱和土强度特性及其边坡稳定性研究一直是岩土工程界的一个重要的课题，许多学者对此进行了研究。大量的统计资料表明：大多数的地质灾害是发生在降雨期间或降雨之后，且呈现出随降雨量增多而增强的规律^[2]。陈南等^[3]分析了红黏土浅层滑坡的破坏机理，提出了合理的评定性方法。陈东霞等^[4]探究了非饱和花岗岩残积土抗剪强度与基质吸力、含水率间的关系，认为基质吸力对含水率的变化敏感性很强，且基质吸力对抗剪强度有很大的贡献，许四法等^[5]在非饱和粉质黏土的研究中也得到了类似的结论。员康峰等^[6]通过非饱和直剪试验探究了膨胀土黏聚力、内摩擦角与基质吸力的关系，并结合有限元分析得出基质吸力对边坡稳定性有较大影响的结论。刘顺青等^[7]探究了粉状和砾状煤系土的水敏性，并结合强度折减法分析了含水率对边坡稳定性的影响规律。张社荣等^[8]结合强度折减法计算分析了瞬态渗流作用下的边坡稳定性，认为低强度长持时的降雨会增加深层土体的孔隙水压力，易使边坡产生深层滑动失稳。李哲等^[9]、楼晓昱等^[10]、马士国等^[11]认为雨水入渗会使边坡表面及其内部孔隙水压力增大，基质吸力逐步丧失而使边坡土体最终达到饱和状态，导致边坡安全系数降低。

然而，由于受土体的结构、土颗粒的矿物成分、应力历史状态等条件的影响，非饱和土的抗剪强度随基质吸力变化的规律对于不同的土体具有较大的差异性^{[4, 5, 12]}，因而不同的土体中基质吸力对边坡稳定的影响也有所不同。目前，关于玄武岩残积土的研究主要集中在其天然状态和饱和状态的物理特性及工程运用上，王参松等^[13]、徐良等^[14]探究了玄武岩残积土的工程及力学特性，张瑞敏等^[1]、余沛等^[15]探究了玄武岩残积土的胀缩性，唐军等^[16]探究了玄武岩残积土的崩解性，而利用非饱和土力学原理与方法来研究其强度特性及边坡稳定性却较少。因此，此次试验研究通过直接剪切试验探究了含水率与抗剪强度及其指标的关系，在此基础上运用强度折减法进行了边坡稳定性分析，并结合试验测试所得土水特征曲线分析了基质吸力对非饱和玄武岩残积土边坡稳定性的影响。

1 强度特性试验及分析 1.1 强度特性试验

(1) 土样基本物理性质

试验土样取自贵州省境内毕威高速公路玄武岩残积土坡，为黄褐色粉质黏土，其物理力学性质如表 1所示。

(2) 直接剪切试验

试验采用重塑土样，干密度为1.29 g/cm³，根据设计含水率27.5%，30%，32.5%，35%，37.5%，40%，42.5%和45%，制备试样并进行直接剪切试验。试验装置为FDJ-20型非饱和直剪仪，采用不固结快剪试验，剪切速率为0.8 mm/min，竖向荷载分别为50，100，200，400 kPa，剪切位移取7 mm。为了尽可能的减小剪切过程中水分的散失，采用不透水的塑料薄膜代替滤纸将试样与透水石隔开，并在剪切盒上围放含水率与试样含水率一致的纱布。剪切试验结果如表 2所示。

(3) 土水特征曲线试验

配制与直剪试验相同干密度、含水率的试样，采用滤纸法(国产双圈牌No.203型定量慢速滤纸)测定其基质吸力，其土水特征曲线如图 1所示。

1.2 试验结果分析

试验为不固结快剪试验，因此获得的抗剪强度为总应力强度，抗剪强度指标为表观黏聚力和表观内摩擦角，都隐含了基质吸力的影响^{[4, 12]}。

(1) 黏聚力与含水率的关系

黏聚力与含水率关系如图 2所示。由图 2可知，黏聚力随含水率的增大呈逐步减小的趋势且具有较强的阶段性。非饱和黄褐色玄武岩残积土的黏聚力主要由土颗粒间的静电引力、基于表面张力所形成的水膜连接力以及土体自身含有游离氧化铁所形成的特殊胶结作用力而组成。其中水膜连接力和胶结作用力与土体含水率密切相关。由于胶结作用力来自于土体自身，只有当土体中的自由水超过界限后，胶结物才会逐步溶解，因而在曲线第1段，随含水率的增大，主要是水膜连接力在减小，黏聚力下降的速度较小；在曲线第2段，随含水率增大，胶结作用开始丧失导致黏聚力减小的速度加快；在曲线的第3段，随含水率继续增大，土体接近饱和状态，水膜连接力基本为零，胶结物的溶解逐渐趋于稳定，此阶段随含水率增加，黏聚力下降的速度逐步减小且最终将趋近于零。

(2) 内摩擦角与含水率的关系

内摩擦角与含水率关系如图 3所示。由图 3可发现，内摩擦角随含水率的增大呈逐渐减小的趋势。在含水率小于32.5%的范围内，随含水率的增大，内摩擦角减小的速率较慢，可认为随含水率的增大，内摩擦角几乎不变；含水率超过32.5%时，随含水率增大，内摩擦角减小速度较快且基本呈线性减小趋势。

内摩擦角的大小反映了土体的摩擦强度，主要由滑动摩擦强度和咬合强度两部分组成。原状玄武岩残积土经风干重塑后，胶结物氧化铁从游离态转变为胶结态和结晶态共存，且此过程具有不可逆性。胶结态游离氧化铁在土体中形成粒团结构，结晶态氧化铁附着于土颗粒表面使其具有较大的粗糙度。在较低的含水率范围内，不同的含水率在相邻的土颗粒间形成的公共水化膜厚度较薄，土体中的水对土颗粒间的润滑作用微弱，故而随含水率的增加，内摩擦角下降的幅度很小。当含水率超过某一界限值后，随含水率的增大，公共水化膜厚度增大，土颗粒间的自由水增多，对土颗粒间润滑作用增大，导致滑动摩擦强度减小；同时随含水率的增大，游离氧化铁的胶结作用减弱，粒团较小的部分被包裹于公共水化膜中，导致土颗粒间的咬合强度减小。因而，在较高含水率范围内，内摩擦角随含水率增大而显著减小。

(3) 含水率对抗剪强度的影响

图 4为抗剪强度与含水率关系曲线，由图 4可知，不同竖向荷载作用下，抗剪强度随含水率变化特征基本一致，呈先不变后大幅减小的趋势。且竖向荷载越大，抗剪强度随含水率增加下降的幅度越大。竖向荷载为50 kPa，含水率从27.5%增加到45%时，抗剪强度由67 kPa下降到18 kPa，降幅为73%；竖向荷载为400 kPa，含水率从27.5%增加到45%时，抗剪强度由254 kPa下降到55 kPa，降幅为78%。对于边坡土体而言，当雨水入渗使深层土体饱和时，会导致抗剪强度丧失过大而诱发边坡失稳滑塌。

(4) 基质吸力对抗剪强度的影响

不同竖向应力下，非饱和玄武岩残积土的抗剪强度与基质吸力的关系曲线如图 5所示。由图 5可知，在各级竖向荷载作用下，基质吸力对抗剪强度的影响大致可以分为3段。基质吸力在0~210 kPa的范围内(阶段Ⅰ)，抗剪强度随基质吸力的增大显著增加；在210~526 kPa的范围内(阶段Ⅱ)，随基质吸力的增大，基质吸力对抗剪强度的贡献相比阶段Ⅰ明显减小；当基质吸力大于526 kPa后(阶段Ⅲ)，抗剪强度随基质吸力的增大几乎不变，这与各类非饱和土中对抗剪强度有影响的基质吸力范围^[17]通常为0~500 kPa的结论相一致。由此说明，试验所获得的抗剪强度隐含了基质吸力的影响，基质吸力在一定范围内对非饱和玄武岩残积土的抗剪强度有明显的提高作用，这种提高作用会随基质吸力的增大而有所下降，即基质吸力对抗剪强度的贡献是有限的。反之，当土体中的含水率增大时，非饱和土体中的基质吸力减小会导致抗剪强度降低。

2 边坡稳定性分析 2.1 工程概况

毕威高速公路(毕节至威宁)沿线穿越玄武岩分布地带，地形地貌复杂，地貌类型变化大，高陡边坡多，坡面上常附有厚度较大的坡积物，导致施工过程中地质灾害频发，严重威胁公路的建设安全。因此，开展玄武岩残积土的边坡稳定分析具有重要的实际价值。

2.2 计算模型建立

以毕威高速公路某右侧滑坡为研究对象，该段玄武岩残积土滑坡的自然坡度约为36°，坡高16 m，线路从斜坡下部通过，如图 6所示。边坡上覆残、坡积耕植土厚小于0.5 m，下部为全风化黄褐色玄武岩残积土，厚度大于24 m，为均质土坡，以此滑坡为原型建立简化后的理想模型，其基本力学参数如表 1所示。郑颖人等^[18]的研究表明弹性模量与泊松比对边坡安全系数的计算结果基本无影响，故可将各含水率下的弹性模量与泊松比统一取为天然含水率状态下的数值。

约束模型左右边界的水平位移，底部边界的水平和竖向位移。计算模型中的土体采用Mohr-Coulomb破坏准则，计算时采用关联流动法则。采用强度折减法对各含水率下的黏聚力和内摩擦角进行折减，求取各含水率下边坡的安全系数。

2.3 边坡安全系数的确定及结果分析

强度折减法其基本实质是材料的c和φ逐渐降低，导致某些单元的应力无法与强度相匹配，或超出了屈服面，不能承受的应力将逐渐转移到周围土体单元中去，当出现连续滑动面(屈服点连成贯通面)后，土体就将失稳。利用Abaqus提供的combine函数绘制场变量FV1与X方向的位移U1关系曲线，如图 7所示，以位移突变点所对应的折减系数为此含水率状态下边坡的安全系数。由图 7可知，对于含水率为42.5%的边坡，当折减系数小于0.61时，随着折减系数增大，边坡水平向的位移值几乎不变；但当折减系数超过0.61时，随折减系数增大，边坡土体水平向位移陡然增大，此时边坡失稳。

边坡失稳后，边坡内部会产生明显的塑性区并逐步形成贯通区，贯通区内外的土体变形相差十分明显，此后位移快速增加，失稳的边坡土体将沿着贯通区滑出，边坡失稳时的塑性应变云图及滑动面云图分别如图 8、图 9所示。

在此模型的基础上，通过强度折减法对含水率分别为27.5%，30%，32.5%，35%，37.5%，40%和45%的边坡进行计算分析，求取各含水率下边坡的安全系数，分析含水率对边坡稳定的影响，含水率与安全系数关系如图 10所示。由图 10可知，随着含水率的增大，非饱和黄褐色玄武岩残积土的边坡安全系数逐渐减小。在较低的含水率范围内(27.5%~30%)，边坡安全系数基本不变，含水率超过30%后，随含水率增大，边坡安全系数几乎成直线下降，当边坡土体含水率大于39.1%时，计算所得的安全系数均小于1，而实际上边坡已发生失稳滑塌。

Fredlund非饱和剪切强度理论表明^[19]，基质吸力是非饱和土抗剪强度的重要组成部分，边坡土体的安全系数与土体的抗剪强度成正比关系，且此次试验为不固结快剪试验，所得的黏聚力和内摩擦角均包含了基质吸力的影响，因此可用边坡的安全系数来间接的反应基质吸力对边坡稳定性的影响。结合滤纸法所得的土水特征曲线(含水率与基质吸力关系曲线)与强度折减法所得各含水率下边坡的安全系数绘制如图 11所示的安全系数与基质吸力的关系图，由图 11可知，基质吸力越大，安全系数越高，且基质吸力对边坡安全系数的影响大致可分为3段。在此次试验设置的条件下，对安全系数与基质吸力关系曲线进行分段线性拟合分析后可初步认为，当基质吸力大于526 kPa时(此值大小与干密度、温度有关，在此次试验设置的条件下为526 kPa)，基质吸力对边坡安全系数基本无影响，此段斜率k₃几乎为0；在小于526 kPa的基质吸力范围内，根据基质吸力对安全系数影响的大小可将基质吸力与安全系数的曲线分为两直线段(k₁>k₂)，这与基质吸力对抗剪强度的影响一致，因此可认为边坡土体在接近饱和状态时，基质吸力丧失的最快，且此时基质吸力的丧失对边坡稳定性的影响最大。

3 讨论

玄武岩残积土具有显著的水敏性，边坡土体含水率的增大会使土体抗剪强度大幅度下降而导致边坡失稳滑塌。毕威高速公路沿线穿越大量的玄武岩残积土质边坡，为保证公路能安全稳定的运行，需采取相应的措施降低边坡土体的雨水入渗量^{[3, 11]}。结合此次试验研究结论及毕威高速公路沿线边坡的特点，给出了以下建议：(1) 修建并完善边坡土体的排水系统，及时疏导边坡雨水。(2) 加强边坡的勘察检测，及时封堵坡体中的裂隙，减小边坡土体的雨水入渗量。(3) 对于裸露的边坡土体，应及时采取挂网喷射混凝土或植物生态护坡等措施，减小雨水冲刷面积。(4) 对于高陡边坡，应及时刷方减重并修建挡土墙以提高边坡的稳定性。

4 结论

(1) 含水率对非饱和玄武岩残积土的抗剪强度及其指标具有显著影响。随含水率的增大，抗剪强度先基本不变后大幅减小，黏聚力呈先小幅度减小后急剧下降的特点，内摩擦角呈先不变后线性减小的特点。

(2) 在526 kPa的基质吸力范围内，基质吸力对抗剪强度的影响呈阶段性减小；当基质吸力超过526 kPa后，基质吸力对抗剪强度基本无影响。

(3) 边坡土体含水率在30%以内，边坡稳定安全系数与含水率关系不大。含水率超过30%以后，随含水率增大，边坡稳定安全系数基本呈线性减小且减小的幅度较大，在边坡土体接近饱和时，安全系数最低，此时边坡最易发生失稳滑塌。

(4) 在约526 kPa基质吸力范围以内，基质吸力的变化对非饱和玄武岩残积土边坡稳性定有较大影响，随基质吸力降低，安全系数阶段性减小。可认为含水率增大，导致边坡土体的基质吸力降低，是边坡安全系数降低的最主要原因之一，玄武岩残积土自身的弱膨胀性、崩解性等对边坡稳定性的影响还有待进一步探究。