平穩(wěn)時間序列模型優(yōu)選及階數(shù)比選研究
——以年最高洪水位預(yù)報為例

宋輝（遼寧廣播電視大學(xué)丹東分校，遼寧 丹東 118000）

?

平穩(wěn)時間序列模型優(yōu)選及階數(shù)比選研究
——以年最高洪水位預(yù)報為例

宋輝
（遼寧廣播電視大學(xué)丹東分校，遼寧 丹東 118000）

摘 要：水位預(yù)報因具有數(shù)據(jù)直觀，復(fù)核方便，容易確定避險區(qū)域等優(yōu)點而逐漸取代復(fù)雜的流量預(yù)報.平穩(wěn)時間序列模型計算簡單，容易操作，是水位預(yù)報的理想方法.本文對比分析自回歸模型、滑動平均模型和自回歸滑動平均模型，以簡捷原則為原則建立了自回歸模型的參數(shù)計算方法和水位預(yù)報模型.并以丹東一水文站25年水位資料為例進(jìn)行了兩個階數(shù)對比，結(jié)果表明3階預(yù)報精度高于5階，比選結(jié)果認(rèn)為一般情況下可只進(jìn)行3階預(yù)報就能滿足精度要求.

關(guān)鍵詞：平穩(wěn)時間序列；自回歸模型；階數(shù)比選；自相關(guān)函數(shù)；水位預(yù)報

1 引言

水資源是生命之源、生產(chǎn)之基、生態(tài)之要.隨著國民經(jīng)濟(jì)的發(fā)展，五位一體戰(zhàn)略的實施，水資源對社會經(jīng)濟(jì)發(fā)展的制約作用逐漸顯現(xiàn).無論是主動的資源節(jié)約制約，還是被動的環(huán)境達(dá)標(biāo)制約，都需要對水資源進(jìn)行準(zhǔn)確的預(yù)測.

隨著水文站網(wǎng)密度的增加，常規(guī)的流量預(yù)測逐漸被高密度的水位預(yù)測所取代.水位預(yù)測在功能上具有數(shù)據(jù)直觀，復(fù)核方便，容易確定避險區(qū)域等優(yōu)點.為提高水位預(yù)報精度，簡化預(yù)報方法，提出簡單、易操作、程序化的水文預(yù)報方法就具有重要的現(xiàn)實意義.

平穩(wěn)時間序列模型計算簡單，過程易懂，應(yīng)用方便，只要有某處的實測資料就可以進(jìn)行預(yù)報［1］.因此，本文在對比分析平穩(wěn)時間序列模型的基礎(chǔ)上，利用自回歸模型（AR），對比分析不同階數(shù)的預(yù)報結(jié)果，以期為水文預(yù)報提供可行性方法.

2 平穩(wěn)時間序列模型

時間序列是指將某種現(xiàn)象某一個統(tǒng)計指標(biāo)在不同時間上的各個數(shù)值，按時間先后順序排列而形成的序列.平穩(wěn)時間序列是指含變量的時間序列，它的前期演變過程的統(tǒng)計相關(guān)規(guī)律在未來的一段時間內(nèi)是不變的.

平穩(wěn)時間序列需滿足均值與時間無關(guān)的常數(shù)，方差是與時間無關(guān)的常數(shù)，協(xié)方差只與時間間隔有關(guān)，與時間無關(guān)的常數(shù).

2.1 自回歸模型（AR）

AR（p）模型（Auto-Regressive Moving Average）是一種平穩(wěn)的自回歸模型.AR模型是一種線性預(yù)測，可由模型推出第N點前面或后面的數(shù)據(jù)（設(shè)推出P點）［2］.

設(shè)X1，X2，…,Xt是平穩(wěn)時間序列，自回歸模型AR（p）：

式中b1,b2,…bp是預(yù)報系數(shù)，也稱回歸系數(shù)，表示序列在時間上的相依程度；P是預(yù)報階數(shù)；Xt,Xt-1,Xt-2,…,Xt-p+1為預(yù)報值Xt+1前1、2、…、p個時間間隔的實測值；εt為外部條件影響值.

2.2 滑動平均模型（MA）

MA（q）模型（moving average model）滑動平均模型，是現(xiàn)代預(yù)報中常用的模型.

設(shè)｛εt｝是一白噪聲序列，時間序列｛Xt｝滿足

