張曉蕾 ,劉榮華 ,劉 啟 ,李昌志 ,郭 良
(1.中國水利水電科學(xué)研究院,北京 100038; 2.水利部防洪抗旱減災(zāi)工程技術(shù)研究中心,北京 100038)
對于山洪災(zāi)害風(fēng)險分析,國內(nèi)外都進行了大量研究。從基于臨界雨量的統(tǒng)計分析[2],到基于指標(biāo)體系的定量分析,到基于遙感數(shù)據(jù)的空間分析,綜合天、空、地信息,計算方法趨于多元化。Wieczorek G F[3]應(yīng)用臨界雨量去預(yù)測山洪災(zāi)害風(fēng)險,并進行預(yù)警。S.Peruccacci等在1980年Nel Caine降雨時間和強度分析的基礎(chǔ)上,發(fā)展了臨界雨量方法,收集全球數(shù)據(jù)庫的2626場場次降雨數(shù)據(jù),進行了統(tǒng)計分析,并進行了區(qū)域劃分。譚徐明等[4]利用近300年的水災(zāi)資料及自然、社會經(jīng)濟基礎(chǔ)數(shù)據(jù),利用GIS技術(shù),完成了區(qū)域洪水風(fēng)險的分析。韓松等[5]構(gòu)建了太湖流域洪水災(zāi)害系統(tǒng)各因素之間的洪水風(fēng)險動因響應(yīng)關(guān)系,判斷了各因素的重要性和不確定程度。樊建勇等[6]研究了小時雨量與山洪災(zāi)害發(fā)生時間及小流域參數(shù)之間的關(guān)系,并運用2009年的場次山洪過程模擬,準(zhǔn)確預(yù)報了山洪風(fēng)險等級。美國水文研究中心HRC(2010)提出了動態(tài)山洪臨界雨量計算模型,并嵌入其國家山洪預(yù)警系統(tǒng)(FFGS)[7]。Totschnig[8]分析了阿爾卑斯山的山洪災(zāi)害風(fēng)險脆弱性,Holub[9]從建筑物的建筑結(jié)構(gòu)、風(fēng)暴降雨特征、冰雪覆蓋等方面,評估了歐洲中部山丘區(qū)建筑物的脆弱性。Guzzetti[10-11]基于1279—2002年的洪水及山洪災(zāi)害損失,分析了意大利北部山區(qū)的因災(zāi)死亡高發(fā)地,季節(jié)高發(fā)期的時空分布特征。Elkhrachy Ismail[12]基于SPOT和SRTM的DEM數(shù)據(jù),分析了沙特阿拉伯地區(qū)Najran市的山洪風(fēng)險,構(gòu)建了徑流、土壤類型、坡度、糙率、河網(wǎng)密度和距主河道距離等指標(biāo)體系,評估了地區(qū)山洪風(fēng)險。
目前,國內(nèi)外對于山洪災(zāi)害風(fēng)險的研究多集中在大尺度,比如山洪災(zāi)害區(qū)劃,或者以柵格為單元進行評估,以自然小流域單元的評估方法較少,而自然流域單元的產(chǎn)匯流是山洪形成的基本條件,為此,基于山丘區(qū)地形地貌和下墊面條件等小流域?qū)傩赃M行山洪風(fēng)險評估是具有一定物理機理的評估方法。
為此,本文基于全國山洪災(zāi)害防治項目形成小流域基礎(chǔ)數(shù)據(jù)集和產(chǎn)匯流特征參數(shù)集,結(jié)合山洪災(zāi)害調(diào)查評價成果,以福建省為例,探討了基于小流域單元的山洪災(zāi)害風(fēng)險因子選取、風(fēng)險評估指標(biāo)體系構(gòu)建和風(fēng)險評估模型建立等問題,并通過歷史山洪事件的對比分析,檢驗了方法的有效性和結(jié)果的合理性。
福建省位于我國的東南沿海,境內(nèi)峰嶺聳峙,丘陵連綿,山丘區(qū)面積占90%,遠高于全國平均水平(66.7%)。復(fù)雜的地形地質(zhì)條件、暴雨多發(fā)的氣候特征、密集的人口分布和人類活動的影響,導(dǎo)致山洪災(zāi)害發(fā)生頻繁,造成巨大的經(jīng)濟損失和大量人員傷亡。2010—2016年,福建省85個區(qū)縣中的74個區(qū)縣實施了全國山洪災(zāi)害防治縣級非工程措施項目建設(shè),涉及68574個自然村(圖1),人口2347萬人,防治區(qū)總面積達9.96萬km2,約占福建省陸地面積的81.55%。其中受山洪威脅縣城74個,受山洪威脅鄉(xiāng)鎮(zhèn)949個,受山洪威脅嚴(yán)重的沿河村落10 922個,受山洪威脅的企事業(yè)單位5478個。
根據(jù)全國山洪災(zāi)害調(diào)查評價成果[13],福建省建國后共發(fā)生歷史山洪災(zāi)害2306場(圖2)。2016年福建泰寧“5·8”大型泥石流災(zāi)害造成1座4層辦公樓、5處工棚推倒掩埋,39名施工人員失聯(lián)。
圖1 福建省山洪災(zāi)害防治村分布
圖2 福建省歷史山洪災(zāi)害分布(1949—2015年)
3.1 風(fēng)險識別與指標(biāo)體系構(gòu)建根據(jù)福建省降雨、下墊面和人類活動的特點,選擇山丘區(qū)降雨特性、小流域特性參數(shù)及下墊面特征、社會經(jīng)濟情況、現(xiàn)狀防洪能力評價結(jié)果、涉水工程及監(jiān)測預(yù)警設(shè)施設(shè)備調(diào)查成果,對指標(biāo)進行歸類,并按照小流域進行數(shù)據(jù)整理和指標(biāo)篩選,得到備選指標(biāo)共54個。主要包括:
(1)降雨指標(biāo)8個:各年型設(shè)計暴雨絕對值(5年一遇、10年一遇、20年一遇、50年一遇、100年一遇)、設(shè)計暴雨變差系數(shù)(1 h、3 h、6 h);
(2)下墊面指標(biāo)9個:土地利用和植被類型、土壤類型、土壤質(zhì)地類型、洪峰模數(shù)(10 min30 mm、10 min35 mm、60 min50 mm)、匯流時間(10 min30 mm、10 min35 mm、60 min50 mm);
(3)社會經(jīng)濟指標(biāo)23個:小流域人口、村落數(shù)、財產(chǎn)類型(Ⅰ類財產(chǎn)、Ⅱ類財產(chǎn)、Ⅲ類財產(chǎn)、Ⅳ類財產(chǎn))、房屋類型統(tǒng)計及占比(Ⅰ類房屋、Ⅱ類房屋、Ⅲ類房屋、Ⅳ類房屋、堅固房屋比例、脆弱房屋比例)、企事業(yè)單位數(shù)、小流域內(nèi)距河網(wǎng)的高程(3 m、5 m、10 m、30 m)及距離因子(50 m、100 m)組合10組;
(4)現(xiàn)狀防洪能力指標(biāo)7個:小于5年一遇的村落數(shù)、5~10年一遇的村落數(shù)、10~20年一遇村落數(shù)、20~50年一遇村落數(shù)、50~100年一遇村落數(shù)、大于100年一遇村落數(shù)、臨界雨量與設(shè)計暴雨比值小于1的村落數(shù);
(5)涉水工程指標(biāo)3個:小流域內(nèi)橋梁數(shù)、路涵數(shù)、塘(堰)壩數(shù);
(6)監(jiān)測預(yù)警配置率指標(biāo)4個:小流域內(nèi)自動監(jiān)測站數(shù)量、簡易雨量站、水位站及無線預(yù)警廣播站數(shù)量。
運用SPSS工具采用主成分分析法對54個備選指標(biāo)進行了降維處理,篩選出核心且獨立的指標(biāo)用于風(fēng)險評估,得到了10個主成分(表1),成分主要集中在自然村落與河網(wǎng)位置關(guān)系,小流域短歷時強降雨特性、匯流時間、洪峰模數(shù)、涉水工程、房屋類型、現(xiàn)狀防洪能力、監(jiān)測預(yù)警設(shè)施設(shè)備等方面。
根據(jù)風(fēng)險分析理論(ADRC,2005)[14],將以上指標(biāo)整合構(gòu)建了山洪災(zāi)害風(fēng)險分析指標(biāo)體系(表1),從危險性、暴露性和易損性三個維度進行指標(biāo)體系構(gòu)建和分析。
風(fēng)險(risk)=危險性(Hazard)×暴露性(Exposure)×易損性(Vulnerability)
危險性:(1)降雨因子:短歷時強降雨是山洪事件的激發(fā)因子,是山洪危險性的初始來源,主成分分析表明,歷時為3 h和6 h的稀遇降雨量及變率是非常突出的主成分,故選擇100年一遇3 h、6 h設(shè)計暴雨值及變差系數(shù)CV作為降雨因子;(2)地形因子:流域形狀、坡度、植被覆蓋、土壤類型及分布、溝道比降等因素,都會對產(chǎn)匯流造成影響,是山洪危險性的主要來源,洪峰模數(shù)及匯流時間是綜合考慮了以上因素分析計算后得到的結(jié)果,在流域尺度可反映下墊面特征及產(chǎn)匯流特性,故選擇這兩項指標(biāo)綜合為危險性的地形因子。
暴露性:山洪災(zāi)害調(diào)查評價對受山洪威脅的沿河村落進行了詳查和分析評估,主成分分析結(jié)果表明了村落與山洪溝的空間關(guān)系對山洪事件影響突出,故選擇受山洪威脅的沿河村落的人口、房屋及居民家庭財產(chǎn)類型為承險體暴露性的具體指標(biāo)。
易損性:表現(xiàn)為在風(fēng)險事件中容易遭受損害的特性及其程度(也可以理解為承險體所具備防災(zāi)力的相對面)。根據(jù)山洪災(zāi)害調(diào)查評價內(nèi)容,主成分分析結(jié)果表明,沿河村落房屋結(jié)構(gòu)及各種監(jiān)測預(yù)警設(shè)施的信息可以大致直接反映易損性的基本情況,故選擇Ⅲ類和Ⅳ類房屋占比例、自動監(jiān)測站和簡易監(jiān)測站的站點密度作為表征易損性的指標(biāo)。
表1 山洪風(fēng)險評估指標(biāo)體系及主成分分析結(jié)果
3.2 權(quán)重設(shè)定權(quán)重設(shè)定時,危險性、暴露性和易損性都是風(fēng)險要素的組成,同等重要,三者權(quán)重均為1/3。在確定二級指標(biāo)權(quán)重時,危險性指標(biāo)中短歷時強降雨特征明顯,降雨賦予權(quán)重較高;暴露性指標(biāo)由小流域人口數(shù)、房屋類型和居民家庭財產(chǎn)類型組成,人口數(shù)權(quán)重相對較高;易損性指標(biāo)包括Ⅲ和Ⅳ類房屋比例、自動監(jiān)測站和簡易監(jiān)測站的密度,考慮到監(jiān)測站點對與應(yīng)急避險轉(zhuǎn)移的重要性,權(quán)重相對較高。各指標(biāo)權(quán)重如表2。
3.3 風(fēng)險等級計算對數(shù)據(jù)進行歸一化處理后,構(gòu)建風(fēng)險立方體模型,將危險性、暴露性、易損性要素疊積。計算方法如式(1):
其中,危險性、暴露性、易損性三要素計算方法如式(2)—式(4):
表2 指標(biāo)權(quán)重設(shè)置表
式中:Hi、Ej、Vk為各指標(biāo)項指標(biāo);m、n、l為指標(biāo)項數(shù)量;m′、n′、l′為指標(biāo)項的各要素數(shù)量;i、j、k、k′為累計求和中間變量;W為要素和指標(biāo)的權(quán)重。
對危險性、暴露性、易損性三要素樣本量的計算值進行降維排序,以每類要素的1/3、2/3處設(shè)為劃分節(jié)點,將三要素計算值分別劃分為低、中、高3個等級風(fēng)險,綜合風(fēng)險需由三個要素組合疊加,即27個等級組合,如圖3所示,具體分類方法如表3:
圖3 風(fēng)險等級閾值
表3 風(fēng)險等級閾值組合
根據(jù)本文提出的風(fēng)險評估模型對福建省山洪災(zāi)害綜合風(fēng)險進行了評價,結(jié)果如圖4所示。圖中紅色為高風(fēng)險區(qū),黃色為中風(fēng)險區(qū),藍色為低風(fēng)險區(qū)。高風(fēng)險區(qū)主要集中在福建沿海臺風(fēng)高發(fā)區(qū)及內(nèi)陸暴雨高值山區(qū),如武夷山脈、玳瑁山、博平嶺等,涉及小流域1151個,流域面積近2萬km2,人口608萬人,占比為19.9%;中風(fēng)險區(qū)在省內(nèi)分布較為廣泛,以人口密集區(qū)為主,涉及小流域3280個,流域面積4.7萬km2,人口1227萬人,占比為40.2%;低風(fēng)險區(qū)涉及小流域4480個,流域面積5.7萬km2,人口1215萬人,占比為39.8%。高、中風(fēng)險區(qū)分布與山洪災(zāi)害防治區(qū)范圍基本一致。圖5是福建省典型山區(qū),即閩江千嶺水庫以上流域的山洪災(zāi)害綜合風(fēng)險分布圖。
圖4 福建省山洪災(zāi)害綜合風(fēng)險分布
圖5 閩江千嶺水庫以上流域山洪防治區(qū)綜合風(fēng)險圖
根據(jù)福建省山洪災(zāi)害調(diào)查評價成果中的歷史山洪災(zāi)害資料計算小流域山洪災(zāi)害事件發(fā)生密度,并與風(fēng)險評估結(jié)果對照,結(jié)果顯示,80%的歷史山洪落在高風(fēng)險區(qū)和中風(fēng)險區(qū),高風(fēng)險區(qū)密度為221場/萬km2,中風(fēng)險區(qū)密度為202場/萬km2,低風(fēng)險區(qū)密度為66場/萬km2。高風(fēng)險區(qū)的山洪災(zāi)害發(fā)生密度是低風(fēng)險區(qū)的3倍,驗證結(jié)果說明基于小流域的山洪災(zāi)害風(fēng)險評估模型指標(biāo)選擇具有較好的代表性,權(quán)重設(shè)定合理,風(fēng)險早期識別結(jié)果可靠。
本文構(gòu)建了山洪災(zāi)害早期識別與定量分析模型,構(gòu)建了指標(biāo)體系和風(fēng)險立方體模型,并以福建省為例,分析了全省的山洪災(zāi)害綜合風(fēng)險,并運用2306歷史山洪場次數(shù)據(jù)對評估結(jié)果進行驗證,結(jié)論如下:(1)基于小流域的山洪風(fēng)險識別和災(zāi)害風(fēng)險分析,可反映山洪災(zāi)害對短時強降雨、流域下墊面與人類活動的響應(yīng)關(guān)系,以自然流域為單元的風(fēng)險評估方法具有一定的物理機制。(2)評估結(jié)果表明福建省山洪災(zāi)害高發(fā)區(qū)集中在沿海地區(qū)受臺風(fēng)暴雨影響的區(qū)域和內(nèi)陸山區(qū)的強對流暴雨高值區(qū),山洪災(zāi)害發(fā)生頻度受自然環(huán)境和人類活動雙重因素影響。評估結(jié)果統(tǒng)計得到高發(fā)區(qū)人口為608萬,面積近2萬km2。(3)通過歷史山洪場次數(shù)據(jù)驗證,80%的歷史山洪落在高風(fēng)險區(qū)和中風(fēng)險區(qū),高風(fēng)險區(qū)的山洪災(zāi)害發(fā)生密度是低風(fēng)險區(qū)的3倍,風(fēng)險識別結(jié)果可靠,風(fēng)險評估模型可較為準(zhǔn)確的反映山洪災(zāi)害風(fēng)險。