劉 歡,朱谷昌,,劉海利,張建國,歐陽華平

(1.中南大學地球科學與信物物理學院,長沙410083;2有色金屬礦產(chǎn)地質(zhì)調(diào)查中心,北京100012)

0 引言

近年來,隨著全球災害性氣候的增多和人類工程活動的加強,崩塌、滑坡和泥石流等地質(zhì)災害活動也趨向頻繁,嚴重影響了國民經(jīng)濟建設(shè)、社會的穩(wěn)定及人民生命財產(chǎn)安全。地質(zhì)災害危險性評價是隨著地質(zhì)災害的日趨嚴重和相關(guān)學科理論與技術(shù)的發(fā)展而興起的,特別是引進了遙感(RS)和地理信息系統(tǒng)(GIS)方法之后,縮短了危險性評價數(shù)據(jù)獲取的周期,提高了數(shù)據(jù)的精度,便于各類信息的集成分析,同時也促進了該領(lǐng)域的迅速發(fā)展[1]。

在昌都縣進行了地質(zhì)災害詳查、收集地質(zhì)和社會經(jīng)濟資料的基礎(chǔ)上,通過RS地質(zhì)災害解譯、統(tǒng)計分析方法確定危險性評價7個主要影響因素,采用層次分析法與 GIS相結(jié)合的方法進行了地質(zhì)災害的危險性評價,并利用 GIS進行快速分析評價及制圖,有很好的借鑒意義。地質(zhì)災害的危險性是指地質(zhì)災害發(fā)生活動的可能性[2],利用RS和 GIS開展昌都縣地質(zhì)災害調(diào)查和評價應(yīng)用研究,對于制定區(qū)域地質(zhì)災害防治規(guī)劃,部署與實施地質(zhì)災害勘察與防治工作具有重要意義。

1 地質(zhì)環(huán)境條件

昌都縣位于橫斷山脈西北部、青藏高原東南部的邊緣地帶,北與青海省接壤,西與類烏齊縣為鄰,南靠八宿縣、察雅縣,東與江達、貢覺縣毗鄰。地理坐標介于東經(jīng) 96°42′—97°58′,北緯 30°44′—32°19′之間,總面積10 652 km2。昌都縣地勢北高南低,東西呈“W”形。高原寒溫帶季風性氣候混雜,夏季溫和多雨,冬季日照充足,年均氣溫3～8℃,多年平均降水量477.7 mm。

昌都縣位于喀喇昆侖—心嶺—昌都板片和羌塘—唐古拉—保山板片,經(jīng)歷了加里東期—喜山期多期構(gòu)造運動的影響,但巖漿活動和變質(zhì)作用微弱。使大片中生代地層發(fā)生構(gòu)造變形的是燕山運動,形成了區(qū)內(nèi)5條較大的斷裂展布和左行雁列的復式背、向斜。出露的地層主要有前泥盆系、石炭系、二疊系、三疊系、侏羅系、第四系等;出露的巖漿巖有安 山巖、安山質(zhì)凝灰?guī)r、英安巖、斑巖等;變質(zhì)巖主要有千枚巖、片巖等。

2 地質(zhì)災害遙感解譯特征分析

根據(jù)對研究區(qū)已有資料的分析與總結(jié),主要地質(zhì)災害類型有滑坡、崩塌和泥石流。根據(jù)對這3種地質(zhì)災害的ETM+影像特征分析,建立了解譯標志(表1)。在重點區(qū)中,采用高分辨率 IKONOS數(shù)據(jù)進行真彩色合成,同時建立解譯標志[3-5](圖1)。根據(jù)解譯標志,采用人機交互方式對研究區(qū)的3種地質(zhì)災害進行了遙感解譯工作。

鑒于遙感信息的多解性,開展了適量的遙感野外驗證工作,對室內(nèi)遙感解譯建立的3種地質(zhì)災害解譯標志進行修改和完善,驗證室內(nèi)解譯的效果(圖2,圖 3)。

