呂蘭州，胡建朋 ，楊吉華

（1.山東省膠東調(diào)水局，山東 濟(jì)南 250100；2.山東農(nóng)業(yè)大學(xué)林學(xué)院，山東 泰安 271018）

山東省山丘區(qū)面積大，山高坡陡，土壤瘠薄，蓄水保土性能差，水土流失嚴(yán)重；土地生產(chǎn)力低而不穩(wěn)，災(zāi)害頻繁發(fā)生，嚴(yán)重制約了區(qū)域經(jīng)濟(jì)發(fā)展和生態(tài)環(huán)境的改善，水土流失問題成為影響生態(tài)、生產(chǎn)和生活安全的重要問題之一。研究?jī)蓚€(gè)不同治理效果的小流域所產(chǎn)生的徑流量差異，對(duì)小流域綜合治理具有重要意義。

1 試驗(yàn)區(qū)概況

曲阜市紅山水土保持試驗(yàn)站位于紅山小流域境內(nèi)，屬典型魯中南低山丘陵區(qū)，流域面積8.88 km2，海拔 548.1 m，溝壑密度 3.8 km/km2，地面坡度 3°～25°，巖石為花崗片麻巖，土壤為棕壤。多年平均氣溫13.7℃，年平均無(wú)霜期192 d，多年平均降水量671.6 mm，降水多集中在7、8、9月份，年平均徑流深246 mm。植被屬溫帶落葉闊葉林，流域內(nèi)喬木樹種以刺槐、黑松、麻櫟、側(cè)柏、毛白楊為主。灌木主要有黃荊、酸棗等。經(jīng)濟(jì)樹種有蘋果、棗、桃、梨和杏。自然生長(zhǎng)的雜草以黃草為主，人工種植的草本植物有小麥、谷子、地瓜、花生、沙打旺、草木樨等。

本試驗(yàn)選擇紅山水保試驗(yàn)站的九仙山、平坡小流域，這兩個(gè)流域除治理效果不同外，其他條件如地形、地質(zhì)、土壤、氣候、流域面積等基本相同或相似，對(duì)比小流域選擇在地面分水線和地下分水線基本重合的閉合流域，溝長(zhǎng)、比降比較接近。流域自然地理特征見表1。

2 試驗(yàn)方法

2.1 降雨量

1）降雨量觀測(cè)：在流域內(nèi)根據(jù)不同的方位和海拔高程布設(shè)雨量觀測(cè)點(diǎn)，采用人工雨量計(jì)、自記雨量記同步觀測(cè)、校準(zhǔn)，據(jù)此計(jì)算日降雨量和降雨強(qiáng)度。

表1 小流域自然地理特征對(duì)比表

2）降雨量和降雨強(qiáng)度分級(jí)：根據(jù)日降雨量（P），共分 6 個(gè)雨量級(jí)，即 0～9.9 mm、10～29.9 mm、30～49.9 mm、50～79.9 mm、80～99.9 mm、＞100 mm；根據(jù)計(jì)算的 30 min最大降雨強(qiáng)度（I30），共分 7 個(gè)雨強(qiáng)級(jí)，即 0～5 mm/h、5～10 mm/h、10～20 mm/h、20～30 mm/h、30～40 mm/h、40～50 mm/h、＞50 mm/h。

2.2 小流域徑流觀測(cè)

在流域出口控制斷面處修建三角形量水堰，小流域流量測(cè)流采用流速儀法，使用旋漿流速儀，水位用自記水位計(jì)，觀測(cè)精度可讀0.01 m，根據(jù)自記水位計(jì)記錄的水位變化曲線，求出各時(shí)段和日徑流量。

2.3 數(shù)據(jù)處理

對(duì)所得的數(shù)據(jù)運(yùn)用SPSS11.5軟件EXCEL2003軟件進(jìn)行統(tǒng)計(jì)分析，建立不同條件下的徑流量模型。

3 結(jié)果分析

3.1 不同治理效果小流域徑流總量變化規(guī)律

1）典型降雨條件下徑流量變化。采用4年的觀測(cè)資料（1982—1985年），從中選出22場(chǎng)典型降雨，在0～9.9、10～29.9、30～49.9、50～79.9、80～99.9、＞100 mm各雨量級(jí)均勻分配，在每個(gè)雨量級(jí)下求出22場(chǎng)典型降雨所產(chǎn)生的平均徑流量見表2。從不同雨量級(jí)看出九仙山小流域的徑流量明顯高于平坡小流域，說明治理效果好的小流域產(chǎn)生的徑流量遠(yuǎn)遠(yuǎn)小于治理差的。從不同治理效果的小流域看出，九仙山小流域在50～79.9 mm雨量段徑流量增長(zhǎng)迅速，但在80～99.9 mm雨量段增長(zhǎng)幅度達(dá)到最大，增長(zhǎng)了7 977 m3/km2，增長(zhǎng)率為244.7%。平坡小流域隨著降雨量的增加，徑流量增幅較小，但是在＞100 mm雨量段徑流量迅速增長(zhǎng)，增長(zhǎng)了10 250 m3/km2，增長(zhǎng)率為117.1%。

