黃祥仕，黃定謀，覃麗梅

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

1 流域概況

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

圖1 百色（三）站水文站網(wǎng)、水利工程分布圖

2 水文站現(xiàn)狀

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

3 預報方案

3.1 新安江模型

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

圖2 百色（三）水文站預報方案結(jié)構(gòu)圖

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

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

3.2 多元回歸法

（1）百色水利樞紐和澄碧河水庫控制的流域面積占百色（三）水文站總流域面積的99.45%，區(qū)間面積僅占全流域面積的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）預報應用。作業(yè)預報時，利用百色水利樞紐和澄碧河水庫下泄流量合成的增量、百色（三）水文站起漲水位、最大下泄流量時百色（三）水文站同時水位等數(shù)據(jù)帶入公式（1），即可推求出百色（三）水文站的相應洪峰水位。

3.3 調(diào)洪演算法

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

水庫水量平衡方程為：

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

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

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

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

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

3.4 實例分析

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

4 結(jié)語

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