闫钇全, 刘 琦,2, 邓大鹏, 廖启迪
(1.同济大学 土木工程学院 地下建筑与工程系, 上海 200092;2.岩土及地下工程教育部重点实验室(同济大学), 上海 200092)
中国西南岩溶地区具有独特的地表、地下二元水文地质结构,因此其水土流失存在特殊的地表流失和地下漏失两方面,贵州省人口密集和不合理的社会经济活动的影响,呈现出大面积岩溶石漠化现象[1-3]。其中,贵州省石漠化面积约37 597 km2,占中国石漠化总面积的28.1%[4-5]。影响西南岩溶地区水土流失和石漠化发展程度的诸多因素包括气候、地理地貌、地质条件和人类活动等[6]。经多年研究成果可以看出,侵蚀性降雨是其中非常重要的影响因素之一。侵蚀性降雨是指能够造成真正意义土壤侵蚀的降雨[7],且短时(1-3小时)强降雨更易引发相关地质灾害发生[8]。杜波等[9]提出小雨量、短历时、大雨强的降雨是岩溶地区最易产流产沙的降雨类型,西南地区雨热同期,降雨时空分布不均,暴雨发生频率高、强度大,且大暴雨和特大暴雨所占比例高[10-11]。这些自然因素都加剧了岩溶地区的水土流失。
目前国内外学者对岩溶地区水土流失机理的研究已有一定成果,研究方法主要有径流小区观测法,核素示踪法、3S技术、室内模拟降雨试验法[12-16],由于野外场地与研究方法的限制,室内试验法针对进一步研究水土流失的机理具有一定优越性,目前已有相关研究集中在基岩裸露率、地下孔隙度、坡面地貌特征、降雨强度等影响因子。雨量与土壤侵蚀联系密切,雨强是一个十分重要的降雨参数,有学者研究指出[7,17],使用平均雨强对侵蚀性降雨进行定义,有一定的错选率,存在放大低强度降雨对土壤的侵蚀作用,而弱化了高强度降雨对土壤侵蚀的影响。在黄土地区存在总降雨量和土壤量相关性低的情况。目前还少有学者进行不同降雨强度交替作用下室内模拟降雨试验,因此本文基于贵州贞丰关岭—花江石漠化治理示范区的野外研究现状,通过模拟降雨强度交替组合方式探究自然界中降雨条件变化对表层岩溶坡面土壤地表流失与地下漏失的影响,该项研究对揭露岩溶地区水土流失的规律及防治石漠化具有理论和现实意义。
研究区位于中国贵州省关岭县与贞丰县相交处的北盘江峡谷花江地段,东经105°36′30″—105°46′30″、北纬25°39′13″—25°41′00″,示范区总面积51.62 km2。多年平均降雨量1 100 mm,降水量时空分布不均。研究区岩性以碳酸盐岩为主。示范区石漠化面积2 527.71 hm2,其中轻度石漠化面积1 497.18 hm2,占59.23%;中度—强度石漠化面积占比40.77%。研究区内中度—强度石漠化较为常见,局部地区已经严重到威胁当地居民的生产生活[18]。
有学者通过研究确定了坡地土壤侵蚀主要由大雨量、中历时、大雨强的降雨造成[19-20]。现对研究区所属贵州省安顺地区气象站1991—2020年连续30 a的逐日降雨数据进行分析,得出研究区共发生降雨5 423次。降雨量资料来源于中国气象数据网(http:∥data.cma.cn)。按气象学的划分方法(24小时降雨量)对研究区的降雨等级进行分级,即大暴雨(p≥100 mm)、暴雨(50≤p<100 mm)、大雨(25≤p<50 mm)、中雨(10≤p<25 mm)、小雨(p<10 mm)进行分级。将5 423次降雨按照降雨量大小统计分析,得出研究区降雨特征见图1。小雨共4 394次,占降雨总次数的81%,累计降雨量8 885 mm,占降雨总量的23%;中雨共609次,占降雨总次数的11%,累计降雨量9 710 mm,占降雨总量的25%;大雨共281次,占降雨总次数的5.2%,累计降雨量9 535 mm,占降雨总量的24.7%;暴雨次数120次,占降雨总次数的2.2%,累计降雨量8 107.3 m m,占降雨总量的21%;大暴雨次数19次,占降雨总次数的0.6%,累计降雨量2 392 mm,占降雨总量的6.3%。
图1 不同等级降雨总量及次数
通过对研究区30年的降雨特征分析得出:降雨等级按照发生频次排序为:小雨>中雨>大雨>暴雨>大暴雨;按照降雨量排序为:中雨>大雨>小雨>暴雨>大暴雨。暴雨及大暴雨主要集中在4—10月份,中雨、大雨、暴雨、大暴雨均呈现雨季降雨频次增加,旱季减少的情况。但这些降雨是否引起侵蚀性降雨还需进一步验证。
为了进一步了解研究区的侵蚀性降雨特征,分析侵蚀性降雨对土壤地表流失/地下漏失的影响。本课题组在研究区选取典型坡面修建径流小区和气象站以获取数据。该气象站每30 min记录一次降雨量数据并自动储存,共统计了2018年4月至2020年9月的降雨数据,以每次降雨后径流小区是否产沙作为发生侵蚀性降雨的依据。随机选择研究区2 h降雨历时的5组侵蚀性降雨量数据分析(图2),以该区2019年5月26日的降雨数据为例,其降雨强度瞬时最高可达34.8 mm/h,但仅0.5 h过后降雨强度变为零,之后又开始降雨。可知研究区的降雨具有瞬时雨强大、降雨在时间上较为集中且不同雨强交替发生的特点。
图2 降雨强度随降雨时长变化
试验土壤样品采自示范区中碳酸盐岩风化形成的红黏土,该区土壤一般都很浅薄,但质地黏重、透水性差,易在降雨季节发生水土流失,其物理性质见表1。进行试验前先对土壤进行筛分,去除石砾,配备成含水率25%的土壤均匀地铺在模型表面并压实,用以模拟示范区表层土。
表1 试验土样物理性质
试验设备为自制的表层岩溶裂隙带土壤地表流失和地下漏失模拟装置(图3)[21],试验设备主要由模型箱和降雨平台组成。模型箱中由切割成型的泡沫、承压铁板和承压铁板上的石灰岩石板和填充碎石组成,该部分用以模拟基岩,之后在基岩上铺设土壤和放置碎石充当土岩界面和裸露岩体。承压铁板上预留出一个长宽为10 cm×2 cm的开口,用以模拟野外裂隙。降雨模拟装置采用针头式降雨方式模拟降雨,降雨均匀系数大于80%,降雨器高度为2 m。降雨强度的变化范围10 mm/h~120 mm/h,调节变化时间<15 s。
依据熊康宁等[22]及本课题组在示范区典型岩溶坡面所建径流小区调查统计结果[23],结合示范区野外调查可知示范区典型坡面的倾角范围为10°~25°,基岩裸露率为60%~85%,将本试验的坡度选取为15°,基岩裸露率为60%,选取4个降雨强度(10, 20, 40, 80 mm/h),试验每0.5 h变化一次降雨强度,共变化4次,一组试验持续2 h。将选取的降雨强度换算为24 h降雨量,按照气象学中依靠降雨量区分降雨等级的方法区分本试验的降雨等级,故对本试验降雨强度的描述为:10 mm/h的降雨为中雨,20 mm/h的降雨为大雨,40 mm/h的降雨为暴雨,80 mm/h的降雨为大暴雨。降雨强度的设定根据所统计的示范区降雨特征及径流小区侵蚀性降雨特征选取,正式试验前先进行一组预试验标定降雨强度,将室内试验次降雨径流深情况与野外实际情况相比较,依据动力相似、产生径流量相似原则[24]确定。试验开始后,每5分钟用采样瓶分别接满土壤地表流失与地下漏失泥水样,用烘干法测得每组数据的土壤地表流失和地下漏失量。本试验的试验条件设置可见表2。
注:1:供水水龙头;2:稳压水箱;3:阀门;4:针孔式降雨器;5:支架;6:土壤层;7:地表碎石;8:石灰岩石板;9:采样瓶;10:水土地下漏失;11:水土地表流失;12:电子秤;13:岩溶裂隙;14:泡沫填充物。
降雨是发生土壤地表流失/地下漏失的重要驱动力,分析两种不同降雨强度交替作用下土壤地表流失、地下漏失的特征(图4),可知在此试验条件下土壤地下漏失占土壤侵蚀总量的63%~91%,平均占比为80%,6组试验均呈现出土壤地下漏失占土壤侵蚀主导地位的现象。组合1土壤地下漏失占土壤侵蚀总量的85%,土壤地下漏失总量约为其他降雨强度交替作用情况下的1.7~2.4倍,相比于其他降雨强度交替作用的情况更容易发生土壤地下漏失。组合6的土壤地表流失占土壤侵蚀总量的37%,土壤地表流失总量约为其他降雨强度交替作用情况下的1.4~4.4倍。在降雨量一定的情况下,组合6有最高的土壤地表流失量,而地下漏失量为同组最低。
表2 试验条件
组合1、组合3、组合6均设置连续降大雨1 h,这3个组合的土壤地表流失量分别为1.97 g,1.82 g,2.83 g,土壤地表流失量高于其它间断性降大雨的3组。结果表明在总降雨量相同的情况下,高降雨强度下的集中降雨更容易导致土壤地表流失的发生。土壤地表流失的发生常常伴随着坡面径流的形成,高降雨强度集中降雨雨水多转化为土壤水,更容易改变土体结构,快速形成坡面径流,发生坡面径流携沙流失的现象。组合2、组合4、组合5即间断性降大雨的情况中土壤地表流失量小于其它3种组合,土壤地下漏失量差异不大,结果表明在总降雨量相同情况中,间断性大雨的不同组合下土壤地下漏量对降雨交替作用的变化响应不敏感。
图4 两种降雨强度交替作用下土壤地表流失/地下漏失特征
分析多种降雨强度交替作用下土壤地表流失、地下漏失总量的特征(图5)可知,在这6种组合情况下,总降雨量相同,土壤地下漏失占土壤侵蚀总量的56%~92%,平均占比为82%,土壤地下漏失仍然为土壤侵蚀的主要形式。降雨量与土壤侵蚀呈正相关关系,按大暴雨、暴雨、大雨、小雨组合下的组合9,组合10具有同组最大的土壤地表流失/地下漏失总量,分别为22.75 g,28.07 g。在组合9和组合10中,组合10具有最高的土壤地表流失量和土壤地下漏失量。在降雨前期,岩溶坡面的土壤经历的短历时大暴雨会快速导致土壤水饱和,降雨强度大于土壤水入渗率,形成地表径流通道和地下漏失通道,之后降雨强度降低但已经形成的土壤侵蚀通道依然存在,降雨难以转化为土壤水,通过已形成的土壤侵蚀通道携带土壤侵蚀。
组合11和组合12模拟了间断性降中雨的情况,在这种情况下土壤地表流失、土壤地下漏失量均较小,两种组合土壤地表流失/地下漏失总量为6.93 g,5.46 g,土壤地表流失量较小,土壤地表流失量占比仅为8%,土壤地下漏失占比92%。土壤地下漏量对降雨强度交替作用的变化响应不敏感。间断性降中雨,降雨强度不能导致土壤水饱和,雨水入渗土壤较为充分,沿着已有的漏失通道,携带土壤漏失。由于雨水大部分入渗土壤,未能形成坡面径流,土壤地表流失量较小。
图5 不同降雨强度作用下土壤地表流失/地下漏失特征
分析不同降雨强度交替作用下土壤地表流失量随时间的变化特征(图6)可知,在降雨全过程中存在土壤地表流失量和降雨强度呈正比的关系。稳态降雨的情况下,随着降雨强度的增高,土壤地表流失量依次为2.12 g,2.37 g,4.83 g,5.92 g,且土壤地表流失速率较为稳定。
组合9,组合10的土壤地表流失特征不同,组合9降雨全过程土壤地表流失量为4.42 g,在大暴雨降雨时段里,土壤地表流失量为1.92 g,占比43%,组合10降雨全过程土壤地表流失量为7.83 g,在大暴雨降雨时段里,土壤地表流失量为3.55 g,占比45%。可以发现,30 min后不同降雨强度下土壤地下漏失速率较为稳定,前30 min的不同降雨强度中土壤漏失速率差异较大,基本呈现土壤漏失速率随降雨强度增大而增大的现象。在土壤水暂未饱和的情况下,较高的降雨强度使得土壤水快速饱和,继续降雨,土壤水难以继续入渗,在表层岩溶坡面形成土壤地表径流,雨水沿着径流通道携带土壤,造成稳态的土壤地表流失,即在图表中体现为高降雨强度具有更高的地表流失初值,之后不同降雨组合的土壤地表流失速率处于稳定状态。组合9降雨情况下的总土壤地表流失量大于稳态大雨,小于稳态暴雨情况中的土壤地表流失。而在90分钟以后,降雨强度增大到大暴雨以后,土壤地表流失速率显著增大,土壤随着形成的坡面径流通道发生流失,组合10后期降雨强度减小,但已经形成的坡面径流通道并未消失,土壤地表流失速率仍维持在相对稳定的水平。
图6 不同降雨强度下土壤地表流失量特征
分析不同降雨强度交替下土壤地下漏失量随时间的变化特征(图7)可知,总体上稳态降雨情况下不同降雨强度的土壤地下漏失量波动不大,在降雨全过程中土壤地下漏失量和降雨强度呈正相关。稳态降雨的情况下,随着降雨强度增大土壤地下漏失量依次为1.64 g,4.15 g,5.68 g,21.35 g。4组稳态降雨的情况中,除大暴雨的情况下,土壤地下漏失速率相对较稳定。
组合9、组合10具有不同的土壤地下漏失特征,土壤地下漏失量分别为18.33 g,20.24 g,组合9随着降雨强度的升高具有越来越大的地下漏失速率,不同降雨阶段下,组合9的地下漏失量均小于80 mm/h降雨强度下的稳态降雨。从总降雨量来考虑,组合9的降雨量为75 mm,40 mm/h下稳态降雨两小时降雨量为80 mm,两组的总降雨量之间相差不大,但组合9降雨情况下最后30 min内,降雨强度增大为之前的2倍,加速了坡面的土壤地表流失/地下漏失能力,分析其原因在大暴雨气候条件下,降雨强度超过一定临界值之后,雨水入渗速率加快,土壤含水率急速增加,土体结构发生改变,更容易形成裂隙流进而加速土壤地下漏失。组合10在降雨初期30 min内,相比于其他降雨情况具有最大的土壤地下漏失量11.08 g,之后随着降雨强度的减小,土壤地下漏失速率也在减小,并在60 min之后土壤地下漏失速率小于80 mm/h的稳态降雨,终态情况下土壤地下漏失量小于80 mm/h的稳态降雨,表明土壤地下漏失量与降雨总量呈正比。
图7 不同降雨强度下土壤地下漏失量特征
(1) 降雨强度与土壤地表流失量/地下漏失量均呈正比关系,但降雨强度的大小交替顺序对其影响存在差异。大雨和中雨组成的6种强度交替降雨情况下,先集中性中雨再集中性大雨产生的土壤地表流失和地下漏失总量最大;先集中性大雨再集中性中雨产生的土壤地表流失量最大,而土壤地下漏失量最小。考虑大暴雨、暴雨、大雨、中雨这4种雨强情况下,按雨强从大到小实施所产生的土壤地表流失和地下漏失总量大于雨强从小到大实施的情况。雨强大的降雨发生在降雨过程的后期更易发生土壤地下漏失。
(2) 短时多次降雨之间的时间间隔对土壤地表流失量和地下漏失量影响显著,集中性大雨比间断性大雨更易引发岩溶坡面土壤地表流失及地下漏失。短时多次降雨之间无雨的时间越长越不易发生水土流失和漏失,间断性强度低的降雨对土壤地表流失和地下漏失的影响远小于高强度集中降雨。要注意间断性中雨对土壤地下漏失的防控。
(3) 降雨强度与时长对土壤地表流失/地下漏失的速率也有影响。降雨前期,土壤地表流失速率与降雨强度呈正相关关系,在降雨中后期其速率不受降雨强度的变化影响;在整个降雨过程中,土壤地下漏失速率均与降雨强度呈正相关关系。因此,在应对短时强降雨极端天气时应该着重加强降雨初期的水土流失/漏失灾害防治。