用混合式遺傳算法進行給水管網(wǎng)現(xiàn)狀分析

1、用混合式遺傳算法進行給水管網(wǎng)現(xiàn)狀分析    摘要：將廣義簡約梯度法（GRG）思想應(yīng)用于傳統(tǒng)遺傳算法，提出一種新的混合式遺傳算法（HGA），此方法具有很強的全局和局部搜索能力，并且無需復(fù)雜的編碼、解碼過程，可以準(zhǔn)確解決多變量、多峰值的給水管網(wǎng)現(xiàn)狀分析問題。采用此算法，絕大多數(shù)的計算值與實測值的誤差可減小到1以下，能準(zhǔn)確反映管網(wǎng)的實際工況。 關(guān)鍵詞：混合式遺傳算法 廣義簡約梯度法 遺傳算法  The Use of Hybrid Genetic Algorithmin Analyzing Status of Water Distribution Ne

2、tworksAbstract: The philosophy of General Reduced Gradient (GRG) method is applied to conventional GeneticAlgorithm (GA) and a new Hybrid Genetic Algorithm (HGA) is proposed. Tl1is Il1ethod has a powerful searchingability globally and locally, with which the issues in analyzing the current statlls o

3、f water distribution networksfeaturing more variables and more peak values can be solved accurately without any complicated coding and de-coding processes. With this method, the differences hetween most calculated values and actually-measured valuescan be reduced to below l% and the actuaI working c

4、onditions of water distribution networks can be reflected ac-curately.Key words: water distribution; network; Hyhrid Genetic Algorithm (HGA); General Reduced Gradient(GRG); Genetic Algorithm (GA) 隨著用水規(guī)模的不斷擴大及管道使用年限的增長，城市給水管網(wǎng)逐漸暴露出一些問題，例如各水源間不合理調(diào)度引起管網(wǎng)壓力過高，造成能量浪費，甚至引起爆管問題；管道淤積結(jié)垢，閥門未正常開啟使管道摩阻增大，管網(wǎng)壓

5、力局部降低；管道銹蝕，或接口處漏水造成水資源浪費，管網(wǎng)水量。水壓不足等問題。為了解這些情況，及時、準(zhǔn)確掌握管網(wǎng)工況，必須進行管網(wǎng)現(xiàn)狀分析，并且其結(jié)果可為管網(wǎng)優(yōu)化調(diào)度、改建、擴建提供可靠依據(jù)。1 給水管網(wǎng)現(xiàn)狀分析常用方法簡介給水管網(wǎng)現(xiàn)狀分析是在水源的供水壓力和供水量已知，通過部分節(jié)點、管段的水壓、流量的實際量測值來推斷節(jié)點流量、管線摩阻及所有的節(jié)點水壓和管段流量。常用方法有現(xiàn)狀平差法、遺傳算法和廣義簡約梯度法。 現(xiàn)狀平差法1是通過實測與經(jīng)驗得出節(jié)點流量和管線摩阻，然后轉(zhuǎn)化為管網(wǎng)平差問題求解節(jié)點水壓與管段流量，再與實測值相比較得知管網(wǎng)工況。由于節(jié)點流量和管段摩阻均很難準(zhǔn)確確定，此方法常導(dǎo)致平差計算

6、結(jié)果與實際值存在較大偏差，以致現(xiàn)狀分析得出錯誤的結(jié)論。王榮和等2提出通過以節(jié)點流量和管線摩阻力控制變量的非線性規(guī)劃來進行管網(wǎng)現(xiàn)狀分析，并以遺傳算法求解，得到較好的效果。但由于遺傳算法中隨機抽取樣本，存在偶然性，有時對部分節(jié)點或管段提供的結(jié)論并不準(zhǔn)確，只能達到80的準(zhǔn)確率，而且對管網(wǎng)末梢的樹狀管段和節(jié)點，也不能正確控制。由于遺傳算法局部收斂能力較差，對于大、中型管網(wǎng)，計算時間很長，且很難逼近最優(yōu)解。而文獻3則嘗試采用解決非線性規(guī)劃的有效方法廣義簡約梯度法來求解，其基本思想是利用簡約梯度構(gòu)造一個使目標(biāo)函數(shù)改善的可行方向，然后沿此方向進行搜索，找出一個更優(yōu)點，從而逐步逼近最優(yōu)解。此方法具有很強的局部

