右江河庫預報調度中百色（三）水文站的洪水預報方法

黃祥仕，黃定謀，覃麗梅

(百色水文中心，廣西 百色 533000)

1 流域概況

右江位于西江最大的一級支流郁江的上游，為珠江流域西江水系，發源于云南省廣南縣，流域地勢從西向東傾斜，多為山地地帶，植被覆蓋較好，主流不易沖刷。主干流的上游為馱娘江，從云南省剝隘鎮境內轉入百色市境內后稱為右江，右江在百色市境內的集水面積為21 903 km2，全長208 km，自然落差362 m。百色（三）水文站位于百色市城區右江河上，屬國家基本水文站，集水面積21 720 km2。右江百色水利樞紐位于百色（三）水文站上游22 km，集水面積19 600 km2，澄碧河水庫位于百色（三）水文站上游7 km 的支流澄碧河上，集水面積2000 km2，右江那吉航運樞紐位于百色（三）水文站下游約41 km（見圖1）。百色（三）水文站屬于典型的上下游都有水庫影響的水文站，河庫預報調度一體化研究非常有必要。

圖1 百色（三）站水文站網、水利工程分布圖

2 水文站現狀

百色（三）水文站是右江上游段的主要控制水文站，其上游百色水利樞紐2006 年開始下閘蓄水，其下游那吉航運樞紐2008年開始蓄水運行，其回水頂托影響到百色（三）水文站斷面，百色（三）水文站斷面就完全處在那吉航運樞紐庫區里。因此2008年之后，百色（三）水文站斷面就失去了天然河道的屬性，其洪水形成完全受百色水利樞紐、那吉航運樞紐、澄碧河水庫共同調度影響，3個水庫影響的因素各不相同，影響的組合比較復雜。1999年修編的百色（三）水文站的天然河道洪水預報方案已經完全不能適用受水庫調度影響的情形，因此在2010年《廣西洪水預報方案》修編中，不再將百色（三）水文站洪水預報方案納入其中。為了響應全區水文推行“預測-預警-預報”的水情分析工作模式，以及推行河庫預報調度一體化工作，將傳統的“哨兵”模式向“偵察兵”模式轉變，需要研究出多種適用于現階段百色（三）水文站的洪水預報方法，準確地預報出百色（三）水文站洪水過程，充分發揮水文的“尖兵”“耳目”作用，為防汛減災決策提供技術支撐。

3 預報方案

3.1 新安江模型

（1）百色（三）水文站的新安江模型洪水預報方案的河道匯流采用MSK法，區間產匯流采用SMS_3和LAG_3法，方案是在中國洪水預報系統平臺進行編制，共設置3 個方案輸入：百色水利樞紐壩下站、澄碧河水庫壩下站流量及區間。百色水利樞紐壩下站、澄碧河水庫壩下站輸入采用馬斯京根河道連續演算法（MSK），區間輸入采用蓄滿產流模型（SMS_3）和滯后演算模型（LAG_3）；雨量站控制權重采用泰森多邊形法。方案計算時段為1 小時，方案輸出類型為水位流量（模型計算結果為流量，水位通過水位流量綜合關系曲線進行反推）。預報方案結構圖見圖2。

圖2 百色（三）水文站預報方案結構圖

（2）預報方案率定選取2008—2014 年資料，所選取的百色（三）水文站水文資料經過整編，資料系列具有一致性、代表性和可靠性，而百色水利樞紐壩下和澄碧河壩下所采用的數據為實測數據，未經過資料整編。方案率定的確定性系數為0.910，達到甲級標準。方案檢驗確定性系數0.782，達到乙級標準。

（3）預報應用。利用中國洪水預報系統平臺進行百色（三）水文站作業預報，選擇需要的預見期和遇見期內的雨量預報數值（可選中國、歐洲或日本），然后進行洪水作業預報，方案能夠預報出預見期內每小時的水位流量過程，根據實際情況可對預報結果進行實時校正，該方案的洪水預報應用效果良好。

3.2 多元回歸法

（1）百色水利樞紐和澄碧河水庫控制的流域面積占百色（三）水文站總流域面積的99.45%，區間面積僅占全流域面積的0.55%，百色水利樞紐和澄碧河水庫的下泄流量合成基本代表了百色（三）水文站的流量，但由于受下游那吉航運樞紐下閘蓄水的頂托影響，百色水利樞紐和澄碧河水庫的下泄流量跟百色（三）水文站的水位無法建立單一的相關關系。在百色水利樞紐和澄碧河水庫的下泄流量較小或變化不大時，百色（三）水文站相應的水位變化不大，百色水利樞紐和澄碧河水庫的下泄流量變化較大時，百色（三）水文站水位才有明顯變化。因此多元回歸法選用百色水利樞紐和澄碧河水庫下泄流量的合成流量作為上游來水量，用后一個時段的上游來水量與前一個時段的上游來水量比值，作為上游來水增量，利用上游來水增量、百色（三）水文站起漲水位、百色（三）水文站同時水位與百色（三）水文站洪峰水位等建立Z峰～（Q增量、Z起、Z同）多元回歸方程。方案采用2009—2020 年實測資料進行分析，建立關系式為：

