李春燕，孟 暉，張若琳，李亞民

(中國地質(zhì)環(huán)境監(jiān)測院,北京 100081)

中國縣域單元地質(zhì)災(zāi)害風(fēng)險評估

李春燕，孟 暉，張若琳，李亞民

(中國地質(zhì)環(huán)境監(jiān)測院,北京 100081)

文中地質(zhì)災(zāi)害包括崩塌、滑坡、泥石流、地面塌陷四類突發(fā)性地質(zhì)災(zāi)害。采用地質(zhì)災(zāi)害危險性與承災(zāi)體易損性分項測算，定性綜合評估的方式實現(xiàn)我國縣域單元地質(zhì)災(zāi)害風(fēng)險評估。將承災(zāi)體易損性分為人口安全易損性與資產(chǎn)易損性，用因災(zāi)死亡人口比與因災(zāi)直接經(jīng)濟(jì)損失比表征?；谝詺v史數(shù)據(jù)分析指導(dǎo)未來預(yù)測的思想，以國土資源部2001—2015年地質(zhì)災(zāi)害災(zāi)情數(shù)據(jù)為樣本，提出了人口安全易損性與資產(chǎn)易損性的分級標(biāo)準(zhǔn)。對我國2869個縣域單元進(jìn)行地質(zhì)災(zāi)害風(fēng)險評價，結(jié)果表明共有216個地質(zhì)災(zāi)害高風(fēng)險單元、643個中風(fēng)險單元和2010個低風(fēng)險單元。地質(zhì)災(zāi)害高風(fēng)險區(qū)集中分布于烏蒙山區(qū)、四川盆地周邊山區(qū)以及云南西部和西藏東南部。

縣域單元；地質(zhì)災(zāi)害；易損性；風(fēng)險評估

我國是世界上地質(zhì)災(zāi)害災(zāi)情最為嚴(yán)重的國家之一。開展縣域單元地質(zhì)災(zāi)害風(fēng)險評估，有助于從宏觀層面對地質(zhì)災(zāi)害形勢進(jìn)行研判，減輕地質(zhì)災(zāi)害造成的損失，進(jìn)而實現(xiàn)地質(zhì)災(zāi)害風(fēng)險管理，服務(wù)于防災(zāi)減災(zāi)和生態(tài)文明建設(shè)。

對地質(zhì)災(zāi)害風(fēng)險評估的關(guān)注與研究于近幾十年興起。1989—2005年，聯(lián)合國召開數(shù)次會議明確防災(zāi)減災(zāi)的必要性與總體思路，研討地質(zhì)災(zāi)害風(fēng)險評估技術(shù)方法，使該領(lǐng)域逐步成為熱點(diǎn)。瑞士、法國等國相繼開展了地質(zhì)災(zāi)害風(fēng)險評估工作[1～5]。在國內(nèi)，張梁等[6]構(gòu)造了地質(zhì)災(zāi)害風(fēng)險評估指標(biāo)體系與模型；朱良峰等[7]基于GIS平臺研發(fā)了地質(zhì)災(zāi)害風(fēng)險評估系統(tǒng)，進(jìn)行了全國地質(zhì)災(zāi)害風(fēng)險試評價；張春山等[8～9]先后對黃河上游地區(qū)崩塌、滑坡、泥石流災(zāi)害危險性與風(fēng)險性進(jìn)行了評價；殷坤龍等[10]從地質(zhì)災(zāi)害預(yù)警的角度出發(fā)，探索了地質(zhì)災(zāi)害風(fēng)險預(yù)測評估方法；張茂省等[11]構(gòu)建了風(fēng)險矩陣評估方法，并在延安寶塔區(qū)進(jìn)行了實踐；吳樹仁等[12～13]重點(diǎn)探索了地質(zhì)災(zāi)害易發(fā)程度、危險性和風(fēng)險評估的工作流程，初步提出定性分析-定量化評價相結(jié)合的地質(zhì)災(zāi)害風(fēng)險評估技術(shù)方法，并根據(jù)張梁、殷坤龍等的全國人口地質(zhì)災(zāi)害傷亡率分級研究成果[14]，基于1995—2005年全國地質(zhì)災(zāi)害數(shù)據(jù)，定量提出我國地質(zhì)災(zāi)害人口傷亡率的最高等級閾值。

