張雅君, 劉全勝
(北京建筑工程學(xué)院城市建設(shè)工程系，北京100044)

　　摘 要：從灰色模型的特點(diǎn)與需水量本身的變化規(guī)律出發(fā)，揭示了GM(1,1)模型的應(yīng)用條件及其在需水量預(yù)測(cè)中的應(yīng)用范圍，以東京都需水量為例進(jìn)行了驗(yàn)證，并對(duì)北京市需水量進(jìn)行了灰色預(yù)測(cè)。?
　　關(guān)鍵詞：需水量；灰色模型；預(yù)測(cè)?
　　中圖分類號(hào)：TU991.31
　　文獻(xiàn)標(biāo)識(shí)碼：C
　　文章編號(hào)：1000-4602(2002)03-0079-03

　　灰色系統(tǒng)理論認(rèn)為，盡管客觀事物或系統(tǒng)表象復(fù)雜、數(shù)據(jù)離亂，但它總是有整體功能和有序的，具有某種內(nèi)在規(guī)律，關(guān)鍵在于怎樣用適當(dāng)?shù)姆椒ㄈネ诰蚝屠盟?］?；疑A(yù)測(cè)模型(GM)通過對(duì)原始數(shù)據(jù)進(jìn)行生成處理，使其呈指數(shù)趨勢(shì)變化，建立指數(shù)微分方程，最終得到預(yù)測(cè)模型。?
　　GM的一般表達(dá)式為GM(h,n)，即對(duì)n個(gè)變量用h階方程表示?；诨疑到y(tǒng)的特點(diǎn)，大多數(shù)預(yù)測(cè)模型采用GM(1,1)及其改進(jìn)模型。
　　由于初始時(shí)序列是否滿足非負(fù)的齊次指數(shù)變化規(guī)律是能否用GM建模的關(guān)鍵，因此建模條件應(yīng)是初始時(shí)序列滿足指數(shù)變化規(guī)律。從數(shù)學(xué)角度看，非負(fù)的指數(shù)序列實(shí)質(zhì)上就是相對(duì)增長率恒定的序列，所以判斷一個(gè)序列是否適用于GM建模的依據(jù)也應(yīng)從序列的相對(duì)增長率入手，基于此提出一種新的相對(duì)增長率判定準(zhǔn)則。 　

1 　理論依據(jù)及應(yīng)用條件

　　灰色模型是一個(gè)指數(shù)模型，這也就決定了它應(yīng)該用在按指數(shù)變化的系統(tǒng)。所以，探討一個(gè)系統(tǒng)是否可用灰色預(yù)測(cè)的方法，實(shí)質(zhì)上是討論該系統(tǒng)是否按指數(shù)規(guī)律變化(同樣可用判定)。?
　　在其他因素不變的條件下，城市需水量增長與人口增長呈線性關(guān)系，而人口的增長方式在一定時(shí)期內(nèi)符合指數(shù)增長規(guī)律^［2］。此外，經(jīng)濟(jì)的發(fā)展在一定時(shí)期內(nèi)表現(xiàn)為等相對(duì)增長率的指數(shù)增長規(guī)律(如目前我國的經(jīng)濟(jì)增長率維持在7%～8%左右)。在不受水資源短缺制約的條件下，經(jīng)濟(jì)的增長必然帶動(dòng)需水量的增長。而用水系統(tǒng)本質(zhì)上并沒有一階正反饋系統(tǒng)特性，需水量增長速率與該時(shí)期的用水量并沒有邏輯上的關(guān)系。但是，基于需水量與人口、經(jīng)濟(jì)之間的密切關(guān)系，在需水量增長不受水資源約束時(shí)，人口、經(jīng)濟(jì)的指數(shù)增長特性將傳遞給需水量，此階段用GM建?？傻玫捷^好的預(yù)測(cè)結(jié)果。?
　　因水資源短缺，使得工業(yè)和生活用水受到約束，各國政府都采取了節(jié)水措施，此時(shí)用水量的增長不再隨人口、經(jīng)濟(jì)等比例增長。一方面經(jīng)濟(jì)、人口的增長使需水量呈上升趨勢(shì)；另一方面由于用水政策、結(jié)構(gòu)調(diào)整、工藝改革等諸多原因使需水量出現(xiàn)下降的趨勢(shì)。所以，此時(shí)的用水量呈波動(dòng)性上升趨勢(shì)(值偏高)，不宜用GM建模預(yù)測(cè)。?
　　未來某些國家可能出現(xiàn)人口下降趨勢(shì)(特別是中國、日本等人口密度高的國家)，而且在以節(jié)約資源為特點(diǎn)的可持續(xù)發(fā)展模式中，產(chǎn)業(yè)結(jié)構(gòu)的調(diào)整、科技的進(jìn)步也勢(shì)必引起需水量的下降。需水量的負(fù)增長在這一時(shí)期內(nèi)也以等相對(duì)增長率的指數(shù)方式減少，因此用GM建模也能得到較好的預(yù)測(cè)結(jié)果。?
　　以日本東京都需水量情況為例，GM模型的應(yīng)用條件得到了很好的驗(yàn)證。東京都水道局負(fù)擔(dān)23個(gè)區(qū)及多摩地區(qū)的供水，東京都1960年—1999年的人口增長情況如圖1所示。其中，20世紀(jì)60年代初到70年代中呈現(xiàn)出快速增長的勢(shì)頭，這一時(shí)期經(jīng)濟(jì)增長率相對(duì)較高，人口、經(jīng)濟(jì)的指數(shù)增長特性傳遞給需水量，因此它呈現(xiàn)出指數(shù)增長趨勢(shì)，用GM建?？傻玫捷^好的預(yù)測(cè)結(jié)果。

