含水率对安溪县花岗岩崩岗土体胀缩特性的影响*
2020-06-23林金石蒋芳市葛宏力连彬彬黎造宇黄炎和
章 智,陈 洁,林金石,蒋芳市,葛宏力,连彬彬,黎造宇,黄炎和
(福建农林大学资源与环境学院,福州 350002)
崩岗是发生在南方花岗岩红壤区,山坡土体在水力和重力复合作用下受到破坏失稳而发生崩塌的一种特殊土壤侵蚀形式和水土流失现象[1-3]。作为崩岗最重要的组成部分之一,崩壁是水力和重力复合作用最明显的部位,也是崩岗失稳崩塌的前提和崩积堆产生的物质来源[4]。崩岗发育的主要驱动力是降雨,土体中水分变化是诱发崩壁坍塌的主要因素。在降雨过程中,雨水通过侧向渗透进入崩壁,使得崩壁土体因含水率增大而发生膨胀;降雨结束后,由于水分快速蒸发,崩壁土体因含水率减小而发生收缩。水分快速进出崩壁土体,使其内部因干湿胀缩交替作用而产生裂隙并持续扩大,从而导致外层土体在水力和重力共同作用下不断崩塌剥落,加速崩岗发育和水土流失发展[5-6]。关于土壤胀缩变形问题,国内外学者已经开展了大量工作[7-10]。相关研究发现膨胀土的无荷膨胀率随初始含水率增加而降低[11-12],膨胀力则相反[11],其膨胀程度与时间[13]、上覆载荷[14]有关。随着含水率的降低,土体线性收缩率明显减小[12]。土体自身失水收缩会引起黏性土内部颗粒的变化,导致其土体强度降低[15],从而改变土体微观结构和内部颗粒排列[16-17],最终影响土体胀缩特性。土壤干湿胀缩过程同样会引起土壤容积和孔隙比的变化。土壤膨胀速度远快于收缩速度[18],且土壤容积变化速率小于含水率变化速率[12]。有研究结果表明对于同一种土体,初始干密度对其缩限几乎没有影响[19],含水率是主要影响因素。
目前,关于土壤胀缩特性的相关研究主要集中于干旱和半干旱地区的膨胀土,对于亚热带季风气候条件下的崩岗土(即花岗岩残积风化土)胀缩特性方面的研究较少。就南方花岗岩崩岗区而言,夏季降雨较多且气候炎热,水分频繁进出崩壁土体,在强烈的干湿胀缩交替作用之下引起崩壁崩塌,从而产生崩岗。在这一过程中,含水率变化是如何影响土体膨胀收缩尚不明晰,因此需要针对含水率对崩岗土体胀缩特性的影响进行研究。由于福建省安溪县龙门镇洋坑村为典型的南亚热带季风气候区,是全省崩岗数量最多的地区,是福建省乃至华南地区受崩岗侵蚀影响最严重的典型区之一。故本研究以福建省安溪县龙门镇洋坑村(118°03′E,24°57′N)的花岗岩崩岗土体为研究对象,以不同梯度的初始含水率为控制变量,开展室内无荷膨胀率和线性收缩率试验,分析无荷膨胀率和线性收缩率随时间变化规律,分别建立无荷膨胀率和线性收缩率与初始含水率之间的拟合关系,比较崩壁不同土层的土壤胀缩特性,进一步认识崩壁失稳发生崩塌的原因,对崩岗侵蚀机理的深入研究具有重要意义。
1 材料与方法
1.1 样品采集
图1 所示为研究区内的典型崩岗,包括崩壁、崩积体、冲积扇和集水坡面等主要组成部分。本研究于2017 年12 月3 日在典型崩岗崩壁上按照土层垂直分布,自下而上分碎屑层、砂土层和红土层进行采样。在同一崩壁上,不同土层的厚度、各采样点的具体位置和深度以及各土层的基本特征详见表 1。在各土层不同深度范围内多点混合取土,根据实际测量得到的各土层厚度,分别在各土层内纵向平均分为3 个具体采样深度范围,每个具体采样深度范围内采集3 个重复样,最终每个采样深度范围内的三个重复样混合为一个混合样,即为此采样深度的待测土样。各土层分别采集混合样3 个(红土层A1~A3,砂土层A4~A6,碎屑层A7~A9),共计 9 个土样,A1~A9 土层深度逐渐加深。将采集到的崩壁不同层次土样自然风干后,取一部分土样通过烘干法(105℃,8 h)测定风干土含水率,剩余土样用碾土棒对大颗粒进行碾压,通过 2 mm筛后装袋密封备用。供试土壤理化性质见表2。
1.2 试验设计
图1 研究区典型崩岗Fig. 1 Typical collapse in the study area
表1 采样层分布信息Table 1 Information of sampling layers
表2 供试花岗岩崩岗土体理化性质Table 2 Physical and chemical properties of the tested soil
本试验控制干密度为1.5 g·cm-3。根据表2 中崩壁不同层次土壤液塑限数据结果所示,各土层塑限极大值为36.85,出现在红土层中,极小值为27.98,出现在碎屑层中。因此初始含水率设置分别为红土层 15%、20%、25%、30%、35%五个梯度,砂土层和碎屑层初始含水率均为15%、20%、25%、30%四个梯度,膨胀和收缩试验各39 组处理,每组试验3 个重复,图中数据均以均值和标准差来体现。膨胀试验为有侧限饱和膨胀,历时24 h,数据采集时间步长为2 s;收缩试验为无侧限无荷收缩,历时96 h,数据采集时间步长为60 s。膨胀和收缩试验室内环境条件控制为室温不超过 30℃,相对湿度50%~60%。试样制备和试验步骤严格遵守《土工试验方法标准》(GB/T50123-1999)和《土工试验规程》(SL237-1999),同时也参照其他试验手册[22]。
无荷膨胀率试验装置和线性收缩率试验装置采用南京泰克奥科技有限公司生产的 TKA-PZY-1 型土壤膨胀仪和TKA-SSY-1 型土壤收缩仪。本试验采用南京泰克奥科技有限公司研发生产的TKA-DAI-8D 型数据采集盒和V2.1 型数据采集处理软件,以数据采集盒为中转设备,通过数据线将数显式位移传感器(测量精度为0.01 mm,测量误差为5%,最大量程为12.70 mm,工作温度为0~40℃)、数据采集盒和电脑端数据采集处理软件连接起来,准备工作完成后在电脑端软件页面点击“开始”即可自动采集试验数据。改变以往同类试验人工记录耗时费力,人工判读对数据精度和误差影响较大的弊端;同时还可以根据试验需要,设置数据采集时间步长,可获得精度较高的持续稳定数据,降低因人工判读而造成的试验误差,并且能够实时监测膨胀和收缩曲线变化。试验所得数据的有效测量精度和误差与后续计算分析精度均以数显式位移传感器的测量精度和误差为依据。
2 结果与讨论
2.1 不同水分条件下崩壁土体膨胀特性
为探究崩壁土体膨胀时程规律,对不同初始含水率条件下崩壁土体(红土层、砂土层和碎屑层)试样进行无荷膨胀率试验,相关试验结果可见无荷膨胀率随时间变化曲线,如图2 所示。从图中可知,崩壁土体无荷膨胀率随时间增加而增大,在24 h 左右吸水量达到饱和,膨胀趋于稳定。试样初始含水率越低,其无荷膨胀率越大,反之则无荷膨胀率越小。对于含水率在 15%~35%的五组不同初始含水率试样,各土层24 h 最终膨胀率总体趋势上均表现为红土层最大,砂土层次之,碎屑层最小。其中红土层在 15%初始含水率时的最终膨胀率最大,其值为11.64%;在35%初始含水率时的最终膨胀率最小,其值为1.13%;红土层的最大膨胀率与最小膨胀率之间相差10.30 倍。砂土层和碎屑层的最大膨胀率和最小膨胀率与红土层表现为同样的规律,均出现在初始含水率梯度的两个极值,但不同于红土层的是其最大膨胀率和最小膨胀率的差距分别为19.02 倍和14.80 倍。此外,由图2还可知,崩壁土体在不同初始含水率条件下的膨胀时程均有明显规律性,可将无荷膨胀率变化曲线分为三个明显的不同阶段,即快速膨胀阶段、缓慢膨胀阶段和膨胀稳定阶段,这与其他研究者的结论一致[13,23],说明崩壁土体的膨胀程度和延吉盆地的黑土一样,主要是由土体吸水增湿的初期所控制。初始含水率越小,土层深度越浅,膨胀时程曲线三阶段越明显,特别是快速膨胀阶段;随着初始含水率的增加和土层深度的加深,膨胀过程曲线渐趋平缓,三阶段膨胀变化程度明显渐渐降低。通过比较发现初始含水率越小,其快速膨胀历时越短,快速膨胀速度越快。
图2 不同初始含水率条件下无荷膨胀率随时间变化曲线Fig. 2 Temporal variation curve of no-loading expansion rate with initial water content
由图2 可知,总体而言大部分膨胀在快速膨胀这一阶段内完成,速度快且明显,大致在膨胀开始后的1 h 以内完成,约占最终膨胀量的75%~85%;缓慢膨胀阶段的膨胀量约占最终膨胀量的 10%~15%,大致在膨胀开始后的1~5 h 时间段;稳定阶段的膨胀量几乎不会发生较大波动,在膨胀曲线上表现为相对平稳,大致在膨胀开始后的5~24 h 时间段,其膨胀量不足最终膨胀量的5%。分析原因可知,崩壁不同层次重塑土壤内部因含水率较低而存在相对较多的内部孔隙。当试样遇水时,崩壁土壤吸力大、吸水较快,其膨胀速率也较快,在这一阶段完成绝大部分膨胀变形。不同于膨胀土矿物成分主要由蒙脱石等亲水性矿物组成,其膨胀变形主要发生在膨胀试验开始后的30 min 内[13],崩岗土壤的黏土矿物成分不含有蒙脱石[24],其膨胀变形主要发生在膨胀试验开始后的1 h 内。随着水分含量的增多,土体颗粒间距增大,孔隙中的空气逐渐被水分排挤出,试样吸水速率减缓,无荷膨胀率增幅趋于平缓直至达到稳定状态,崩岗土壤快速膨胀时间大于膨胀土,其最终膨胀率小于膨胀土。
不同土层初始含水率与无荷膨胀率的拟合关系曲线如图3 所示,随着初始含水率的增大,红土层、砂土层和碎屑层的无荷膨胀率均呈现出减小的趋势,反之,其无荷膨胀率均表现为增大的趋势。随着初始含水率的变化,红土层的无荷膨胀率变幅约为10%,砂土层和碎屑层的无荷膨胀率变幅均小于红土层,分别约为 7.50%和 6.50%。通过对无荷膨胀率和初始含水率的关系进行回归分析发现,不同于膨胀土的无荷膨胀率与初始含水率存在着良好的线性关系[11,13,23],初始含水率对崩壁不同土层无荷膨胀率的影响存在着明显的指数递减关系。这是因为膨胀土由于其矿物成分主要由蒙脱石和伊利石等亲水性黏土矿物组成,膨胀潜势远大于崩岗土壤,另一方面则是由于其他研究者的膨胀时间数倍于本试验的膨胀时间,所以在含水率和膨胀率的拟合关系上表现出差异性。红土层、砂土层和碎屑层的拟合方程的决定系数R2均大于0.96,其拟合所得回归方程均符合以下函数关系:
式中,δe为无荷膨胀率,%;ω为初始含水率,% ;a、b和c均为回归方程参数,与土壤初始干密度和土层性质有关。
图3 不同土层初始含水率与无荷膨胀率的拟合关系曲线Fig. 3 Fitting relationship between initial water content and no-loading expansion rate relative to soil layer
2.2 不同水分条件下崩壁土体收缩特性
通过对不同初始含水率条件下的崩壁土体进行线性收缩率试验,研究崩壁土体收缩时程规律,相关试验结果见线性收缩率随时间变化曲线,如图 4所示。从图中可知,崩壁土体的线性收缩率随时间的增加而增大,最终趋于某一稳定值。试样初始含水率越高,其线性收缩率越大,反之则线性收缩率越小。对于含水率在 15%~35%的五组不同初始含水率试样,各土层的96 h 最终线性收缩率整体趋势上均表现为红土层最大,砂土层次之,碎屑层最小。其中红土层在35%初始含水率时的最终线性收缩率最大,其值为 5.53%;在 15%初始含水率时的最终线性收缩率最小,其值为 0.23%;红土层的最大线性收缩率与最小线性收缩率之间相差 24.04 倍。砂土层和碎屑层的最大线性收缩率和最小线性收缩率与红土层表现为同样的规律,均出现在初始含水率范围的两个极值,但不同于红土层的是其最大线性收缩率和最小线性收缩率的差距分别为7 倍和2.80倍。此外,由图4 还可知,崩壁土体在不同初始含水率条件下的收缩时程同样具有明显的规律性,可将线性收缩率变化曲线分为三个明显不同的阶段,即快速收缩阶段、缓慢收缩阶段和收缩稳定阶段,这与其他研究者的结论一致[19,24],说明崩壁土体的收缩程度和膨胀土一样,主要是由土体失水减湿的初期所控制。初始含水率越大,土层深度越浅,收缩时程曲线三阶段越明显,特别是快速收缩阶段;随着初始含水率的减小和土层深度的加深,收缩时程曲线渐趋平缓,收缩变化三阶段明显程度逐渐降低。通过比较发现初始含水率越大,其快速收缩历时越短,快速收缩速度越快。
由图4 可知大部分收缩在快速收缩这一阶段内完成,速度快且明显,大致在收缩开始后5 h 内完成,约占最终收缩量的85%;缓慢收缩阶段的收缩量约占最终收缩量的 10%,大致对应图中 5~18 h时间段;稳定阶段的收缩量几乎不会发生明显波动,在收缩曲线上表现为平直线,大致对应图中18~96 h 时间段,其收缩量约占最终收缩量的5%。分析原因可知,初始阶段土体蒸发损失的是自由水,而随着土体中的水分不断蒸发散失,自由水全部损失后结合水也开始逐步散失直至损失殆尽,在这一过程中,水膜逐渐变薄,双电层厚度因此而减小;同时由于水分散失产生的负孔隙压力迫使土粒之间不断靠拢距离逐渐缩小,颗粒间的斥力渐大以至于超过负空隙压力,此时土体不再收缩并趋于稳定[23]。
图4 不同初始含水率条件下线性收缩率随时间变化曲线Fig. 4 Temporal variation curve of linear shrinkage rate with initial water content
不同土层初始含水率与线性收缩率的拟合关系曲线如图5 所示。随着初始含水率的增大,红土层、砂土层和碎屑层的线性收缩率均呈现出增大的趋势,反之,其线性收缩率均表现为减小的趋势。随着初始含水率的变化,红土层的线性收缩率变幅约为5%,砂土层和碎屑层的线性收缩率变幅均小于红土层,分别约为2%和0.65%。通过对线性收缩率和初始含水率的关系进行回归分析发现,初始含水率对崩壁不同土层线性收缩率的影响存在着明显的指数递增关系,与章李坚等[24]研究结果一致。红土层、砂土层和碎屑层的拟合方程的决定系数R2均大于0.96,其拟合所得回归方程均符合以下函数关系:
图5 不同土层初始含水率与线性收缩率的拟合关系曲线Fig. 5 Fitting relationship between initial water content and linear shrinkage rate relative to soil layer
式中,δsl为线性收缩率,%;ω为初始含水率,% ;a、b和c均为回归方程参数,与土壤初始干密度和土层性质有关。
2.3 不同水分条件下崩壁土体胀缩特性差异
从上述膨胀试验和收缩试验的结果分析可知,膨胀时程曲线和收缩时程曲线在整体上具有明显的相似性。首先,膨胀时程曲线和收缩时程曲线的总体变化趋势基本上一致,均表现为膨胀率和收缩率随着胀缩时间的增大而增大。其次,在胀缩过程特征上,膨胀和收缩均可分为明显相似的3 个阶段,均经历了从快速胀缩的斜直线至缓慢胀缩的外凸弧线最后至胀缩稳定的平直线的胀缩变化过程;但膨胀和收缩两个相反过程各对应阶段的历时及其胀缩变形程度具有较大差异,表现为膨胀速率远快于收缩速率,最终膨胀量远大于最终收缩量。
崩壁不同土层最终膨胀率和最终收缩率随初始含水率的变化曲线见图 6。从图中可知,随着初始含水率的增大,无荷膨胀率随之减小,线性收缩率随之增大,具体表现为红土层的最终膨胀率由11.64%减小至1.13%,砂土层的最终膨胀率由8.37%减小至 0.44%,碎屑层的最终膨胀率由 7.40%减小至 0.50%;红土层的最终收缩率由 0.23%增加至5.53%,砂土层的最终收缩率由0.33%增加至2.30%,碎屑层的最终收缩率由 0.35%增加至 0.98%。在初始含水率变化范围内,不同土层的最大无荷膨胀率均表现为红土层最大,砂土层次之,碎屑层最小的趋势,最大线性收缩率也呈现出一致的趋势。对同一土层而言,膨胀变化幅度均较收缩变化幅度大。其中红土层的膨胀变化幅度为10.51%,收缩变化幅度为 5.30%;砂土层的膨胀变化幅度为 7.93%,收缩变化幅度为 1.97%;碎屑层的膨胀变化幅度为6.90%,收缩变化幅度为0.63%。相同土层的膨胀变化幅度和收缩变化幅度同样均为红土层最大,砂土层次之,碎屑层最小;膨胀变幅远大于收缩变幅,红土层的膨胀变幅为其收缩变幅的1.98 倍,砂土层和碎屑层的膨胀变幅分别为其收缩变幅的 4.03 倍和10.95 倍。通过对比不同土层的胀缩性能发现,膨胀和收缩之间存在着发展不同步、胀缩变形不均衡,表现为红土层的膨胀和收缩变化幅度最大,分别较砂土层高 2.58%和3.33%,较碎屑层高3.61%和4.67%。说明随着土层深度的加深,表层红土的胀缩性能最大,砂土层次之,碎屑层最小,这与在野外实际观察到的现象是一致的。这种不可逆的胀缩现象可能是造成土体产生裂隙进而引起崩壁坍塌的原因。
图6 不同初始含水率条件下各土层最终胀缩率对比图Fig. 6 Comparison between soil layers in final expansion and shrinkage rate relative to initial water content
有研究认为,土壤胀缩变形是一对作用相反、能够互逆的过程。从前述膨胀时程曲线和收缩时程曲线中可知,虽然膨胀曲线和收缩曲线在整体上具有一定的相似性,但是其膨胀和收缩过程并不是一个完全可逆的过程。数学意义上膨胀过程的逆过程为收缩过程,即膨胀曲线与收缩曲线为中心对称曲线。从图6 中同理可知,红土层、砂土层和碎屑层所对应的随初始含水率变化的胀缩对比曲线均呈相互交叉状而非中心对称曲线,这从数学意义上说明膨胀与收缩是一个不可逆的过程。另一个方面,在初始含水率变化范围内,红土层的最大胀缩差为6.11%,最大膨胀率为最大收缩率的 2.10 倍;砂土层的最大胀缩差为6.07%,最大膨胀率为最大收缩率的3.64 倍;碎屑层的最大胀缩差为6.42%,最大膨胀率为最大收缩率的 7.55 倍。由此可说明,膨胀和收缩是不可逆的过程,同时也是不同步和不均衡的。
3 结 论
本试验通过室内无荷膨胀率和线性收缩率试验对安溪县崩岗区崩壁土体胀缩性能进行研究,总体上崩壁土体膨胀和收缩过程具有明显的相似性,均可大致划分为斜直线快速阶段-外凸弧线缓慢阶段-平直线稳定阶段三个不同过程。不同土层的无荷膨胀率均随初始含水率增大而减小,线性收缩率则相反。初始含水率与无荷膨胀率和线性收缩率之间均具有良好的指数关系,其中初始含水率与崩壁不同层次土壤无荷膨胀率之间存在着明显的指数递减关系。不同土层最大无荷膨胀率和最大线性收缩率均表现为红土层最大,砂土层次之,碎屑层最小,说明随着土层深度的增大,崩壁土体的胀缩性能逐渐降低。这种不可逆的胀缩现象是造成土体产生裂隙进而引起崩壁坍塌的原因。本试验以崩壁土体为对象,由于只考虑含水率变化对土壤胀缩特性的影响,仅从一个侧面揭示了崩岗土体的胀缩特性,因此后续还要在此基础上结合土体结构、胶体含量、矿物成分、上部载荷和干密度等内外影响因素进一步研究崩岗土体胀缩特性。