何素明　譚喬鳳

摘要：以漓江三期補水工程體系為研究對象，對漓江實時補水優(yōu)化調(diào)度進行研究。首先基于馬斯京根流量演算法建立了供水水庫群放水和控制斷面需水的定量響應(yīng)關(guān)系，在此基礎(chǔ)上提出考慮河道水流演進滯時和生態(tài)需水要求的漓江水庫群聯(lián)合補水實時優(yōu)化調(diào)度模型及相應(yīng)的調(diào)度規(guī)則，并通過約束處理技術(shù)使模型轉(zhuǎn)化為能通過線性規(guī)劃求解的問題。實時優(yōu)化調(diào)度模型在漓江流域的應(yīng)用結(jié)果表明：斧子口、川江、小溶江、青獅潭、五里峽、思安江六庫聯(lián)合調(diào)度，可以使特枯年份的供水日保證率達到94%以上。研究對維持漓江河道健康、提高通航保證率以及旅游經(jīng)濟效益具有重要的意義。

關(guān)鍵詞：實時調(diào)度；優(yōu)化調(diào)度；馬斯京根；線性規(guī)劃算法；調(diào)度規(guī)則

中圖分類號：TV213文獻標志碼：A文章編號：

16721683（2018）04009806

Realtime waterssupply optimal operation in Li River Basin

HE Suming1，TAN Qiaofeng2，LEI Xiaohui3，LI Hongliang1，WANG Xu3，YANG Mingxiang3，ZHANG Peilun3

（

1.Guangxi Water and Power Design Institute，Nanning 530023，China；

2.College of Water Resource & Hydropower，Sichuan University，Chengdu 610065，China；

3.China Institute of Water Resources and Hydropower Research，Beijing 100038，China

）

Abstract：

Choosing the thirdphase watertransfer project as the research object，we studied the realtime watertransfer optimal operation in Li River Basin.Firstly，we built the quantitative response relationship between the release of the watersupply reservoirs and water demand of the control section based on Muskingum′s flow algorithm.Then we proposed a realtime watertransfer optimal operation model and corresponding operating rules for the reservoirs group on Li River，which took into consideration the water flow evolution delay and the ecological water demand.Through constrainthandling techniques，we transformed the model into a question which could be solved by linear programming.The application results showed that the joint operation of the six reservoirs on Li River （including Fuzikou，Chuanjiang，Xiaorongjiang，Qingshitan，Wulixia and Si′anjiang reservoirs）could raise the watersupply daily assurance rate to above 94% in especially dry years.This study has important significance for maintaining a healthy river channel，increasing navigation assurance rate，and improving economic benefits of tourism on Li River.

Key words：

realtime operation；optimal operation；Muskingum；linear programming；operating rules

水庫實時調(diào)度是一個“預(yù)報、決策、實施、再預(yù)報、再決策、再實施”的滾動向前過程，不僅要顧及徑流的長期變化規(guī)律，[HJ1.85mm]又要適應(yīng)當前水情變化[1]。利用優(yōu)化技術(shù)指導(dǎo)水庫實時調(diào)度，一直是水資源管理者嘗試突破的關(guān)鍵科學(xué)問題[234]。為了提高漓江流域供水效益：于曉[5]根據(jù)是否考慮森林蓄水工程對漓江補水的影響，對斧子口、川江和小溶江水庫的月供水調(diào)度方案進行了研究；李平[6]采用HJHAT法確定了漓江流域的環(huán)境流量，并根據(jù)區(qū)域內(nèi)耗水情況以及漓江補水目標，設(shè)置青獅潭、斧子口、川江和小溶江四庫聯(lián)合長期調(diào)度方案；李若男等[7]通過建立漓江流域水環(huán)境模型，模擬了青獅潭水庫補水及非補水模式下，枯水期其下游河道內(nèi)水環(huán)境變化情況。這些研究對漓江供水效益的提高具有重要的意義，但是大都集中于中長時間尺度及規(guī)劃調(diào)度層面。

