喀斯特槽谷區(qū)坡面集中水流曼寧系數(shù)變化特征

曾榮昌,李天陽,何丙輝

(西南大學(xué) 資源環(huán)境學(xué)院 三峽庫區(qū)生態(tài)環(huán)境教育部重點實驗室，重慶 400715)

坡面水流水動力學(xué)變化對坡面產(chǎn)流產(chǎn)沙特性有決定作用，常用的水流阻力參數(shù)包括Darcy-Weisbach阻力系數(shù)f和曼寧系數(shù)n，其中曼寧系數(shù)n常用于溝渠流的阻力分析[1]。了解水流阻力影響因素及其變化規(guī)律是認(rèn)識坡面土壤侵蝕過程和動力學(xué)機(jī)制的基礎(chǔ)，有助于構(gòu)建坡面土壤侵蝕預(yù)報模型[2-3]。關(guān)于水流阻力的研究中，張科利等[4]研究表明受溝道沖刷形態(tài)的影響，n的變化與細(xì)溝中的徑流量及坡度的陡緩有關(guān)。湯珊珊等[5]研究發(fā)現(xiàn)相比較于其他水動力學(xué)參數(shù)，相對水深(h)曼寧系數(shù)n/h可以更好預(yù)測水力侵蝕過程中輸沙量和產(chǎn)流。Xia等[6]研究發(fā)現(xiàn)隨著巖石覆蓋率的增加n和f均隨之增大。同時，前人研究發(fā)現(xiàn)水流阻力與其他水力學(xué)參數(shù)之間存在一定的數(shù)學(xué)關(guān)系[7-10]。如拜亞茹等[7]研究發(fā)現(xiàn)糙率系數(shù)n與水深h的關(guān)系受流態(tài)的影響。李朝棟等[10]研究表明n與雷諾數(shù)Re和弗汝德數(shù)Fr之間存在著冪函數(shù)關(guān)系。以上研究多關(guān)注于土壤結(jié)構(gòu)單一、空間結(jié)構(gòu)為單層的非喀斯特區(qū)域，而對于具有特殊二元地質(zhì)結(jié)構(gòu)的喀斯特地區(qū)的相關(guān)研究還不多見。西南喀斯特槽谷區(qū)緊密式箱型背斜/向斜構(gòu)造、巖層傾角大、典型順/逆向邊坡地質(zhì)構(gòu)造發(fā)育充分等特征造成不同傾角巖層出露地面，使地表裸露巖石沿坡向與斜坡坡面形成不同夾角[11]，導(dǎo)致其存在著不同的產(chǎn)流產(chǎn)沙方式和地表徑流水力學(xué)特性[12]，且在降雨過程中易匯聚局部徑流成集中流[13]，造成劇烈的土壤侵蝕。目前，已有較多針對喀斯特槽谷區(qū)土壤侵蝕的研究，但這些研究多集中于水土流失方式[14-15]，產(chǎn)流產(chǎn)沙變化[16-17]等方面，從水動力學(xué)角度闡述喀斯特槽谷區(qū)土壤侵蝕特征的研究還不多，特別是關(guān)于水流曼寧系數(shù)變化特征的研究更是少見。為此，本文通過模擬喀斯特槽谷區(qū)不同巖石與斜坡夾角條件坡面和地下裂隙雙層空間結(jié)構(gòu)，采用室內(nèi)土槽放水沖刷試驗，研究在限定裂隙下，不同坡度、流量和巖石與斜坡坡面夾角條件下的坡面集中流曼寧系數(shù)變化特征，揭示了喀斯特槽谷區(qū)坡面集中流曼寧系數(shù)隨時間變化規(guī)律及其與其他水力學(xué)參數(shù)之間的關(guān)系，以期為建立喀斯特槽谷區(qū)土壤侵蝕預(yù)報模型提供科學(xué)依據(jù)。

1 試驗方法

1.1 供試儀器

