基于NDVI的汶川大地震前后北川縣次生地質(zhì)災害區(qū)植被破壞評估

梁超，趙廷寧?，史常青，高志海，周進，張曉龍，薛天貴，肖輝全

(1.北京林業(yè)大學，水土保持與荒漠化防治教育部重點實驗室，100083，北京;2.中國林業(yè)科學研究院資源信息研究所，100091，北京;3.四川省北川縣林業(yè)局，622750，四川北川)

基于NDVI的汶川大地震前后北川縣次生地質(zhì)災害區(qū)植被破壞評估

梁超1，趙廷寧1?，史常青1，高志海2，周進3，張曉龍3，薛天貴3，肖輝全3

(1.北京林業(yè)大學，水土保持與荒漠化防治教育部重點實驗室，100083，北京;2.中國林業(yè)科學研究院資源信息研究所，100091，北京;3.四川省北川縣林業(yè)局，622750，四川北川)

摘要：選取北川縣“5·12”大地震前后2個時相的遙感數(shù)據(jù)，利用歸一化植被指數(shù)(NDVI)提取2期影像的植被信息，以此為基礎(chǔ)，反演植被覆蓋度，采用破壞指數(shù)DDI(the damage degree index)表示3種災害類型區(qū)域震后植被破壞情況。結(jié)果表明：1)研究區(qū)內(nèi)共解譯滑坡103處，崩塌122處，泥石流10處，災害面積共計17.5 km2，震前85%以上區(qū)域處于中植被覆蓋度以上級別，震后中植被覆蓋度以上級別土地面積減少8.01 km2，占災害區(qū)域總面積的45.77%;2)滑坡區(qū)域植被破壞程度相對較低，中度及重度破壞的面積占總災害面積的68.66%，崩塌區(qū)域植被破壞程度次之，中度及重度破壞的面積占總災害面積的88.15%，泥石流區(qū)域的植被破壞最徹底，中度及重度破壞的面積占總災害面積的99.74%;3)植被破壞與海拔、坡度有一定關(guān)系，破壞較嚴重的地區(qū)主要集中在海拔611～1 543m、坡度25°～45°范圍內(nèi)，植被破壞主要集中在重度破壞這一級別，并隨海拔、坡度增加而增加，植被破壞與坡向關(guān)系不大。

關(guān)鍵詞：地質(zhì)災害;歸一化植被指數(shù)(NDVI);植被破壞指數(shù);北川縣

2008年5月12日14時28分，中國四川省發(fā)生了里氏8.0級的汶川大地震，地震在汶川、北川等縣誘發(fā)了大量滑坡、崩塌、泥石流、塌陷以及堰塞湖等次生地質(zhì)災害，造成了嚴重的生態(tài)環(huán)境破壞。北川縣在地震前就是地質(zhì)災害易發(fā)區(qū)，“5·12”地震在北川縣更是誘發(fā)了大規(guī)模的次生地質(zhì)災害，震后遙感解譯滑坡與崩塌就達284處［1］，滑坡和崩塌產(chǎn)生的堆積物在強降雨作用下又進一步誘發(fā)了泥石流［2］，導致大量耕地、林地被掩埋毀壞，自然植被也遭到不同程度的破壞，山地生態(tài)環(huán)境極度退化［3］。

目前對北川縣地震災害的研究主要集中在地質(zhì)災害體數(shù)量、分布、形成條件及誘發(fā)因素、堰塞湖動態(tài)監(jiān)測等幾個方面［4-9］，對地震災害導致的植被破壞情況的研究主要集中在汶川縣［10］、都江堰市［11-12］等地，針對北川縣的研究較少［3］。為了準確評估北川縣災情，為災后植被恢復提供科學依據(jù)，需要了解植被變化情況。筆者利用北川縣地震前后的高分辨率遙感數(shù)據(jù)，通過計算NDVI，獲取地震前后次生地質(zhì)災害發(fā)生區(qū)域植被覆蓋變化等數(shù)據(jù)，為植被重建規(guī)劃提供參考依據(jù)。

1 研究區(qū)概況

北川羌族自治縣位于四川盆地西北邊緣山區(qū)。地理位置介于 E103°44'～104°42'和 N31°14'～ 32°14'之間。全縣東西長91.8 km，南北寬73.8 km，縣域總面積3 084 km2。

