摘要：為了解坡面薄層水流在坡度、單寬流量及水流流床等外界條件發(fā)生變化時，其水流水動力特證參數(shù)的變化規(guī)律，從水力學(xué)的角度出發(fā)，在室內(nèi)開展阻力試驗，分別對3°、6°、9°、12°、15°等5種不同坡度、不同級配人工沙粒鋪設(shè)的粗糙床面展開水流水力學(xué)試驗研究，并對試驗得出的有關(guān)數(shù)據(jù)進行分析，初步探明了坡面水流流態(tài)為過渡流區(qū)的延伸，得出坡面薄層水流阻力系數(shù)、雷諾數(shù)、佛汝德數(shù)與單寬流量的變化規(guī)律及相互之間的關(guān)系。

關(guān)鍵詞：坡面流；流態(tài)；雷諾數(shù)；佛汝德數(shù)

中圖分類號：TV698.1 文獻標(biāo)識碼：A 文章編號：0439-8114（2013）22-5455-04

坡面水流主要指的是大氣降水或冰山融雪后產(chǎn)生的水流在自身重力的作用下沿著坡面向下運動的薄層水流，是形成河道水流的主要部分，有時也稱為片流或漫流。坡面水流的流態(tài)和徑流常被用來表示薄層坡面水流的動力學(xué)特征[1]。由于坡面水流的流態(tài)歸屬于坡面薄層水流的研究領(lǐng)域，因此與土壤、植被環(huán)境、降雨及其沖刷的坡溝細(xì)度、密度等參數(shù)有著十分密切的關(guān)系[2-4]，其中的相關(guān)參數(shù)不同會使得坡面薄層水流的流態(tài)發(fā)生較大改變。

一些研究人員認(rèn)為坡面地表的徑流既是地表徑流的匯集演進過程，也決定坡面土壤侵蝕發(fā)生的過程和侵蝕量[5]。當(dāng)水流侵蝕能力超過土壤抗蝕性后就可能發(fā)生跌坎并產(chǎn)生細(xì)溝，形成坡面細(xì)溝流，最終可出現(xiàn)淺溝和切溝流[6]。坡面薄層水流所流過的土壤表面的粗糙度對坡面水流流態(tài)變化也會有所影響，且這種影響引起的水流動力學(xué)相關(guān)的特征變化也比較明顯?，F(xiàn)在也用糙率表征水流床面的粗糙度，它也是坡面水流水動力學(xué)特性的直接影響因素，糙率又稱為有效糙率[7，8]。

當(dāng)前大部分研究是針對坡面水流在沒有植被覆蓋的裸土表面上展開[3]，隨著坡面粗糙程度不同引起的水流水力學(xué)特性的研究相對較少。很多研究者認(rèn)為粗糙床面主要是通過增加水流流過的阻力使得包含薄面水流的近壁水流流區(qū)發(fā)生變化[9]。本研究通過對不同粗糙流床展開試驗，對試驗數(shù)據(jù)進行分析，從而得出粗糙度不同的土壤坡面水流的水流特性，以期為坡面水流的侵蝕研究以及坡面水流水動力學(xué)研究提供理論基礎(chǔ)。

1 試驗方法

1.1 試驗設(shè)計

實際中坡面上的水流大多數(shù)是流床發(fā)生變化的非恒定水流，其流動過程中的有關(guān)流態(tài)參數(shù)也不斷發(fā)生變化，要直接進行比較全面的觀測非常困難，因此可以對其進行定量描述。本試驗將坡面薄層水流看作二維水流，采用專門的坡面流設(shè)計定床阻力試驗，沿水力方向設(shè)置5個觀測斷面，對水深和流速進行觀測記錄。

1）試驗用的水槽為5 mm厚的有機玻璃制作，長600 cm，寬60 cm，深20 cm，其底板的坡度可以進行調(diào)節(jié)，變化范圍在0°～15°。

在供水箱設(shè)置閥門和試驗水槽連接，加上自來水的閥門，共設(shè)置5個調(diào)流閥門，穩(wěn)流板2道，讓水流比較平緩地流進試驗水槽。考慮到薄面水流的水力學(xué)特性受邊界條件影響較大，為使得試驗數(shù)據(jù)更為精確，需要對出口的流速進行控制，保證水流平緩穩(wěn)定，因此在出口水箱中設(shè)置穩(wěn)流板和溢流板。

