尹　璐，閆慶武，卞正富

(中國礦業(yè)大學環(huán)境與測繪學院，江蘇 徐州　221116)

?

基于RUSLE模型的六盤水市土壤侵蝕評價

尹璐，閆慶武①，卞正富

(中國礦業(yè)大學環(huán)境與測繪學院，江蘇 徐州221116)

摘要：礦產(chǎn)資源的開采容易導致礦區(qū)土地資源破壞及生態(tài)環(huán)境惡化，為探索煤礦開采對土壤侵蝕的影響，分析礦區(qū)土地侵蝕狀況。以貴州省西南部煤礦城市六盤水市為研究對象，基于土壤質(zhì)地、地形、氣候、土地覆蓋和土地利用等數(shù)據(jù)，運用RUSLE模型，結合GIS空間分析方法，研究六盤水市土壤侵蝕現(xiàn)狀及空間分布特征，著重分析煤礦區(qū)土壤侵蝕分布特征。結果表明：(1)六盤水市土壤侵蝕以微度、中度侵蝕為主，侵蝕嚴重區(qū)域主要集中在水城縣東部和南部，六枝特區(qū)東部以及盤縣中部。(2)土壤侵蝕主要發(fā)生在人類負向干擾活動強烈的林地、草地、旱地、煤礦區(qū)和裸地。(3)六盤水市南部煤礦區(qū)侵蝕較嚴重，以中度侵蝕為主，主要原因是采煤造成的地質(zhì)災害使土壤侵蝕抑制因素作用減弱。(4)私營煤礦區(qū)土壤侵蝕程度比國有煤礦區(qū)嚴重。私營煤礦區(qū)應在煤炭開采過程中加強保護與治理。

關鍵詞：GIS；水土流失；RULSE模型；煤礦區(qū)

水土流失是一個全球性的土地退化問題,是當今全球變化的重要表現(xiàn)形式[1]。我國水土流失與生態(tài)安全綜合科學考察活動研究表明,我國現(xiàn)有土壤侵蝕面積為356.92萬km2,其中水力侵蝕面積為161.22萬km2,風力侵蝕面積為195.70萬km2[2]。定量評價區(qū)域土壤侵蝕量、土壤侵蝕強度及其空間分布特征,對水土保持和耕地保護具有重要意義[3]。就煤礦區(qū)而言,隨著煤炭開采活動的不斷深入,對土地資源的破壞日益加劇,煤礦的采、運、堆等過程使得植被遭到嚴重破壞,地表裸露,土壤表層松散性加大,喪失了原地表土壤的抗蝕力,造成土地劣化、貧瘠化和干旱化,不僅使耕地急劇減少,還會造成地表水源和地下含水層水源的漏失,引起礦區(qū)水土流失[4]。

目前,土壤侵蝕研究應用最廣泛的模型是修正通用水土流失方程(revisied universal soil loss equation,RUSLE),該模型也被成熟地運用到我國土壤侵蝕調(diào)查工作中。劉寶元等[5]借鑒通用土壤流失方程(universal soil loss equation,USLE)的成果經(jīng)驗,建立適用于中國陡坡的土壤侵蝕預報模型(Chinese soil loss equation,CSLE);肖洋等[6]通過查表法,構建不同植被覆蓋度和不同土地覆蓋類型下的地表植被覆蓋因子C值圖層,估算了重慶地區(qū)土壤侵蝕量、土壤保持量和土壤保持價值;田鵬等[7]以RUSLE模型為基礎,增加流域輸沙能力與淤地壩攔沙效率模塊,并在黃土高原粗沙多沙區(qū)皇甫川流域進行驗證。筆者在GIS和RS技術的支持下,結合修訂的RUSLE模型對研究區(qū)域及研究區(qū)內(nèi)煤礦區(qū)土壤侵蝕現(xiàn)狀進行定量分析,為促進煤礦區(qū)土地復墾、土壤改善和土地資源利用率的提高提供科學依據(jù)。

1研究區(qū)概況

