基于多情景的巫山縣山洪災害危險評價

林孝松, 崔夢瑞, 牟鳳云, 余 情, 徐 州

(1.重慶交通大學 建筑與城市規劃學院, 重慶 400074; 2.重慶交通大學 土木工程學院, 重慶 400074)

山洪災害是山丘區經濟社會發展面臨的突出問題，每年均會造成重大人員傷亡和巨大財產損失。由于我國山丘區山洪孕災、致災因子復雜多變，空間異質性強，山洪形成機理不確定性顯著，人類活動不利影響加劇，山丘區建設與發展將長期面臨緊迫的山洪災害防災減災重任。當前國內外科技工作者在山洪災害孕災致災、風險評估以及災害預警等方面取得了較豐碩成果。Garambois等[1]基于面向過程水文模型對降雨觸發山洪進行了深入研究；杜俊等[2]分析基本因子和一般因子對不同類型山洪災害的影響程度；蔡維英等[3]基于分布式中小流域山洪災害模擬模型分析了各因素對山洪形成影響；Azmeri等[4]對印度尼西亞典型山區小流域山洪災害進行危險區識別研究；林孝松等[5]對跳石鎮山洪災害危險進行了綜合分區；Zhang等[6]對中國霞浦流域進行山洪風險制圖研究；黃國如等[7]對瑤安小流域山洪災害進行風險評估；岳琦等[8]完成了閩江上游山洪災害風險區劃；趙瑜等[9]對黑龍江省各區縣山洪災害風險進行等級區劃；劉少軍等[10]應用DEA模型對海南島山洪災害危險性大小進行排序，根據成災效率進行山洪災害危險性的等級劃分；葉金印等[11]提出采用降雨徑流相關法進行產流計算，采用Nash匯流模型進行匯流計算的山洪預警方法。當前以及將來全球氣候變化將使高強度降雨發生頻率、強度和區域分布變得更加復雜，所造成的山洪災害也將更為嚴重。

以重慶市巫山縣山洪災害為研究對象，針對山區縣域山洪災害形成特點及尺度差異，構建由微地貌類型、地形坡位、地形指數、土地利用綜合程度和地表匯流累積量等組成的山洪災害靜態、動態危險評價指標；在GIS軟件支持下，獲取研究區各危險評價靜態指標30 m分辨率數據，同時按照48 h降雨量為50，100，150，200，250 mm這5種強度值，考慮不同土地利用和土壤類型的下滲和產水特性，模擬獲取研究區地表匯流累積量；采用綜合指數法計算得到各網格單元山洪災害危險綜合指數并進行綜合分區；利用GIS區域匯總功能，將網格評價結果平均匯總到巫山縣各鄉鎮區和行政村，以利于各行政管理部門應急響應和實際應用。

1 研究區及數據

1.1 研究區概況

巫山縣位于重慶市東北部，地處三峽庫區腹心，位于東經109°33′—110°11′，北緯30°45′—23°28′，幅員面積2 958 km2，全縣轄25個鄉鎮和9個林、農場，2016年末常住人口45.55萬人，地區生產總值101.79億元。地形以深谷、中山和低山為主，海拔高度為75～2 690 m(圖1)。氣候屬亞熱帶季風性濕潤氣候，常年降水量1 041～1 600 mm，最大月降雨量達300 mm。2014年8月30日—9月2日，巫山縣普降暴雨，最大降雨量達384.6 mm，高強度降雨導致山洪暴發，95 616人受災，793間房屋倒塌，6 187 hm2農作物受災，36處道路中斷，直接經濟損失超過2.4億元。1978—2016年巫山縣歷史山洪災害點共770處(圖2)。

圖1 巫山縣數字高程模型

1.2 研究數據

研究所采用的數據主要包括巫山縣2015年土地利用類型圖、巫山縣1∶5萬30 m分辨率數字高程模型數據、巫山縣土壤類型圖、包含有鄉鎮(街道)和行政村(居委會)的2015年行政區劃圖、巫山縣歷史山洪災害點數據、2017年巫山縣統計年鑒等。

2 評價指標與數據獲取