隨著漓江三期補水工程體系逐漸形成，漓江水庫調(diào)度研究需要逐漸從規(guī)劃設(shè)計階段向?qū)崟r管理階段過渡。實時補水調(diào)度中，由于補水距離長，補水斷面不能對水庫放水作出即時響應(yīng)。受水流演進滯時影響，各補水水庫的放水經(jīng)過河道水流演進到達需水斷面的流量可能遠大于需求，而使水庫水位降低，不利于后期補水效益發(fā)揮；也可能遠小于需求，無法滿足沿經(jīng)河段的生態(tài)需水及漓江河段補水需求，而帶來生態(tài)破壞以及旅游、航運等經(jīng)濟損失。因此，為了科學(xué)合理地利用各水庫的補水潛力，迫切需要建立一個能考慮河道水流演進滯時及生態(tài)需水要求的漓江水庫群聯(lián)合補水實時優(yōu)化調(diào)度模型。

1漓江補水工程體系

漓江補水分三期實施[10]：一期補水工程以青獅潭水庫作為補水水源，通過減少灌溉用水、調(diào)整水庫調(diào)度等方式向漓江補水，設(shè)計補水使枯水期流量達到30 m3/s；二期補水工程五里峽水庫和思安江水庫與青獅潭水庫聯(lián)合運用，設(shè)計補水使枯水期流量達到42 m3/s；三期補水工程需結(jié)合斧子口、小溶江、川江水庫，與一、二期補水工程聯(lián)合向漓江補水，以使枯水期流量達到60 m3/s。各補水水庫的主要特性參數(shù)見表1。

漓江補水對象具體指桂林陽朔河段共86 km的河段，以桂林水文站作為漓江補水的控制斷面。斧子口、川江、小溶江、青獅潭和五里峽水庫位于桂林斷面以上，思安江水庫位于桂林斷面以下的支流上，其網(wǎng)絡(luò)概化見圖1。圖中，斧子口、川江、小溶江和青獅潭水庫放水至桂林斷面的演進河段按照河道上下游關(guān)系編號。以斧子口水庫為例，其放水要到達桂林斷面，需要經(jīng)過斧子口-大溶江河段（河段1）、大溶江-上支鋪河段（河段3）、上支鋪-匯合口河段（河段5）、匯合口-桂林河段（河段8）。五里峽水庫為跨流域引水工程，通過干渠連接到漓江流域補水，而思安江水庫位于桂林斷面以下，桂林斷面對于這兩個水庫的補水響應(yīng)通過設(shè)計時采用的補水有效系數(shù)法表示，即有效補水量等于水庫出庫流量乘以補水有效系數(shù)。五里峽和思安江水庫的補水有效系數(shù)分別為085，065。

2.1優(yōu)化目標

以漓江三期補水工程體系為研究對象，綜合考慮河道調(diào)蓄作用、庫容和桂林斷面需水等約束，以斧子口、小溶江、川江、青獅潭、五里峽、思安江水庫群的總放水最小化為目標建立實時補水優(yōu)化調(diào)度模型。目標函數(shù)為：

min∑[DD（]6[]i=1[DD）]Ri，t[JY]（1）

式中：Ri，t表示水庫i在t時段的放水流量。該目標函數(shù)既能使水庫以最小的放水量滿足當前時段的流量需求，又能保證水庫群的總蓄水量盡可能大以保證后期供水，能兼顧當前供水和未來供水，有利于供水保證率的提高。

2.2約束條件

2.3約束處理

六庫的放水需要一定的河道調(diào)蓄作用才能到達桂林斷面。河道水流演進包括水力學(xué)方法[810]和水文學(xué)方法兩種[1115]。水力學(xué)方法嚴格基于圣維南方程組，能得到精細的演進過程，但是其對河道沿程斷面數(shù)據(jù)要求較高；水文學(xué)方法以馬斯京根法為代表，通過對圣維南方程組進行概化以簡化計算過程。由于漓江流域六個水庫至桂林斷面的數(shù)據(jù)不支撐建立水動力學(xué)模型，本文采用馬斯京根法進行河道水流演進。但馬斯京根流量演進滯時特性[1617]使得桂林斷面對于水庫放水不能得到即時響應(yīng)，使實時補水調(diào)度變得尤為復(fù)雜。如果能建立考慮馬斯京根流量演進滯時的水庫出庫和桂林斷面需水的響應(yīng)關(guān)系，問題能大大簡化。下面將通過對約束條件進行處理，實現(xiàn)這一目標。