六盤水市位于貴州省西部,云貴高原一、二級臺地斜坡上,地理位置為25°18′～26°50′ N,104°18′～105°53′ E,現(xiàn)轄六枝特區(qū)、盤縣、水城縣和鐘山區(qū)4個縣級行政區(qū),總面積為9 914 km2,占全省總面積的5.63%,2013年常住人口為285.9萬,GDP總量為876.95億元,人均GDP全省排名第2。六盤水市屬亞熱帶高原性季風氣候區(qū),氣候溫和,夏無酷暑,冬無嚴寒,年平均氣溫為13～14 ℃,被稱為中國涼都,境內(nèi)降水充沛,但時空分布不均,年降水量表現(xiàn)為南部多北部偏少,東部多西部偏少,旱澇災害頻繁,多年平均降水量為1 244 mm。六盤水市地處長江水系和珠江水系的分水嶺地區(qū)。市境礦藏儲量豐富,礦種較多,煤儲量居貴州省之首,素有“江南煤都”之譽,是國家“西電東送”的重點電源基地之一。2013年5月發(fā)布的《貴州省水土保持公告》顯示,2010年貴州全省水土流失面積為55 269.40 km2,占土地總面積的31.37%,全省水土流失狀況呈現(xiàn)西北重東南輕的格局,六盤水市地處生態(tài)環(huán)境脆弱的喀斯特地區(qū),水土流失治理難度大。

2研究方法

2.1數(shù)據(jù)資料

為開展六盤水市水土流失評價,需收集以下6類數(shù)據(jù):

(1)土壤質(zhì)地數(shù)據(jù)。一部分來源于中國科學院資源環(huán)境科學數(shù)據(jù)中心,數(shù)據(jù)包括砂粒、粉粒和黏粒的組成;另一部分來源于六盤水市國土部門收集的土壤類型數(shù)據(jù)及土壤采樣點數(shù)據(jù),通過插值生成的研究區(qū)域土壤有機碳分布圖。

(2)地形數(shù)據(jù)。運用ASTER GDEM 30 m分辨率的數(shù)字高程模型(DEM),利用ArcGIS的表面分析方法計算坡度和坡長。

(3)氣候數(shù)據(jù)。通過中國氣象數(shù)據(jù)網(wǎng)獲取六盤水市及其周圍15個氣象站點2004—2014年各月平均降水量,并用克里金插值方法生成30 m分辨率的多年平均降水量分布圖和多年平均各月降水量分布圖。

(4)土地覆蓋數(shù)據(jù)。利用2013年6月16日Landsat 8 OLI遙感影像,經(jīng)拼接、裁剪、配準和大氣校正等預處理后進行數(shù)據(jù)反演,得到研究區(qū)域植被覆蓋度分布圖。

(5)土地利用數(shù)據(jù)。采用國土部門收集的2010年第2次全國土地利用矢量調(diào)查數(shù)據(jù),比例尺為1∶25萬。對土地利用數(shù)據(jù)進行柵格化,與計算出的水土流失等級圖進行疊加分析,得到不同土地利用類型的水土流失分布數(shù)據(jù)。

(6)六盤水市煤礦區(qū)數(shù)據(jù)。采用國土部門收集的《六盤水市礦產(chǎn)資源總體規(guī)劃(2011—2015)》中的煤礦區(qū)拐點坐標生成礦區(qū)位置矢量圖。將六盤水市分為南、北和東3個煤礦區(qū)域。將煤炭生產(chǎn)能力大于120萬t·a-1的國有企業(yè)煤礦區(qū)作為國有煤礦區(qū)土壤侵蝕現(xiàn)狀的重點研究區(qū)域,將與上述國有煤礦研究區(qū)域有著相近環(huán)境和地理位置、生產(chǎn)能力為30萬～60萬t·a-1的私營煤礦區(qū)作為私營煤礦區(qū)土壤侵蝕現(xiàn)狀的重點研究區(qū)域。

將柵格數(shù)據(jù)統(tǒng)一為30 m×30 m大小的柵格單元,統(tǒng)一采用WGS-1984-UTM投影。

2.2土壤侵蝕模型

采用修正的RUSLE模型對六盤水市土壤侵蝕量進行估算。RUSLE是對USLE加以改進發(fā)展而來的經(jīng)驗土壤侵蝕預報模型,得到國內(nèi)外的廣泛研究和應用,目前已有學者驗證了RUSLE模型在中國南方喀斯特地貌的應用是可行的[8-9],因此選用RUSLE模型作為研究土壤侵蝕的定量模型,其表達式[10]為

A=R×K×L×S×C×P。

(1)

式(1)中,A為單位面積上時間和空間的平均土壤流失量,t·km-2·a-1;R為降雨侵蝕力因子,MJ·mm·hm-2·h-1·a-1;K為土壤可蝕性因子,t·h·MJ-1·mm-1;L為坡長因子;S為坡度因子，(°);C為地表植被覆蓋因子;P為水土保持措施因子。

2.3研究區(qū)RUSLE因子計算

2.3.1降雨侵蝕力因子(R)估算