式中：Z峰為百色（三）水文站洪峰水位，m；Q增量為百色水利樞紐和澄碧河水庫的相應下泄流量之和的增量，m3/s；Z起為百色（三）水文站起漲水位，m；Z同為最大下泄流量時百色（三）水文站同時水位，m。

（2）本方案選用2009—2020年期間百色水利樞紐和澄碧河水庫下泄流量之和大于660 m3/s 時，百色（三）站不同的水位流量級共31 場洪水資料點子進行編制，方案合格率為65%。采用2018年6場洪水資料點子來檢驗，檢驗合格率為67%。方案精度等級屬于丙級，可用于參考預報。

（3）預報應用。作業預報時，利用百色水利樞紐和澄碧河水庫下泄流量合成的增量、百色（三）水文站起漲水位、最大下泄流量時百色（三）水文站同時水位等數據帶入公式（1），即可推求出百色（三）水文站的相應洪峰水位。

3.3 調洪演算法

由于百色（三）水文站受那吉航運樞紐頂托，百色（三）水文站位于庫區內，而百色（三）水文站與那吉航運樞紐壩首距離41 km，在不同的水位級有不同的庫面水位落差。因此利用那吉航運樞紐當前庫水位、庫容曲線、出庫流量和入庫流量（百色水利樞紐+澄碧河水庫+區間的合成流量）演算那吉航運樞紐庫水位未來變化過程，然后加上庫面落差即求得百色（三）水文站水位變化過程。

水庫水量平衡方程為：

式中：Q平均為時段內那吉航運樞紐入庫流量的平均值，m3/s；q平均為時段內那吉航運樞紐出庫流量的平均值，m3/s；ΔV為時段內那吉航運樞紐蓄水量的變化，億m3；A平均為時段內那吉航運樞紐庫面面積平均值，m2；P、E分別為時段內那吉航運樞紐庫面積的降雨量及水面蒸發量，mm；Δt為計算時段時長，h。

百色（三）水文站斷面與那吉航運樞紐壩首之間的庫面水位落差公式為：

式中：Z百為百色（三）水文站水位，m；Z那為同一時間那吉航運樞紐庫水位，m；S為百色（三）水文站與那吉航運樞紐的不同水位級平均庫面水位落差，m。

由于庫面比例不大，同時計算時段較短（1 小時），則公式（2）的（P-E）A平均項可不計。

預報應用：首先，按合成流量預報方法預報那吉航運樞紐逐時入庫流量，進而求得各時段內的平均入庫流量；其次，由已知的當前庫水位、入庫流量過程和水庫調度計劃來推算蓄水量變化過程，然后根據庫容曲線推求那吉航運樞紐庫水位變化過程；最后根據百色（三）水文站與那吉航運樞紐的平均庫面落差推求百色（三）水文站水位變化過程。在實際應用中，已經將各種數據存于計算機內，輸入當前庫水位和水庫調度計劃即可自動計算產生結果。根據調洪演算法預報的百色（三）水文站水位結果準確度相對較高。

3.4 實例分析

以2018年8月8日洪水為例進行分析。2018年8月4～7日，受季風槽影響，百色市出現持續性強降雨天氣過程。8 月6 日08 時，百色水利樞紐庫水位已經超汛限水位0.28 m，而當時上游入庫流量為1950 m3/s，出庫流量是四臺機組滿發共702 m3/s，庫水位仍持續上漲。根據預報，8月8日入庫流量將達到峰值約3500 m3/s，為了確保安全，百色水利樞紐決定于8月7日15時開閘泄洪，泄洪流量（不含發電流量）暫按1000 m3/s以內控制。8月6日20時，根據百色水利樞紐、那吉航運樞紐與澄碧河水庫調度計劃，利用新安江模型、多元回歸法、調洪演算法分別分析計算右江百色水文站的洪峰水位，結果見表1。從本場洪水的預報結果與實際來看，新安江模型預見期最長，調洪演算法預報水位結果更接近于實際，從歷年多場洪水的預報結果來看，大都如此。

4 結語

由于百色（三）水文站受上下游水庫調度共同影響，百色水利樞紐、那吉航運樞紐、澄碧河水庫3個水庫影響的因素各不相同，影響的組合比較復雜，因此方案在實際應用中存在一定的局限性：一是百色水利樞紐下泄流量傳播到百色（三）水文站斷面只有3 h，實測數據分析預見期很短，如果要延長預見期，則需要利用百色水利樞紐及澄碧河水庫的未來調度計劃參與計算分析；二是預報結果受下游那吉航運樞紐調度的影響因素較大，計算分析出預報結果后，那吉航運樞紐調度計劃有變化（加大或減小出庫流量）時，需要及時調整預報結果；三是多元回歸法適用于上游水庫出庫流量變化較大的洪水，當上游水庫出庫流量相對穩定時，方案的適用性不強；四是方案所用的資料系列不長，僅為10年左右，今后需要收集更多的資料系列加入分析，不斷修訂完善預報方案，提高預報精度?；谝陨显?，在實際的洪水預報作業中，要多種方法一起進行綜合分析預報，發布預報結果后，要隨時關注有關聯的水庫調度情況，一旦有調度變化要及時調整預報結果，這樣才能變被動為主動，更好地發揮河庫預報調度在防汛抗旱中的參謀作用。