王曉輝，鄭優(yōu)男　(.安徽省環(huán)境科學(xué)研究院，安徽合肥 300；.合肥工業(yè)大學(xué)資源與環(huán)境工程學(xué)院，安徽合肥 30009)

?

基于RS和GIS技術(shù)的水土流失動態(tài)監(jiān)測與分析——以淮南礦區(qū)為例

王曉輝1，鄭優(yōu)男2(1.安徽省環(huán)境科學(xué)研究院，安徽合肥 230022；2.合肥工業(yè)大學(xué)資源與環(huán)境工程學(xué)院，安徽合肥 230009)

摘要為了及時、準確地掌握水土流失現(xiàn)狀，了解水土流失變化規(guī)律，為區(qū)域水土流失防治確定科學(xué)的思路和目標，基于通用土壤流失方程(USLE)，應(yīng)用RS和GIS技術(shù)對淮南礦區(qū)水土流失進行動態(tài)監(jiān)測與分析。結(jié)果表明，1986~2012年，淮南礦區(qū)水土流失中度侵蝕和輕度侵蝕區(qū)面積比例減少，微度侵蝕區(qū)比例增加，礦區(qū)水土流失狀況總體上趨向變好的態(tài)勢。研究結(jié)果有助于指導(dǎo)淮南礦區(qū)水土流失防治及其生態(tài)系統(tǒng)的恢復(fù)與重建。

關(guān)鍵詞淮南礦區(qū)；水土流失；土壤流失方程

Dynamic Monitoring and Analysis of Water and Soil Erosion Based on RS and GIS Technology—Taking Huainan Mining Area as an Example

WANG Xiao-hui1, ZHENG You-nan2(1. Anhui Institute of Environmental Science, Hefei, Anhui 230022; 2. School of Nature Resources and Environmental Engineering, Hefei University of Technology, Hefei, Anhui 230009)

AbstractIn order to timely and accurately grasp water and soil erosion situation, understand the variation law, provide scientific thought and objective of regional water and soil erosion control, based on USLE, RS and GIS technology was applied to monitor and analyze water and soil erosion in Huainan mining area. The results showed that during 1986-2012, proportion of moderate erosion and mild erosion area decrease, micro erosion area proportion increase. The results are helpful to guide the prevention and control of soil erosion in Huainan mining area and the restoration and reconstruction of ecological system.

Key wordsHuainan mining area; Water and soil loss; USLE

水土流失是我國生態(tài)環(huán)境的主要破壞因素之一。開展水土流失動態(tài)監(jiān)測與分析，及時掌握水土流失變化特點，對加強水土流失的預(yù)防、保護、監(jiān)督和管理，以及有效減少和控制水土流失的發(fā)生和發(fā)展具有重要的意義。水土流失動態(tài)監(jiān)測與分析主要是指對區(qū)域水土流失的面積、程度和分布狀況實施長期、連續(xù)的監(jiān)測，需要能及時獲取現(xiàn)勢性的資料并對現(xiàn)有資料進行快速處理分析的新技術(shù)手段。衛(wèi)星遙感技術(shù)(RS)的快速更新信息能力和地理信息系統(tǒng)(GIS)強大的數(shù)據(jù)處理和分析能力，已成為近年來水土保持領(lǐng)域倍受關(guān)注的研究熱點之一。

我國礦區(qū)在資源開采過程中伴隨著棄土廢渣的排放、地表塌陷、破壞地下水資源的平衡等現(xiàn)象，導(dǎo)致地表植被的破壞，使得地表裸露，在降雨影響下，會造成嚴重的水土流失。但目前國內(nèi)外針對于礦區(qū)的研究主要集中于礦區(qū)的生態(tài)環(huán)境質(zhì)量評價、污染治理和污染現(xiàn)狀評價及預(yù)警等[1-4]，且研究技術(shù)日益趨于成熟，相比之下礦區(qū)水土流失方面的研究略顯薄弱[5-8]。筆者利用RS和GIS技術(shù)對淮南煤礦區(qū)多時期水土流失進行動態(tài)監(jiān)測和分級評估，以便能及時掌握礦區(qū)水土流失變化情況，從而為有效防治和治理提供科學(xué)支撐。