R反映了降雨對土壤的潛在侵蝕能力,與降雨量、降雨強度和降雨歷時等有關。采用WISCHMEIER等[11]提出的R值簡便算法,即利用多年月降水量建立的年均降雨侵蝕力因子的經(jīng)驗公式:

(2)

式(2)中,p為多年平均降水量,mm;pi為各月平均降水量,mm。

2.3.2土壤可蝕性因子(K)估算

(3)

式(3)中,w砂為砂粒(粒徑為0.1～2 mm)含量,%;w粉為粉粒(粒徑為0.002～<0.1 mm)含量,%,w黏為黏粒(粒徑<0.002 mm)含量,%;w碳為土壤有機碳含量,%。N=1-w砂/100。由EPIC模型計算出的KEPIC值與實測值具有良好的線性關系,可通過修正轉換公式對KEPIC進行修訂,以用于我國的土壤可蝕性估算[13]。修訂公式為

K=-0.013 83+0.515 75KEPIC。

(4)

利用式(4)得到每個像元單位的K值,用面積加權法求得土壤亞類K值(表1)。

表1六盤水市各土壤類型的可蝕型因子(K)取值

Table 1Valuation of erodible factor(K) of various types of soils in Liupanshui City

土壤類型K值土壤類型K值粗骨土0.110漂洗黃壤0.138黑色石灰土0.136石灰(巖)土0.135紅壤0.116水稻土0.131紅色石灰土0.128酸性石質(zhì)土0.105黃紅壤0.140酸性紫色土0.130黃壤0.126巖石0.128黃壤性土0.135中性紫色土0.119黃色石灰土0.122紫色土0.124黃棕壤0.130棕壤0.132黃棕壤性土0.112棕色石灰土0.141灰潮土0.122

2.3.3坡長(L)和坡度(S)因子估算

坡長和坡度因子反映了地形地貌對降雨侵蝕動力的加速影響,表示在降雨強度、植被覆蓋度等條件相同的情況下,坡度越大,坡長越長,徑流速度越大,對土壤的沖刷能力越強[14]。使用Hickey和Van Remortel構建的方法計算L和S因子,計算公式[15]為

(5)

β=(sinθ/0.089 6)/[3(sinθ)0.8+0.56],

(6)

m=β/(1+β),

(7)

(8)

式(5)～(8)中,λ為像元坡長,m;m為坡長指數(shù)；θ為坡度。

2.3.4地表植被覆蓋因子(C)估算

植被覆蓋多少對土壤有著重要影響,C是土壤侵蝕模型中的抑制因子,C值與植被覆蓋度(VFC,CVF)有著較強的相關性。用Landsat 8 OLI遙感影像計算六盤水市歸一化植被指數(shù)(normalized differential vegetation index,NDVI,INDV)值,計算公式為

(9)

式(9)中,RNI為近紅外波段(0.7～1.1 μm)像元亮度值;R為紅光波段(0.4～<0.7 μm)像元亮度值。

選取INDV值累積概率為5%左右數(shù)為區(qū)域內(nèi)最小INDV值(INDV,min),選取INDV值累積概率為95%左右數(shù)為最大INDV值(INDV,max),求得INDV,min=0.288 868,INDV,max=0.829 726。利用INDV值計算六盤水市CVF,計算公式為

利用建立的CVF與C因子之間的關系計算C值,計算公式[16]為

C=1-CVF。

2.3.5水土保持措施因子(P)估算

P是土壤侵蝕模型中的抑制因子。當P=0時,表示采取水土保持措施后不發(fā)生土壤侵蝕;當P=1時,表示未采取任何水土保持措施[17]。根據(jù)六盤水市土地利用現(xiàn)狀和實地調(diào)查情況,結合前人的研究成果[7,18-20],確定水田、旱地、園地、有林地、其他林地、灌木林地、草地、建設用地、采礦用地、水域和裸地的P值分別為0.01、0.4、0.4、1、0.9、1、0.8、0、0.5、0和1,并計算得到P值分布圖(圖1)。

圖1　六盤水市各土壤侵蝕因子計算結果

根據(jù)SL 190—2007《土壤侵蝕分類分級標準》[21],六盤水市屬于水力侵蝕區(qū)的西南土石山區(qū),該類型侵蝕區(qū)的容許土壤流失量為500 t·km-2·a-1,應用RUSLE模型計算六盤水市土壤侵蝕強度,并將其分為6級。

3結果與分析

3.1土壤侵蝕等級評價

