開發建設項目新增水土流失的試驗研究方法

藺明華 藺 青 蘇佳園 田 博 陳小科

(黃河水利委員會晉陜蒙接壤地區水土保持監督局)

資源開發和基礎設施建設，是人類社會發展過程中必須進行的社會活動，新增水土流失也是客觀必然結果。近30 a來，隨著我國經濟社會發展的加速推進，開發建設項目大規模上馬，人為造成的水土流失日益嚴重。如何對開發建設項目已造成的水土流失進行調查與分析評價，對即將可能產生的水土流失進行科學預測，以便采取針對性措施，確保項目建設與生態環境建設雙贏，是急需研究解決的一個重大問題。為此，水利部將“黃河中游地區開發建設項目新增水土流失預測研究”課題列為水利技術開發基金項目，在開發建設項目比較集中的烏蘭木倫河流域進行了專項試驗研究。根據研究成果，本文僅對開發建設項目新增水土流失的試驗研究方法及其成果的應用做一簡要介紹。

1 開發建設項目新增水土流失特點

開發建設項目實施過程中，破壞了地表原有植被和土壤結構，產生了大量的棄土棄渣，為新增水土流失提供了動力基礎與物質基礎。同時，形成了較多的人工堆墊地面，這些新的下墊面，其土壤的機械物理組成與原生地面相比發生了很大變化，必然導致土壤抗蝕力的變化，進而引起土壤侵蝕的變化。一般而言，因新下墊面抗蝕力的下降，其侵蝕強度必然大于原生地面，產生新增水土流失。

開發建設項目在建設與生產過程中新增水土流失的主要表現形式:一是因廢棄的土、巖石或其混合物任意堆放未采取水土保持措施而所產生的水土流失，其侵蝕方式有重力侵蝕(瀉溜、崩塌、泥石流等)、水蝕和風蝕，以重力侵蝕和水蝕為主要方式;二是因破壞地表土壤結構和地面植被，開挖巖土使土壤基巖裸露而形成新的地表面(土質面或巖質面)，因其抗蝕力下降，比原地表侵蝕強度加大，其侵蝕方式有重力侵蝕、風力侵蝕和水力侵蝕等，一般以水蝕為主;三是根據設計使移動后的巖土按一定密實度在指定位置有序地堆放，但因堆積體表面未采取水保措施或水保措施尚未發揮效能而產生水土流失，其侵蝕方式有水蝕、風蝕等，一般以水蝕為主;四是因施行地下挖、采而導致地表下沉、地面裂縫、地下水的運動結構破壞和地表植被枯萎、死亡，所引起的穴陷、穴蝕和加重風蝕、水蝕等。

這四種新的水土流失表現形式，第一種我們稱之為棄土棄渣的水土流失，第二種稱之為裸露地貌的水土流失，第三種稱之為堆墊地貌的水土流失，第四種稱之為裂陷地貌的水土流失。根據調查，在開發建設過程中，第四種破壞地貌范圍不大，其新增水土流失量微小且侵蝕機理復雜。我們將裸露地貌和堆墊地貌統稱為“擾動地面”。所以，開發建設項目新增水土流失的主要來源地是棄土棄渣和擾動地面。

2 試驗研究內容與技術路線

黃河中游地區的土壤侵蝕多以水蝕為主，其影響因子主要有氣候、土壤、植被、地形和人類活動等5類。對開發建設形成的新下墊面而言，植被因子可不予考慮，所以開發建設形成的新下墊面的侵蝕影響因子主要為氣候、土壤、地形3類。開發建設新增水土流失量即為各類新下墊面因降雨侵蝕比原生地面增加的侵蝕量。本項研究的目的，就是通過試驗研究建立包括原生地面在內的各類下墊面以土壤侵蝕模數為目標量，以氣候因子、土壤因子、地形因子為預報因子的土壤侵蝕模型。

為達上述目的需開展試驗研究的主要內容:①典型區域開發建設棄土棄渣和人為擾動地面調查;②典型區域降雨資料收集和降雨特性分析;③典型區域水保、水利措施減沙量與河道沖淤量調查分析;④進行各類下墊面天然徑流小區和降雨量觀測;⑤在各類下墊面上開展人工降雨和放水沖刷試驗;⑥根據上述試驗結果建立不同下墊面的入滲方程，確定各種下墊面在不同頻率降雨作用下的水沙關系;⑦依據土壤侵蝕原理對典型區域所有開發建設項目的新增水土流失量進行分析計算，同時應用水文法和水保法分析同一典型區域所有開發建設項目的新增流失量，并對前述分析結果進行驗證;⑧開發建設項目區土壤抗蝕性研究;⑨建立用于開發建設項目新增水土流失量分析計算及預測的數學模型、土壤侵蝕系數及流失系數等。

3 試驗研究方法及其成果

3.1 天然降雨侵蝕試驗

