西藏東南部山洪災害過程水文動力模擬和臨界雨量

林志強， 尼瑪吉， 黃志誠

(1.西藏自治區氣候中心, 西藏 拉薩 850000; 2.西藏自治區氣象局信息網絡中心, 西藏 拉薩 850000)

西藏東南部山洪災害過程水文動力模擬和臨界雨量

林志強1， 尼瑪吉1， 黃志誠2

(1.西藏自治區氣候中心, 西藏 拉薩 850000; 2.西藏自治區氣象局信息網絡中心, 西藏 拉薩 850000)

[目的] 研究西藏山洪致災臨界雨量確定方法，為西藏地區水文氣象預警提供基礎數據資料，以服務于高原山洪災害防御工作。[方法] 采用水文動力模式Floodarea模型對西藏東南部的一次山洪過程進行模擬研究，并利用氣象資料逐時降水量進行了淹沒高度計算，得到了山洪致災臨界雨量。通過實地考察獲取西藏自治區林芝市巴宜區2015年8月一次山洪災害的基本參數資料，與模型模擬研究結果進行對比分析和模型驗證。[結果] (1) Floodarea對暴雨誘發山洪過程的模擬結果較好，結合多源降水數據能更準確地模擬西藏山區山洪暴發的過程； (2) 通過調整雨量情景設定可推算山洪不同淹沒水深的臨界雨量，得到較為準確的洪水淹沒范圍和降雨量—淹沒深度關系； (3) 應用降雨量—淹沒深度關系計算得到巴宜區永久河山洪溝的4個淹沒深度災害山洪等級(0.1，0.6，1.2，1.8 m)的6 h累計臨界雨量分別為33，55.7，75.4，91.9 mm。[結論] 研究結果表明Floodarea軟件適用于西藏高原水文觀測資料匱乏的復雜地形山區，能較準確地重現山洪災害過程和確定山洪臨界雨量，可為復雜地形山區山洪防治和預警提供參考。

Floodarea模型; 西藏林芝; 山洪; 臨界雨量

文獻參數： 林志強, 尼瑪吉, 黃志誠.西藏東南部山洪災害過程水文動力模擬和臨界雨量[J].水土保持通報，2017,37(1)：183-187.DOI:10.13961/j.cnki.stbctb.2017.01.033； Lin Zhiqiang, Nimaji, Huang Zhicheng, et al. Hydrological dynamics simulation and critical rainfall for flash flood in Southeastern Tibet[J]. Bulletin of Soil and Water Conservation, 2017,37(1)：183-187.DOI:10.13961/j.cnki.stbctb.2017.01.033

山洪是發生在山區溝谷或河道地形的突發性洪水，多為短歷時強降雨造成，具有持續時間短、能量集中、破壞性大，難以預防等特點，山洪及其誘發的泥石流和滑坡等次生災害常造成人員傷亡，沖毀房屋、田地、道路橋梁，甚至能導致水壩潰決，是國民經濟和人民生命財產安全的重大威脅[1-3]。西藏地區地處中國西南邊陲，是青藏高原的主體，向有“亞洲水塔”之美稱，中國的主要大河長江、黃河、怒江、金沙江都起源于西藏，雅魯藏布江、獅泉河、象泉河、瀾滄江等河流還流向南亞、東南亞國家，這些河流和淡水資源是下游地區生態和經濟的生命線，為數以億計的下游人民提供生存和發展所需的核心資源[4]，由于西藏地區地形復雜、地層巖性多變、構造運動頻繁，地形陡峻、切割破碎，生態環境脆弱，山洪災害是影響西藏經濟發展、基礎工程建設和人民生命財產安全的主要氣象災害之一[5]。在氣候變暖的背景下，西藏地區極端降水事件發生頻繁[6]，加之冰川退化，山洪災害的發生越來越頻繁，影響越來越嚴重。如何減輕山洪災害的影響，保障西藏生態安全屏障建設是西藏氣象工作者責無旁貸的重任。

