馮衛(wèi)，唐亞明，馬紅娜，宿曉虹，陳新建

(1.長(zhǎng)安大學(xué)地質(zhì)工程與測(cè)繪學(xué)院，陜西 西安 710064；2.自然資源部黃土地質(zhì)災(zāi)害重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710054；3.中國(guó)地質(zhì)調(diào)查局西安地質(zhì)調(diào)查中心，陜西 西安 710054)

當(dāng)今社會(huì)經(jīng)濟(jì)發(fā)展和人口數(shù)量急劇增長(zhǎng)，城市內(nèi)的各類資源環(huán)境矛盾日益突出，自然災(zāi)害對(duì)城市造成的影響也越來(lái)越大，并逐漸成為當(dāng)今人類社會(huì)面臨的最主要問題之一(閃淳昌，2005)。許多學(xué)者針對(duì)各單項(xiàng)災(zāi)種開展了一系列深入的風(fēng)險(xiǎn)評(píng)價(jià)和風(fēng)險(xiǎn)管理研究(Fell R，2008；吳樹仁，2012；唐亞明，2015)。隨著當(dāng)今世界地質(zhì)災(zāi)害、洪水災(zāi)害和地震災(zāi)害等多種自然災(zāi)害頻發(fā)，災(zāi)害群聚和群發(fā)現(xiàn)象導(dǎo)致災(zāi)情加重，人們?cè)絹?lái)越關(guān)注如何降低多災(zāi)種自然災(zāi)害的風(fēng)險(xiǎn)。

近幾十年來(lái)，國(guó)內(nèi)外學(xué)者針對(duì)多災(zāi)種形成的災(zāi)害風(fēng)險(xiǎn)開展了一定程度的研究。美國(guó)聯(lián)邦應(yīng)急管理局(FEMA)與國(guó)家建筑科學(xué)研究院(NIBS)聯(lián)合開發(fā)了基于GIS平臺(tái)的多災(zāi)種風(fēng)險(xiǎn)評(píng)估軟件包HAZUS-MH，對(duì)美國(guó)各個(gè)行政區(qū)開展了洪水災(zāi)害、地震災(zāi)害和颶風(fēng)災(zāi)害的綜合風(fēng)險(xiǎn)評(píng)估(FEMA，2004)。歐洲委員會(huì)聯(lián)合研究中心提出了JRC綜合風(fēng)險(xiǎn)評(píng)估法，并對(duì)歐洲多個(gè)國(guó)家開展了天氣災(zāi)害的綜合風(fēng)險(xiǎn)評(píng)估(Lavalle C et al.，2005)。史培軍對(duì)災(zāi)害系統(tǒng)的性質(zhì)、動(dòng)力學(xué)機(jī)制、綜合減災(zāi)范式以及災(zāi)害科學(xué)體系等進(jìn)行了一系列探討，構(gòu)建了由災(zāi)害科學(xué)、應(yīng)急技術(shù)和風(fēng)險(xiǎn)管理共同組成的“災(zāi)害風(fēng)險(xiǎn)科學(xué)”學(xué)科體系(史培軍，1996，2009)。葛全勝等在系統(tǒng)論述國(guó)內(nèi)外自然災(zāi)害風(fēng)險(xiǎn)研究理論方法的基礎(chǔ)上進(jìn)行了地級(jí)市、縣、村三級(jí)尺度的自然災(zāi)害綜合風(fēng)險(xiǎn)評(píng)估研究(葛全勝，2008)。蓋程程等為了解決多災(zāi)種耦合的風(fēng)險(xiǎn)評(píng)估問題，提出了針對(duì)自然災(zāi)害和事故災(zāi)難耦合的風(fēng)險(xiǎn)評(píng)估方法(蓋程程，2011)。明曉東等在厘定多災(zāi)種以及多災(zāi)種風(fēng)險(xiǎn)基本概念的基礎(chǔ)上，系統(tǒng)綜述了國(guó)內(nèi)外有較大影響的多災(zāi)種風(fēng)險(xiǎn)評(píng)估方法及其分類體系(明曉東，2013)。倪曉嬌等以長(zhǎng)白山地區(qū)為研究區(qū)，采用極差法、層次分析法和綜合指數(shù)法進(jìn)行了基于多災(zāi)種自然災(zāi)害風(fēng)險(xiǎn)水平的長(zhǎng)白山地區(qū)生態(tài)安全綜合評(píng)價(jià)(倪曉嬌，2014)。

