張升堂, 梁 博, 張 楷

(1.山東科技大學(xué) 地球科學(xué)與工程學(xué)院, 山東 青島 266590; 2.陜西省水文水資源勘測(cè)局, 陜西 西安710068)

植被分布對(duì)地表糙率影響

張升堂1, 梁 博1, 張 楷2

(1.山東科技大學(xué) 地球科學(xué)與工程學(xué)院, 山東 青島 266590; 2.陜西省水文水資源勘測(cè)局, 陜西 西安710068)

[目的] 分析在植被淹沒(méi)或非淹沒(méi)狀態(tài)不同植被覆蓋密度影響下，地表糙率的變化規(guī)律與特征，為進(jìn)一步研究植被對(duì)地表糙率的影響提供參考。[方法] 試驗(yàn)以塑料棒模擬植被分布，模擬出不同密度，分別進(jìn)行放水試驗(yàn)。[結(jié)果] 在植被淹沒(méi)與非淹沒(méi)狀態(tài)下地表糙率的變化不同，在非淹沒(méi)狀態(tài)時(shí)，地表糙率隨著平均水深的增加呈減小的趨勢(shì)，在淹沒(méi)水深小時(shí)，地表糙率隨著平均水深的增加而先增大后減小，而當(dāng)淹沒(méi)水深很大時(shí)，地表糙率趨向于一個(gè)穩(wěn)定值。[結(jié)論] 地表糙率取值隨植被覆蓋密度的增大而增大；同一下墊面情況下，不同的水流方向，地表糙率取值不相同。

地表糙率; 植被密度; 淹沒(méi)植被; 非淹沒(méi)植被

地表糙率對(duì)地表徑流的流速、滲透、土壤侵蝕等都具有很大的影響作用，研究地表糙率的變化規(guī)律對(duì)于估算洪峰流量以及水土保持措施決策等具有重要意義。

地表糙率的影響因素很多，降雨、土壤類型、植被覆蓋等因素都對(duì)地表糙率有很大的影響，因此地表糙率的取值復(fù)雜多變，難以確定。

在現(xiàn)實(shí)流域下墊面經(jīng)常有雜草、灌木等，或耕種的作物、草被、樹(shù)木等植被覆蓋于地表，這些植被其形狀、剛?cè)嵝?、密度、以及分布形式等在很大程度上影響了水流的水力學(xué)參數(shù)，增大了水流阻力，加大了地表糙率，降低了地表過(guò)流能力，同時(shí)也起到了防止水土流失的作用。

因此研究植被分布對(duì)地表糙率的影響對(duì)于探討地表徑流對(duì)地表沖刷能力、坡面水沙運(yùn)動(dòng)規(guī)律以及調(diào)節(jié)洪水、削減洪峰等具有重要意義。所以植被分布對(duì)地表糙率影響問(wèn)題的研究具有重要的理論價(jià)值以及實(shí)際應(yīng)用價(jià)值。

近幾十年來(lái)，國(guó)內(nèi)外諸多學(xué)者在植被對(duì)地表糙率的影響方面進(jìn)行了研究[1-5]。

Hsieh T.[6]和Li R. M.等[7]，分別用圓柱體模擬植被，研究矩形渠道中其對(duì)水流的影響過(guò)程，指出樹(shù)的密度與樹(shù)的排列方式對(duì)水流速度的影響是不同的，水流阻力也會(huì)隨圓柱布置密度的增大而增大。Dunn C.等[8]在矩形渠道中，對(duì)淹沒(méi)狀態(tài)下的剛性植物以及柔性圓棒進(jìn)行試驗(yàn)，發(fā)現(xiàn)水流阻力也是隨植物密度的增加而增加。閆旭峰等[9]通過(guò)試驗(yàn)研究得出：由于植被阻水及水流混摻效應(yīng)，Darcy—Weisbach阻力系數(shù)基本隨雷諾數(shù)增大而增大，且植被剛度及密度越大，增大趨勢(shì)越顯著。Kim J.等[10]通過(guò)植被部分淹沒(méi)研究坡面流曼寧糙率系數(shù)變化，指出坡面流糙率隨植被部分淹沒(méi)情況變化而變化，將顯著影響坡面流路徑和坡面侵蝕狀況。拉尼亞[11]通過(guò)對(duì)人工模擬的菖蒲進(jìn)行水槽試驗(yàn)，研究了4種分布密度對(duì)水流阻力及水流結(jié)構(gòu)的影響。研究表明，植被密度的增大將導(dǎo)致過(guò)水?dāng)嗝娴脑龃蠛妥枇Φ脑黾印?/p>

