楊建鵬,韓 晨,王 龍,王景麗
(1.西北大學地質學系,陜西 西安710069;2.渭南師范學院資源環境系,陜西 渭南714000;3.西南大學地理科學學院,重慶404000)
地理信息系統技術在滑坡及泥石流災害評估中的應用
楊建鵬1,韓 晨2,王 龍3,王景麗1
(1.西北大學地質學系,陜西 西安710069;2.渭南師范學院資源環境系,陜西 渭南714000;3.西南大學地理科學學院,重慶404000)
地理信息系統已廣泛應用于礦產資源勘探、管理,土地資源管理與地表覆被變化,環境保護,城市規劃與建設,政府部門行政與決策等諸多方面。以韓城為例探討地理信息系統在滑坡、泥石流等地質災害中的諸多應用。利用AacGIS圖層疊加及圖層要素綜合計算得出研究區韓城大斷裂以東的黃河河道階地及河谷沖積區為潛在的滑坡、泥石流高發區。
地理信息系統;地質災害;預測評估
目前大部分地區地質災害防治工作由于經濟技術等因素的限制仍處于艱難的被動防御狀態。其中一個重要原因就是缺乏一種以先進技術為依托從而建立起系統化的實用有效的地質災害防治系統。由于地質災害發生的地質水文背景與巖土的質地結構、降雨量、地表植被覆蓋率以及其它外在因素錯綜交錯且在時間上又有較大的不確定性。所以地質災害確有其難以預測的一方面,但地災部門調查手段和管理落后往往造成對地質災害數量、性質、規模、危害程度等心中無數。常常難以兼顧或時常處于消極的被動防御和善后處理。這樣無疑對廣大人民群眾生命財產和國家經濟發展造成不可估量的損失。因此,建立一套基于GIS技術的地質災害監測、評估和防御系統意義重大,這樣各地質災害決策管理專家就可以通過GIS提供的全方位實時的區域空間信息,全面掌握可能發生地質災害的位置、時間、性質、規模、危害程度、數量并通過合適的模型匹配分析從而得到相應的整治措施,防患于未然。
地理信息系統(GIS)是用于輸入、存儲、查詢、分析和顯示地理數據的計算機系統。是介于信息科學、空間科學和地球科學之間的交叉科學,同時也是一門新興的技術學科。它是計算機科學、遙感學、信息技術與現代地學理論和方法的有機結合。此外地理信息系統是基于數據庫系統、地理信息的空間分析、模擬和地圖可視化的計算機系統。像其他任何信息技術一樣,GIS可分為如下四個部分[1]:計算機系統、GIS軟件、智囊和設施。
韓城市位于關中平原東北隅,距省會西安240余公里,北依宜川,西鄰黃龍,南接合陽,東隔黃河與山西省河津、鄉寧、萬榮等縣市相望。全市總面積1 621 km2,耕地42萬畝,地形地貌為“七山一水二分田”[2]。由于本區位于關中平原像陜北黃土高原過度的黃土臺塬地帶,所以黃土地貌諸如:塬、梁、峁及流水沖刷溝壑區等典型的黃土地貌遍布。由于黃土固有的疏松、黏性較小、垂直節理的巖土特性、本區特大降水多集中余夏秋之交及地表植被覆蓋率較小。所以本區發生滑坡、泥石流或山崩的可能性極大。屬地區地質災害重點監測區之一。
構造方面,研究區位于祁(連山)呂(梁山)賀(蘭山)山字形構造的前弧東翼與新華夏構造體系第三沉降帶的復合部位。以東北~西南向的山前大斷裂(即韓城大斷層)為界,東南面屬汾渭地塹沉降區,西北面則屬于鄂爾多斯臺向斜的陜北高原邊緣褶皺區。構造復雜,礦產豐富,存在著發生中強度以上地震的條件。本區地層,總體上可分兩個單元。以東北~西南向的大斷層為界,東南部為該正斷層的上盤,系下沉地層,為一套厚層第四紀沉積地層;西北部為相對上升的侵蝕地層,以古生代巖層為主,上覆黃土或坡積、殘積物[3](如圖1)。
三維地形圖,它是遙感、地理信息系統、三維仿真等高新技術的結合。它是以三維電子地圖數據庫為基礎,按照一定比例對現實世界或其中一部分的一個或多個方面的三維、抽象的描述。三維地形圖根據高程的不同賦予不同顏色值,制作出三維地形圖,用以表達不同的地形起伏情況。地理信息系統中主要采用等高線法、地貌暈渲法(一種旨在模擬太陽光與地表要素相互作用下的地形容貌)、分層設色法(用不同顏色表示不同高度帶已明顯表現出地貌的起伏特征)、透視法(地形的三位視圖)。其受控于四個參數:觀察方位、觀察角度、觀察距離、豎向比例尺。ArcGIS中的3-D Analyst擴展模塊提供了用于觀察參數操作的圖形界面。

圖1 研究區地層巖性空間分布
數字高程模型(DEM)作為地理數據庫中最為重要的空間信息資料和賴以進行地形分析的核心數據系統,不僅可以提取地形起伏度、坡度和坡向等地形參數,同時能夠成為定量分析地形地貌和地質災害評價的基礎。格網數字高程模型,即格網DEM,它是把DEM覆蓋區劃分成規則格網,每個網格大小和形狀都相同,用相應矩陣元素的行列號來實現網格點的二維地理空間定位,第三維為特性值,可以是高程或屬性[4]。

