基于GIS的多年土壤侵蝕對安基山水庫淤積的影響

鐘 磊，周 迎，劉 暢，程公德(. 濟南市水利建筑勘測設計研究院，濟南 5004；．河海大學水利水電學院，南京 0098)

0 引 言

土壤侵蝕是土壤及其母質在外營力作用下被破壞、剝蝕、搬運和沉積的過程，其輸出程度受人為因素和自然因素的綜合影響[1]。土壤侵蝕引起的水土流失不僅使土壤養分顆粒流失和土地退化，還使流域上游大量泥沙向下遷移造成水庫和河道淤積，以致影響水庫的正常使用[2,3]。因此，對土壤侵蝕的定量檢測和水庫淤積的估算是有效治理水土流失和保護水庫的前提。

GIS(Geographic Information System)具有的數據提取、管理和空間分析功能，在多尺度、大范圍集水地以及各級管理區域水土保持預測模型中應用極其廣泛[4]，如游松財等[5]在GIS技術支持下，運用通用土壤流失方程USLE(Universal Soil Loss Equation)對灌溪鄉土壤侵蝕進行定量估計，計算結果與實測數據有良好的相關性；周偉峰等[1]利用RS(Remote Sensing)和GIS技術在通用土壤流失方程的框架上建立區域土壤侵蝕模型，對北京市密云水庫土壤侵蝕量進行定量估算。針對水庫淤積計算問題，很多學者進行了深入研究。如焦恩東等[6]對萬兩河水庫進行水庫泥沙淤積分析計算，包括淤積程度判斷、淤積形態判別、淤積量計算和壩前淤積高程計算，并提出防治水庫淤積措施；王彥芳等[7]基于泥沙隨機沉降的原理，結合汾河靜樂水文站的歷年實測資料，建立汾河水庫泥沙淤積的隨機數學模型，并應用模型對汾河水庫庫區泥沙淤積發展情況進行模擬預測。然而以上研究并未針對流域上游土壤侵蝕和下游水庫淤積情況進行探討，不能反映多年土壤侵蝕變化對下游水庫淤積影響的嚴重程度。因此，在防治水庫泥沙淤積工作中，研究多年土壤侵蝕變化對下游水庫淤積的影響十分迫切。

本文研究的是在GIS技術支持下多年土壤侵蝕對水庫泥沙淤積的影響，以安基山水庫為例，分析2007-2011年5 a間流域上游土壤侵蝕對安基山水庫淤積情況。利用RS遙感影像數據獲取流域內的植被覆蓋度和土地利用分類信息，收集降雨資料和土壤數據；在流域DEM數據基礎上利用ArcGIS表面分析功能提取流域內的坡度和坡長等地形數據；采用通用土壤流失方程USLE建立流域土壤侵蝕定量檢測模型，利用ArcGIS空間分析技術對安基山水庫集水區域連續5 a間的土壤侵蝕量進行定量估算；根據多年土壤侵蝕量分析安基山水庫攔沙和淤沙量，以期為水土流失和水庫泥沙淤積治理提供指導作用。

1 研究區概況及數據準備

1.1 研究區概況

安基山水庫位于七鄉河上游，是南京市江寧區的村鎮飲用水源地(119°3′～119°7′E、32°4′～32°8′N)，流域面積為15.74 km2，庫區與句容市亭子鎮接壤。水庫壩頂長350 m，壩頂高程53 m，壩頂寬6 m。水庫總庫容618 萬m3，屬于小(一)型水庫。流域內有安基山村、孟塘、葉庵、鹿山、袁家咀、韓家邊、螺絲沖、、華山村等村。流域上游為丘陵地形，水庫附近地勢較平緩，流域海拔范圍23～403 m。流域范圍內以耕地和林地為主，其中林地占大部分。多年平均降雨量1 047 mm，降雨多集中在6-9月，暴雨一般多集中在7-8月。

1.2 數據準備

研究中用到的數據包括該地區的數字高程模型數據(DEM)、Landsat7遙感影像數據(2007-2011年)、降雨資料(2007-2011年)、土壤數據。其中DEM數據用于提取地形信息[見圖1(a)]，根據匯水區域分水嶺利用ArcGIS軟件中Arc Hydro Tools功能確定安基山水庫集水流域范圍[8][見圖1(b)]。

