检索项 检索词
  土壤学报  2021, Vol. 58 Issue (5): 1179-1189  DOI: 10.11766/trxb202002140046
0

引用本文  

曾荣昌, 张玉启, 何丙辉, 等. 喀斯特槽谷区岩层与坡面夹角对坡面集中流水力学特性的影响. 土壤学报, 2021, 58(5): 1179-1189.
ZENG Rongchang, ZHANG Yuqi, HE Binghui, et al. Hydraulic Characteristics of Concentrated Flow Relative to Angle between Rock Strata and Slope in the Karst Trough Valley Area. Acta Pedologica Sinica, 2021, 58(5): 1179-1189.

基金项目

国家重点研发计划课题(2016YFC0502303)和国家自然科学基金项目(41771312)资助

通讯作者Corresponding author

何丙辉, E-mail: hebinghui@swu.edu.cn

作者简介

曾荣昌(1995-), 男, 河南商城人, 硕士研究生, 研究方向为城市水土保持。E-mail: sjhpforever01@163.com
喀斯特槽谷区岩层与坡面夹角对坡面集中流水力学特性的影响
曾荣昌, 张玉启, 何丙辉, 李天阳, 曾成    
西南大学资源环境学院/三峡库区生态环境教育部重点实验室, 重庆 400715
摘要:喀斯特槽谷区地表出露的岩层与坡面形成不同夹角,塑造了复杂的微地貌,显著改变坡面集中水流路径,使得水流特性发生重大变化。目前,不同岩层与坡面夹角下集中流水力学特性的变化特征还不清楚。通过室内放水冲刷模拟试验,研究了3个坡度(10°、15°和20°)、3个冲刷流量(5、7.5、10 L·min-1)和6个岩层与坡面夹角(30°、60°、90°、120°、150°和180°(0°))组合条件下喀斯特槽谷区地表集中流的水力学特性。结果表明:集中水流雷诺数(Re)随冲刷历时变化趋势较复杂,试验条件下,Re随流量的增大而增大,随坡度变化规律不明显,变化范围为517~3 343。Darcy-Weisbach阻力系数(f)随冲刷历时表现出增大的变化趋势,其随坡度增大而增大,随流量增大而减小,变化范围为0.62~5.70。Re和岩层与坡面夹角关系均不显著,f和岩层与坡面夹角间的关系在各流量和坡度组合下表现不同。fRe之间的耦合关系受岩层与坡面夹角的显著影响,当岩层与坡面夹角小于90°时,f-Re关系用对数方程描述比较好;夹角大于等于90°时,则用幂函数方程描述更好。试验条件下,当Re < 1 791时,fRe关系不显著;当Re>1 791时,fRe则呈显著正相关。本研究可为喀斯特槽谷区土壤侵蚀预报模型的建立提供依据。
关键词喀斯特    岩层与坡面夹角    地下裂隙    雷诺数    阻力系数    集中流    
Hydraulic Characteristics of Concentrated Flow Relative to Angle between Rock Strata and Slope in the Karst Trough Valley Area
ZENG Rongchang, ZHANG Yuqi, HE Binghui, LI Tianyang, ZENG Cheng    
College of Resources and Environment, Southwest University, Key Laboratory of Eco-environments in Three Gorges Reservoir Region, Ministry of Education, Chongqing 400715, China
Abstract: 【Objective】Outcrops of rock strata in the Karst Trough Valley Area often form angles with slopes, thus making micro landforms complicated and altering significantly paths of concentrated flows, and consequently bringing about great changes in characteristics of the flows. So, it's important to study hydraulic characteristics of the concentrated flows relative to angle between rock strata and slope.【Method】To that end, a flume experiment, designed to have combinations of three slope gradients (10°, 15° and 20°), three flow rates (5, 7.5 and 10 L·min-1) and six angles between the rock strata and slope (30°, 60°, 90°, 120°, 150° and 180°(0°)), was carried out.【Result】Results show that hydraulic characteristics of a concentrated flow varied significantly with the angle between the rock stratum and slope. Reynolds number (Re) of the concentrated flow changed with duration of the scouring in a complicated trend, but exhibited a significant trend of increasing with rising flow rate, and an insignificant one of changing with slope gradient. And it varied in the range of 517~3 343 in the experiment. Darcy-Weisbach friction coefficient (f) of the concentrated flow increased with rising slope gradient and with scouring going on, but decreased with rising flow rate. And it varied in the range of 0.62~5.70 in the experiment. The relationship of Re with angles between rock strata and the slope was not significant, but that of f varied with combinations of slope gradient and flow rate. The coupling relationship between f and Re was influenced significantly by the angle between the rock stratum and the slope. It could be better described with a logarithmic equation when the angle between the rock stratum and the slope was narrower than 90°, and with a power function equation when the angle between the rock stratum and the slope was wider than or equal to 90°. In the experiment, when Re < 1 791, the relationship between f and Re was not apparent, but when Re>1 791, it was significantly positive.【Conclusion】To sum up, f and Re and their relationship vary with the angle between the rock stratum and the slope. The findings in this experiment may provide strong data support for establishing soil erosion prediction models in the Karst Trough Valley Area.
Key words: Karst    Angles between rock strata and slope    Underground fissure    Reynolds number    The resistance coefficient    Concentrated flow    