2）為保證試驗的準(zhǔn)確性，試驗在人工粗糙流床進行，選取坡度為3°、6°、9°、12°、15°[10]，同時在流床底部用油漆黏上自然級配粒徑均勻的沙粒，形成沙粒糙度。

3）由于研究坡面薄層水流的主要應(yīng)用是為水土流失提供理論支持，因此流量不應(yīng)該超出實際降雨量的范疇。為了能達到試驗的要求，試驗水流的單寬流量分別設(shè)置為5.0、7.5、10.0、12.5、15.0、20.0、25.0、30.0、40.0、60.0、90.0 L/min。

4）鑒于進行試驗的水槽中的水流比較深，使得水流的水動力學(xué)特征參數(shù)很容易因為外界條件的變化而發(fā)生變化，因此在試驗水槽的制作以及后期的安放等環(huán)節(jié)的質(zhì)量都做出了相對比較規(guī)范的操作要求，主要是為了保證所有的試驗材料為整體裝置，避免接縫對水流的干擾，使床面均勻平滑，無明顯局部擾動產(chǎn)生。

1.2 測定方法

在試驗過程中，要在水槽的坡面水流達到平穩(wěn)后反復(fù)準(zhǔn)確地測量單位時間內(nèi)的單寬流量、流速、水深，同時注意氣泡對試驗的影響。在試驗中，用染色劑法[11]分別在水槽進口1 m以下和水槽末端0.3～5.0 m的位置觀測流速，水深測量在粗糙床面采用自上而下的方法，水深計算方法見公式（1）。

h=-H′+H0 （1）

式（1）中，h為實際中采用的有效水深（mm）；H′是測針的讀數(shù)（mm）；H0是試驗采用的有機玻璃的厚度，本研究定位為5 mm。

本研究中量測流速采用染色劑法，測定的區(qū)域位于距水槽末端0.6 m以上的2.0～3.0 m范圍以內(nèi)；水層較深（>2.0 cm）時采用染色劑法，即在測定時將KMnO4溶液用塑料管注入水流中，記錄染色水流流經(jīng)測定區(qū)域所用的時間，再用測定區(qū)域長度除以時間即得水流表面速度[11，12]。水深使用精度為0.01 mm的重慶水文儀器制造廠生產(chǎn)的SX40-1型水位測針測定。

2 結(jié)果與分析

2.1 水流的平均流速

通常，謝才公式用來表征水流的層流，而曼寧公式用來表征水流的紊流，它們都把流速看成是水深與坡度之間的冪函數(shù)[13]?？墒呛髞淼难芯孔C明，在水流流動形成的細(xì)小的溝中，水流的速度與坡度幾乎沒有關(guān)系，只與流量之間形成函數(shù)關(guān)系，這種情況可以看成是由于水流流過的細(xì)小溝中的阻力引起的[14，15]。隨著坡面坡度不斷變大，水流流速也相應(yīng)增加，具有的能量也不斷增加，使得水流的侵蝕能力和流床的阻力也相應(yīng)變大。當(dāng)流量變大時，水流的平均流速以冪函數(shù)的形式變大；在流量較小時，坡面坡度對水流的流速產(chǎn)生的影響比較小，隨著流量的不斷變大，坡面坡度對水流流量的影響也越來越大。在坡度變大的時候，水流平均流速相應(yīng)變大，但流量不同時，坡度引起的水流流速的變化也不相同，在流量比較大的時候，坡度的影響也相應(yīng)較大。

由于坡面水流流態(tài)確定比較困難，使得其水力計算也較為困難。1945年，Horton[16]假定坡面流是紊流與層流流態(tài)相互交錯的混合流，從而可以認(rèn)為：

式（2）中，q為單寬流量；Kf為綜合系數(shù)，表示坡面阻力；S0為坡面坡度；h為有效水深；M1、M2為水流流區(qū)決定的指數(shù)，一般情況下，層流情況下M2為3，紊流情況下M2為5/3。式中，單寬流量、坡度、有效水深均為可測量要素，因此以指數(shù)M1為主要研究對象，并以此來確定水流流態(tài)。

