王 鑫 東

(遼寧省江河流域管理局，沈陽 110003)

0 引 言

當(dāng)前，水資源日趨緊張，對于區(qū)域水資源供需平衡而言，農(nóng)業(yè)灌溉用水占據(jù)較為重要的比重。目前，在最嚴(yán)格的水資源管理制度下，需要對農(nóng)業(yè)灌溉用水進(jìn)行規(guī)劃和設(shè)計，保證水資源的可持續(xù)利用。而合理規(guī)劃農(nóng)業(yè)灌溉用水的前提是需要對灌溉用水進(jìn)行預(yù)測，在預(yù)測的基礎(chǔ)上進(jìn)行合理規(guī)劃。對于農(nóng)業(yè)灌溉用水預(yù)測，主要分為兩種方式，一種是基于長系列的數(shù)據(jù)進(jìn)行統(tǒng)計分析，第二種是采用數(shù)學(xué)模型方法進(jìn)行預(yù)測分析，第一種方法由于區(qū)域灌溉用水?dāng)?shù)據(jù)較少，在區(qū)域農(nóng)業(yè)灌溉用水預(yù)測研究較少。第二種方法由于較為明晰的數(shù)學(xué)原理，被國內(nèi)學(xué)者運用于農(nóng)業(yè)灌溉用水預(yù)測的研究中，并取得一定研究成果[1-5]，但在這些模型中，運用加權(quán)馬爾科夫鏈模型進(jìn)行農(nóng)業(yè)灌溉用水預(yù)測進(jìn)行研究較少，而加權(quán)馬爾科夫鏈模型可考慮中間過程狀態(tài)的影響，在水資源其他領(lǐng)域得到一定程度運用，效果較好[6-8]，為此，本文引入加權(quán)馬爾科夫鏈模型，對農(nóng)業(yè)灌區(qū)灌溉用水進(jìn)行預(yù)測和分析，研究成果可對區(qū)域農(nóng)業(yè)灌溉用水預(yù)測提供參考價值。

1 加權(quán)馬爾科夫鏈模型原理

加權(quán)馬爾科夫鏈先計算一組隨機變量序列{x,t∈T}為隨機生成的一組隨機序列，在隨機序列中：