以往的研究成果豐富，評價方法多采用多元統(tǒng)計分析法、定性權(quán)重打分法、GIS系統(tǒng)疊加分析等方法，神經(jīng)網(wǎng)絡(luò)法也得到了初步應(yīng)用[15～19]；研究區(qū)則以小區(qū)域，大比例尺居多；評價結(jié)果則均為空間化結(jié)果。受資料來源的可靠性和積累程度的限制，以往的研究難以精確反映風(fēng)險的區(qū)域差異。不論何種風(fēng)險性評估方法，易損性定量評價都是地質(zhì)災(zāi)害風(fēng)險評估的難點(diǎn)之一，目前尚無統(tǒng)一的指標(biāo)和方法。

本文基于國土資源部地質(zhì)災(zāi)害報告制度所取得的長時間序列災(zāi)情數(shù)據(jù)，以概率統(tǒng)計法為基礎(chǔ)，探討以縣域為單元易損性量化評價方法及其分級標(biāo)準(zhǔn)，在地質(zhì)災(zāi)害危險性和易損性分項測算的基礎(chǔ)上綜合分析，從而實現(xiàn)區(qū)域地質(zhì)災(zāi)害風(fēng)險性評估。通過災(zāi)害調(diào)查數(shù)據(jù)和災(zāi)情統(tǒng)計數(shù)據(jù)的不斷更新，可進(jìn)一步實現(xiàn)風(fēng)險動態(tài)評估。

1 數(shù)據(jù)來源與算法

1.1 數(shù)據(jù)來源

地質(zhì)災(zāi)害調(diào)查數(shù)據(jù)來源于全國地質(zhì)災(zāi)害調(diào)查數(shù)據(jù)庫；基礎(chǔ)地質(zhì)條件數(shù)據(jù)主要包括地形起伏度圖、地貌圖、地質(zhì)圖、第四紀(jì)地質(zhì)圖、活動斷裂圖、工程地質(zhì)圖、水文地質(zhì)圖、可溶巖類型圖、土壤侵蝕圖等；地質(zhì)災(zāi)害點(diǎn)災(zāi)情數(shù)據(jù)來源于國土資源部2001—2015年地質(zhì)災(zāi)害速報制度和月報制度報送的數(shù)據(jù)；人口與經(jīng)濟(jì)數(shù)據(jù)來源于2001—2015年社會經(jīng)濟(jì)統(tǒng)計年鑒；縣級行政區(qū)劃數(shù)據(jù)采用國家地理信息中心制作的2013年縣級行政區(qū)劃數(shù)據(jù)。

1.2 算法

1.2.1 地質(zhì)災(zāi)害危險性評價

地質(zhì)災(zāi)害危險性評價時，為保障結(jié)果的合理性，應(yīng)首先進(jìn)行自然單元評價，再將評價結(jié)果轉(zhuǎn)化為縣域單元。

地質(zhì)災(zāi)害危險性的評價方法比較多，常見的方法有模糊綜合評價法、證據(jù)權(quán)法、多元回歸分析法、人工神經(jīng)網(wǎng)絡(luò)法、層次分析法、灰色關(guān)聯(lián)分析法等等。這些方法基于的理論不同、含義不同，對評價數(shù)據(jù)的需求不同，各有優(yōu)點(diǎn)和不足。本文選擇證據(jù)權(quán)法進(jìn)行地質(zhì)災(zāi)害危險性評價。

證據(jù)權(quán)法是一種由點(diǎn)及面的定量化評價方法，早期主要應(yīng)用于礦產(chǎn)資源預(yù)測等相關(guān)領(lǐng)域，近年來才被引入地質(zhì)災(zāi)害危險性評價中[16]，具有易于理解、操作簡便、受主觀因素影響小等特點(diǎn)。通過將已經(jīng)發(fā)災(zāi)的滑坡、崩塌、泥石流、地面塌陷災(zāi)害點(diǎn)與地形坡度、暴雨次數(shù)等評價指標(biāo)空間疊加，進(jìn)行相關(guān)性分析，求得每個指標(biāo)與滑坡、崩塌、泥石流、地面塌陷的相關(guān)程度，基于相關(guān)程度賦予各指標(biāo)不同的權(quán)重，進(jìn)行疊加計算，得到崩塌、滑坡、泥石流危險性指數(shù)圖和地面塌陷危險性指數(shù)圖，根據(jù)危險性指數(shù)由高到低將危險性劃分為高、中、低、極低4個等級。將崩塌、滑坡、泥石流危險性評價分級圖和地面塌陷危險性評價分級圖按照“就高原則”疊加，形成自然單元地質(zhì)災(zāi)害危險性評價結(jié)果。