上式為q階平穩(wěn)滑動模型，式中εt為外部白噪聲；β1,β2, βq是預(yù)報系數(shù)；q是預(yù)報階數(shù)；Xt為預(yù)報值.

2.3 自回歸滑動平均模型（ARMA）

ARMA（p,q）模型（auto regressive moving average model）自回歸滑動平均模型，集中AR模型法與MA模型法優(yōu)點［2］［3］.

設(shè)｛Xt｝是時間序列，｛εt｝是白噪聲序列，兩序列滿足平穩(wěn)條件，則

式中符號同上文，其中a1,…,ap與b1,…,bp都是預(yù)報系數(shù).

對比分析上述三個模型，以適用性及簡捷性為判斷條件，本文選擇自回歸模型為預(yù)報模型，該模型原理簡單，參數(shù)模擬不復(fù)雜，預(yù)報簡捷.

3 自回歸模型求解

本文根據(jù)數(shù)據(jù)特點，選用AR（p）模型進(jìn)行預(yù)報.AR（P）模型參數(shù)一個是預(yù)報階數(shù)，二是回歸系數(shù).

3.1 預(yù)報階數(shù)

預(yù)報階數(shù)P采用經(jīng)驗公式法，

且P多取3階或5階.實際預(yù)報中，選取幾種不同階數(shù)分別進(jìn)行預(yù)報并比較預(yù)報值，選擇吻合最好的做為應(yīng)取階數(shù).

3.2 回歸系數(shù)

回歸系數(shù)一般按三步分別計算.

首先計算不同時間間隔的自相關(guān)函數(shù)，不同時間間隔即不同階數(shù).一個單位時間間隔的自相關(guān)函數(shù)計算公式為：τ個單位的自相關(guān)函數(shù)為R（τ），計算公式為：

其次，為比較不同序列，將自相關(guān)函數(shù)標(biāo)準(zhǔn)化，即用自相關(guān)函數(shù)值除以方差：

公式5至公式8中，ΔX是序列的距平；δX2是序列方差；b1,b2,b3是預(yù)報系數(shù)；γ（1）,γ（2）,γ（3）是標(biāo)準(zhǔn)化的自相關(guān)函數(shù).

4 實例應(yīng)用

4.1 實例簡介

本文以丹東地區(qū)某中小河流入干流控制水文站1990-2014年合計25年系列年最高洪水位資料進(jìn)行預(yù)報.

該河流位于丹東暴雨中心區(qū)，年平均降雨量1100mm.流域面積40km2，流域植被保持良好，下墊面條件基本平穩(wěn).1990年以來，流域內(nèi)無大的水利工程建設(shè)，大的取水、排水設(shè)施.最高洪水位主要受天然降水強度、降水量、流域產(chǎn)匯流等水文要素影響.

4.2 平穩(wěn)性檢驗本文對系列資料的平穩(wěn)性檢查，分為定性和定量檢查.定性檢查根據(jù)流域水文特性及最高洪水位物理機理，最高洪水位為隨機變量，不受時間影響.

定量分析主要計算不同時間的均值和方差的變化.分別計算總均值、總方差，1990-2004年合計15年的均值、方差；2000-2014年合計15的均值、方差.均值分別為70.35、70.58、70.57；方差分別為3.19、3.42、3.25.從統(tǒng)計結(jié)果看，25年最高洪水位時間序列為寬平穩(wěn)序列.

4.3 預(yù)報模型

預(yù)報模型建立的依據(jù)是認(rèn)為預(yù)報值與前期的觀測值有關(guān)［4］［5］［6］.因為平穩(wěn)時間序列的自相關(guān)函數(shù)不隨時間變化，最高水位自身前后關(guān)系的線性方程：

式中，Ht+1是未來預(yù)測的最高水位值；ΔHt+1是未來一個時間間隔的預(yù)報距平值；ΔHt,…,ΔHt-p+1為預(yù)報年的上1年、…、p年的最高水位值距平值；b1,b2,…,bp是預(yù)報系數(shù)；P是預(yù)報階數(shù).

4.4 預(yù)報計算及評價

根據(jù)水文站最高洪水位資料,參數(shù)計算見表1.

25年系列資料，預(yù)報階數(shù)為3或5，為了對比預(yù)報精度，本文對兩個階數(shù)分別進(jìn)行計算.標(biāo)準(zhǔn)化的自相關(guān)函數(shù)值計算結(jié)果見表2.