表1 地質(zhì)災害遙感解譯標志Table 1 Remote sensing interpretation marks of geological hazards

圖1 遙感解譯標志Fig.1 Remote sensing interpretation marks

圖2 崩塌點的遙感解譯圖像和實地驗證照片F(xiàn)ig.2 The remote sensing interpretation of slumping and the on-the-spot check photo

圖3 研究區(qū)遙感調(diào)查路線圖Fig.3 Sinvestgation line for remote sensing interpretation in the study area

遙感信息提取和統(tǒng)計中應(yīng)用MAPGIS軟件,以實現(xiàn)屬性數(shù)據(jù)的管理和成果圖件的制作。具體實施步驟:①在PCI軟件中進行遙感影像處理[6],并轉(zhuǎn)為Geotiff格式;②利用 MAPGIS中“圖像分析”工具將Geotiff轉(zhuǎn)換為MAPGIS默認的圖像文件(MSI)格式;③在MAPGIS軟件中建立地質(zhì)災害的屬性信息字段,進行地質(zhì)災害的初步提取[7-9](表2,圖4)。

3 地質(zhì)災害危險性評價指標

地質(zhì)災害的危險性是指存在地質(zhì)災害的風險區(qū)發(fā)生不同程度地質(zhì)災害活動的可能性,危險性評價即為通過對歷史地質(zhì)災害活動程度及可能引起地質(zhì)災害活動的各種因素進行綜合的分析研究,確定地質(zhì)災害可能的危險區(qū)位置范圍、概率等。地質(zhì)災害的危險性指標有很多,根據(jù)研究區(qū)的特點,本文分別選取了7個因素作為危險性評價因子,其權(quán)重結(jié)合前人的經(jīng)驗與該區(qū)的實際情況利用A HP模型計算而確定(表3),具體影像如下:

圖4 研究區(qū)解譯地質(zhì)災害點分布圖Fig.4 Map showing geological hazard distribution in the study area

表2 地質(zhì)災害分類統(tǒng)計Table 2 Statistics of geological hazards

表3 昌都縣地質(zhì)災害評價要素Table 3 Evaluation elements for the geological hazards

(1)地質(zhì)構(gòu)造。主要是指褶皺和斷層。當?shù)貙?產(chǎn)狀與坡向一致時,容易沿層面發(fā)生滑動形成滑坡,特別是地層傾角小于坡角時,滑坡災害更容易發(fā)生。再者,控制區(qū)域性滑坡發(fā)育、分布最基本的構(gòu)造因素就是構(gòu)造線走向,當河谷走向與構(gòu)造線方向垂直時,岸坡比較穩(wěn)定,較少分布滑坡與崩塌;反之,當河谷走向與構(gòu)造線方向一致時,處于向斜河谷的順向坡是產(chǎn)生滑坡的最敏感部位,滑坡的分布多沿斷裂帶成帶狀分布。

(2)地層巖性?；屡c地層巖性的關(guān)系非常密切,軟弱的地層巖性是產(chǎn)生滑坡的物質(zhì)基礎(chǔ),對滑坡影響最為主要的是巖性。同時,已有滑坡體的巖性特征、下伏巖體的巖性特征等也顯著影響著滑坡的類型及其演化趨勢。

(3)地形地貌特征。坡度與滑坡有一定的關(guān)系,一般來說,隨著坡度的增加,包括重力在內(nèi)的剪切力增大,相應(yīng)的滑坡發(fā)生的概率也會增大。在坡度較大的部位,崩塌作用會為滑坡提供大量的物質(zhì),從而為滑坡創(chuàng)造條件?；路植寂c高程也有一定的關(guān)系,坡向的影響主要表現(xiàn)為山坡的小氣候和水熱比的規(guī)律性差異。陽坡巖體風化破碎,易發(fā)生基巖崩滑和泥石流;陰坡土層厚,易發(fā)生淺層塌滑。