致災臨界雨量閾值是山洪預報預警的關鍵指標之一[7]：當小流域范圍內的降雨量達到或超過一定量級，形成的山洪流量超出河道的安全泄洪能力，該雨量即為致災臨界面雨量。結合準確的區域定量雨量預報和致災臨界面雨量，氣象部門即可向社會和決策部門提供較為準確的山洪預警服務。基于精細DEM數據的洪水動態水文模擬[8-10]可以較為準確地分析山洪的發生和發展過程，模擬不同雨量量級的淹沒范圍和水深，是一種成熟的致災臨界雨量研究方案，已在中國不同區域得到了應用：張明達等[11]利用Floodarea模型分析了云南宣威2012年“7·12”山洪災害淹沒高度，并得到不同等級山洪災害臨界面雨量；姬興杰等[12]模擬了豫西山區無水文資料的暴雨誘發山洪過程，通過雨量情景設定方法確定了洛河上游地區4種雨量分配方式下的致災臨界雨量；張磊等[13]針對缺乏水文資料但有歷史洪水淹沒記錄的山區小流域進行山洪風險雨量計算，實現山洪災害精細化風險評估；謝五三等[14]利用Floodarea模型模擬了安徽大通河流域的暴雨淹沒情景，模擬值和實際的淹沒過程較為吻合；姜智懷等[15]分析了江西省曹水流域降水導致山洪過程的推進路線、淹沒范圍和淹沒水深，表明Floodarea模型嵌入河道柵格方法可用于推算曹水流域致災臨界雨量；文明章等[16]利用Floodarea再現了福建上清溪流域無水文資料山區的山洪過程，通過模擬結果水深和雨量相關關系推算流域受災隱患點不同等級淹沒水深的臨界雨量；葉麗梅等[17]依據蓄滿產流平衡原理和暴雨洪澇淹沒模型對通城縣一次強降水造成的洪澇淹沒水深和范圍進行了模擬。最近，劉義花等[18]對青海省羊智溝的洪水動態模擬表明Floodarea能適用于青藏高原這樣的復雜地形地區的山洪模擬，尤其是Floodarea模型可以適用于無水文觀測資料地區山洪分析，這尤其適用于地廣人稀、觀測資料匱乏的西藏地區，本文擬采用該模型通過對西藏東南部山洪災害過程分析，探討西藏地區山洪致災臨界雨量研究，為水文氣象預警提供基礎數據資料，為防御山洪災害提供參考。

1 災害過程概況

本研究選取的山洪溝位于西藏自治區林芝市巴宜區，山洪發生災害點位于94°21′38.69″E，緯度：29°36′07.61″N，海拔高度2 984 m。林芝地處西藏東南部，南接喜馬拉雅山脈，北靠念青唐古拉山脈，東連橫斷山脈，位于雅魯藏布江下游地區，主要有帕隆藏布和尼洋河，高山峽谷相間，多陡峻峽谷，侵蝕強烈，是地質災害的多發區[19]。林芝屬于高山溫濕季風氣候，印度洋暖濕氣流在夏季西南季風驅動下沿雅魯藏布江河谷北上[20]，形成較為豐沛的降水，植被生長茂盛，有“雪域江南”之譽，也是廣受游客喜愛的旅游目的地。林芝山洪災害多由強降水引發，多發于5—9月的汛期，山高坡陡，山洪匯流快、成災迅猛。由于林芝地處川藏交界，是入藏出藏重要交通干道，山洪及其次生地質災害往往造成交通中斷，甚至是嚴重的交通事故，是旅游事業發展的重大威脅。

2015年8月17—19日，在西太副高、孟加拉暖濕氣流輸送和西風槽東移的綜合影響下，西藏中東部地區大部持續強降水，其中林芝市多地有大到暴雨過程，其中林芝、波密和米林3個氣象站的日降水量打破了有氣象站記錄以來的日降水量極值，此次過程按照林志強等[6]的環流分型，屬于典型的西太副高西伸型。強降水造成318國道一鋼便橋被山洪沖毀，致交通中斷達20余天；山洪引發泥石流災害，導致地處永久河邊永久村的部分道路、橋梁、交通工具和房屋、店鋪等不同程度損毀，所幸未造成人員傷亡。此次山洪及其次生災害過程共造成直接經濟損失900多萬元，間接經濟損失達2 000多萬元。與其他山洪災害相比，由于災害發生地點距聚居區較近，有道路施工人員見證過程，為西藏東南部山區的山洪災害分析提供了寶貴的資料。

