徐洪福,袁一星,蘭宏娟
（哈爾濱工業(yè)大學(xué)給水排水系統(tǒng)研究室）

　　摘要: 周期性指數(shù)平滑法是時間序列分析方法中的一種。以天津市給水管網(wǎng)系統(tǒng)為例，進(jìn)行日用水量預(yù)測。通過對預(yù)測結(jié)果進(jìn)行殘差分析，說明這種日用水量預(yù)測模型是可行的。
　　關(guān)鍵詞: 日用水量預(yù)測 時間序列分析 周期性指數(shù)平滑法 水平因子 趨勢因子 周期因子

Forcasting Daily Water Demands of Period Index Smooth Method
Xu Hongfu,Yuan Yixing,Wei Baocheng
(Workgroup of Water & Wastewater System, Harbin Institute of Technology)
　　Abstract: Period index smooth method is an approach of time-series analyses. The study applies this method in the daily demand forecast of Tianjin distribution networks. The error analysis of the forecast result demonstrates that this model is effective.
　　Keywords: daily demand forecast; time-series analysis; period index smooth method; level factor; tendency factor; period factor.

1． 引言

　　城市用水量預(yù)測在城市建設(shè)規(guī)劃﹑供水系統(tǒng)的優(yōu)化調(diào)度中具有重要的作用，它的準(zhǔn)確程度直接影響到供水系統(tǒng)調(diào)度決策的可靠性及實用性。
　　用水量預(yù)測按預(yù)測方法可分為兩類：回歸分析方法和時間序列分析方法。回歸分析也稱為解釋性預(yù)測，它通過觀測系統(tǒng)輸出的結(jié)果，并對系統(tǒng)的輸入數(shù)據(jù)加以分析﹑比較﹑來尋找兩者是否存在因果關(guān)系，回歸模型對長期預(yù)測來說是一種有效的方法；時間序列分析與回歸分析有明顯的不同之處，它把系統(tǒng)看成一個“暗盒”，只依賴于歷史觀測數(shù)據(jù)模式，用適當(dāng)?shù)臄?shù)理統(tǒng)計方法對這個時間序列加以解釋，確定它的數(shù)據(jù)模式，然后選用預(yù)測模型進(jìn)行預(yù)測。

2． 周期性指數(shù)平滑法

　　許多時間序列的變化情況和季節(jié)因素有關(guān)，或者說呈現(xiàn)周期性的變化規(guī)律。有的不僅含有周期性的變化因素而且還有線性增長或減少的趨勢。對于這種時間序列，前幾種指數(shù)平滑法基本上是無效的，而要采取更為高級的指數(shù)平滑方法。它的基本原理是要把這種時間序列分解成三個部分，第一部分是水平因子；第二部分是趨勢因子1；第三部分是周期因子。先把這三部分從時間序列中分離出來，然后再合起來進(jìn)行預(yù)測。這就存在一個周期長度的問題。周期長度要通過自相關(guān)分析來確定。
　　在應(yīng)用周期性指數(shù)平滑法進(jìn)行預(yù)測時，必須事先獲得前兩個周期內(nèi)每一時期的觀測值。設(shè)時間序列的周期長度，其計算過程如下，共分為九個步驟。
　　第一步：
　　分別計算前兩個周期每期的平均數(shù)。
　　第一個周期的平均數(shù)V₁：V1=(x₁+x₂+∧+x_l)/l······（1）
　　第二個周期的平均數(shù)V₂：V2=(x_1+l+x_2+l+∧+x_2l)/l······（2）
　　第二步：
　　計算兩個周期內(nèi)平均每個時期的增量B：B=(V₂-V₁)/l······ （3）
　　第三步：
　　計算初始指數(shù)平滑值S：S= V₂+(l-1)B/2······（4）
　　第四步：
　　分別計算前兩個周期內(nèi)每一時期的季節(jié)因子
　　第一個周期內(nèi)每一個時期的季節(jié)因子C_t＇：

······（5）

　　其中：
　　當(dāng)t=l時，m=1；當(dāng)t=2時，m=2；.....當(dāng)t=l時，m=l。
　　第二個周期內(nèi)每一個時期的季節(jié)因子C_t＇：

······（6）

　　其中：
　　當(dāng)t=l+1時，m=1；當(dāng)t=l+2時，m=2；當(dāng)t=2l時，m=l。
　　第五步：
　　計算前兩個周期中平均每個時期的季節(jié)因子C_t”：

　　C_t”=(C_t-l＇+C_t＇)/2 ······（7）
　　第六步：
　　將季節(jié)因子正態(tài)化
　　按第五步計算出來的l個平均季節(jié)因子之和可能不等于周期l，因此需要按比例縮小或擴(kuò)大，使它們之和等于l。
　　先計算這l個平均季節(jié)因子之和l＇：

　　l＇=C_l+1＂+C_l+2＂+∧+C_2l＂······（8）

　　再計算正態(tài)化以后的季節(jié)因子C_t：