目前，陜西省內(nèi)各縣市的多災(zāi)種風(fēng)險(xiǎn)研究尚處于空白狀態(tài)。筆者以陜西省咸陽(yáng)市為例，將咸陽(yáng)市地質(zhì)災(zāi)害、洪水災(zāi)害和地震災(zāi)害結(jié)合起來(lái)，開展多災(zāi)種自然災(zāi)害風(fēng)險(xiǎn)綜合評(píng)價(jià)，以期為咸陽(yáng)市應(yīng)急管理和國(guó)土空間規(guī)劃提供一定的理論指導(dǎo)。

1 研究區(qū)概況

咸陽(yáng)市位于陜西省中部，東與銅川市、渭南市為鄰，西與寶雞市接壤，北同甘肅省慶陽(yáng)市、平?jīng)鍪信B，南接西安市，轄9縣、2區(qū)、2市，面積為10 196.35 km2。咸陽(yáng)市屬大陸性季風(fēng)氣候，多年平均降水量為518.27 mm，最大年降水量為812.33 mm，最小年降水量為369.42 mm，渭河從咸陽(yáng)市南界自西向東流過，涇河自西北向東南流經(jīng)全境。全市以橫亙于中部的嵯峨山、九峻山和五峰山為界分為地貌景觀迥異的南北2個(gè)區(qū)。北部總體為黃土高原南緣，除轉(zhuǎn)角、石門一帶為子午嶺南端石質(zhì)低山區(qū)外，主要為黃土殘塬梁峁溝壑區(qū)；南部為渭河平原。區(qū)內(nèi)地層巖性由老到新依次為寒武系、奧陶系、石炭系、二疊系、三疊系、侏羅系、白堊系、古近系、新近系和第四系。在區(qū)域地質(zhì)構(gòu)造上屬于中朝準(zhǔn)地臺(tái)西南部的二級(jí)地質(zhì)構(gòu)造單元渭河斷陷帶內(nèi)，區(qū)域地震構(gòu)造上屬于華北地震區(qū)。

現(xiàn)階段對(duì)咸陽(yáng)市影響較為重大的災(zāi)害類型包括地質(zhì)災(zāi)害、洪水災(zāi)害和地震災(zāi)害3種類型。地質(zhì)災(zāi)害以崩塌、滑坡、泥石流為主；洪水災(zāi)害以山前山洪和平原區(qū)內(nèi)澇為主；地震災(zāi)害以地裂縫、地表斷層、砂土液化、黃土震陷為主(圖1)。筆者從這3類災(zāi)害入手，開展相應(yīng)的自然災(zāi)害綜合風(fēng)險(xiǎn)評(píng)價(jià)。

2 總體思路和研究方法

2.1 總體思路

首先，以柵格單元作為評(píng)價(jià)單元，兼顧咸陽(yáng)市已有的各類自然災(zāi)害發(fā)育情況，針對(duì)每類災(zāi)種分別采用層次分析法(analytic hierarchy process，簡(jiǎn)稱AHP)進(jìn)行計(jì)算，求出各類災(zāi)種的危險(xiǎn)性分析結(jié)果。其次，在上述各類危險(xiǎn)性分析結(jié)果的基礎(chǔ)上，再次采用層次分析法進(jìn)行計(jì)算，得出多災(zāi)種危險(xiǎn)性綜合評(píng)價(jià)結(jié)果。最后，在多災(zāi)種危險(xiǎn)性綜合分析結(jié)果的基礎(chǔ)上，對(duì)危害后果(人員、財(cái)產(chǎn)及易損性)進(jìn)行量化，將危險(xiǎn)性量化結(jié)果與危害后果量化結(jié)果進(jìn)行柵格乘積運(yùn)算，進(jìn)而得出多災(zāi)種自然災(zāi)害綜合風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果。