通過RUSLE模型計算研究區(qū)土壤侵蝕強度并劃分等級(表2),研究區(qū)超過容許土壤流失量(500 t·km-2·a-1)的土地面積達6 710.06 km2,占全市侵蝕總面積的68.11%,其中,微度、輕度侵蝕面積為7 874.50 km2,占比達79.93%,但其土壤侵蝕量強度占比僅為42.13%;輕度、中度和強烈土壤侵蝕量強度較大,占土壤侵蝕量總強度的85.68%。從圖2～3可知,六盤水市土壤侵蝕嚴重的地區(qū)位于水城縣東部和南部,六枝特區(qū)東部以及盤縣中部,這

些地區(qū)地勢起伏較大,降雨較多,植被對降雨的攔截力度不夠是其受侵蝕的主要因素。

表2研究區(qū)域土壤侵蝕強度分級統(tǒng)計表

Table 2Statistics and grading of soil erosion intensities in the study area

侵蝕等級侵蝕模數(shù)/(t·km-2·a-1)侵蝕面積/km2侵蝕面積占比/%侵蝕量強度占比/%微度<5003141.7231.893.71輕度500～<25004732.7848.0438.41中度2500～<50001408.1414.2931.05強烈5000～<8000411.224.1816.21極強烈8000～<15000147.781.509.45劇烈≥1500010.140.101.17

3.2土壤侵蝕與土地利用的相關分析

不同的土地利用類型對區(qū)域土壤侵蝕的發(fā)生面積和強度特征有著重要影響[22]。六盤水市不同土地利用類型土壤侵蝕強度分布見表3。

圖2　六盤水市土壤侵蝕強度等級分布

由表3可知,水域和建設用地屬于侵蝕極少區(qū)域,水田屬于不易發(fā)生侵蝕區(qū),主要原因是水田多分布在地勢平坦、水熱條件好的地區(qū),且人為管理措施強;園地水土流失較輕,主要為微度和輕度流失,主要原因是人為保護措施較好;六盤水市土壤侵蝕主要發(fā)生在草地、旱地、林地、裸地和煤礦區(qū),其中,灌木林地、有林地和其他林地侵蝕特征較接近,主要為輕度、微度和中度侵蝕,分別占各地類面積的88.80%、91.77%和89.22%;旱地超過土壤侵蝕容許模數(shù)的侵蝕面積占旱地總面積的67.76%,主要發(fā)生輕度侵蝕,占旱地總面積的57.28%;草地超過土壤侵蝕容許模數(shù)的侵蝕面積占比達80.46%,主要發(fā)生輕度和中度侵蝕,分別占草地總面積的43.56%和22.72%;煤礦區(qū)超過土壤侵蝕容許模數(shù)的侵蝕面積占比達73.95%,主要發(fā)生微度和輕度侵蝕,分別占煤礦區(qū)總面積的26.05%和58.19%;裸地是極易發(fā)生土壤侵蝕的土地類型,侵蝕等級以輕度、中度和強烈為主,占裸地總面積的84.25%。

圖3　六盤水市礦區(qū)分布

表3各土地利用類型土壤侵蝕面積分布特征

Table 3Distribution of soil erosion areas relative to land use

地類面積比例/%微度侵蝕輕度侵蝕中度侵蝕強烈侵蝕極強烈侵蝕劇烈侵蝕面積/km2比例/%面積/km2比例/%面積/km2比例/%面積/km2比例/%面積/km2比例/%面積/km2比例/%草地8.57164.4419.54366.5343.56191.1222.7285.0210.1031.513.742.760.34建設用地1.36122.3099.730.250.210.060.050.020.010.000.000.000.00旱地40.261279.1332.242272.7057.28384.879.7029.580.751.080.030.020.00裸地3.4836.8210.73153.2544.6796.3928.0939.4311.4916.064.681.130.34煤礦區(qū)6.89177.6026.05396.6758.1995.1313.9511.451.680.870.130.000.00水田2.97289.0899.750.510.180.140.050.040.010.020.010.000.00水域0.2826.3199.830.030.100.010.050.000.010.000.010.000.00有林地15.45506.7333.25652.0442.77240.1315.7589.565.8834.032.231.890.12其他林地3.6679.0621.93165.3745.8877.1721.4128.337.8610.072.790.480.13灌木林地16.83442.5226.70708.0342.71321.5019.39127.597.7054.123.273.870.23園地0.2511.0245.2111.6347.711.516.200.210.850.010.030.000.00