1研究區(qū)概況及數(shù)據(jù)來源

1.1研究區(qū)概況淮南礦區(qū)地處安徽省中北部，淮河中游，位于116°21′21″～117°11′59″E、32°32′45″～ 33°00′24″ N，海拔16.5～240 m，自西北向東南，坡降1/10 000。礦區(qū)地跨淮河兩岸，東西長約70 km，南北傾斜寬約25 km，總面積約1 571 km2。整個礦區(qū)位于江淮丘陵與黃淮平原的交界處，地貌類型兼有平原和丘陵的特點，地形較為復(fù)雜，山丘、崗地、平原、湖洼、河流湖庫兼而有之。淮河由西向東橫穿礦區(qū)，北岸為地勢平坦的淮北平原，南岸為丘陵地區(qū)。礦區(qū)煤炭儲量豐富，總儲量約占安徽省總儲量的74%，占華東地區(qū)煤炭儲量的45%，被列為全國 13 個億噸煤炭生產(chǎn)基地和6 個煤電基地之一，是我國東部和南部地區(qū)資源最好、儲量最大的整裝煤田。

1.2數(shù)據(jù)來源該研究采用RS和GIS技術(shù)，對20世紀80、90年代和近期3個時段的淮南礦區(qū)GIS數(shù)據(jù)進行疊加分析。3個時段空間數(shù)據(jù)分別來自于20世紀80年代中期淮南礦區(qū)數(shù)字地圖(環(huán)保系統(tǒng)生態(tài)環(huán)境現(xiàn)狀調(diào)查對比數(shù)據(jù))、20世紀90年代末數(shù)字地圖(環(huán)保系統(tǒng)生態(tài)環(huán)境現(xiàn)狀調(diào)查數(shù)據(jù))和2012年淮南礦區(qū)SPOT6的1.5 m全色及6 m多光譜捆綁數(shù)據(jù)(通過遙感解譯、現(xiàn)場校核得到相匹配的空間數(shù)據(jù))。

2研究方法

2.1土壤流失方程及監(jiān)測流程水土流失監(jiān)測手段包括常規(guī)地面調(diào)查法和定性遙感法，而定性遙感法以其快速、簡單、宏觀等特點而在國內(nèi)外得到廣泛應(yīng)用，并建立了眾多的遙感侵蝕監(jiān)測模型。其中，從數(shù)據(jù)適用性和應(yīng)用范圍而言，以Wischmeier等[9]于20世紀60年代提出通用土壤流失方程(USLE)應(yīng)用為佳。國內(nèi)外學(xué)者做了大量的水土流失研究，包括對方程中各因子的計算和各因子參數(shù)的本地化，使該方程的適用范圍和監(jiān)測精度得到進一步提高。由于該方程的經(jīng)典性和因子的解釋具有物理意義，仍是目前監(jiān)測土壤侵蝕最為廣泛使用的方法。該方程是將水土流失的6個因子用連乘的形式組成，表達式為：

A=f×R×K×LS×C×P

(1)

式中,A為年土壤流失量(t/km2·a)；f為單位轉(zhuǎn)換系數(shù)，f=224.2；R為降雨和徑流因子；K為土壤可蝕性因子；LS為坡度坡長因子；C為植被與經(jīng)營管理因子；P為水土保持因子。

運用ArcGIS的空間數(shù)據(jù)管理分析功能，建立研究區(qū)的數(shù)字高程模型(DEM)以及土地利用現(xiàn)狀圖、土壤類型圖、植被覆蓋分布圖等矢量數(shù)據(jù)圖，對其屬性數(shù)據(jù)進行相應(yīng)的數(shù)據(jù)編碼操作，并將其柵格化(Grid格式)，求得各參數(shù)值的因子圖。再根據(jù)USLE方程的形式，將各因子值相乘，獲得研究區(qū)土壤侵蝕強度等級圖，其工作流程見圖1。