(4)植被。植被覆蓋率與滑坡有一定的關(guān)系,植被覆蓋率越高,發(fā)生滑坡等地質(zhì)災害的幾率就越小。

(5)降雨。在具備滑坡的地形、地質(zhì)條件下,降雨在很大程度上是滑坡發(fā)生的誘發(fā)因素,特別是暴雨和綿綿細雨的危害性更大。

(6)水系。在水系較為發(fā)達的區(qū)域,尤其是河流湖泊附近,發(fā)生泥石流滑坡等地質(zhì)災害的可能性較大。

(7)人類活動。隨著城市規(guī)模的擴大,近年來由于城鎮(zhèn)建設(shè),三峽移民搬遷、交通設(shè)施和通訊設(shè)施的建設(shè)、礦山開采等活動,對山體斜坡的改造力度不斷加強,破壞作用亦日趨嚴重,常造成地質(zhì)災害發(fā)生。

4 地質(zhì)災害危險性評價方法流程

采用層次分析法與 GIS相結(jié)合的方法計算研究區(qū)內(nèi)評價單元的地質(zhì)災害危險性指數(shù)[10-16],從而得到研究區(qū)危險性區(qū)劃圖。危險性指數(shù)計算模型的表達式為:

式中,Wj為第j單元的危險性指數(shù);θi為控制地質(zhì)災害危險程度的i類因素的作用權(quán)重;Qi為控制地質(zhì)災害危險程度的i類因素的評分;n為評價因子的個數(shù)。

選取7種評價因子進行危險性評價(表3),其權(quán)重結(jié)合前人的經(jīng)驗與該區(qū)的實際情況利用AHP模型計算而確定,然后利用MAPGIS空間分析平臺,按1 km×1 km的評價單元進行數(shù)據(jù)處理,利用公式(1)通過地圖代數(shù)運算,在MAPGIS空間分析平臺上按1 km×1 km評價單元的最終得分值,在空間分析中提取各單元的地質(zhì)災害危險性指數(shù),并對其進行離散數(shù)據(jù)網(wǎng)格化后,利用網(wǎng)格化后數(shù)據(jù)得出地質(zhì)災害危險性綜合等值線圖及危險性區(qū)劃圖[9,17](表 4,圖 5,圖 6) 。

表4 地質(zhì)災害危險性的等級劃分Table 4 Order division of the geological hazards

圖5 昌都縣地質(zhì)災害危險性分區(qū)圖Fig.5 Geological hazard risk zone division into

4 地質(zhì)災害解譯結(jié)果分析

根據(jù)上述評價分析,將昌都縣危險性區(qū)劃分為:1個地質(zhì)災害高危險區(qū),6個中危險區(qū),14個低危險區(qū),2個無危險區(qū)。

(1)地質(zhì)災害高危險區(qū)。地質(zhì)災害高危險區(qū)占昌都縣總面積的0.42%,主要分布在昌都縣城一帶,縣城內(nèi)滑坡、崩塌、泥石流較多,人文活動頻繁,同時214國道和317國道都在該區(qū)內(nèi)。區(qū)內(nèi)為高山地貌,海拔在 3 100～4 500 m,相對切割深度達1 000 m以上,地形坡度多大于35°,局部地段近于直立,為瀾滄江窄谷區(qū),多為“V”型;構(gòu)造上受俄洛橋斷裂帶與額哎普斷裂帶所控制;出露地層為侏羅系泥質(zhì)粉砂巖、砂巖、泥巖夾中厚-厚層泥質(zhì)粉砂巖以及第四系松散堆積層。

根據(jù)崩塌、滑坡、泥石流地質(zhì)災害點與研究區(qū)危險性區(qū)劃圖疊加分析,大部分地質(zhì)災害都發(fā)生在本區(qū)。危險較大的主要是滑坡和崩塌,危及城鎮(zhèn)重要設(shè)施、交通干線及居民的安全。據(jù)實地調(diào)查,目前這些危害嚴重的地質(zhì)災害點大部分都已得到治理,但還有少部分正在治理當中。