通過對不同土地利用類型土壤侵蝕狀況的定量分析可知,研究區(qū)土壤侵蝕主要發(fā)生在人類負向干擾活動強烈的林地、草地、旱地、煤礦區(qū)和裸地,草地和裸地是最易發(fā)生強烈以上等級侵蝕的土地利用類型,因此要加強對六盤水市草地和裸地的防護和治理工作,尤其是要防止草地和裸地石漠化演變,一旦石漠化將不易恢復。林地多處于坡度大的地區(qū),遇到暴雨沖刷時,也會引起輕度和中度土壤侵蝕。旱地和煤礦區(qū)屬于人類影響作用強烈的地區(qū),農(nóng)牧區(qū)和工礦活動區(qū)內(nèi)人類不合理的生產(chǎn)方式是導致侵蝕發(fā)生的主要原因,因此要重視農(nóng)牧區(qū)和工礦區(qū)的水土保持工作,如采取等高耕作、平衡施肥和邊開采邊治理等措施。

3.3煤礦區(qū)土壤侵蝕的空間分布

3.3.1煤礦區(qū)分區(qū)土壤侵蝕

為了反映六盤水市煤礦區(qū)土壤侵蝕狀況,將六盤水市分為南、北和東3個煤礦區(qū)域(圖3～4),以鐘山區(qū)和水城縣134個煤礦為北區(qū),以盤縣101個煤礦為南區(qū),以六枝特區(qū)40個煤礦為東區(qū)(表4～5)。

圖4　南、北和東煤礦區(qū)土壤侵蝕強度分級

從煤礦區(qū)總體情況來看,煤礦區(qū)極強烈和劇烈侵蝕不明顯,以微度、輕度和中度侵蝕面積占比較大,達98.17%,煤礦區(qū)中度侵蝕面積比例(表4)比六盤水市中度侵蝕面積比例(表3)高7.54百分點；從土壤侵蝕量強度來看,煤礦區(qū)中度侵蝕造成的土壤侵蝕量強度比例(表5)比六盤水市中度侵蝕量強度比例(表3)高15.45百分點,說明煤礦區(qū)與土壤侵蝕空間相關性較大。從煤礦分區(qū)來看,南部煤礦區(qū)土壤侵蝕較嚴重,中度及其以上土壤侵蝕量強度達50.33%,遠大于北部和東部煤礦區(qū),主要原因為盤縣是六盤水市煤炭產(chǎn)業(yè)發(fā)展的最重要地區(qū),大型煤礦有80%位于盤縣,煤炭開采造成的崩塌、滑坡和泥石流等大型地質(zhì)災害,使得煤礦區(qū)植被覆蓋和人為的土壤保持等抑制侵蝕因素作用變小。

表4煤礦區(qū)分區(qū)土壤侵蝕面積比例

Table 4Proportion of soil erosion in area in different coal mining areas

煤礦區(qū)不同等級侵蝕面積比例/%微度輕度中度強烈極強烈劇烈北部32.0958.908.390.560.050.00南部19.1958.4219.392.780.220.00東部40.4952.696.500.320.010.00全部26.0950.2521.831.700.140.00

表5煤礦區(qū)分區(qū)土壤侵蝕量強度比例

Table 5Percentages of soil erosions different in intensity in different coal mining areas

煤礦區(qū)不同等級侵蝕量強度比例/%微度輕度中度強烈極強烈劇烈北部6.9964.7324.892.980.410.00南部2.9946.6839.279.841.210.01東部9.0466.8622.081.966.510.00全部4.5940.7846.507.230.890.00

3.3.2典型煤礦區(qū)土壤侵蝕

六盤水市國有煤礦區(qū)土壤侵蝕狀況(圖2～3)顯示,六盤水市大型國有企業(yè)煤礦區(qū)共8個,其中發(fā)耳煤礦和玉舍煤礦位于水城縣,其他6個煤礦均位于盤縣,由表6～7可知,大型國有煤礦區(qū)主要為輕度和中度侵蝕,侵蝕模數(shù)大于5 000 t·km-2·a-1侵蝕面積比例和侵蝕量強度比例均較小,無劇烈侵蝕。

表6國有煤礦區(qū)土壤侵蝕面積比例

Table 6Proportion of soil erosion in area in the state-owned coal mining areas

侵蝕等級侵蝕模數(shù)/(t·km-2·a-1)侵蝕面積比例/%發(fā)耳煤礦玉舍煤礦火鋪煤礦土城煤礦松河煤礦響水煤礦月亮田煤礦金佳煤礦微度<50012.7667.7812.4114.5218.6422.7417.7838.09輕度500～<250067.0930.1568.0958.9558.3162.2058.8447.19中度2500～<500018.241.7218.4922.3419.6914.0019.6912.44強烈5000～<80001.770.281.013.783.061.043.372.21極強烈8000～150000.140.080.010.420.300.020.320.07