2.2參數(shù)因子的計算

2.2.1降雨侵蝕因子R值。降雨侵蝕力因子是評價降雨引起土壤分離和搬運的潛在動力指標，難以直接測定，大多用降雨參數(shù)，如雨強、雨量等來估算降雨侵蝕力。該研究認為采用福建農(nóng)業(yè)大學(xué)提出的R值計算公式較為合適，計算公式為：

R=-13.86+0.179 2×P

(2)

式中，R為全年的降雨侵蝕力；P為流域年降雨量(mm)。

通過淮南區(qū)氣象站的降雨量資料1986、1999和2012年3個年份區(qū)域降雨量分別為773、785和797 mm)，整個研究區(qū)每年的R值只有一個，即R1986=124.66，R1999=126.81，R2012=128.96。

2.2.2土壤可蝕性因子K值。K因子反映了土壤對侵蝕的敏感性及降水所產(chǎn)生的徑流量與徑流速率的大小。K值的大小與土壤質(zhì)地、土壤有機質(zhì)含量有較高的相關(guān)性。根據(jù)二次土壤普查成果，查找Wischmeier等建立的土壤可蝕性K值諾謨圖，近似確定出研究區(qū)不同土壤類型的可蝕性因子K值(表1)。

表1　淮南礦區(qū)土壤可蝕性因子K值估算結(jié)果

將研究區(qū)現(xiàn)有17類土壤類型經(jīng)ArcGIS軟件的“數(shù)據(jù)融合”命令處理并編輯轉(zhuǎn)化為9類易于K因子賦值的土壤質(zhì)地類型，按照表1將其增加一列屬性K并賦相應(yīng)數(shù)值，最后利用 “矢量數(shù)據(jù)柵格化”命令生成K值因子?xùn)鸥駡D(圖2)。

2.2.3坡度坡長因子LS值。區(qū)域范圍內(nèi)坡長L與坡度S組成的復(fù)合參數(shù)LS，它代表地形條件變化產(chǎn)生侵蝕的主要水力因素。該研究充分利用ArcGIS軟件的表面分析功能，在已生成的研究區(qū)數(shù)字高程模型(DEM)基礎(chǔ)上，利用軟件表面分析功能中的“坡度”分析命令得到礦區(qū)坡度S值因子圖。由于區(qū)域DEM的柵格大小為10 m×10 m，故具體到一個柵格(圖3)，易得：

L=10/Cosθ

(3)

式中,θ為坡度角。

因此，可通過軟件的“柵格計算器”得到坡長L值因子圖。當坡長和坡度確定以后LS的計算借鑒施為光[10]四川清平水庫流域非點源污染負荷計算公式：

(4)

同樣，利用“柵格計算器”計算式(4)便可得到研究區(qū)坡長坡度LS值因子?xùn)鸥駡D(圖4)。

2.2.4植被覆蓋因子C值。植被覆蓋因子反映了植被對地表的保護作用。在完全沒有植被保護的裸露地面C值定為1，完全被植被覆蓋的地面C值可取0。顯然，C的取值范圍為0~1。其值大小取決于具體的植被覆蓋、輪作順序及管理措施的綜合作用等。也就是說，C值的取值主要與植被覆蓋和土地利用類型有關(guān)。根據(jù)淮南礦區(qū)輪作制度的調(diào)查和相關(guān)研究報道，利用植被分布矢量圖、輪作方式及覆蓋度更能反映C因子值。根據(jù)植被分布圖中作物種類及覆蓋度情況確定相應(yīng)C值(見表2)，建立C因子的屬性連接，并利用“矢量數(shù)據(jù)柵格化”命令生成不同年度C值因子圖(以2012年為例，見圖5)。

表2　淮南礦區(qū)C因子取值

