老撾無資料地區小流域設計洪水計算

賈春強　劉尚

【摘 要】老撾大多數小流域水文實測資料匱乏，無法滿足工程設計洪水計算。針對這一問題，本文以老撾Nam Pae河小流域設計洪水計算為例，采用STRM90m的DEM資料，提出了基于Arcgis軟件和林平一公式的小流域設計洪水計算方法，有效解決了林平一公式在計算小流域設計洪水時參數難以獲取的問題，大大提高了該公式的實用性，使其能夠廣泛應用于老撾無資料小流域地區的設計洪水計算中。

【關鍵詞】無資料地區；設計洪水；林平一公式；DEM；Arcgis

1 引言

老撾水資源豐富，其中流域面積在100km2以下的河流更是分布廣泛，小河流之上修建小型水庫、灌溉工程等對地區工業及農業都具有重要意義，而水利建筑物的設計需要依據小流域某種頻率的設計洪水[1]。老撾社會經濟發展較慢，科學技術水平不高，小河流上通常不設置水文站點，無法提供設計所需水文數據。林平一法計算小流域設計洪水在無資料地區廣泛應用，但由于公式中匯流參數較多且不易獲取，給水文工作者帶來較大困難，隨著DEM 數據的不斷完善，應用Arcgis 軟件計算的流域水文參數已日趨精確。鑒于此，本文以老撾Nam Pae河小流域設計洪水計算為例，提出了基于Arcgis軟件和林平一公式的小流域設計洪水計算方法，以期為老撾無資料地區設計洪水計算提供一種行之有效的解決途徑。

2 計算原理

林平一法在進行流域匯流的研究中，把流域的組成假定為具有一條坡度均勻的理想單干型的水道，其兩側為整齊而對稱的坡面，即設想為矩形的模型，該法主要適用于150km2以下的小流域洪水計算。用此種成因法理論來推算流域匯流過程的核心是幾個匯流參數的確定[2]。DEM數據通過有序數值陣列形式描述地面實體高程，能夠定量反映地貌結構和水文過程等空間變化特征[3]，可利用Arcgis軟件對DEM數據進行分析，計算林平一公式中的匯流參數。

2.1 林平一法

2.2 ArcGIS水文分析

基于DEM數據的Arcgis水文分析的主要內容是提取地表水流徑流模型的水流方向、匯流累積量、河流網絡以及對研究區的流域進行分割等。

2.2.1 無洼地DEM生成

DEM 是比較光滑的地形表面模型，但由于DEM 誤差以及一些真實地形或特殊地形的影響，使得DEM 表面存在一些凹陷的區域。在進行水流方向計算時，由于這些區域的存在，往往得到不合理的甚至錯誤的水流方向。因此，在進行水流方向的計算之前，應該首先對原始DEM 數據進行洼地填充，得到無洼地的DEM。洼地填充的基本過程是先利用水流方向數據計算出 DEM 數據中的洼地區域，并計算洼地深度，然后，依據這些洼地深度設定填充閾值進行洼地填充。

2.2.2 水流方向的確定

水流方向是指水流離開每一個柵格單元時的指向。地表徑流在流域空間內總是從地勢高處向地勢低處流動，最后經流域出口排出流域。為了準確地劃定流域界線，首先要確定水流在每個柵格單元格內的流動方向。在ArcGIS 中通過將中心柵格的8 個鄰域柵格編碼，水流方向便可以以其中的某一值來確定。

2.2.3 匯流累積量計算

在地表徑流模擬過程中，匯流累積量是基于水流方向數據計算得到的。首先基于無洼地DEM 生成水流方向圖，利用該數據，采用Hydrology 工具集中的Flow Accumulation 工具計算出匯流累積量數據。

2.2.4 水流長度計算

水流長度指地面上一點沿水流方向到流向起點（或終點）間的最大地面距離在水平面上的投影長度。它分為順流計算及溯流計算兩種，可利用Hydrology 工具集中的FlowLength 工具實現，其中計算方向分別選擇順流計算Downstream 或溯流計算Upstream。

2.2.5河網的提取