自然單元評價結(jié)果向縣域單元轉(zhuǎn)化時采用最大面積法，即以縣域面積內(nèi)面積最大的評價等級作為該縣域的危險性評價等級。再結(jié)合已有資料、根據(jù)專家經(jīng)驗對評價結(jié)果進(jìn)行校驗修訂，最終確定縣域地質(zhì)災(zāi)害危險性等級。

1.2.2 承災(zāi)體易損性評價

承災(zāi)體易損性是指地質(zhì)災(zāi)害發(fā)生時可能對承災(zāi)體造成的破壞、損失程度。地質(zhì)災(zāi)害造成的損失主要有直接和間接經(jīng)濟(jì)損失、人員傷亡和資源環(huán)境破壞等。由于人口、財產(chǎn)和環(huán)境的性質(zhì)不同，且考慮到災(zāi)情統(tǒng)計數(shù)據(jù)的內(nèi)容，將承災(zāi)體易損性分為人口安全易損性與資產(chǎn)易損性兩部分。

人口安全易損性是指因地質(zhì)災(zāi)害造成人口死亡的最大可能程度。鑒于本次評價是以縣域為單元的宏觀尺度評價，可采用縣年均地質(zhì)災(zāi)害最大可能死亡人口占縣年均總?cè)丝诘谋龋匆驗?zāi)死亡人口比來表達(dá)。因災(zāi)死亡人口比越大，人口安全易損程度越高。

資產(chǎn)易損性是指縣域內(nèi)地質(zhì)災(zāi)害造成直接經(jīng)濟(jì)損失的最大可能程度。采用縣年均地質(zhì)災(zāi)害直接經(jīng)濟(jì)損失占縣年均經(jīng)濟(jì)總量的比，即因災(zāi)直接經(jīng)濟(jì)損失比來表達(dá)。其比值越大，資產(chǎn)易損程度越高。

由于縣年均地質(zhì)災(zāi)害最大可能死亡人口與可能直接經(jīng)濟(jì)損失難以預(yù)測，基于以歷史數(shù)據(jù)分析指導(dǎo)未來預(yù)測的思想，本文采用地質(zhì)災(zāi)害實際造成的縣年均地質(zhì)災(zāi)害死亡人口與年均地質(zhì)災(zāi)害直接經(jīng)濟(jì)損失代替計算。

以GIS空間分析統(tǒng)計功能為基礎(chǔ)，根據(jù)國土資源部2001—2015年地質(zhì)災(zāi)害災(zāi)情統(tǒng)計數(shù)據(jù)，以縣為單元統(tǒng)計2001—2015年縣域內(nèi)所有地質(zhì)災(zāi)害點(diǎn)共造成的死亡人口總數(shù)與直接經(jīng)濟(jì)損失數(shù)，計算年均值；2001—2015年期間，僅在2010年開展過一次人口普查，數(shù)據(jù)確實。因此，在實際計算因災(zāi)死亡人口比時，采用2010年人口普查的縣域總?cè)丝跀?shù)據(jù)作為縣年均總?cè)丝跀?shù)。

通過以全國縣域“因災(zāi)死亡人口比”、“因災(zāi)直接經(jīng)濟(jì)損失比”數(shù)據(jù)為樣本進(jìn)行測算分級，結(jié)合20世紀(jì)地質(zhì)災(zāi)害災(zāi)情數(shù)據(jù)和前人的研究成果[13～14]，提出人口安全易損性評價分級標(biāo)準(zhǔn)(表1)、資產(chǎn)易損性評價分級標(biāo)準(zhǔn)(表2)。

表1 人口安全易損性評價分級

表2 資產(chǎn)易損性評價分級

將縣域地質(zhì)災(zāi)害人口安全易損性評價分區(qū)結(jié)果與資產(chǎn)易損性評價分區(qū)結(jié)果，按照“就高原則”進(jìn)行疊加，形成縣域單元承災(zāi)體易損性評價分區(qū)結(jié)果。

1.2.3 地質(zhì)災(zāi)害風(fēng)險定性綜合評估

按照地質(zhì)災(zāi)害風(fēng)險定性綜合評估方法，根據(jù)地質(zhì)災(zāi)害危險性4個等級、承災(zāi)體易損性4個等級的組合特征，將地質(zhì)災(zāi)害風(fēng)險劃分為高風(fēng)險、中風(fēng)險、低風(fēng)險3個等級(表3)。為了保證地質(zhì)災(zāi)害風(fēng)險性評估結(jié)果的可靠性，應(yīng)通過統(tǒng)計分析、野外實地驗證或同行檢查對評估結(jié)果進(jìn)行必要的檢驗和完善。