選擇北川縣舊縣城至擂鼓鎮(zhèn)地區(qū)作為研究區(qū)，該區(qū)域位于汶川地震高烈度區(qū)的北川縣湔江河谷地帶(圖1)，是典型的地震引起的次生地質(zhì)災害發(fā)育地段，面積90 km2。研究區(qū)處于北川—映秀斷裂帶上，該斷裂帶是“5·12”汶川地震的發(fā)震斷裂帶，是龍門山斷裂的中央斷裂帶，在北川縣境內(nèi)沿擂鼓—舊縣城—鄧家—陳家壩延伸。研究區(qū)出露的地層巖性主要為志留系茂縣群的板巖、千枚巖以及寒武系清平組的泥灰?guī)r、砂板巖［13］。該類巖石風化強烈，節(jié)理裂隙發(fā)育，大部分以殘坡積碎塊石土出露在大于25°斜坡上，為滑坡、崩塌、泥石流等次生地質(zhì)災害的形成提供了物質(zhì)基礎(chǔ)［14］。

研究區(qū)屬于山間河谷地貌，山地海拔一般在700～1 500m。由于河流兩岸侵蝕下切，岸坡多呈25°～45°的斜坡。北川縣全年平均降水量1 399.1mm，年最大降水量2 175.3mm(1967年)，日最大降雨量323.4mm(1977年)，降雨主要集中在5—9月，地表侵蝕作用強烈 ，水土流失嚴重［15］。

圖1 研究區(qū)位置圖Fig．1 Location of the study area

2 研究方法

2.1 數(shù)據(jù)來源

研究使用的數(shù)據(jù)主要為北川縣高分辨率遙感數(shù)據(jù)和該縣1∶5萬地形圖?！?·12”地震前遙感影像為SPOT-5，該影像包括4個多光譜波段(分辨率為10m)和1個全色波段(分辨率為2.5m)，覆蓋范圍為北川縣全境;震后遙感影像為Formosat-2，該影像也包含4個多光譜波段(分辨率為8m)和1個全色波段(分辨率為2m)，覆蓋范圍為北川境內(nèi)斷裂帶及周邊地震災害嚴重地區(qū)。

2.2 遙感影像預處理

1)幾何精校正及投影轉(zhuǎn)換。在ERDAS9.2中，利用北川縣震前1∶5萬地形圖，均勻選取10個控制點(GCP)，采用三次多項式對震前2007年夏季SPOT-5影像進行幾何精校正。以校正后的SPOT-5影像為標準，同樣選取10個控制點(GCP)，采用三次多項式對震后2008年夏季Formosat-2影像進行校正［16］，校正的同時對 SPOT-5及 Formosat-2影像作投影轉(zhuǎn)換，根據(jù)北川縣經(jīng)緯度范圍將影像統(tǒng)一到WGS84 坐標系，UTM 投影［17］。由于“5·12”地震造成大量房屋倒塌、道路毀壞以及次生地質(zhì)災害導致的堰塞湖，給校正帶來了一定困難，在選取控制點時，盡量在震害較輕的區(qū)域選點。

2)影像拼接及剪裁。利用ERDAS中的Mosaic工具將地震前、后各一景影像中涉及到北川縣的SPOT-5及Formosat-2影像分別拼接成震前、震后2幅影像，并利用北川縣邊界矢量圖分別對影像進行剪裁。

3)影像融合。采用主成分融合方法，分別對SPOT-5和Formosat-2影像進行全色影像和多光譜影像的融合處理，融合后的SPOT-5影像分辨率2.5m，F(xiàn)ormosat-2影像分辨率2m。融合后的影像既能保留多光譜影像的高光譜分辨率，還能保留全色影像的高空間分辨率［18］。

2.3 遙感影像解譯

在對研究區(qū)地震次生災害野外調(diào)查的基礎(chǔ)上，結(jié)合其他相關(guān)解譯研究成果［19］，利用ArcGIS對震后Formosat-2影像進行次生地質(zhì)災害解譯。災害類型分為滑坡、崩塌、泥石流3種，采用目視解譯的方法，勾繪出3種災害類型的分布信息。采用波段1(R)、2(G)、3(B)不同的波段組合方式，得到最接近真實色彩的Formosat-2影像。

2.4 植被破壞計算

根據(jù)解譯結(jié)果，采用差分植被指數(shù)NDVI分別提取3種災害類型范圍內(nèi)及震前相同區(qū)域的植被指數(shù)，計算公式為

式中：INDVI為差分植被指數(shù)NDVI值;B2為影像第2波段;B3為影像第3波段。

