段煥豐　 俞海寧　 俞國平

（同濟大學(xué) 環(huán)境科學(xué)與工程學(xué)院，上海 200092）

摘要：對城市用水需求量的預(yù)測，是配水系統(tǒng)為其操作運行系統(tǒng)準備和執(zhí)行計劃在線控制的主要任務(wù)。本文運用了時間序列法和BP神經(jīng)網(wǎng)絡(luò)法分別對配水系統(tǒng)的短期和中長期用水量進行預(yù)測，并引入了遺傳算法對這兩種預(yù)測方法進行了優(yōu)化改進。從算例結(jié)果比較來看，優(yōu)化后的預(yù)測方法具有良好的有效性和實用性。

關(guān)鍵詞：遺傳算法　 時間序列法　 BP神經(jīng)網(wǎng)絡(luò)　 用水量預(yù)測

Study on Improved Forecasting Models of Urban Water Demands

Duan Huan-feng，Yu Hai-ning，Yu Guo-ping
(School of Environmental Science and Engineering, Tongji University, Shanghai 200092)

Abstract：Forecasting for urban water demands is the main task for on-line control of a water distribution system to prepare and execute a plan for operating the system. This paper forecasts the short-term and medium or long-term water demands by Time-series Analysis and BP Neural Networks respectively, and introduces Genetic Algorithms to improve and optimize these two forecasting methods. The cases study shows that the improved methods can obtain favorable validity and practicability.
Key words： Genetic Algorithms；Time-series Analysis；BP Neural Networks；Forecasting for Water Demands

　　對配水系統(tǒng)未來用水量的預(yù)測是進行水規(guī)劃的一項重要的任務(wù)。一般來說，用水量預(yù)測分為兩種：一種是短期預(yù)測，以滿足操作系統(tǒng)在線實時模擬和給水系統(tǒng)的優(yōu)化調(diào)度；另一種是對水量的中長期預(yù)測，用以進行水資源規(guī)劃和城市整體設(shè)計規(guī)劃。目前，用以預(yù)測的方法很多，科學(xué)合理地選擇預(yù)測模型應(yīng)首先考慮預(yù)測的時間間隔，不用時間間隔應(yīng)選用不同的預(yù)測模型。根據(jù)目前多年來各種預(yù)測積累的實踐經(jīng)驗和各種預(yù)測方法的特點，對短期預(yù)測一般用到時間序列法，對中長期預(yù)測一般用到BP神經(jīng)網(wǎng)絡(luò)法或灰色模型法（本文只就前者討論，后者可用同樣方法討論）。但是，預(yù)測與實際必然存在一定的誤差，如何減少誤差，以使得預(yù)測結(jié)果盡量準確，是目前急需解決難題。本文針對傳統(tǒng)的預(yù)測方法，引入遺傳算法對其進行參數(shù)優(yōu)化，得到了兩種改進預(yù)測方法。

1、短期預(yù)測優(yōu)化模型

1.1時間序列預(yù)測法
　　時間序列法是根據(jù)已有的時間序列數(shù)據(jù)，通過統(tǒng)計分析的方法建立和已有數(shù)據(jù)盡量擬合的數(shù)學(xué)模型，如自回歸模型AR(p)序列、滑動平均模型MA(q)序列以及回歸與滑動平均混合模型ARMA(p,q)序列，并將得到數(shù)學(xué)模型根據(jù)遞推預(yù)報法對未來的變化趨勢進行預(yù)測。其一般的預(yù)測流程步驟如下：
　　第一步：取樣。就是通過觀測得到的樣本序列Q₁,…, Q_n。
　　第二步：數(shù)據(jù)的預(yù)處理。
　?。?）　 數(shù)據(jù)序列平穩(wěn)性檢驗
　　（2）　 均值歸零處理。即作變換
　　　　 ，
　　其中，
　　第三步：計算出{W_i}的樣本自相關(guān)函數(shù)，樣本自協(xié)方差函數(shù)和樣本偏相關(guān)函數(shù)。
　　第四步：模型識別。利用和來判斷模型的類別（AR(p) 、MA(q)、 ARMA(p,q)）和階數(shù)（p和q值）。
　　第五步：模型的參數(shù)估計。根據(jù)第四步確立的模型類別和階數(shù)，應(yīng)用統(tǒng)計分析的方法估計出模型的參數(shù)。
　　第六步：寫出模型方程。先寫出適合Wt的模型方程，再用代入到模型方程中去。
　　第七步：用水量預(yù)報。根據(jù)已有的模型方程對未來的用水量進行預(yù)測。
1.2改進時間序列法
　　改進時間序列法就是對用遺傳算法上述預(yù)測中第五步估計得到的參數(shù)進行優(yōu)化，使得得到的模型方程更準確，預(yù)測結(jié)果更有效。其優(yōu)化的思想方法為：利用得到的參數(shù)構(gòu)造變化區(qū)間，以時間序列模型中的離散白噪聲的均方誤差為目標函數(shù)和各參數(shù)為變量，再利用遺傳算法的搜索優(yōu)化性能，在區(qū)間范圍內(nèi)找出更佳的參數(shù)解。
　　優(yōu)化目標函數(shù)：

