基于模型預測控制的泵站群聯合調度研究

張少愷，龍 巖，管 一，李 樹

（1.河北工程大學水利水電學院，河北 邯鄲 056038；2.河北工程大學河北省智慧水利重點實驗室，河北 邯鄲 056038；3.通遼市水利事業發展中心，內蒙古 通遼 028000）

受全球氣候變暖的影響和城市化的快速發展，由極端惡劣天氣引發的城市洪澇災害愈發嚴重，受災范圍愈發廣泛[1]。泵站群作為福州市排澇的主要工程措施，在汛期發揮重大作用[2]。本文以福州市汛期承擔主要排水功能的晉安河流域為研究對象，從安全、經濟層面考慮構建泵站群聯合防洪優化調度模型，同時引入模型預測控制（MPC：Model Predictive Control）策略，通過優化算法求解，提出福州市晉安河流域汛期調度方案。

1 研究區域

1.1 流域狀況

晉安河發源于福州市北峰，橫跨福州市晉安區南北，是福州市最長的城市內河。晉安河全長6.7 km，流域面積69.7 km2。晉安河下游為光明港，最后注入閩江。

1.2 水利工程現狀

晉安河上游為琴亭湖，湖水日常水位保持在4.50 m 左右。琴亭湖周邊地面高程均在7.80 m 以上，當湖水位在7.80 m 以下時，周邊不會產生內澇。琴亭湖泵站水泵數量為5 臺，總裝機3000 kW，總設計流量25 m3/s，平均設計流量5 m3/s，設計揚程0～7.0 m。光明港下游的魁岐排澇站承擔主要排水任務，水泵數量為3 臺，總裝機6600 kW，總設計流量120 m3/s，平均設計流量40 m3/s，設計凈揚程2.6 m。琴亭湖、魁岐泵站群調度概化，如圖1所示。

圖1 琴亭湖、魁岐泵站群調度概化

2 模型預測控制策略

為了更好地解決模型優化過程中的失配、畸變、干擾等因素引起的不確定性[3]，引入基于MPC 的滾動優化策略。MPC 是一種滾動優化控制算法，主要包括模型預測、滾動優化、反饋校正3 個環節。MPC滾動優化策略，如圖2所示。

圖2 MPC滾動優化策略示意

圖2中，Np為預測時域，Nc為控制時域。在t時，利用預測時域內的預測數據，在控制時域內進行多目標優化，得到控制時域內最優的控制序列；但在當前t時段，只執行第一個控制序列的指令。令t=t+1，重復進行上一過程，直至完成整個時間軸的優化。

3 河道水動力模型構建

為實現對晉安河的洪澇過程模擬，需獲得琴亭湖上游來水和琴亭湖至光明港的區間來水，建立基于SWMM 的晉安河河道水力學模型。SWMM 模型，如圖3所示。

圖3 SWMM 模型

通過查閱相關文獻，以三八泵站和魁岐排澇站為例，模擬2021年福州市的2 次降雨，采用Nash-Sutcliffe 效率系數法，對SWMM 模型參數進行率定與驗證。模擬水位過程線與實測水位過程線形狀比較吻合，水位峰值及峰值出現時間的偏差均較小。三八泵站和魁岐排澇站6月29日降雨事件率定結果如圖4—5所示，其納什效率系數分別為0.53、0.92；8月8 日降雨事件驗證結果如圖6—7 所示，其納什效率系數分別為0.68、0.55。

圖4 三八泵站率定結果

圖5 魁岐排澇站率定結果

圖6 三八泵站驗證結果

圖7 魁岐排澇站驗證結果

4 基于模型預測控制的泵站群調度模型構建

4.1 預測模型

本文在泵站群優化調度模型內建立基于水量平衡原理的預測模型，選擇時間間隔為1 h 進行采樣，對未來8 h的水位進行預測。預測模型為：

式中：Vi,t為i水庫在t時段的庫容值（m3）；Vi,t+1為i水庫在t+ 1 時段的庫容值（m3）；Ii,t為i水庫在t時段的入庫流量（m3/s）；Qi,t為i水庫在t時段的出庫流量（m3/s）；Vi,0為i水庫初始庫容值（m3）；b0i為i水庫初始水位值（m）；Vi,T為i水庫T時段末庫容值（m3）；bT i為i水庫T時段末水位值（m）;t為調度時段( )

h ,t=1,2,…,T;T為調度總時段( )h ，取8 h；ΔT表示單位時間長度為1 h；i=1代表琴亭湖,i=2代表光明港區域。

4.2 梯級泵站優化調度模型

4.2.1 目標函數

（1）調度優化目標1。琴亭湖超警風險最小。其目標函數為：

式中：F1為調度周期內琴亭湖各時段距死水位的最小值（m）；Zx為琴亭湖汛限水位（m）；Z1,t為t時段琴亭湖水位（m）；Zs為琴亭湖死水位（m）；t、T含義及取值同上。

（2）調度優化目標2。光明港峰值水位最低。其目標函數為：

式中：F2為調度周期內光明港各時段峰值水位的最小值（m）；Z2,t為t時段光明港的水位（m）；t的含義及取值同上。

（3）調度優化目標3。調度運行期間抽水費用最少。其目標函數為：

式中：F3為調度周期內總抽水費用（元）；k為電費單價（元/kW·h）；N為泵站數量，n=1 代表琴亭湖泵站，n=2 代表魁岐排澇站；ρ為水的密度（g/cm3）；g為重力加速度（m/s2），取9.8 m/s2；Hn,t為第n個泵站在t時段的揚程（m）；Qn,t為第n個泵站在t時段的過泵流量（m3/s）；η為水泵損失效率;ΔT、t、T含義及取值同上。