2.1 評價指標選取

國內外相關研究者針對不同的研究尺度，選用不同的評價指標開展山洪災害危險評價。在大尺度方面，岳琦等[8]選用多年降水均值、高程、坡度、土壤類型、糙率、穩定下滲率、最長匯流路徑長度、最長匯流路徑比降、植被覆蓋度9個指標，進行了閩江上游山洪災害風險區劃；劉少軍等[10]選用地形因子、水系因子、暴雨日指數、滑坡和泥石流密度、綜合災度作為海南島山洪災害危險性評價指標；唐川等[12]采用地形坡度、暴雨天數、河網緩沖區、標準面積洪峰流量、泥石流分布密度和綜合災度6個指標，完成了紅河流域山洪災害危險評價圖；杜俊等[13]選用暴雨綜合指標、地形起伏度、土地利用產匯流能力、巖性軟硬程度4個指標對長江流域山洪災害危險進行評估。在小尺度方面，黃國如等[7]選用最大3 d降雨、高程、坡度、土壤類型等危險指標對瑤安小流域山洪災害進行風險評估；高展等[14]采用不同等級降水量及頻次、高程、高程標準差等指標對黃山風景區暴雨山洪災害進行風險區劃研究；左倩云等[15]采用年均降雨量、地形坡度、起伏度、植被覆蓋度、路網密度、居民點密度、歷史災害點密度、匯流累積量、暴雨強度、河網密度10個指標對接龍鎮山洪災害危險進行綜合評價。

圖2 巫山縣山洪災害點分布

綜合上述研究來看，基本上是從氣象、地形地貌、水文、人類活動、歷史災害情況等方面遴選相關指標進行山洪災害的危險評價，不足之處在于所有指標均是靜態的歷史性數據，未考慮降雨尤其是暴雨這一指標的動態性對山洪災害形成的影響。論文根據巫山縣山地丘陵地形特點以及研究區尺度差異特征，從地形地貌、水動力條件、人類活動強度等方面選用微地貌類型、地形坡位、地形指數、土地利用綜合程度等作為山洪危險評價的靜態指標，從不同降雨強度下洪水匯集差異角度選用地表匯流累積量作為山洪危險評價的動態指標。

2.2 微地貌類型與地形坡位

利用研究區30 m分辨率DEM數據，采用地形坡位指數(Topographic Position Index，TPI)法[16]可獲取巫山縣微地貌類型和地形坡位類型數據。首先，利用巫山縣30 m×30 m柵格DEM數據，采用90 m×90 m尺度矩形鄰域計算得到巫山縣TPI指數；然后，按照Weiss劃分標準，采用90 m×90 m和330 m×330 m尺度矩形鄰域提取到巫山縣10種微地貌類型(圖3)；最后，采用標準偏差(SD)方法對TPI值進行分等，并結合坡度值大小(表1)劃分得到巫山縣地形坡位類型(圖4)。

圖3 巫山縣微地貌類型分布局部示意圖

圖4 巫山縣地形坡位分布局部示意圖

地形坡位山谷谷地山坡下部平緩坡地陡峻坡地山坡上部山脊TPI標準偏差值≤-1(-1,-0.5](-0.5,0.5](-0.5,0.5](0.5,1]>1坡度——≤5°>5°——

2.3 地形指數

地形指數(Topographic Index，TI)是綜合考慮區域地形坡度和流域水文特性的一個綜合指標。其空間分布可用于確定土壤飽和缺水量的空間分布情況以及產流區空間位置和范圍。計算公式為：

TI=ln(α/tanβ)

(1)

式中：TI為地形指數；α為匯流面積；β為坡度。

地形指數TI值越大，則表征土壤越容易達到飽和而產流，從而有利于發生山洪災害。故可利用該值作為暴雨后地表產流并由此引發山洪的綜合評價指標。首先利用巫山縣DEM數據得到其坡度圖；然后利用GIS水文分析功能得到巫山縣各網格坡面匯流面積；最后基于公式(1)利用柵格計算器可得到巫山縣地形指數分布數據。

2.4 土地利用綜合程度

