魏光輝

(新疆塔里木河流域管理局，新疆 庫爾勒841000)

合理開發(fā)和利用好地下水，是一個政策問題，也是一個可持續(xù)發(fā)展問題。需水量預(yù)測方法主要分為:時間序列法、結(jié)構(gòu)分析法和系統(tǒng)分析法。幾種最為典型的方法包括:ARMA方法，回歸分析法，指標(biāo)分析法，灰色預(yù)測法，人工神經(jīng)網(wǎng)絡(luò)方法，系統(tǒng)動力學(xué)方法。各種需水量預(yù)測方法都有其自身的優(yōu)點和不足，而需水量預(yù)測就是結(jié)合預(yù)測的目的、特點，結(jié)合用水量變化規(guī)律，合理地選擇一種或幾種預(yù)測方法，并收集所需的數(shù)據(jù)進(jìn)行預(yù)測。

作為灰色系統(tǒng)理論重要內(nèi)容之一的 GM(1，1)模型［1－2］，其應(yīng)用價值在越來越多的領(lǐng)域中得到體現(xiàn)［3－4］。GM(1，1)模型的突出特點是:建模過程簡單，模型表達(dá)式簡潔，便于求解，應(yīng)用廣泛。在發(fā)展系數(shù)a的絕對值較小時(當(dāng)時間間隔很小、序列數(shù)據(jù)變化平緩時，如0(－a≤0.3)，模擬值精度較高。但當(dāng)發(fā)展系數(shù)a的絕對值較大(如－a＞0.5)時，模型偏差較大，無法用于中長期預(yù)測，甚至不宜作短期預(yù)測。一些學(xué)者對此進(jìn)行了改進(jìn)，得到了比原GM(1，1)模型模擬精度高和適應(yīng)性更強(qiáng)的新模型［5－8］。在此基礎(chǔ)上，羅黨等人［9］經(jīng)過分析GM(1，1)模型產(chǎn)生模擬誤差的原因，從構(gòu)造背景值公式入手，優(yōu)化GM(1，1)模型，既保持了原GM(1，1)模型的優(yōu)點，又使優(yōu)化GM(1，1)模型適用于各種發(fā)展系數(shù)的情形，尤其是當(dāng)發(fā)展系數(shù)絕對值較大時也可用于中長期預(yù)測。本文利用優(yōu)化的GM(1，1)模型，對某水庫周邊地區(qū)地下水埋深變化進(jìn)行了預(yù)測，通過數(shù)據(jù)的分析處理、函數(shù)擬合預(yù)測，對于未來水庫周邊地區(qū)地下水資源的可持續(xù)利用提供參考。

1 優(yōu)化GM(1，1)模型建模過程

設(shè)x(0)為非負(fù)準(zhǔn)光滑序列，x(0)={x(0)(1)，x(0)(2)，…，x(0)(n)}

對其進(jìn)行一次累加生成，得到生成序列:

其中當(dāng) x(1)(k)=x(1)(k－1)時，z(1)(k)=z(1)(k－1)

則:

(1)灰色微分方程x(0)(k)+αzx(1)(k)=b的最小二乘估計參數(shù)滿足:

(2)灰色微分方程x(0)(k)+αz(1)(k)=b的白化方程

(3)灰色微分方程x(0)(k)+αzx(1)(k)=b的時間響應(yīng)式為:

(4)還原值為:

(5)灰色預(yù)測檢驗

①殘差檢驗

殘差檢驗就是計算相對誤差，以殘差的大小來判斷模型的好壞。

殘差:Δ(k)=x(0)(k)－0)(k)

GM(1，1)的建模精度:p=1－ε

則ε(k)越小越好，p越大越好，一般要求ε(k)＜20%，p＞80%;最好是 ε(k)＜10%，p＞90%。x(0)為原始數(shù)列，^X(0)是由式(5)、(6)得到的預(yù)測數(shù)據(jù)列。

②后驗差檢驗

式中:S1為原始數(shù)列 x(0)的均方差;S2為殘差序列{Δ(k)}的均方差;C越小，模型越好。

表1 檢驗指標(biāo)等級標(biāo)準(zhǔn)

2 應(yīng)用實例

2.1 模型計算

本文以新疆某平原水庫周邊地區(qū)地下水水位觀測資料為依據(jù)，選擇該水庫2003年至2009年共計7年的地下水水位觀測數(shù)據(jù)，見表2。

表2 某平原水庫周邊地區(qū)歷年平均地下水位

依據(jù)灰色模型原理建立GM(1，1)模型，并求得模型參數(shù):a= －0.098 502，b=2.947 962，最終得模型計算式為:

2.2 模型檢驗

(1)殘差檢驗。經(jīng)檢驗，相對誤差值為－4% ～5%(見表3)。

其殘差平均值ε=2.3% ＜10%，平均精度 p=97.7% ＞95%，模型擬合精度較高，模型判為優(yōu)。

(2)后驗差檢驗。經(jīng)計算，后驗差比值 C=0.1323＜0.35，小誤差概率P=1，模型級別為好。

表3 模型擬合效果分析

2.3 用水量預(yù)測

經(jīng)過檢驗的模型符合精度要求后，可用于外推預(yù)測。2010～2014年該水庫周邊地下水水位預(yù)測結(jié)果見表4。

表4 2010～2014年地下水水位預(yù)測

3 結(jié)論

(1)優(yōu)化GM(1，1)模型建模過程簡單，需要的數(shù)據(jù)也比較少，擬合精度較高，對于地下水埋深預(yù)測，具有很強(qiáng)的適用性。

(2)通過建立數(shù)學(xué)模型來反映地下水的動態(tài)變化規(guī)律，對于今后該地區(qū)的持續(xù)、合理的開采和利用地下水資源起到積極作用，對于農(nóng)業(yè)及水資源的可持續(xù)發(fā)展具有重大意義。

［1］鄧聚龍.灰色系統(tǒng)理論教程［M］.武漢:華中理工大學(xué)出版社，1990.

［2］劉思峰，郭天榜，黨慧國，等.灰色系統(tǒng)理論及其應(yīng)用［M］.北京:科學(xué)出版社，1999.

［3］劉思峰，孫書安.應(yīng)用數(shù)學(xué)［M］.鄭州:河南科技出版社，1993.

［4］劉思峰，鄧聚龍.GM(1，1)模型的適用范圍［J］.系統(tǒng)工程理論與實踐，2000，(5):121 －124.

［5］王義鬧，劉光珍，劉開第.GM(1，1)的一種逐步優(yōu)化直接建模法［J］.系統(tǒng)工程理論與實踐，2000，(9):99 －104.

［6］潭冠軍.GM(1，1)模型的背景值構(gòu)造方法和應(yīng)用［J］.系統(tǒng)工程理論與實踐，2000，(4):98 －103.

［7］呂林正，吳文江.灰色模型 GM(1，1)優(yōu)化探討［J］.系統(tǒng)工程理論與實踐，2001，(8):92 －96.

［8］吉培榮，黃巍松，胡翔勇.灰色預(yù)測模型特性的研究［J］.系統(tǒng)工程理論與實踐，2001，(9):105 －108.

［9］羅黨，劉思峰，黨耀國.灰色模型 GM(1，1)優(yōu)化［J］.中國工程科學(xué)，2003，(8):50 －53.