使用植被覆蓋度fc反映研究區(qū)內(nèi)植被變化情況［20-21］，計算公式為

式中：fc為植被覆蓋度;INDVIS為研究區(qū)無植被覆蓋區(qū)域像元NDVI值;INDVIV為研究區(qū)純植被覆蓋區(qū)域像元NDVI值。

依據(jù)《開發(fā)建設項目水土保持技術(shù)規(guī)范》［22］中植被覆蓋分級要求，將地震前后不同災害類型區(qū)域植被覆蓋度分為5個等級：低(fc＜0.1)、中低(0.1≤fc＜0.3)、中(0.3≤fc＜0.5)、中高(0.5≤fc＜0.7)、高(fc≥0.7)。

采用破壞指數(shù)DDI(the damage degree index)即研究區(qū)地震前后不同災害類型區(qū)域NDVI的差值來表示3種災害類型區(qū)域震后植被破壞情況［12］，計算公式為

式中：IDDI為植被破壞指數(shù)值;INDVI1為震前NDVI值;INDVI2為震后NDVI值。

將植被破壞指數(shù)分為3個等級：輕度破壞(IDDI≤0.05)、中度破壞(0.05 ＜IDDI≤0.15)、重度破壞(IDDI＞0.15)。

3 結(jié)果與分析

3.1 解譯結(jié)果驗證

解譯結(jié)果會存在誤差，主要是因為：1)研究區(qū)內(nèi)存在的采石場、裸露的廢棄地等會對解譯時判斷有一定影響;2)研究區(qū)內(nèi)大部分地質(zhì)災害活動強烈，但有個別潛在災害體在影像上較難識別;3)2期遙感影像的分辨率、清晰度不同會帶來一定的解譯誤差。

為了驗證解譯類型的正確性并分析上述可能存在的誤差，選擇北川縣擂鼓鎮(zhèn)及舊縣城魏家溝流域進行解譯結(jié)果驗證。主要進行現(xiàn)場調(diào)查驗證災害類型和使用GPS測量災害體面積2項工作。共選取10處災害區(qū)域作為驗證數(shù)據(jù)，驗證結(jié)果見表1。對比現(xiàn)場調(diào)查和解譯結(jié)果可以看出，有3處災害區(qū)域面積精度在95%以上，4處災害區(qū)域面積精度在90%～95%之間，2處災害區(qū)域面積精度在80% ～90%之間，只有1處災害區(qū)域面積精度較低。另外，在現(xiàn)場調(diào)查過程中，發(fā)現(xiàn)2處崩塌屬于漏判，主要是因為面積小，在遙感影像上很難識別。

由于研究區(qū)屬于震害嚴重地區(qū)，該區(qū)域誘發(fā)的次生地質(zhì)災害規(guī)模大、數(shù)量多，且絕大多數(shù)災害造成了較強的破壞，在遙感影像上較易識別。同時，通過現(xiàn)場調(diào)查與解譯結(jié)果對比驗證可知，調(diào)查結(jié)果與解譯結(jié)果有很好的擬合性，因此，利用遙感影像可以較為準確地確定不同次生地質(zhì)災害的類型及邊界。

3.2 地震前后植被覆蓋度等級狀況

通過對北川縣地震前后遙感影像的處理，得到地震前后滑坡、崩塌、泥石流3種災害類型區(qū)域植被覆蓋度等級狀況，結(jié)果見表2?？梢钥闯觯喊l(fā)生地質(zhì)災害的區(qū)域震前植被覆蓋度較高，地震后，植被覆蓋度大于30%的土地面積急劇減少，覆蓋度小于30%的土地面積有所增加。其中高覆蓋度等級的土地面積減少最多，達到4.17 km2，低覆蓋度等級的土地面積增加最多，達到了4.22 km2，說明地震引起的地質(zhì)災害對植被造成了嚴重的破壞。

表1 遙感解譯結(jié)果的野外調(diào)查驗證Tab．1 Comparison and verification for the results of RS intepretation and site investigation

發(fā)生滑坡的區(qū)域震前植被覆蓋度在中植被覆蓋度以上等級的面積占86.8%，震后降到了41.2%，中低植被覆蓋度及低植被覆蓋度面積比例由13.2%升到了58.8%;發(fā)生崩塌的區(qū)域震前植被覆蓋度在中植被覆蓋度以上等級的面積占96.4%，震后降到了31.2%，中低植被覆蓋度及低植被覆蓋度面積比例由3.6%升到了68.8%;發(fā)生泥石流的區(qū)域震前植被覆蓋度在中植被覆蓋度以上等級的面積占93.1%，震后降到了72.9%，中低植被覆蓋度及低植被覆蓋度面積比例由6.9%升到了27.1%。