2.2.5水土保持因子P值。水土保持因子P又稱侵蝕防治措施因子，是指采用專門措施后的土壤流失量與順坡種植時的土壤流失量的比值，通常的控制措施有等高耕作等。坡面侵蝕量試驗研究結(jié)果表明，如以自然植被和坡耕地P因子為1，則灌木叢為0.8，農(nóng)耕地為0.35。將前面整理好的6類流域土地利用類型屬性項與P值記錄項建立鏈接，將P值按表3賦值給相應(yīng)土地利用類型并轉(zhuǎn)化成柵格圖形式，得到淮南礦區(qū)不同年度水土保持P值因子圖(以2012年為例，見圖6)。

3結(jié)果與分析

根據(jù)式(1)，以及上述所計算的各因子圖層，在ArcGIS 中使用柵格計算器，將各因子圖層相乘，即可得到淮南礦區(qū)水土流失分布圖(圖7～9)。又根據(jù)水利部土壤侵蝕面蝕分級標準[11]，把研究區(qū)水土流失狀況按水土流失強度分為6級，即微度(小于200 t·km-2·a-1)、輕度(200～2 500 t·km-2·a-1)、中度(2 500～5 000 t·km-2·a-1)、強度(5 000～8 000 t·km-2·a-1)、極強度(8 000～15 000 t·km-2·a-1)和劇烈(大于15 000 t·km-2·a-1)。由此得到淮南礦區(qū)不同年份水土流失程度各等級面積(表3)。

表3　淮南礦區(qū)年際間水土流失程度各等級面積變化

由表3可以看出，各研究區(qū)水土流失程度等級區(qū)域面積中以輕度侵蝕面積為最，其次是微度侵蝕，中度侵蝕面積最少。其中，微度侵蝕強度面積由大到小依次為2012、1999、1986年；輕度侵蝕強度面積由大到小依次為1986、1999、2012年；而中度侵蝕強度面積也以2012年為最小?？梢?，2000年以后，淮南礦區(qū)整體水土流失程度持續(xù)減小，這與其不斷加強生態(tài)保護與恢復(fù)，大力開展植樹造林活動是一致的。

1986、1999、2012年由3個時期的土壤流失分布可以看出，淮南礦區(qū)土地流失主要以輕度侵蝕為主，但1986~2012年間輕度侵蝕的面積在逐年遞減，相反的，微度侵蝕的面積在逐年增加；而中度侵蝕的面積在1999~2012年也有著顯著的縮減，同時其所處位置也從火石山、金家?guī)X一帶轉(zhuǎn)移到了平阿山附近。

根據(jù)上述所得數(shù)據(jù)計算淮南煤礦3個時期的土地退化指數(shù)，按照2006 年生態(tài)環(huán)境狀況評價技術(shù)規(guī)范[12]，計算公式為：

(5)

式中，X為土地退化指數(shù)；Di為各侵蝕類型土地面積；Wi為各侵蝕類型的權(quán)重，輕度為0.05、中度為0.25、重度為0.7；D為研究區(qū)內(nèi)總面積。

運用公式(5)，計算得到淮南煤礦3個時期的土地退化指數(shù)分別為X1986=0.042，X1999=0.041，X2012=0.034，從1986~2012年逐漸減小，所以可得出淮南礦區(qū)3個時期土地退化指數(shù)在減少，說明研究時段內(nèi)淮南礦區(qū)水土流失情況得到了很好的控制和改善。

4結(jié)論

該研究根據(jù)淮南礦區(qū)1986、1999和2012年3期的遙感影像數(shù)據(jù)，并結(jié)合數(shù)字地圖等數(shù)據(jù)，以RS和GIS技術(shù)為支持，生成淮南礦區(qū)1986、1999和2012年的水土流失數(shù)據(jù)并加以分析，得出以下結(jié)論：

(1)淮南礦區(qū)3個時期內(nèi)的水土流失都以輕度侵蝕為主，1986年、1999年和2012年輕度侵蝕的面積分別占總面積的83.03%、80.92%、68.36%。