在大柳塔選擇4種下墊面布設了11個2 m×5 m的天然徑流小區和雨量站，獲得2001、2002兩年間的5次有效降雨徑流資料，這些小區基本包括了開發建設過程中產生的下墊面類型。研究過程中，還參考借鑒了綏德水土保持科學試驗站在大柳塔獲得的1992－1994年間7個標準徑流小區的5次有效降雨徑流資料。這些資料將作為分析天然降雨條件下土壤入滲方程、降雨—徑流—泥沙關系和土壤抗沖性指標的主要依據。

3.1.1 天然降雨侵蝕試驗觀測內容與方法

(1)降雨量觀測。利用自記雨量計和標準雨量筒同時觀測降雨量及過程。

(2)前期土壤含水率。在每次降雨前分別取樣測定各小區的土壤含水量。

(3)土壤干容重。在每次降雨前分別取樣測定各小區的土壤干容重。

(4)洪水總量與產沙量。當一次產流降雨結束，及時測出徑流池的洪水總量，同時取渾水樣分析含沙量，進而分析計算次降雨情況下的產沙量和清水總量，以此作為建立土壤入滲方程及降雨—徑流—泥沙關系的依據。

3.1.2 徑流小區土壤入滲方程的推求

為了克服天然降雨徑流小區未觀測到徑流過程和產沙過程的不足，研究采用試算法與累計下滲曲線扣損法相結合來推求徑流小區的入滲率曲線，即根據試驗地的土壤類型，假定一條土壤入滲能力曲線。本項研究選取土壤入滲曲線形式為Horton型曲線:

式中:f——土壤入滲率，mm/min;

t——入滲時間，min;

a、b——系數;

k——指數。

根據每場降雨過程的實測雨強、土壤前期含水量以及假定的入滲曲線，利用累計下滲曲線扣損法，逐時段進行產流計算，得到總產流量。如果利用假定的入滲曲線所得到的總產流量與實測產流量之差滿足精度要求，則假定的入滲曲線就是要求的徑流小區土壤入滲曲線，否則，另假定一條土壤入滲曲線重新進行產流量計算，直至理論產流量與實測產流量之差滿足精度要求為止。

按照上述方法，編寫了C語言程序，依據實測資料，分別推求出11個徑流小區5場產流降雨條件下的Horton入滲方程(見表1)。

表1 不同下墊面小區天然降雨土壤入滲方程

3.2 人工降雨入滲試驗

3.2.1 試驗布設及觀測方法

(1)試驗布設。人工降雨入滲試驗的主要目的是建立各種下墊面下的土壤入滲方程，為推求不同頻率降雨條件下的降雨過程和徑流過程奠定基礎。為達此目的，在野外布設了5類下墊面11種坡度的34個2 m×1 m人工降雨入滲試驗小區，利用人工降雨設備按不同的降雨強度和降雨歷時實施了34場人工降雨試驗。降雨器由6個下噴式噴頭組成，噴頭距地面的距離不小于6 m。

(2)觀測內容與方法。在每一個小區上進行一場雨強衡定的人工降雨試驗，設計降雨強度和降雨歷時根據對試驗區域降雨特性分析結果分別確定為:1.0－2.5 mm/min和30－60 min之間。正式試驗前，先進行雨強率定試驗，使雨強接近設計值，且小區上的雨強均勻度在85%以上。試驗開始后，利用均勻布設于小區內的承雨器和接于小區出口集流槽上的徑流桶，每間隔2 min同時觀測一次降雨量、徑流量和徑流含沙量，直至實驗結束。

3.2.2 資料分析與試驗結果

根據觀測到的各時段徑流量及其含沙量，推算出該小區逐時段的清水徑流量和凈雨量;再根據各時段降雨量、凈雨量即可計算出逐時段的入滲量，并可繪制出該小區在一場人工降雨試驗過程中的土壤入滲過程線，進而擬合成Horton入滲方程。

在34場人工降雨試驗中，有部分試驗的計算結果因數據零亂不能應用;當年棄土棄渣小區試驗過程中，因其滲透性非常強，以致于降雨強度調到2.45 mm/min時仍然未產流，因而這類下墊面沒有試驗結果。對有效試驗結果經參數綜合匯總后得到各類下墊面的土壤入滲方程(見表2)。

表2 不同下墊面人工降雨土壤入滲擬合方程

3.3 降雨入滲試驗成果匯總

采用天然降雨和人工降雨兩種試驗方法得到的各類下墊面入滲方程，存在一定的差異。根據實驗條件和觀測方法的優劣勢，對兩種試驗獲得的成果進行綜合分析匯總，最終確定試驗區各類下墊面的入滲方程見表3。同時，將入滲方程f=a+bekt對時間積分，并以t=0時F=0為條件代入，得到降雨入滲能力方程:

式中:F——降雨入滲能力，mm。

(2)式和前面的(1)式即是降雨產流計算模型，可作為計算典型區域開發建設項目新增水土流失量的主要技術依據。

表3 不同下墊面土壤入滲方程匯總表

3.4 放水沖刷試驗