試驗土槽為自制的可調(diào)孔(裂)隙度的變坡鋼槽(見圖1)，土槽長5.0 m，寬1.0 m，深0.2 m。通過錯位調(diào)節(jié)土槽底部0.2 m處兩塊可調(diào)節(jié)的帶圓孔鋼槽地板來模擬地下裂隙，地下裂隙度為地下裂隙水平投影面積占鋼槽底板總面積的百分比[12]，可在0%～5%范圍內(nèi)調(diào)節(jié)。土槽上、下端分別設(shè)置地表徑流、地下裂隙流集流槽，沖刷水流從置于土槽頂端的穩(wěn)流水箱(長1.0 m，寬0.2 m，深0.2 m)出流，試驗沖刷長度為4.8 m。前期預(yù)試驗表明，土槽寬度為1.0 m的條件下，集中流在坡面形成的細(xì)溝能夠自由發(fā)育不受邊界限制，土槽尺寸滿足試驗要求。

注：圖中坡度α=10°，15°，20°；巖石與坡面夾角β=30°，60°，90°，120°，150°，180°(0°)。

1.2 供試材料處理

試驗土壤于2018年6月采自重慶市北碚區(qū)雞公山(106°27′19″E，29°47′41″N)喀斯特槽谷區(qū)碳酸鹽發(fā)育的石灰性土坡耕地土壤，土壤粒徑組成及容重詳見表1。試驗前將土壤風(fēng)干并過10 mm篩，計算填土量后按上層15 cm(容重1.2 g/cm3)，下層5 cm(容重1.3 g /cm3)分層填裝土壤于試驗土槽內(nèi)。土槽邊界處涂抹凡士林，填土?xí)r壓實耙平土壤，兩邊填高以減小邊緣效應(yīng)。填土完成后，為模擬喀斯特槽谷區(qū)野外實際條件，采用坡度傾斜儀將9塊直徑大于25 cm的不規(guī)則碳酸鹽石塊按一定巖石裸露率及設(shè)計夾角角度隨機(jī)排列裝于土槽內(nèi)，壓實巖石與土壤接觸面。

表1 試驗用土壤機(jī)械組成及容重

1.3 試驗方案

室內(nèi)模擬試驗于2019年4月在西南大學(xué)后山人工降雨大廳內(nèi)進(jìn)行。根據(jù)野外調(diào)查重慶市酉陽喀斯特槽谷區(qū)典型基巖裸露坡面，本試驗設(shè)置巖石裸露率為15%，地下裂隙度為2%(輕度發(fā)育)，巖石與坡面夾角設(shè)置為30°，60°，90°，120°，150°和180°。夾角為30°和60°的巖石傾向與坡向一致；夾角120°和150°的巖石傾向與坡向相反；夾角為90°和180°時，巖石分別與坡面垂直和平行。根據(jù)研究區(qū)典型坡度和雨強(qiáng)，試驗坡度設(shè)計為10°，15°和20°，流量設(shè)計為5，7.5，10 L/min。所有試驗均重復(fù)2次，共計108場試驗。

試驗開始前，調(diào)整地下裂隙度達(dá)到設(shè)計水平2%。使用側(cè)噴式降雨器降小雨至地下產(chǎn)生均勻不間斷漏滲且地表有蓄流時，停止降雨，使土槽中土壤達(dá)到飽和，放置48 h。試驗開始時，將坡度及流量均調(diào)整至設(shè)計水平。待出流穩(wěn)定后，調(diào)整水箱與土槽過渡段左右高度使水流集中在一側(cè)出流。試驗沖刷歷時為18 min，在前6 min每隔1 min分別于距底端1，2.5，3.5 m處設(shè)置的測定斷面測定流速、水深及水寬，后12 min每隔1.5 min在相同位置測定流速、水深及水寬。采用染色法測定水流流速，測距為0.5 m；用薄鋼尺量測水深及水寬，用溫度計測量水流溫度用于計算水運動黏滯系數(shù)。試驗結(jié)束后，根據(jù)沖刷溝道深度更換表層土壤，將表面抹平打毛后為下一場試驗作準(zhǔn)備。