遙感影像成像時間選擇在6-9月為最佳，并用流域范圍對遙感數據進行裁剪，進而應用監督分類的方法處理遙感影像。如2011年6月27日研究區Landsat遙感影像數據經圖像預處理、影像增強、監督分類、人工校核等，獲取該時期集水流域內農田、水域、村莊、裸地、林地等土地利用類型[見圖1(c)]。土壤數據是全國第二次土壤普查土壤類型及其理化性質分析資料。

圖1 研究區數據

2 研究方法

2.1 模型建立

通用土壤流失方程(USLE)是眾多土壤侵蝕預報模型應用最為廣泛的模型[9,10]，它具有形式簡單、參數易求、利于計算、適用于資料短缺的地區等特點。USLE模型綜合考慮人為因素和自然因素的影響，其核心內容是R、K、L、S、C和P6個因子的計算，具體表達式如下：

A=fRKLSCP

(1)

式中：A為年土壤侵蝕量，t/km2；L為坡長因子；S為坡度因子；R為降雨侵蝕因子；K為土壤可蝕性因子；C為植被覆蓋度(管理)因子；P為水土保持因子；f為單位換算系數，若A取用英美單位，f=1，若取我國常用的t/(km2·a)，f=224.2。

2.2 因子計算

2.2.1R因子計算

降雨侵蝕因子是雨滴擊濺產生削蝕現象和雨水徑流形成地表沖涮作用的綜合評價指標，其計算是采用吳素業[11]研究提出的適合我國長江流域地區降雨侵蝕因子的簡單計算方法，以此根據安基山流域多年降雨資料計算降雨侵蝕因子。由于流域面積小，降雨分布比較均勻，全流域每年只采用一個R值，并利用ArcGIS軟件轉化為柵格圖像。R由下式計算：

(2)

式中：Pi為月降雨量，mm；170.2為轉化為R的美制單位換算系數。

2.2.2K因子計算

土壤可蝕性因子K采用EPIC模型中的計算方法[12]，在ArcGIS軟件中由土壤類型分布進行賦值，再轉化為土壤K值得到柵格數據。計算如下：

(3)

式中：SAN為砂粒(0.05～2 mm)含量百分數；SIL為粉粒(0.002～0.05 mm)含量百分數；CLA為黏粒(<0.002 mm)含量百分數；C為有機碳含量百分數；SN=1-SAN/100。

2.2.3S、L因子計算

坡度和坡長因子是反映地形特征對土壤侵蝕的影響程度，其計算是利用ArcGIS軟件對研究區DEM數據處理得到。坡度因子和坡長因子計算公式[13]如下：

(6)

式中：θ為坡度，(°)；λ為坡長，m；m為坡度坡長系數；β為細溝和細溝間侵蝕的比率。

2.2.4C、P因子計算

植被覆蓋度因子C表示植物葉面生長覆蓋地面情況，可抑制土壤侵蝕發生，根據馬超飛等[14]建立的植被覆蓋度與C值數學關系式求得：

(7)

式中：C為植被覆蓋度因子；f為植被覆蓋度，可通過遙感數據歸一化植被指數求得。

水土保持因子P是指采用特定水土保持措施后的土壤侵蝕量與相應無實施順坡耕地土壤侵蝕量比值，P大小范圍為0～1，0表示未發生土壤侵蝕地區，1表示未采取任何保持措施地區。根據相關研究數據[12,15]知，自然植被區、裸地、農村居民區和建筑用地均為1，水域為0，水田為0.15，旱田為0.35。由該集水流域內土地利用類型，利用ArcGIS對C、P因子賦值計算，得到柵格化圖層。

2.3 多年侵蝕量統計

基于6大因子計算得到的柵格因子專題圖層，運用GIS技術空間分析功能進行疊加運算，得到2007-2011年安基山水庫流域土壤侵蝕模數G，并依據《土壤侵蝕分類分級標準》(SL190-2007)[16]對計算結果進行分級統計，分級級別為：微度[0～500 t/(km2·a)]；輕度[500～2 500 t/(km2·a)]；中度[2 500～5 000 t/(km2·a)]；強度[5 000～8 000 t/(km2·a)]。安基山水庫流域各點土壤侵蝕模數不同，故根據土壤侵蝕模數柵格圖層屬性值和像元大小，計算安基山水庫流域每年向下游水庫輸沙量W(入庫泥沙量)。

3 水庫淤積