　　C_t=l· C_t＂/l＇ ······（9）
　　共計算出l個正態(tài)化以后的季節(jié)因子C_t，它們之和必然為l。
　　第七步：
　　對第三個周期內(nèi)每一個時期做初步預(yù)測

　　F_t+m=(S+mB)C_{t-l+m······（10）}

　　第八步：
　　當(dāng)?shù)谌齻€周期的第一個時期的觀測值得到(x_t=x_2l+1)，就可以用一組確定的平滑常數(shù)α,β,γ的數(shù)值來修正指數(shù)平滑值、趨勢和季節(jié)因子，修正公式如下：

　　這樣可以重新預(yù)測第三個周期內(nèi)其余(l-1)個時期的數(shù)值
　　F_t+m=(S_t+mB)C_t-l+m ······（12）

　　修正以后的指數(shù)平滑預(yù)測值比修正前的預(yù)測值要準(zhǔn)確。
　　第九步：
　　以后，每次當(dāng)獲得前t時期的觀測值x_t時，就可以用以下公式分別計算單指數(shù)平滑值、趨勢和季節(jié)因子：

······（13）

　　對(t+m)時期的預(yù)測值
　　F_t+m=(St+mBt)C_{t-l+m······（14）}每當(dāng)計算完一個周期，得到l個季節(jié)因子以后，就要按第六步的方法，把它們重新加以正態(tài)化。

3． 工程應(yīng)用

　　本文利用上述周期性平滑指數(shù)法對天津市日用水量建模預(yù)測。對1999年7月7日~1999年9月28日天津市日用水量記錄資料可由圖 和表 中看到。經(jīng)自相關(guān)分析，天津市日用水量序列自相關(guān)系數(shù)列于表1中，對應(yīng)變化曲線見圖1。從表1和圖1中可以看出此時間序列表現(xiàn)出很強的周期性，所以非常適于建立周期性指數(shù)平滑模型。本文作者經(jīng)過實際調(diào)查得出，對于每一個星期，休息日與工作日用水量存在顯著的差異，是導(dǎo)致其自相關(guān)系數(shù)出現(xiàn)周期性的直接原因。而且，從圖1中我們也可以發(fā)現(xiàn)，自相關(guān)系數(shù)每隔六個出現(xiàn)一次“峰值”，所以可以選定此序列的周期為7，即L=7。

利用上述算法編制計算機程序，搜索得到一組最佳平滑常數(shù)為α=0.71；β=0.53；γ=0.03，。所以

　　預(yù)測模型為：

　　利用上述模型對天津市日用水量進(jìn)行預(yù)測，預(yù)測結(jié)果如圖2所示。預(yù)測殘差序列變化曲線如圖3所示。

4．預(yù)測精度分析

　　殘差序列中最大殘差值為-121.54千噸，最大百分比誤差為-9.4%，對于工程精度來說，可以接受。對殘差序列進(jìn)行自相關(guān)檢驗，取前30個自相關(guān)系數(shù)，如表2所示，對應(yīng)的殘差自相關(guān)系數(shù)分布曲線如圖4所示。由表2和圖4可知，前30個自相關(guān)系數(shù)全部落在95%置信區(qū)間[-0.234，+0.234]內(nèi)，也就是說有95%的置信度認(rèn)為所有的自相關(guān)系數(shù)與零沒有顯著差異。由此可以斷定這個殘差序列是一個隨機時間序列。
　　對前30個自相關(guān)系數(shù)求Q值，得Q=13.61，與附表中df=m-1=29，置信度為95%時所對應(yīng)x₂₉²的值，即x₂₉²=42.55比較Q29²，，同樣說明有95%置信度認(rèn)為這30個自相關(guān)系數(shù)中與零沒有顯著性差異，預(yù)測殘差序列為隨機性誤差。所以用周期性平滑指數(shù)法預(yù)測日用水量不存在必然性誤差，預(yù)測模型是成功的。

　　參考文獻(xiàn):
　　[1] 趙洪賓等 城市供水管網(wǎng)系統(tǒng)節(jié)能技術(shù)課題技術(shù)報告，哈爾濱建筑大學(xué)。
　　[2] 徐洪福 城市用水量預(yù)測理論與方法的研究，碩士學(xué)位論文。
　　[3] 王勇領(lǐng) 《預(yù)測計算方法》，科學(xué)出版社。
　　[4] 楊位欽等 《時間序列分析與動態(tài)數(shù)據(jù)建?！?，北京理工大學(xué)出版社，1988。
　　[5] 姚庭寶等 《Turbo Pascal 7.0程序設(shè)計及Turbo Vision 使用大全》，電子工業(yè)出版社。

作者簡介：
　　徐洪福 博士 哈爾濱工業(yè)大學(xué)市政環(huán)境工程學(xué)院
　　通 訊 處： 哈爾濱工業(yè)大學(xué)（二區(qū)）624#信箱 150090
　　0451-6282281
　　Email: xxhf@0451.com