(2)地質(zhì)災害中危險區(qū)。地質(zhì)災害中危險區(qū)占全縣總面積的3.43%,主要分布在卡若鎮(zhèn)、俄洛鎮(zhèn)、拉多鄉(xiāng)、嘎馬鄉(xiāng)、柴維鄉(xiāng)、約巴鄉(xiāng)等地,分布區(qū)人類活動較多。區(qū)內(nèi)主要為高山、極高山地貌,相對切割深度 1 000～2 000 m,地形坡度 30°～45°,局部地段近直立,河(溝)谷呈“V”型;主要的構(gòu)造有額哎普斷裂帶和妥壩斷裂帶、桑多—吉塘斷裂帶;主要出露侏羅系、三疊系、二疊系、石炭系地層,巖性以泥巖、頁巖、安山巖、灰?guī)r和礫巖為主。

地質(zhì)災害中危險區(qū)的主要地質(zhì)災害亦為滑坡、崩塌,主要危害到317國道、縣內(nèi)交通公路和附近居民的安全,需要對其范圍采取地質(zhì)災害防治措施。

(3)地質(zhì)災害低危險區(qū)。地質(zhì)災害低危險區(qū)占全縣總面積的18.56%,分布在妥壩鄉(xiāng)、日通鄉(xiāng)、如意鄉(xiāng)、柴維鄉(xiāng)、拉多鄉(xiāng)等地,人類活動較少。區(qū)內(nèi)為高山地貌,切割深度在1 000 m左右,地形坡度為30°～40°,局部較陡;構(gòu)造上屬于額哎普斷裂帶與妥壩斷裂帶;出露地層有侏羅系、三疊系、侏羅系,巖性為砂巖、泥巖、粉砂巖、頁巖;植被覆蓋率為中等。

地質(zhì)災害低危險區(qū)是已有地質(zhì)災害較為分散的區(qū)域,災害發(fā)生較少。地質(zhì)災害點主要沿公路沿線分布,主要是危害公路的暢通和附近居民的安全。

(4)無地質(zhì)災害危險區(qū)。無地質(zhì)災害危險區(qū)占全縣總面積的77.59%,分布于瀾滄江、昂曲、扎曲、色曲兩側(cè)高山、極高山地區(qū),基本上無人類活動。主要為高山地貌,海拔一般在4 000 m以上。

5 結(jié)論與討論

(1)充分利用遙感技術(shù)、已有的地質(zhì)調(diào)查資料和實地野外資料結(jié)合,建立了研究區(qū)ETM+衛(wèi)星數(shù)據(jù)的崩塌、滑坡、泥石流地質(zhì)災害解譯標志。本次研究的過程中,主要使用的是 ETM+衛(wèi)星數(shù)據(jù),但在昌都縣城關(guān)鎮(zhèn)地區(qū)選用了高分辨率的IKONOS遙感衛(wèi)星數(shù)據(jù),提高了災害發(fā)生重點區(qū)的解譯精度,也有助于本區(qū)地質(zhì)災害產(chǎn)生因素的分析。

(2)根據(jù)對重點區(qū)城關(guān)鎮(zhèn)一帶的實地調(diào)查,發(fā)現(xiàn)城關(guān)鎮(zhèn)大部分崩塌、滑坡的發(fā)生與人類活動關(guān)系密切,大多數(shù)的崩塌、滑坡都產(chǎn)生于道路、河流沿線附近。

(3)研究過程中,利用 GIS管理和空間分析功能,統(tǒng)一了昌都縣地質(zhì)-地形等基礎(chǔ)資料、已有地質(zhì)災害點調(diào)查資料和利用遙感技術(shù)提取的各類信息資料的投影參數(shù)和投影類型,建立了地質(zhì)災害資料信息庫,為研究區(qū)地質(zhì)災害監(jiān)測分析和預防規(guī)劃提供了科學依據(jù)。

