朱雪梅

（柳州市水文水資源局，廣西 柳州 545001）

無資料地區(qū)中小河流洪水預報方案構(gòu)建

朱雪梅

（柳州市水文水資源局，廣西 柳州 545001）

針對中小河流地區(qū)新建水文站缺乏實測及調(diào)查的水文資料，無法建立預報方案，水文預警預報工作限于被動局面的問題。以廣西富樂水文站為例，采用瞬時單位線法構(gòu)建洪水預報方案，經(jīng)成果驗證：該方法對中洪以上洪水的預報成果較理想，可應用于其它新建水文站中小河流地區(qū)的洪水預報，對于提高中小河流防洪減災能力具有現(xiàn)實意義。

無資料地區(qū)；中小河流；洪水預報

柳州市的中小河流點多面廣、源短流急、洪水暴漲暴落，許多中小河流沿岸的縣城、重要集鎮(zhèn)防洪標準普遍偏低，有的甚至沒有設(shè)防。近年來，受異常天氣系統(tǒng)、中小河流不合理開發(fā)等因素影響，中小河流暴雨洪水頻發(fā)，成災嚴重。2011年，柳州水文局開始實施中小河流水文監(jiān)測系統(tǒng)建設(shè)，截止2015年底，已完成8個水文站，26個水位站，90個雨量站的建設(shè)并投入使用中。此中小河流水文監(jiān)測系統(tǒng)建設(shè)的完成，填補了柳州中小河流水文信息在防汛決策中的空白，大大提高了中小河流的水文預警預報能力。但新建水文站中小河流地區(qū)，由于缺乏實測及調(diào)查的水文資料，無法建立預報方案，水文預警預報工作限于被動的局面。為此，筆者擬利用瞬時單位線法構(gòu)建無資料地區(qū)中小河流洪水預報方案，供缺乏水文資料的中小河流防洪減災決策參考。

1 瞬時單位線法

1.1 定義

單位時間內(nèi)流域上均勻分布的單位凈雨量在流域出口斷面處形成的地面徑流量過程線，簡稱單位線，記為UH。單位凈雨量通常取雨深10mm；單位歷時可以是一個時段如1、3、6h等，也可以是瞬時，即凈雨歷時趨于無限小的情況。相應于前者的單位線稱為時段單位線，相應于后者的稱為瞬時單位線。

瞬時單位線記為IUH。應用較廣的是納須的IUH模型，納須把流域匯流的調(diào)蓄作用模擬為n個相同的線性串聯(lián)水庫，推導出的IUH的理論公式是：

式中：u（O，t）——t時刻瞬時單位線縱高；

?！が敽瘮?shù)；

n——相當于線性水庫的個數(shù)（調(diào)節(jié)次數(shù)或調(diào)節(jié)系數(shù)）；

K——自然庫線性蓄洪方程的匯流歷時（反映流域匯流時間參數(shù)或調(diào)蓄系數(shù)）；

e——自然對數(shù)的底；

t——時刻。

由以上公式可知，決定瞬時單位線的參數(shù)為n，K兩參數(shù)，只要n，K確定，瞬時單位線的值即可求出。

1.2 凈雨過程的推求

瞬時單位線的計算核心為凈雨過程的推求。凈雨過程的確定采用初損和平均扣損法。時段降雨量的計算采用算術(shù)平均法。一般把起漲點以前的降雨作為初損Io扣除，以后的降雨則考慮水量平衡，減去平均入滲。如果部分雨強小于f時則采用多次扣損，直至水量平衡為止。f為產(chǎn)流期平均入滲率，系根據(jù)瞬時單位線法匯流計算的要求進行分析，求取時把起漲點以前的降雨視為初損Io扣除，漲點以后的降雨，視雨強大小，當雨強大于f時稱為產(chǎn)流期降雨，當雨強小于f時稱為不產(chǎn)流降雨。

1.3 資料準備

推求一個特定河流的匯流單位線需要準備所分析河流的流域面積，河道比降等資料。流域面積：可根據(jù)預報橫斷面所在位置，可在1/50 000地形圖上勾繪，用電子數(shù)字求積儀量計或在CAD上勾畫出流域面積。目前河湖普查已經(jīng)有詳細成果，可直接應用，有時量算出斷面以下部分，只需用全面積減去下游部分即可。每個斷面至少量計3次，各次誤差在5％以內(nèi)時，取其算術(shù)平均值為最后采用值，單位為km2。

河道比降：是反映河道高程沿河長平均變化況情，是指自從分析斷面起至分水嶺止河道干流平均坡降（見圖1）。根據(jù)沿程地形變化的特征點高程，按下式采用加權(quán)平均法計算：