（1）約束條件（5）和（6）處理。

從馬斯京根演進公式可以看出，桂林斷面t時段的流量與各水庫t時段的放水有關(guān)，還與各個水庫從出庫演進至桂林斷面過程中，涉及到的各個演進斷面t-1時段的流量狀態(tài)有關(guān)。桂林斷面的流量由各水庫的演進流量和區(qū)間流量共同構(gòu)成，約束條件（5）和（6）可以轉(zhuǎn)化為下式：

∑[DD（]6[]i=1[DD）]Ai×Ri，t+Bi，t-1+QJGL，t≥60[JY]（9）

式中：QJGL，t為時段t各水庫至桂林斷面的區(qū)間徑流；

Ai為各水庫t時段的放水量經(jīng)馬斯京根分段演進不斷累積相乘的常數(shù)項系數(shù)，該系數(shù)只與水庫i放水流經(jīng)河段j的馬斯京根參數(shù)C0有關(guān)，具體的計算公式可表示為：

Ai=∏[DD（][]j∈Φ（i）[DD）]Cj，0[JY]（10）

式中：Φ（i）為水庫i放水演進至桂林斷面所經(jīng)過的河段集合，Cj，0為河段j的馬斯京根演進系數(shù)C0；

Bi，t-1為馬斯京根演進涉及到的t-1時段的流量演進至桂林斷面的常數(shù)項，該常數(shù)項不僅與水庫i放水流經(jīng)河段j的馬斯京根參數(shù)C1、C2有關(guān)，還與各個演進斷面t-1時段的流量狀態(tài)有關(guān)，具體的計算公式可以表示為：

Bi，t-1=∑[DD（][]j∈Φ（i）[DD）][（Cj，1×Ij，t-1+Cj，2×Qj，t-1）×∏[DD（][]k∈Φ（j）[DD）]Ck，0][JY]（11）

式中：Ij，t-1，Qj，t-1分別為河道j時段t-1上斷面和下斷面流量；Φ（j）為河段j演進至桂林斷面涉及到下游河段的集合。Φ（i）為水庫i演進至桂林斷面涉及到的河段的集合。j和k為河段編碼。由于t-1時段的流量狀態(tài)已知，所以Bi，t-1為常數(shù)項。

（2）約束條件（2）和（3）處理。

在實時補水調(diào)度中，最高庫容約束可轉(zhuǎn)化為最小放水量約束，最小庫容約束可轉(zhuǎn)化為最大放水量約束。最終的流量約束為：

[JB（]Ri，t[TXX-]≤Ri，t≤Ri，t[TX-]

Ri，t[TXX-]=max（Qi，in，t-（Vi，t，max-Vi，t）/Δt，Ri，t，min）

Ri，t[TX-]=minQi，in，t-（Vi，t，min-Vi，t）/Δt，Ri，t，max）[JB）][JY]（12）

考慮到各河段的最小生態(tài)環(huán)境流量需求，Ri，t，min取各個水庫的生態(tài)環(huán)境基流。

（3）水庫水量平衡約束和水位庫容關(guān)系曲線，用于計算本時段的調(diào)度末狀態(tài)，也即下一時段的調(diào)度初始狀態(tài)，不作處理。

2.4模型求解

約束條件經(jīng)過處理之后，得到實時補水優(yōu)化調(diào)度模型為：[HJ1.8mm]

目標函數(shù)。

（1） 輸入各水庫和各水庫至桂林斷面區(qū)間的實時徑流預(yù)報信息，假定各水庫按最大允許放水流量放水，若演進至桂林斷面的流量小于60 m3/s，則各水庫按最大允許放水流量放水；否則，轉(zhuǎn)入（2）。