表3 地質(zhì)災(zāi)害風(fēng)險定性綜合評估

2 全國地質(zhì)災(zāi)害風(fēng)險評估

2.1 評價單元

本次評價的基本單元包括市轄區(qū)、縣、自治縣、縣級市。根據(jù)采用的2013年行政區(qū)劃，全國陸域共有2869個評價單元(不包括臺灣省和南海諸島)。

2.2 評價與結(jié)果分析

2.2.1 地質(zhì)災(zāi)害危險性評價

地質(zhì)災(zāi)害危險性評價選取地形起伏度、地貌類型、工程地質(zhì)巖組、活動斷裂距離、多年平均降水量、年平均暴雨日數(shù)(天)、土地利用程度、礦山分布密度、崩塌滑坡泥石流發(fā)育密度、地面塌陷發(fā)育密度等指標(biāo)。將各評價指標(biāo)空間分布圖進(jìn)行1 km×1 km柵格化處理，按照評價流程形成自然單元地質(zhì)災(zāi)害危險性評價分區(qū)結(jié)果。將自然單元評價結(jié)果轉(zhuǎn)化，以縣域單元的形式表達(dá)，得到全國縣域單元地質(zhì)災(zāi)害危險性評價分區(qū)圖(圖1)。

圖1 地質(zhì)災(zāi)害危險性評價分區(qū)圖Fig.1 Zoning map of geological hazard assessment

據(jù)圖1，我國共有地質(zhì)災(zāi)害高危險單元453個、中危險單元1106個、低危險單元746個、極低危險單元564個。地質(zhì)災(zāi)害危險性高的區(qū)域主要連片分布于云南、四川、重慶、貴州大部以及西藏西南部、陜西南部、甘肅東部的山地丘陵區(qū)，遼寧、河北、江西等省有零星分布。這些地區(qū)地形切割劇烈，構(gòu)造發(fā)育強(qiáng)烈，與我國地質(zhì)環(huán)境背景條件相對應(yīng)。

2.2.2 承災(zāi)體易損性評價

2001—2015年地質(zhì)災(zāi)害共在798個縣(市、區(qū))造成了人口死亡，因災(zāi)死亡人口總數(shù)為8 100人(圖2、表4)。人口安全高易損單元25個、中易損單元324個、低易損單元422個、極低易損單元2 098個。人口安全高易損單元零星分布于云南、四川、甘肅、陜西等省，與地質(zhì)災(zāi)害危險性并不完全重合，說明危險性與實際的人口傷亡并沒有高度的相關(guān)性。同時，可以發(fā)現(xiàn)人口安全高易損單元基本分布在人口密度比較低的區(qū)域，說明人口集中的地區(qū)地質(zhì)環(huán)境相對安全或者防災(zāi)減災(zāi)措施到位，發(fā)災(zāi)較少，即使發(fā)災(zāi)也不易造成人口死亡；同時，人口安全高易損單元的人口基數(shù)低，導(dǎo)致這些單元的因災(zāi)死亡人口比偏高。甘肅省舟曲縣由于2010年發(fā)生特大型泥石流的緣故，死亡人數(shù)是所有評價單元中最高的，共計1 451人。因災(zāi)死亡人口比同樣是舟曲縣最高，為7.322 29。

圖2 地質(zhì)災(zāi)害人口安全易損性評價分區(qū)圖Fig.2 Zoning map of population safety vulnerability assessment

省(市)名市(州)名縣(區(qū))名縣(區(qū))人口總數(shù)(萬人)2001—2015年因災(zāi)死亡人口總數(shù)(人)因災(zāi)死亡人口比(人/萬人)人口安全易損性等級甘肅省甘南藏族自治州舟曲縣13.2114517.32229高云南省怒江傈僳族自治州貢山獨(dú)龍族怒族自治縣3.79891.56577高云南省玉溪市新平彝族傣族自治縣28.53680.15887高西藏自治區(qū)拉薩市墨竹工卡縣4.47660.98491高四川省甘孜藏族自治州丹巴縣5.97650.72590高四川省綿陽市北川羌族自治縣19.71630.21308高重慶市—武隆縣35.10620.11775高陜西省榆林市子洲縣17.40540.20691高四川省甘孜藏族自治州康定縣13.01530.27150高四川省雅安市石棉縣12.36430.23193高四川省阿壩藏族羌族自治州汶川縣10.08360.23816高貴州省黔東南苗族侗族自治州三穗縣15.57360.15411高浙江省麗水市慶元縣14.15300.14130高福建省三明市建寧縣12.00300.16670高陜西省安康市嵐皋縣15.42280.12109高山西省呂梁市中陽縣14.14260.12261高云南省怒江傈僳族自治州福貢縣9.86260.17577高遼寧省本溪市南芬區(qū)8.29190.15276高陜西省延安市甘泉縣7.72190.16410高四川省甘孜藏族自治州九龍縣6.21140.15022高新疆自治區(qū)昌吉回族自治州木壘哈薩克自治縣6.57140.14202高四川省雅安市寶興縣5.61130.15460高四川省阿壩藏族羌族自治州小金縣7.77130.11150高甘肅省張掖市肅南裕固族自治縣3.37130.25753高陜西省安康市鎮(zhèn)坪縣5.1090.11773高