2 資料和方法

2.1 災害調查

由于研究區域地形起伏較大，不同的研究地點將導致淹沒深度變化，因此本文將以垮塌鋼便橋和永久河附近為致災臨界雨量的淹沒深度參考區，由于河道無防護堤，淹沒水深為基于河道測量的水深。由于山洪暴發過程時間短，災害強，親歷者并不是專業氣象、水文工作人員，因此很多信息并不完整，以山洪災害發生時受洪水淹沒的橋體設施和附近的樹木、施工簡易設施等淹沒痕跡進行實地測量，并進行相應的實地災害調查，準確記錄淹沒點的位置，包括經度、緯度、海拔高度、最大淹沒深度等信息。

附近村民和搶險施工人員和資料查詢，并通過多方信息相互映證，得到盡可能準確的淹沒過程信息，由于無法確定淹沒水深隨時間的變化，因此采取多點測量災害過程中淹沒跡線獲取最大淹沒深度的方法，與模擬時段得到的最大淹沒深度進行對比(表1)。

表1 山洪淹沒深度調查信息

2.2 氣象和地理信息數據

本文采用的降水數據為西藏自治區氣象局信息網絡中心提供的自動氣象站逐小時降水觀測數據，為彌補氣象觀測站不足，采用CMORPH(CPC MORPHing technique)多衛星觀測融合降水數據作為降水觀測的補充，以往的研究表明[10]，結合CMORPH數據和站點觀測能更準確地反映降水的強度和分布情況。

地形數據采用ASTER GDEM V2全球數字高程數據，空間分辨率達30 m，由日本METI和美國NASA聯合研制，于2015年1月6日免費面向公眾分發，數據由地理空間數據云(www.gscloud.cn)提供下載；土地利用數據采用歐洲空間局(ESA，http:∥due.esrin.esa.int/page_globcover.php)2009年全球陸地覆蓋數據集，分辨率為300 m，在ArcGIS平臺上通過最鄰近法插值到與DEM相一致的30 m。按照張洪江等[21]的試驗結果，由土地利用類型得到流域Manning水力糙度系數。

2.3 Floodarea模型

Floodarea是德國Geomer公司開發的水文動力學洪水淹沒模型，內嵌于ArcGIS平臺，被廣泛應用于洪水演進模擬、洪水動態風險評估等，國內的應用表明了其在不同地區山洪淹沒模擬的適用性[11-18]。Floodarea洪水過程同時考慮柵格周圍八個單元，水流量由Manning-Stricker公式計算，坡度由單元最低水位和最高高程差異決定，對每個單元都進行計算[22]，相鄰單元的水流寬度被認為是相等的；位于對角線的單元，以不同的長度算法來計算；水流方向柵格間坡度決定。Floodarea每個時相運行過程的相應淹沒范圍和水深都以柵格形式存儲，直觀地呈現洪水演進動態過程。

3 結果分析

3.1 山洪災害淹沒過程模擬

圖1給出了Floodarea模擬的2015年8月19日08時至20日08時(北京時，下同)暴雨過程鋼便橋附近的淹沒過程，同時給出了距山洪暴發區域最近的林芝自動氣象站的逐小時降水監測數據。降水過程具有顯著的日變化特征[23]，早晨降水較弱；隨著地面增暖加強，地面加熱形成對流，午后14時降水強度開始增大，到18時降水強度最大，達到10.1 mm/h；較強降水持續到夜間20時，隨后開始減弱；凌晨04時又開始形成較弱的降水。鋼便橋的高度約為2 m左右，根據調查鋼便橋約在19時40分垮塌，此時模擬深度達2.5 m左右，沒過鋼便橋約0.5 m，較降水峰值滯后2 h左右。模擬的淹沒水深較降水強度有所滯后，在20日0時達到最大，接近3 m，這較降水峰值滯后達5 h左右，較強降水時段滯后約4 h，這可能與流域集水過程有關，經過進一步檢查降水區域分布情況，發現在流域上游地區的自動站監測得到在20—23時上游有較強的降水，這綜合形成了淹沒深度的峰值；隨著降水的減弱，淹沒深度逐漸減小，到20日04時的降水后淹沒深度稍有回升。從山洪淹沒過程的模擬來看，Floodarea很好地反映了降水強度和分布的影響，與災情調查的結果吻合較好。