不同強度的人類活動直接或間接影響山洪災害的形成或強度大小。人類活動對地表形態造成一定的破壞作用，減弱對雨水的滲透和截留作用，土地利用程度不同其對地表破壞程度也會有所差異。由此可用土地利用綜合程度來表征人類活動強度大小以及其對山洪災害形成的影響。按照人類利用程度的大小差異，將城鎮村及工礦用地、交通用地等界定為高度利用等級，耕地、園地為中度利用等級，林地、草地、水域及水利設施用地為低度利用等級，其他用地為未利用等級。將巫山縣土地利用類型矢量數據按照利用等級標準進行柵格轉化即可得到研究區土地利用綜合程度柵格數據。

2.5 地表匯流累積量

結合巫山縣夏半年多年平均降雨情況，按照48 h降雨量50，100，150，200，250 mm這5種降雨強度模擬計算地表匯流累積量。首先基于研究區DEM數據，利用ArcGIS水文分析功能進行洼地填充、水流方向計算、匯流累積量(柵格數)模擬；然后根據巫山縣2015年各土地利用類型(圖5)和土壤類型(圖6)的下滲與產水特性，按照表2中系數賦予各土地利用類型和土壤類型相應的徑流系數值，并將兩個圖層進行相乘疊加得到巫山縣各柵格單元徑流系數；最后按照5種降雨強度值，基于以網格數量為單位的匯流累積量圖層，通過單位換算得到以m3為單位的地表匯流累積量，并以巫山縣徑流系數分布柵格數據為權重進行修正，最終得到不同降雨強度下巫山縣地表匯流累積量模擬分布數據。從模擬結果來看，隨著降雨量的增大，研究區各個網格單元地表匯流累積量均有不同程度增多。

圖5 巫山縣土地利用類型

圖6 巫山縣土壤類型

土地利用類型徑流系數值土地利用類型徑流系數值土地利用類型徑流系數值土壤類型徑流系數值建制鎮1果園0.60坑塘水面1紫色土0.88村莊0.99茶園0.65水工建筑用地1黃棕壤0.90采礦用地0.99其他園地0.70水庫水面1黃壤0.91風景名勝及特殊用地0.96有林地0.45河流水面1潮土0.92公路用地1灌木林地0.55裸地0.90石灰巖土0.94水田0.95其他林地0.60設施農用地0.95裸巖0.95灌溉水田0.90其他草地0.70沙地0.65水稻土0.96旱地0.75內陸灘涂0.95

3 結果與分析

巫山縣山洪災害危險評價5個指標中，微地貌類型、地形坡位、土地利用綜合程度為定性指標，地形指數和地表匯流累積量為定量指標。為進行綜合量化評價，將研究區山洪災害危險等級劃分為高、中、低危險和基本安全4個級別，并依據表3標準對各指標進行分等級標準化賦值處理。

表3 巫山縣山洪災害危險評價指標量化表

采用層次分析法對評價指標權重進行計算，得到微地貌類型、地形坡位、地形指數、土地利用綜合程度和地表匯流累積量的權重值依次為0.178 5，0.111 0，0.290 7，0.075 1，0.344 7。

基于評價指標標準量化柵格數據，利用柵格計算器功能計算得到不同降雨強度下巫山縣山洪災害危險指數綜合值(表4)，并按照閾值1.2，1.8，2.4進行等級劃分(圖7—8)；最后利用區域匯總功能可得到不同降雨強度下巫山縣各鄉鎮和各行政村的平均危險評價值，其中巫山縣48 h降雨量為250 mm時各鄉鎮和各行政村的危險評價結果如圖9—10所示。

表4 不同降雨強度下巫山縣山洪災害危險評價

由表4可知，巫山縣山洪災害以低危險和中危險為主，分別占總面積的43.4%和31.0%，而中危險和高危險面積占總面積的39%。在48 h內，降雨量由50 mm增加到250 mm的情況下，基本安全區面積基本不變，僅減少0.412 0 km2；低危險區面積稍有降低，降低了2.890 3 km2；中危險區面積變化較明顯，減少了10.425 1 km2；而高危險區面積增加較顯著，面積增加了13.727 5 km2。由此可預見，在降雨量不斷增加的情況下，巫山縣山洪災害高危險區的面積在不斷增加，防治任務加劇。