（2） 假定各水庫按最小允許放水流量放水，若演進至桂林斷面的流量大于60 m3/s，此時各水庫按最小放水流量放水即可滿足當前時段桂林斷面的流量需求?？紤]到馬斯京根流量演進的滯時特性，若當前時段各水庫放水太小，河道中的流量衰減劇烈，下一時段需要很大的放水才能滿足桂林斷面的需水要求，會導(dǎo)致各水庫出庫和桂林斷面的流量劇烈變化。因此，除了五里峽水庫外，給定其余水庫均按10 m3/s 放水。五里峽水庫按其渠道最大輸水能力8 m3/s 放水。當（1）和（2）均不能滿足時，各水庫放水存在優(yōu)化空間，轉(zhuǎn)入（3）。

（3）通過線性規(guī)劃模型得到最優(yōu)放水決策。

4調(diào)度結(jié)果

現(xiàn)收集到桂林水文站特枯年份1961年枯水期9月1日至次年2月28日的逐日徑流資料，各水庫入庫徑流資料和各庫至桂林斷面的區(qū)間徑流資料通過流域面積比法求解。采用實時優(yōu)化調(diào)度模型得到各水庫的調(diào)度結(jié)果見表2。

從調(diào)度結(jié)果表可以看出，桂林斷面的實時流量分為以下三種情況。

（1） 桂林斷面流量為60 m3/s。表示滿足進入實時優(yōu)化調(diào)度模型的條件，可通過線性規(guī)劃模型得到各個水庫的最優(yōu)放水量。

（2）桂林斷面的流量大于60 m3/s，如時段5和時段7。這主要是因為梯級水庫的前一刻為了滿足60 m3/s的流量要求而使水庫加大放水，但由于馬斯京根的后效性，當前時刻桂林斷面才對前一時段的放水作出響應(yīng)，使得水庫按最小放水流量放水演進至桂林斷面的流量仍然大于60 m3/s。馬斯京根的后效性使得桂林斷面的流量存在波動，當前時刻滿足桂林斷面60 m3/s的需求，可能導(dǎo)致下一時刻桂林斷面的流量大于60 m3/s。

（3） 桂林斷面流量小于60 m3/s，如時段170、時段173等。此時水庫已經(jīng)按最大允許放水量下泄，仍然不能滿足60 m3/s的流量需求。以缺水最嚴重的時段173為例：小溶江、青獅潭已經(jīng)達到死水位限制，只能按出入庫平衡進行放水；思安江和川江水庫的上一時段來水大于放水，有一定的庫存，該時段按照來水加庫存流量進行放水；斧子口和五里峽水庫已經(jīng)按最大允許放水量和最大渠道過水能力進行放水。各個水庫均按照自身的最大供水能力供水，但是桂林斷面的流量仍不足60 m3/s。

圖3為各庫的出庫流量過程。從調(diào)度過程可以看出青獅潭、小溶江、思安江、川江、斧子口的集中補水時段分別為1至153、1至84、77至120、84至159和153至180。這五個水庫距離桂林斷面的距離分別為32 km、40 km、 45 km、57 km與58 km，因此滿足“就近補水”原則。值得注意的是：雖然圖3顯示各庫流量存在很強的波動，但是調(diào)度時段為1日，不存在閘門的頻繁啟閉，且流量變幅已經(jīng)考慮在最大和最小放水流量約束中，各庫的流量變幅控制在可接受的范圍內(nèi)。

圖4為桂林斷面的流量過程，可以看出六庫聯(lián)合調(diào)度絕大多數(shù)情況下能使桂林斷面達到60 m3/s以上的補水要求。其中，流量剛好等于60 m3/s ，大于60 m3/s以及小于60 m3/s的時段總數(shù)分別為134、37、10。有7403 %的時段能使桂林斷面流量剛好維持在60 m3/s，滿足60 m3/s的供水日保證率為9448 %，超過設(shè)計階段要求的補水目標（9008 %）。另外，從圖4還可以看出，缺水時段主要集中在時段155-175。這主要是由于該時段范圍內(nèi)水庫來水量?。ㄒ姳?）且余留水量已經(jīng)不多，受到最小庫容約束，使得最大允許放水流量大大減小，水庫按最大放水流量決策仍然不能滿足桂林斷面的流量需求，如表2中時段170、時段173等。時段176以后各庫的來水明顯增加，水庫余留水量增多，又能保證桂林斷面的流量需求。