1.4 數(shù)據(jù)處理及分析

曼寧系數(shù)n計算公式為：

n=V-1·S1/2·R2/3

(1)

式中：V為水流平均流速(m/s)；R為水力半徑(m)，R=A/P，A為過水?dāng)嗝婷娣e(m2)；P為濕周(m)，P=2h+b，h為流深(m)；b為徑流寬(m)；S為水力坡降，S=sinθ，θ為試驗坡度(°)。

徑流雷諾數(shù)Re是判斷水流流態(tài)的重要指標(biāo)，當(dāng)Re<900時，水流為層流；9002 000時水流處于紊流狀態(tài)[18]。其計算公式為：

Re=VR/v

(2)

(ν=0.017 75×10-4/(1+0.033 7t+0.000 221t2)

式中：t為水溫(℃)。

弗汝德數(shù)Fr計算公式為：

(3)

式中：g為重力加速度(m/s2)。

1.5 數(shù)據(jù)分析方法

采用Excel 2016進(jìn)行試驗數(shù)據(jù)的整理。利用SPSS 24.0進(jìn)行數(shù)據(jù)分析：采用單因素方差(one-way ANOVA)對不同巖石與坡面夾角條件的水流曼寧系數(shù)和相對水深曼寧系數(shù)進(jìn)行分析，用Duncan 法進(jìn)行多重比較；對n與h，Re，F(xiàn)r以及n/h與Re，F(xiàn)r進(jìn)行回歸分析并檢驗方程的顯著性。本文顯著水平設(shè)置為：p<0.05為關(guān)系顯著，p<0.01為關(guān)系極顯著。使用Origin 2018及AutoCAD 2018作圖。

2 結(jié)果與分析

2.1 巖石與坡面夾角對曼寧系數(shù)n的影響

試驗條件下曼寧系數(shù)n隨沖刷歷時的變化如圖2所示。由圖2可知，n在試驗前6 min內(nèi)迅速增大，之后隨著沖刷歷時的延長呈緩慢增長或趨于穩(wěn)定的趨勢。這是因為喀斯特槽谷區(qū)地表裸露巖石與坡面形成不同夾角，在試驗條件下，前6 min內(nèi)集中流迅速沖刷出溝道，水流經(jīng)過巖石時由于阻擋和巖土界面發(fā)生改變，水深和流速隨之變化，混亂度增大，易發(fā)生劇烈侵蝕，巖石下部產(chǎn)生較深跌水，溝道糙度增加，這與在野外所觀察到的現(xiàn)象一致[13]。同時溝道中徑流實際水力坡降隨著溝道深度增大而不斷增大，巖石對水流阻擋作用增強(qiáng)，流速減小，使得n明顯增大[19-20]。隨著沖刷繼續(xù)進(jìn)行，溝道加深，在地下裂隙影響下地下漏滲增加[21]，地表集中流動能減小，水流阻力較大，沖刷能力減小，溝道形態(tài)變化程度較小，n呈緩慢增長或趨于穩(wěn)定[22]。試驗條件下，n隨流量的增大而減小。這可能是由于流量變化使得相對糙率和溝道形態(tài)發(fā)生變化，其中糙率減小的程度大于形態(tài)變化引起的增大程度，二者綜合作用使得n隨流量的增大而減小[4]。試驗條件下，水流曼寧系數(shù)變化于0.036～0.116，這與傅斌等[23]研究結(jié)果相比較小，而相對于田風(fēng)霞等[24]研究結(jié)果較大。原因與喀斯特槽谷區(qū)下墊面條件及坡面集中水流特性有關(guān)?？λ固夭酃葏^(qū)坡面巖石改變了地表微地形，增大了地表糙度，同時地下孔隙也減少了地表徑流能量，使得該區(qū)域水流曼寧系數(shù)相對于非喀斯特區(qū)域較大；而坡面集中水流具有流速和水深較大，水流能量較集中等特征，這與在降雨條件下形成的坡面流不同。以上因素使得本試驗研究結(jié)果相較于前人有所不同。

圖2 坡面集中水流曼寧系數(shù)n隨沖刷歷時的變化

不同巖石與坡面夾角對n的影響如圖3所示。夾角角度分別為30°，60°，90°，120°，150°和180°時，n平均值分別為0.070，0.064，0.072，0.068，0.075，0.077。由圖3可知，巖石與坡面夾角小于等于90°時n隨夾角角度增大而先減小后增大；巖石與坡面夾角大于90°時，n隨夾角角度增大而增大。這種現(xiàn)象可能和與坡面形成不同夾角的巖石影響坡面水流運動及其能量變化進(jìn)而影響到溝道形態(tài)發(fā)育有關(guān)。巖石與坡面夾角小于等于90°時，巖石傾向與坡向相同，此時巖石對坡面水流有促進(jìn)作用[25]，水流流速相對較大，在夾角30°表現(xiàn)較為明顯。試驗觀察到此時巖石更易影響水流集中狀態(tài)，細(xì)溝形態(tài)復(fù)雜，水流濕周增大。由于巖石阻力較小，水流能量較大，水流巖石周圍易產(chǎn)生跌水和壅水，溝道糙度增大，n反而較大。而隨著夾角增大，巖石分散水流能力減弱，但阻擋作用也增大，流速減小，水流阻力增大，n隨之增大。巖石與坡面夾角大于90°時，隨著夾角角度增大，巖石維持水流處于集中狀態(tài)，阻擋作用增強(qiáng)。水深有所增大的同時，流速減小較多，使得水流總體能量隨之減小。水流路徑在巖石影響下更易彎曲，產(chǎn)生側(cè)切，流程增加，水流阻力和過水?dāng)嗝娲植诙入S之增大[18]，水流克服巖石阻擋做功增加，水流流速減小，n增大。巖石與坡面夾角為180°時徑流路徑最長，在濕周和巖石阻擋能力較大的綜合作用下，水流阻力和過水?dāng)嗝娌诙容^大，此時n為最大[10]。對試驗條件下巖石與坡面夾角(x)與n進(jìn)行擬合得到：