表2 不同雨量級(jí)下的平均徑流量 m3/km2

2）降雨因子與徑流量回歸方程。為建立在不同治理效果下的小流域降雨因子與徑流量的經(jīng)驗(yàn)?zāi)Ｐ?，?duì)降雨量和I30與徑流量進(jìn)行逐步回歸分析，建立多元線性回歸方程，從表3可以看出降雨量和I30與徑流量有較好的線性關(guān)系，復(fù)相關(guān)系數(shù)隨著流域治理程度的提高而逐漸降低，表現(xiàn)出流域產(chǎn)流的復(fù)雜性和多變性。

表3 不同治理效果下降雨因子與徑流量的回歸方程

表3方程中：P為降雨量，I30為30 min最大雨強(qiáng)，V為徑流量，V=aP+bI30+c為多元回歸方程。

3.2 徑流系數(shù)變化規(guī)律

1）典型降雨條件下徑流系數(shù)變化。徑流系數(shù)的大小與降雨量、植被、前期土壤濕度、土壤理化性質(zhì)、地形地貌等因素密切相關(guān)，在小流域地形地貌特征相似的前提下，對(duì)兩個(gè)小流域的比較見表4。從50～79.9 mm雨量段開始九仙山、平坡小流域徑流系數(shù)都有較大增長(zhǎng)，在這一雨量段九仙山流域增長(zhǎng)了10.375，而平坡小流域增長(zhǎng)了7.65，從不同雨量級(jí)看九仙山小流域的徑流系數(shù)明顯高于平坡小流域，說明治理程度的提高可以明顯減少?gòu)搅髁?，使得徑流系?shù)明顯降低。

表4 不同雨量級(jí)下的平均徑流系數(shù)

2）降雨量與徑流系數(shù)的回歸方程。

對(duì)降雨量和徑流系數(shù)進(jìn)行回歸分析，發(fā)現(xiàn)徑流系數(shù)Y和降雨量P呈很好的線性關(guān)系，見表5。相關(guān)系數(shù)都在0.9以上，而且九仙山流域的相關(guān)系數(shù)大于平坡流域，對(duì)不同治理效果下的流域徑流系數(shù)在0.01水平上進(jìn)行方差分析，看出兩個(gè)流域的徑流系數(shù)差異性呈顯著水平（F=8.379＞F0.01=7.28；P=0.021），說明小流域治理效果的好壞導(dǎo)致了流域徑流系數(shù)的差異。

表5 不同治理效果下降雨量與徑流系數(shù)的回歸方程

表5方程中：P為降雨量，Y為徑流系數(shù)，Y=aP+b為直線回歸方程。

4 結(jié)語(yǔ)

1）九仙山、平坡流域不同雨量級(jí)下平均徑流量分析表明治理效果差流域的徑流量明顯大于治理效果好的流域的徑流量，對(duì)降雨量和I30與徑流量進(jìn)行相關(guān)分析，其相關(guān)程度均達(dá)到極顯著水平，其相關(guān)程度的排序?yàn)榻涤炅浚綢30，降雨量和I30與徑流量的相關(guān)性排序?yàn)槠狡铝饔颍季畔缮搅饔?。?duì)降雨量、I30與徑流量進(jìn)行逐步回歸，發(fā)現(xiàn)降雨量、I30與徑流量呈多元線性關(guān)系，其中九仙山流域?yàn)閂=109.897P+135.725I30-3596.809，復(fù)相關(guān)系數(shù)R=0.948；平坡流域?yàn)閂=95.568P+70.469I30-2758.068，復(fù)相關(guān)系數(shù)R=0.934。

2）九仙山、平坡流域不同雨量級(jí)下平均徑流系數(shù)分析表明治理效果差的流域的徑流系數(shù)明顯大于治理效果好的流域的徑流系數(shù)，對(duì)降雨量和徑流系數(shù)進(jìn)行回歸分析，發(fā)現(xiàn)降雨量和徑流系數(shù)呈線性關(guān)系，其中九仙山流域?yàn)閅=0.265P-0.1084，R2=0.815；平坡流域?yàn)閅=0.2116P-2.4754，R2=0.7969。