土壤顆粒在降雨作用下隨著地表徑流形成的匯流進行遷移，當進入水庫后，水深增加，水流速度減小，水流挾沙能力降低，導致水流中的部分懸移質和推移質泥沙在庫區內沉淀，使水庫不斷淤積，進而影響水庫使用壽命和污染水庫[6]。

3.1 水庫攔沙率

水庫攔沙率β表示土壤顆粒淤積在水庫庫底的沙量與進入水庫的沙量比值[3]，與水庫庫容和徑流量的比值有關，計算公式為：

(9)

式中：β為水庫攔沙率，%；T為年入庫水量，m3；f為流域年降雨量，mm；αi為各種下墊面徑流系數，由《室外排水設計規范》(GB50014-2006)[17]查得；Ai為各種下墊面面積，m2。

3.2 水庫淤積量

對長期蓄水情況下的中小型水庫來說，通過計算水庫淤沙量分析滯留在水庫庫底的泥沙情況，以便確定水庫的壽命或使用年限，可以根據水庫攔沙率法計算[3]：

S=Wβ

(10)

式中：S為水庫淤積泥沙量，以質量計算，t；β為水庫攔沙率；W為入庫泥沙量，t。

4 結果與分析

4.1 侵蝕強度分析

根據安基山水庫上游土壤侵蝕監測結果，統計土壤侵蝕強度分布狀況，2007-2011年土壤侵蝕量分級結果見圖2，5 a土壤侵蝕統計見表1。

4.1.1土壤侵蝕區域分析

由圖2可知，侵蝕級別所占面積關系是微度>輕度>中度>強度>，強度侵蝕主要發生在流域中下部，中度侵蝕分布區域緊鄰強度侵蝕發生區附近，輕度侵蝕主要是在流域中部和下部，其余均是微度侵蝕，無明顯侵蝕現象。

分析知，強度侵蝕主要發生在裸地和地勢起伏大的區域，其原因為裸地地面植被覆蓋度較小和未采取水土保持措施；地勢起伏大的區域地面坡度坡長較大，地表徑流速度和水流長度增加，致使坡度因子和坡長因子較大。中度和輕度侵蝕主要發生在人類活動頻繁地區和丘陵山區，其原因是人類活動頻繁地區土地利用類型發生了改變，丘陵山區山高坡陡等地形條件。微度侵蝕所占比例均在80%以上，該區植被覆蓋較好，大部分區域是以林地為主，人類活動干擾小。

表1 2007-2011年安基山水庫上游土壤侵蝕變化統計表

4.1.2土壤侵蝕年際分析

由圖2和表1可知，2007-2011年5 a中2007年土壤侵蝕最為嚴重，強度侵蝕所占面積百分比為0.35%，中度侵蝕所占面積百分比為3.26%，輸沙量最多；2008年土壤侵蝕最弱，未發生強度侵蝕，中度侵蝕所占面積百分比僅為0.15%；2009-2011年土壤侵蝕情況大致相同，其中2010年土壤侵蝕情況稍弱。

降雨條件和人為擾動是影響2007-2011年土壤侵蝕的變化主要因素，其中降雨是水土流失的載體，2007年各月降雨強度相差不大，但其降雨侵蝕力R明顯大于其他年份，以致2007年安基山水庫上游輸沙量最多和強度侵蝕面積較大。2008年降雨量較少，導致降雨侵蝕力下降，土壤侵蝕強度最小。人為擾動因素主要影響植被覆蓋度和水土保持控制因子，發生在人類活動頻繁的區域，但相鄰年間植被變化不明顯。

4.2 水庫淤積分析

根據安基山水庫流域降雨資料和土壤侵蝕結果計算水庫攔沙率、輸沙量和淤積量，結果見表2。2008年水庫攔沙率最大，其他年份攔沙率差別不大。2007年水庫的輸沙量明顯多于其他年份，2008年輸沙量最小，僅為2007年輸沙量的43%。2007年和2011年水庫泥沙淤積量最多，2008年最少，僅為2007年淤積量的49%。降雨是影響水庫泥沙淤積量的重要因素，包括降雨侵蝕因子和年入庫水量，致使攔沙率和輸沙量不同，綜合影響水庫泥沙淤積量。人類活動的干擾也會增大土壤侵蝕強度，增加向下游水庫的泥沙輸送量。

表2 2007-2011年安基山水庫淤積情況

5 結 語