坡面集中水流可剥离、搬运大量泥沙,是坡面土壤侵蚀发生的主要动力源[1]。不同于一般明渠流,坡面集中流受降雨特性、坡度及床面状况等因素的影响,水流结构相对复杂,使准确描述其水力学特性变得更加困难[2]。水力学参数中雷诺数(Re)和Darcy-Weisbach阻力系数(f)等是表征径流水力学特性的主要指标[3-4]。目前,国内外关于各类下垫面情况下坡面水流Ref的耦合变化特征研究较多[5-8],但由于试验条件不同,二者间的数学关系还未获得统一认识[9-13]。张光辉[14]发现坡面薄层水流Ref呈良好的幂函数负相关。李占斌等[15]发现黄土区陡坡条件下Ref之间存在显著的指数关系。另有研究则认为f-Re随坡度的变化表现为不同的数学关系[3, 16]。因此,还需要更多实测数据来阐明不同条件下f-Re的耦合关系,并揭示引发二者关系变化的内在机理。

西南喀斯特槽谷区是中国南方喀斯特分布面积最大的“老、少、边、山、穷”地区,其地表岩层裸露,地下裂隙发育丰富,形成了特殊的“二元”空间结构。与其他喀斯特地貌不同,该区紧密式箱型背斜/向斜构造、岩层倾角大、典型顺/逆层边坡地质构造发育充分等特征使得地表裸露的岩层与坡面之间具有不同夹角[17],导致其存在着不同的产流产沙方式和地表径流水力学特性[18]。出露地表的岩层在降雨过程中易改变径流路径,使局部径流汇聚成集中流[19],造成剧烈的土壤侵蚀。目前,针对喀斯特槽谷区坡面水土流失方式[20-23],产流产沙变化[24-26]等方面已有大量研究,但对造成侵蚀的水动力学机制认识还不深入。对该区不同岩层与坡面夹角条件下集中流水力学特性,特别是对f-Re的耦合关系的影响研究还很少见。因此,本文通过模拟喀斯特槽谷区“二元”空间结构,采用室内土槽放水冲刷试验,研究在限定孔隙下,不同坡度、流量和岩层与坡面夹角条件下坡面集中流的流态、阻力等水力学特性随冲刷历时的变化规律,并探究岩层与坡面夹角影响下f-Re间的耦合关系,以期为建立和完善喀斯特槽谷区土壤侵蚀预报模型提供依据。

1 材料与方法 1.1 供试材料