{T={0,1,2,…,K}

(1)

對于任意時間變量t≥0以及隨機變量的狀態(tài)值j,i0,i1,i2,i3,…,it-1,it∈E。都代表

P={Xt+1=j|X0=i0,X1=i1,…,Xi=it}

(2)

則表示{x,t,∈T}為馬爾科夫鏈。通過馬爾科夫鏈的表達(dá)式可以看出，馬爾科夫鏈表示的系統(tǒng)含義是t+s時刻變量的一個狀態(tài)值只和前個t時刻的狀態(tài)相關(guān)，而和t時刻以前的狀態(tài)無任何關(guān)聯(lián)。

傳統(tǒng)的馬爾科夫鏈未能考慮各個變量因子之間的權(quán)重系數(shù)，在收斂計算時存在收斂異常的情況。為此本文引入改進(jìn)的馬爾科夫鏈，改進(jìn)的馬爾科夫鏈主要是通過建立各個因子不同的轉(zhuǎn)移概率權(quán)重，利用標(biāo)準(zhǔn)化的各因子的自相關(guān)系數(shù)作為轉(zhuǎn)移概率權(quán)重，將變量的某一個時刻的狀態(tài)轉(zhuǎn)移概率進(jìn)行加權(quán)，以加權(quán)后的概率作為變量在某個時刻的預(yù)測值。從而預(yù)測該變量標(biāo)準(zhǔn)化因子所處在的變量值。改進(jìn)的馬爾科夫鏈變量的徑流預(yù)測的主要步驟為：

(1)通過建立各個變量因子的不同等級的標(biāo)準(zhǔn)以及各因子變化轉(zhuǎn)移的概率系數(shù)。采用標(biāo)準(zhǔn)方差法進(jìn)行不同等級的標(biāo)準(zhǔn)化。通過這種方法來確定徑流預(yù)測不同時刻的狀態(tài)值。因為需要計算不同等級下的各因子的變化轉(zhuǎn)移的權(quán)重系數(shù)。需要對各個因子的概率權(quán)重系數(shù)進(jìn)行相關(guān)計算。

(2)各因子不同相關(guān)系數(shù)及標(biāo)準(zhǔn)。通過對各個影響因子進(jìn)行相關(guān)系數(shù)的計算，相關(guān)系數(shù)計算公式為：

(3)

對各個變量因子進(jìn)行標(biāo)準(zhǔn)化處理，標(biāo)準(zhǔn)化處理公式為：

(4)

式中：m代表的物理意義為模型預(yù)測計算的最大因子數(shù)。

(3)結(jié)合前面兩個步驟中不同時刻的計算值作為初始的狀態(tài)值，對應(yīng)各個變量的轉(zhuǎn)移概率加權(quán)權(quán)重系數(shù)，可以推算出下個時刻不同的狀態(tài)變量值，其中K表示的模型計算時段的步長(本文選用年、月)兩個時間步長。

(4)將處在同一個時間變量的狀態(tài)預(yù)測值和各個因子的權(quán)重進(jìn)行疊加計算，可以表示某個變量在該狀態(tài)下的概率預(yù)測值，表達(dá)式為：

(5)

該時刻的計算最大概率值表示預(yù)測時刻對應(yīng)下的狀態(tài)值。

(5)基于模糊數(shù)學(xué)的理論，將各個狀態(tài)變量不同權(quán)重系數(shù)可以得到各因子的權(quán)重矩陣wi，則可以得到年徑流的預(yù)測值，表達(dá)式為：

(6)

式中：T和B分別表示的變量狀態(tài)的上下邊界。

(6)結(jié)合改進(jìn)的馬爾科夫鏈的游歷和方程穩(wěn)定的特點，可以計算不同狀態(tài)的變量概率值。

2 模擬結(jié)果分析

2.1 年灌溉水量預(yù)測結(jié)果對比分析

基于加權(quán)前后的馬爾科夫鏈模型對遼寧某灌區(qū)2000-2010年農(nóng)業(yè)灌溉用水進(jìn)行預(yù)測，并和調(diào)查區(qū)域灌溉水量進(jìn)行對比分析。模擬結(jié)果見表1和圖1。

表1 加權(quán)前后馬爾科夫鏈模型模擬在區(qū)域年農(nóng)業(yè)灌溉用水預(yù)測精度對比Tab.1 The comparative simulation irrigations between the improved Markov chain model and Markov chain model in the year scale

圖1 加權(quán)前后馬爾科夫鏈模型的農(nóng)業(yè)灌溉用水預(yù)測與調(diào)查灌溉水量年相關(guān)圖Fig.1 The compares of the forecast accuracy by Improved Markov chain model and Markov chain model in the year scale

從表1中可以看出，加權(quán)后的馬爾科夫鏈模型模擬的灌區(qū)2000-2010年灌溉用水預(yù)測值和調(diào)查的灌溉用水預(yù)測值相對誤差在-6.74%～9.80%之間，相對誤差的均值為7.83%，而傳統(tǒng)馬爾科夫鏈模型預(yù)測值和調(diào)查的灌溉用水量之間的相對誤差在9.85%～15.73%之間，相對誤差均值為13.6%，可見，加權(quán)后的馬爾科夫鏈模型農(nóng)業(yè)灌溉用水預(yù)測精度較好，明顯好于未加權(quán)的馬爾科夫鏈模型，加權(quán)后的馬爾科夫鏈模型預(yù)測相對誤差小于10%，具有較好的精度。從圖1中可以看出，加權(quán)后的馬爾科夫鏈模型農(nóng)業(yè)灌溉用水值和調(diào)查值具有較好的相關(guān)系數(shù)，相關(guān)系數(shù)達(dá)到0.7以上，綜上，可見，加權(quán)后的馬爾科夫鏈模型可用于區(qū)域農(nóng)業(yè)灌溉用水預(yù)測。

2.2 月農(nóng)業(yè)灌溉水量預(yù)測結(jié)果對比分析

農(nóng)業(yè)灌溉用水在各個月份的值均不同，為此，基于加權(quán)前后的馬爾科夫鏈模型分別對遼寧某灌區(qū)2000-2010年農(nóng)業(yè)灌溉用水進(jìn)行預(yù)測，并和調(diào)查區(qū)域灌溉水量進(jìn)行對比分析。模擬結(jié)果見表2和圖2。

表2 加權(quán)前后馬爾科夫鏈模型模擬在區(qū)域年農(nóng)業(yè)灌溉用水預(yù)測精度對比Tab.2 The comparative simulation irrigations between the improved Markov chain model and Markov chain model in the month scale

圖2 加權(quán)前后馬爾科夫鏈模型的農(nóng)業(yè)灌溉用水預(yù)測與調(diào)查灌溉水量月相關(guān)圖Fig.2 The compares of the forecast accuracy by Improved Markov chain model and Markov chain model in the month scale

從表2中可以看出，加權(quán)后的馬爾科夫鏈模型模擬的灌區(qū)2000-2010年月尺度灌溉用水預(yù)測值和調(diào)查的灌溉用水預(yù)測值相對誤差在-14.02%～24.07%之間，相對誤差的均值為17.9%，而傳統(tǒng)馬爾科夫鏈模型預(yù)測值和調(diào)查的灌溉用水量之間的相對誤差在-22.92%～38.24%之間，相對誤差均值為30.1%，可見，加權(quán)后的馬爾科夫鏈模型月尺度農(nóng)業(yè)灌溉用水預(yù)測精度明顯好于未加權(quán)的馬爾科夫鏈模型，但是在月尺度上，加權(quán)后的馬爾科夫鏈模型預(yù)測相對誤差明顯高于年尺度模擬的相對誤差，這主要是因為區(qū)域灌溉水量主要為調(diào)查數(shù)據(jù)，月尺度調(diào)查數(shù)據(jù)的精度低于年尺度調(diào)查精度，因此預(yù)測相對誤差較大。從圖2中可以看出，在月尺度上，加權(quán)后的馬爾科夫鏈模型農(nóng)業(yè)灌溉用水值和調(diào)查值具有較好的相關(guān)系數(shù)，相關(guān)系數(shù)達(dá)到0.86以上，好于年尺度相關(guān)系數(shù)，這就是因為，加權(quán)后的馬爾科夫鏈模型可考慮中間過程，因此在月尺度上的相關(guān)系數(shù)較好。

3 結(jié) 語

本文應(yīng)用加權(quán)的馬爾科夫鏈模型進(jìn)行農(nóng)業(yè)灌溉用水進(jìn)行預(yù)測，并和傳統(tǒng)馬爾科夫鏈模型預(yù)測結(jié)果進(jìn)行對比，研究取得以下結(jié)論。

(1)加權(quán)的馬爾科夫鏈模型可用于區(qū)域農(nóng)業(yè)灌溉用水預(yù)測，在年尺度和月尺度模擬精度符合規(guī)范要求，預(yù)測精度好于傳統(tǒng)馬爾科夫鏈模型。

(2)加權(quán)的馬爾科夫鏈模型可考慮變量中間過程的狀態(tài)，因此，在月尺度預(yù)測的灌溉用水和調(diào)查灌溉水量相關(guān)系數(shù)好于年尺度預(yù)測的相關(guān)系數(shù)。

□

[1] 萬玉文,蘇 超,方 崇. 灰色預(yù)測模型在廣西達(dá)開水庫灌區(qū)農(nóng)業(yè)灌溉用水管理中的應(yīng)用[J]. 節(jié)水灌溉,2012,(9):58-60,66.

[2] 張 丹. 基于BP人工神經(jīng)網(wǎng)絡(luò)的水稻需水量預(yù)測及灌區(qū)灌溉制度研究[D]. 沈陽：沈陽農(nóng)業(yè)大學(xué),2007.

[3] 唐延芳. 灌溉用水量預(yù)測方法的研究[D]. 沈陽：沈陽農(nóng)業(yè)大學(xué),2007.

[4] 馬勝利. 基于灰色理論與模糊控制的智能灌溉系統(tǒng)研究[D]. 西安：陜西科技大學(xué),2013.

[5] 陳文軒. 多目標(biāo)模糊決策方法在鐵山水庫用水調(diào)度管理中的應(yīng)用研究[D]. 長沙：中南大學(xué),2009.

[6] 茹正亮,楊芝艷,朱文剛,等. 加權(quán)馬爾科夫AR-GARCH-GED模型在降水量中的預(yù)測[J]. 系統(tǒng)工程,2013,(12):98-102.

[7] 秦 松. 基于加權(quán)的滑動平均-馬爾科夫預(yù)測模型及其應(yīng)用[J]. 水資源與水工程學(xué)報,2013,(1):185-188.

[8] 羅 坤,陶麗平,黃銀華. 加權(quán)馬爾科夫的灰色殘差修正模型在酸性降水pH值預(yù)測中的應(yīng)用[J]. 科技資訊,2013,24:73-74.

[9] 錢 會,李培月,王 濤. 基于滑動平均-加權(quán)馬爾科夫鏈的寧夏石嘴山市年降雨量預(yù)測[J]. 華北水利水電學(xué)院學(xué)報,2010,(1):6-9.