周耀華，張 濤，郭國先，王天巍

（1武漢市水土保持監(jiān)測站，武漢430010；2華中農(nóng)業(yè)大學資源與環(huán)境學院，武漢430070；3武漢市黃陂區(qū)蔡店水保站，武漢430323）

水土流失預測是生產(chǎn)建設項目水土保持方案中的重要內(nèi)容之一，是水土保持措施布局和設計的基礎，也是項目區(qū)水土保持動態(tài)監(jiān)測與分析評價的一項重要內(nèi)容［1］。目前，大部分水土流失預測主要依據(jù)RULSE模型來估算。野外對于水土流失背景值的測定多采用徑流小區(qū)法，但是由于地形條件的限制，常常找不到合適的位置布設徑流小區(qū)，因此，水土流失背景值的數(shù)據(jù)很難監(jiān)測到。所以，在一個特定的區(qū)域內(nèi)，水土流失監(jiān)測徑流場的設立變得至關重要。徑流場監(jiān)測數(shù)據(jù)的合理利用成為快速、準確預測某區(qū)域土壤流失量的有效途徑。同時該預測方法的研究將有利于探索區(qū)域水土流失規(guī)律，尋求水土流失防治途徑和方法［2］。

水土流失監(jiān)測設立的坡面徑流場是建立在自然坡面上的小集水區(qū)域。一般面積在1hm2左右，比自然徑流小區(qū)大［3］。大型坡面徑流場由于坡面寬闊，徑流易于集中，降雨產(chǎn)流后，坡面徑流形態(tài)除少部分為薄層片流外，大多為股流，侵蝕形態(tài)既有片蝕和細溝侵蝕，也有淺溝侵蝕、谷坡侵蝕及重力侵蝕，其測定的土壤流失量代表了自然坡面的實際土壤流失量［4］。這類小區(qū)反映了下墊面異質(zhì)性相互作用后產(chǎn)生的土壤侵蝕狀況。坡面徑流場更貼近自然狀態(tài)，監(jiān)測結果可以直接用于實際工作需要。鄭粉莉、唐克麗等也研究得出大型徑流場所監(jiān)測的泥沙數(shù)據(jù)和小流域的監(jiān)測數(shù)據(jù)基本一致，而徑流數(shù)據(jù)卻相差較大［4］。所以徑流場所監(jiān)測的數(shù)據(jù)可以作為特定區(qū)域生產(chǎn)建設項目水土保持方案編制中非耕地水土流失背景值預測的類比依據(jù)，預測該區(qū)域的土壤流失量。生產(chǎn)建設項目水土保持監(jiān)測中計算的水土流失背景值是土壤流失量，指土壤及其母質(zhì)在侵蝕營力作用下，產(chǎn)生位移并通過某一觀察斷面的泥沙數(shù)量。用土壤流失量推演侵蝕量時，一般只需要除以泥沙輸移比，就能求得流域的土壤侵蝕量［5］。

本文以武漢市黃陂區(qū)為研究區(qū)域，該區(qū)于1986年在北部低山區(qū)建立野佛溝徑流場，1987年開始開展降雨、徑流泥沙監(jiān)測試驗。通過對該徑流場多年來的監(jiān)測結果進行分析，模擬出預測土壤流失量的模型，以期為該區(qū)水土流失及生產(chǎn)建設項目背景值的預測提供依據(jù)。

1 數(shù)據(jù)與方法

1．1 研究區(qū)概況

武漢市黃陂區(qū)位于武漢市中心城區(qū)北部，地理坐標為東經(jīng)114°09′—114°37′，北緯30°41′—31°22′。東臨新洲區(qū)和紅安縣，北與大悟縣為界，西與孝感市和孝昌縣接壤。土地總面積2 257．20km2，總人口110萬，是武漢市面積最大和人口最多的一個區(qū)。全區(qū)境內(nèi)地勢北高南低，以丘陵平原為主，兼有低山、洼地等地貌類型。屬濕潤的亞熱帶季風氣候，年平均氣溫16．3℃，年降雨量1 100～1 200mm，是武漢市水土流失重點治理區(qū)。

野佛溝徑流場位于黃陂區(qū)蔡店鄉(xiāng)，以花崗巖、片麻巖發(fā)育的土壤為主，試驗區(qū)垂直投影長度277m，寬64．1m，封閉徑流場總面積1．775hm2，平均坡度25°。監(jiān)測開始時林草覆蓋率31．7%，目前林草覆蓋率上升為71．5%。該區(qū)建有巴歇爾槽監(jiān)測站，巴歇爾槽長5m，高1m，喉導寬0．5m，下設60°三角量水堰，堰口海拔高程147．4m。

1．2 數(shù)據(jù)應用范圍