　　以1978年—1986年的用水量數(shù)據(jù)建立灰色模型，結(jié)果見表1。?

表1 東京都1978年—1986年用水量數(shù)據(jù)擬合結(jié)果 年份 實(shí)際值(10⁴m³/a) 擬合值(104m3/a) 相對(duì)誤差(%) 平均相對(duì)誤差(%) C值 1978 1776 1776.0 0.00 0.90 0.5719 1979 1728 1706.3 1.25 1980 1665 1704.5 -2.37 1981 1694 1702.6 -0.51 1982 1706 1700.8 0.31 1983 1728 1698.9 1.68 1984 1710 1697.1 0.75 1985 1690 1695.3 -0.31 1986 1678 1693.4 -0.92 　注： 實(shí)際值＞擬合值時(shí)，相對(duì)誤差為正值；實(shí)際值＜擬合值時(shí)，相對(duì)誤差為負(fù)值。

　　從表1的C值結(jié)果看，也證實(shí)了用水量波動(dòng)上升期間不宜采用GM建模。由人口緩慢增長、傳統(tǒng)家庭形式的變遷、經(jīng)濟(jì)緩慢增長等現(xiàn)象推測(cè)，未來較長時(shí)期內(nèi)東京都需水量增加緩慢，其預(yù)測(cè)模型如下：

　　

2 　北京市需水量的灰色預(yù)測(cè)

　　以北京市1989年—1999年的用水量數(shù)據(jù)(見表2)建立GM模型并進(jìn)行預(yù)測(cè)。

表2 北京市用水量數(shù)據(jù)　10⁴m³/a? 年份 工業(yè)用水量 生活用水量 總用水量 1989 454105.3 51428.7 96834 1990 40939.8 52852.8 93792.6 1991 38327 56619.5 94946.5 1992 41575 58376.6 99951.6 1993 38251.7 62571.9 100823.6 1994 41431.2 67295.2 108726.4 1995 39437.3 70540.5 109977.8 1996 39824 69860.5 109684.5 1997 41622.1 75573.5 117195.6 1998 35644.8 78731.9 114376.7 1999 37354.7 82715.5 120070.2

　　工業(yè)、生活、總用水量的相對(duì)增長率k及其方差S的計(jì)算結(jié)果如表3所示。

表3 ?k及S計(jì)算結(jié)果 項(xiàng)目 工業(yè)用水量的k 生活用水量的k 總用水量的k 1 -0.098 0.028 -0.031 2 -0.064 0.071 0.012 3 -0.085 0.031 0.053 4 -0.080 0.072 0.009 5 0.083 0.075 0.078 6 -0.048 0.048 0.012 7 0.010 -0.010 -0.003 8 0.045 0.082 0.068 9 -0.145 0.042 -0.024 10 0.048 0.051 0.050 -0.0164 0.0490 0.0224 s 0.0770 0.0266 0.0359 4.70 0.54 1.604

　　工業(yè)用水量的GM(1,1)模型如下：

　　

　　生活用水量的GM(1,1)模型如下：

　　

　　總用水量的GM(1,1)模型如下：

　　

　　擬合及預(yù)測(cè)結(jié)果如表4所示。

