丁麗，李再興，李學(xué)問，張瑞香

(河南省地質(zhì)礦產(chǎn)勘查開發(fā)局測繪地理信息院，鄭州　450006)

?

基于層次分析法的登封市地質(zhì)災(zāi)害危險性評價

丁麗，李再興，李學(xué)問，張瑞香

(河南省地質(zhì)礦產(chǎn)勘查開發(fā)局測繪地理信息院，鄭州450006)

摘要：文章針對登封市地質(zhì)環(huán)境條件、發(fā)育特征等特性，選取了災(zāi)害點密度、巖土體類型、斷裂發(fā)育程度、地貌類型、坡度、降雨量、人類工程活動等7個地質(zhì)災(zāi)害危險性評價因子，通過層次分析法確定評價因子權(quán)重，并采用專家系統(tǒng)法賦值，建立了該市地質(zhì)災(zāi)害危險性評價模型。通過劃分評價單元及利用危險性評價模型，基于Arcgis空間分析功能將該市地質(zhì)災(zāi)害危險性分為低危險區(qū)、中危險區(qū)和高危險區(qū)。研究成果對相似地區(qū)地質(zhì)災(zāi)害危險性的劃分也提供了一定的借鑒意義。

關(guān)鍵詞：層次分析法；地質(zhì)災(zāi)害；危險性分區(qū)；危險性評價

地質(zhì)災(zāi)害是自然和人為因素導(dǎo)致地質(zhì)環(huán)境或地質(zhì)體發(fā)生變化，也給人類和社會造成危害的災(zāi)害事件，如滑坡、泥石流、崩塌、地面塌陷、地面沉降等，其主要體現(xiàn)在對居民點、交通、水利水電工程、礦山等造成的危害。目前被用于區(qū)域環(huán)境質(zhì)量評價的方法已有多種，在對地質(zhì)災(zāi)害危險性評價的方法的選擇上一直困擾著廣大地質(zhì)工作者。隨著數(shù)理知識在各學(xué)科中的應(yīng)用，相互穿插，特別是GIS與數(shù)學(xué)模型結(jié)合的評價方法，使得地質(zhì)災(zāi)害危險性評價有了得心應(yīng)手的方法，當(dāng)影響因素多、分布復(fù)雜時，層次分析理論使地質(zhì)災(zāi)害理論評價更加清晰化、層次化、定量化，且較為成熟，在大區(qū)域以及其他領(lǐng)域中的應(yīng)用也尤為方便、簡捷[1]。

文中結(jié)合現(xiàn)有地質(zhì)災(zāi)害信息，擬在查明登封市地質(zhì)災(zāi)害發(fā)育概況及其類型的基礎(chǔ)上，闡明登封市地質(zhì)災(zāi)害的發(fā)育特征、分布規(guī)律，分析登封市地質(zhì)災(zāi)害的形成條件，運用GIS、層次分析法實現(xiàn)登封市地質(zhì)災(zāi)害危險性評價。

1登封市地質(zhì)災(zāi)害發(fā)育概況

登封市位于河南省中西部，穎河上游，北與鞏義、偃師相連，東與新密市相接，南鄰汝州、禹州，西與伊川接壤。區(qū)內(nèi)東西長55 km，南北寬36 km。登封市交通便利，境內(nèi)有鄭少高速公路、少洛高速公路、207國道、237省道等。鄉(xiāng)村公路已基本網(wǎng)絡(luò)化，交通十分方便。

登封市地處嵩山地區(qū)，地勢南部和北部高，中部低緩，境內(nèi)主要以低山、構(gòu)造剝蝕丘陵為主。巖性以石英巖、石英砂巖、灰?guī)r為主，第四系零星分布在境內(nèi)近山前地帶、山區(qū)較大溝谷以及穎河河谷，隨著人類工程活動的增加及降雨的陣法性影響，使該市地質(zhì)災(zāi)害頻頻發(fā)生。據(jù)走訪調(diào)查發(fā)現(xiàn)，該市地災(zāi)類型有滑坡、崩塌、地面塌陷和泥石流4類共計289處，其中滑坡74處，占總數(shù)的25.61 %；崩塌169處，占總數(shù)的58.47%；地面塌陷37處，占總數(shù)的12.80%；泥石流9處，占總數(shù)的3.11%(圖1)。