黃文沛[12]通過(guò)室內(nèi)水槽分別進(jìn)行了植被密度試驗(yàn)與植被格局試驗(yàn)，其中格局試驗(yàn)主要為改變植被段長(zhǎng)度；分析得到隨著植被密度的增大，曼寧糙率系數(shù)明顯增大，植被段長(zhǎng)度的變化對(duì)曼寧糙率系數(shù)的影響不明顯。姬昌輝等[13]采用概化水槽試驗(yàn)研究了在不同植物間距、水流條件下，含淹沒(méi)植被明渠水位、曼寧糙率系數(shù)的變化特征，試驗(yàn)結(jié)果表明水流條件相同的情況下，隨著植物排列間距的減小，糙率系數(shù)明顯增大。房春艷等[14]通過(guò)室內(nèi)變坡水槽試驗(yàn)研究了復(fù)式河槽中灘地植被對(duì)水流阻力的影響，認(rèn)為影響灘地曼寧糙率系數(shù)的主要因素為植被密度、水深、淹沒(méi)與非淹沒(méi)植被。

雖然目前關(guān)于植被分布對(duì)地表糙率影響的研究已有很多，但多集中在覆蓋密度或植被淹沒(méi)方面。本文通過(guò)室內(nèi)水槽放水模擬試驗(yàn)，研究分析在植被淹沒(méi)與非淹沒(méi)狀態(tài)下，不同植被覆蓋密度影響下，地表糙率的變化規(guī)律與特征，并且首次通過(guò)試驗(yàn)分析相同的植被排列方式下不同水流方向?qū)Φ乇聿诼实挠绊?，以期為進(jìn)一步研究植被對(duì)地表糙率的影響提供參考。

1 試驗(yàn)設(shè)計(jì)與試驗(yàn)數(shù)據(jù)

1.1 試驗(yàn)裝置

試驗(yàn)裝置主要由水箱、矩形水槽、量水堰等組成。采用底寬0.6 m，側(cè)壁高0.9 m，長(zhǎng)度13 m的矩形水槽，水槽分上游平水段、試驗(yàn)鋪設(shè)段、下游量水段。下游量水段末端設(shè)有矩形薄壁量水堰觀測(cè)流量。水槽由水箱定水頭供水，在水槽與水箱連接處設(shè)有閥門控制流量。在矩形水槽末端設(shè)有回流池，水流流經(jīng)水槽后進(jìn)入回流池，然后由水泵將回流池中的水抽入水箱中，用以實(shí)現(xiàn)水流循環(huán)，節(jié)約用水。

1.2 試驗(yàn)下墊面條件

天然植物不管是柔性的還是剛性的，由于其柔韌度、植株和植被葉子形狀各異等植被生物特性，導(dǎo)致對(duì)直接模擬天然植被的特性十分困難，而且不易固定，不易用來(lái)在變坡水槽中進(jìn)行試驗(yàn)測(cè)量，目前多數(shù)這方面的研究都是將植被進(jìn)行簡(jiǎn)化[15-19]。因此本試驗(yàn)概化地表植被類型，以塑料棒模擬地表植被，塑料棒高度為0.15 m。為了可以更好地分析植被不同分布方式對(duì)地表糙率的影響，試驗(yàn)設(shè)計(jì)了4種不同的玻璃底板用以模擬植被的不同分布。取2塊有機(jī)玻璃板，分別在玻璃板上按與水流方向呈90°和45°夾角有規(guī)律的鉆孔，孔間縱、橫間距均為60 mm，模擬植物以60 mm×60 mm和60 mm×120 mm這2種固定間距“種植”于有機(jī)玻璃板孔內(nèi)，模擬出不同的植被覆蓋密度。