圖1 河道比降計算示意圖

2 應用實例

2.1 富樂水文站簡介

富樂水文站屬省級重要水文站，是珠江流域西江水系浪溪河的唯一控制站，位于柳州市融安縣大將鄉(xiāng)富樂村浪溪河河段左岸，地理位置為東經(jīng)109°35′，北緯 25°17′，集水面積 623km2，平均坡降2.62‰。浪溪河，融江支流，因河床陡，灘多水急浪大而得名。發(fā)源于融安縣板欖鎮(zhèn)東嶺村的翁古頂，西南流經(jīng)融安縣內(nèi)的板欖、大將、長安等鄉(xiāng)鎮(zhèn)，在長安鎮(zhèn)江口村匯入融江。浪溪河沿河兩岸為山嶺連綿的低山丘陵地區(qū)。廣福頂山脈呈人字形向北邊、東邊地界展開，為高程800m以上的中山地帶，南側(cè)和西側(cè)逐漸下降為高程500m左右的山丘陵，浪溪河位于西側(cè)中下游，地勢降到高程200m以下，南方暖濕氣流，因地勢抬升，易于成云致雨。流域溫暖潮濕，森林植被茂盛。

2.2 匯流單位線計算

查《廣西暴雨徑流查算圖表》[1]可知浪溪河富樂水文站：產(chǎn)流分區(qū)屬第2區(qū)，匯流分區(qū)屬第一（2）區(qū)但在第一（1）區(qū)附近，雨型分區(qū)屬第1區(qū)，下墊面分區(qū)屬丘陵區(qū)，退水分區(qū)屬2區(qū)，富樂水文站所在地區(qū)為非巖溶地區(qū)，瞬時單位線主要參數(shù)m1穩(wěn)、n采用非巖溶地區(qū)計算公式，即

匯流分區(qū)第一（1）區(qū)：

因富樂水文站處在匯流分區(qū)第一（2）區(qū)但在第一（1）區(qū)附近，則m1穩(wěn)取第一（1）區(qū)和第一（2）區(qū)的平均值，即m1穩(wěn)=（7.34+5.95）/2=6.65，n采用匯流分區(qū)第一（2）區(qū)公式計算，即

式中：m1穩(wěn)——單位線匯流參數(shù)；

n——相當于線性水庫的個數(shù)（調(diào)節(jié)次數(shù)或調(diào)節(jié)系數(shù)）；

K——自然庫線性蓄洪方程的匯流歷時（反映流域匯流時間參數(shù)或調(diào)蓄系數(shù)）；

C——匯流速度；

F——流域面積，km2；

J——河道比降；

△t——時段長。

S曲線可查《廣西壯族自治區(qū)暴雨徑流查算圖表》[1]附錄表5：先確定t/K，然后根據(jù)n從表中讀出S（t）值，將S（t）錯后一個時段相減，即得無因次單位線u（△t、t），再乘以F/△t×3.6即得流量為單位的時段單位線q（△t，t），富樂水文站匯流單位線見圖2。

圖2 富樂水文站匯流單位線

2.3 凈雨計算（以“2013.6”暴雨洪水為例）

2.3.1 降雨過程摘錄

2013年6月8日5時～6月10日11時，浪溪河流域普降中到大雨，局部暴雨到大暴雨，暴雨中心位于浪溪河流域的下游，累積最大降雨量大坡站232.5mm，強降雨主要集中在6月9日5時-15時，是富樂站“2013.6”暴雨洪水的造峰雨。取浪溪河富樂站以上雅瑤、板欖、沙木田、中村、車平、東嶺、里當、富樂共8個雨量站的算術(shù)平均作為浪溪河面平均降雨量。浪溪河流域6月9日5時-15時降雨過程見表1。

表1 浪溪河流域各站點6月9日5時-15時降雨過程摘錄

2.3.2 時段凈雨過程

依據(jù)富樂水文站所在地屬于第2產(chǎn)流區(qū)，查降雨徑流相關(guān)特征參數(shù)綜合表得Wm=100mm，K=0.93。土壤前期含水量W0=KWm=0.93×100=93mm，則土壤前期缺水量I0=100-93=7mm。根據(jù)浪溪河流域下墊面條件，取產(chǎn)流期平均入滲率為7mm/h。則時段凈雨量=時段降雨量-前期缺水量-產(chǎn)流期平均入滲率。富樂水文站“2013.6”強降雨時段凈雨過程見表2。

表2 富樂水文站“2013.6”強降雨時段（6月9日）凈雨過程計算 mm

2.4 地表徑流過程計算