圖1 林芝巴宜流域永久河鋼便橋逐時淹沒水深模擬結果

3.2 災害調查與模擬對比

對山洪暴發區最大淹沒水深的調查結果與Floodarea模擬結果進行對比，圖2給出所有調查點的對比結果。模擬結果與實地調查得到的最大水深誤差在-0.2～2.9 m之間，其中5個調查點的模擬結果誤差小于0.2 m，8個調查點模擬結果誤差小于0.5 m；2個調查點的誤差大于1 m；大部分的模擬結果的淹沒深度大于調查結果，河床上的模擬結果較差，偏高幅度最大。

總體而言，Floodarea的模擬結果效果較好地反映了此次山洪過程的影響，淹沒范圍與實際發生山洪災害分布基本一致，準確地反映了地勢低地區的永久村、橋梁公路與其相鄰的渥羅幫嘎村的淹沒差異情況，表明在精細的地理信息基礎上，Floodarea可以較準確地反映洪水在地勢低洼區積水匯集和水流演進的過程。

3.3 降雨量—淹沒深度關系和臨界雨量

綜合考慮流域山洪災害影響的區域和模擬結果，地形、海拔高度、與河溝的距離遠近，居民聚居點、公路、橋梁等設施和社會影響等因素，選取鋼便橋、柏油公路和永河村3個地點為預警(山洪隱患)點，通過逐步增加面雨量的方式動態調整降水情景，分別模擬不同面雨量下的淹沒深度和淹沒水深；當隱患點的淹沒水深達到山洪災害等級時，即可得出預警點的山洪致災臨界雨量。

圖2 山洪災害調查點淹沒水深模擬與觀測對比

根據西藏東南部地區短時強降水的特征，以6 h為模擬時間長度，模擬降水時間步長為1 h，按照6 h累積雨量10～100 mm的降水強度分別進行模擬，以5 mm的增量逐漸調整降水強度，得到不同設定雨量強度下的淹沒模擬結果，并提取預警點的最大淹沒水深。圖3給出得到的3個預警點的淹沒水深和降雨量的關系，淹沒水深隨降雨量增加而增大，但兩者并不成線性關系，淹沒深度在降雨量較小時一直保持較低的淹沒深度，當6 h降雨量達到40～55 mm以上時，淹沒深度隨降雨量增加開始呈現線性增加趨勢；淹沒深度與降雨量的這種關系與日常業務中對山洪的認識是一致的，與山洪暴發的集水過程、地形和上游地面產流過程有關，這種降雨量—淹沒深度關系的非線性特征也是山洪多由短時強降水造成的原因。

圖3 林芝巴宜區永久河流域降雨量—淹沒深度關系

根據預警點得到的降雨量—淹沒水深關系，當預警點的淹沒水深分別達到0.1，0.6，1.2和1.8 m時，其對應的雨量即為該預警點的1,2,3和4級山洪災害致災臨界雨量(表2)。

3.4 臨界雨量山洪淹沒模擬

以地勢較低、受災影響較大的柏油公路得到的致災臨界雨量作為流域山洪致災臨界雨量，把1—4級山洪臨界雨量輸入Floodarea模型作為驅動雨量，計算得到不同山洪災害等級臨界雨量對應的淹沒范圍和深度。隨著臨界雨量的增大，洪水泛濫的范圍擴大，淹沒水深增大，在業務中，結合精細化分布的居民點和交通數據與淹沒柵格數據相疊加，能夠直觀地了解山洪災害的影響程度。以圖形和柵格數據的方式提供給流域地區氣象部門，根據氣象觀測數據和降水預報可得到流域山洪精細化預警。本研究還以相同的方法得到了流域4級山洪臨界雨量的1,3,12和24 h累積降雨量及其對應淹沒柵格數據，為不同時段山洪預警提供參考。

