王星薈，孟凡奇

(1.中國地質(zhì)大學(武漢)環(huán)境學院，湖北 武漢　430074；2.山東省地質(zhì)環(huán)境監(jiān)測總站，山東 濟南　250014)

可拓學在濟南市地質(zhì)災(zāi)害易發(fā)性評價中的應(yīng)用

王星薈1，孟凡奇2

(1.中國地質(zhì)大學(武漢)環(huán)境學院，湖北 武漢430074；2.山東省地質(zhì)環(huán)境監(jiān)測總站，山東 濟南250014)

介紹了可拓學的基本理論，并應(yīng)用于地質(zhì)災(zāi)害易發(fā)性劃分。選取了巖土體類型、地形坡度、地質(zhì)構(gòu)造、已有災(zāi)害點數(shù)量、人類工程活動強度、多年平均降雨量為地質(zhì)災(zāi)害易發(fā)性評價指標，運用層次分析法獲得了評價指標的權(quán)重值，建立了基于可拓學理論的地質(zhì)災(zāi)害易發(fā)性評價模型。以濟南市為例，對可拓學模型的評價結(jié)果進行驗證，該方法獲得結(jié)果與實際情況吻合程度高，證明這一方法是切實可行的。

可拓學；地質(zhì)災(zāi)害；易發(fā)性評價；濟南

引文格式：王星薈，孟凡奇.可拓學在濟南市地質(zhì)災(zāi)害易發(fā)性評價中的應(yīng)用[J].山東國土資源，2015,31(4)：67-71. WANG Xinghui，MENG Fanqi. Application of Extension Method of Probability Evaluation of Geological Hazarads in Jinan City[J].Shandong Land and Resources, 2015,31(4):67-71.

0　引言

我國的地質(zhì)災(zāi)害種類多、發(fā)災(zāi)頻率高、威脅人口多、造成經(jīng)濟損失大。地質(zhì)災(zāi)害易發(fā)性評價是地質(zhì)災(zāi)害防治工作的要點，其評價結(jié)果可為地方政府的地質(zhì)災(zāi)害防治規(guī)劃、預(yù)警預(yù)報、經(jīng)濟發(fā)展規(guī)劃等工作提供重要的科學依據(jù)。國內(nèi)外學者多采用專家經(jīng)驗法、打分法、綜合危險性指數(shù)法、信息量法等方法進行研究，并取得了豐碩的成果[1-5]。但是傳統(tǒng)方法受研究人員的專業(yè)知識水平不同造成對研究區(qū)認知程度差異的影響，往往選用的評價模型運算精度不高，造成評價結(jié)果與實際吻合程度較低的情況。因此建立一種既能通過將評價因子數(shù)據(jù)分級來幫助研究者認知地質(zhì)災(zāi)害發(fā)育程度，又能提高評價精度的地質(zhì)災(zāi)害模型十分必要。

可拓學是以廣東工業(yè)大學的蔡文研究員為首的中國學者創(chuàng)立的新學科，該方法在泥石流危險度評價、邊坡穩(wěn)定性、以及工業(yè)、醫(yī)學經(jīng)濟等領(lǐng)域應(yīng)用十分廣泛[6]。該文將可拓學運用到地質(zhì)災(zāi)害易發(fā)性評價中，采取了層次分析法對易發(fā)性影響因子進行權(quán)重計算，構(gòu)建了基于可拓學的地質(zhì)災(zāi)害易發(fā)性評價模型，并以濟南市地質(zhì)災(zāi)害為例進行驗證。

1　可拓方法

1.1物元理論

可拓學理論是以事物、特征及事物關(guān)于該特征的量值三者所組成三元組，記作R(事物、特征、量值)[7]，可將諸多同征物元用簡單的方法表示出來：

(1)

式中：R為n維物元；Ni為事物的名稱；ai為特征；vi為特征ai所對應(yīng)的量值，合稱為物元的三要素。

1.2確定經(jīng)典域

(2)