(2)淮南礦區(qū)研究期間輕度和中度侵蝕的面積在逐年減少，而微度侵蝕的面積在遞增，2012年礦區(qū)中度侵蝕的面積已減少至0.6km2，中度侵蝕對的面積基本已達消失。由此可得淮南礦區(qū)3個時期內(nèi)水土流失程度在減輕，這與其加強生態(tài)保護和大力植樹造林的方針是密不可分的。

(3)淮南礦區(qū)3個時期的土地退化指數(shù)在逐年遞減，2012年淮南礦區(qū)的土地退化指數(shù)降低至0.034，說明礦區(qū)的水土流失情況得到了有效的控制和改善。

所以為了綜合治理礦區(qū)的水土流失情況，應(yīng)加強礦區(qū)的生態(tài)保護與恢復(fù)，大力開展退耕還林活動，復(fù)墾塌陷區(qū)土地為林地、草地等；同時對礦區(qū)進行水土流失動態(tài)監(jiān)測是必不可少的，根據(jù)監(jiān)測數(shù)據(jù)更好地進行下一步治理，更有效地控制礦區(qū)水土流失，改善礦區(qū)生態(tài)環(huán)境。

參考文獻

[1] 郭金陵，於世為，熊慧，等.基于AHP的煤炭礦區(qū)生態(tài)環(huán)境質(zhì)量模糊綜合評價[J].科技管理研究，2011(15)：50-53.

[2] 徐嘉興，李 鋼，陳國良，等.土地復(fù)墾礦區(qū)的景觀生態(tài)質(zhì)量變化[J].農(nóng)業(yè)工程學(xué)報，2013，29(1)：232-239.

[3]VATALISKI，KALIAMPAKOSDC.Anoverallindexofenvironmentalqualityincoalminingareasandenergyfacilities[J].Environmanage，2006，38：1031-1045.

[4] 李保杰，顧和和，紀亞洲.礦區(qū)土地復(fù)墾景觀格局變化和生態(tài)效應(yīng)[J].農(nóng)業(yè)工程學(xué)報，2012，28(3)：251-256.

[5] 劉會玉，王充，林振山，等.基于RS和GIS的滇池流域水土流失動態(tài)監(jiān)測[J].南京師大學(xué)報(自然科學(xué)版)，2012，35(2)：120-124.

[6] 張明陽，王克林，陳洪松.基于RS和GIS的喀斯特區(qū)域水土流失動態(tài)監(jiān)測與分析：以廣西環(huán)江縣為例[J].資源科學(xué)，2007，29(3)：125-131.

[7] 李曉松，姬翠翠，曾源，等.基于遙感和GIS的水土流失動態(tài)監(jiān)測：以河北省赤城縣為例[J].生態(tài)學(xué)雜志，2009，28(9)：1723-1729.

[8] 劉文龍，王堅，趙小平.利用GIS進行忻州窯礦區(qū)水土流失評估[J].測繪通報，2014(5)：107-109．

[9]WISCHMEIERWH，SMITHDD．Predictingrainfallerosionlossesfromcroplandeastoftherockymountains：Guideforselectionofpracticesforsoilandwaterconservation[M]．Washington，DC．：AgriculturalResearchService，USDepartmentofAgriculture，1965．

[10] 施為光.四川省清平水庫流域非點源污染負荷計算[J].重慶環(huán)境科學(xué)，2000，22(2)：33-36.

[11] 水利部水土保持司.土壤侵蝕分類分級標準：SL 190-96[S].北京：中國水利水電出版社，1997.

收稿日期2015-10-28

作者簡介王曉輝(1972- )，男，安徽合肥人，正高級工程師，碩士，從事生態(tài)評價與規(guī)劃研究。

基金項目淮南礦業(yè)集團科研項目[HNKY JT JS (2013)31]。

中圖分類號S 127

文獻標識碼A

文章編號0517-6611(2015)33-373-04