野佛溝徑流場共進行了23a的徑流和泥沙監(jiān)測，根據(jù)《水利水電工程水文計算規(guī)范》SL278－2002，“徑流頻率計算依據(jù)的資料系列應在30a以上”，所以降雨徑流需要繼續(xù)監(jiān)測8a以上或插補延長，才能作為設計依據(jù)。懸移質(zhì)泥沙“統(tǒng)計系列長度定為20a”，推移質(zhì)泥沙統(tǒng)計系列長度“不宜少于10a”，野佛溝徑流場監(jiān)測站泥沙監(jiān)測資料超過20a，可將統(tǒng)計泥沙特征值作為黃陂北部低山區(qū)河流、水庫、塘堰泥沙淤積預測，土壤流失量可以作為生產(chǎn)建設項目水土保持方案編制中非耕地水土流失背景值預測的類比依據(jù)。

1．3 研究方法

黃陂區(qū)北部低山區(qū)非耕地水土流失背景值預測時運用了我國相關學者改進后的、適用于我國水土流失評價的RUSLE模型［6］，其表達式如下：

式中：A——土壤年流失量；f——使A代表我國單位量綱［t／（km2·a）］的綜合轉(zhuǎn)換常數(shù)；R——降雨侵蝕力因子；K——土壤可侵蝕性因子；LS——坡長坡度因子的乘積；CP——植被與經(jīng)營管理因子、水土保持措施因子的乘積；LS、CP為無量量綱單位。

在各參數(shù)計算過程中，降雨侵蝕力因子R、土壤可侵蝕性因子K的確定非常困難，所以在同一個水土流失類型區(qū)內(nèi)，近似認為R、K不變，僅考慮坡長、坡度、植被與經(jīng)營管理因子的變化。本文根據(jù)生產(chǎn)建設項目水土流失預測和通用土壤流失方程RUSLE的特點，利用類比法［7］推導土壤流失預測模型。

2 結果與分析

2．1 野佛溝徑流場監(jiān)測結果及分析

野佛溝徑流場主要監(jiān)測了降雨、徑流、泥沙和土壤流失量。其中1997年因進水喉道嚴重下沉，只監(jiān)測了降雨量，未監(jiān)測土壤流失量。監(jiān)測結果詳見表1。

表1 野佛溝徑流場水土流失監(jiān)測結果

根據(jù)以上數(shù)據(jù)通過計算得出徑流場內(nèi)多年平均徑流系數(shù)μ1＝0．402，平均徑流系數(shù)標準差σ1＝0．114，變異系數(shù)Cv1＝σ1／μ1＝0．283。多年平均土壤流失量μ2＝186t／（km2·a），土壤流失量標準差σ2＝391，變異系數(shù)Cv2＝σ2／μ2＝2．10。土壤流失量標準差、變異系數(shù)大，這可能是由于土壤流失量受影響因子的影響后波動較大，所以變化比較明顯。

隨著年降雨量的增大，徑流深度也在增大，而降雨量與土壤流失量的相關性不明顯（圖1），這是因為土壤流失量受環(huán)境因素的影響遠比徑流深度的影響復雜［8］，同時也說明土壤流失量與降雨強度密切相關。1987年土壤流失量明顯高于其他年份，經(jīng)調(diào)查，主要原因是1986年建設巴歇爾槽監(jiān)測站時，在徑流場內(nèi)采用魚鱗坑植樹擾動地表，導致1987年流失量明顯提高，該年監(jiān)測受人為活動的影響，理論上說應該剔除該年監(jiān)測結果，但是考慮到該類型區(qū)其他小流域都可能會受到不同程度的人為活動影響，為了保持應用該結果的一致性所以不剔除該年監(jiān)測結果。但是2001年、2006年、2009年的土壤流失量明顯偏低，為了保證與實際侵蝕情況的一致性，在計算過程中剔除了這三年監(jiān)測的土壤流失量。

圖1 降雨、徑流和土壤流失量變化規(guī)律

2．2 水土流失預測模型的建立

以徑流場多年平均土壤流失量為類比依據(jù)，采用比例估算法預測黃陂區(qū)北部低山區(qū)土壤流失背景值，即：

項目區(qū)土壤流失量＝186×LS×CP／L0S0×

2．2．1 模型參數(shù)的計算

（1）坡長坡度因子LS的確定。坡長坡度因子反映的是地形地貌特征對土壤侵蝕的影響［9］，LS的確定根據(jù)《水土保持監(jiān)測技術》推薦的適用于南方紅壤區(qū)的公式：

式中：h——相對高差；β——地面平均坡度角（°）。

野佛溝徑流場的相對高差：

野佛溝徑流場的坡長坡度因子：

（2）植被與經(jīng)營管理因子、水土保持措施因子的確定。植被與經(jīng)營管理因子、水土保持因子CP表征自然植被和人為水保措施通過截留、改良土壤結構、改變地表粗糙度等對土壤侵蝕產(chǎn)生的抑制作用［10］。C、P取值范圍為0～1。地面完全裸露時C值為1，當?shù)孛姹Ｗo良好時C值為0．001；未采取任何水保措施的地區(qū)P值為1，水土保持很好、無侵蝕風險的地區(qū)P值為0．001。C值的確定綜合參考植被覆蓋度（表2）與土地利用類型數(shù)據(jù)（表3），P值的確定綜合參考耕地類型數(shù)據(jù)和地形坡度數(shù)據(jù)（表4）。

