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

?

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

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

摘要為了及時(shí)、準(zhǔn)確地掌握水土流失現(xiàn)狀，了解水土流失變化規(guī)律，為區(qū)域水土流失防治確定科學(xué)的思路和目標(biāo)，基于通用土壤流失方程(USLE)，應(yīng)用RS和GIS技術(shù)對(duì)淮南礦區(qū)水土流失進(jìn)行動(dòng)態(tài)監(jiān)測(cè)與分析。結(jié)果表明，1986~2012年，淮南礦區(qū)水土流失中度侵蝕和輕度侵蝕區(qū)面積比例減少，微度侵蝕區(qū)比例增加，礦區(qū)水土流失狀況總體上趨向變好的態(tài)勢(shì)。研究結(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

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

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

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

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ū)地跨淮河兩岸，東西長(zhǎng)約70 km，南北傾斜寬約25 km，總面積約1 571 km2。整個(gè)礦區(qū)位于江淮丘陵與黃淮平原的交界處，地貌類型兼有平原和丘陵的特點(diǎn)，地形較為復(fù)雜，山丘、崗地、平原、湖洼、河流湖庫(kù)兼而有之?；春佑晌飨驏|橫穿礦區(qū)，北岸為地勢(shì)平坦的淮北平原，南岸為丘陵地區(qū)。礦區(qū)煤炭?jī)?chǔ)量豐富，總儲(chǔ)量約占安徽省總儲(chǔ)量的74%，占華東地區(qū)煤炭?jī)?chǔ)量的45%，被列為全國(guó) 13 個(gè)億噸煤炭生產(chǎn)基地和6 個(gè)煤電基地之一，是我國(guó)東部和南部地區(qū)資源最好、儲(chǔ)量最大的整裝煤田。

1.2數(shù)據(jù)來(lái)源該研究采用RS和GIS技術(shù)，對(duì)20世紀(jì)80、90年代和近期3個(gè)時(shí)段的淮南礦區(qū)GIS數(shù)據(jù)進(jìn)行疊加分析。3個(gè)時(shí)段空間數(shù)據(jù)分別來(lái)自于20世紀(jì)80年代中期淮南礦區(qū)數(shù)字地圖(環(huán)保系統(tǒng)生態(tài)環(huán)境現(xiàn)狀調(diào)查對(duì)比數(shù)據(jù))、20世紀(jì)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ù)(通過(guò)遙感解譯、現(xiàn)場(chǎng)校核得到相匹配的空間數(shù)據(jù))。

2研究方法

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

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

(1)

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

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

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

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

R=-13.86+0.179 2×P

(2)

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

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

2.2.2土壤可蝕性因子K值。K因子反映了土壤對(duì)侵蝕的敏感性及降水所產(chǎn)生的徑流量與徑流速率的大小。K值的大小與土壤質(zhì)地、土壤有機(jī)質(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坡度坡長(zhǎng)因子LS值。區(qū)域范圍內(nèi)坡長(zhǎng)L與坡度S組成的復(fù)合參數(shù)LS，它代表地形條件變化產(chǎn)生侵蝕的主要水力因素。該研究充分利用ArcGIS軟件的表面分析功能，在已生成的研究區(qū)數(shù)字高程模型(DEM)基礎(chǔ)上，利用軟件表面分析功能中的“坡度”分析命令得到礦區(qū)坡度S值因子圖。由于區(qū)域DEM的柵格大小為10 m×10 m，故具體到一個(gè)柵格(圖3)，易得：

L=10/Cosθ

(3)

式中,θ為坡度角。

因此，可通過(guò)軟件的“柵格計(jì)算器”得到坡長(zhǎng)L值因子圖。當(dāng)坡長(zhǎng)和坡度確定以后LS的計(jì)算借鑒施為光[10]四川清平水庫(kù)流域非點(diǎn)源污染負(fù)荷計(jì)算公式：

(4)

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

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

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

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

3結(jié)果與分析

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

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

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

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

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

(5)

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

運(yùn)用公式(5)，計(jì)算得到淮南煤礦3個(gè)時(shí)期的土地退化指數(shù)分別為X1986=0.042，X1999=0.041，X2012=0.034，從1986~2012年逐漸減小，所以可得出淮南礦區(qū)3個(gè)時(shí)期土地退化指數(shù)在減少，說(shuō)明研究時(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個(gè)時(shí)期內(nèi)的水土流失都以輕度侵蝕為主，1986年、1999年和2012年輕度侵蝕的面積分別占總面積的83.03%、80.92%、68.36%。

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

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

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

參考文獻(xiàn)

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

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

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

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

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

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

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

[8] 劉文龍，王堅(jiān)，趙小平.利用GIS進(jìn)行忻州窯礦區(qū)水土流失評(píng)估[J].測(cè)繪通報(bào)，2014(5)：107-109．

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

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

[11] 水利部水土保持司.土壤侵蝕分類分級(jí)標(biāo)準(zhǔn)：SL 190-96[S].北京：中國(guó)水利水電出版社，1997.

收稿日期2015-10-28

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

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

中圖分類號(hào)S 127

文獻(xiàn)標(biāo)識(shí)碼A

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