圖1登封市地質(zhì)災(zāi)害點餅狀圖

近年來，隨著經(jīng)濟(jì)的發(fā)展，人口增長，在城鎮(zhèn)建設(shè)、交通及旅游景區(qū)建設(shè)，特別是礦業(yè)開發(fā)、群眾削坡建窯建房等工程活動中，誘發(fā)的地質(zhì)災(zāi)害事件呈上升趨勢，因此，對區(qū)內(nèi)的災(zāi)害治理要求十分迫切。同時，區(qū)內(nèi)地質(zhì)災(zāi)害危險性分區(qū)對其地質(zhì)災(zāi)害的防治工作也具有重要指導(dǎo)意義。目前，層次分析法是地質(zhì)災(zāi)害危險性分區(qū)評價中應(yīng)用最為廣泛的方法之一[2]。

2評價方法及評價過程

層次分析法(analytic hierarchy process,AHP)是美國運籌學(xué)家Saaty教授于20世紀(jì)70年代初提出的。該方法是一種定性與定量分析的多目標(biāo)決策分析方法，能夠?qū)Q策者的經(jīng)驗判斷給予量化，廣泛應(yīng)用于目標(biāo)結(jié)構(gòu)復(fù)雜且缺乏必要數(shù)據(jù)情況下的分析與決策[4]。其主要思想是首先按照指標(biāo)間的相互關(guān)聯(lián)影響以及隸屬關(guān)系，將因子依不同層次聚集組合，形成一個多層次的分析結(jié)構(gòu)模型；其次根據(jù)對礦山客觀現(xiàn)象的主觀判斷，就每一層因子的相對重要性給予量化描述；最后，利用數(shù)學(xué)方法，確定每一層次全部因子相對重要次序的數(shù)值[3]。

2.1評價因子分析與量化

2.1.1評價因子選取

登封市地質(zhì)災(zāi)害主要以滑坡和崩塌為主，影響其危險性的因子有很多，根據(jù)形成的地質(zhì)環(huán)境條件、發(fā)育特征等方面分析，最后篩選出災(zāi)害點密度(x1)、巖土體類型(x2)、斷裂發(fā)育程度(x3)、地貌類型(x4)、坡度(x5)、降雨量(x6)、人類工程活動(x7)等7個地質(zhì)災(zāi)害危險性評價因子。

2.1.2構(gòu)造判斷矩陣

為了量化各判別因子的權(quán)重，判別因子之間需形成一個兩兩對比的矩陣，采用Saaty T.L提出的1～9比較標(biāo)度法(表1)來表示因子之間的影響強弱。通過相關(guān)工程的實際情況，確定各判別因子之間的相對重要性并賦與相應(yīng)的分值，最后得到所需的判斷矩陣(表2)。

2.1.3評價因子權(quán)重系數(shù)確定

采用層次分析法來確定本市地質(zhì)災(zāi)害評價因子的權(quán)重。

可分4個步驟：

表1　1～9標(biāo)度的意義表

表2　判斷矩陣U

第4步： 一致性檢驗針對階數(shù)大于2階的矩陣，只有當(dāng)隨機(jī)一致性比值CR≤0.1時判讀矩陣才具有一致性。隨機(jī)一致性比值的計算公式如下：

其中，公式中λmax為矩陣最大特征值；n為矩陣階數(shù)；R.I為平均一致性指標(biāo)(根據(jù)階數(shù)n查表3可得)。當(dāng)CR<0.1，認(rèn)為判斷矩陣的一致性是可以接受的，表明權(quán)重(表4)計算合理。

2.1.4評價因子分級賦值

地質(zhì)災(zāi)害危險性評價因子的分級賦值是根據(jù)本市實際情況，依據(jù)專家系統(tǒng)法收集調(diào)查獲取的，賦值越高，表示該等級對本市地質(zhì)災(zāi)害危險性的貢獻(xiàn)率越大[4]。具體分級賦值見表5。

2.2評價的模型建立