表7國有煤礦區(qū)土壤侵蝕量強度比例

Table 7Percentages of soil erosions different in intensity in the state-owned mining areas

侵蝕等級侵蝕模數(shù)/(t·km-2·a-1)侵蝕量強度比例/%發(fā)耳煤礦玉舍煤礦火鋪煤礦土城煤礦松河煤礦響水煤礦月亮田煤礦金佳煤礦微度<5002.7320.732.411.912.534.603.216.68輕度500～<250052.9062.6857.2143.7346.2957.2444.5947.02中度2500～<500037.2511.6836.8040.4438.9333.7138.7935.01強烈5000～<80006.323.533.5511.9110.694.3111.7810.75極強烈8000～150000.811.390.032.011.560.141.630.53

在8個私營煤礦區(qū)中,保興和攀枝花煤礦區(qū)位于水城縣,其他6個私營煤礦均位于盤縣,由表8～9可知,私營煤礦區(qū)侵蝕面積以輕度、中度為主,侵蝕量強度以輕度、中度和強烈為主,強烈侵蝕較明顯,有劇烈侵蝕出現(xiàn)。

綜上所述,私營煤礦區(qū)土壤侵蝕程度要比國有企業(yè)煤礦區(qū)嚴重,主要表現(xiàn)為強烈、劇烈等級侵蝕面積占比增大,就強烈及其以上等級的平均土壤侵蝕量強度占整個研究區(qū)域侵蝕量總強度而言,私營煤礦區(qū)達30%以上,國有煤礦區(qū)為9%,除個別自然因素以外,造成這種現(xiàn)象的原因是相對于私營煤礦國有企業(yè)煤礦的開采保護措施力度更大,開采后復墾工作更加到位。

表8私營煤礦區(qū)土壤侵蝕面積比例

Table 8Proportion of soil erosion in area in the private-owned mining areas

侵蝕等級侵蝕模數(shù)/(t·km-2·a-1)侵蝕面積比例/%祥興煤礦紅果煤礦攀枝花煤礦保興煤礦梓木嘎煤礦上紙廠煤礦森林煤礦昌興煤礦微度<50010.2229.434.968.253.6822.316.4110.22輕度500～<250045.1049.8961.7862.9249.0446.3544.9856.86中度2500～<500031.2513.1728.3923.6833.3718.3335.8925.64強烈5000～<800010.595.264.434.3010.709.2310.714.73極強烈8000～<150002.742.250.440.833.213.722.022.55劇烈≥150000.090.000.000.020.000.070.000.00

表9私營煤礦區(qū)土壤侵蝕量強度比例

Table 9Percentages of soil erosions different in intensity in the private-owned coal mining areas

侵蝕等級侵蝕模數(shù)/(t·km-2·a-1)侵蝕量強度比例/%祥興煤礦紅果煤礦攀枝花煤礦保興煤礦梓木嘎煤礦上紙廠煤礦森林煤礦昌興煤礦微度<5000.943.760.781.360.452.110.711.45輕度500～<250024.0736.8942.0543.9726.2227.0322.9435.94中度2500～<500040.8426.5343.2438.3939.8329.1645.8237.73強烈5000～<800023.8320.2111.9612.4023.1125.2223.3013.43極強烈8000～<150009.7712.621.963.7510.3916.047.2311.46劇烈≥150000.550.000.000.130.000.440.000.00

4討論

由于RUSLE模型是經(jīng)驗統(tǒng)計模型,難以保持很高的精度,因此在運用該模型進行土壤侵蝕估算時模型因子計算的本地化尤為重要[23]。由于條件限制,該研究在這方面做得還欠缺。此外,模型中的地表植被覆蓋因子和水土保持措施因子目前尚無統(tǒng)一的計算方法,在一定程度上影響了模型精度。