地質(zhì)災(zāi)害共在1685個縣(市、區(qū))造成了直接經(jīng)濟(jì)損失(圖3、表5)，資產(chǎn)高易損單元26個、中易損單元136個、低易損單元395個、極低易損單元2 412個。較于人口安全高易損單元的分散分布，資產(chǎn)高易損單元分布相對集中于成都-蘭州-西安區(qū)域，除受地形因素影響外，還受到龍門山斷裂、鮮水河斷裂、祁連山北緣斷裂的控制。數(shù)據(jù)特征與人口安全易損性類似，資產(chǎn)高易損單元分布區(qū)域GDP密度較低。直接經(jīng)濟(jì)損失最高的縣為遼寧省新賓滿族自治縣。據(jù)統(tǒng)計數(shù)據(jù)，2001—2015年該縣因地質(zhì)災(zāi)害共造成直接經(jīng)濟(jì)損失1 002 377萬元，由于縣GDP總量為760 189.44萬元，故因災(zāi)直接經(jīng)濟(jì)損失比為8.790 59，居全國第三；西藏自治區(qū)洛扎縣因災(zāi)直接經(jīng)濟(jì)損失7 863.48萬元，遠(yuǎn)低于新賓滿族自治縣，但洛扎縣年均GDP總量僅2 281.72萬元，所以因災(zāi)直接經(jīng)濟(jì)比達(dá)到22.98，為全國最高。

將人口安全易損性與資產(chǎn)易損性按照“就高原則”疊加，形成全國縣域單元地質(zhì)災(zāi)害承災(zāi)體易損性評價分區(qū)圖(圖4)。承災(zāi)體高易損單元49個、中易損單元380個、低易損單元541個、極低易損單元1 899個。

2.2.3 地質(zhì)災(zāi)害風(fēng)險定性綜合評估

采用地質(zhì)災(zāi)害風(fēng)險定性綜合評估方法(表4)，將地質(zhì)災(zāi)害危險性與承災(zāi)體易損性集成分析，形成全國縣域單元地質(zhì)災(zāi)害風(fēng)險評估結(jié)果(圖5)。全國共有216個地質(zhì)災(zāi)害高風(fēng)險單元、643個中風(fēng)險單元和2 010個低風(fēng)險單元。地質(zhì)災(zāi)害高風(fēng)險區(qū)集中分布于烏蒙山區(qū)、四川盆地周邊山地區(qū)、甘肅東部以及云南西部和西藏東南部。與我國地質(zhì)環(huán)境背景條件比較吻合。

圖3 地質(zhì)災(zāi)害資產(chǎn)易損性評價分區(qū)圖Fig.3 Zoning map of population property vulnerability assessment