n=7.3×10-7x2-9.3×10-5x+0.071

(R2=0.626，p=0.228)

(4)

根據(jù)公式(4)可以看出，二者關(guān)系不顯著。

注：不同小寫字母表示不同巖石與坡面夾角間差異顯著性(p<0.05)。下同。

2.2 曼寧系數(shù)n與水深h的關(guān)系

試驗條件下n與水深h的關(guān)系如表2所示。通過對不同坡度、流量和巖石與坡面夾角條件下的n與h進(jìn)行函數(shù)擬合發(fā)現(xiàn)，坡度為10°，流量為10 L/min，夾角為150°時的決定系數(shù)最大。由表2可知，坡度為15°，流量為7.5 L/min，夾角為30°和坡度為20°，流量為5 L/min、夾角為90°條件下，n隨h的增大而減小，但二者關(guān)系均不顯著。在其他條件下，n隨h的增大而增大，且存在著顯著的冪函數(shù)正相關(guān)關(guān)系，這與Fraga等[26]研究結(jié)果不一致。這可能是因為本試驗為集中出流，水流在試驗前期快速沖刷出溝道，隨著細(xì)溝不斷加深，侵蝕溝內(nèi)產(chǎn)生了一定間距的跌水，在巖石影響下集中水流沿溝道能量分布不具有連續(xù)性，存在水深和侵蝕突變點[27]，水流所塑造的溝道形態(tài)也更復(fù)雜。巖石本身增強(qiáng)了水流本身的擾動，使得水深增大的同時增強(qiáng)水流阻力作用[28]。巖石與坡面夾角30°條件下n與h擬合效果較差，決定系數(shù)R2變化于0.025～0.801；夾角為60°，90°，120°，150°和180°時，二者擬合效果較好，相關(guān)性顯著。這與巖石與坡面夾角30°條件下巖石更易分散和改變徑流路徑，在地下裂隙和巖石影響下沖刷過程和溝道發(fā)育更為復(fù)雜有關(guān)。