六盤水市土壤侵蝕主要發(fā)生在人類負向干擾活動強烈的林地、草地、旱地、煤礦區(qū)和裸地,草地和裸地是最易發(fā)生強烈以上侵蝕的土地利用類型,目前已成為六盤水市可持續(xù)發(fā)展的制約因素。因此,應該采取有效的保護和管理措施來防止生態(tài)環(huán)境惡化。建議加強對六盤水市草地和裸地的防護和治理工作,尤其要防止草地、裸地的石漠化演變。林地多處于坡度大的地區(qū),遇到暴雨沖刷時,也會引起輕度和中度水土流失,要加強對高山區(qū)域的生態(tài)保護工作,護坡護林護草,有效保護好高山地區(qū)草地和林地。旱地和煤礦區(qū)屬于人類影響作用強烈的地區(qū),農(nóng)牧區(qū)和工礦活動區(qū)人類不合理的生產(chǎn)方式是導致土壤侵蝕發(fā)生的主要原因,因此要重視農(nóng)牧區(qū)和工礦區(qū)的水土保持工作,如采取等高耕作、平衡施肥和邊開采邊治理等措施。

針對煤礦區(qū)土壤侵蝕問題要重點加強對礦區(qū)土地復墾與生態(tài)重建研究,包括:土地復墾與生態(tài)重建的基礎理論問題;礦區(qū)土地復墾與生態(tài)重建實用技術;煤礦區(qū)域土地與生態(tài)環(huán)境管理的綜合研究;新技術的應用研究[24]。

所選用的衛(wèi)星影像分辨率為30 m×30 m,而研究區(qū)域面積較小,在處理精度上不夠精確,又因受數(shù)據(jù)的限制,主要采用小比例尺的土壤質(zhì)地數(shù)據(jù),因此后續(xù)研究中要提高數(shù)據(jù)源的精確性,且可選擇多時相數(shù)據(jù)來反映研究區(qū)域土壤侵蝕情況的時空演變特征。

5結論

六盤水市土壤侵蝕以微度、輕度侵蝕為主,侵蝕嚴重的地區(qū)位于水城縣東部和南部,六枝特區(qū)東部以及盤縣中部。六盤水市土壤侵蝕主要發(fā)生在林地、草地、旱地、煤礦區(qū)和裸地,其中草地和裸地是最易發(fā)生強烈以上侵蝕的土地利用類型。六盤水市煤礦區(qū)土壤侵蝕以微度、輕度和中度為主,中度侵蝕造成的水土流失量最多,南部煤礦區(qū)土壤侵蝕比北部和東部嚴重,要重點加強南部煤礦區(qū)的礦區(qū)修復工作。煤礦開采時要做到邊開采邊治理,開采后要進行環(huán)境評估來科學地進行土地復墾工作。六盤水市私營煤礦區(qū)土壤侵蝕程度比國有煤礦區(qū)嚴重,要通過政策約束和政策激勵,使私營煤礦區(qū)在開采和治理過程中注重水土保持。

參考文獻：

[1]OLDEMAN L R.The Global Extent of Soil Degradation[R].Netherlands:[s. n.],1991.

[2]鄂竟平.中國水土流失與生態(tài)安全綜合科學考察總結報告[J].中國水土保持,2008(12):3-7.

[3]蘭敏.基于GIS的秦巴山區(qū)土壤侵蝕空間特征分析：以陜西省寧強縣為例[J].干旱地區(qū)農(nóng)業(yè)研究,2012,30(2):196-200.

[4]朱麗,秦富倉.露天煤礦開采項目水土流失量預測：以內(nèi)蒙古錫林郭勒盟勝利礦區(qū)一號露天煤礦為例[J].水土保持通報,2008,28(4):111-115.

[5]劉寶元,謝云,張科利,等.土壤侵蝕預報模型[M].北京:中國科學技術出版社,2001:217-223.

[6]肖洋,歐陽志云,徐衛(wèi)華,等.基于GIS重慶土壤侵蝕及土壤保持分析[J].生態(tài)學報,2015,35(21):7130-7138.

[7]田鵬,趙廣舉,穆興民,等.基于改進RUSLE模型的皇甫川流域土壤侵蝕產(chǎn)沙模擬研究[J].資源科學,2015,37(4):832-840.

[8]龍?zhí)煊?喬敦,安強,等.基于GIS和RULSE的三峽庫區(qū)土壤侵蝕量估算分析[J].灌溉排水學報,2012,31(2):33-37.

[9]王堯,蔡運龍,潘懋.貴州省烏江流域土壤侵蝕模擬：基于GIS、RUSLE和ANN技術的研究[J].中國地質(zhì),2014,41(5):1735-1747.

[10]RENARD K G,FOSTER G R,WEESIES G A,etal.Predicting Soil Erosion by Water:A Guide to Conservation Planning With the Revised Universal Soil Loss Equation (RUSLE)[K]∥USDA Agricultural Handbook.[s. l.]:[s. n.],1997.

[11]WISCHMEIER W H,SMITH D D.Predicting Rainfal Erotion Loss Cropland East of the Rocky Mountains[K]∥USDA Agricultural Handbook.[s. l.]:[s. n.],1965.