2 結(jié)果與分析

圖1給出了4種不同玻璃底板類型下，平均水深h與地表糙率n的關(guān)系圖。試驗(yàn)結(jié)果表明，隨著平均水深的增大，地表糙率取值呈減小的趨勢(shì)，并且在平均水深較小時(shí)，地表糙率取值變化幅度較大；但隨著平均水深逐漸增大，地表糙率取值變化幅度越來(lái)越小，并逐漸趨向一個(gè)穩(wěn)定值。由此可見(jiàn)，在植被覆蓋下水深對(duì)地表糙率有很大的影響作用。這與李勉等[20]通過(guò)放水沖刷試驗(yàn)所得結(jié)論相似。

圖1 4種不同底板類型條件下平均水深與糙率系數(shù)的關(guān)系

2.1 植被淹沒(méi)狀態(tài)對(duì)地表糙率的影響

根據(jù)圖1可以看出，在平均水深小于0.15 m時(shí)，地表糙率取值變化幅度很大，隨著平均水深的逐漸增大，其取值變化幅度逐漸減小。

由于模擬植被高度為0.15 m，因此將試驗(yàn)過(guò)程按水深分為3個(gè)區(qū)域，第1區(qū)域?yàn)槠骄钚∮?.15 m時(shí)，第2區(qū)域?yàn)槠骄钐幱?.15～0.20 m時(shí)，第3區(qū)域?yàn)槠骄畲笥?.20 m時(shí)。

由圖1可知，在第1區(qū)域，水流平均水深小于植被高度，為非淹沒(méi)狀態(tài)，此時(shí)平均水深很小，植被對(duì)地表水流的影響作用很大，隨著平均水深的增大，地表糙率取值呈減小趨勢(shì)，并且變化幅度較大。而在第2區(qū)域，植被處于淹沒(méi)狀態(tài)，此時(shí)植被剛剛被淹沒(méi)，但是淹沒(méi)水深很小。通過(guò)試驗(yàn)現(xiàn)場(chǎng)觀察發(fā)現(xiàn)當(dāng)水流剛剛淹沒(méi)模擬植被時(shí)，隨著水流流動(dòng)植被頂端會(huì)出現(xiàn)高頻率的擺動(dòng)，這將對(duì)地表糙率取值產(chǎn)生影響，因此地表糙率的取值規(guī)律也更加復(fù)雜，通過(guò)圖1分析得出地表糙率系數(shù)在第2區(qū)域基本呈先增大后減小的趨勢(shì)。對(duì)于第3區(qū)域，水流平均水深比植被高度高很多，處于完全淹沒(méi)狀態(tài)，并且淹沒(méi)水深逐漸增大，此時(shí)植被頂端不再擺動(dòng)，而是根據(jù)水流流向彎折。隨著淹沒(méi)水深的增大，淹沒(méi)植被對(duì)水流的影響逐漸減小，甚至可以忽略不計(jì)，因此地表糙率系數(shù)趨向于一個(gè)穩(wěn)定值。Chow V. T.[21]指出只有水位低于某一特定水位，植被才對(duì)水流有顯著影響，因此當(dāng)確定一個(gè)比河道中植被還要高得多的洪水位時(shí)，可以用一個(gè)固定的糙率系數(shù)。

因此可以得出，植被在非淹沒(méi)和淹沒(méi)狀態(tài)下，地表糙率取值變化規(guī)律并不相同，但總的來(lái)說(shuō)，非淹沒(méi)狀態(tài)時(shí)植被對(duì)地表糙率取值影響較大，隨著淹沒(méi)水深的增加，植被對(duì)地表糙率取值影響逐漸減小甚至可以忽略不計(jì)。