5結(jié)論

本文針對漓江三期補水工程體系，建立了考慮馬斯京根演進滯時和生態(tài)環(huán)境基流的實時補水優(yōu)化調(diào)度模型，并提出了相應(yīng)的實時補水調(diào)度規(guī)則。通過約束處理技術(shù)將模型轉(zhuǎn)化成能用線性規(guī)劃求解的問題。模型在漓江實時補水優(yōu)化調(diào)度中的應(yīng)用表明：

（1） 實時補水優(yōu)化調(diào)度模型可以根據(jù)桂林斷面流量情況實時確定水庫的最優(yōu)調(diào)度方案，模型具有可行性和一定的優(yōu)越性；

（2） 為了提高補水效益，除跨流域引水工程五里峽水庫以外，其他水庫應(yīng)按距離桂林斷面的距離，按“就近原則”安排補水順序；

（3） 實時優(yōu)化調(diào)度模型能使特枯年份桂林斷面流量有7403 %的時段剛好維持在60 m3/s，滿足需水要求的供水保證率為944 %。因此水庫用盡可能少的放水量滿足桂林斷面的流量需求，有利于漓江供水工程長期效益的發(fā)揮。

研究屬于規(guī)劃設(shè)計階段的內(nèi)容，在制定實時調(diào)度規(guī)則的時候沒有考慮徑流預(yù)報的不確定性，規(guī)則實際應(yīng)用效果跟流域徑流預(yù)報水平息息相關(guān)。提高漓江流域預(yù)報水平，再輔以科學(xué)的調(diào)度手段，漓江流域的缺水狀況可得到有效的緩解。

參考文獻（References）：

[1]邱林，陳守煜.水電站水庫實時優(yōu)化調(diào)度模型及其應(yīng)用[J].水利學(xué)報，1997（3）：7477.（QIU L，CHEN S Y.Optimum model for realtime hydroelectric power station reservoir operation[J].Journal of Hydraulic Engineering，1997（3）：7477.（in Chinese））

[2]謝崇寶，袁宏源，郭元裕.水庫實時優(yōu)化調(diào)度模糊隨機模型[J].水電能源科學(xué)，1994（3）：194199.（XIE C B，YUAN H Y，GUO Y Y.The fuzzy stochastic model for a reservoir time optimal operation[J].Water Resources and Power，1994（3）：194199.（in Chinese））

[3]李梅，劉俊萍，黃強，等.水庫實時優(yōu)化調(diào)度余留庫容的云決策方法研究[J].西北農(nóng)林科技大學(xué)學(xué)報（自然科學(xué)版），2007，35（3）：238244.（LI M，LIU J P，HUANG Q.Cloud decisionmaking method for the remaining capacity of reservoir in realtime optimal operation[J].Journal of North West A & F University （Natural Science ），2007，35（3）：238244.（in Chinese））

[4]魏明華，萬芳，楊鑫，等.水庫群實時供水優(yōu)化調(diào)度研究[J].人民黃河，2013，35（8）：5659.（WEI M H，WAN F，YANG X.Study on the realtime multistorage system water supply optimization dispatching[J].Yellow River，2013，35（8）：5659.（in Chinese））

[5]于曉.變化環(huán)境下漓江流域上游水庫群防汛抗旱調(diào)度研究[D].桂林：桂林理工大學(xué)，2011.（YU X.The study of drought and flood prevention of the upstream reservoir group in Li River under the changing environment[D].Guilin：Guilin University of Technology，2011.（in Chinese））

[6]李平.基于河流環(huán)境流量的漓江上游水庫優(yōu)化調(diào)度[D].桂林：桂林理工大學(xué)，2009.（LI Ping.Reservoir optimal operation in the upstream of Li River based on river environmental discharge[D].Guilin：Guilin University of Technology，2009.（in Chinese））