表2 林芝巴宜永久河流域山洪臨界雨量閾值

4 討論與結論

(1) Floodarea水動力模型可以較好地模擬西藏復雜地形地區的山洪淹沒過程，與災害普查的山洪影響范圍和淹沒深度較為一致，在無水文觀測資料的山區小流域山洪分析具有較高應用價值；

(2) 通過逐步增加面雨量值得方法，以模擬不同強度雨量情景，可用于確定山洪不同淹沒水深的雨量，建立降雨量—淹沒水深關系，進而推算不同山洪災害等級的致災臨界雨量；利用這個方法確定了3個山洪預警點的臨界雨量；

(3) 林芝巴宜區永久河山洪溝的0.1，0.6，1.2，1.8 m淹沒高度災害等級的6 h臨界雨量分別為33，55.7，75.4和91.9 mm。

從實地考察和衛星圖片上可以發現研究流域地區有多條泥石流溝，災情報告也強調了山洪誘發的泥石流災害，山洪和泥石流災害共同作用加大了山洪研究的難度；由于山洪的暴發過程較快，匯水和集水的過程時間短，Floodarea模型并未考慮水流的下滲問題，山洪臨界雨量計算時沒有考慮前期降雨和土壤含水狀況；在利用Floodarea進行淹沒模擬過程中發現，普查的GPS數據和海拔高度與DEM模型有一定的差異，這些因素都可能給淹沒模擬帶來一定的誤差，這是在服務中利用淹沒區域模擬結果需要加以注意的。在沒有山洪淹沒過程的準確資料情況下，采用多地點考察最大淹沒深度的方式進行驗證，可為西藏地區這樣幅員遼闊、人口稀少，山洪災害過程記錄匱乏的地區提供一種可行的災害普查方法和途徑；還需要再加強山洪災害風險普查工作，為西藏山洪分析和預警工作的開展提供更詳盡、準確的基礎資料，改善水文模型的率定效果和驗證；西藏地區的氣象觀測站點稀疏，山洪過程集水范圍較小，在業務應用中需要綜合多種資料得到盡可能精細的降水分布和面雨量以更準確評估山洪災害的影響。

[1] 劉少軍,張京紅,張明潔,等.海南島山洪災害風險區劃研究[J].水土保持研究,2013,20(5):165-169.

[2] 周忠發,李波,鄒長慧.望謨縣特大山洪泥石流災后重建分區與功能定位[J].水土保持通報,2013,33(1):174-177.

[3] 王凱博,高保玲,李洪濱,等.博州山洪災害防治對策探討[J].水土保持研究,2007,14(1):289-294.

[4] 鐘祥浩,劉淑珍,王小丹,等.西藏高原生態安全研究[J].山地學報,2010,28(1):1-10.

[5] 劉光軒.中國氣象災害大典：西藏卷[M].北京:氣象出版社,2008.

[6] 林志強,德慶,文勝軍,等.西藏高原汛期大到暴雨的時空分布和環流特征[J].暴雨災害,2014,33(1):73-79.

[7] 葉金印,李致家,劉靜,等.山洪災害氣象風險預警指標確定方法研究[J].暴雨災害,2016,35(1):25-30.

[8] 蘇布達,姜彤,郭業友,等.基于GIS柵格數據的洪水風險動態模擬模型及其應用[J].河海大學學報:自然科學版,2005,33(4):370-374.

[9] 葛小平,許有鵬,張琪,等. GIS支持下的洪水淹沒范圍模擬[J].水科學進展,2002,13(2):456-460.

[10] 林志強,洪健昌,尼瑪吉,等.基于HBV模型的尼洋曲流域洪水致災臨界面雨量研究[J].水土保持通報,2016,36(4):22-26.

[11] 張明達,李蒙,戴叢蕊,等.基于FloodArea模型的云南山洪淹沒模擬研究[J].災害學,2016,31(1):78-892.

[12] 姬興杰,李鳳秀,朱業玉,等.河南省洛河上游暴雨山洪致災臨界面雨量的確定[J].氣象與環境學報,2015,31(6):43-50.