2.2 植被覆蓋對(duì)地表糙率的影響

2.2.1 不同覆蓋密度對(duì)地表糙率取值影響 植被覆蓋密度是影響地表糙率取值的重要因素之一。本試驗(yàn)設(shè)計(jì)了不同水流方向下2種不同植被覆蓋密度下墊面，根據(jù)試驗(yàn)數(shù)據(jù)做出了不同植被覆蓋密度下水深與地表糙率的關(guān)系圖(圖2—3)。

圖2 θ為90°，45°條件下不同密度水深與地表糙率的關(guān)系

圖3 60 mm×60 mm和60 mm×120 mm密度，不同水流方向條件下地表糙率與水深關(guān)系

由圖2可以看出，不同的植被覆蓋密度下，其地表糙率取值是不相同的。在θ為90°時(shí)，從整體上可以看出，在平均水深較小時(shí)，不同植被覆蓋密度下地表糙率取值并不具有明顯規(guī)律性，只有在平均水深較大時(shí)，地表糙率取值才隨著植被覆蓋密度增加而增大，在前3個(gè)測(cè)量點(diǎn)時(shí)，植被覆蓋密度大時(shí)，地表糙率取值并不一定大。但是當(dāng)θ為45°時(shí)，可以看出隨著試驗(yàn)的進(jìn)行，不同植被覆蓋密度下地表糙率取值有著明確的規(guī)律性，即植被覆蓋密度大時(shí)，地表糙率取值大。并且在θ為90°時(shí)，只有前3個(gè)測(cè)量點(diǎn)不具有規(guī)律性，從第3個(gè)測(cè)量點(diǎn)之后也是呈植被覆蓋密度大，地表糙率取值大的規(guī)律。由于試驗(yàn)條件等原因，并沒(méi)有進(jìn)行重復(fù)試驗(yàn)，因此，圖2中的前3個(gè)測(cè)量點(diǎn)可能是在試驗(yàn)過(guò)程中或者測(cè)量時(shí)產(chǎn)生了誤差，其具體原因還有待進(jìn)一步分析。但是總的來(lái)說(shuō)，通過(guò)分析圖2可以看出，植被覆蓋密度對(duì)地表糙率取值有很大的影響，不同的植被覆蓋密度下地表糙率取值基本呈覆蓋密度大、地表糙率取值大的變化規(guī)律。這與肖培青等[22]、潘成忠等[23]、田風(fēng)霞等[24]所得出的結(jié)論一致。

2.2.2 相同覆蓋密度、不同水流方向?qū)Φ乇聿诼实挠绊?影響地表糙率取值變化的因素很多，在相同的植被覆蓋密度下，不同的水流方向?qū)Φ乇聿诼嗜≈狄灿泻艽蟮挠绊懽饔茫鶕?jù)試驗(yàn)數(shù)據(jù)做出了相同密度、不同水流方向下，地表糙率與平均水深關(guān)系圖(圖3)。

在現(xiàn)實(shí)流域面，坡面水流一般漫流于坡面，且受坡面地表的影響使其流動(dòng)方向多向。目前關(guān)于地表糙率的研究一般認(rèn)為在同一地表情況下其地表糙率都取同一定值。Engman E.T.[25]曾由野外小區(qū)實(shí)測(cè)洪水資料，經(jīng)迭代適線方法得出了不同地表情況下的Manning糙率系數(shù)，例如遭受侵蝕的裸露黏壤地表上n=0.02，天然牧場(chǎng)n=0.13，刈割放牧地n=0.10，藍(lán)草草地n=0.45。但是張升堂等[26]首次提出了“矢量糙率理論”，即糙率的取值應(yīng)該是隨匯流方向的不同而不同。由圖3可以看出，在相同的植被覆蓋密度、相同的平均水深下，水流流向不同其地表糙率取值是不相同的。也就是說(shuō)，在同一下墊面情況下，坡面水流由于流向的不同而具有不同的地表糙率取值，從而進(jìn)一步驗(yàn)證了張升堂所提出的“矢量糙率理論”。