省(市)名市(州)名縣(區(qū))名縣(區(qū))年均GDP總量(萬元)2001—2015年因災(zāi)直接經(jīng)濟(jì)損失數(shù)(萬元)因災(zāi)直接經(jīng)濟(jì)損失比(萬/百萬)資產(chǎn)易損性等級西藏自治區(qū)山南地區(qū)洛扎縣2281.727863.4822.98000高甘肅省隴南市文縣132340.62178502.408.99207高遼寧省撫順市新賓滿族自治縣760189.4410023778.79059高云南省怒江傈僳族自治州貢山獨(dú)龍族怒族自治縣38631.0935686.106.15844高云南省怒江傈僳族自治州福貢縣33998.7131054.066.08926高西藏自治區(qū)昌都地區(qū)八宿縣6082.084089.204.48224高甘肅省隴南市武都區(qū)479098.07237245.863.30128高四川省阿壩藏族羌族自治州理縣81004.2634890.532.87150高四川省阿壩藏族羌族自治州松潘縣81842.0433998.332.76943高四川省甘孜藏族自治州丹巴縣64755.1426804.022.75953高西藏自治區(qū)日喀則地區(qū)亞東縣23214.318171.402.34666高西藏自治區(qū)林芝地區(qū)察隅縣29214.0910143.652.31478高四川省阿壩藏族羌族自治州馬爾康縣99106.9831169.702.09670高四川省阿壩藏族羌族自治州小金縣62945.9219045.262.01710高西藏自治區(qū)昌都地區(qū)邊壩縣10370.203051.701.96184高西藏自治區(qū)昌都地區(qū)左貢縣11410.803277.701.91497高四川省甘孜藏族自治州瀘定縣82196.9422531.101.82741高重慶市-武隆縣594881.631438421.61200高甘肅省甘南藏族自治州夏河縣86285.51208401.61016高四川省阿壩藏族羌族自治州茂縣139977.9533568.151.59874高甘肅省天水市麥積區(qū)839998.741996771.58474高甘肅省天水市清水縣191355.6244706.701.55754高四川省巴中市南江縣542636.31116532.391.43168高四川省阿壩藏族羌族自治州汶川縣301320.9961563.591.36208高甘肅省隴南市康縣94068.0018888.601.33865高甘肅省隴南市宕昌縣97215.44146861.00711高

圖4 地質(zhì)災(zāi)害承災(zāi)體易損性評價分區(qū)圖Fig.4 Zoning map of hazard-affected bodies separately vulnerability assessment

圖5 地質(zhì)災(zāi)害風(fēng)險程度評價分區(qū)圖Fig.5 Zoning map of geological hazard risk assessment

3 結(jié)論與建議

(1)本文將承災(zāi)體易損性分為人口安全易損性和資產(chǎn)易損性，分別用因災(zāi)死亡人口和因災(zāi)直接經(jīng)濟(jì)損失比量化表達(dá)。以國土資源部權(quán)威數(shù)據(jù)為基礎(chǔ)，通過測算，分別給出了二者分級標(biāo)準(zhǔn)以及疊加原則，再結(jié)合地質(zhì)災(zāi)害危險性進(jìn)行定性綜合評估，形成以縣為單元的區(qū)域地質(zhì)災(zāi)害風(fēng)險評估方法。此方法簡便快捷、數(shù)據(jù)連續(xù)有保障。

(2)對于地質(zhì)災(zāi)害高風(fēng)險縣(區(qū))，在規(guī)劃和土地利用前，應(yīng)考慮地質(zhì)災(zāi)害對土地利用適宜性的限制，進(jìn)行詳細(xì)的地質(zhì)災(zāi)害危險性評估，重點(diǎn)區(qū)需詳細(xì)勘查地質(zhì)條件，注重事前降低地質(zhì)災(zāi)害風(fēng)險，使區(qū)域經(jīng)濟(jì)與地質(zhì)環(huán)境可持續(xù)發(fā)展。需強(qiáng)調(diào)，地質(zhì)災(zāi)害低風(fēng)險區(qū)是基于歷史災(zāi)情判斷其地質(zhì)災(zāi)害風(fēng)險相對低，并不表示絕對的安全，在開展社會經(jīng)濟(jì)活動前也應(yīng)按需求進(jìn)行適度的調(diào)查，必要時采取一定的防災(zāi)措施。

(3)對于承災(zāi)體易損性的評價，是以歷史災(zāi)情數(shù)據(jù)作為分析樣本，基于歷史數(shù)據(jù)對未來預(yù)測具有指導(dǎo)意義的思想。但如何據(jù)此實現(xiàn)更加精細(xì)化的預(yù)測，尚需更深入研究。通過數(shù)據(jù)的不斷更新來實現(xiàn)評價的動態(tài)化或許是一種可行方法。

(4)以縣域單元進(jìn)行地質(zhì)災(zāi)害風(fēng)險評估，其評價結(jié)果是每個縣域為同一個風(fēng)險等級，默認(rèn)地質(zhì)災(zāi)害危險性、承災(zāi)體易損性等在該評價單元內(nèi)是均勻分布的。但實際上，縣域面積變化很大，其內(nèi)地形地貌不同，地質(zhì)災(zāi)害、人口、財產(chǎn)、土地等并非均勻分布。因此，本評估方法適用于具有較長時間災(zāi)情記錄的宏觀空間尺度的風(fēng)險分析。如何解決縣域內(nèi)部易損性差異的問題，提升地質(zhì)災(zāi)害風(fēng)險評估的精確性是需要進(jìn)一步研究的。

