序列二次規(guī)劃法在多水源管網(wǎng)優(yōu)化調(diào)度中-應用研究

1、序列二次規(guī)劃法在多水源管網(wǎng)優(yōu)化調(diào)度中-應用研究摘要：供水管網(wǎng)優(yōu)化調(diào)度的一級優(yōu)化是一個非線性優(yōu)化問題，本文通過分析管網(wǎng)的水力關(guān)系，對管網(wǎng)水力關(guān)系進行合理的線性化，使目標函數(shù)和約束條件顯式化，將問題轉(zhuǎn)化為序列二次規(guī)劃問題。在求解二次規(guī)劃問題中，考慮到大部分節(jié)點水頭的約束是非作用約束，利用線性化結(jié)果，將非作用約束從約束集中剔除，同時將齒行法的思想和水力學上的基本概念相結(jié)合，提出了一種適合本問題的修正齒行法，將二次規(guī)劃結(jié)果拉回到原約束面，保證了解的可行性。最后還初步分析了優(yōu)化計算的計算量。數(shù)值試驗表明本文的方法計算量小、效率高，結(jié)果可靠。 關(guān)鍵詞：多水源管網(wǎng) 優(yōu)化 序列二次規(guī)劃 齒行法 1 問題的提出

2、 城市供水管網(wǎng)是城市的生命線之一。這一復雜的網(wǎng)絡系統(tǒng)，主要通過幾個供水泵站為城市血液提供能量，送至城市的各個角落。通過對供水泵站的優(yōu)化調(diào)度，可以降低企業(yè)的制水成本、使管網(wǎng)的供水壓力分布更合理，根據(jù)初步估計，對于一個日供水量為10萬噸的自來水公司，如果供水揚程降低1m，每年可以節(jié)電15萬kWh；由于管道系統(tǒng)的滲漏與水頭有關(guān)，降低供水水頭也可以在一定程度上減少管網(wǎng)的滲漏；供水水頭的降低還可以減少爆管的風險，這對于管網(wǎng)的管理有更深刻的意義。因此管網(wǎng)合理調(diào)度研究一直是供水企業(yè)一個重要課題，同時也是一個難題。 供水管網(wǎng)運行的合理調(diào)度可以用一個最優(yōu)化問題來描述，管網(wǎng)的水力方程組是一組非線性方程，各水源水泵

3、的開啟狀態(tài)作為離散變量，因此這是一個混合變量的非線性最優(yōu)化問題。由于離散變量與連續(xù)變量的同時存在，求解極為不方便，最常用的方法是將該問題分作兩級進行優(yōu)化：一級優(yōu)化是針對管網(wǎng)而言，目的在于求各水源的最佳供水量或最佳供水揚程；二級優(yōu)化是在一級最優(yōu)化的基礎(chǔ)上，根據(jù)水源的具體情況，確定滿意的水泵開啟方案和水泵的調(diào)速比。采用以上方法可以在一定程度上降低求解的困難，但一級優(yōu)化也是一個非線性的優(yōu)化問題，求解起來相當麻煩，目前國內(nèi)外最常用的方法是廣義簡約梯度法。廣義簡約梯度法雖屬較優(yōu)秀的約束非線性規(guī)劃算法，根據(jù)作者在以往其他優(yōu)化應用方面的研究，其重分析次數(shù)相當多。在本優(yōu)化問題中水力計算是計算量的主體部分，由數(shù)

4、值試驗的經(jīng)驗知，在目前中等配置的微機上完成一個2000個左右節(jié)點的供水管網(wǎng)，一次水力計算需要10s左右，如果采用廣義簡約梯度法，需要反復迭代計算，花費的時間是相當可觀的。由此可見采用廣義簡約梯度法實現(xiàn)管網(wǎng)的在線優(yōu)化調(diào)度存在較大難度。Tucker最優(yōu)化條件所形成的非線性方程進行迭代計算，而這一迭代過程恰好可以用求解一相應的二次規(guī)劃問題替代，故原問題的求解過程轉(zhuǎn)化為求解一個二次規(guī)劃的序列。其中二次規(guī)劃問題的二次目標函數(shù)是原問題Lagrange函數(shù)的二次展開式，包含了目標與約束函數(shù)的二次信息。通常其二階導數(shù)矩陣由變尺度的思想通過先前迭代點的梯度信息逐步生成。序列二次規(guī)劃法綜合利用了K條件、變尺度、線

5、性及二次近似等有效手段，在理論上是一種比廣義簡約梯度法優(yōu)秀的算法1，但是它的迭代序列通常從不可行域逐步逼近可行域，需要在極限情況下才能完全達到約束要求，這顯然不利于盡快獲得可行的較優(yōu)解，故約束條件的妥善處理非常重要，本文將結(jié)構(gòu)優(yōu)化中齒行法的思想和水力學的基本概念相結(jié)合，提出了一種新的算法，可以方便地將迭代中的非可行點拉回到約束界面上，獲得了較高的計算效率，有助于實現(xiàn)管網(wǎng)的在線優(yōu)化調(diào)度。2 供水優(yōu)化調(diào)度一級優(yōu)化的數(shù)學模型 管網(wǎng)的運行調(diào)度一般以經(jīng)濟性作為目標函數(shù)，與水源的供水量、供水水頭有關(guān)，據(jù)此可以建立供水管網(wǎng)的目標函數(shù)：minFG(Qs,Hs)(1)式中：FG為各水源的制水成本和供水的動力費用

