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

水利工程論文-序列二次規(guī)劃法在多水源管網(wǎng)優(yōu)化調(diào)度中的應用研究摘要：供水管網(wǎng)優(yōu)化調(diào)度的一級優(yōu)化是一個非線性優(yōu)化問題，本文通過分析管網(wǎng)的水力關系，對管網(wǎng)水力關系進行合理的線性化，使目標函數(shù)和約束條件顯式化，將問題轉(zhuǎn)化為序列二次規(guī)劃問題。在求解二次規(guī)劃問題中，考慮到大部分節(jié)點水頭的約束是非作用約束，利用線性化結(jié)果，將非作用約束從約束集中剔除，同時將齒行法的思想和水力學上的基本概念相結(jié)合，提出了一種適合本問題的修正齒行法，將二次規(guī)劃結(jié)果拉回到原約束面，保證了解的可行性。最后還初步分析了優(yōu)化計算的計算量。數(shù)值試驗表明本文的方法計算量小、效率高，結(jié)果可靠。關鍵詞：多水源管網(wǎng)優(yōu)化序列二次規(guī)劃齒行法1問題的提出城市供水管網(wǎng)是城市的生命線之一。這一復雜的網(wǎng)絡系統(tǒng)，主要通過幾個供水泵站為城市血液提供能量，送至城市的各個角落。通過對供水泵站的優(yōu)化調(diào)度，可以降低企業(yè)的制水成本、使管網(wǎng)的供水壓力分布更合理，根據(jù)初步估計，對于一個日供水量為10萬噸的自來水公司，如果供水揚程降低1m，每年可以節(jié)電15萬kWh；由于管道系統(tǒng)的滲漏與水頭有關，降低供水水頭也可以在一定程度上減少管網(wǎng)的滲漏；供水水頭的降低還可以減少爆管的風險，這對于管網(wǎng)的管理有更深刻的意義。因此管網(wǎng)合理調(diào)度研究一直是供水企業(yè)一個重要課題，同時也是一個難題。供水管網(wǎng)運行的合理調(diào)度可以用一個最優(yōu)化問題來描述，管網(wǎng)的水力方程組是一組非線性方程，各水源水泵的開啟狀態(tài)作為離散變量，因此這是一個混合變量的非線性最優(yōu)化問題。由于離散變量與連續(xù)變量的同時存在，求解極為不方便，最常用的方法是將該問題分作兩級進行優(yōu)化：一級優(yōu)化是針對管網(wǎng)而言，目的在于求各水源的最佳供水量或最佳供水揚程；二級優(yōu)化是在一級最優(yōu)化的基礎上，根據(jù)水源的具體情況，確定滿意的水泵開啟方案和水泵的調(diào)速比。采用以上方法可以在一定程度上降低求解的困難，但一級優(yōu)化也是一個非線性的優(yōu)化問題，求解起來相當麻煩，目前國內(nèi)外最常用的方法是廣義簡約梯度法。廣義簡約梯度法雖屬較優(yōu)秀的約束非線性規(guī)劃算法，根據(jù)作者在以往其他優(yōu)化應用方面的研究，其重分析次數(shù)相當多。在本優(yōu)化問題中水力計算是計算量的主體部分，由數(shù)值試驗的經(jīng)驗知，在目前中等配置的微機上完成一個2000個左右節(jié)點的供水管網(wǎng)，一次水力計算需要10s左右，如果采用廣義簡約梯度法，需要反復迭代計算，花費的時間是相當可觀的。由此可見采用廣義簡約梯度法實現(xiàn)管網(wǎng)的在線優(yōu)化調(diào)度存在較大難度。本文針對一級優(yōu)化問題，采用序列二次規(guī)劃法進行求解。M.J.D.Powell所給出的序列二次規(guī)劃法實質(zhì)上是運用KuhnTucker最優(yōu)化條件所形成的非線性方程進行迭代計算，而這一迭代過程恰好可以用求解一相應的二次規(guī)劃問題替代，故原問題的求解過程轉(zhuǎn)化為求解一個二次規(guī)劃的序列。其中二次規(guī)劃問題的二次目標函數(shù)是原問題Lagrange函數(shù)的二次展開式，包含了目標與約束函數(shù)的二次信息。通常其二階導數(shù)矩陣由變尺度的思想通過先前迭代點的梯度信息逐步生成。序列二次規(guī)劃法綜合利用了K一種比廣義簡約梯度法優(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ù)，與水源的供水量、供水水頭有關，據(jù)此可以建立供水管網(wǎng)的目標函數(shù)：minFG(Qs,Hs)(1)式中：FG為各水源的制水成本和供水的動力費用；Qs、Hs為各水源的供水量和供水水頭。供水調(diào)度的主要約束條件有：管網(wǎng)的水力關系，各水源的水量和水壓的約束，管網(wǎng)中各節(jié)點的最小服務水頭。這些約束條件分別表示如下：管網(wǎng)水力關系F(Hs,HN,QN)=0(2)各水源的供水水頭約束HsminHsHsmax(3)各水源的供水量約束Qsmin(Hs)QsQsmax(Hs)(4)管網(wǎng)各節(jié)點服務水頭約束(5)其中：Hsmax、Hsmin分為水源的最大、最小供水水頭；Qsmin(Hs)、Qsmax(Hs)分為水源的最大、最小供水能力，通常水源的供水量的能力與供水水頭有關。HN為管網(wǎng)中各節(jié)點的服務水頭；HNmax、HNmin分為管網(wǎng)中各節(jié)點的最大、最小服務水頭；QN為管網(wǎng)中各節(jié)點的節(jié)點流量。3模型的求解模型求解主要有2個難點：(1)約束條件太多，一個中等復雜的城市管網(wǎng)可能會有上千個約束；(2)目標函數(shù)中各變量隱式相關水源的供水水頭Hs和供水水量Qs隱式相關。如果能對以上兩個方面進行適當?shù)奶幚恚梢源蟠蟮慕档碗y度，提高求解效率。針對以上兩點，本文從管網(wǎng)的水力條件出發(fā)，提出了一套求解方法：在一定負荷N下，將管網(wǎng)的水力計算公式(2)在H0處作一階泰勒展開有：(6),稱為敏度矩陣。如果用哈真-威廉公式表示管道的能量損失，用矩陣A、B可以分別表示為,，管網(wǎng)的水力學公式可以用式(7)表達。A和B僅與管網(wǎng)中管道的水力坡度有關。當任一水源的供水水頭發(fā)生變化，由于管網(wǎng)自身的調(diào)節(jié)作用，每根管道的水力坡度的變化幅度要比節(jié)點水頭變化小得多，A、B的變化都比較小。管網(wǎng)的水力計算公式(2)在H0附近可以線性近似為式(7)，且方程有足夠的精度(算例的數(shù)值計算結(jié)果參見附錄)。(7)在文獻2中已證明B是正定對稱矩陣，其逆矩陣存在。令則(8)矩陣C的分量ci,j反映了第j個水源對節(jié)點i的影響，矩陣C也稱為影響矩陣。如果管網(wǎng)中所有水源的供水水頭同步上升h，即HNh1,h2,hsT，相當于管網(wǎng)的參考水位提高了h。由式(8)知管網(wǎng)中任一點的水頭上升的水位，因此其中矩陣C的行向量的各分量之和必等于1，各管段的水力坡度不變。如果各水源的供水水頭和節(jié)點流量已知，可求得管網(wǎng)中的各節(jié)點的水頭，同樣可以求出各水源的供水量。水源泵站供水的動力能耗可以表示為QsHs/(為水源效率)，供水的動力費用與耗能成正比。水源供水量在H0附近可以線性近似為：Qs=LsHs。水源的制水0費用(除動力費用)可以表示成RsLsHs(Rs表示各水源的單位制水成本)，因此目標函數(shù)在H0處可以近似用水源水頭的二次函數(shù)表示如下：(9)各水源的供水量約束在H0處可以線性近似表示：(10)由于在管網(wǎng)中往往只是一部分的最不利節(jié)點違反約束，只要最不利的節(jié)點滿足了服務水頭的要求，其他節(jié)點也滿足了要求，因此可以將最不利的一些節(jié)點與水源節(jié)點的水頭關系從式(9)中的影響矩陣C中抽取出來，表示成矩陣G，管網(wǎng)節(jié)點水頭的約束方程(5)可以簡化表示如下：(11)其中為最不利節(jié)點的集合。通過上述方法，一級優(yōu)化模型在H0附近可以近似表示為線性約束的二次規(guī)劃問題：(12a)s.t.HsminHsHsmax(12b)Qsmin(H0s)+KsminHsLsHsQsmax(H0s)+KsmaxHs(12c)(12d)在原優(yōu)化問題中，各水源的供水量、管網(wǎng)中節(jié)點的水頭是水源供水水頭的函數(shù)，是隱式關系，求解起來非常不方便。通過把管網(wǎng)水力關系式(2)線性化，消去原目標函數(shù)(1)中的變量水源供水量Qs，可以把目標函數(shù)表示僅含水源水頭變量的形式，將管網(wǎng)中各節(jié)點的水頭HN表示成水源的供水水頭Hs的線性函數(shù)，只取其中最不利一部分作為每次優(yōu)化計算的約束條件，這樣大大地減少了約束條件。如果管網(wǎng)有上千個節(jié)點，只要保證最不利的10%左右節(jié)點滿足服務水頭約束，就能基本上保證每次優(yōu)化計算結(jié)果不會離約束邊界太遠，同時優(yōu)化計算的計算量成倍的減少。由于采用了線性近似的方法簡化約束條件和目標函數(shù)，采用二次規(guī)劃法(QP法)優(yōu)化之后會導致結(jié)果越過實際約束邊界，其中主要是最不利點不滿足管網(wǎng)最小服務水頭的要求。在結(jié)構(gòu)優(yōu)化設計中經(jīng)常采用齒行法進行優(yōu)化迭代，其基本思想是在每次優(yōu)化迭代后，通過射線步(即將所有設計變量以同一倍數(shù)放大或縮小)將結(jié)果拉到最嚴格的約束邊界上。根據(jù)管網(wǎng)水力學，所有的水源的供水水頭同時都提高或降低相同的水位，使管網(wǎng)的最不利點的水位恰好處于約束邊界上，不會改變各個水源的供水關系。利用這一特性，可以構(gòu)造一修正的射線步，能夠方便的將中間優(yōu)化迭代點拉回到約束界面上(見圖1)。由于在這一修正的射線步中，每個水源提高的水位相同，因此變化后的值