河網提取采用地表徑流漫流模型計算。首先是在無洼地DEM上利用最大坡降的方法得到每一個柵格的水流方向，然后利用水流方向柵格數據計算出每一個柵格在水流方向上累積的柵格數，即匯流累積量，所得到的匯流累積量則代表在一個柵格位置上有多少個柵格的水流方向流經該柵格；假設每一個柵格處攜帶一份水流，那么柵格的匯流累積量則代表著該柵格的水流量。基于上述思想，當匯流量達到一定值的時候，就會產生地表水流，那么所有那些匯流量大于那個臨界數值的柵格就是潛在的水流路徑，由這些水流路徑構成的網絡，就是河網。

2.2.6流域劃分

在河網提取時通過設置不同的閾值可以得到不同尺度的河網，然后利用watershed 工具生成集水區，即可以得到不同細化的集水區域，通過對比可以得到適合的流域劃分結果[4]。

3 實例分析

老撾川壙省某硅廠取水口位于Nam Pae上，該河流域面積約90km2，河道上無水文觀測站，屬于典型無資料小流域地區，本研究以該硅廠取水口設計洪水計算為例。

3.1 流域劃分及地形參數確定

根據STRM90m的數字高程資料（DEM資料），采用Acrgis軟件的ArcHydro工具進行流域劃分，得到Nam Pae流域的水系圖及Nam Pae流域地形圖，分別見圖1、圖2。

在Arcgis軟件中計算出Nam Pae流域的面積為88.9km2，其中廠址取水口處斷面控制流域面積為40.39km2，主河道平均坡度Jc為9.4‰，河長L為8.663km。Nam Pae流域氣候溫熱濕潤，河槽中水草蔓生，河道兩側及坡面上植被茂密，根據小流域河道匯流糙率系數表，匯流糙率系數Nc取0.1，流域內土壤主要為黏土、肥黏土質土壤、山地草甸土，根據下滲系數值表，穩定下滲率μ取0.12 mm/min。

3.2 暴雨參數確定

根據收集到的川壙氣象站1982～2015年最大一日降水量資料，采用Gumbel分布進行適線，分析得出不同重現期所對應的一日降水量，最大24h降雨量H24由最大一日降雨量乘以系數1.1得到，雨力Sp根據公式6求出，暴雨衰減指數n取0.8，見表1。

3.3 暴雨洪水計算結果

根據林平一法計算公式求出Nam Pae河流硅廠取水口處不同重現期對應暴雨洪水流量，見表2。

3.4 洪水水位計算

采用曼寧公式進行計算分析，該處河床為礫石、卵石組成，底坡尚均勻，床面不平整，河槽部分糙率取0.033，左側河灘長有稠密的植被，糙率取0.180，右側河灘3/4地帶長滿茂密的雜草、灌木，糙率取0.130，實測該處河道比降為0.004，根據實測水水口處河道橫斷面斷面采用曼寧公式求出Nam Pae河流硅廠取水口處斷面不同重現期洪水流量對應的洪水位

4 結語

本研究以老撾川壙省Nam Pea流域的DEM 為基礎數據，采用Arcgis軟件提取出該流域的水文特征，如控制流域面積、主河道平均坡度，河長等參數，有效解決了林平一公式在計算小流域設計洪水時參數難以獲取的問題，大大提高了該公式的實用性，使其能夠廣泛應用于老撾無資料小流域地區的設計洪水計算中。

參考文獻：

[1] 綦晶，薛蓮.無資料地區特小流域設計洪水計算方法分析[J]，黑龍江水利科技.2017；45（10）：41-2，50.

[2] 晁銳.兩種常用小流域洪水計算方法的靈敏度分析 [D]：長安大學；2014.

[3] 畢曉玲，李小娟，胡卓瑋，賴晗.DEM網格尺寸對地形因子精度的影響分析[J]，測繪科學.2012；37（6）：150-2，88.

[4] 程崢，李永勝，高微微.基于ArcGIS的DEM流域劃分[J].地下水，2011，33（06）：128-130.

作者簡介：

賈春強（1988-），男，河南省人，碩士研究生。主要從事工程水文氣象研究。E-mail：125320825@qq.com

（作者單位：1.機械工業勘察設計研究院有限公司；2.中國機械設備工程股份有限公司）