放水沖刷試驗的目的，是確定不同下墊面在降雨產流條件下的水沙關系和建立新增水土流失量預測方法體系所需的參數。由于條件所限，人工降雨小區試驗結果不能真實反映自然坡面的水沙關系，實際情況是小區上部降雨產生的徑流對小區也產生侵蝕，坡面越長，徑流量越大，徑流的侵蝕能力越大，坡面水沙關系越接近實際情況。為解決這個問題，需進行坡面模擬沖刷試驗，而坡面模擬沖刷試驗由于條件的限制，不能進行全坡長試驗，為真實反映各類下墊面在全坡長降雨沖刷條件下的水沙關系，小區放水模擬沖刷試驗流量設計時，考慮了各類下墊面的實際匯水面積問題。

3.4.1 試驗布設

模擬沖刷試驗小區規格為10.0 m×1.0 m，共布設13類下墊面小區73個。小區的周圍埋設1.0 m×0.3 m的鋼板，出口處用三角堰將小區的徑流收集到徑流桶中，以便進行測定和分析。根據各類下墊面小區降雨試驗得到的凈雨過程及小區上部匯水面積，確定的放水沖刷試驗流量分別為5、10、15、20、25(L/min)，模擬沖刷試驗時間一般在15－40 min之間。

3.4.2 試驗觀測項目

沖刷試驗前，需測定各下墊面土壤的容重、土壤顆粒級配、土壤前期含水量，并且采用體積法率定放水流量。

試驗開始后，每間隔1 min測量一次小區出口徑流桶承接到的渾水徑流的體積、重量，并將徑流桶靜置一段時間后，倒掉桶中上層清水，測定桶中泥沙層的重量，在泥沙層中部取部分泥沙樣，利用酒精烘干法測定泥沙中的含水量，據此反推桶中的泥沙量，多次測定取其平均值。

3.4.3 試驗結果

對每個小區放水沖刷過程得到的各時段徑流、泥沙數據進行分析，可得到各個時段的清水量(L)和泥沙量(kg)，再將每一類多個小區的試驗結果進行綜合，最后建立起該類小區放水沖刷條件下的水沙關系:

式中:y——產沙量，kg;

x——徑流量，L;

a、b——系數。

試驗結果見表4。

表4 不同下墊面水沙關系試驗成果表

4 試驗研究成果應用實例

應用上述試驗獲得的各類下墊面的降雨入滲方程和水沙關系方程，對開發建設初期的神東礦區所在區域—烏蘭木倫河流域的新增水土流失量進行了分析研究。

4.1 基本資料調查分析方法

烏蘭木倫河是黃河一級支流窟野河的上游，流域面積3 849 km2。在神東礦區大開發初期的1986－1998年的13 a間，共開工建設了包括工礦、煤礦、公路、鐵路、化工、建材、城鎮建設等項目208個。此外，還調查了農村民宅建設項目。這些項目總計擾動土地面積 5 392.93 hm2，產生棄土棄渣量13 681.53萬 t。

為了分析該時段開發建設項目的新增流失量，對該流域在1986－1998年間共發生的44次產洪降雨特性、每一場降雨中流域內21個雨量站的前期影響雨量進行了調查分析;對流域內所有開發建設項目的基本情況進行了調查，調查了各個項目棄土棄渣的數量、時間、堆棄位置、幾何形態、顆粒級配、棄土棄渣坡面上游來水面積，并調查了各個項目擾動地面的面積、數量、擾動時間、坡度、擾動形態及其上游來水面積;同時，將208個項目和農村民宅逐年建設面積，分配在該流域的21個雨量站控制的泰森多邊形內。此外，還在開發建設項目比較集中的晉陜豫“金三角”地區、陜西黃陵礦區及晉陜蒙接壤區，調查了棄土棄渣和各類擾動地面的植被及地表抗蝕性恢復時間;最后，對前述試驗獲得的降雨入滲方程和徑流泥沙關系方程，以開發建設產生各類下墊面的地面坡度和下墊面形成年份為參數進行了細化。

4.2 新增流失量分析方法

首先，根據超滲產流原理和初損扣損法，利用前面建立起來的降雨入滲方程，分析計算某個項目的每一種下墊面及其原生地面在逐次產洪條件下的產流量;其次，依據該種下墊面的水沙關系方程，分析計算其在逐次產洪降雨條件下的產沙量;該下墊面的產沙量與原生地面產沙量之差，即是該下墊面在本次產洪降雨條件下的新增流失量;最后，對該項目各類下墊面新增流失量求和，即獲得本次產洪降雨條件下該項目的新增流失量。至此，分析得到了某一個項目在一次產洪降雨條件下的新增流失量，依此方法，可分析獲得烏蘭木倫河流域所有開發建設項目在研究時段(1986－1998年)內的新增流失量。

“黃河中游地區開發建設項目新增水土流失預測研究”課題在研究過程中，還采用水文法分析計算了烏蘭木倫河流域開發建設項目在同一研究時段內的新增流失量，對上述試驗研究成果進行了驗證，并利用試驗研究資料建立起用于分析、預測開發建設項目新增水土流失量的數學模型等方法體系，這里不再詳述。

本文僅對開發建設項目新增水土流失試驗研究的內容、方法及利用其成果對典型區域新增水土流失量進行分析的方法做了簡要介紹，不足之處，敬請提出寶貴意見。