推薦一種泥石流坡面侵蝕物源量的計算方法

溫清茂，鐘東

推薦一種泥石流坡面侵蝕物源量的計算方法

溫清茂，鐘東

（四川省地礦局九一五水文地質工程地質隊，四川 眉山 620020）

在收集、研究前人研究成果的基礎上，應用土壤水蝕模型和流域泥沙輸移比初步建立了泥石流坡面侵蝕物源量計算模型，同時列舉了降雨侵蝕力因子、土壤可蝕性因子、地形因子及植被覆蓋因子等水蝕因子的標準算法和經驗算法，具有較強的理論依據，對泥石流坡面侵蝕物源量計算具有借鑒意義。

泥石流；坡面侵蝕物源；土壤水力侵蝕；計算方法

坡面侵蝕物源是泥石流形成過程中物源條件的重要組成部分，在泥石流形成過程中占有重要地位，特別是對以坡面侵蝕為主的泥石流的形成起著至關重要的作用。泥石流坡面侵蝕物源量系指參與泥石流活動的坡面侵蝕物源動儲量，傳統估算方法主要有工程經驗法和侵蝕強度法，工程經驗法即按所占靜儲量的比重估算，給定參與泥石流活動的坡面松散物質量的百分比[1]，或依據其它類似工程進行類比，按平均侵蝕深度估算；侵蝕強度法即依據工程地域、植被覆蓋程度和土地類型等定性判斷其侵蝕強度等級，再按侵蝕強度分級標準的區間值給定其平均侵蝕模數或平均侵蝕厚度[2]，這也是《泥石流災害防治工程勘查規范》（T /CAGHP 006-2018）推薦的算法[3]。上述兩種方法雖然計算簡單，但受人為因素影響很大，估算結果往往與實際情況懸殊，計算精度很差。事實上，泥石流坡面物源的啟動多依靠雨滴濺蝕和坡面流沖刷而完成，其與當地降雨強度、土層性質、地形條件、植被發育程度和人類活動等因素密切相關，應屬于土壤水力侵蝕理論的范疇。為此，在綜合考慮降雨特征、土壤特性、地形條件、植被覆蓋程度、水保措施和耕作方式的基礎上，應用土壤水力侵蝕模型，推薦一種泥石流坡面侵蝕物源量的計算方法，旨在拋磚引玉。

1 坡面侵蝕物源量計算模型

坡面侵蝕系指坡面土壤在水動力（雨水及其坡面流）作用下發生的土壤水蝕現象，主要有面蝕（雨滴和地表徑流作用的結果，有濺蝕、片蝕和細溝蝕）和溝蝕（坡面徑流作用的結果，往往形成淺溝和沖溝），其侵蝕模數可采用土壤水蝕模型表征。《泥石流災害防治工程勘查規范》（T /CAGHP 006-2018）要求：坡面侵蝕物源量按侵蝕分區計算工程有效期內的侵蝕總量，但因其粒度較小，易被常年洪水帶走，僅部分可計為泥石流動儲量，故坡面侵蝕物源量應按常年洪水可輸移的粒徑所占比例予以折減[3]。據此，泥石流坡面侵蝕物源量的計算模型可由式（1）、式（2）[4，5]表征：

式中：為坡面侵蝕物源量（m3）；V為工程設計年期（a）；為泥沙輸移比；為泥沙干密度（t/m3）；MF為第計算分區的年平均水蝕模數（t/hm2·a）、面積（hm2）；為年平均土壤水蝕模數（t/hm2·a）；為降雨侵蝕力因子（MJ·mm/hm2·h·a）；為土壤可蝕性因子（t·h/MJ·mm）；為坡長坡度因子；為生物措施因子；為工程措施因子；為耕作措施因子。

2 土壤水蝕因子計算方法

2.1 降雨侵蝕力因子R

1）標準計算方法（規范推薦算法）：降雨侵蝕力是降雨引起土壤侵蝕的潛在能力，是降雨特性的函數，反映了雨滴對土壤顆粒的擊濺分離和降雨形成徑流對土壤沖刷的綜合作用。標準的計算方法為一次降雨總動能E與該次降雨的最大30?min雨強I30的乘積，采用剔除日降雨量小于12mm的非侵蝕性降雨后的多年降雨過程資料進行計算。計算公式[2、3]為：

2）等值線圖內插法：第一次全國水利普查之水土保持情況普查根據上述標準計算方法，分別編繪了各省（市、區）及全國多年平均降雨侵蝕力等值線圖[6，7]，可根據等值線圖內插取值。

表1 岷江上游和西藏高原主要土壤可蝕性K值 單位：t·h/MJ·mm