[1] Cascini L, Bonnard C, Corominas J,etal. Landslide hazard and risk zoning for urban planning and development[C].Balkema:CRC Press Inc,2005.

[2] Hungr O, Fell R, Couture R,etal. Landslide Risk Management: Proceedings of the International Conference on Landslide Risk Management[C]. Vancouver: CRC Press Inc,2005.

[3] Fell R, Corominas J, Bonnard C,etal. Guidelines for landslide susceptibility, hazard and risk zoning for land use planning[J]. Engineering Geology, 2008, 102(3/4):85-98.

[4] Leroi E, Bonnard C, Fell R,etal. Risk assessment and management[M]// Weinheim: wiley-vchverlaggmbh & Co. kgaa, 2005:693-694.

[5] Jibson R W, Harp E L, Michael J A. A method for producing digital probabilistic seismic landslide hazard maps[J]. Engineering Geology, 2000, 58(3/4):271-289.

[6] 張梁，張建軍.地質(zhì)災(zāi)害風(fēng)險區(qū)劃理論與方法[J].地質(zhì)災(zāi)害與環(huán)境保護(hù)，2000，11(4):323-328.[ZHANG L,ZHANG J J. The theory and method of risk zonation of geo-hazard[J]. Journal of Geologiccal Hazard and Environment Preservation,2000,11(4):323-328.(in Chinese)]

[7] 朱良峰，殷坤龍，張梁，等.GIS支持下的地質(zhì)災(zāi)害風(fēng)險分析[J].長江科學(xué)院院報，2002(5)：42-45.[ZHU L F,YIN K L, ZHANG L,etal. Risk analysis system of geo-hazard supported by GIS[J]. Joural of Yangze River Scientific Research Institute, 2002(5): 42-45.(in Chinese)]

[8] 張春山，張業(yè)成，馬寅生.黃河上游地區(qū)崩塌、滑坡、泥石流地質(zhì)災(zāi)害區(qū)域危險性評價[J].地質(zhì)力學(xué)學(xué)報，2003，9(2):143-153.[ZHANG C S, ZHANG Y C, MA Y S. Regional dangerous on the geological hazards of collapse, landslide and debris flows in the upper reaches of the Yellow River[J]. Journal of Geomechanics, 2003，9(2): 143-153.(in Chinese)]

[9] 張春山,張業(yè)成,馬寅生,等. 區(qū)域地質(zhì)災(zāi)害風(fēng)險評價要素權(quán)值計算方法及應(yīng)用——以黃河上游地區(qū)地質(zhì)災(zāi)害風(fēng)險評價為例[J]. 水文地質(zhì)工程地質(zhì),2006, 33(6):84-88. [ZHANG C S, ZHANG Y C, MA Y S,etal. Calculation method and application of the right-weighty value on geological hazard in region[J]. Hydrogeology & Engineering Geology,2006,33(6):84-88.(in Chinese)]

[10] 殷坤龍，陳麗霞，張桂榮.區(qū)域滑坡災(zāi)害預(yù)測預(yù)警與風(fēng)險評價[J].地學(xué)前緣，2007，14(6):85-97.[YIN K L, CHEN L X, ZHANG G R. Regional landslide hazard warning and risk assessment[J].Earth Science Frontiers,2007, 14(6):85-97.(in Chinese)]

[11] 張茂省，唐亞明.地質(zhì)災(zāi)害風(fēng)險調(diào)查的方法與實踐[J].地質(zhì)通報，2008，08:1205-1216.[ZHANG M S, TANG Y M. Risk investigation method and practice of geohazards[J]. Geological Bulletin of China, 2008, 27(8):1205-1216.(in Chinese)]

[12] 吳樹仁,石菊松,張春山,等. 地質(zhì)災(zāi)害風(fēng)險評估技術(shù)指南初論[J]. 地質(zhì)通報,2009,28(8):995-1005.[WU S R, SHI J S, ZHANG C S,etal. Preliminary discussion on technical guideline for geohazard risk assessment[J]. GeologicalBulletin of China, 2009, 28(8):995-1005.(in Chinese)]

[13] 吳樹仁，石菊松，張春山，等.滑坡風(fēng)險評估理論與技術(shù)[M]. 北京: 科學(xué)出版社，2012.[WU S R, SHI J S, ZHANG C S,etal. Theory and Technique of Landslide Risk Assessment[M]. Beijing: Science Press, 2012.(in Chinese)]