（4）調度優化目標4。閘門啟閉次數最少。泵站的檢修次數和檢修成本與泵站機組在使用過程中的啟閉有關。閘門啟閉次數的目標函數為：

式中：F4為調度期間內泵站機組啟閉次數之和；L(n,t)為n泵站的t時段與上一時段泵站機組相比變化的臺數；n= 1 代表琴亭湖泵站，n= 2 代表魁岐排澇站；t、T含義及取值同上。

（5）調度優化目標5。琴亭湖調度時段末水位最低。其目標函數為：

式中：F5為琴亭湖調度周期內各時段末最低水位（m）；Z1,t為t時段琴亭湖水位（m）；t的含義及取值同上。

（6）調度優化目標6。光明港調度時段末水位最低。其目標函數為：

式中：F6為光明港調度周期內各時段末最低水位（m）；Z2,t為t時段光明港水位（m）；t的含義及取值同上。

4.2.2 約束條件

（1）水位約束。其約束條件公式為：

式中：Zn,t為第n個泵站在t時段的上游水位( m );、分別為第n個泵站的上游水位上、下限約束（m），為0，為7.0 m，為0，為5.0 m；n= 1 代表琴亭湖泵站，n= 2 代表魁岐排澇站；t的含義及取值同上。

（2）流量約束。其約束條件公式為：

（3）水泵揚程約束。其約束條件公式為：

4.2.3 模型求解

由于琴亭湖、魁岐泵站群聯合優化調度是多目標優化問題，所以采用構造懲罰函數的方法求解。將琴亭湖調度時段末水位低于目標水位、光明港峰值水位最低和調度時段末水位低于目標水位、泵站啟閉次數最少作為懲罰函數，融入抽排水費用最低目標函數中，基于JAVA 語言通過Eclipse 軟件編寫優化調度程序，運用人工蜂群算法計算。

4.3 MPC設計

本文將模型預測應用于泵站群優化調度模型中，在已知當前時刻琴亭湖、光明港的水位、水位庫容曲線和未來一段時間的入庫流量前提下，進行晉安河流域泵站群未來最優調控，使得琴亭湖、光明港充分發揮調蓄能力。MPC優化調度方案步驟，如圖8所示。

圖8 MPC優化調度方案步驟

結合圖8，將MPC 應用于對未來時段晉安河防洪調度的最優控制決策，具體步驟如下。

步驟1：首先進行人工蜂群算法的相關參數設置。①初始種群規模：本文種群規模取值為200。種群規模SN包含采蜜蜂和觀察蜂，各取種群數量的1/2。②一般食物源為種群的1/2。③設置limit限制，超過此限制沒有更新的采蜜蜂變為觀察蜂。④設置maxcycle為100。

步驟2：數據輸入。設定采樣時間為1 h 和預測步長為8 h，總調度時長為24 h。輸入k時段未來8 h的琴亭湖上游來水、光明港區間來水。

步驟3：在已知k時刻琴亭湖、光明港初始水位并已知未來8 h入庫流量的前提下，將晉安河泵站群起始數據代入預測模型。

步驟4：利用人工蜂群算法計算懲罰函數，求解出使性能指標最小的最優控制序列，輸出最優控制序列中的首個控制指令進行晉安河泵站群優化調度。

步驟5：將MPC 防洪優化調度最優控制序列結果代入河道動力模型中，獲取k+ 1 時刻初始水位，進行水位反饋校正。

步驟6：進入下一時刻最優控制序列的求解。

5 實例分析

選取100 a 一遇設計降雨序列，作為雨量邊界，輸入河道動力模型獲取琴亭湖、光明港入流過程。琴亭湖主要考慮上游來水，光明港主要考慮琴亭湖至光明港的區間來水以及魁岐排澇站泵、閘排水方式的影響。當福州城區晉安河水位低于閩江潮位時，將沿江的閘門全部關閉，若超過警戒水位，此時開啟魁岐排澇站泵站機組，降低晉安河水位；當福州城區內河水位高于閩江潮位時，將魁岐泵站機組關閉，僅采用魁岐水閘排水，將閘內水位降到警戒水位以下。琴亭湖、光明港區間來水過程如圖9所示，閩江潮位如圖10所示。

圖9 琴亭湖、光明港區間來水

圖10 閩江潮位過程

在100 a一遇的設計洪水情況下，結合所建水力學模型和實際防洪排澇需求，得出以下結論。

（1）從經濟運行角度分析，MPC 優化調度運行費用為2.36萬元，規則調度運行費用為2.48萬元，運行費下降了4.8%；從設備安全及維護角度分析，MPC優化調度的泵站啟閉次數為43，規則調度的泵站啟閉次數為46。規則調度、MPC 優化調度泵站啟閉情況對比，詳見表1。

表1 規則調度、MPC優化調度泵站啟閉情況對比

（2）從防洪效果上看，MPC 優化調度方案琴亭湖、光明港峰值水位分別為6.15、5.621 m，規則調度方案琴亭湖、光明港峰值水位分別為7.11、5.732 m，洪峰水位分別降低了13.5%、1.93%。琴亭湖、光明港100 a一遇調度，如圖11—12所示。

圖11 琴亭湖100 a一遇調度

圖12 光明港100 a一遇調度

（3）從實際工程角度出發，因MPC 優化調度采用滾動策略實時計算，可以應對突發情況，相比其他依靠既定規則的優化調度，實用性大幅提升。

6 結語

本文基于實際工程需求，構建了基于模型預測控制的泵站群聯合調度模型，為福州市江北城區應對100 a一遇的典型洪水場景提供了優化調度方案，可作為福州市汛期實時排澇調度的參考依據。