基于水文模型及水動力模型的山洪臨界雨量研究論文
引言
山洪是指山區(qū)溪溝中發(fā)生的暴漲洪水，影響面積小，具有突發(fā)、水量集中、流速大、沖刷破壞力強的特點。山洪及其誘發(fā)的泥石流、滑坡，常造成人員傷亡、設(shè)施毀壞，對人民生命財產(chǎn)造成嚴重危害。日趨嚴重的山洪災(zāi)害已得到了各級政府的高度重視，這對氣象部門山洪災(zāi)害風(fēng)險評估業(yè)務(wù)提出了迫切的需求。
山洪災(zāi)害風(fēng)險評估的關(guān)鍵有兩點，一是致災(zāi)臨界雨量的確定，另一個是風(fēng)險評估。致災(zāi)臨界雨量的分析計算與確定是山洪災(zāi)害研究的重要基礎(chǔ)。對于無資料或資料比較缺乏的地方，臨界雨量的分析方法主要采用內(nèi)插法、比擬法、山洪災(zāi)害事例調(diào)查法、災(zāi)害與降水頻率分析法等，這些方法多從雨量本身或災(zāi)情方面出發(fā)，沒有考慮水文過程。美國水文研究中心研發(fā)的山洪指導(dǎo)（flashfloodguidance,FFG）法基于動態(tài)臨界雨量的概念，通過水文模型分析計算，反推出流域出口洪峰流量達到預(yù)警流量所需的降水量，該方法兼顧了運行效率和水文過程，適用于有水文資料的流域。山洪災(zāi)害風(fēng)險評估的另一個關(guān)鍵點是風(fēng)險評估，目前的風(fēng)險評估多為定性的結(jié)果，或是基于歷史數(shù)據(jù)得出的“靜態(tài)”的評估結(jié)果，無法“動態(tài)”的進行風(fēng)險評估，利用GIS技術(shù)則可以解決這一問題，獲取洪水淹沒范圍，進而結(jié)合承災(zāi)體屬性實現(xiàn)山洪災(zāi)害的動態(tài)定量評估。
國外學(xué)者已嘗試將降水預(yù)報與水文模型、水力學(xué)模型結(jié)合應(yīng)用于洪水預(yù)報及早期預(yù)警、洪災(zāi)風(fēng)險評估中，國內(nèi)有學(xué)者亦嘗試了在較大的流域上建立從降水預(yù)報到洪水預(yù)報的模型。本文以閩江上游的寧化漁潭境內(nèi)的一個小流域作為研究對象，首先利用流域面雨量和流量率定該流域的TopModel模型參數(shù)，根據(jù)實地調(diào)查結(jié)果確定不同風(fēng)險等級的臨界流量，然后利用已率定的模型參數(shù)確定臨界流量對應(yīng)的臨界面雨量，同時將臨界水位輸入FloodArea淹沒模型得到淹沒水深和面積，并利用GIS平臺對不同隱患點進行風(fēng)險評估。
1方法簡介
1.1“TopModel”水文模型
1979年，Beven和Kirkby提出了以地形為基礎(chǔ)的TopModel半分布式流域水文模型（TopgraphybasedhydrologicalModel）。Topmodel是以地形為基礎(chǔ)的模型，基于DEM推求地形指數(shù)來反映流域水文響應(yīng)特性，模擬變動產(chǎn)流面積，即流域總徑流是飽和坡面流和壤中流之和。地形指數(shù)的大小和分布是影響模型產(chǎn)匯流計算的重要因子之一，模型中采用Quinn提出的多流向法，將網(wǎng)格單元的地形指數(shù)看作是隨機變量，經(jīng)過統(tǒng)計分析得到地形指數(shù)的頻率密度分布，地形指數(shù)相同的網(wǎng)格具有相同的水文響應(yīng)，用“地形指數(shù)—面積分布函數(shù)”來描述水文特性的空間不均勻性。Topmodel的匯流計算主要應(yīng)用坡面徑流滯時函數(shù)和河道演算函數(shù)，匯流計算主要采用單位線的方法，河道演算采用河道平均洪峰波速的方法，常采用簡單的常波速洪水演算方法。本文采用蘭開斯特大學(xué)環(huán)境與生物科學(xué)研究所開發(fā)的TopModel95.02版本，共有7個參數(shù)，分別是:Szm為土壤下滲率呈指數(shù)衰減的速率;T0為土壤剛達飽和時的傳導(dǎo)度;Td為重力排水的時間滯時參數(shù);Srmax為根帶最大蓄水能力;Sr0為根帶土壤的初始缺水量;RV為地表坡面匯流的有效速度;CHV為主河道匯流的有效速度。
1.2“FloodArea”淹沒模型
基于GIS柵格數(shù)據(jù)，德國Geomer公司開發(fā)了內(nèi)嵌于ArcGIS平臺的擴展模塊FloodArea，專門用于洪水演進模擬與動態(tài)風(fēng)險制圖。FloodArea采用ArcGRID數(shù)據(jù)格式，基于數(shù)字高程模型進行水文—水動力學(xué)建模，充分考慮了地形坡度和不同地表覆蓋形態(tài)下地面糙度對洪水演進形態(tài)的影響；洪水以給定水位、給定流量和給定面雨量三種方式進入模型，并可根據(jù)水文過程線進行實時調(diào)整，可視化表達流向、流速和淹沒水深等水文要素的時空物理場，為洪水淹沒風(fēng)險動態(tài)制圖提供了有效工具。
2典型個例分析
2.1山洪流域基本情況
本文選擇的流域位于福建省三明市寧化縣境內(nèi)，流域內(nèi)部地勢相對平坦，流域周圍被高山環(huán)繞，整個流域海拔為349～1100m。流域內(nèi)有水茜溪和東溪兩條河流，并在流域西南部出口處匯合。整個流域面積632km2，流域涉及到水茜、泉上、湖村、中沙4個鄉(xiāng)鎮(zhèn)的部分或全部區(qū)域。
2.2資料
2.2.1水位、流量資料
流域內(nèi)有一個漁潭水文觀測站，觀測逐日水位、流量，因資料有限，選擇2010年5月19～23日、2010年6月13～27日兩個過程的資料進行水文模型參數(shù)率定，選擇2012年5月13日的過程作為驗證數(shù)據(jù)。
2.2.2面雨量資料
2.3臨界雨量的計算
2.3.1“TopModel”模型的參數(shù)率定
通過2010年5~6月兩次暴雨過程的逐時面雨量和流量數(shù)據(jù)，根據(jù)目標函數(shù)調(diào)整模式的7個參數(shù)，使模擬流量盡可能接近輸入流量來確定參數(shù)大小。通過模式的參數(shù)率定，使兩次暴雨過程的確定性系數(shù)均達到0.94以上。
2.3.2致災(zāi)臨界雨量確定
根據(jù)實地調(diào)查，漁潭斷面壩高326m，大部分水田位置高出壩高1m。2010年6.18過程，漁潭發(fā)生了小洪水，漫溝1.75m并淹沒水田，調(diào)查得知這次洪水最高水位327.75m，在歷史上屬中等程度。根據(jù)調(diào)查結(jié)果，分別將漫溝、淹沒水田、淹沒居民點的水位作為不同風(fēng)險等級的`臨界水位，即326m為低風(fēng)險的臨界水位，327m和328m分別為中、高風(fēng)險的水位。
2.4山洪災(zāi)害風(fēng)險評估
取漁潭流域中落差相對較小且兩岸居民相對集中的主河道，對河道進行柵格化，河道高程代表河道基面高度，根據(jù)水位過程線的高度進行“河道漫頂式”模擬，得到漁潭流域下游河道兩岸不同風(fēng)險等級的淹沒情況。根據(jù)得到的淹沒情況圖，利用GIS的疊置分析功能，即可實現(xiàn)對不同風(fēng)險等級淹沒范圍內(nèi)居民點、道路等承載體影響程度的評估。
2.5效果檢驗
2012年前汛期期間，漁潭流域發(fā)生了一次低風(fēng)險山洪過程（圖5），發(fā)生低風(fēng)險預(yù)警開始時間是13日10時，此時水位正好達到低風(fēng)險水位326m，預(yù)警開始時間準確；12時水位上漲至預(yù)先設(shè)定的中風(fēng)險等級而面雨量仍然在低風(fēng)險級別，可見原先設(shè)定的中風(fēng)險水位偏低，因此在中風(fēng)險臨界雨量不變的情況下，上調(diào)中風(fēng)險水位至328m，高風(fēng)險水位相應(yīng)調(diào)整為329m，調(diào)整結(jié)果與實際情況是否符合，尚需進一步驗證。
3結(jié)論
本文用典型個例分析方法，以位于福建省三明市寧化縣境內(nèi)的漁潭流域為例研究具備完整水文觀測資料的山洪流域災(zāi)害風(fēng)險評估方法。具體步驟如下：首先利用雨量和流量資料率定流域的TopModel模型參數(shù)，把根據(jù)實地考察結(jié)果確定的致災(zāi)臨界流量和率定的模型參數(shù)回代模擬得到各級致災(zāi)臨界雨量；同時結(jié)合FloodArea淹沒模型，獲取不同臨界水位下的淹沒范圍，結(jié)合GIS技術(shù)確定淹沒范圍的影響程度，從而實現(xiàn)山洪災(zāi)害風(fēng)險的動態(tài)評估。
通過分析檢驗，得到以下初步結(jié)論：
（1）在具有完整水文觀測資料的流域，采用水文模型模擬，能夠很好地體現(xiàn)山洪災(zāi)害發(fā)生過程的雨洪關(guān)系，并得到不同風(fēng)險等級下的臨界雨量。這對于山洪災(zāi)害的防御有著重要作用，只要預(yù)報出流域內(nèi)降水量達到或超過風(fēng)險臨界雨量，就可以發(fā)布不同等級的山洪災(zāi)害預(yù)警。
（2）結(jié)合淹沒模型，可獲取流域不同風(fēng)險等級下可能淹沒范圍、淹沒水深，為可能發(fā)生山洪災(zāi)害區(qū)域進行人員遷移、災(zāi)害防御提供指導(dǎo)性建議。
（3）根據(jù)2012年前汛期強降水過程的檢驗，漁潭流域山洪事件得到準確的預(yù)警，但是尚需在今后的應(yīng)用中就臨界水位調(diào)整和提高預(yù)警時效做進一步研究。
【基于水文模型及水動力模型的山洪臨界雨量研究論文】相關(guān)文章：
人工定額的系統(tǒng)動力學(xué)模型研究論文11-05
基于蛛網(wǎng)模型的集團國際戰(zhàn)略化研究論文05-31
周邊績效模型研究論文12-01
基于Watterson模型的相關(guān)信道下Goldencode性能研究論文11-03
基于CMMI的軟件開發(fā)模型應(yīng)用研究論文11-04
基于模型的軟件測試技術(shù)探析論文05-03