　　其中，E(ξ_t²)為離散白噪聲均方誤差，φ_j和θ_j為參數(shù)變量。

　　其優(yōu)化步驟如下：
　　步驟1：構(gòu)造變量區(qū)間。設(shè)有上述時間序列法得到的某一參數(shù)變量為c_j，則它的變化區(qū)間可構(gòu)造為[a_j，b_j]，其中：
　　a_j =c_j-d，b_j =c_j+d.　 且d為一正常數(shù)。
　　步驟2：參數(shù)編碼。采用二進制編碼，設(shè)編碼長度為e=10，群體規(guī)模為200。
　　步驟3：生成初始父代個體即初始解并評價適應(yīng)度。若適應(yīng)度滿足要求，則結(jié)束。否則，繼續(xù)。
　　步驟4：選擇（采用比例選擇法）、雜交（pc=1.0）、變異（pm=0.10）操作。
　　步驟5：進化迭代。有步驟5得到的子代個體作為新的父代個體，轉(zhuǎn)步驟3。
　　步驟6：加速循環(huán)。把第一次和第二次進化迭代產(chǎn)生的優(yōu)秀個體的參數(shù)變化范圍作為新的變化區(qū)間，算法轉(zhuǎn)入步驟2，如此往復(fù)循環(huán)，直到給定的加速次數(shù)。
1.3應(yīng)用
　　以國內(nèi)某市1991年6月份某段時間的時用水量為例，分別用時間序列法和改進時間序列法來預(yù)測未來24小時的時用水量（由于篇幅限制，原始數(shù)據(jù)省略），其得到預(yù)測模型方程分別為：

　　時間序列法模型：X_t=7966-0.17X_t-1+0.25X_t-2+ξ_t　　

　　改進時間序列法模型：X_t=7983-0.175X_t-1+0.257X_t-2+ξ_t　　 　

　　預(yù)測結(jié)果對比如下表1：

表1 預(yù)測結(jié)果及其誤差比較（部分）

　　從上表可以看出，改進時間序列預(yù)測方法比一般的時間序列方法得到的結(jié)果更加準確有效，同時也看出，這類的預(yù)測方法在預(yù)測短時期的用水量是比較準確，誤差較小，而越往后誤差越大。

2、中長期預(yù)測優(yōu)化模型

　　對城市用水量的中長期預(yù)測，BP網(wǎng)絡(luò)算法使用較為普遍，并且從實踐來看，結(jié)果相對較為理想。BP神經(jīng)網(wǎng)絡(luò)是以樣本訓(xùn)練學(xué)習(xí)而不是用程序指令來完成某一特定的任務(wù)，一般由輸入層、隱層和輸出層三層拓撲結(jié)構(gòu)組成，用來解決模式識別、預(yù)測及擬合等問題。利用BP神經(jīng)網(wǎng)絡(luò)解決用水量預(yù)測問題的模型如下：
2.1 BP神經(jīng)網(wǎng)絡(luò)預(yù)測用水量模型
　　該模型以已有的實際用水量序列以及氣溫、季節(jié)、節(jié)假日等影響因素作為輸入量，通過學(xué)習(xí)訓(xùn)練得到預(yù)測模式。具體如下：

　　其中，Q_I為最終預(yù)測用水量，Q_j為實測系統(tǒng)漏失量，

