阎思宇李斌斌于坤霞丛佩娟戴 宁
(1.西安理工大学 省部共建西北旱区生态水利国家重点实验室,陕西 西安710048;2.水利部 水土保持监测中心,北京100053)
受全球气候变化影响,中国暴雨发生频率逐年增加,受暴雨影响,土壤结构遭到破坏、肥力下降,建筑设施遭到损毁[1-2],严重可威胁民众生命财产安全[3],不仅如此,暴雨引发的洪水灾害会导致地区发生严重的水土流失现象。生产道路是农田与非农田系统重要的基础设施,在暴雨条件下是重要的泥沙来源[4-5],国内外诸多专家学者对生产道路侵蚀发生原因进行了研究探讨,表明道路自身土壤特性及路面状况会对道路侵蚀产生影响。Ziegler等[6]认为道路下垫面条件、道路边坡、道路沟渠及侵蚀沟渠等4个道路特征是影响生产道路侵蚀的主要因素,并在研究中指出路面透水性会影响径流系数,径流系数越大,即使小雨强也可能产生径流,对道路造成侵蚀[7];MacDonald等[8]总结了道路土壤性质及结构、透水性、压实度等是影响道路侵蚀强弱的主要因素;张强等[9]分析了黄土区土质路与被草土质路侵蚀状况,研究表明被草土质路侵蚀产沙量远低于土质路面;曹世雄等[10]在黄土丘陵区土质路面开展实验,表明种草不仅能美化环境,起到生态修复功效,还能有效防止道路水土流失的发生。
虽然国内外研究表明了影响道路侵蚀的主要因素,但不同区域、不同气候条件、下垫面条件下,道路侵蚀量及道路侵蚀防治研究仍不足。Wang等[11]和Zhang等[4]对暴雨后道路侵蚀量进行了试算,分析了产生道路侵蚀的原因及驱动因素,并提出水土保持措施能有效减轻道路侵蚀。水土保持措施具有改善土壤结构、蓄水保水、农田防护等基础功能[12-14],也能够在防洪减灾中发挥重要作用[15],尤其是在中国水土流失最严重的黄土高原地区[16],多名学者[17-20]研究表明,水土保持措施能够有效减少黄土高原水土流失现象。为此,调查暴雨后水土保持措施发挥的效益具有重要的实践和理论意义。
2022年7月13—16日,宁夏回族自治区西吉县出现明显降雨过程,马莲乡、什字乡出现局地强降水,有96站累计降雨量超过50 mm,有6站累计降雨量超过100 mm,在唐家河流域出现了淤地坝溃坝等安全事故,造成了严重的水土流失灾害。为深入认识暴雨条件下土壤侵蚀规律与影响,了解此次暴雨水土流失状况、土壤侵蚀强度及水土保持工程损毁状况,评价暴雨条件下水土保持各项措施安全性和有效性。2022年7月25日至29日,水利部水土保持监测中心组织开展了宁夏西吉县“7·13”暴雨水土保持综合调查,对宁夏西吉县唐家河流域土壤侵蚀状况、水土流失危害和水土保持措施效果进行全面调查和分析。本文选取西吉县唐家河流域五十岔小流域生产道路侵蚀状况进行分析,对其生产道路侵蚀开展了实地调查和测量,并采用断面法计算了五十岔小流域生产道路侵蚀模数,分析了生产道路侵蚀发生原因。
1.1 研究区概况
五十岔小流域位于宁夏西吉县东南将台乡李咀村,面积为1.07 km2,属于大陆性季风气候,位于半干旱区,多年平均气温为5.81℃,多年平均降水量为460 mm,多年平均径流深为42 mm,年径流量为1.21×104m3,年输沙量为4.58×105t。
1.2 研究方法
小流域道路侵蚀调查工作主要分为内业及外业工作,内业包括降雨前影像资料及降雨后影像资料的收集、解译和抽样工作。外业工作主要是对道路进行实地勘测。对于道路调查及计算方法均参考以往暴雨调查经验和专家学者方法[10-11],采用断面法进行调查和侵蚀强度计算。
1.2.1 实地调查方法 在暴雨过后及时采用无人机对调查小流域进行影像拍摄工作,获取基于无人机的正射影像图(DOM)和精度为8 cm的数字地面模型(DSM),并生成三维遥感影像。通过对影像资料的解译,提取小流域中所有道路的空间分布情况,从而得到整个小流域内道路总长为7.81 km,其中生产道路长度为7.62 km,占总道路长度的97.61%。按照道路主要土地管理措施(包含林地、草地和梯田)、高程和流域内分布情况作为抽样依据,从图中抽样选出11条道路,抽样道路总长为1.11 km,占流域道路总长的14.26%。流域道路分布情况及抽样道路分布情况如图1所示。调查时以10 m作为断面间隔,对抽样道路划分n个断面,并测量各断面侵蚀沟的宽、深情况,如道路划分断面出现不足10 m间隔区间时,则舍弃该区间。分别在道路起始断面、结束断面和中间随机断面进行取样,取样体积为100 cm3,用以测量道路土壤容重。
图1 五十岔小流域道路及抽样道路分布
1.2.2 道路侵蚀强度计算方法 道路侵蚀强度计算
采用断面法进行计算,计算方法为:
式中:k表示抽样道路编号;n表示k道路的断面数;分别表示抽样道路侵蚀沟的平均宽度和平均深度;wi,hi分别表示断面侵蚀沟的宽度和深度;Qk表示k道路的侵蚀模数;l表示断面区间长,本文取10 m;Ak表示道路k的面积;Lk,Wk分别表示道路k的长度和宽度;ρ表示土壤容重,经取样测量得到五十岔小流域道路土壤容重平均为1.43 g/cm3。同时采用无人机影像得到的三维影像数据对流域内道路侵蚀进行测量,以此得到的流域道路侵蚀总量作为实际道路侵蚀总量。
2.1 道路损毁情况
现场调查中发现,五十岔小流域生产道路侵蚀严重,大部分生产道路受到了不同程度的冲刷,部分道路完全损毁,仅有少部分横向道路及覆被达90%以上的纵向道路断面未受侵蚀。本次道路抽样调查总长为1.11 km,共调查断面118个,各抽样道路断面损毁统计状况详见表1。
表1 五十岔小流域抽样道路损毁状况统计
由表2统计可知,调查抽样道路断面总计118个,63个断面存在不同程度的侵蚀沟,占总断面数的53.39%,其中如W1,W2,W11号道路各断面均有切沟出现,且切沟长度较长。根据道路损毁状况进行分级,分级标准详见表2。对抽样调查道路断面区间侵蚀强度进行分级(表3)。
表2 道路损毁等级分级标准
表3 五十岔小流域抽样道路侵蚀强度分级
调查发现抽样道路侵蚀主要集中在W1,W2,W4,W11号道路,部分断面现场情况如图2所示。现场调查中发现,抽样道路中裸露黄土路面受到侵蚀状况较为严重,如图2a—2f所示,其中图2c—2f展示了道路转弯处形成的特大侵蚀沟,这些侵蚀沟深度达到3~4 m,对道路造成了严重破坏;图2g—2h展示了覆被道路的侵蚀状况,这些道路由于覆被率较高,基本不会形成侵蚀沟。
图2 五十岔小流域部分调查断面道路侵蚀现状
2.2 道路侵蚀量计算及分析
抽样道路侵蚀状况计算结果详见表4。依据表4的推算成果,根据断面法计算抽样道路侵蚀总量为161.86 t,在此基础上推算小流域的道路侵蚀总量为1 163.06 t,基于小流域的道路侵蚀模数为1 054.49 t/km2;利用无人机航拍影像建立的三维影像进行道路侵蚀浅沟和切沟的提取和计算,可以得到抽样道路的侵蚀总量为200.80 t,流域道路侵蚀总量为1 409.49 t,通过计算,可以得到流域内的道路侵蚀模数为1 277.91 t/km2。通过计算发现,现场断面法量测与三维影像提取计算的结果存在偏差。以抽样道路侵蚀总量为例,现场断面法量测道路侵蚀总量较三维影像提取计算道路侵蚀总量低了19.39%,主要是由于断面法以断面近似替代区间侵蚀状况,对侵蚀沟也是近似看作规则图形进行量测计算,若两断面间侵蚀状况发生显著变化,如断面区间侵蚀状况减弱,则计算结果可能偏大,断面区间侵蚀增大,则计算结果可能偏小,而三维影像能够较精确地量算出侵蚀沟体积,因此造成了两种方法间存在偏差、传统方法推算成果偏小的现象。
表4 五十岔小流域抽样道路侵蚀总量结果
2.3 道路侵蚀强度影响因素
通过上文侵蚀量分析,各抽样道路侵蚀总量各不相同,且在实际调查过程中所见各条道路的侵蚀状况各有差异。道路侵蚀主要受到自然因素和人类活动共同作用[21],五十岔小流域道路侵蚀的影响因素可归结为降雨因素、道路特征和人类活动3部分。
2.3.1 降雨因素 2022年7月13—16日,西吉县受强对流天气影响出现明显降雨过程,其中96站累计降雨量超过50 mm,6站累计降雨量超过100 mm,其中7月15日降雨最为集中,具有短历时、雨强大的特点。距五十岔小流域最近雨量站为什字水库站,其最大1 h降雨、最大3 h降雨、最大6 h降雨分别为43.6,91.6和115.8 mm,分别为超10 a一遇、超50 a一遇和超100 a一遇暴雨。在如此暴雨条件下,暴雨在短时间内迅速转换为地表径流,对路面造成冲刷,严重时出现侵蚀沟,是道路侵蚀发生的前置条件。
2.3.2 道路特征 许多专家学者研究发现道路是否覆被、土壤质地、走向、坡度等特征因素对道路侵蚀有着显著影响[7,21-23]。五十岔小流域道路侵蚀均发生在未硬化的生产道路路面,道路两侧土地管理措施均为梯田,且道路均无排水措施,抽样道路调查情况详见表5所示。
表5 五十岔小流域抽样道路路面及土地利用情况
(1)道路形态与土壤性质。王天巍等[22]发现路面形态通过影响路面径流的水力学特征进而影响道路侵蚀,其中凹型路面产沙最多。通过调查发现,五十岔小流域道路两侧均为梯田,因此大部分道路呈现两边高中间低的凹型形态,这种汇聚型路面形态会使降雨径流集中在低洼处,水量和水流动能呈线性或指数倍增加[24],又由于道路平均比降较大,使得降雨径流流速加快,侵蚀能力更强,而调查发现土壤构成主要为湿陷性黄土,其在天然状态下有一定承载力,但受水冲刷后,其强度变低、承载力降低,结构易受到破坏,发生沉陷,综合两种原因,使得五十岔小流域道路发生侵蚀。
(2)道路比降。抽样调查的11条道路平均比降为0.12,统计调查抽样道路侵蚀沟发生数与道路比降关系如图3所示。随道路比降增加,侵蚀沟发生个数呈先增后减的趋势,在比降为0.15~0.2区间时侵蚀沟发生数最多,占侵蚀沟发生总数的32.08%,0.2~0.25比降区间次之,占侵蚀沟发生总数的22.64%。图3也反映了侵蚀量与道路比降的关系,各比降区,侵蚀量为5 t以上的侵蚀沟主要集中发生于0.1~0.2比降区间,这些侵蚀沟主要出现在W1,W2,W4和W11号道路,通过前文分析,这些道路是位于流域上游的顺径流方向道路,周围植被覆盖率较差,比降普遍较大,因此受雨水冲刷后易出现严重的侵蚀沟。
图3 五十岔小流域道路侵蚀与比降关系
(3)道路覆被。土壤表层植被能够减弱径流动能,有效起到减弱水蚀的作用[25-26]。依据表2统计调查道路覆被情况结果显示,仅有W5,W6和W9号道路路面有覆被。由表4可知,相较于植被覆盖较差的W1,W2号道路,这些道路的侵蚀总量均较小,表明道路覆被和环境覆被能有效减少道路侵蚀的发生。尽管植被能有效减弱水蚀作用,但部分覆被较好的道路仍发生了较严重的侵蚀。部分道路覆被和环境覆被情况如图4所示。其中W5号道路覆被主要在路中,道路两侧仍为裸露地表(如图4所示),在道路两侧形成了侵蚀沟,而W9号道路因路面植被覆盖率高,几乎不发生侵蚀。表明植被覆盖位置和盖度不同,也会对道路侵蚀产生不同的影响,盖度较低或植被未覆盖的部分易发生道路侵蚀。
图4 五十岔小流域道路覆被情况
2.3.3 人类活动因素 五十岔小流域生产道路侵蚀人类活动影响主要体现在路旁土地管理措施和道路修建工程。调查发现,五十岔小流域生产道路两旁均为梯田,虽然梯田作为水土保持措施中重要的保水保土措施,能够有效削减流域水土流失现象,但此次调查发现,侵蚀量较大的道路旁的梯田基本处于未耕或育苗稀少的状态,使得土壤大面积裸露,而田面土层内均采用了覆膜措施(如图5a—5b所示),极大地增加了梯田田面的过流能力,使得这些梯田形成了天然的汇流过流平台,水流在重力作用下朝梯田田埂汇集,根据设计标准,田埂应达到30 cm高,但通过对梯田损毁情况的调查发现,田埂高度均未达标,实际高度不足20 cm(如图5c所示)。因此,水流极易冲毁田埂,汇集于道路上,通过前文分析,道路呈凹型路面,其水流不仅来自上游降雨径流,还汇集了来自梯田的汇水,因此,导致道路侵蚀更易发生,侵蚀状况急剧。五十岔小流域部分生产道路修建时有转弯路面,调查中发现这些路面发生的侵蚀状况比平直路面更为严重(如图5d所示)。路面径流在重力作用下顺平直路面汇流,当路面发生转弯时,径流会对转弯外侧路面产生冲刷后再顺路面而下,因此,路面转弯外侧发生的侵蚀状况往往较为严重。通过调查发现,五十岔小流域生产道路均未修建排水措施(如图5e所示),生产道路过流面均为道路路面。图5f为调查过程中发现较严重的道路的侵蚀,侵蚀已掏空路基,其中发现道路埋管,这部分道路内部土质较为松散;而通过观察现场进行的道路整修作业发现,整修过程中仅对侵蚀路面进行覆土填平,并没有采用压实工具对路面进行压实整平,未压实路面存在浮土,渗透能力较低,即便雨强很小时,也会在短时间内产生路面径流进而引发路面侵蚀[27]。
图5 五十岔小流域道路侵蚀人类活动影响因素
许多专家学者对农村土质道路进行了产沙模拟试验,表明土质道路受道路形态、坡度和压实度等影响,入渗率较低,易发生侵蚀[22-23]。五十岔小流域生产道路均为未硬化的黄土道路,路面整体呈现凹型,均未修建排水措施,部分路段有地下埋管,道路压实度较低,加之道路两侧梯田均未布设截排水设施,并且田面在耕作前铺设了地膜,地膜影响了水量下渗,加剧产汇流过程,梯田汇水很容易冲毁田埂顺边坡汇入道路,综合上述条件,五十岔小流域部分生产道路发生了较为严重的侵蚀现象。
根据侵蚀量计算结果,由于断面法以断面近似替代区间侵蚀状况,因此造成断面法计算结果与三维影像测量结果存在一定偏差。在实际调查过程中,受人力、物力及时间限制,无法对全流域所有道路侵蚀沟进行详细测量,采用抽样调查和断面法进行测量能节省人力、物力资源,并能在短时间内对抽样样本进行详细测量,以此来表征五十岔小流域生产道路侵蚀状况。在抽样过程中,应全方位考虑土地管理措施、分布情况,尽可能使抽样道路覆盖流域内道路所包含的特征,才不会使断面法进行侵蚀量计算时偏差过大,使得计算结果更接近实际侵蚀量,但断面法仍存在不足,若断面区间内侵蚀沟发育发生明显变化,则以断面表征区间内侵蚀状况则会产生较大误差。
五十岔小流域道路侵蚀受到降雨因素、道路特征和人类活动三方面共同作用,受到大雨强作用,地表径流急剧增大,小流域内生产道路是主要的汇流路径,同时,生产道路因其自身未硬化、坡度较大且呈凹型的路面形态,无排水、有地下埋管且未压实的不完善的修建工程和不达标的梯田修建措施,使得生产道路汇水增加,流速增大,集中成股能力增强,导致径流侵蚀力增大,对路面造成的侵蚀增大,致使小流域生产道路部分道路形成侵蚀量达5 t以上的侵蚀沟,而侵蚀沟在形态上与凹型路面相似,道路上的侵蚀沟作用于排水措施相似,径流汇聚其中,加剧了侵蚀沟的发育,最终造成了严重的水土流失现象。
通过现场的调查发现,良好的水土保持措施和植被覆盖能有效减少侵蚀现象的发生。五十岔小流域上下游、山顶和山脚植被覆盖状况、梯田内种植状况均不相同,上游梯田内作物长势较好,植被覆盖度较高,山顶采取了退耕还林还草措施,因此位于小流域上游和山顶的道路侵蚀量普遍较小,甚至不发生侵蚀;而小流域下游梯田种植密度较为稀疏,部分梯田处于未耕作状态,山脚植被覆盖率较低,因此,小流域下游及山脚处的生产道路侵蚀量较大。
(1)五十岔小流域生产道路侵蚀主要受到强降雨、生产道路自身特性和人类活动三方面因素影响,使得小流域内生产道路侵蚀状况各异,其中位于小流域下游、顺径流方向且植被盖度较低的道路侵蚀状况较为严重。
(2)基于断面法得到的流域道路侵蚀总量为1 163.06 t,基于流域道路侵蚀模数为1 054.49 t/km2,通过三维影像测量的流域道路侵蚀总量为1 409.49 t,基于流域的道路侵蚀模数为1 277.91 t/km2,调查断面损毁率达53.39%。
(3)本次调查结果反映出实施水土保持措施是必要且有效的,而且要进一步完善截排水措施,良好的水土保持措施和植被覆盖能够有效减少侵蚀现象发生。结合本次暴雨调查成果及本文的不足之处,提出以下展望:①高精度影像资料的覆盖。实地调查受限于时间紧、人力和物力资源短缺等因素,不能在短时间内对全流域全部道路展开详细测量调查,因此,高精度影像资料的拍摄和分析能够弥补实地调查的不足。通过高精度影像资料,可以清晰提取流域内各种要素,对生产道路侵蚀状况进行全方位分析,以弥补断面法产生的计算误差。②暴雨过后及时开展暴雨水土保持调查是有必要且必须的。调查能够及时发现暴雨造成的损毁情况,反馈水土保持措施发挥的功效。调查结果能对区域水土流失治理和优化调度提供参考,使得区域水土保持工作更加完善、高效。③在未来研究中,将进一步进行建模分析,采用不同雨强、坡度、土地利用、植被盖度等因素与道路侵蚀量建立回归模型及机器学习模型,探究各因子对于道路侵蚀的贡献度,为改善道路侵蚀,完善水土保持措施管理制度提供参考。
猜你喜欢梯田暴雨水土保持“80年未遇暴雨”袭首尔环球时报(2022-08-10)2022-08-10水土保持探新路 三十九年写春秋中国水土保持(2022年7期)2022-07-14暴雨小猕猴学习画刊(2022年10期)2022-01-01当暴雨突临疯狂英语·初中天地(2021年8期)2021-11-20《水土保持通报》征稿简则水土保持通报(2021年1期)2021-04-16水土保持河南水利年鉴(2020年0期)2020-06-09《山区修梯田》艺术品鉴(2019年12期)2020-01-18秀美梯田乡村地理(2018年2期)2018-09-19梯田之恋民族音乐(2017年4期)2017-09-22暴雨袭击支点(2017年8期)2017-08-22