3 結(jié) 論

(1) 在植被非淹沒(méi)狀態(tài)下，地表糙率隨著平均水深的增加而減??；在淹沒(méi)狀態(tài)下，當(dāng)淹沒(méi)水深很小時(shí)，地表糙率隨著平均水深的增加而先增大后減??；隨著淹沒(méi)水深越來(lái)越大，地表糙率趨向于一個(gè)穩(wěn)定值。

(2) 植被覆蓋密度對(duì)地表糙率取值有很大的影響，不同的植被覆蓋密度下地表糙率取值基本呈覆蓋密度大、地表糙率取值大的變化規(guī)律。

(3) 影響地表糙率取值變化的因素很多，在相同的植被覆蓋密度下，不同的水流方向?qū)Φ乇聿诼嗜≈狄灿泻艽蟮挠绊懽饔?。即在同一下墊面情況下，坡面水流由于流向的不同而具有不同的地表糙率取值。

[1] Wu Fuchun, Shen Wenhseis, Chou Yiju. Variation of roughness coefficients for unsubmerged and submerged vegetation[J]. Journal of Hydraulic Engineering, 1999,125(9):934-942.

[2] Stone B M, Shen H T. Hydraulic resistance of flow in channels with cylindrical roughness[J]. Journal of Hydraulic Engineering, 2002,128(5):500-506.

[3] James C S, Birkhead A L, Jordanova A A, et al. Flow resistance of emergent vegetation [J]. Journal of Hydraulic Research, 2004,42(4):390-398.

[4] 拾兵,王川源,尹則高,等.淹沒(méi)植物對(duì)河道糙率的影響[J].中國(guó)海洋大學(xué)學(xué)報(bào),2009,39(2):295-298.

[5] 曾玉紅,槐文信,張健,等.非淹沒(méi)剛性植被流動(dòng)阻力研究[J].水利學(xué)報(bào),2011,42(7):834-838.

[6] Hsieh T. Resistance of cylinder piers in open-channel flow[J]. Journal of Hydraulics Division, 1964,90(1):161-173.

[7] Li R M, Shen H W. Effect of tall vegetation on flow and sediment[J]. Journal of Hydraulics Division, 1973,99(5):793-814.

[8] Dunn C, Lopez F, Garcia M. Mean flow and turbulence in a laboratory channel with simulated vegetation[M]. Hydrosystems Laboratory, Department of Civil Engineering, University of Lllinois at Urbana-Champaign,1996.

[9] 閆旭峰,周蘇芬,黃爾,等.植被條件下坡面薄層水流動(dòng)力學(xué)特性試驗(yàn)研究[J].四川大學(xué)學(xué)報(bào):工程科學(xué)版,2012,44(2):26-30.

[10] Kim J, Ivanov V Y, Katopodes N D, Hydraulic resistance to overland flow on surfaces with partially submerged vegetation[J], Water Resources Research, 2012,48(10):1-19.

[11] 拉尼亞.非淹沒(méi)植被水流特性研究[D].南京:河海大學(xué),2005.

[12] 黃文沛.非淹沒(méi)剛性植被對(duì)水流阻力和水流能量特性影響的研究[D].北京：華北電力大學(xué),2013.

[13] 姬昌輝,洪大林,丁瑞,等.含淹沒(méi)植被明渠水位及糙率變化試驗(yàn)研究[J].水利水運(yùn)工程學(xué)報(bào),2013,(1):60-65.

[14] 房春艷,羅憲.灘地植被化復(fù)式河槽的水流阻力特性試驗(yàn)[J].重慶交通大學(xué)學(xué)報(bào)：自然科學(xué)版,2013,32(4):668-672.

[15] 閆靜,唐洪武,田志軍,等.植物對(duì)明渠流速分布影響的試驗(yàn)研究[J].水利水運(yùn)工程學(xué)報(bào),2011(4):138-142.