3）簡易計算方法：標準計算方法采用EI30指標度量，需要長序列降雨過程資料且計算繁瑣、費時，為此，很多學者提出了采用常規降雨資料（年降雨量、月降雨量）計算多年平均降雨侵蝕力的簡易算法。伍育鵬等[8]通過對國內外具有代表性的9種簡易算法的計算結果與降雨侵蝕力指標E10進行比較，較適用于我國的降雨侵蝕力簡易算法有：

a）馬志尊（海河流域太行山）：適用于除西北地區以外的廣大地區，相對偏差4.7%~93.7%。

式中：為年降雨侵蝕力（100?ft·sht·in/acre·h），乘以17.02轉化為MJ·mm/hm2·h·a；P為多年月平均降雨量（mm）；為多年平均降雨量（mm）。

b）YU等（澳大利亞南部地區）：適用于我國大部分地區，相對誤差0.3%~193%，平均44%。

式中：為年降雨侵蝕力（MJ·mm/hm2·h·a）。

2.2 土壤可蝕性因子K

1）標準計算方法：土壤可蝕性因子表征土壤被沖蝕的難易程度，反映土壤遭受外營力剝蝕和搬運的敏感性，是影響土壤水蝕的內在因素。一般采用EPIC模型中Williams等提出的計算方法[5，9]：

式中：為土壤可蝕性因子（ton·h/100ft·tonf·in），乘以0.1318轉換為t·h/MJ·mm；a為砂粒含量（%）；i為粉粒含量（%）；l為黏粒含量（%）；為有機碳含量/%。

2）查圖法：第一次全國水利普查之水土保持情況普查根據標準計算方法，利用第二次全國土壤普查資料分別編繪了各省（市、區）土壤可蝕性空間分布圖，可查詢對應土壤的可蝕性K值。

3）類比法: 多年來，國內眾多學者利用不同地區的小區資料對土壤可蝕性進行了分析計算，積累了大量的數據，可根據土壤類型類比取值。如岷江上游[9]和西藏高原[10]地區的K值空間分布如表1。

表2 全國分省份土壤可蝕性K值 單位：t·h/(MJ·mm)

4）平均值法：2013年，梁音等[11]在收集、整理了全國31個省（市、區）17558個土壤樣品理化數據的基礎上，采用EPIC模型中Williams提出的K值公式計算，形成全國水蝕區土壤可蝕性K值分布圖，研究表明：全國水蝕區K值變化范圍為0.0004~0.0828t·h/MJ·mm，平均K 值為0.033t·h/MJ·mm（表2），呈現出由北向南逐漸減小的趨勢，黃土區K 值最大，紫色土和紅壤區北部的K 值居中，最小的分布在青藏高原北側。

2.3 坡長坡度因子LS

1)標準計算方法:坡度坡長因子即地形因子。隨坡度的增大，土壤顆粒的抗滑力減小，水流動能增大，水土流失量隨之增大，坡面越長，流失量也越大。地形因子LS表示在其它條件均相同的情況下，某一給定坡長（22.13?m）和坡度（9%）的坡面上，土壤流失量與標準徑流小區典型坡面土壤流失量的比值，它是土壤水蝕動力的加速因子。其計算公式[5]為：

表3 基于地貌及環境因子的泥沙輸移比計算公式

式中：為坡長因子；為坡度因子；λ為坡長（m）；為坡度（o）。

2)經驗公式估算方法: 我國部分地區坡度坡長因子LS經驗公式主要有[12]：

a）黃土高原（汪忠善）：

b）東北黑土地區（張憲奎）：

c）閩東南地區（黃炎和）：

4）南方紅壤地區（楊艷生）：

式中：為坡高（m）。

2.4 生物措施因子B

生物措施因子即植被覆蓋因子。植被具有截持降雨、分散雨量、降低雨滴沖擊力、減緩徑流、降低流速、提高滲透速度、保土固土等生態功能，對控制土壤流失起著決定性的作用，受植被、作物種植順序、栽培措施、降雨分布等眾多因素制約，一般采取如下方法求算。

1)標準計算方法: 利用植被覆蓋度結合土地利用和降雨侵蝕力季節分布，根據各半月時段的土壤流失比例和各半月時段降雨侵蝕力占全年的比率，計算生物措施因子[5]：

2）經驗公式估算法：蔡崇法等[13]利用徑流小區人工降雨和部分天然降雨的觀測結果，通過計算坡面產沙量與植被覆蓋度的相關關系，建立了植被覆蓋度與生物措施因子間的數學關系：

式中：為生物措施因子；為植被平均覆蓋度（%）。