6、；Qs、Hs為各水源的供水量和供水水頭。 供水調(diào)度的主要約束條件有：管網(wǎng)的水力關(guān)系，各水源的水量和水壓的約束，管網(wǎng)中各節(jié)點的最小服務水頭。這些約束條件分別表示如下： 管網(wǎng)水力關(guān)系F(Hs,HN,QN) =0 (2)2007-04-23 各水源的供水水頭約束HsminHsHsmax (3) 各水源的供水量約束Qsmin(Hs)QsQsmax(Hs) (4) 管網(wǎng)各節(jié)點服務水頭約束(6)稱為敏度矩陣。如果用哈真-威廉公式表示管道的能量損失，用矩陣A、B可以分別表示為 , ，管網(wǎng)的水力學公式可以用式(7)表 達。A和B僅與管網(wǎng)中管道的水力坡度有關(guān)。當任一水源的供水水頭發(fā)生變化，由于管網(wǎng)自身的調(diào)節(jié)作用

7、，每根管道的水力坡度的變化幅度要比節(jié)點水頭變化小得多，2007-04-23A、B的變化都比較小。管網(wǎng)的水力計算公式(2)在H0附近可以線性近似為式(7)，且方程有足夠的精度(算例的數(shù)值計算結(jié)果參見附錄)。則，因此其中矩陣C的行向量的各分量之和必等于1，各管段的水力坡度不變。如果各水源的供水水頭和節(jié)點流量已知，可求得管網(wǎng)中的各節(jié)點的水頭，同樣可以求出各水源的供水量。水源泵站供水的動力能耗可以表示為QsHs/(為水源效率)，供水的動力費用與耗能成正比。水源供水量在H0附近可以線性近似為：Qs=LsHs。水源的制水0費用(除動力費用)可以表示成RsLsHs(Rs表示各水源的單位制水成本)，因此目標函

8、數(shù)在H0處可以近似用水源水頭的二次函數(shù)表示如下： (10) 由于在管網(wǎng)中往往只是一部分的最不利節(jié)點違反約束，只要最不利的節(jié)點滿足了服務水頭的要求，其他節(jié)點也滿足了要求，因此可以將最不利的一些節(jié)點與水源節(jié)點的水頭關(guān)系從式(9)中的影響矩陣C中抽取出來，表示成矩陣G2007-04-23，管網(wǎng)節(jié)點水頭的約束方程(5)可以簡化表示如下： (12a)s.t . HsminHsHsmax (12b)Qsmin(H0s) +KsminHsLsHsQsmax(H0s) +KsmaxHs (12c)4 算 例 2007-04-23 為了檢驗上述模型的可行性和可靠性，本文采用Fortran語言編制了優(yōu)化計算程序，

9、算例中管網(wǎng)基本形狀如圖3所示。管網(wǎng)有兩個水源，12個節(jié)點，19根管道，管網(wǎng)節(jié)點的基本數(shù)據(jù)見表1。 在計算中，不失一般性，以泵站供水的能耗作為目標函數(shù)，即 ，與理論分析有一點差別，主要是由計算中截斷誤差引起。 同時本文還比較了計算了在不同的比例負荷作用下，(即各節(jié)點的流量同時乘以一個相同的比例系數(shù))，此時敏度矩陣C矩陣仍然保持基本不變，矩陣假設基本合理。向量E的各系數(shù)與比例系數(shù)的1852次方成正比，數(shù)字計算中間結(jié)果見附表2。 附表 不同狀態(tài)下C、E矩陣的變化情況狀態(tài)1水源1供水水頭：100.227m水源2供水水頭： 97.227m狀態(tài)2水源2007-04-231供水水頭： 98.402m水源2供

10、水水頭：101.131m123456789101112Ci,10.904480.903260.668420.040950.945000.895360.667750.371150.918210.865290.732650.66456Ci,20.035440.096660.331500.958970.054910.104560.332170.628770.081710.134630.267260.33534Ei-4.3374-9.9901-22.4775-2.0184-6.8753-11.7309-24.0387-17.4159-11.4742-16.0033-26.5912-27.2211Ci,

11、10.966630.909430.694080.041450.948330.901980.693260.377060.923070.873920.075140.69799Ci,20.033290.090480.377060.958370.051580.097930.306630.622930.076840.125980.248470.30189Ei-4.3379-9.9913-22.4826-2.0183-6.8760-11.7323-24.0439-17.4166-11.4752-16.0051-26.5951-27.2264附表2 不同比例負荷狀態(tài)下C、E矩陣的變化情況 狀態(tài)1水源1供水水頭：100.227m，水源2供水水頭：97.227m，負荷比1/Q0=1.2，且(Q1/Q0)1.8521.4017狀態(tài)2水源1供水水頭：100.227m，水源2供水水頭：97.227m，負荷比2/Q0=1.5，且(Q2/Q0)1.8522.1189123456789101112Ci,10.964820.904250.