圖1 咸陽(yáng)市地質(zhì)環(huán)境與災(zāi)害分布簡(jiǎn)圖Fig.1 Sketch map of geological environment and disasters distribution in Xianyang city

2.2 層次分析法基本原理

層次分析法是比較常用的對(duì)各因子、指標(biāo)進(jìn)行賦權(quán)的方法之一(鄧雪等，2012)。該方法把復(fù)雜問題按照主次或支配關(guān)系分組而形成有序的遞階層次結(jié)構(gòu)，利用數(shù)學(xué)方法確定每一層中所有元素的相對(duì)重要性權(quán)值，最后通過排序結(jié)果來(lái)確定最終權(quán)重，是一種解決多目標(biāo)復(fù)雜問題的定性與定量相結(jié)合的決策分析方法。它包括：①建立遞階層級(jí)結(jié)構(gòu)。②構(gòu)造兩兩判斷矩陣。③由判斷矩陣計(jì)算被比較元素對(duì)于該準(zhǔn)則的相對(duì)權(quán)重，并進(jìn)行一致性檢驗(yàn)(表1)。

表1 兩兩比較矩陣Bij的取值及含義表Tab.1 The value and meaning of the pairwise comparison matrix Bij

以A表示準(zhǔn)則，Bi、Bj(i,j=1,2,…,n)表示因素，Bij表示Bi對(duì)Bj的相對(duì)重要性，則由Bij可組成一個(gè)兩兩判斷矩陣P=(Bij)n×n(Bij的取值及含義如表1)。對(duì)矩陣的每一行向量求幾何均值并歸一化，所得特征向量Wi即為各指標(biāo)權(quán)重。對(duì)判斷矩陣還需進(jìn)行一致性檢驗(yàn)，一致性指標(biāo)的計(jì)算公式為：

(1)

式中：λmax為判斷矩陣的最大特征值，計(jì)算公式為：

(2)

式中：(PW)i為向量PW的第i個(gè)元素；Wi為正規(guī)化處理的W的對(duì)應(yīng)元素。

最后檢驗(yàn)時(shí)使用一致性比例進(jìn)行判定，其公式為：

(3)

式中：RI為1～9階矩陣的平均隨機(jī)一致性指標(biāo)。若CR<0.1，認(rèn)為判斷矩陣的一致性是可以接受的；若CR>0.1，則應(yīng)對(duì)矩陣做適當(dāng)修正。

3 多災(zāi)種自然災(zāi)害危險(xiǎn)性分析

3.1 建立評(píng)價(jià)指標(biāo)模型

選取斜坡坡度、坡向、坡高、工程地質(zhì)巖組、地質(zhì)災(zāi)害發(fā)育程度和距公路距離6項(xiàng)主要因素作為地質(zhì)災(zāi)害的危險(xiǎn)性評(píng)價(jià)指標(biāo)；選取斜坡坡度、坡向、高程、工程地質(zhì)巖組、距水系距離和年平均降水量作為洪水災(zāi)害的危險(xiǎn)性評(píng)價(jià)指標(biāo)；選取斜坡坡度、工程地質(zhì)巖組、黃土震陷程度、砂土覆蓋區(qū)域、地震動(dòng)峰值加速度和距斷裂帶距離作為地震災(zāi)害的危險(xiǎn)性評(píng)價(jià)指標(biāo)。分別構(gòu)造各災(zāi)種的判斷矩陣，計(jì)算出最大特征值和對(duì)應(yīng)的特征向量，并進(jìn)行判斷矩陣一致性檢驗(yàn)。結(jié)果顯示均滿足一致性檢驗(yàn)，可作為計(jì)算用的權(quán)向量，各災(zāi)種的危險(xiǎn)程度量化評(píng)分表見表2、表3、表4。

3.2 各單災(zāi)種危險(xiǎn)性分析

按照各類災(zāi)害危險(xiǎn)性劃分的依據(jù)及要求，以“區(qū)內(nèi)相似、區(qū)間相異”為原則，并結(jié)合每個(gè)評(píng)價(jià)單元危險(xiǎn)性計(jì)算結(jié)果，將相應(yīng)的研究區(qū)劃分為高危險(xiǎn)區(qū)、中危險(xiǎn)區(qū)、低危險(xiǎn)區(qū)和極低危險(xiǎn)區(qū)4類(圖2、圖3、圖4)。