本研究根據通用土壤流失方程USLE利用遙感影像數據、降雨資料、土壤信息和地形資料提取相關影響因子，應用GIS技術實現安基山水庫上游2007-2011年多年土壤侵蝕變化檢測，進而分析上游水土流失對水庫的泥沙淤積影響，為流域土壤侵蝕和水庫淤積治理提供技術支持和參考。其中流域強度侵蝕主要發生在裸地和地勢起伏大的流域中下部；中度和輕度侵蝕主要發生在人類活動頻繁地區和丘陵山區；微度侵蝕則以植被覆蓋較好的林地為主，所占比例均在80%以上；5 a中2007年土壤侵蝕最為嚴重，2008年土壤侵蝕最弱。在水庫淤積方面，2008年水庫攔沙率最大；2007年水庫的輸沙量明顯多于其他年份，2008年輸沙量最小。2007年和2011年水庫泥沙淤積量最多，2008年最少。

根據土壤侵蝕動態檢測分析和統計結果知，流域土壤侵蝕治理對減少水庫淤積至關重要。通過植樹造林、退耕還林、修建截流溝和前置庫、改造坡耕地等措施，減少地表徑流速度，提高地表抗蝕能力，以減少泥沙輸送量，保護水庫蓄水功能和延長使用年限。

因數據資料的精度限制和計算方法的選取不同，其計算結果會與真實值有所偏差，但不影響流域土壤侵蝕的分布趨勢，在土壤侵蝕預測方面仍有指導作用。在今后的應用中，盡可能使用高精度數據并選擇合理的區域性計算方法，以此提高土壤侵蝕預測準確度。

□

[1] 周為峰，吳炳方. 基于遙感和GIS的密云水庫上游土壤侵蝕定量估算[J]. 農業工程學報，2005：21(10)：46-50.

[2] 韓其為，楊小慶. 我國水庫泥沙淤積研究綜述[J]. 中國水利水電科學研究院學報，2003：1(3)：169-175.

[3] 宋曉紅，付 慧，趙 敏. 淺談水庫泥沙淤積量計算的方法[J]. 水利科技與經濟，2005：11(3)：145-149.

[4] 繆馳遠，何丙輝，陳曉燕. “3S”技術在水土保持預測預報模型中的應用[J]. 西南農業大學學報，2004，2(4)：1-4.

[5] 游松財，李 文. GIS支持下的土壤侵蝕量估算——以江西省泰和縣灌溪鄉為例[J]. 自然資源學報，2000，14(2)：19-24.

[6] 焦恩東，于德萬. 水庫泥沙淤積分析計算及防治措施[J]. 吉林水利，2009，(1)：64-66.

[7] 王彥芳，馮民權，尹宏偉，等. 汾河水庫泥沙淤積數值模擬[J]. 水資源與工程學報，2011，22(1)：101-104.

[8] 朱海玲，楊曉輝，張學培. 基于DEM的密云水庫上游流域特征提取與分析[J]. 中國水土保持科學， 2013，11(3)：66-71.

[9] 張光輝. 土壤侵蝕模型研究現狀與展望[J]. 水科學進展，2002，13(3)：389-393.

[10] 謝 云，林 燕，張 巖. 通用土壤流失方程的發展與應用[J]. 地理科學進展，2003，22(3)：279-287.

[11] 吳素業. 安徽大別山區降雨侵蝕力簡化算法與時空分布規律[J]. 中國水土保持，1994，(4)：12-13.

[12] 徐小華，雷 聲，張秀平. 三峽庫區土壤侵蝕定量遙感監測[J]. 江西水利科學，2009，35(4)：253-265.

[13] 張宏鳴，楊勤科，劉晴蕊，等. 基于GIS的區域坡度坡長因子提取算法[J]. 計算機工程，2010，36(9)：246-248.

[14] 馬超飛，馬建文. USLE模型中植被覆蓋因子的遙感數據定量估算[J]. 水土保持通報，2001，21(4)： 6-9.

[15] 齊述華，蔣梅鑫，于秀波. 基于遙感和USLE模型評價1995-2005年江西土壤侵蝕[J]. 中國環境科學，2011, 31(7)：1 197-1 203.

[16] SL190-2007，土壤侵蝕分類分級標準[S].

[17] GB50014-2006，室外排水設計規范[S].

[18] 劉旦旦，王 健，尹武君. 天然降雨對黃土坡地土壤侵蝕和養分流失的影響[J]. 節水灌溉，2011,(8):17-20.