[12]WILLIAMS J R.The Erosion-Productivity Impact Calculator(EPIC) Model:A Case History[J].Royal Society of London Philosophical Transactions,1990,329(1255):421-428.

[13]張科利,彭文英,楊紅麗.中國土壤可蝕性值及其估算[J].土壤學報,2007,44(1):7-13.

[14]鄧輝.基于遙感和GIS的瀘定縣生態(tài)地質(zhì)環(huán)境質(zhì)量評價[D].成都：成都理工大學,2011.

[15]楊潔,郭曉敏,宋月君,等.江西紅壤坡地柑橘園生態(tài)水文特征及水土保持效益[J].應用生態(tài)學報,2012,23(2):468-474.

[16]GUTMAN G,IGNATOV A.The Derivation of the Green Vegetation Fraction From NOAA/AVHRR Data for Use in Numerical Weather Prediction Models[J].International Journal of Remote Sensing,1998,19(8):1533-1543.

[17]鄧輝,何政偉,陳曄,等.基于GIS和RUSLE模型的山地環(huán)境水土流失空間特征定量分析:以四川瀘定縣為例[J].地球與環(huán)境,2013,41(6):669-679.

[18]陸建忠,陳曉玲,李輝,等.基于GIS/RS和USLE鄱陽湖流域土壤侵蝕變化[J].農(nóng)業(yè)工程學報,2011,27(2):337-344.

[19]高峰,華璀,盧遠,等.基于GIS和USLE的欽江流域土壤侵蝕評估[J].水土保持研究,2014,21(1):18-22,28.

[20]李奎,岳大鵬,劉鵬,等.基于GIS與RUSLE的榆林市土壤侵蝕空間分布研究[J].水土保持通報,2014,34(6):172-178.

[21]SL 190—2007,土壤侵蝕分類分級標準[S].

[22]倪晉仁,李英奎.基于土地利用結構變化的水土流失動態(tài)評估[J].地理學報,2001,56(5):611-621.

[23]陳貴廷.基于“3S”技術的呼倫貝爾土壤侵蝕研究[D].蘭州：甘肅農(nóng)業(yè)大學,2009.

[24]卞正富.我國煤礦區(qū)土地復墾與生態(tài)重建研究[J].資源·產(chǎn)業(yè),2005,7(2):18-24.

(責任編輯： 李祥敏)

Evaluation of Soil Erosion of Liupanshui City Based on Revised Universal Soil Loss Equation(RUSLE).

YIN Lu， YAN Qing-wu， BIAN Zheng-fu

(School of Environment Science and Spatial Informatics, China University of Mining and Technology, Xuzhou 221116, China)

Abstract:Mining mineral resources will easily lead to destruction of land resources and deterioration of eco-environment in mining areas. In order to explore effect of coal mining on soil erosion and evaluate the situation of soil erosion in Liupanshui, a major coal mining based city in Southwest Guizhou Province, based on related soil texture, terrain, meteorology, land cover and land use data, status quo and spatial distribution of soil erosion in the city was studied and characterized, with emphasis on distribution of soil erosion in the mining area using the RUSLE model and GIS spatial analysis method. Results show that the soil erosion in Liupanshui is mostly mild and moderate, with severe erosion distributed only in the east and south of Shuicheng County, the east of Liuzhi Special Zone and the central of Panxian County; that soil erosion mainly occurs in woodlands, meadows, dry lands, mining area and bare land where human disturbance is intense; that the mining area in the south of Liupanshui is quite severe in soil erosion as compared with others, but yet dominated with moderate erosion, mainly because mining activities triggers geological disasters which in turn lower the effect of soil erosion controlling factors; and that the situation of private-owned coal minings is much worse than that of state-owned ones. It is, therefore, important for private-owned coal minings to intensify protection and management of the land resources in the process of coal mining.

Key words:GIS;soil erosion;RULSE;coal mining area

收稿日期：2015-11-08

基金項目：國家科技基礎性工作專項(2014FY110800);江蘇高校優(yōu)勢學科建設工程基金(PAPD)

中圖分類號：S157；X87

文獻標志碼：A

文章編號：1673-4831(2016)03-0389-08

DOI：10.11934/j.issn.1673-4831.2016.03.009

作者簡介：尹璐(1990—),女,山東濟寧人,碩士生,主要研究方向為礦區(qū)土地退化空間分析。E-mail: yinlu112233@163.com.

① 通信作者E-mail: 765427936@qq.com