表2 咸陽(yáng)市地質(zhì)災(zāi)害危險(xiǎn)程度量化評(píng)分表Tab.2 Quantitative score table for geological disaster in Xianyang city

表3 咸陽(yáng)市洪水災(zāi)害危險(xiǎn)程度量化評(píng)分表Tab.3 Quantitative score table for flood disaster in Xianyang city

表4 咸陽(yáng)市地震災(zāi)害危險(xiǎn)程度量化評(píng)分表Tab.4 Quantitative score table for earthquake disaster in Xianyang city

圖2 咸陽(yáng)市地質(zhì)災(zāi)害危險(xiǎn)性分區(qū)評(píng)價(jià)圖Fig.2 Evaluation map of geological disaster analysis in Xianyang city

圖3 咸陽(yáng)市洪水災(zāi)害危險(xiǎn)性分區(qū)評(píng)價(jià)圖Fig.3 Evaluation map of flood disaster analysis in Xianyang city

從整體上看，地質(zhì)災(zāi)害高危險(xiǎn)區(qū)主要分布于北部的中低山丘陵區(qū)，尤其是在彬州市、長(zhǎng)武縣、旬邑縣和淳化縣為較集中，其他區(qū)域相對(duì)較低；洪水災(zāi)害高危險(xiǎn)區(qū)主要分布于渭河河流兩側(cè)武功縣至咸陽(yáng)市城區(qū)段沿線以及涇河河流兩側(cè)區(qū)域，其他區(qū)域相對(duì)較低；地震災(zāi)害高危險(xiǎn)區(qū)主要分布于渭河斷裂帶沿線(橫穿武功縣、興平市、秦都區(qū)及渭城區(qū))、扶風(fēng)-禮泉-富平斷裂帶局部段(禮泉縣縣城至三原縣城區(qū)沿線)、岐山-乾縣-蒲城斷裂帶局部段(乾縣縣城至淳化縣南段)以及口鎮(zhèn)-關(guān)山斷裂帶局部段(永壽縣東側(cè)至三原縣北側(cè)段)，其他區(qū)域相對(duì)較低。

3.3 多災(zāi)種危險(xiǎn)性綜合分析

在上述地質(zhì)災(zāi)害危險(xiǎn)性分析結(jié)果、洪水災(zāi)害危險(xiǎn)性分析結(jié)果和地震災(zāi)害危險(xiǎn)性分析結(jié)果的基礎(chǔ)上，采用層次分析法進(jìn)行計(jì)算，進(jìn)而得出多災(zāi)種危險(xiǎn)性綜合評(píng)價(jià)結(jié)果(圖5)。根據(jù)層次分析法計(jì)算結(jié)果，咸陽(yáng)市重大災(zāi)害危險(xiǎn)性分區(qū)的3個(gè)指標(biāo)權(quán)重(地質(zhì)災(zāi)害危險(xiǎn)性分區(qū)、洪水災(zāi)害危險(xiǎn)性分區(qū)和地震災(zāi)害危險(xiǎn)性分區(qū))依次為0.539、0.297、0.164。從3個(gè)指標(biāo)權(quán)重中可以看出，地質(zhì)災(zāi)害危險(xiǎn)性分區(qū)是影響咸陽(yáng)市重大災(zāi)害危險(xiǎn)性分區(qū)的最重要因素，其次為洪水災(zāi)害危險(xiǎn)性分區(qū)，最后是地震災(zāi)害危險(xiǎn)性分區(qū)。

圖4 咸陽(yáng)市地震災(zāi)害危險(xiǎn)性分區(qū)評(píng)價(jià)圖Fig.4 Evaluation map of earthquake disaster analysis in Xianyang city

圖5 咸陽(yáng)市多災(zāi)種自然災(zāi)害危險(xiǎn)性分區(qū)評(píng)價(jià)圖Fig.5 Evaluation map of multi natural disaster risk assessment in Xianyang city

4 多災(zāi)種自然災(zāi)害綜合風(fēng)險(xiǎn)評(píng)價(jià)