式中：N0j表示地質(zhì)災(zāi)害易發(fā)性評價的第j(j=1，2，3，…，m)個等級；ci(i=1，2，3，…，n)表示地質(zhì)災(zāi)害易發(fā)性評價等級N0j的第j個評價因子；V0ji=為N0j關(guān)于因素ci所確定的量值范圍，即地質(zhì)災(zāi)害易發(fā)性評價各等級關(guān)于對應(yīng)的評價因子所取的數(shù)據(jù)范圍——經(jīng)典域[8]。

1.3確定節(jié)域

(3)

式中：P為地質(zhì)災(zāi)害易發(fā)性評價等級的全體；Vpi=為P關(guān)于因素ci所取量值的范圍，即P的節(jié)域。

1.4確定待評物元

對于待評p，把所收集到的地理信息數(shù)據(jù)或分析結(jié)果用物元表示，即得到待評物元R。

(4)

式中：p為待評的目標區(qū)域，vi為p關(guān)于評價ci的量值，即待評區(qū)域的具體數(shù)據(jù)。

1.5確定權(quán)系數(shù)

該文采用層次分析法來進行權(quán)重計算，可按下面4個步驟進行：建立遞階層次結(jié)構(gòu)模型；構(gòu)造出各層次中的所有判斷矩陣；層次單排序及一致性檢驗；層次總排序及一致性檢驗。其核心內(nèi)容是通過對因子進行兩兩比較建立成對判斷矩陣，即每次取兩個因子xi和xj，以aij表示xi和xj對Z的影響大小之比，其全部比較結(jié)果用矩陣A=(aij)n×n表示。

1.6確定評價指標關(guān)于各等級的關(guān)聯(lián)度

各單項評價指標vi關(guān)于各等級j的關(guān)聯(lián)度為：

(5)

1.7待評物元關(guān)于各類別等級的關(guān)聯(lián)度

待評目標p關(guān)于評價等級j的關(guān)聯(lián)度為

(6)

1.8可拓評判等級

(7)

(8)

則稱j*為p的級別變量特征值，從j*中可看出偏向另一類的程度。

2　濟南市地質(zhì)災(zāi)害易發(fā)性評價

2.1濟南市地質(zhì)災(zāi)害概況

濟南市是山東省地質(zhì)災(zāi)害較發(fā)育的城市之一，地貌整體上呈南高北低之勢，南部地處泰山余脈為低山丘陵區(qū)，巖性以泰山巖群變質(zhì)巖和寒武-奧陶紀灰?guī)r為主，北部為黃河沖積平原區(qū)。據(jù)統(tǒng)計，截至2012年底，濟南市共有地質(zhì)災(zāi)害點230處，其中崩塌133處，滑坡32處，泥石流10處，地面塌陷54處，地裂縫1處，其中崩塌、滑坡、泥石流災(zāi)害點主要分布在濟南市南部的低山丘陵區(qū)，地面塌陷、地裂縫主要分布在地下礦產(chǎn)開采區(qū)和第四系覆蓋層較薄，且地下水位波動較大的隱伏灰?guī)r區(qū)。

2.2評價思路

根據(jù)《縣市地質(zhì)災(zāi)害調(diào)查與區(qū)劃基本要求》將地質(zhì)災(zāi)害易發(fā)程度劃分為高易發(fā)區(qū)、中易發(fā)區(qū)、低易發(fā)區(qū)和不易發(fā)區(qū)，濟南市底圖按照1km×1km單元格進行劃分，每一個單元格視作一個待評物元，參考前人研究成果[9-10]，將每個單元格對應(yīng)的評價指標進行賦值，運用層次分析法對評價指標進行賦權(quán)，構(gòu)建可拓學模型進行分析評價，確定其對應(yīng)的等級。

2.3評價指標確定

根據(jù)濟南市地質(zhì)災(zāi)害實地調(diào)查中發(fā)現(xiàn)，該市地質(zhì)災(zāi)害發(fā)育程度存在受地形地貌、地質(zhì)構(gòu)造、地層巖性影響較大，多由大氣降雨所引發(fā)，且在削坡修路建房的地區(qū)尤甚的情況，參考前人在地質(zhì)災(zāi)害易發(fā)性評價時所總結(jié)的經(jīng)驗，選取了巖土體類型、地形坡度、地質(zhì)構(gòu)造、已有災(zāi)害點數(shù)量、人類工程活動強度、多年平均降雨量共6個地質(zhì)災(zāi)害易發(fā)性評價指標。由于各個評價指標的數(shù)值和單位不同，且有的指標是非定量的，因此必須對所有指標進行量化處理，將定性和定量描述的指標有機的結(jié)合起來，指標量化值見表1。量化后將評價指標歸一化處理，使數(shù)值滿足0～1之間，再參考其他文獻所劃定分類標準[11-13]，劃分標準見表2。