通過對試驗數(shù)據(jù)進行分析，對水流的單寬流量（q）與平均流速（V）進行擬合，得出水流的平均流速和單寬流量成冪函數(shù)關(guān)系擬合情況最好，其擬合方程見表1。

2.2 水流的平均水深

通過試驗，得出不同坡度條件下平均水流深度隨單寬流量的變化規(guī)律。由表2可知，當(dāng)單寬流量逐漸增大的時候，坡面水流的平均水深也逐漸增大。但是，隨著坡度變大，水深卻變小。在水流流量比較小的時候，平均水深的差異相對很小，而在流量變大的情況下坡度的影響也隨之增加。另外還發(fā)現(xiàn)當(dāng)坡度由小變大的過程中，平均水深之間的變化卻逐漸變小。

對水流單寬流量（q）與平均水深（H）進行擬合得到擬合方程，具體見表3。由表3可以看出，單寬流量和平均水深之間的冪函數(shù)關(guān)系擬合較好。

2.3 傳統(tǒng)雷諾數(shù)的變化

采用水力學(xué)的計算方法，通過公式（3）對水流流態(tài)雷諾數(shù)（Re）確定坡面流的流態(tài)，使用Darcy-Weisbach公式計算阻力系數(shù)，用公式（4）表示。水流的粘滯系數(shù)依據(jù)測定的水溫，通過公式（5）計算。

Re=hV/v （3）

式（3）中，Re表示二元流雷諾數(shù)；V表示水流平均流速（m/s）；v表示運動粘滯系數(shù)（m2/s）；h表示有效水深（m），把坡面水流看成二元流，水力半徑數(shù)值上等于有效水深。

？姿=8ghJ/V2 （4）

式（4）中，？姿表示阻力系數(shù)；g為重力加速度（g=9.8 m/s2）；J為水力坡度，對于均勻流，J=i，i為水槽底坡度。

ν=0.017 75/（1+0.033 7t+0.000 22t2） （5）

式（5）中，ν表示粘滯系數(shù)；t為水溫（℃）。

通過對人工粗糙床面情況下試驗數(shù)據(jù)的分析，得到3°、6°、9°、12°、15°不同坡度下不同單寬流量與相應(yīng)雷諾數(shù)的關(guān)系。

對試驗數(shù)據(jù)的分析得出，在流床相同以及坡度不變的情況下，流量增加，雷諾數(shù)也增加，且雷諾數(shù)在單寬流量較大時的變化更為明顯。

經(jīng)過對相關(guān)試驗數(shù)據(jù)分析，擬合得到雷諾數(shù)與單寬流量之間的關(guān)系式，通過對試驗數(shù)據(jù)的分析，得到雷諾數(shù)與單寬流量之間的冪函數(shù)關(guān)系擬合最好，其擬合方程見表4。

2.4 阻力系數(shù)

研究表明，坡面水流也存在著紊流、層流和過渡流流態(tài)的區(qū)別。在此次試驗中雷諾數(shù)與坡面流阻力在降雨均勻及坡度恒定的條件下也有著比較明顯的分區(qū)規(guī)律。通過對本試驗數(shù)據(jù)進行分析，得到了阻力系數(shù)λ的數(shù)值，具體見表5。

2.5 坡面水流水動力學(xué)特性參數(shù)的關(guān)系

2.5.1 阻力系數(shù)與雷諾數(shù)的關(guān)系 通過對5個不同坡度的11個不同單寬流量的55組試驗數(shù)據(jù)進行分析研究發(fā)現(xiàn)，雷諾數(shù)增大，阻力系數(shù)減小。其線性方程為公式（6）。

λ=24/Re （6）

依據(jù)普遍的水力學(xué)雷諾數(shù)層流的界限值，試驗中表現(xiàn)出了層流特性，即λ與Re之間的關(guān)系在進入層流區(qū)后呈線性關(guān)系。利用傳統(tǒng)雷諾數(shù)判斷出3°、6°、9°、12°、15°時點據(jù)都在直線上方，而且點據(jù)大部分位于過渡區(qū)，只有少部分位于層流區(qū)和紊流區(qū)。

2.5.2 佛汝德數(shù)與單寬流量的關(guān)系 佛汝得數(shù)（Fr）為流速與波速的比，常用公式（7）計算。