登封市的地質(zhì)災(zāi)害危險性評價方法主要采用綜合指數(shù)評價方法[5]，由綜合評價指數(shù)計算獲取，其計算公式為：

式中,W為地質(zhì)災(zāi)害危險性綜合指數(shù)；ai為因子權(quán)重；Fi為因子的賦值。

因此，可得出本區(qū)地質(zhì)災(zāi)害危險性評價的數(shù)學(xué)模型為：

W=0.19F1+0.18F2+0.12F3+0.13F4+0.12F5+0.11F6+0.15F7

表3　平均一致性指標(biāo)RI[5]

表4　評價因子權(quán)重值表

表5　地質(zhì)災(zāi)害評價因子賦值表

3登封市地質(zhì)災(zāi)害危險性綜合分區(qū)

3.1評價單元劃分

斜坡單元是滑坡、崩塌等地質(zhì)災(zāi)害發(fā)育的基本單元，并且在各類控制或影響因素中，河流和溝谷的發(fā)育階段對滑坡、崩塌的形成具有明顯的控制作用，且登封市的地質(zhì)災(zāi)害類型以崩塌、滑坡為主，因此本文采用基于幼年期溝谷劃分的斜坡單元作為評價單元，可以與地質(zhì)環(huán)境條件緊密聯(lián)系，綜合體現(xiàn)各類控制或影響因素的作用，使評價結(jié)果更貼近于實際[6]。本文借助GIS水文分析工具從1∶50 000的地形數(shù)據(jù)中提取劃分斜坡單元，根據(jù)登封市1∶50 000地形圖生成DEM，將幼年期溝谷中的三級支流坳溝、沖溝劃分為1 667個單元。

3.2危險性綜合分區(qū)

在前述評價因子和評價單元劃分的基礎(chǔ)上，將各因子區(qū)特征屬性添加到單元網(wǎng)格中，使區(qū)域內(nèi)所有網(wǎng)格帶有所有因子的特征屬性；利用Arcgis的統(tǒng)計和分析功能，提取每個評價單元的各項指標(biāo)值，對各評價因子進(jìn)行屬性疊加，利用公式(1)計算各評價單元的危險性指標(biāo)值，根據(jù)指標(biāo)值突變點將區(qū)域劃分為低危險區(qū)、中危險區(qū)和高危險區(qū)3個等級(圖2)。在定量分析的基礎(chǔ)上，結(jié)合野外實地調(diào)查，進(jìn)一步勾繪出危險程度分區(qū)圖(圖3、表6)。

圖2登封市評價單元劃分圖

圖3　地質(zhì)災(zāi)害危險性分區(qū)圖