表2 不同植被覆蓋條件下的C值

表3 不同土地利用條件下的C值

表4 不同坡度和水保條件下的P值

通過計算得出監(jiān)測期野佛溝徑流場平均林草覆蓋率為51．6%，所以其植被與經(jīng)營管理因子取有林地C0＝0．11，水土保持措施因子取等高壟作P0＝0．3。

將以野佛溝徑流場計算出的各個參數(shù)值代入，即得出黃陂區(qū)北部低山區(qū)水土流失背景值的預測模型：

項目區(qū)土壤流失量

2．2．2 預測模型應用范圍 不同區(qū)域地形地貌、氣象、植被類型、土壤、擾動地表組成物質(zhì)及形態(tài)、侵蝕類型等各個方面都會有所不同。所以在類比法選擇時，要充分考慮水土流失主要因子的相似性、可比性［11］。該預測模型只適用于武漢黃陂區(qū)北部低山區(qū)，同一個水土流失類型區(qū)內(nèi)，近似認為R、K不變，所以在該區(qū)域內(nèi)依據(jù)上述方法，合理確定預測小流域的坡長坡度因子LS，植被與經(jīng)營管理因子、水土保持因子CP就能估算出項目區(qū)水土流失背景值。

2．3 水土流失預測模型的驗證

根據(jù)建立的預測模型，在武漢黃陂區(qū)北部低山區(qū)選取不同土地利用狀況下的片區(qū)進行對比，對比情況如表5所示。

表5 不同土地利用方式下土壤流失量預測值與實測值比較 t／（km2·a）

表5表明，在不同的土地利用條件下，應用建立的土壤流失預測模型得出的土壤流失預測值與實測值十分接近，由此可見由小區(qū)觀測數(shù)據(jù)研究得出的土壤流失預測模型在武漢黃陂區(qū)北部低山區(qū)是適用的。

3 結 論

在水土流失重點治理區(qū)建立的大型坡面徑流場不僅可以研究流域內(nèi)不同地形部位的泥沙來源和產(chǎn)沙部位，而且可以研究整個坡面土壤侵蝕過程及土壤侵蝕與泥沙搬運的關系，建立坡面土壤流失預報方程［4］。野佛溝徑流場23a實測土壤流失量特征值即可用于黃陂北部低山區(qū)河流、水庫、塘堰泥沙淤積預測，同時建立的武漢市黃陂北部低山區(qū)水土流失背景值預測模型可以預測該區(qū)非耕地水土流失背景值，需要強調(diào)的是，預測得出的土壤流失量是土壤侵蝕量的一部分。該區(qū)降雨、土壤抗蝕性基本相同，僅需考慮坡長坡度因子、植被與經(jīng)營管理因子、水土保持因子的變化對土壤流失量的影響，建立的模型計算方法科學，是從土壤流失監(jiān)測到土壤流失預測研究的重大突破，為區(qū)域水土流失背景值的預測提供了一種新的思路。

［1］ 聶國輝，許建新，葉永棋，等．Cs－137在浙江省水土流失背景值監(jiān)測中的應用［J］．浙江水利科技，2008（3）：22－26．

［2］ 張志萍，陳智漢，王志雄．黃土高塬溝壑區(qū)典型流域水土流失監(jiān)測發(fā)展歷程及對策［J］．中國水土保持科學，2011，9（4）：41－45．

［3］ 趙愛軍，馬力剛．湖北省水土保持監(jiān)測點建設實踐［J］．水土保持通報，2009，29（2）：117－182．

［4］ 鄭粉莉，唐克麗，白紅英．標準小區(qū)和大型坡面徑流場徑流泥沙監(jiān)測方法分析［J］．人民黃河，1994，17（7）：19－21．

［5］ 景可，焦菊英，李林育，等．中國南方紅壤丘陵區(qū)土壤侵蝕量與泥沙輸移比研究［J］．水利學報，2010，41（12）：1476－1482．

［6］ 陳云明，劉國彬，鄭粉莉，等．RUSLE侵蝕模型的應用及進展［J］．水土保持研究，2004，11（4）：80－83．

［7］ 曲德雙．類比法預測開發(fā)建設項目水土流失幾個問題的探討［J］．黑龍江水利科技，2008，36（5）：26－27．

［8］ 江淼華．閩北裸露坡地水土流失與降雨量的關系研究［J］．梧州學院學報，2010，20（1）：7－11．

［9］ 楊艷生．區(qū)域性土壤流失預測方程的初步研究［J］．土壤學報，1990，27（1）：73－79．

［10］ 王秀艷，郭兵，姜琳．基于 USLE、GIS、RS的流域土壤侵蝕研究進展［J］．亞熱帶水土保持，2012，24（1）：42－48．

［11］ 潘獻鴻，李海林，高榮．類比法預測建設項目水土流失量［J］．水土保持應用技術，2009（6）：14－16．