圖2 研究區坡度專題圖
如圖3所示:本文將研究區坡度以十度為間隔分為九個層次,分別以九種不同的漸變色表示??梢悦黠@看出坡度大小多集中在45°~85°之間。而坡向則被分為東、南、西、北、東北、東南、西北、西南、八個方向并分別用八種不同的漸變色表示;研究區坡度以韓城大斷裂為界西部中低山區坡度(降)較小,斷裂以東的黃河河道沖積區由于處于斷層上盤,所以沉降較大,且表面幾乎均被第四紀風積黃土及河流沖積物覆蓋(圖2),屬水土流失嚴重區。
如圖3研究區坡向則大多呈東-東南方向,表現出從西北向東南高度逐漸降低的趨勢,直到黃河河道處降為水平。而區域居民聚居區(韓城市區)及公路正好處于山坡谷底低凹處及主要的順坡向區,這無疑處于滑坡泥石流等地質災害等的高度威脅區。此外,在這種區域坡度較大的情況下極容易產生重力作用下的巖土塌陷及巖土的蠕動或滑坡陷現象,尤其是在遇上季節性暴雨時無疑是滑坡、泥石流等地質災害的高度危險區。

圖3 研究區東-東南坡向(深色區)專題圖
根據陜西省地情數據庫記載研究區滑坡、泥石流等地質災害在一年中不同月份發生的次數百分率統計,如圖可以看出當降雨量達到一定(40 mm)時發生滑坡、泥石流等地質災害的頻率顯著增大,所以降雨的集中出現是影響研究區滑坡、泥石流等地質災害發生的主要影響因素。此外,空間分布上的特點為在東部夏季風暖濕氣流自東向西推進過程中,受西部山區地形的抬升所形成的地形雨影響,研究區內降水量不均衡,其分布趨勢是由西北向東南遞減,西部最大降雨量達 657.0 mm,東南部最少為 468.4 mm[2]。所以降水多集中于西部山地(海拔較高區)加之斷裂帶附近坡度的急劇下降也使主要人口聚居區成為發生滑坡泥石流等地質災害的潛在高發區。
1)地表覆被度的提取
由于植物的光譜特征可使其在遙感影像上有效的與其他地物相區別。同時不同植物各有其自身的光譜特征,從而成為區分研究區地表植被覆蓋類型、長勢等[5]。如利用Landsat TM光學遙感影像解譯提取植被區域、計算歸一化植被指數(NDVI)[6]并進而反演植被覆蓋率來反映地表植被的覆蓋情況。譬如,通過等密度模型可以從遙感影像中反演植被覆蓋度:

其中NDVI為像元的歸一化植被指數,NDVImax、NDVIsoil分別是遙感影像上裸土的最大和裸土的NDVI值。

圖4 降雨量與地質災害率關系(引自林斌)

圖5 研究區植被覆蓋度(深色)遙感圖像
2)地表歸一化植被指數NDVI的提取
歸一化植被指數[6]是從Landsat TM影像獲得。
計算式為:

其中r3、r4分別為 Landsat TM第三(紅光)、第四(近紅外)通道的反射率。

圖6 研究區NDVI(深色)遙感圖像
如圖6(兩條斜線區)可以看出從西部山地到斷裂帶以東植被覆蓋度明顯減小,表明斷裂帶及其以東為表層剝蝕區植被覆蓋度較低。此外,如圖7(斜線區)所示研究區植被覆蓋率(深色)較低且由于斷裂帶的存在,地形地貌復雜、剝蝕作用強烈從而極易發生滑坡、泥石流等地質災害;此外,水土流失也非常嚴重。

圖7 綜合預測地質災害高發區(深色區)
綜上,通過對研究區坡度、坡向、地層巖性、植被覆蓋度、植被歸一化指數及降雨量等因素的綜合分析,利用AacGIS圖層疊加及圖層要素綜合計算得出研究區韓城大斷裂以東的黃河河道階地及河谷沖積區為潛在的滑坡、泥石流高發區(圖7斜線區)。
本文簡要介紹了地理信息系統的主要概念、相關工具;同時以韓城為實例簡要分析了滑坡、泥石流等地質災害發生的潛在的地質構造背景、地形地貌因素、巖土質地、結構因素、地表覆被因素,植被覆被情況,降雨因素等;得出韓城存在發生滑坡、泥石流等地質災害的潛在的高發區域主要位于韓城大斷裂以東的河谷地帶,即第四紀風積黃土和河道沖積物廣泛沉積的黃土臺塬及河谷川道地帶,而這種地理單元正是人口、公路和工業廣泛分布區(如圖8),值得重視。
[1]陳健飛,Kang-tsung Chang等.地理信息系統導論[M].北京:科學出版社.2003,1-3.
[2]陜西省地情資料庫[DB].
[3]陜西省區域地質志[M].北京:地質出版社.1989.520-540.
[4]湯國安,張勇,劉詠梅等.不同比例尺DEM提取地面坡度研究.水土保持通報[J].2001,21(1):53-56.
[5]梅安新,秦其明等遙感導論[M].北京:高等教育出版社.2001,225-247.
[6]陳云浩,李曉兵,陳晉等.1983-1992年中國陸地植被NDVI演變特征的變化向量分析.遙感學報[J].2002,6(1):12-17.
P642.2
B
1004-1184(2012)04-0181-03
2012-05-25
楊建鵬(1987-),男,陜西富平人,在讀碩士研究生,主攻方向:礦物學、巖石學、礦床學。