從巫山縣各鄉鎮山洪災害危險評價結果來看(圖9)，34個鄉鎮(林、農場)各危險等級面積占各自面積的比重分別為10.91%～27.87%(基本安全)，27.87%～52.64%(低危險)，25.47%～35.25%(中危險)和5.38%～10.32%(高危險)；中、高危險兩者總和比例超過40%的鄉鎮有鄧家、抱龍、金坪、龍溪、平河、竹賢、廟堂、當陽、巫峽和官陽10個鄉鎮。此外，按照自然斷點分類法將各鄉鎮危險綜合值分為4個等級，危險指數最高的第四等級鄉鎮主要集中在研究區的南部和西部，共11個；而危險等級較低的第一和第二等級主要為研究區內植被相對較好的林農場。

圖7 巫山縣山洪災害危險評價局部示意圖(降雨量50 mm)

從巫山縣各行政村山洪災害危險分布來看(圖10)，巫山縣316個行政村(居委會)中，危險指數值大于1.80的有23個，1.71～1.80的有148個，1.61～1.70的有140個，小于1.60的有5個，同時各行政村山洪災害危險評價結果與歷史災害分布情況吻合相對較好。

圖9 巫山縣各鄉鎮山洪災害危險評價

圖10 巫山縣各行政村山洪災害危險評價

4 討論與結論

(1) 在對各評價指標進行分級賦值量化的基礎上，采用多因素綜合指數法計算得到各評價單元山洪災害危險綜合指數值并據此將山洪災害危險劃分為高危險、中危險、低危險和基本安全4個等級；從總體來看，巫山縣山洪災害以低、中危險為主，占比分別為43.4%和31.0%，而中、高危險面積占比為39%。

(2) 在研究區48 h降雨量由50 mm增加到250 mm的情況下，基本安全區面積基本不變，低危險區面積稍有降低，中危險區面積變化較明顯，面積減少10.425 1 km2，而高危險區面積則增加較顯著，面積增加13.727 5 km2。由此可見，采用動態的地表匯流累積量既可進行動態評價，也可進一步開展不同降雨強度下的山洪災害空間預警研究。

(3) 從各鄉鎮山洪評價結果來看，中、高危險兩者總和比例超過40%的鄉鎮達10個，且危險指數值相對較高的鄉鎮主要集中在研究區的南部和西部；從各行政村山洪危險分布來看，危險指數值大于1.7的行政村占總數的54.11%，且評價結果與歷史災害分布情況吻合相對較好，可為研究區山洪災害防治提供借鑒。

論文研究成果可為巫山縣山洪災害防災減災提供科學依據，但也存在以下需要進一步研究之處。

首先，地表匯流累積量指標的賦值量化問題。將地表匯流累積量由通常的網格數或面積變為具體匯水量，可以模擬不同降雨條件下的地表匯流情況，并由此進行山洪災害危險的動態評價。但在對該指標進行分級賦值時是按照研究區河流整體特征確定的固定數額，未充分考慮河流上下游間河道寬度的差異，故按照此分級標準會使得對上游的評價結果偏輕，而對下游地區的評價會相應偏重，使得評價結果與實際情況存在一定偏差。

其次，地表匯流累積量的精確性問題。在進行不同降雨強度下地表匯流累積量的模擬時，充分考慮了不同地表下墊面尤其是不同土地利用類型和土壤類型對洪水形成的差異影響，但在確定各土地和土壤類型的徑流系數時主要借鑒專家經驗和少量試驗參數確定，其結果也會與實際情況存在一定誤差。

再次，山洪災害形成關鍵指標作用的均化問題。山洪災害的形成受各種因素的綜合影響，是一個復雜的綜合過程。雖然在進行危險評價時盡量選擇綜合性和關鍵性指標，但采用綜合指數方法通過線性組合計算山洪災害危險指數值，在一定程度上平滑了部分關鍵性指標所起的主要作用。