表2 坡面集中水流曼寧系數(shù)n與水深h的關(guān)系

2.3 曼寧系數(shù)n與雷諾數(shù)Re、弗汝德數(shù)Fr的關(guān)系

在試驗條件下，對n與Re分別進(jìn)行直線、對數(shù)、冪函數(shù)擬合，結(jié)果如表3所示。由表3可知，總體上，直線方程的決定系數(shù)R2比對數(shù)方程和冪函數(shù)方程要大，直線方程更適合用來描述n與Re的關(guān)系，這與Emmett等[29]研究結(jié)果一致。n與Re之間關(guān)系受坡度、流量和巖石與坡面夾角改變的影響，不呈現(xiàn)特定的正負(fù)相關(guān)函數(shù)關(guān)系。這可能是隨著沖刷進(jìn)行在地下裂隙、溝道糙度以及巖石阻擋的作用下流速和水深綜合作用的結(jié)果。直線方程的R2普遍不高，其中最大值出現(xiàn)在坡度15°，流量7.5 L/min，夾角為180°情況下，為0.882，最小值出現(xiàn)在坡度10°，流量5 L/min，夾角為90°情況下，為0.001。說明在本試驗條件下，二者關(guān)系不顯著(p>0.05)，雷諾數(shù)Re不適合用來描述曼寧系數(shù)n，這可能是受巖石與坡面夾角和喀斯特槽谷區(qū)特殊二元地質(zhì)結(jié)構(gòu)影響，使得試驗過程中水深和流速變化復(fù)雜，Re變化規(guī)律不明顯，從而影響上述關(guān)系。王仁新等[28]的研究也獲得類似結(jié)果。

表3 曼寧系數(shù)n與雷諾數(shù)Re之間的關(guān)系

對各巖石與坡面夾角條件下n與Fr關(guān)系進(jìn)行回歸分析，得到關(guān)系式如表4所示。

表4 各巖石與坡面夾角條件下曼寧系數(shù)n與弗汝德數(shù)Fr的關(guān)系

由表4可知夾角為60°時二者關(guān)系為極顯著(p<0.01)，夾角為90°時二者關(guān)系顯著(p<0.05)，其他條件下不顯著(p>0.05)。說明巖石與坡面夾角影響著二者關(guān)系。n與Fr總關(guān)系式為

(5)

由公式(5)可知，在總體試驗條件下，n與Fr呈極顯著的冪函數(shù)負(fù)相關(guān)關(guān)系。這可能是由于坡面徑流受到的巖石阻擋作用和溝道形態(tài)阻力增大，水流受到的阻滯作用越強(qiáng)，水流將越緩[30]，這與郝好鑫等[31]研究結(jié)果一致。

2.4 巖石與坡面夾角對相對水深曼寧系數(shù)n/h的影響

由不同巖石與坡面夾角對n/h的影響(見圖4)可知，6個夾角角度條件下的平均n/h呈先增大后減小的變化趨勢。巖石與坡面夾角小于等于90°時，各夾角之間表現(xiàn)出顯著性差異，n/h隨著夾角角度增大而增大。這可能是由于隨著巖石與坡面夾角角度增加，巖石對水流產(chǎn)生的阻擋作用及相應(yīng)溝道形態(tài)阻力的增大大于巖石對水流影響而造成的水深變化。Yang等[32]研究結(jié)果表明，n/h能夠較好地預(yù)測徑流能量，且與滾波速度呈反比。本試驗在巖石與坡面夾角小于等于90°的條件下，隨著夾角減小，巖石阻擋作用減小，巖石處由于水流被抬高以及水流流經(jīng)的溝道界面糙度突變，可能更易產(chǎn)生滾波。當(dāng)巖石與坡面夾角大于90°時，各夾角之間沒有顯著差異，而夾角90°與120°之間則表現(xiàn)出了顯著差異，說明試驗條件下，當(dāng)巖石與坡面夾角大于等于90°時存在臨界夾角，角度在90°～120°之間。對試驗條件下巖石與坡面夾角(x)與n/h進(jìn)行擬合得到：