危險性分區(qū)區(qū)亞區(qū)分布位置面積/km2分布特征高危險區(qū)ⅢⅢ1君召鄉(xiāng)-石道鄉(xiāng)-大金店鎮(zhèn)-嵩陽辦-少林辦-唐莊鄉(xiāng)(Ⅲ1)亞區(qū)149.49區(qū)內(nèi)共發(fā)育地質(zhì)災(zāi)害點38處,其中滑坡10處,泥石流8處,不穩(wěn)定斜坡12處,崩塌8處;該區(qū)域景區(qū)建設(shè)活動強烈,流動人口密度較大,分布有少林寺、初祖庵、塔林、少室闕、嵩陽書院、會善寺等文物建筑單體和建筑群,人類工程活動較強烈,景區(qū)建設(shè)、修路切坡等工程活動較多,極易發(fā)生災(zāi)害,危害嚴(yán)重,潛在地質(zhì)災(zāi)害也發(fā)育,對古建筑群危害大Ⅲ2送表鄉(xiāng)-白坪鄉(xiāng)-徐莊鎮(zhèn)亞區(qū)224.99區(qū)內(nèi)共發(fā)育地質(zhì)災(zāi)害點84處,其中地面塌陷10處,滑坡18處,不穩(wěn)定斜坡11處,崩塌45處;該區(qū)人口密集,經(jīng)濟(jì)發(fā)展較快;經(jīng)濟(jì)主要以礦業(yè)開采為主,煤礦企業(yè)約有二十幾家,石料開采隨處可見,地貌景觀破壞強烈,因煤礦開采形成的采空區(qū)逐年增大,形成采空區(qū)面積較大,地面出現(xiàn)塌陷,房屋出現(xiàn)裂縫,影響居民生活;G207兩側(cè)經(jīng)?？梢娦⌒捅浪?告城鎮(zhèn)-大冶鎮(zhèn)-宣化鎮(zhèn)亞區(qū)130.37區(qū)內(nèi)共發(fā)育地質(zhì)災(zāi)害點41處,其中地面塌陷13處,滑坡7處,不穩(wěn)定斜坡4處,崩塌17處;該區(qū)分布有觀星臺古建筑,人類工程活動強烈,永登高速、臨白鐵路、省道237由此通過,區(qū)內(nèi)大冶鎮(zhèn)、告城鎮(zhèn)經(jīng)濟(jì)以礦業(yè)開采為主,地質(zhì)地貌景觀破壞嚴(yán)重,地面塌陷、崩塌、滑坡較為發(fā)育,地質(zhì)災(zāi)害如果發(fā)生影響嚴(yán)重,影響居民生活、道路交通,同時這些地質(zhì)災(zāi)害的發(fā)生直接或間接威脅著古建筑群的安全中危險區(qū)ⅡⅡ1穎陽鎮(zhèn)-君召鄉(xiāng)亞區(qū)140.29區(qū)內(nèi)共發(fā)育地質(zhì)災(zāi)害點9處,其中地面塌陷3處,滑坡2處,不穩(wěn)定斜坡1處,崩塌3處;該區(qū)域涉及地貌單元主要為低山、剝蝕丘陵區(qū)。黃土覆蓋丘陵區(qū)人口密集、公路和鐵路等人類活動強烈,地表為Qp3黃土狀粉土、黃土狀粉質(zhì)粘土,下伏基巖,基巖風(fēng)化強烈,穩(wěn)定性較差,在外力作用下易于形成滑坡,如果發(fā)生地質(zhì)災(zāi)害,危害較大,危險性較高;低山區(qū)上覆黃土狀粉質(zhì)粘土較薄,下伏基巖出露較淺,基巖產(chǎn)狀與坡體傾向一致,在外力作用下易于發(fā)生滑坡,但人口較少,危害性雖高,但危害較小Ⅱ2唐莊鄉(xiāng)-盧店鎮(zhèn)-大冶鎮(zhèn)亞區(qū)152.52區(qū)內(nèi)共發(fā)育地質(zhì)災(zāi)害點39處,其中地面塌陷4處,滑坡12處,不穩(wěn)定斜坡3處,崩塌19處;該區(qū)人類工程活動強烈,鄭盧高速、237省道由此通過,居民主要沿路而居,人口密集,地質(zhì)地貌景觀破壞嚴(yán)重,崩塌、滑坡較為發(fā)育,地質(zhì)災(zāi)害如果發(fā)生影響較嚴(yán)重,影響居民生活、道路交通低危險區(qū)ⅠⅠ1君召鄉(xiāng)、石道鄉(xiāng)、大金店鎮(zhèn)、東華鎮(zhèn)、中岳辦、告城鎮(zhèn)420.38區(qū)內(nèi)共發(fā)育地質(zhì)災(zāi)害點78處,其中地面塌陷7處,滑坡14處,不穩(wěn)定斜坡11處,崩塌43處;中西部地勢相對平坦,起伏較小,地表為黃土覆蓋,人口密度大,交通網(wǎng)絡(luò)發(fā)達(dá),鄭盧高速由此通過,局部有切坡建房修路造成邊坡失穩(wěn),該區(qū)分布有中岳廟、啟母闕、太室闕古建筑群,如發(fā)生災(zāi)害對古建筑群發(fā)生威脅;北部山區(qū)氣候濕潤,植被覆蓋率較高,大部分基巖出露,主要以農(nóng)田耕作為主,人類工程活動微弱,但人口稀少,居住分散,地質(zhì)災(zāi)害發(fā)生造成的危害較小

4結(jié)論