Fr=■=■ （7）

式（7）中，w為波速（m/s）。

通過試驗得到不同單寬流量下的佛汝德數(shù)，具體見表6。由表6可知，在相同床面狀況、等坡度情況下，當(dāng)單寬流量變大時，佛汝德數(shù)隨之增加，表示隨著單寬流量的增加，水流流態(tài)也在發(fā)生改變。

3 結(jié)論

1）根據(jù)坡面水流的形成起因，依據(jù)水力學(xué)原理，結(jié)合當(dāng)前的試驗條件，設(shè)計了當(dāng)坡度及床面等外界條件發(fā)生變化的情況下，坡面水流侵蝕能力的試驗方法。

2）通過對試驗數(shù)據(jù)的分析，探明坡面流通常被認(rèn)為是在層流區(qū)、過渡流區(qū)和阻力平方區(qū)的三區(qū)和過渡流區(qū)的延伸。

3）試驗中發(fā)現(xiàn)，坡面水流基本表現(xiàn)出層流的特性，阻力系數(shù)與傳統(tǒng)水力學(xué)雷諾數(shù)之間的關(guān)系在進入層流區(qū)后表現(xiàn)為線性關(guān)系。

4）在床面相同、坡度不變的條件下，佛汝德數(shù)隨單寬流量的增加而增加，表明單寬流量增加時，水流流態(tài)同時也在發(fā)生著變化。

參考文獻：

[1] 王玉強，王英華，張寬地.基于模式搜索算法的Van Genuchten方程參數(shù)優(yōu)化[J].中國農(nóng)村水利水電，2010（9）：8-10，14.

[2] 潘成忠，上官周平.牧草對坡面侵蝕動力參數(shù)的影響[J].水利學(xué)報，2005，36（3）：371-377.

[3] 郭雨華，趙廷寧，孫保平，等.草地坡面水動力學(xué)特性及其阻延地表徑流機制研究[J].水土保持研究，2006，13（4）：264-267.

[4] 沙際德，白清俊.粘性土坡面細(xì)溝流的水力特性試驗研究[J].泥沙研究，2001（6）：39-44.

[5] 姚文藝.坡面流阻力規(guī)律實驗研究[J].泥沙研究，1996（1）：74-81.

[6] 胡世雄，靳長興.坡面流與坡面侵蝕動力過程研究的最新進展[J].地理研究，1998，17（3）：326-335.

[7] 沈 冰， 李懷恩 沈 晉，等.坡面降雨漫流過程中有效糙率的實驗研究[J].水利學(xué)報，1994（10）：61-68.

[8] ENGMAN E T. Roughness coefficient for routing surface runoff[J]. Journal of Irrigation and Drainage Engineering，1986，112（1）：39-54.

[9] 劉春晶，李丹勛，曲兆松，等.過渡區(qū)動床明渠流的流速分布[J].水科學(xué)進展，2006，17（1）：49-54.

[10] 敬向鋒，呂宏興，潘成忠，等.侵蝕性坡面流流態(tài)的試驗研究[J].灌溉排水學(xué)報，2008，27（2）：82-85.

[11] 敬向鋒，呂宏興，潘成忠，等.坡面薄層水流流態(tài)判定方法的初步探討[J].農(nóng)業(yè)工程學(xué)報，2007，23（5）：56-61.

[12] 敬向鋒.坡面薄層水流水動力學(xué)規(guī)律的試驗研究[D].陜西楊凌：西北農(nóng)林科技大學(xué)，2007.

[13] 吳持恭.水力學(xué)[M].第二版.北京：高等教育出版社，1998.

[14] GOVERS G. Relationships between discharge， velocity， and flow area for rills eroding loose， non-layered materials[J].Earth Surface Processes and Landforms，1992，17（5）：515-528.

[15] NEARING M， SIMANTON R， NORTON D， et al. Soil erosion by surface water flow on a stony， semiarid hillslope[J].Earth Surface Processes and Landforms，1999，24（8）：677-686.

[16] HORTON R E. Erosional development of streams and their drainage basin：Hydrophysical approach to quantitative morphology[J]. Bull Geol Soc AM，1945，56：275-370.

（責(zé)任編輯 呂海霞）