[14] 張梁，殷坤龍，雷明堂，等.全國地質(zhì)災(zāi)害風(fēng)險區(qū)劃研究報告[R].廊坊：中國國土資源經(jīng)濟(jì)研究院，2002.[ZHANGL, YIN K L, LEI M T,etal. Research Report on the risk zoning of geological hazards of China[R]. Langfang: Chinese Academy of Land&Resource Economic, 2002(in Chinese)]

[15] 向喜瓊，黃潤秋.基于GIS的人工神經(jīng)網(wǎng)絡(luò)模型在地質(zhì)災(zāi)害危險性區(qū)劃中的應(yīng)用[J].中國地質(zhì)災(zāi)害與防治學(xué)報，2000，11(3):26-30.[XIANG X Q, HUANG R Q. Application of GIS-based artificial Neural Networks on assessment of geo-hazards risk[J]. The Chinese Journal of Geological Hazard and Control, 2000,11(3):26-30.(in Chinese)]

[16] 范強(qiáng),巨能攀,向喜瓊,等. 證據(jù)權(quán)法在區(qū)域滑坡危險性評價中的應(yīng)用——以貴州省為例[J]. 工程地質(zhì)學(xué)報,2014,22(3):474-481. [FAN Q, JU N P, XIANG X Q,etal. Landslides hazards assessment with weights of evidence—a case study in Guizhou, China[J]. Journal of Engineering Geology, 2014,22(3):474-481.(in Chinese)]

[17] 程凌鵬,楊冰,劉傳正. 區(qū)域地質(zhì)災(zāi)害風(fēng)險評價研究述評[J]. 水文地質(zhì)工程地質(zhì),2001,28(3):75-78.[CHENG L P, YANG B, LIU C Z. Review on risk assessment of regional geological hazards[J]. Hydrogeology & Engineering Geology, 2001,28(3):75-78.(in Chinese)]

[18] 張超,陳艷,張宇飛,等.基于多元線性回歸模型的云南昭通地質(zhì)災(zāi)害易發(fā)性評價[J].水文地質(zhì)工程地質(zhì),2016,43(3):159-163. [ZHANG C, CHEN Y, ZHANG Y F,etal. Geohazard susceptibility evaluation in Zhaotong on Yunnan based on the multivariate linear regression model[J].Hydrogeology & Engineering Geology,2016,43(4):159-163. (in Chinese)]

[19] 石愛軍,馬娟,齊安文,等.物聯(lián)網(wǎng)技術(shù)在突發(fā)地質(zhì)災(zāi)害應(yīng)急響應(yīng)中的應(yīng)用研究[J].水文地質(zhì)工程地質(zhì),2014,41(5):148-152. [SHI A J, MA J, QI A W,etal. Research on application of the Internet of things(IOT) platform of geological disaster emergency response[J].Hydrogeology & Engineering Geology, 2014,41(5):148-152. (in Chinese)]

Risk assessment of geo-hazard of China in county unit

LI Chunyan, MENG Hui, ZHANG Ruolin, LI Yamin

(ChinaInstituteofGeo-environmentalMonitoring,Beijing100081)

In this paper, geo-hazards include landslide, collapse, debris flow and ground collapse. By calculating the risk of geo-hazard and vulnerability of hazard-affected bodies separately, we evaluated the geo-hazard risk of county territory unit by qualitative comprehensive assessment. First, we separated the vulnerability of hazard-affected bodies into population safety vulnerability and property vulnerability, then we use the rate of dead population and rate of direct economic loss that caused by geo-hazard as representation as index. Based on historical data analysis guiding the future prediction, combing with the database of geo-hazard situation during 2001—2005 from Ministry of Land and Resources, we proposed the grading standard of population safety vulnerability and property vulnerability. Based on the geological hazard risk evaluation from the 2869 county territory units, 216 units are in the high-risk, 643 units are in middle-risk and 2010 are in low-risk. The high-risk areas are dominated in Wumeng mountainous, Mountainous around Sichuan Basin, Western Yunnan and Southeast part of Tibet.

county unit; geo-hazard; vulnerability; risk assessment

10.16030/j.cnki.issn.1000-3665.2017.02.24

2016-10-08;

2017-01-11

中國地質(zhì)調(diào)查項目“地質(zhì)環(huán)境承載能力評價與監(jiān)測示范研究”(12120115050701)

李春燕(1987-)，男，碩士，工程師，主要從事地質(zhì)環(huán)境綜合研究工作。E-mail：licy@mail.cigem.gov.cn

孟暉(1962-)，女，本科，教授級高級工程師，主要從事環(huán)境地質(zhì)研究。E-mail：mengh@mail.cigem.gov.cn

P642.2

A

1000-3665(2017)02-0160-07