表4 數(shù)據(jù)擬合值及2000年—2002年預(yù)測(cè)結(jié)果 項(xiàng) 目 工業(yè)用水量 生活用水量 總用水量 實(shí)際值(10⁴m³/a) 擬合值(10⁴m³/a) 相對(duì)誤差(%) 實(shí)際值(10⁴m³/a) 擬合值(104m3/a) 相對(duì)誤差(%) 實(shí)際值(10⁴m³/a) 擬合值(104m3/a) 相對(duì)誤差(%) 1989年 45405.3 45405.3 0.00 51428.7 51428.7 0.00 96834.0 96834.0 0.00 1990年 40939.8 40746.1 0.47 52852.8 53920.0 -2.02 93792.6 94252.9 -0.49 1991年 38327.0 40450.4 -5.54 56619.5 56562.8 0.10 94946.5 96817.1 -2.03 1992年 41575.0 40156.8 3.41 58376.6 59335.1 -1.64 99951.6 99562.4 0.40 1993年 38251.7 39865.4 -4.22 62571.9 62243.3 0.53 100823.6 102328.4 -1.49 1994年 41431.2 39576.1 4.48 67295.2 65294.1 2.97 108726.4 105171.1 3.27 1995年 39437.3 39288.9 0.38 70540.5 68494.3 2.90 109977.8 108092.9 1.71 1996年 39824.0 39003.7 2.06 69860.5 71851.5 -2.85 109684.5 111095.8 -1.29 1997年 41622.1 38720.7 6.97 75573.5 75373.2 0.27 117195.6 114182.1 2.57 1998年 35644.8 38439.7 -7.84 78731.9 79067.5 -0.43 114376.7 117354.1 -2.80 1999年 37354.7 38160.7 -2.16 82715.5 85942.9 -0.27 120070.2 120614.3 -0.45 2000年 37883.8 87008.2 123965.1 2001年 37608.8 91272.8 127409.0 2002年 37335.9 95746.4 130948.4 平均相對(duì)誤差（%） 3.412 1.271 1.500 C值 0.5066 0.0602 0.1236

　　由表3、4中數(shù)據(jù)可見，工業(yè)、生活用水量及總用水量中工業(yè)用水量的=4.7為最大(C=0.5066)，預(yù)測(cè)結(jié)果不可靠。究其原因，一方面北京市工業(yè)產(chǎn)業(yè)結(jié)構(gòu)的調(diào)整、節(jié)水措施的加強(qiáng)使工業(yè)用水量呈減少趨勢(shì)；另一方面工業(yè)的高速發(fā)展使其用水仍不斷增加，因此表現(xiàn)出極強(qiáng)的波動(dòng)性，不宜用灰色方法預(yù)測(cè)。生活用水量的=0.54為最小(C=0.0602)，可判定結(jié)果可靠。這是由于生活水平的提高、城市人口增長緩慢，生活用水量仍保持較為平穩(wěn)的指數(shù)增長趨勢(shì)所致?？傆盟吭诠I(yè)、生活用水雙重因素的影響下表現(xiàn)出一定的波動(dòng)性，但用GM進(jìn)行較短周期的預(yù)測(cè)可得到較好的效果(C=0.1236)，若周期過長則可靠度降低。

3 　結(jié)論

　　GM(1,1)模型本身是一個(gè)指數(shù)模型，其預(yù)測(cè)效果在很大程度上取決于原始數(shù)據(jù)的特點(diǎn)。用水量在某些時(shí)期內(nèi)表現(xiàn)出一定的指數(shù)變化規(guī)律，適用于GM(1,1)模型；而在某些時(shí)期則表現(xiàn)出較強(qiáng)的波動(dòng)性，不宜用GM(1,1)模型。
　　此外，一般按指數(shù)規(guī)律變化的狀態(tài)變量的增長常數(shù)隨時(shí)間變化，所以GM(1,1)模型的預(yù)測(cè)周期不宜太長，否則將帶來很大的系統(tǒng)誤差。

參考文獻(xiàn)：

　?。?］聶相田,邱林,朱普生,等.水資源可持續(xù)利用管理不確定性分析方法及應(yīng)用［M］.鄭州:黃河水利出版社，1999.
　?。?］王其藩.系統(tǒng)動(dòng)力學(xué)［M］.北京:清華大學(xué)出版社，1988.

---

　　電　　話：(010)68322517(O)?
　　E-mail:uced@mail.bicea.net.cn
　　收稿日期：2001-07-10?
　　修回日期：2001-08-08