[13] 張磊,王文,文明章,等.基于FloodArea模型的山洪災害精細化預警方法研究[J].復旦學報：自然科學版,2015,54(3):282-287.

[14] 謝五三,田紅,盧燕宇.基于FloodArea模型的大通河流域暴雨洪澇災害風險評估[J].暴雨災害,2015,34(4):384-387.

[15] 姜智懷,張毅之,蔡哲,等.基于河道柵格的山洪災害淹沒模擬[J].氣象,2014,40(8):1013-1018.

[16] 文明章,林昕,游立軍,等.山洪災害風險雨量評估方法研究[J].氣象,2013,39(10):1325-1330.

[17] 葉麗梅,周月華,李蘭,等.通城縣一次暴雨洪澇淹沒個例的模擬與檢驗[J].氣象,2013,39(6):699-703.

[18] 劉義花,魯延榮,周強,等.基于GIS柵格數據的青海省羊智溝洪水動態模擬[J].中國農業大學學報,2015,20(3):169-174.

[19] 胡桂勝,陳寧生,鄧虎.基于GIS的西藏林芝地區泥石流易發與危險區分析[J].水土保持研究,2012,19(3):195-199.

[20] 林志強,唐叔乙,何曉紅,等.西藏高原汛期水汽輸送特征與降水異常[J].氣象,2011,37(8):984-990.

[21] 張洪江,北原曜,謝明曙,等.晉西幾種不同土地利用狀況下糙率系數的研究[J].北京林業大學學報,1994,16(4):86-92.

[22] Geomer. Floodarea and Floodarea HPC: ArcGIS-extension for calculating flooded areas: User Manual V10.0[M/OL]. Heidelberg: Geomer GmbH,2011[2016-05-01]. http:∥www. geomer. de/fileadmin/templates/main/res/downloads/floodarea10_manual_englisch. pdf

[23] 林志強,假拉,羅骕翾,等.西藏高原閃電特性時空分布特征[J].氣象科技,2012,40(6):1002-1006.

Hydrological Dynamics Simulation and Critical Rainfall for Flash Flood in Southeastern Tibet

LIN Zhiqiang1, Nimaji1, HUANG Zhicheng2

(1.ClimateCenterofTibetAutonomousRegion,Lhasa850000,China; 2.MeteorologyInformationNetworkCenter,MeteorologicalBureauoftheTibetanAutonomousRegion,Lhasa850000,China)

[Objective] In this paper, we aimed to determine the critical rainfall of flash flood in Tibet, in order to provide the basic data for hydrological and meteorological warning in Tibet and serve the flash flood prevention. [Methods] We simulated the process of a flood in Southeastern Tibet by hydrodynamic model Floodarea, and calculated the submergence depth by using meteorological data of hourly precipitation, and we then obtained the critical rainfall of flash flood. To test the simulation result, we conducted a field survey to obtain the basic parameter of a flash flood disaster in August 2015 in Bayi Nyingchi region. [Results] The Floodarea model could be well used in the simulation of rainfall induced mountain flood process. The critical rainfall of flood submerged in different water depth could be calculated by adjusting the rainfall scenarios. According to the relationships between rainfall amount and submergence depth, the critical 6 hours accumulated rainfall amounts in four mountain flood ditches with different submergence depths(0.1, 0.5, 1.2, 1.8 m) were 33, 55.7, 75.4 and 91.9 mm, respectively. [Conclusion] Floodarea model software can be applied in the complex terrain mountainous areas that lack of hydrological observation data, in Tibetan Plateau. It can accurately simulate the flash flood disaster process and determine the critical rainfall. This study can provide reference for flood prevention and early warning in complex terrain mountainous areas.

Floodarea model; Nyingchi prefecture in Tibet; flash flood; critical rainfall

2016-05-05

2016-06-11

第三次青藏高原資助項目專題“西藏高原災害性天氣分析和預報方法研究”(GYHY201406001); 中國氣象局業務建設項目“西藏暴雨洪澇災害風險評估”

林志強(1982—),男(漢族),福建省漳州市人,碩士,工程師,主要從事西藏氣候分析工作。E-mail:linzq82@gmail.com。

A

1000-288X(2017)01-0183-05

P339