3.3 震后植被破壞情況

統(tǒng)計的研究區(qū)滑坡、崩塌、泥石流3種災害類型發(fā)生區(qū)域植被破壞情況見表3?？梢钥闯觯?種災害類型區(qū)域內(nèi)，植被破壞狀況均以重度破壞的面積最多，其中，滑坡區(qū)域重度破壞的面積略少，僅占滑坡災害面積的55.53%，崩塌和泥石流區(qū)域重度破壞的面積較大，分別占其災害面積的80.20%和92.40%，滑坡和崩塌區(qū)域中，還有部分面積植被屬于輕度破壞，分別占其災害面積的31.35%和12.03%。

其原因主要是由于震前這些區(qū)域植被覆蓋很高，地震誘發(fā)的地質(zhì)災害對區(qū)域內(nèi)的植被造成了不同程度的破壞。在現(xiàn)場實地調(diào)查中也發(fā)現(xiàn)：泥石流破壞力最強，發(fā)生該類型災害時對地表植被的破壞是毀滅性的;崩塌對植被有一定的破壞，如崩落的石塊會掩埋植被、砸斷樹木等;滑坡發(fā)生時，不同的滑動速度會對植被造成不同程度的破壞，滑動速度慢的滑坡其滑坡體上的樹木不會完全遭到破壞，會隨土體滑動而歪斜，在滑坡停止滑動后依然能夠存活，形成馬刀樹(醉林)。因此，發(fā)生泥石流的區(qū)域基本不會有植被存活，發(fā)生崩塌的區(qū)域大部分植被被破壞，可能會有少量植被存活，發(fā)生滑坡的區(qū)域根據(jù)滑坡速度等的不同，會有部分植被存活。

表3 3種災害類型區(qū)域植被破壞狀況Tab．3 Vegetation damage in landslides，collapses andmudslides disaster

3.4 植被破壞與海拔的關(guān)系

在ArcGIS中，通過疊加研究區(qū)植被破壞指數(shù)數(shù)據(jù)和DEM數(shù)據(jù)，確定植被破壞與海拔的關(guān)系，結(jié)果見表4。可以看出，植被破壞主要集中在海拔611～1 543m范圍，破壞面積14.42 km2，占破壞總面積的82.40%。各海拔范圍內(nèi)，重度破壞一等級與其他等級相比占有比例最大，均占55%以上，同時，該破壞等級面積占的比例在海拔611～1543m范圍內(nèi)隨海拔增加而增加，當海拔大于1 543m時，植被破壞面積有所減小。說明植被破壞與海拔有一定關(guān)系，但并非線性關(guān)系。

表4 植被破壞與海拔的關(guān)系Tab．4 Relationship between vegetation damage and elevation

3.5 植被破壞與坡度的關(guān)系

疊加研究區(qū)植被破壞數(shù)據(jù)和坡度數(shù)據(jù)，確定植被破壞與海拔的關(guān)系，結(jié)果見表5?？梢钥闯?，植被破壞與坡度有一定關(guān)系，植被破壞主要集中在坡度25°～45°范圍，破壞面積 11.61 km2，占破壞總面積的66.30%。植被破壞面積隨坡度的增加而增加，在坡度達到35°～45°時破壞面積最大，坡度大于45°時，破壞面積隨坡度增加的趨勢有所降低。各坡度范圍中，重度破壞這一等級與其他等級相比占的比例最大，約60%以上，同時，該破壞等級的面積在坡度25°～45°范圍內(nèi)隨坡度的增加而增加，當坡度大于45°時，植被破壞面積有所減小。說明植被破壞與坡度也并非線性關(guān)系。

3.6 植被破壞與坡向的關(guān)系

疊加研究區(qū)植被破壞數(shù)據(jù)和坡向數(shù)據(jù)，確定植被破壞與坡向的關(guān)系，結(jié)果見表6。可以看出，9個等級的坡向范圍中，各個等級的植被破壞面積及所占的比例大致相當，植被破壞與坡向無明顯關(guān)系。各坡向范圍中，重度破壞等級與其他等級相比均占的比例最大，在30%以上。