該方法在國內已被較多學者所采用[14]。從式（22）可以看出，植被覆蓋度≥78.3%時其B因子為0，因此植被覆蓋度在80%以上的地區可不考慮土壤水蝕（坡面侵蝕）作用。

2.5 工程措施因子E

對于泥石流坡面侵蝕物源區，一般為自然斜坡、農用地，取E值=1[4]。

2.6 耕作措施因子T

耕作措施因子即水保措施因子，指在特定水保措施的土壤流失量與起伏地耕作的相應土壤流失之比。自然斜坡和坡耕地（順坡耕作）的T值取1.0，橫坡耕作（包括等高耕作、等高帶狀種植、修梯田等）的T值取0.5[4]。

3 流域泥沙輸移比估算方法

流域泥沙輸移比表征流域輸沙總量與流域產沙量的比值[15]，二者之差即為該流域范圍內的泥沙沉積量。影響泥沙輸移比的因素很多，一般可概括為四類：一是地貌及環境因子，如流域面積、形態、坡度坡長、溝道特征等；二是降水及水文因素，如降雨強度、降雨歷時等；三是侵蝕土壤的粒徑與土壤質地結構；四是植物覆被和土地利用狀況等[16]。為方便計算，本文僅列舉考慮地貌及環境因子的泥沙輸移比估算方法（表3）。

此外，也可根據相關學者已有研究總結的部分流域泥沙輸移比分布進行類比確定。謝旺成[16]根據國內外泥沙輸移比（單因子和多因子的影響）計算的實例，總結了主要流域泥沙輸移比數值的變化范圍，可參考利用。

4 應用實例

4.1 研究區概況

以四川省甘孜州鄉城縣定波鄉萬絨溝泥石流為例。其流域形態呈基本對稱的扇形，流域面積22.85?km2，主溝長7.04?km，溝口高程3 009?m，相對高差1 323?m，平均縱坡降172?‰。溝谷形態總體呈“V”型，岸坡一般35?°~45?°，多懸崖峭壁，間布較平緩的臺地，山脊渾圓。流域內森林植被發育密集，以山地針葉林和高山灌叢草甸帶為主，植被覆蓋率一般在80?%以上，水土流失較輕微，但近河谷地帶植被破壞較嚴重，植被覆蓋率一般小于40?%，土壤類型為灌叢草甸土，土壤層厚一般1~1.5m，以面蝕和溝蝕為主，侵蝕強烈。據此，圈劃坡面侵蝕物源點4處，分別編號為P1、P2、P3、P4。

4.2 影響因子取值

4.2.1 水蝕因子取值

1）降雨侵蝕力因子R：萬絨溝屬大陸性季風高原型氣候，降雨量少但集中，干濕季節分明，多年平均降雨量447.6?mm；選用式（8）計算，R=1 749?MJ·mm/hm2·h·a。

2）土壤可蝕性因子K：依據表1選取亞高山灌叢草甸土K=0.043?t·h/MJ·mm。

3）坡度坡長因子LS：選用式（17）計算。

4）生物措施因子B：選用式（20）計算。

5）工程措施因子E和耕作措施因子T：為自然斜坡，取E=1，T=1。

4.2.2 其它因子選取

工程設計年限按20?a，泥沙干密度[4]取1.35t/m3。流域泥沙輸移比選用式（25）計算，SDR=0.287。

4.3 計算結果

表4 坡面侵蝕物源量計算結果

根據式（1）、式（2）計算，其結果列于表4。可見，萬絨溝坡面侵蝕物源的侵蝕模數為5491~9577?t/km2·a，依據《土壤侵蝕分類分級標準》（SL190-2007）[4]劃分其侵蝕強度分屬強烈——極強烈等級，計算結果與定性判別的侵蝕強度基本吻合；在工程設計年限內，坡面侵蝕物源量48 847?m3。

5 結語

1）本文在收集研究國家相關規范和前人研究成果的基礎上，應用土壤水力侵蝕理論及流域泥沙輸移比初步建立了泥石流坡面侵蝕物源量的計算模型，各相關因子的計算和選取在采用我國通用計算方法的基礎上還推選了國內外相關學者較簡便的算法和其適用性，具有計算方便、精度較高的特點，有可能成為評價泥石流坡面侵蝕物源量的統一算法。

2）土壤水蝕因子的降雨侵蝕力R因子、土壤可蝕性K因子的標準算法源于第一次全國水利普查之水土保持情況普查培訓教材，已廣泛應用于土壤侵蝕調查，應充分利用其調查成果，采取等值線圖或因子圖內插，其計算精度會更高。