[7]李若男，陳求穩(wěn)，蔡德所，等.漓江枯水期水庫補水對下游水環(huán)境的影響[J].水利學(xué)報，2010，39（1）：716.（LI R N，CHEN Q W，CAI D S.The impacts of flow recharge in dry season on the downstream aquatic environment of the Lijiang River[J].Journal of Hydraulic Engineering，2010，39（1）：716.（in Chinese））

[8]葛守西，程海云，李玉榮.水動力學(xué)模型卡爾曼濾波實時校正技術(shù)[J].水利學(xué)報，2005，36（6）：687693.（ GE S X，CHEN H Y，LI Y R.Real time updating of hydrodynamic model by using Kalman filter[J].Journal of Hydraulic Engineering，2005，36（6）：687693.（in Chinese））

[9]張小峰，穆錦斌，袁晶.一維非恒定水動力學(xué)模型的實時洪水預(yù)報[J].水動力學(xué)研究與進展，2005，20（3）：400404.（ZHANG X F，MU J B，YUAN J.Realtime forecasting method based on 1 D unsteady channel flow[J].Journal of Hydrodynamics.2005，20（3）：400404.（in Chinese））

[10]劉海嬌，苑希民，楊敏，等.基于二維水動力學(xué)模型的洪水演進數(shù)值模擬[J].南水北調(diào)與水利科技，2013（s1）：8991.（LIU H J，YUAN X M，YANG M.Numerical modeling of flood routing based on twodimensional hydrodynamic model[J].South to North Water Transfers and Water Science & Technology.2013（s1）：8991.（in Chinese））

[11]袁曉輝，張雙全，張勇傳，等.非線性馬斯京根模型參數(shù)率定的新方法[J].水利學(xué)報，2001（5）：7781.（YUAN X H，ZhANG S Q，ZHANG Y C.Parameter estimation of nonlinear Muskingum model using mixed genetic algorithm[J].Journal of Hydraulic Engineering，2001（5）：7781.（in Chinese））

[12]孫美云，劉俊，左君，等.分段馬斯京根法在黃河龍潼河段洪水預(yù)報中應(yīng)用[J].水電能源科學(xué)，2012（8）：5052.（SUN M Y，LIU J，ZUO J.Application of Muskingum segmentation algorithm in flood forecasting in LongTong reach of the Yellow River[J].Water Resources and Power.2012（8）：5052.（in Chinese））

[13]黃國如，胡和平，尹大凱.馬斯京根康吉洪水演算方法的穩(wěn)定性分析[J].水科學(xué)進展，2001，12（2）：206209.（HUANG G R，HU H P，YI D K.Stability condition analysis of MuskingumCunge flood routing method[J].Advances in Water Science.2001，12（2）：206209.（in Chinese））

[14]張碩輔.水庫洪水演算方法探討[J].水利水電技術(shù)，1991（3）：49.（ZHANG S F.Discussion of reservoir flood routing method[J].Water Resources and Hydropower Engineering，1991（3）：49.（in Chinese））

[15]許拯民.馬斯京根法應(yīng)用研討[J].華北水利水電大學(xué)學(xué)報（自然科學(xué)版），2005，26（4）：13.（XU Z M.Application of Muskingum flood rooting method[J].Journal of North China Institute of Water Conservancy and Hydroelectric Power.2005，26（4）：13.（in Chinese））

[16] 趙人俊.馬斯京根法——河道洪水演算的線性有限差解[J].河海大學(xué)學(xué)報（自然科學(xué)版），1979（1）：4658.（ZHAO R J.Muskingum method Linear finite difference solution of river flood routing[J].Journal of Hohai University （Natural Sciences），1979（1）：4658.（in Chinese））

[17] 張旺，劉攀，李澤君，等.基于EDA與GA算法的水庫群短期防洪優(yōu)化調(diào)度研究[J].水資源研究，2015，04（5）：395403.（ZHANG W，LIU P，LI Z J，et al.The study of reservoir group shortterm flood prevention optimal regulation based on EDA and GA algorithms[J].Water resources research，2015，04（5）：395403.（in Chinese））

[18] 黃平.最優(yōu)化理論與方法[M].北京：清華大學(xué)出版社，2009.（HUANG P.Theory and methods of optimization[M].Beijing：Tsinghua University Press，2009.（in Chinese））