另外,由于事物是在不斷變化的,評價區(qū)劃的高危險區(qū)與低危險區(qū)隨著時間的推移及其他條件的變化可以相互轉(zhuǎn)化。高危險區(qū)如果政府防災意識增強、興建抗災建筑、居民點的合理布局等措施有可能降為中危險區(qū)或低危險區(qū);同樣,在低危險區(qū)內(nèi)如果 實施過度的人文-經(jīng)濟開發(fā)活動,危害的程度也會隨之加大。

[1] 齊信,唐川,鐵永波,等.基于 GIS技術(shù)的汶川地震誘發(fā)地質(zhì)災害危險性評價——以四川省北川縣為例[J].成都理工大學學報(自然科學版),2010,37(2):160-167.

[2] 張春山,何淑軍,辛鵬,等.陜西省寶雞市渭濱區(qū)地質(zhì)災害風險評價[J].地質(zhì)通報,2009,28(8):1053-1063.

[3] Vijith H,Madhu G.Western Ghat’s of Kerala(India)through frequency ratio and GIS[J].Environmental Geology,2008,55(7):1397-1405.

[4] 李遠華,姜琦剛.基于遙感調(diào)查與GIS分析的林芝地區(qū)地質(zhì)災害評價[J].國土資源遙感,2006(2):57-60.

[5] 謝韜,何政偉,黃民奇,等.RS與 GIS技術(shù)在庫區(qū)地質(zhì)災害調(diào)查中的應(yīng)用研究[J].測繪科學,2007,32(1):130-134.

[6] 趙文吉.PCI圖像處理教程[M].北京:中國環(huán)境科學出版社,2007:5-14.

[7] 高改萍,楊建宏.GIS在地質(zhì)災害研究中的應(yīng)用[J].人民長江,2003,34(6):32-33.

[8] 魏從玲.基于GIS的巫山縣滑坡地質(zhì)災害危險性評價(碩士學位論文)[D].重慶:西南大學,2008.

[9] 魏平新,湯連生,張建國,等.基于 GIS的廣東省滑坡災害區(qū)劃研究[J].水文地質(zhì)工程地質(zhì),2005(4):6-9.

[10] Jiang Xiao bo,Cui Peng.Hazards assessment of regional debris flows based on geographic information science[J].Wuhan University Journal of Natural Sciences,2007,12(4):651-656.

[11] Piyoosh Rautela,Ramesh Chandra Lakhera.Landslide risk analysis between Giri and Tons Rivers in Himachal Himalaya(India)[J].The International Institute for Aerial Survey and Earth Sciences,2000,2(3):153-160.

[12] 高克昌,崔鵬,趙純勇,等.基于地理信息系統(tǒng)和信息量模型的滑坡危險性評價——以重慶萬州為例[J].巖石力學與工程學報,2006,32(5):991-996.

[13] 謝韜,何政偉,黃民奇,等.RS與 GIS技術(shù)在庫區(qū)地質(zhì)災害調(diào)查中的應(yīng)用研究[J].測繪科學,2007,32(1):130-134.

[14] 徐小飛,馬東濤,孟國才,等.滇藏公路類烏齊—俄洛橋段地質(zhì)災害危險性評估及分區(qū)[J].中國地質(zhì)災害與防治學報,2007(9):20-23.

[15] 張春山,吳滿路,張業(yè)成.地質(zhì)災害評價方法及展望[J].自然災害學報,2003,12(1):96-102.

[16] 張春山,張亞成,張立海,等.中國崩塌、滑坡、泥石流災害危險性評價[J].地質(zhì)力學學報,2004,10(1):27-32.

[17] 張梁,張業(yè)成,胡景江,等.中國地質(zhì)災害時空分布特征與形成條件[J].第四紀研究,2000,20(6):559-566.