　　 T_j⁵為s層第j天輸出的預(yù)測用水量修正值系數(shù)，

　　π_j⁵為s層第j天輸出的預(yù)測用水量，

　　 且

　　　　 Q_rj為實測系統(tǒng)漏失量閾值，

　　W_j⁵為氣溫及節(jié)假日等影響因素權(quán)重值，

S.t.　　　　　

　 其中，E為最終輸出系統(tǒng)誤差，_ε為系統(tǒng)誤差允許值，

　　　　 e(k)為第k個樣本系統(tǒng)誤差，n為樣本對數(shù)，

　　　 且：

　　t_I^k為第k個樣本最終輸出期望值，m為網(wǎng)絡(luò)節(jié)點數(shù)， Q_I為第k個樣本最終預(yù)測用水量。

　　BP神經(jīng)網(wǎng)絡(luò)算法一般程序流程圖如下：

2.2 優(yōu)化改進模型
　　如同改進時間序列模型一樣，引入遺傳算法對BP網(wǎng)絡(luò)的參數(shù)優(yōu)化，這些參數(shù)包括各個權(quán)值w_ij和各個閾值θ_j，使得網(wǎng)絡(luò)全局誤差函數(shù)極小化。其參數(shù)優(yōu)化步驟同1.2中步驟1~6。其優(yōu)化目標函數(shù)為：

　 其中，n為樣本對數(shù)，m為網(wǎng)絡(luò)節(jié)點數(shù)，Q_j為實測量，Q_ij為預(yù)測量。

2.3 應(yīng)用
　　下面采用國外某城市1991年上半年部分日用水量數(shù)據(jù)，來預(yù)測該年7月份的用水量。預(yù)測的BP網(wǎng)絡(luò)結(jié)構(gòu)模型中，網(wǎng)絡(luò)節(jié)點數(shù)m=40，隱層數(shù)s=1，節(jié)點作用函數(shù)采用Sigmiod函數(shù)yi=1/(1+e^-x_i)，并設(shè)最大循環(huán)次數(shù)為5000次，允許系統(tǒng)誤差設(shè)為0.01。分別用BP網(wǎng)絡(luò)和改進BP網(wǎng)絡(luò)來預(yù)測7月份的用水量（由于篇幅限制，原始數(shù)據(jù)省略），其得到預(yù)測結(jié)果如下表2，并且改進前后輸出的全局網(wǎng)絡(luò)誤差分別為0.0081和0.0019。

表2 預(yù)測結(jié)果及誤差比較（部分）

從上表的預(yù)測結(jié)果可以看出，利用BP神經(jīng)網(wǎng)絡(luò)法對中長期用水量的預(yù)測較為合適，并且優(yōu)化改進后的BP網(wǎng)絡(luò)預(yù)測更加準確有效。

3、結(jié)語

通過上述實例仿真對比結(jié)果可以看出，通過遺傳算法優(yōu)化改進的時間序列預(yù)測模型和BP神經(jīng)網(wǎng)絡(luò)預(yù)測模型在預(yù)測精度上有很大的提高，使得預(yù)測結(jié)果更加準確，也證明了該方法的有效性和實用性；而且，這種方法也同樣可以推廣應(yīng)用到其它的預(yù)測模型中（如灰色理論預(yù)測模型）。至于遺傳算法作為優(yōu)化方法，其自身的參數(shù)確定具有一定的經(jīng)驗性，仍需進一步的探究。

參考文獻：

[1] 金菊良，丁晶. 遺傳算法及其在水科學(xué)中的應(yīng)用. 成都，四川大學(xué)出版社，2000，68~88
[2] 周明，孫樹棟. 遺傳算法原理及應(yīng)用. 北京，國防工業(yè)出版社，1999，78~89
[3] 袁一星，張杰等. 城市用水量中長期預(yù)測模型的研究. 給水排水，2004，30（6）：102~105
[4] 俞峰，陶建科. 城市給水系統(tǒng)時用水量預(yù)測方法的探討.城市公用事業(yè)，2003，17(6)：23~26

作者簡介：段煥豐，男，同濟大學(xué)環(huán)境學(xué)院，021－65987044 ，duanhuanfeng@sohu.com

通訊地址：上海市同濟大學(xué)環(huán)境科學(xué)與工程學(xué)院　　　　 郵編：200092