風(fēng)險(xiǎn)評(píng)價(jià)是在上述多災(zāi)種危險(xiǎn)性綜合分析的基礎(chǔ)上，對(duì)危害后果(人員、財(cái)產(chǎn)及易損性)進(jìn)行量化(表5)，將危險(xiǎn)性量化結(jié)果與危害后果量化結(jié)果(危害后果=人口數(shù)量或財(cái)產(chǎn)價(jià)值×易損率)進(jìn)行柵格乘積運(yùn)算，進(jìn)而得出多災(zāi)種自然災(zāi)害綜合風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果(圖6)。

圖6 咸陽(yáng)市多災(zāi)種自然災(zāi)害風(fēng)險(xiǎn)分區(qū)評(píng)價(jià)圖Fig.6 Evaluation map of multi-hazard natural disaster analysis in Xianyang city

表5 承災(zāi)體類型、價(jià)值及易損性估算值表Tab.5 Estimate values of hazard bearing body type,value and vulnerability

綜合分析認(rèn)為，咸陽(yáng)市自然災(zāi)害高風(fēng)險(xiǎn)區(qū)所占面積很小，零星分布于秦都區(qū)城區(qū)東南側(cè)、武功縣西北側(cè)、乾縣縣城區(qū)、三原縣城區(qū)東側(cè)、淳化縣城區(qū)、彬州市城區(qū)以及旬邑縣城區(qū)，總面積約65.57 km2，占全區(qū)面積的0.64%；自然災(zāi)害中風(fēng)險(xiǎn)區(qū)分布面積較廣，主要分布于武功縣城區(qū)至渭城區(qū)城區(qū)沿線地帶、永壽縣城區(qū)東側(cè)、三原縣城區(qū)外圍、淳化縣城區(qū)外圍、彬州市中部區(qū)域、旬邑縣城區(qū)外圍以及長(zhǎng)武縣南側(cè)，總面積約1 974.42 km2，占全區(qū)面積的19.40%；自然災(zāi)害低風(fēng)險(xiǎn)區(qū)分布面積最廣，遍布于各個(gè)區(qū)(縣)，總面積約6 576.48 km2，占全區(qū)面積的64.63%；自然災(zāi)害極低風(fēng)險(xiǎn)區(qū)分布面積較大，主要分布在彬州市、旬邑縣、禮泉縣、乾縣、武功縣城區(qū)至秦都區(qū)城區(qū)沿線北側(cè)區(qū)域以及渭河河谷區(qū)局部地段，總面積約1 558.42 km2，占全區(qū)面積的15.32%。

5 結(jié)論

(1)咸陽(yáng)市主要的災(zāi)害類型包括地質(zhì)災(zāi)害、洪水災(zāi)害和地震災(zāi)害3種。地質(zhì)災(zāi)害是影響咸陽(yáng)市綜合風(fēng)險(xiǎn)分區(qū)的最重要因素，其次為洪水災(zāi)害，最后是地震災(zāi)害。

(2)咸陽(yáng)市多災(zāi)種自然災(zāi)害風(fēng)險(xiǎn)分區(qū)呈現(xiàn)出明顯的空間差異性，整體上渭河河谷沿線地段和西北側(cè)的中低山丘陵區(qū)風(fēng)險(xiǎn)程度以中風(fēng)險(xiǎn)和高風(fēng)險(xiǎn)為主，主要集中在秦都區(qū)城區(qū)東南側(cè)、武功縣西北側(cè)、乾縣縣城區(qū)、三原縣城區(qū)東側(cè)、淳化縣城區(qū)、彬州市城區(qū)以及旬邑縣城區(qū)，其他區(qū)域?yàn)?zāi)害風(fēng)險(xiǎn)程度以低風(fēng)險(xiǎn)和極低風(fēng)險(xiǎn)為主。

(3)首次對(duì)陜西省咸陽(yáng)市開展了多災(zāi)種自然災(zāi)害綜合風(fēng)險(xiǎn)評(píng)價(jià)，評(píng)價(jià)結(jié)果可以為咸陽(yáng)市應(yīng)急管理和國(guó)土空間規(guī)劃帶來(lái)一定的借鑒意義和指導(dǎo)作用。