圖4 不同巖石與坡面夾角下n/h的變化

n/h=-2×10-4x2+0.033x+4.013

(R2=0.675，p=0.185)

(6)

根據(jù)公式(6)可以看出，二者關(guān)系不顯著。

2.5 相對水深曼寧系數(shù)n/h與Re，F(xiàn)r的關(guān)系

不同巖石與坡面夾角條件下，n/h與Re關(guān)系如圖5所示。

圖5 不同巖石與坡面夾角的n/h與Re關(guān)系

由圖5可知，n/h與Re均呈冪函數(shù)負(fù)相關(guān)關(guān)系。巖石與坡面夾角為90°和120°時二者關(guān)系為極顯著(p<0.01)，R2分別為0.830，0.729；巖石與坡面夾角為60°，150°和180°時二者關(guān)系顯著(p<0.05)，R2分別為0.608，0.483，0.486；巖石與坡面夾角30°時二者關(guān)系則不顯著(p>0.05)。這可能是因為夾角30°時巖石促進(jìn)了水流的流動，在水流能量較大的條件下更易加重侵蝕，促進(jìn)崩塌和細(xì)溝流路的發(fā)育，使徑流路徑更為復(fù)雜，沖刷過程中溝道侵蝕形態(tài)變化的不確定性增加，減弱了二者關(guān)系[33]。

不同巖石與坡面夾角的n/h與Fr的關(guān)系如圖6所示。由圖6可知，n/h與Fr關(guān)系不顯著(p>0.05)。Guo等[34]研究結(jié)果表明，n/h和雷諾數(shù)Re均可較好預(yù)測徑流波動，同時其能夠更好預(yù)測泥沙輸移，且n/h與Fr關(guān)系不顯著，這與本試驗結(jié)果一致。

圖6 不同巖石與坡面夾角的n/h與Fr關(guān)系

3 結(jié) 論

(1) 巖石與坡面夾角對喀斯特槽谷區(qū)坡面集中流曼寧系數(shù)n影響顯著。沖刷過程中，n在試驗前6 min內(nèi)迅速增大，后隨著沖刷歷時的延長呈緩慢增長或趨于穩(wěn)定的趨勢。n隨流量的增大而減小，不同巖石與坡面夾角的n值大小順序為：n180°>n150°>n90°>n30°>n120°>n60°。試驗條件下，n變化范圍為0.036～0.116。

(2)n與水深h關(guān)系受巖石與坡面夾角影響，當(dāng)巖石與坡面夾角大于30°時，n與水深h顯著相關(guān)，夾角為30°時，二者關(guān)系不顯著；試驗條件下，n與Re相關(guān)性不顯著，與Fr則呈顯著冪函數(shù)負(fù)相關(guān)。

(3) 相對水深曼寧系數(shù)n/h隨巖石與坡面夾角增大呈先增后減的趨勢，當(dāng)巖石與坡面夾角大于等于90°時存在臨界角度，臨界角度在90°～120°之間。巖石與坡面夾角大于30°時，n/h與Re均呈冪函數(shù)負(fù)相關(guān)，夾角為30°時，二者關(guān)系不顯著；不同巖石與坡面夾角的n/h與Fr關(guān)系均不顯著?？λ固夭酃葏^(qū)巖石與坡面夾角及地下裂隙結(jié)構(gòu)復(fù)雜多變，室內(nèi)模擬試驗條件與野外實際情況可能有所偏差，同時野外條件下的土壤異質(zhì)性及地表植被等因素均會對集中水流曼寧系數(shù)產(chǎn)生影響，因此還需結(jié)合野外喀斯特槽谷區(qū)徑流小區(qū)試驗對本試驗結(jié)果進(jìn)行進(jìn)一步驗證。