7、收斂能力，但全局搜索能力欠佳，對于多峰值的規(guī)劃問題，很容易陷入局部最優(yōu)解。本文將提出一種集廣義簡約梯度法和遺傳算法于一體的混合式遺傳算法，針對現(xiàn)狀分析問題目標(biāo)函數(shù)。約束條件均可微且較易求得的特點，在遺傳算法全局搜索的基礎(chǔ)上加以沿負(fù)梯度方向的變異，以提高其局部搜索能力，很快達到全局最優(yōu)解。2 數(shù)學(xué)模型的建立依據(jù)在滿足管網(wǎng)水力條件約束下，在允許的調(diào)幅范圍內(nèi)，通過對節(jié)點流量q，管段過水能力系數(shù)r進行調(diào)整，使測壓點水壓和測流管段流量的計算值與實測值之差降至最小的原則，建立數(shù)學(xué)模型3。3 數(shù)學(xué)模型的求解1，4-53.1 初始群體的產(chǎn)生樣本染色體為 MN維實向量，由 M個節(jié)點流量和N個管段過水能力系數(shù)組成

8、。初始種群的產(chǎn)生是在各自的界限值范圍內(nèi)隨機選取n組。以下的交叉。變異過程將直接采用控制變量的參數(shù)值，而無需編碼。解碼。另外，本方法樣本規(guī)模n可適當(dāng)減小，一般可取傳統(tǒng)遺傳算法的一半左右。3.2 樣本染色體交叉對于目標(biāo)函數(shù)值較小的樣本，則將其作為父代染色體按算術(shù)組合進行交叉得到子代染色體。父代染色體向量從第k代染色體向量集的子集中隨機抽取產(chǎn)生。設(shè)分別為1（k）和2（k），則子代染色體3（k+1）=×1（k）+（1-）×2（k）。其中 可由式= Z1（ Z1 Z2）確定（其中 Z為染色體向量所對應(yīng)的目標(biāo)函數(shù)值）。3.3 樣本染色體變異對于目標(biāo)函數(shù)值較大的樣本，則需進行變異。其中目

9、標(biāo)函數(shù)值超過一定限度的，按照初始群體產(chǎn)生的方法隨機選取新的個體，以確保樣本的多樣性，避免陷入局部收斂，稱之為隨機變異。另外一部分則沿目標(biāo)函數(shù)的負(fù)梯度方向變異，即：q(k+1)=q(k)+(k)·Pq(k)，r(k+1)=r(k)+(k)·Pr(k)其中(k)為第k步的步長，Pq(k)，Pr(k)為q，r的搜索方向，即目標(biāo)函數(shù)的負(fù)梯度方向。3.4 交叉、變異概率的動態(tài)調(diào)整在迭代之初，樣本以隨機變異為主，目的是進行充分的全局搜索以達到最優(yōu)解的鄰域，并加以少量的交叉和沿負(fù)梯度方向變異操作。在以后的過程中將逐步加強交叉和沿負(fù)梯度方向變異操作，以加快收斂速度，逐步逼近最優(yōu)解。本文采用

10、一次函數(shù)來動態(tài)確定交叉和兩種變異的概率。3.5 終止法則依經(jīng)驗采用最大迭代步數(shù)GENMAX，即迭代步數(shù)達到規(guī)定值時即停止計算。也可采用限制精度來終止計算，即當(dāng)q(k)-時終止計算。4 算例驗證筆者分別采用廣義簡約梯度法（GRG）、傳統(tǒng)遺傳算法（GA）以及前述混合式遺傳算法（HGA）進行了軟件編制，并對如圖1簡單管網(wǎng)進行計算驗證。對于該管網(wǎng)，假設(shè)實測值如表1所示。根據(jù)此值可計算出各節(jié)點流量和管段過水能力系數(shù)的準(zhǔn)確值，如表2所示。表1 假設(shè)實測值計算參數(shù)計算點實測值節(jié)點水壓/m節(jié)點164.1696節(jié)點262.2248管段流量/（L·s-1）管段10.1587管段20.0587管段30.0413表2 本文方法計算結(jié)果計算參數(shù)計算點準(zhǔn)確值計算值誤差/%節(jié)點流量/（L·s-1）節(jié)點10.10000.10000.00節(jié)點20.10000.10000.00管段過水能力系數(shù)管段10.04620.04610.22管段20.04100.04181.95管段30.01090.01090.00由圖2、圖3可以看出，利用傳統(tǒng)遺傳算法計算盡管可以在很短時間內(nèi)收斂到最優(yōu)解的鄰域，但之后卻很難達到最優(yōu)解，迭代200次后目標(biāo)函數(shù)值為4.7419。而本文所提出的混合式遺傳算法則具有很強的局部搜索能力，迭代20次后目