表5 植被破壞與坡度的關(guān)系Tab．5 Relationship between vegetation damage and slope

表6 植被破壞與坡向的關(guān)系Tab．6 Relationship between vegetation damage and aspect

4 結(jié)論

1)研究區(qū)內(nèi)共解譯滑坡103處，崩塌122處，泥石流10處，災害面積共計17.5 km2。研究區(qū)震前植被覆蓋較高，85%以上的區(qū)域處于中植被覆蓋度以上級別，地震后中植被覆蓋度以上級別土地面積減少8.01 km2，占災害區(qū)域總面積的45.77%，說明汶川地震對研究區(qū)內(nèi)植被造成了嚴重破壞。

2)3種不同類型災害對植被破壞程度不同，發(fā)生滑坡的區(qū)域植被破壞程度相對較低，中度及重度破壞的面積占總災害面積的68.66%，發(fā)生崩塌的區(qū)域植被破壞程度次之，中度及重度破壞的面積占總災害面積的88.15%，發(fā)生泥石流的區(qū)域植被破壞最嚴重，中度及重度破壞的面積占總災害面積的99.74%。

3)植被破壞與海拔有較為密切的關(guān)系，破壞較嚴重的地區(qū)主要集中在海拔611～1 543m范圍內(nèi)，植被破壞主要集中在重度破壞這一級別，且在海拔611～1 543m范圍內(nèi)隨海拔增加而增加，當海拔大于1 543m時，植被破壞程度相對減輕。

4)植被破壞與坡度也有一定關(guān)系，破壞較嚴重的地區(qū)主要集中在坡度25°～45°范圍內(nèi)，植被破壞主要集中在重度破壞這一級別，且在坡度25°～45°范圍內(nèi)隨坡度增加而增加，當坡度大于45°時，植被破壞程度相對減輕。

5)植被破壞與坡向關(guān)系不大。

5 參考文獻

［1］國家減災委員會-科學技術(shù)部抗震救災專家組.汶川地震災害綜合分析與評估［M］.北京：科學技術(shù)出版社，2008：90-120

［2］朱忠禮，劉希林，尚志海.四川省北川縣趙家溝地震次生泥石流形成條件分析［J］.中國地質(zhì)災害與防治學報，2010，21(4)：57-62

［3］范建容，張建強，田兵偉，等.汶川地震次生災害毀壞耕地的遙感快速評估方法：以北川縣唐家山地區(qū)為例［J］.遙感學報，2008，12(6)：917-924

［4］王猛，王軍，江煜，等.汶川地震災區(qū)北川縣地震地質(zhì)災害遙感調(diào)查與分析［C］.中國遙感應用協(xié)會2010年會暨區(qū)域遙感發(fā)展與產(chǎn)業(yè)搞成論壇論文集，2010：249-255

［5］黃庭，張志，谷延群，等.基于遙感和GIS技術(shù)的北川縣地震次生地質(zhì)災害分布特征［J］.遙感學報，2009，13(1)：177-182

［6］黃潤秋，李為樂.“5·12”汶川大地震觸發(fā)地質(zhì)災害的發(fā)育分布規(guī)律研究［J］.巖石力學與工程學報，2008，27(12)：2585-2592

［7］程思思，甘甫平，劉欽甫，等.CBERS-2B星數(shù)據(jù)在“5·12”地震滑坡調(diào)查中的應用：以唐家山至北川縣城堵江河道滑坡為例［J］.國土資源遙感，2009(1)：93-96

［8］王根龍，張軍慧，劉紅帥.汶川地震北川縣城地質(zhì)災害調(diào)查與初步分析［J］.中國地質(zhì)災害與防治學報，2009，20(3)：47-51

［9］王志一，甘甫平，李賢慶，等.CBERS-2B星數(shù)據(jù)在“5·12”地震災害監(jiān)測中的應用：以唐家山堰塞湖動態(tài)監(jiān)測為例［J］.國土資源遙感，2009(1)：89-92

［10］王芳，楊武年，劉漢湖，等.基于混合像元線性分解的5·12地震后汶川縣由地質(zhì)災害導致的植被覆蓋度變化研究［C］.International Conference on Remote Sensing，2010：310-314

［11］倪忠云，何政偉，趙銀兵，等.汶川地震前后都江堰植被蓋度變化的遙感研究［J］.水土保持研究，2009，16(4)：45-48

［12］Yang Cunjian，Ren Xiaolan，Huang He.The vegetation damage assessment of the Wenchuan earthquake ofmay 2008 using remote sensing and GIS［J］.Nature Hazards，2012，62：45-55