3）目前，我國的泥沙輸移比SDR研究區域主要集中在黃河流域，長江流域、珠江流域、淮河流域、山地丘陵地區等地區或流域的研究相對較少。由于各流域侵蝕條件和輸移條件不同，計算結果也不盡相同，所獲得的計算公式往往也受地域局限，難以推廣。在應用時還應考慮多因素對流域泥沙輸移比的影響。

[1] 宋志, 陳澤碩. 四川石棉簡槽溝泥石流發育特征與減災防災初探[J]. 地質災害與環境保護, 2015, 26(4): 12-17.

[2] 丁保艷, 李征. 鹿鳴鄉松樹林泥石流特征與防治[J]. 建材世界, 2017, 38(3): 103-105.

[3] 中國地質災害防治工程行業協會. T/CAGHP 006-2018 泥石流災害防治工程勘查規范[S]. 武漢: 中國地質大學出版社, 2018.

[4] 中華人民共和國水利部. SL 190-2007 土壤侵蝕分類分級標準[S]. 北京: 中國水利水電出版社, 2008.

[5] 國務院第一次全國水利普查領導小組辦公室. 第一次全國水利普查培訓教材之六: 水土保持情況普查[M]. 北京: 中國水利水電出版社, 2010.

[6] 周琳, 許武成, 尹義星. 近61年四川省降雨侵蝕力的時空變化趨勢分析[J]. 水土保持研究, 2018, 25(4): 1005-3409.

[7] 劉斌濤, 陶和平, 宋春風, 等. 1960—2009年中國降雨侵蝕力的時空變化趨勢[J]. 地理研究, 2013, 32(2): 245-256.

[8] 伍育鵬, 謝云, 章文波. 國內外降雨侵蝕力簡易計算方法的比較[J]. 水土保持學報, 2001, 15(3): 31-34

[9] 楊孟, 李秀珍, 胡遠滿,等. 利用SEDD模型模擬岷江上游小流域的年產沙量[J]. 應用生態學報, 2007, 18(8): 1758-1764.

[10] 王小丹, 鐘祥浩, 王建平. 西藏高原土壤可蝕性及其空間分布規律初步研究[J]. 干旱區地理, 2004, 27(3): 343-346.

[11] 梁音, 劉憲春, 曹龍熹, 等. 中國水蝕區土壤可蝕性K值計算與宏觀分布[J]. 中國水土保持, 2013, 10: 35-40.

[12] 王萬忠, 焦菊英. 中國的土壤侵蝕因子定量評價研究[J]. 水土保持通報, 1996, 16(5): 1-20.

[13] 蔡崇法, 丁樹文, 等. 應用USLE模型與地理信息系統IDRISI預測小流域土壤侵蝕量的研究[J]. 水土保持學報, 2000, 14(2): 19-24.

[14] 吳昌廣, 李生, 任華東, 等. USLE/RUSLE模型中植被覆蓋管理因子的遙感定量估算研究進展[J]. 應用生態學報, 2012, 23(6): 1728-1732.

[15] 中華人民共和國國家質量監督檢驗檢疫總局, 中國國家標準化管理委員會. GB/T 20465-2006 水土保持術語[S]. 北京: 中國標準出版社. 2006.

[16] 謝旺成, 李天宏. 流域泥沙輸移比研究進展[J]. 北京大學學報(自然科學版), 2012, 48(4): 685-694.

[17] 吳傲, 李天宏, 韓鵬. 黃河流域泥沙輸移比與流域面積的關系[J]. 泥沙研究, 2014, 1: 61-67.

[18] 張曉明, 曹文洪, 武思宏, 等. 泥沙輸移比尺度依存及分形特征[J]. 水利學報, 2013, 44(10): 1225-1232.

[19] 高燕, 張延玲, 焦劍, 等. 松花江流域不同空間尺度典型流域泥沙輸移比及其影響因素[J]. 中國水土保持科學, 2016, 14(1): 21-27.

A Method for Calculating the Amount of Slope Erosion Material Source of Debris Flow Recommended

WEN Qing-mao ZHONG Dong

(No. 915 Hydrogeological and Engineering Geological Team, BGEEMRSP, Meishan, Sichuan 620010)

A model for calculating the amount of erosion material on the slope of debris flow is established by means of soil water erosion model and sediment delivery ratio on the basis of previous achievements. At the same time, standard algorithm and empirical algorithm for sediment delivery ratio and water erosion factors such as rainfall erosive power factor, soil erodibility factor, terrain factor and vegetation cover factor are enumerated.

debris flow;slope erosion material source; soil water erosion;sediment delivery ratio

2018-07-02

溫清茂（1965-），男，四川蓬溪人，高級工程師，主要從事水工環、地質災害方面的工作。

P642. 23

A

1006-0995（2019）02-0289-05

10.3969/j.issn.1006-0995.2019.02.023