试验土壤采自重庆市北碚区鸡公山(106°27'19" E,29°47'41" N)碳酸盐发育的石灰性土坡耕地(0~30 cm),土壤粒径组成及容重见表 1。试验前将土壤风干并过10 mm筛。为模拟喀斯特槽谷区不同岩层与坡面夹角条件,确保尽可能地接近野外实际情况,选择9块直径大于25 cm的不规则碳酸盐石块随机排列在试验钢槽中,调试并测量出露地表的面积,确保试验过程中岩石裸露率不变。因岩石体积一定,故岩土接触面面积保持固定。

表 1 试验用土壤粒径及容重 Table 1 Particle size distribution and bulk density of the test soil

试验土槽为自制的可调孔(裂)隙度的变坡钢槽(图 1),地下裂隙度为地下裂隙水平投影面积占钢槽底板总面积的百分比[24]。其计算公式[27]为:

$ \text {地下孔隙度}=S_{\text {孔}} / S_{\text {底}} $ (1)
$ S_{\text {孔 }}=\left(\pi R^{2} \arcsin \frac{1}{R}\right) / 90-R L^{2} $ (2)
图 1 试验过程(a)及装置(b) Fig. 1 Flowchart(a)and apparatus(b)of the experiment

式中,S为圆孔交错重合部分面积,m2S为土槽底板的面积,5 m2;π为圆周率,3.14;R为底板圆孔半径,0.025 m;L为两块底板圆孔重合区域的最大弦长,m。

试验的地下裂隙度可在0~5%范围内调节,土槽长5.0 m、宽1.0 m、深0.2 m,土槽上、下端分别设置地表、地下裂隙流集流槽,用于收集地表和地下径流泥沙样。前期预试验结果表明,在试验土槽宽度为1 m时,集中流及细沟发育均未达小区边界,集中流在坡面形成的细沟能够自由发育不受边界限制,故冲刷土槽长5.0 m、宽1.0 m、深0.2 m的尺寸满足试验要求。

1.2 研究方法

室内模拟试验在西南大学后山人工降雨大厅内进行。根据野外调查重庆市酉阳喀斯特槽谷区典型基岩裸露坡面,本试验设置岩石裸露率为15%,地下裂隙度为2%(轻度发育)。结合实地调查岩层与坡面夹角的大小,岩层与坡面夹角设置为30°、60°、90°、120°、150°和180°(0°),其中岩层与坡面夹角180°(0°)为对照组。岩层与坡面夹角为30°和60°时,岩石倾向与坡向一致;夹角为120°和150°时,岩石倾向与坡向相反;夹角为90°和180°(0°)时,岩石分别与坡面垂直和平行(图 2)。根据研究区典型坡度和雨强,试验坡度设计为10°、15°和20°,冲刷流量分别为5、7.5、10 L·min–1。按照岩层与坡面夹角、坡度和流量完全组合的方法进行试验,每个组合试验重复两次,共计108场试验。

图 2 岩层与坡面夹角及地下裂隙示意 Fig. 2 Diagram of angles between rock strata and slopes and underground fissure

利用钢制变坡土槽(长5 m、宽1 m、深0.5 m)开展模拟试验。试验中在水槽底部0.2 m处设置两块可调节的带圆孔钢槽地板,通过错位调节钢槽圆孔实现地下裂隙模拟。土槽上端安装有稳流水箱(长1 m、宽0.2 m、深0.2 m),用于冲刷时放水流量的控制。为保证水箱与土槽紧密连接,试验过程中不漏水,水箱放置在土槽内且用木板和纱布填充缝隙,因此实际冲刷长度为4.8 m(图 1)。

试验开始前,调节钢槽地板圆孔错位,使其裂隙达到设计值2%。试验填土厚度为20 cm,按上层15 cm(容重1.2 g·cm–3),下层5 cm(容重1.3 g·cm–3)计算填土量后分层填装。在土槽两边内壁处涂抹凡士林,两侧填土稍高于其他部分以减小边缘效应。填土过程中使用特制木板压实耙平土壤。填土完成后将直径大于25 cm的石灰岩块石按岩石裸露率15%及相应固定岩土接触面积随机排列装于土槽内。采用坡度倾斜仪将夹角调整到设计水平后,压实岩石与土壤接触面,保证与野外实际情况一致。使用侧喷式降雨器降小雨,观察到地下产生均匀不间断漏渗且地表有蓄流时停止降雨,使其达到饱和,放置48 h。

试验开始时,调整变坡水槽至设计坡度。打开出水阀门,开始放水。使用量筒对放水流量进行率定,保证放水流量与设计流量相对误差≤2.5%。待流量达到设计值且出流稳定后,调整水箱与土槽过渡段左右高度使水流集中在一侧出流,保证水流出流集中。试验中,分别在距底端1 m、2.5 m、3.5 m处设置测定断面,测定流速、水深及水宽。前期预试验表明,冲刷历时为18 min时产流产沙基本稳定,故本文冲刷时间设置为18 min。在前6分钟每分钟测定流速、水深及水宽,后12分钟每隔1.5 min测定流速、水深及水宽。采用染色法测定流速,测距为0.5 m;用薄钢尺量测水深及水宽,用温度计测量水流温度用于计算水运动黏滞系数。每场冲刷试验结束后,根据冲刷沟道深度更换表层土壤,将表面抹平打毛后为下一场试验作准备。

1.3 数据处理

染色法测得的流速为水流表面流速,根据表面流速计算所得水流流态,将测得的流速分别乘0.67、0.70和0.80的修正系数[11],从而得到水流的平均流速。本试验水流平均流速变化范围为0.17~0.46 m·s–1,标准差变化于0.0001~0.15;水流平均水深变化范围为0.50~2.40 cm,标准差变化于0.005~0.83。径流雷诺数Re是判断水流流态的重要指标,当Re < 900时,水流为层流;900 < Re < 2000时,水流为过渡流;Re > 2000时水流处于紊流状态[7],其计算公式为:

$Re = VR/v$ (3)

式中,V为水流平均流速,m·s–1R为水力半径,m;R=A/PA为过水断面面积,m2P为湿周,m;P=2h+bh为径流深,m;b为径流宽,m;ν为水运动黏滞系数,m2·s–1ν=0.01775×10–4/(1+0.0337 t+0.000221 t2);t为水温,℃。

Darcy-Weisbach阻力系数f计算公式为:

$f = 8gRS/{V^2}$ (4)

式中,g为重力加速度,m·s–2S为水力坡降,S=sinθθ为试验坡度,°。

1.4 数据分析方法

采用Excel2016进行试验数据的整理,采用SPSS24.0对Ref之间关系进行线性、对数和幂函数拟合。本文显著水平设置为:P < 0.05为关系显著,P < 0.01为关系极显著。使用Origin2018及AutoCAD2018进行作图。

2 结果与讨论 2.1 不同岩层与坡面夹角条件下径流雷诺数Re随冲刷历时的变化

不同流量、坡度和岩层与坡面夹角条件下径流雷诺数Re随冲刷时间的变化过程如图 3所示。在试验条件下,Re变化范围为517~3 343,与李占斌等[15]在陡坡条件下放水冲刷的研究结果相似。集中水流在1 min之内快速冲刷出细沟,使得每场试验Re变化范围不大。流量为5 L·min–1条件下,试验时段内各组合的径流流态为过渡流。随着坡度和流量的增大,水流重力势能和动能增加,流速增大,径流紊动性增加,径流处于以紊流为主过渡流为辅的交替状态。当流量和坡度均达到最大时,不同岩层与坡面夹角条件下的地表径流均处于紊流状态。Re随着流量的增加而增大,但随坡度变化规律不明显,这与杨大明等[28]在室内放水冲刷黄绵土的研究结果类似。这是由于单宽流量越大,单位时间通过某一断面的径流量越多,水深和流速越大,Re也越大[4]

图 3 雷诺数随冲刷时间的变化 Fig. 3 Variation of Reynolds number with duration of scouring

图 3可知,岩层与坡面夹角对Re随时间变化过程的影响较复杂,不同夹角条件下的径流雷诺数对坡度和流量的响应规律不同。10°坡面、5 L·min–1流量时,各岩层与坡面夹角下坡面集中流Re随冲刷历时总体表现出增大的趋势,这与张乐涛等[29]研究结果相似。试验前5 min内沟道形态阻力增大,随着沟道深度增加,岩石对水流阻挡作用增大,流速减小,但水力半径增大变化大于流速减小的变化,使得Re增加并始终保持发展或稳定态势[6]。岩层与坡面夹角为180°时,10°坡面各流量条件下,Re随时间变化波动性增加,变化范围为803~2 295。在其他坡度条件下其Re则随冲刷历时波动减小,说明坡度对于夹角180°条件下Re随冲刷历时变化的影响较大。岩层与坡面夹角为30°时,10°坡面、7.5 L·min–1和10 L·min–1流量下与15°坡面、5 L·min–1流量以及20°坡面全流量条件下,Re均随冲刷历时总体呈减小趋势,这是由于夹角为30°时,岩石对地表水流阻挡作用最小[30],水流受惯性力影响较大。随着坡度和流量的增大,径流能量增加,侵蚀能力增强,沟道深度及侵蚀形态复杂程度增大,沟道糙度变大,此时地下裂隙和沟道形态对地表径流雷诺数影响较大[31],流速减小快于径流深增加,Re随时间变化逐渐减小。在其他坡度和流量条件下Re随冲刷历时变化不一,规律性不强,表明岩层与坡面夹角30°的Re随冲刷历时变化趋势受坡度和流量的双重影响。在岩层与坡面夹角60°条件下,15°和20°坡面、10 L·min–1流量的Re随时间变化表现为先增大后减小,最大值出现在4~6 min内,这与李天阳等[32]研究结果类似。这是水深和流速变化综合作用的结果,在冲刷过程前期流速较大且水力半径快速增加,Re增大。在6 min以后,跌坎被冲开,壅水减少,沟道加深促进了地下漏渗,沟道中水力半径有减小的趋势,水流能量减少,流速减小较快,二者综合作用导致Re减小。10°坡面、7.5 L·min–1流量和15°坡面、5 L·min–1流量条件下,其Re随冲刷历时总体呈增大趋势,而在10°坡面、10 L·min–1流量,15°坡面、7.5 L·min–1流量与20°坡面下的5 L·min–1及7.5 L·min–1流量条件时,其Re随冲刷历时呈减小趋势。这可能是冲刷水流能量不同导致了不同的沟道侵蚀强度,在岩层与坡面夹角和地下裂隙影响下流速和水深变化不一致所造成的。岩层与坡面夹角为90°时,15°坡面、10 L·min–1流量和20°坡面、7.5 L·min–1流量的Re随冲刷历时先增大后减小,其原因与夹角60°条件下一致。15°和20°两个坡面下5 L·min–1流量和20°坡面、10 L·min–1流量的Re则随冲刷历时逐渐减小。在15°和20°两个坡面、10 L·min–1流量条件下,岩层与坡面夹角为120°、150°和180°的Re随冲刷历时变化规律较为明显,均逐渐减小,且在前4 min内减少较快。这可能是因为当夹角大于60°时,岩石倾向与坡面水流流向不同,岩石对能量较大水流阻挡效应明显。前4 min内随着冲刷进行,岩石对水流阻挡作用增强,水流黏滞力占主导作用,水深增大;但随着沟道下切,地下裂隙也使得地下漏失量增多,沟道水流能量减小,在沟道形态阻力和岩石阻挡综合作用下,流速减小快于水力半径增加,Re呈现减小趋势。在其他坡度和流量条件下,夹角120°、150°和180°的Re随冲刷历时变化规律不一,说明冲刷水流能量大小影响着不同岩层与坡面夹角的Re随冲刷历时变化趋势。

相较于岩层与坡面夹角30°和60°,当夹角大于60°时,试验过程中Re的变化更为稳定。Re的突变常伴随水流冲刷部位侵蚀作用的增强[33],这说明岩层与坡面夹角为30°和60°时水土流失方向主要是沿坡面顺流而下,随着岩层与坡面夹角增大,岩石有利于控制地表侵蚀发生的变化过程,从而减小地表侵蚀的强度,此时地表水流在地下裂隙存在的情况下,主要通过地下漏渗方式带走土壤颗粒,这与Cruz Nunes等[34]研究结果类似。

2.2 不同岩层与坡面夹角条件下径流阻力系数f随冲刷历时的变化

试验条件下阻力系数随时间变化如图 4所示。细沟侵蚀过程中水流所受阻力与流量和坡度密切相关。本试验结果显示f随坡度增大而增大,随流量增大而减小,这与白玉洁等[35]研究结果相反,与申楠等[36]研究结果类似。这可能与在地下裂隙和岩层与坡面夹角的影响下,流量和坡度变化引起径流深、相对糙度D/h以及沟道形态的变化有关[4]。坡度增大,水流势能转化为动能的能量增加,水流冲刷能力增强,沟道形态阻力增大,沟道深度增加导致地下漏渗增多,地表径流能量减小,增大了径流阻力。当流量增大时,水深增大,相对糙率减小,虽然径流冲刷改变沟道形态产生形态阻力,但流量由小变大时引起的糙率减小的程度大于形态变化引起的阻力增大的程度[3],因此f随流量的增大而减小。

图 4 阻力系数随冲刷时间的变化 Fig. 4 Variation of Darcy-Weisbach friction coefficient with duration of scouring

试验过程中f值介于0.62~5.70,这与张翔等[37]研究结果有所不同,这可能是由于本试验设置的坡度、流量等条件以及下垫面特性、土壤粒径及其排列与前人不同而造成的。随着冲刷时间的延长,各条件下的f总体表现出增大的趋势,这与前人研究结果类似[38]。试验刚开始时,集中水流仅受到颗粒阻力和岩石阻挡作用,来自细沟床面的摩擦阻力及其产生的水流紊动作用小,阻力较小。随后集中水流快速切出细沟,沟道形态阻力产生,f增大。在冲刷过程中,由于岩石阻挡作用,岩土界面处流速和水深发生突变,水流紊动性增强,加剧水流沿流程的横断面与纵降比不规则性,对沟床的局部冲刷和剪切作用增强[33, 39]。试验中观察到径流流经的岩石上部和侧部出现较深壅水及跌水,该处侵蚀程度较大,沟道形态阻力增加,这与王恒松等[19]在喀斯特野外地区观察到的侵蚀现象一致。水流经过岩石时发生明显的侧切现象,流程增加,阻力随之增大[7]。且随着侵蚀沟道深度不断增大,岩石阻挡作用进一步增强,水流路径起伏变大,水力坡度逐渐增大,沟道摩擦力对水流的阻滞作用增强,水流在运动过程中克服阻力所做的功也越大[40]。同时地下裂隙的影响程度不断增大,地下漏渗增加,地表水流能量减小,综合作用使得f增大或保持平稳趋势。由于各条件下水流能量、岩层与坡面夹角不同,使得水流携沙能力、土壤侵蚀程度和岩石阻挡作用不同[13],进而产生不同的细沟形态,在相对糙率D/h变化和冲刷后期水流对沟道的磨蚀作用影响下,各试验过程中f的变化不尽相同[4]

与岩层与坡面夹角小于90°时的f变化过程相比,夹角为120°和150°时f变化过程更为稳定。说明岩层与坡面夹角大于60°时,岩石对于水流有较强的阻挡作用,此时岩石有利于控制径流冲刷过程平稳。上述研究结果表明在喀斯特槽谷区,沟道水流阻力不仅受坡度、流量等水力条件和侵蚀形态本身的影响,还与地下裂隙和岩层与坡面夹角有关。

2.3 Ref和岩层与坡面夹角的关系

本试验研究结果表明各坡度和流量组合条件下,Re和岩层与坡面夹角关系不显著(P > 0.05),而f和岩层与坡面夹角间的关系在各流量和坡度组合下表现不同。在20°坡面、10 L·min–1流量条件下,f与岩层与坡面夹角(x)间的作用关系可表示为:

$ R^{2}=0.039 \quad P=0.706 \quad\left(\text { 含 } 30^{\circ} \text { 夹角数据 }\right) $ (5)
$ R^{2}=0.923 \quad P=0.009 \quad\left(\text { 无 } 30^{\circ} \text { 夹角数据 }\right) $ (6)

式中,f为阻力系数;x为岩层与坡面夹角,°。

由式(5)和式(6)可知,水流能量较大时,岩层与坡面夹角30°条件对于f和岩层与坡面夹角关系影响显著。结合图 4分析,这和在夹角30°时冲刷过程中的岩石对于大能量水流阻挡作用较小有关。此时集中水流冲刷能力强,试验中观察到由于能量较大,水流更易直接流经岩石光滑表面,冲击岩石下部土体,产生深度较大的跌水,剧烈的冲刷作用使得沟道形态阻力较大。同时沟道被快速下切出较大的深度,地下漏渗增加,这与严友进等[27]研究结果类似。冲刷过程中发现岩层与坡面夹角150°条件下地下产流速率从3 min的0.52 L·min–1增大至18 min的1.18 L·min–1,而夹角30°的地下产流速率从3 min的0.3 L·min–1增大至18 min的2.06 L·min–1。地表冲刷流量减小,水流流速减小,f增大。以上综合因素作用使得水流能量较大时,岩层与坡面夹角30°的f相对于其他夹角较大,从而影响着f和岩层与坡面夹角关系。

同时由式(6)可知,无夹角30°影响时,20°坡面、10 L·min–1流量条件下f和岩层与坡面夹角呈显著的幂函数正相关关系。这可能是由于岩石的存在重塑了坡面构造,随着岩层与坡面夹角的增大,岩石对大能量的水流阻挡作用也逐渐增大,不同夹角条件下的岩石通过影响沟道中集中流水力坡降从而影响坡度的实际作用效果,在水力作用下产生不同侵蚀形态和径流深,进而影响到径流阻力[41]

2.4 fRe的关系

不同岩层与坡面夹角下fRe的拟合关系见表 2。Roels[13]研究了细沟出现前后的f-Re关系,发现二者的关系可以用对数方程式描述。在本试验条件下,岩层与坡面夹角小于90°时,对数方程可较好描述fRe之间的关系;夹角大于等于90°时,幂函数方程决定系数R2较大。岩层与坡面夹角180°情况下的方程决定系数R2最小。

表 2 不同岩层与坡面夹角条件下f-Re之间的关系 Table 2 Relationships between f and Re relative to angle between the rock stratum and the slope

本研究中,岩层与坡面夹角为30°时,fRe呈正相关。当岩层与坡面夹角逐渐增大,两者呈负相关,表明存在临界岩层与坡面夹角值,使得Re增加,f开始减小。张科利[3]研究发现受坡度陡缓影响,沟道侵蚀程度及形态会有很大差异,从而使得坡度小于10°时,fRe呈指数函数负相关;而在坡面为12°和15°时,二者呈指数函数正相关,表明坡度对水流阻力的变化趋势有着决定性的影响作用。本文中,相同水力条件下,不同夹角的岩石对水流能量消减作用不同,在地下裂隙影响下造成不同的沟道冲刷形态和侵蚀程度,表明岩层与坡面夹角显著影响f-Re关系,这与甘凤玲等[30]研究结果一致。试验条件下,当Re < 1 791时,fRe关系不显著,当Re > 1 791时,fRe则呈显著正相关(图 5)。Shen和Li[42]在光滑下垫面条件下通过试验得出Re > 2 000时,fRe的增大而减小。这与本研究结果不同,可能是由于本试验为集中水流动床冲刷试验,水力条件不一致,同时特殊二元结构和岩层与坡面夹角存在改变了坡体构造,使得在紊流状态下,水流冲刷作用导致地下漏渗增多、岩石阻挡作用以及沟道形态阻力增大,fRe的增大而增大。本文研究了喀斯特槽谷区不同岩层与坡面夹角条件下集中流f-Re的变化耦合关系,但基于该区域特殊二元地质结构特性,不同地下裂隙度对于地表集中流水力学特性的影响作用还需进一步研究。

图 5 试验条件下fRe间的变化关系 Fig. 5 Regression relationships between f and Re under the experiment conditions
3 结论

试验条件下,喀斯特槽谷区坡面集中水流Re变化范围为517~3 343,其随流量的增加而增大,随坡度变化规律则不明显。岩层与坡面夹角对Re随时间变化过程的影响较复杂。10°坡面、5 L·min–1流量下,各岩层与坡面夹角的Re随冲刷历时均表现出增大趋势。在其他条件下,则出现Re随冲刷进行而减小的现象。f值介于0.62~5.70,随坡度增大而增大,随流量增大而减小。各试验条件下的f随冲刷历时呈增大的变化趋势。各坡度和流量组合条件下,Re和岩层与坡面夹角关系均不显著,f和岩层与坡面夹角间的关系在各流量和坡度组合下表现不同。fRe的耦合关系受岩层与坡面夹角的显著影响。岩层与坡面夹角小于90°时,可用对数方程式来描述二者关系;夹角大于等于90°时,用幂函数方程描述二者关系更好。当Re < 1 791时,fRe关系不显著,当Re > 1 791时,fRe则呈显著正相关。本研究可丰富土壤侵蚀水动力学理论,并为建立和完善喀斯特槽谷区土壤侵蚀预报模型提供依据。

参考文献
[1]
Poesen J, Nachtergaele J, Verstraeten G, et al. Gully erosion and environmental change: Importance and research needs[J]. Catena, 2003, 50(2/4): 91-133. (0)
[2]
Sha J D, Jiang Y J. Attempt for expounding basic dynamic characteristics of very shallow flow on preliminary eco-erosion slopes (In Chinese)[J]. Journal of Soil and Water Conservation, 1995, 9(4): 29-35. DOI:10.3321/j.issn:1009-2242.1995.04.009 [沙际德, 蒋允静. 试论初生态侵蚀性坡面薄层水流的基本动力特性[J]. 水土保持学报, 1995, 9(4): 29-35.] (0)
[3]
Zhang K L. Study of flow resistance law in rill erosion on loess slope (In Chinese)[J]. Yellow River, 1998, 20(8): 13-15, 46. [张科利. 黄土坡面细沟侵蚀中的水流阻力规律研究[J]. 人民黄河, 1998, 20(8): 13-15, 46.] (0)
[4]
Zheng L Y, Li Z B, Li P. Experimental study on hydraulic characteristics of runoff on steep slope in loess area (In Chinese)[J]. Journal of Hydraulic Engineering, 2004, 35(5): 46-51. DOI:10.3321/j.issn:0559-9350.2004.05.008 [郑良勇, 李占斌, 李鹏. 黄土区陡坡径流水动力学特性试验研究[J]. 水利学报, 2004, 35(5): 46-51.] (0)
[5]
Xia L, Song X Y, Fu N, et al. Effects of rock fragment cover on hydrological processes under rainfall simulation in a semi-arid region of China[J]. Hydrological Processes, 2018, 32(6): 792-804. DOI:10.1002/hyp.11455 (0)
[6]
Li M, Yao W Y, Yang J F, et al. Experimental study on the effect of grass cover on the overland flow pattern in the hillslope-gully side erosion system (In Chinese)[J]. Journal of Basic Science and Engineering, 2009, 17(4): 513-523. DOI:10.3969/j.issn.1005-0930.2009.04.003 [李勉, 姚文艺, 杨剑锋, 等. 草被覆盖对坡面流流态影响的人工模拟试验研究[J]. 应用基础与工程科学学报, 2009, 17(4): 513-523.] (0)
[7]
Yao W Y. Experiment study on hydraulic resistance laws of overland sheet flow (In Chinese)[J]. Journal of Sediment Research, 1996(1): 74-82. [姚文艺. 坡面流阻力规律试验研究[J]. 泥沙研究, 1996(1): 74-82.] (0)
[8]
Li P, Li Z B, Zheng L Y. Experimental study on dynamic mechanism of critical erosion on loess steep slopes (In Chinese)[J]. Journal of Sediment Research, 2008(1): 17-20. DOI:10.3321/j.issn:0468-155X.2008.01.003 [李鹏, 李占斌, 郑良勇. 黄土陡坡土壤侵蚀临界动力机制试验研究[J]. 泥沙研究, 2008(1): 17-20.] (0)
[9]
Li Y H, Niu Y B, Wang Z Z, et al. Hydrodynamic parameters and their relationships of runoff over engineering accumulation slope (In Chinese)[J]. Transactions of the Chinese Society of Agricultural Engineering, 2015, 31(22): 83-88. DOI:10.11975/j.issn.1002-6819.2015.22.012 [李永红, 牛耀彬, 王正中, 等. 工程堆积体坡面径流水动力学参数及其相互关系[J]. 农业工程学报, 2015, 31(22): 83-88.] (0)
[10]
Wang W L, Lei A L, Lin Z B, et al. Experimental study on soil erosion dynamic mechanism of slope flow in loess hilly region (In Chinese)[J]. Journal of Hydraulic Engineering, 2003, 34(9): 66-70. DOI:10.3321/j.issn:0559-9350.2003.09.012 [王文龙, 雷阿林, 李占斌, 等. 黄土丘陵区坡面薄层水流侵蚀动力机制实验研究[J]. 水利学报, 2003, 34(9): 66-70.] (0)
[11]
Abrahams A D, Parsons A J. Resistance to overland flow on desert pavement and its implications for sediment transport modeling[J]. Water Resources Research, 1991, 27(8): 1827-1836. DOI:10.1029/91WR01010 (0)
[12]
Nearing M A, Norton L D, Bulgakov D A, et al. Hydraulics and erosion in eroding rills[J]. Water Resources Research, 1997, 33(4): 865-876. DOI:10.1029/97WR00013 (0)
[13]
Roels J M. Flow resistance in concentrated overland flow on rough slope surfaces[J]. Earth Surface Processes and Landforms, 1984, 9(6): 541-551. DOI:10.1002/esp.3290090608 (0)
[14]
Zhang G H. Study on hydraulic properties of shallow flow (In Chinese)[J]. Advances in Water Science, 2002, 13(2): 159-165. DOI:10.3321/j.issn:1001-6791.2002.02.005 [张光辉. 坡面薄层流水动力学特性的实验研究[J]. 水科学进展, 2002, 13(2): 159-165.] (0)
[15]
Li Z B, Qin B S, Kang W, et al. Indoor experimental studies on hydrodynamic characteristics of runoff in rill erosion procession steep slope (In Chinese)[J]. Transactions of the Chinese Society of Agricultural Engineering, 2008, 24(6): 64-68. DOI:10.3321/j.issn:1002-6819.2008.06.013 [李占斌, 秦百顺, 亢伟, 等. 陡坡面发育的细沟水动力学特性室内试验研究[J]. 农业工程学报, 2008, 24(6): 64-68.] (0)
[16]
Abrahams A D, Li G, Parsons A J. Rill hydraulics on a semiarid hillslope, southern Arizona[J]. Earth Surface Processes and Landforms, 1996, 21(1): 35-47. DOI:10.1002/(SICI)1096-9837(199601)21:1<35::AID-ESP539>3.0.CO;2-T (0)
[17]
Jiang Y J, Liu X M, He S Y, et al. Research and development of comprehensive rehabilitation measures for land rocky desertification in Karst trough valley area (In Chinese)[J]. Acta Ecologica Sinica, 2016, 36(22): 7092-7097. [蒋勇军, 刘秀明, 何师意, 等. 喀斯特槽谷区土地石漠化与综合治理技术研发[J]. 生态学报, 2016, 36(22): 7092-7097.] (0)
[18]
Gan F L, He B H, Qin Z Y. Hydrological response and soil detachment rate from dip/anti-dip slopes as a function of rock strata dip in Karst valley revealed by rainfall simulations[J]. Journal of Hydrology, 2020, 581: 124416. DOI:10.1016/j.jhydrol.2019.124416 (0)
[19]
Wang H S, Xiong K N, Zhang F M. Mechanism study on effects of terrain on soil erosion of Karst slope (In Chinese)[J]. Bulletin of Soil and Water Conservation, 2015, 35(4): 1-7. [王恒松, 熊康宁, 张芳美. 地形因子对喀斯特坡面水土流失影响的机理研究[J]. 水土保持通报, 2015, 35(4): 1-7.] (0)
[20]
Zhang X B, Wang S J, Cao J H, et al. Characteristics of water loss and soil erosion and some scientific problems on Karst rocky desertification in southwest China Karst area (In Chinese)[J]. Carsologica Sinica, 2010, 29(3): 274-279. DOI:10.3969/j.issn.1001-4810.2010.03.009 [张信宝, 王世杰, 曹建华, 等. 西南喀斯特山地水土流失特点及有关石漠化的几个科学问题[J]. 中国岩溶, 2010, 29(3): 274-279.] (0)
[21]
Peng X D, Dai Q H, Li C L. Research progress on the process and mechanism of soil water loss or leakage on slope in Southwest Karst of China (In Chinese)[J]. Journal of Soil and Water Conservation, 2017, 31(5): 1-8. [彭旭东, 戴全厚, 李昌兰. 中国西南喀斯特坡地水土流失/漏失过程与机理研究进展[J]. 水土保持学报, 2017, 31(5): 1-8.] (0)
[22]
Akdim B. Karst landscape and hydrology in Morocco: Research trends and perspectives[J]. Environmental Earth Sciences, 2015, 74(1): 251-265. DOI:10.1007/s12665-015-4254-5 (0)
[23]
Dai Q H, Peng X D, Wang P J, et al. Surface erosion and underground leakage of yellow soil on slopes in Karst regions of southwest China[J]. Land Degradation & Development, 2018, 29(8): 2438-2448. (0)
[24]
Peng X D, Dai Q H, Li C L, et al. Effect of simulated rainfall intensities and underground pore fissure degrees on soil nutrient loss from slope farmlands in Karst region (In Chinese)[J]. Transactions of the Chinese Society of Agricultural Engineering, 2017, 33(2): 131-140. [彭旭东, 戴全厚, 李昌兰, 等. 模拟雨强和地下裂隙对喀斯特地区坡耕地养分流失的影响[J]. 农业工程学报, 2017, 33(2): 131-140.] (0)
[25]
Dai Q H, Peng X D, Zhao L S, et al. Effects of underground pore fissures on soil erosion and sediment yield on Karst slopes[J]. Land Degradation & Development, 2017, 28(7): 1922-1932. (0)
[26]
Peng X D, Dai Q H, Yang Z, et al. Sediment yield of surface and underground erosion in the process of rocky desertification of Karst area (In Chinese)[J]. Acta Pedologica Sinica, 2016, 53(5): 1237-1248. [彭旭东, 戴全厚, 杨智, 等. 喀斯特山地石漠化过程中地表地下侵蚀产沙特征[J]. 土壤学报, 2016, 53(5): 1237-1248.] (0)
[27]
Yan Y J, Dai Q H, Fu W B, et al. Runoff and sediment production processes on a Karst bare slope (In Chinese)[J]. Acta Ecologica Sinica, 2017, 37(6): 2067-2079. [严友进, 戴全厚, 伏文兵, 等. 喀斯特裸坡产流产沙过程试验研究[J]. 生态学报, 2017, 37(6): 2067-2079.] (0)
[28]
Yang D M, Gao P L, Liu X Y, et al. Hydraulic characteristics of slope land under rill erosion (In Chinese)[J]. Journal of Hydroelectric Engineering, 2017, 36(10): 65-73. DOI:10.11660/slfdxb.20171007 [杨大明, 高佩玲, 刘小媛, 等. 细沟侵蚀下的坡面水动力学特性研究[J]. 水力发电学报, 2017, 36(10): 65-73.] (0)
[29]
Zhang L T, Gao Z L, Tian H W. Hydro-dynamic characteristics of steep slope runoff of engineering accumulation under simulated runoff conditions (In Chinese)[J]. Journal of Soil and Water Conservation, 2013, 27(4): 34-38. [张乐涛, 高照良, 田红卫. 工程堆积体陡坡坡面径流水动力学特性[J]. 水土保持学报, 2013, 27(4): 34-38.] (0)
[30]
Gan F L, He B H, Qin Z Y. Effect of bedding slope on hydrodynamic parameters in typical Karst Valley (In Chinese)[J]. Acta Pedologica Sinica, 2019, 56(4): 825-837. [甘凤玲, 何丙辉, 覃自阳. 喀斯特槽谷区的顺/逆层坡面对水动力学参数的影响[J]. 土壤学报, 2019, 56(4): 825-837.] (0)
[31]
Wang L S, Cai Q G, Cai C F, et al. Morphological changes of rill on loess slope and its relationship with flow velocity (In Chinese)[J]. Transactions of the Chinese Society of Agricultural Engineering, 2014, 30(11): 110-117. DOI:10.3969/j.issn.1002-6819.2014.11.014 [王龙生, 蔡强国, 蔡崇法, 等. 黄土坡面细沟形态变化及其与流速之间的关系[J]. 农业工程学报, 2014, 30(11): 110-117.] (0)
[32]
Li T Y, He B H, Lei T W, et al. Indoor experiments on eroding rills hydrodynamics characteristics of earth-rock mixed landslide slope in Wenchuan Earthquake area (In Chinese)[J]. Journal of Hydraulic Engineering, 2014, 45(8): 892-902. [李天阳, 何丙辉, 雷廷武, 等. 汶川震区滑坡堆积体土石混合坡面细沟水动力学特征室内试验[J]. 水利学报, 2014, 45(8): 892-902.] (0)
[33]
Ding W F, Li Z B, Lu K X, et al. The elementary study of the reason of rill erosion on slope (In Chinese)[J]. Acta Pedologica Sinica, 2003, 40(6): 822-828. DOI:10.3321/j.issn:0564-3929.2003.06.004 [丁文峰, 李占斌, 鲁克新, 等. 坡面细沟发生临界水动力条件初探[J]. 土壤学报, 2003, 40(6): 822-828.] (0)
[34]
Cruz Nunes F, Delunel R, Schlunegger F, et al. Bedrock bedding, landsliding and erosional budgets in the Central European Alps[J]. Terra Nova, 2015, 27(5): 370-378. DOI:10.1111/ter.12169 (0)
[35]
Bai Y J, Zhang F B, Yang M Y, et al. Variation of hydraulic parameters of shallow flow on steep loess slope (In Chinese)[J]. Acta Pedologica Sinica, 2018, 55(3): 641-649. [白玉洁, 张风宝, 杨明义, 等. 急陡黄土坡面薄层水流水力学参数变化特征[J]. 土壤学报, 2018, 55(3): 641-649.] (0)
[36]
Shen N, Wang Z L, Liu J E, et al. Experimental study of rill flow resistance on loess hillslope (In Chinese)[J]. Bulletin of Soil and Water Conservation, 2012, 32(6): 34-37. [申楠, 王占礼, 刘俊娥, 等. 黄土坡面细沟水流阻力的试验研究[J]. 水土保持通报, 2012, 32(6): 34-37.] (0)
[37]
Zhang X, Gao Z L, Yuan X H, et al. Study on hydraulic parameters of rill flow on slopes of engineering accumulation (In Chinese)[J]. Journal of Sediment Research, 2016(4): 34-40. [张翔, 高照良, 袁雪红, 等. 工程堆积体坡面细沟流水力学参数特性研究[J]. 泥沙研究, 2016(4): 34-40.] (0)
[38]
Li M, Yao W Y, Chen J N, et al. Experimental study on runoff resistance of hilly slopegullied surface with grass coverage (In Chinese)[J]. Journal of Hydraulic Engineering, 2007, 38(1): 112-119. DOI:10.3321/j.issn:0559-9350.2007.01.017 [李勉, 姚文艺, 陈江南, 等. 草被覆盖下坡面-沟坡系统坡面流阻力变化特征试验研究[J]. 水利学报, 2007, 38(1): 112-119.] (0)
[39]
余新晓, 张晓明, 李建劳, 等. 土壤侵蚀过程与机制[J]. 北京: 科学出版社, 2009.
Yu X X, Zhang X M, Li J L, et al. Process and mechanism of soil erosion[J]. Beijing: Science Press, 2009. (0)
[40]
Liang X L, Zhao L S, Wu J, et al. Simulation of response law for soil surface roughness and hydraulics parameters of runoff (In Chinese)[J]. Transactions of the Chinese Society of Agricultural Engineering, 2014, 30(19): 123-131. DOI:10.3969/j.issn.1002-6819.2014.19.015 [梁心蓝, 赵龙山, 吴佳, 等. 地表糙度与径流水力学参数响应规律模拟[J]. 农业工程学报, 2014, 30(19): 123-131.] (0)
[41]
Zhang K L. Hydrodynamic characteristics of rill flow on loess slopes (In Chinese)[J]. Journal of Sediment Research, 1999(1): 57-62. [张科利. 黄土坡面发育的细沟水动力学特征的研究[J]. 泥沙研究, 1999(1): 57-62.] (0)
[42]
Shen H W, Li R M. Rainfall effect on sheet flow over smooth surface[J]. Journal of the Hydraulics Division, 1973, 99(5): 771-792. DOI:10.1061/JYCEAJ.0003646 (0)