表1 最高水位自回歸模型參數(shù)計算表

序號 年份 水位 △Ht△Ht2△Ht△Ht+1△Ht△Ht+2△Ht△Ht+3△Ht△Ht+4△Ht△Ht+517 2006 68.67 -1.68 2.816 -0.809 2.497 -0.121 2.094 5.534 18 2007 70.83 0.48 0.232 -0.717 0.035 -0.602 -1.590 -0.549 19 2008 68.86 -1.49 2.214 -0.107 1.857 4.907 1.693 -3.068 20 2009 70.42 0.07 0.005 -0.090 -0.237 -0.082 0.148 0.158 21 2010 69.1 -1.25 1.558 4.116 1.420 -2.573 -2.736 22 2011 67.05 -3.30 10.877 3.753 -6.800 -7.229 23 2012 69.21 -1.14 1.295 -2.347 -2.494 24 2013 72.41 2.06 4.252 4.520 25 2014 72.54 2.19 4.805合計 1758.7 79.658 11.120 -1.138 -20.733 -17.478 6.197平均 70.35 3.186 0.445 -0.046 -0.829 -0.699 0.248

表2 標(biāo)準(zhǔn)化的自相關(guān)函數(shù)值

4.4.1 3階回歸預(yù)報

預(yù)報系數(shù)的線性方程組為：

根據(jù)EXCELL中MINVERSE,MMULT函數(shù)，計算得：b1為0.135495，b2為0.002916，b3為-0.258751.把計算數(shù)據(jù)代入公式9和公式10，得：

4.4.2 5階回歸預(yù)報

如上，代入數(shù)據(jù)，計算得到：b1為0.115146，b2為0.034907，b3為-0.237259，b4為-0.174399，b5為133096.把計算數(shù)據(jù)代入公式9和公式10，得：

4.4.3 預(yù)報評價

該站2015年最高洪水位為71.13米，其中3階預(yù)測值偏小，5階預(yù)測值偏大.3階預(yù)測誤差為0.26%，5階預(yù)報誤差為0.31%；預(yù)報結(jié)果說明平穩(wěn)時間序列預(yù)報精度較高,3階預(yù)報就能滿足要求，一般可不進(jìn)行5階預(yù)報，可減少工作量，計算更簡捷.

5 結(jié)論

1.通過對比分析，采用自回歸模型進(jìn)行水位預(yù)報，確定了預(yù)報階數(shù)和回歸系數(shù)的計算方法及步驟.

2.案例預(yù)報結(jié)果表明，自回歸模型適合水位預(yù)報，預(yù)報精度高，3階預(yù)報精度好于5階預(yù)報精度，一般可只采用3階預(yù)報，可減少工作量，計算更簡捷.

3.平穩(wěn)時間序列較簡單，但系列數(shù)據(jù)必須滿足平穩(wěn)性，在推廣應(yīng)用中，需對系列數(shù)據(jù)進(jìn)行定性、定量評價分析.

參考文獻(xiàn)：

〔1〕張美英，何杰.時間序列預(yù)測模型研究綜述［J］.數(shù)學(xué)的實踐與認(rèn)識，2011（18）：189-195.

〔2〕趙春艷.平滑轉(zhuǎn)換自回歸模型的平穩(wěn)性問題研究［J］.數(shù)量經(jīng)濟(jì)技術(shù)經(jīng)濟(jì)研究，2012（01）：152-160.

〔3〕林慧君，徐榮聰.組合ARMA與SVR模型的時間序列預(yù)測［J］.計算機與現(xiàn)代化，2009（08）：19-22.

〔4〕郭春霞，諸云強，孫偉.不同時間尺度、季節(jié)的氣溫數(shù)據(jù)空間平穩(wěn)特征及其對插值結(jié)果的影響 ［J］.地理研究，2015 （09）：1675-1684.

〔5〕崔磊，遲道才，曲霞.基于小波消噪的平穩(wěn)時間序列分析方法在降雨量預(yù)測中的應(yīng)用 ［J］.中國農(nóng)村水利水電，2010 （10）：39-43.

〔6〕吳東杰，王金生，滕彥國.小波分解與變換法預(yù)測地下水位動態(tài)［J］.水利學(xué)報，2004（05）：39-45.

中圖分類號：O211

文獻(xiàn)標(biāo)識碼：A

文章編號：1673-260X（2016）06-0005-03

收稿日期：2016-03-22