(1) 層次分析法對影響登封市危險性的相關(guān)指標(biāo)進(jìn)行定量化分析，同時，根據(jù)當(dāng)?shù)貙嶋H情況，調(diào)整了評價分級標(biāo)準(zhǔn)并重新分配了指標(biāo)權(quán)重[2]，然后將數(shù)據(jù)代入典型災(zāi)害評價模型，進(jìn)行地質(zhì)災(zāi)害危險性評價，評價結(jié)果與實際情況相符。研究表明，基于層次分析法所建立的評價模型能較好地應(yīng)用于登封市地質(zhì)災(zāi)害危險性評價中，該評價方法是實用和有效的。

(2) 從地質(zhì)災(zāi)害規(guī)模上看，南北中低山區(qū)、東部采礦區(qū)以及中部沖溝兩側(cè)地質(zhì)災(zāi)害規(guī)模較大，以滑坡、崩塌、地面塌陷為主，這些地質(zhì)災(zāi)害給當(dāng)?shù)鼐用裨斐芍苯咏?jīng)濟(jì)損失約27 011萬元。隨著經(jīng)濟(jì)社會的迅猛發(fā)展，人類工程活動不斷增強，無論是在深度上還是在廣度上都將日益加劇，對自然環(huán)境的破壞也相應(yīng)日趨嚴(yán)重，從而誘發(fā)許多危害極嚴(yán)重的地質(zhì)災(zāi)害。

參考文獻(xiàn)

[1]章正高.采用層次分析法對地質(zhì)災(zāi)害進(jìn)行危險性評價——以漢陰縣地質(zhì)災(zāi)害為例.西安科技大學(xué)學(xué)報,2014,34(1)：87-91．

[2]王小江,孫書勤,盧鵬飛.層次分析法在重慶某縣地質(zhì)災(zāi)害危險性評價中的應(yīng)用[J].地球與環(huán)境,2014,42(3):419-423．

[3]羅娟,陳守余.礦山環(huán)境質(zhì)量評價指標(biāo)體系及層次分析法評價[J].安全與環(huán)境工程,2005,12(1):9-12．

[4]王哲,易發(fā)成.基于層次分析法的綿陽市地質(zhì)災(zāi)害易發(fā)性評價.自然災(zāi)害學(xué)報,2009,18(1):14-23．

[5]尚慧，倪萬魁，等.斜坡單元劃分在彭陽縣地質(zhì)災(zāi)害危險性區(qū)劃中的應(yīng)用.中國水土保持SWCC,2011,3:48-50.

[6]霍艾迪，張駿，等.地質(zhì)災(zāi)害易發(fā)性評價單元劃分方法——以陜西省黃陵縣為例.吉林大學(xué)學(xué)報(地球科學(xué)版),2011,41(2):523-535.

AHP-BASED ASSESSMENT OF GEOLOGICAL HAZARDS RISK IN DENGFENG CITY

DING Li,LI Zai-xing,LI Xue-wen,ZHANG Rui-xiang

(Mapping and Geoinformation of Henan Provincial Bureau of Geo-exploration and Mineral Development，Zhengzhou450006,China)

Abstract:On the basis of the geological environment and development characteristics in a county,seven indicators for evaluating the geological disaster fatalness around the county,density,rock soil mass type,fault developmental degree,geomorphic type,gradient,rainfall and engineering activities are selected.Using the analysis hierarchy process to ascertain the weight of index,and adopting the export system to mark for the indicators,and then,the eveluation model for the geological disaster fatalness along the city is put forward.At last,the area is divided into three risk regions:the most risk region,the medium risk region,and low risk region based on Arcgis spatial analysis,and the results of the study provide some certain referential significance for similar areas.

Key words：AHP;geological disaster;division of disaster fatalness;evaluation of geological disaster fatalness

文章編號：1006-4362(2016)02-0092-05

收稿日期：2015-12-15改回日期：2016-02-23

基金項目：河南省2013年度1∶5萬地質(zhì)災(zāi)害調(diào)查項目(執(zhí)行項目編號：19)

中圖分類號：X43;X141

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

作者簡介：丁麗(1984-)，女，甘肅白銀人，漢族，碩士，河南省地質(zhì)礦產(chǎn)勘查開發(fā)局測繪地理信息院工程師，從事環(huán)境地質(zhì)、地質(zhì)災(zāi)害防治方面的研究。E-mail:83381975@qq.com