表1　地質(zhì)災(zāi)害易發(fā)性分區(qū)可拓方法評價指標量值

表2　地質(zhì)災(zāi)害易發(fā)評價各指標分類標準

2.4權(quán)重確定

該次研究中利用層次分析法來確定各因子的權(quán)重，根據(jù)層次分析法的基本步驟，將評價指標分為基礎(chǔ)因子、響應(yīng)因子和誘發(fā)因子，建立濟南市地質(zhì)災(zāi)害易發(fā)性各評價指標的遞階模型(圖1)。

圖1　地質(zhì)災(zāi)害影響因子遞階層次模型

經(jīng)計算并通過一致性檢驗，C1，C2，C3，C4，C5，C6的權(quán)重分別為0.285，0.082，0.203，0.140，0.194，0.096。

2.5確定經(jīng)典域及節(jié)域

根據(jù)表2數(shù)據(jù)確定經(jīng)典域R01，R02，R03，R04和節(jié)域Rp。

(1)經(jīng)典域

(9)

(10)

(11)

(12)

(2)節(jié)域

(13)

2.6可拓學評價

將濟南市按照1km×1km劃分的每一個單元格作為一個待評單元，對其每一項指標進行量化賦值，根據(jù)表2、表3和所有待評物元及經(jīng)典域、節(jié)域、層次分析法確定權(quán)重系數(shù)，進行編程計算，獲得每個單元的地質(zhì)災(zāi)害易發(fā)程度可拓學評價結(jié)果。例如某一單元格的數(shù)據(jù)為：

(14)

經(jīng)判定，R1為Ⅲ級，對應(yīng)的地質(zhì)災(zāi)害易發(fā)級別為中易發(fā)區(qū)。通過可拓學模型計算研究區(qū)內(nèi)的所有單元格的易發(fā)等級，從而獲得濟南市地質(zhì)災(zāi)害易發(fā)性分區(qū)評價圖(圖2)。

圖2　濟南市地質(zhì)災(zāi)害易發(fā)性分區(qū)圖

2.7應(yīng)用分析

濟南市南部為低山丘陵區(qū)，具有地質(zhì)構(gòu)造較發(fā)育、災(zāi)害分布密度密集、地形高差大、災(zāi)害發(fā)育程度較大的特點，與地質(zhì)災(zāi)害中、低易發(fā)區(qū)評價結(jié)果相符合。濟南北部地區(qū)地形平坦，受地形地貌影響，崩塌滑坡泥石流不發(fā)育，雖有地下采礦和地下水集中開采活動，但形成的地面塌陷、地裂縫災(zāi)害點少，與地質(zhì)災(zāi)害低發(fā)育和不易發(fā)區(qū)評價結(jié)果吻合程度較高。地質(zhì)災(zāi)害易發(fā)性評價結(jié)果和實地發(fā)育情況吻合程度高，證明了選取的評價因子和權(quán)重的賦值是合理的，也對應(yīng)了濟南地質(zhì)災(zāi)害整體上發(fā)育程度不高，災(zāi)害易發(fā)性南強北弱的特點。

值得注意的是，受地理位置的影響，其他地區(qū)在選取地質(zhì)災(zāi)害易發(fā)性評價因子和權(quán)重賦值時，應(yīng)充分考慮當?shù)貙嶋H情況進行適當調(diào)整因子數(shù)量和權(quán)重大小。

3　結(jié)論