[16] Noarayanan L, Murali K, Sundar V. Manning’s ‘n’ co-efficient for flexible emergent vegetation in tandem configuration[J]. Journal of Hydro-environment Research, 2012,6(1):51-62.

[17] Tanaka N, Yagisawa J. Flow structures and sedimentation characteristics around clump-type vegetation[J]. Journal of Hydro-environment Research, 2010,4(1):15-25.

[18] 倪漢根,顧峰峰.濕地非淹沒(méi)蘆葦水流阻力的試驗(yàn)研究[J].水動(dòng)力學(xué)研究與進(jìn)展:A輯,2005,20(2):167-173.

[19] 吳福生.含植物明渠水動(dòng)力特性研究[D].南京:南京水利科學(xué)研究院, 2009.

[20] 李勉,姚文藝,陳江南,等.草被覆蓋下坡面—溝坡系統(tǒng)坡面流阻力變化特征試驗(yàn)研究[J].水利學(xué)報(bào),2007,38(1):112-119.

[21] Chow V T. Open-channel hydraulics[M]. New York: McGraw-Hill, 1959.

[22] 肖培青,姚文藝,李莉,等.植被影響下坡面流阻力變化特征研究[J].泥沙研究,2013(3):1-5.

[23] 潘成忠,上官周平.降雨和坡度對(duì)坡面流水動(dòng)力學(xué)參數(shù)的影響[J].應(yīng)用基礎(chǔ)與工程科學(xué)報(bào),2009,17(6):843-851.

[24] 田風(fēng)霞,劉剛,鄭世清,等.草本植物對(duì)土質(zhì)路面徑流水動(dòng)力學(xué)特征及水沙過(guò)程的影響[J].農(nóng)業(yè)工程學(xué)報(bào),2009,25(10):25-29.

[25] Engman E T. Roughness coefficients for routing surface runoff[J]. Journal of Irrigation and Drainage Engineering, 1986,112(1):39-53.

[26] 張升堂,康紹忠.基于矢量糙率的柵格單元流量分配模型[J].水利學(xué)報(bào),2006,36(11):1326-1330.

Effects of Vegetation Distribution on Earth Surface Roughness

ZHANG Shengtang1, LIANG Bo1, ZHANG Kai2

(1.CollegeofGeologicalScience&Engineering,ShandongUniversityofScienceandTechnology,Qingdao,Shandong266590,China;2.ShaanxiSurveyBureauofHydrologyandWaterResources,Xi’an,Shaanxi710068,China)

[Objective] The change rules and characteristics of the surface roughness was analyzed under the station of submerged and unsubmerged and the different vegetation density in order to provide the reference for studying the effect of vegetation on the surface roughness. [Methods] The plastic rods, as a surrogate for vegetation distribution, were used to simulate two different densities in the flushing experiments. [Results] The earth surface roughness changes were different between the submerged vegetation and unsubmerged vegetation. Under the unsubmerged condition, the earth surface roughness decreased with the increase of the average water depth.While submerged depth was low, the earth surface roughness firstly increased and then decreased with the increase of the average water depth. When the submerged depth was very deep, the earth surface roughness tends to be a stable value. [Conclusion] The earth surface roughness value increases along with the vegetation density. Under the same underlying surface, the surface roughness value changes with different flow direction.

surface roughness; vegetation density; submerged vegetation; unsubmerged vegetation

2014-07-11

2014-09-16

國(guó)家自然科學(xué)基金項(xiàng)目“現(xiàn)代流域分布式坡面匯流模擬的矢量糙率理論” (40971021)， “考慮阻力作用的分布式坡面匯流多流向流量分配算法” (41471025)， 山東省自然科學(xué)基金項(xiàng)目“坡面匯流阻力空間變異及其影響研究”(ZR2014DM004)

張升堂(1970—),男(漢族),陜西省鳳翔縣人,博士(后),副教授,主要從事水文水資源方面研究。 E-mail:zst0077@163.com。

A

1000-288X(2015)05-0045-04

S715.3