2.5 地下徑流過程計算[1]

（1）地下徑流過程的起漲部分采用三角形分配，退水部分則根據(jù)退水曲線公式推算。地下水峰值置于計算得的地表徑流基本結(jié)束附近；

（2）富樂水文站所在地退水分區(qū)屬2區(qū)，查退水指數(shù)a1、a2分區(qū)綜合表，結(jié)合富樂站歷年分析資料，取a1＝0.35，a2=0.16，t′=24.7×a1^（-1.06）=24.7×0.35^（-1.06）=75.2h=3.1d。

2.6 洪水過程計算

洪水過程=地表徑流過程+地下徑流過程。“2013.6”暴雨洪水過程計算成果與實際洪水過程對比見圖3。

圖3 “2013.6”暴雨洪水過程對比圖

從洪水過程計算成果可知，富樂水文站2013年6月9日17時出現(xiàn)洪峰流量為741m3/s，由演算的洪峰流量741m3/s查Z～Q曲線得洪峰水位為156.68m；實際洪水過程是2013年6月9日16：15出現(xiàn)洪峰水位156.78m，相應流量為765m3/s。洪峰流量僅相差24.0m3/s，洪峰水位僅相差10cm，洪水預見期2h，預報成果較好。從圖3可可看出：洪水預報的漲水過程與實際漲水過程比較吻合；洪水預報的峰后洪量偏小較多，原因是6月9日16時至6月10日11時，還有27.9mm的累計降雨過程，因此實際退水過程比推算的退水過程緩慢，峰后洪量偏小。

2.7 預報成果的評定

根據(jù)《水文情報預報規(guī)范》SL250-2000規(guī)定[2]，河道流量（水位）預報以預見期內(nèi)變幅的20％作為許可誤差。當流量許可誤差小于實測值的5％時，取流量實測值的5％作為許可誤差。當水位許可誤差小于實則洪峰流量的5％所相應的水位幅度值或小于0.1m時，則以該值作為許可誤差。峰現(xiàn)時間預報許可誤差是以預報根據(jù)時間至實測洪峰出現(xiàn)時間之間時距的30％計算，當許可誤差小于3h或一個計算時段長時，則以3h或一個時段長為許可誤差。

富樂站的洪水預報成果按上述規(guī)定進行評定。從洪水過程計算成果可知，本次預報洪峰流量預報誤差24m3/s，允許誤差是38.3m3/s；洪峰水位預報誤差是0.1m，允許誤差是0.71m；峰現(xiàn)誤差是0.75h，允許誤差是3h，預報成果合格。

2.8 預報方案檢驗

方案檢驗采用2011～2015年共3場中洪以上洪水進行檢驗，檢驗期洪峰水位合格3場，合格率100％；峰現(xiàn)時間合格3場，合格率100％；預見期2～5h，可用于洪水預報作業(yè)。

3 結(jié)語

瞬時單位線法主要是由暴雨過程推求流量過程，僅適用于流域面積小于1 000km2的中小河流域的洪水計算，是解決無實測水文資料地區(qū)洪水計算的有效方法，可應用于其它新建水文站中小河流地區(qū)的洪水預報，在進行預報時應結(jié)合降雨數(shù)值預報、根據(jù)前期降雨及流域下墊面調(diào)查分析土壤最大蓄水量、確定初損與后損，才能提高洪水預報的精度與延長預見期。新建水文站的中小河流域地區(qū)，應加強測流工作，盡快建立水位～流量關(guān)系曲線，才能預報出水位過程。

[1] 廣西水文總站.廣西壯族自治區(qū)暴雨徑流查算圖表[R].南寧：廣西水文總站，1984.

[2]SL250-2000，水文情報預報規(guī)范[S].

（責任編輯：周 群）

Set up flood forecasting scheme for middle or small river basin without observation data

ZHU Xue-mei
（Hydrology and Water Resources Bureau of Liuzhou,Liuzhou 545001,China）

For a new hydrological station in a middle or small river basin,it is difficult to set up a forecasting scheme due to shortage of observation data.At the Fule Hydrological Station in Guangxi,instantaneous unit hydrograph method had been adopted to set up the flood forecasting scheme.It was verified by actual measurement that this scheme could lead to rather satisfying forecasting results in case of moderate or bigger flood and should be popularized to the new hydrological stations in middle or small river basins,so as to improve the flood control capability.

Region without observation data;middle and small rivers;flood forecasting

P338

B

1003-1510（2016）04-007-04

2016-01-09

朱雪梅（1980-），女，廣西武宣人，柳州市水文水資源局工程師，學士，從事水文情報預報工作。