［13］唐川，齊信，丁軍，等.汶川地震高烈度區(qū)暴雨滑坡活動的遙感動態(tài)分析［J］.地球科學：中國地質(zhì)大學學報，2010，35(2)：317-323

［14］孫瑜.基于遙感技術(shù)的5·12震后北川縣山地災害分布規(guī)律研究［D］.北京：中國地質(zhì)大學，2010：13-15

［15］北川縣志編纂委員會.北川縣志［M］.北京：方志出版社，1996：100-107

［16］黨安榮，王曉棟.ERDAS IMAGINE遙感圖像處理方法［M］.北京：清華大學出版社，2004：69-80

［17］吳秀芹，張洪巖.ArcGIS 9地理信息系統(tǒng)應用于實踐［M］.北京：清華大學出版社，2008：138-145

［18］趙英時.遙感應用分析原理與方法［M］.北京：科學出版社，2008：254-272

［19］王平，王殿琦.飛向汶川：5·12地震災害航空遙感調(diào)查［M］.北京：科學出版社，2009：68-111

［20］賈寶全，邱爾發(fā)，張紅旗.基于歸一化植被指數(shù)的西安市域植被變化［J］.林業(yè)科學，2012，48(10)：6-12

［21］劉淼，秦大庸，劉家宏，等.基于NDVI的山西省植被覆蓋度變化研究［J］.人民黃河，2009，31(5)：17-18

［22］GB 50433—2008開發(fā)建設項目水土保持技術(shù)規(guī)范［S］.北京：中國計劃出版社，2008：16-17

Vegetation damage assessment of Beichuan County after the earthquake inmay 2008 Based on NDVI

Liang Chao1，Zhao Tingning1，Shi Changqing1，Gao Zhihai2，Zhou Jin3，Zhang Xiaolong3，Xue Tiangui3，Xiao Huiquan3

(1.Key Lab.of Soil＆ Water Conservation and Desertification Combating of theministry of Education，Beijing Forestry University，100083，Beijing，China;2.Research Institute of Forest Resource Information Techniques，Chinese Academy of Forestry，100091，Beijing，China;3.Beichuan County Forestry Bureau of Sichuan Province，622750，Beichuan，Sichuan，China)

Abstract：This study selected two phases remote sensing images before and after Wenchuan Earthquake onmay 12th，2008 in Beichuan County.The two images before and after the earthquake was obtained to extract vegetation covering information using normalization difference vegetation index(NDVI)，and NDVI was used to estimate vegetation coverage.Then the damage degree index(DDI)of the vegetation was defined，indicating the difference of the normalized difference vegetation index before and after earthquake.The results showed that：1)103 landslides，122 collapses and 10mudslides were interpreted from Formosat-2 image taken after Wenchuan Earthquake，and the area of disaster was 17.5 km2.The area ofmedium vegetation coverage status was 85%before the earthquake，and it reduced by 8.01 km2after the earthquake，which accounted for 45.77%of the whole area of disaster.2)The vegetation damage was relative less serious in landslides regions，wheremoderate and severe damage vegetation occupied 68.66%of the whole area.The vegetation damage wasmodest in collapses regions，wheremoderate and severe damage vegetation occupied 88.15%of the whole area.Most vegetation damage occurred inmudslides regions，wheremoderate and severe damage vegetation occupied 99.74%of the whole area.3)The vegetation damage was closely related to the elevation and slope grade，which was prone to occur at the elevation of 611-1 543m and slope grade of 25°-45°regions.The vegetation damage was worse，most of which were severe damage with a decreasing trend along with elevation and slope grade increasing，but weakly related to the aspect factor.

Key words：geological disasters;normalization difference vegetation index(NDVI);vegetation damage degree index;Beichuan County

TP79

A

1672-3007(2013)04-0086-07

2013-03-11

2013-05-09

國家林業(yè)公益性項目“四川地震災區(qū)災后植被恢復及可持續(xù)發(fā)展關(guān)鍵技術(shù)研究與示范”(201104109)

梁超(1988—)，男，碩士研究生。主要研究方向：工程綠化與水土保持。E-mail：chao.lyc@gmail.com

簡介：趙廷寧(1962—)，男，教授，博士生導師。主要研究方向：工程綠化與水土保持。E-mail：zhtning@bjfu.edu.cn

(責任編輯：宋如華)