該文利用可拓學評價理論，將地質(zhì)災(zāi)害易發(fā)性評價因子進行了合理分級，并提高了研究者對評價區(qū)地質(zhì)災(zāi)害發(fā)育情況的認識程度。建立了基于可拓學的地質(zhì)災(zāi)害易發(fā)性評價模型，并運用于濟南地區(qū)。研究結(jié)果表明,該方法確定的評價等級真實地反映了該區(qū)地質(zhì)災(zāi)害易發(fā)程度,評價精度較高,應(yīng)用于地質(zhì)災(zāi)害易發(fā)性評價是可行的、可靠的。

[1]張有全，宮輝力，李巧剛.基于GIS的北京市延慶縣地質(zhì)災(zāi)害易發(fā)性區(qū)域劃分[J].中國地質(zhì)災(zāi)害與防治學報，2006,17(4)：17-20.

[2]高治群，薛傳東，尹飛，等.基于GIS的信息量法及其地質(zhì)災(zāi)害易發(fā)性評價中應(yīng)用-以滇中晉寧縣為例[J]．2010,46(6)：1112-1116.

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

[4]阮沈勇，黃潤秋．基于GIS的信息量法模型在地質(zhì)災(zāi)害危險性區(qū)劃中的應(yīng)用[J]．成都理工學院學報，2001，28(1)：89-92．

[5]霍艾迪，張俊，盧玉東，等.地質(zhì)災(zāi)害易發(fā)性評價單元劃分方法——以陜西省黃陵縣為例[J]．吉林大學學報(地球科學版)，2011，41(2)：58-64.

[6]蔡文,楊春燕,林偉初.可拓工程方法[M].北京:科學出版社,2000.

[7]蔡文,楊春燕,王光華.一門新的交叉學科——可拓學[J].中國科學基金,2004,(5):268-272.

[8]蔡文.可拓論及其應(yīng)用[J].科學通報,1999,44(7):673-682.

[9]王學，張阻陸，張超，等.山東省地質(zhì)環(huán)境質(zhì)量評價研究[J].山東國土資源，2011，27(8)：14-17.

[10]張以晨,佴磊,孟凡奇，等.基于最優(yōu)組合賦權(quán)理論的可拓學評價模型在地質(zhì)災(zāi)害易發(fā)區(qū)劃分中的應(yīng)用[J].吉林大學學報(地球科學版)，2011，41(4)：1110-1115.

[11]尚敏,陳劍平,王征亮，等.向家河大橋庫岸地質(zhì)災(zāi)害危險性分區(qū)的可拓學評價[J].巖土力學,2007,28(11):2445-2450.

[12]謝曉娟．地質(zhì)環(huán)境條件及誘發(fā)因素在地質(zhì)災(zāi)害易發(fā)分區(qū)中定量化評價初探[J]．探礦工程，2003，(增刊)：93-95．

[13]呂勇，李高良，張貴，等.彝良縣地質(zhì)災(zāi)害分布特征及防治建議[J].中國巖溶，2009，28(2)：219-224.

ApplicationofExtensionMethodofProbabilityEvaluationofGeologicalHazaradsinJinanCity

WANGXinghui，MENGFanqi

(1.EnvironmentCollegeofChinaUniversityofGeosciences,HubeiWuhan430074，China；2.ShandongMonitoringCenterofGeologicalEnvironment,ShandongJinan250014,China)

Basictheoryofextensionmethodisintroducedinthispaper.Itcanbeappliedingeologicalhazardprobabilityofoccurrenceassessment.Selectingrocksoilmasstypes,hillslope,geologicalstructure,thenumberofdisaster,intensityofhumanengineeringactivityandmulti-yearaverageprecipitationastheevaluationindexs,usingthemethodofanalytichierarchyprocessforweightvalue,susceptibilityevaluationmodelhasbeenestablishedbasedonextensiontheoryofgeologicaldisaster.TakingJi'nancityasanexample,theextensionevaluationmodelisusedtovalidatetheresults.Theresultsobtainedbyusingthismethodcanagreewiththeactualsituationhighly.Itisprovedthatthismethodisfeasible.

extensiontheory；geologicalhazard；evaluationofprobabilityofoccurrence；Jinancity

2014-06-18；

2014-09-29；編輯：曹麗麗

王星薈(1994—)，女，江蘇高郵人，在校學生，地下水科學與工程專業(yè)；E-mail：619637319@qq.com

TD854.6

B