基于灰色神經(jīng)網(wǎng)絡-加權馬爾可夫鏈的大壩變形監(jiān)控模型及預報研究

何 啟，戴 波(1.河海大學 水利水電學院，南京 210098；2.河海大學水文水資源與水利工程科學國家重點實驗室，南京 210098)

1 研究背景

近些年建成或在建的大壩規(guī)模越來越大，如二灘、向家壩、錦屏一級、溪洛渡等[1]。因此，建立大壩安全自動化監(jiān)測系統(tǒng)，分析變形規(guī)律及時進行預報至關重要。許多學者對此進行了研究：楊杰等[2]利用誤差逆?zhèn)鞑ニ惴ǖ腂P神經(jīng)網(wǎng)絡模型對大壩位移監(jiān)測數(shù)據(jù)進行擬合及預報研究；何鮮峰等[3]對實測數(shù)據(jù)與確定性分量預測結果的殘差序列進行分析，把混沌分量預報模型和分形預測模型疊加構成混合預測模型；閆濱等[4]把單純形法嵌入遺傳算法并與神經(jīng)網(wǎng)絡有機結合，構成遺傳單純形神經(jīng)網(wǎng)絡算法模型；何政翔等[5]把模糊聚類和多元逐步回歸法進行結合對大壩變形資料進行分析并作出預報；此外還有多尺度小波分析的自回歸預測[6]，小波消噪與BP神經(jīng)網(wǎng)絡[7]結合等預測方法。

大壩監(jiān)測系統(tǒng)既含已知信息，又有確定、未確定的信息，因此可以把大壩看成一個灰色系統(tǒng)?；疑碚撚糜陬A測的模型主要是GM(1,1)，它可對時間短、數(shù)量少、波動性不大的數(shù)據(jù)序列建模，預測效果良好。BP神經(jīng)網(wǎng)絡在較大量訓練樣本數(shù)據(jù)的基礎上才可以準確地反映非線性函數(shù)的映射關系，這在大壩蓄水初期是較難實現(xiàn)的。而馬爾可夫鏈預測理論是通過狀態(tài)之間的轉移概率來推測系統(tǒng)未來發(fā)展變化的，適用于數(shù)據(jù)隨機波動性較大的情況[8]。本文結合3者的優(yōu)點，構建灰色神經(jīng)網(wǎng)絡-加權馬爾可夫鏈預測模型，首先利用灰色神經(jīng)網(wǎng)絡對少量樣本數(shù)據(jù)進行非線性擬合，揭示大壩水平位移的整體變化趨勢，然后將其相對誤差作為隨機波動數(shù)據(jù)序列，利用加權馬爾可夫鏈確定相對誤差的轉移概率來修正模型預測值，提高預測精度[9]。

2 灰色神經(jīng)網(wǎng)絡-加權馬爾可夫鏈預測模型的構建

由于原始數(shù)據(jù)數(shù)量的有限性，數(shù)據(jù)本身的波動性和無序性，欲將灰色模型預測精度控制在一定范圍內比較困難。神經(jīng)網(wǎng)絡需要大量數(shù)據(jù)才能取得良好的預測效果，故引入灰色理論構建灰色神經(jīng)網(wǎng)絡模型可減少對樣本數(shù)量的依賴。雖然灰色神經(jīng)網(wǎng)絡模型在小樣本數(shù)據(jù)預測中優(yōu)于傳統(tǒng)GM模型，但一些預測值的相對誤差較大，難以滿足預測要求[10]，而馬氏鏈適用于數(shù)據(jù)波動性較大的預測問題。

2.1 GM(1，1)模型

記原始序列X(0)={x(0)(k1),x(0)(k2),…,x(0)(kn)}，x(0)≥0,k=1,2,3,…,n，X(1)為X(0)對原始數(shù)據(jù)作一次累加：

(1)

GM(1，1)模型為：

(3)

k=1,2,…,n

(4)

GM(1，1)模型預測值與實際值的相對誤差為：

(5)

式(5)反應的是擬合值與原始值的偏離程度。

2.2 灰色神經(jīng)網(wǎng)絡模型

n個參數(shù)的灰色神經(jīng)網(wǎng)絡模型的微分方程表達式為：

(6)

式中：y2,…,yn為系統(tǒng)輸入?yún)?shù)；y1為系統(tǒng)輸出參數(shù)；a,b1,b2,…,bn-1為微分方程系數(shù)。

式(9)的時間響應式為：

(7)

令：

式(10)可以作如下轉化：

(8)

將轉化后的公式用擴展的BP神經(jīng)網(wǎng)絡進行映射得到n個輸入值、1個輸出值。網(wǎng)絡拓撲結構如下圖1所示。

圖1 灰色神經(jīng)網(wǎng)絡拓撲結構Fig.1 Grey neural network topology

圖1中，t為輸入?yún)?shù)序號；y2(t),…,yn(t)為網(wǎng)絡輸入?yún)?shù)；ω21,ω22,…,ω2n，ω31,ω32,…,ω3n為網(wǎng)絡權值；y1為網(wǎng)絡預測值。

2.3 加權馬爾可夫鏈預測的基本步驟

根據(jù)實測資料建立訓練樣本，利用灰色神經(jīng)網(wǎng)絡可擬合出一條非線性函數(shù)曲線，確定大壩變形的趨勢，然后利用馬爾可夫鏈的“無后效”性分析相對誤差的波動規(guī)律，修正預測值。對幾種傳統(tǒng)馬爾可夫鏈預測方法的對比研究表明:加權馬爾可夫鏈預測理論在數(shù)據(jù)序列隨機波動性較大的預測中精度最高[11]。

根據(jù)馬爾可夫鏈理論，首先計算灰色神經(jīng)網(wǎng)絡模型ε(k)的轉移概率矩陣，然后對同一狀態(tài)下的概率進行加權求和，其最大值對應的狀態(tài)即為加權馬爾可夫鏈的預測狀態(tài)，并依此對灰色神經(jīng)網(wǎng)絡模型結果進行修正[12]。具體步驟如下。

(2)按(1)所建立的分級標準，確定各擬合值相對誤差所對應的狀態(tài)。

(3)“馬氏性”檢驗。通常用χ2統(tǒng)計量來檢驗隨機離散序列的馬氏鏈。設所討論的序列包含m個可能的狀態(tài)，(fij)ij∈E代表狀態(tài)轉移頻數(shù)矩陣，將(fij)ij∈E的各列之和除以各行各列的總和得到“邊際概率”，記為：

(9)

(4)計算各階自相關系數(shù)。即:

(10)

式中：rk為第k階自相關系數(shù)；xl為第l時段的相對誤差；n表示相對誤差序列長度。

對各階相關系數(shù)規(guī)范化：

(11)

作為各種步長的馬爾可夫鏈的權重。

p(K)=(p(K)ij)ij∈E

(12)

(6)以離預測值最近的K個相對誤差值為初始狀態(tài)，根據(jù)其相應的轉移概率矩陣即可預測該時段相對誤差的狀態(tài)概率p(K)i,i∈E,K為步長，K=1,2,…,m。

(7)確定預測狀態(tài)。對同一狀態(tài)的各預測概率進行加權求和，得到相對誤差的轉移概率為：

(13)

max{pi,i∈E}所對應的狀態(tài)為該時段相對誤差的預測狀態(tài)。

(8)利用狀態(tài)特征值結合線性插值的方法對預測值進行推斷[13]， 先計算狀態(tài)特征值：

(14)

(15)

也可用上限法：

x)

(17)

(9)對后續(xù)時間預測，將前一時間預測值加入原序列中重新構置馬爾可夫鏈，重復步驟(1)～(8)可進行預測。

3 實例分析

某水電站大壩為碾壓混凝土重力壩，位于福建省境內，該工程除發(fā)電外，在防洪、航運、水產(chǎn)養(yǎng)殖等方面也發(fā)揮著效益。水電站建成后，是閩西南最大水電站，也是省網(wǎng)主要電源之一，可以擔任省網(wǎng)調峰、調頻任務，同時對增加網(wǎng)內有調節(jié)能力的水電比重發(fā)揮作用，改善了省網(wǎng)電源北重南輕布局。

3.1 預測因子的選擇與處理

為加快學習的速度，可對訓練樣本數(shù)據(jù)按如下公式進行歸一化處理：

X′=0.1+0.8 (X-Xmin)/(Xmax-Xmin)

(18)

式中：Xmax、Xmin分別表示每組樣本數(shù)據(jù)的最大、最小值。

3.2 灰色神經(jīng)網(wǎng)絡的大壩變形預測結果

根據(jù)4號壩段PL7正垂線 方向(上下游方向為正)位移監(jiān)測資料建立灰色神經(jīng)網(wǎng)絡模型。選取蓄水初期2003年11月14日至2003年12月20日的實測數(shù)據(jù)作為訓練樣本用于擬合計算，共37組，取12月20-26日6組數(shù)據(jù)作為預測樣本，評價模型的預測性能。

根據(jù)原數(shù)據(jù)序列建立灰色神經(jīng)網(wǎng)絡模型并運用Matlab程序訓練、預測，預測結果如表1所示。

表1 水平位移預測結果與實測值的相對誤差Tab.1 Prediction of horizontal displacement and the relative error of the measured value

從表1中可以看出，灰色神經(jīng)網(wǎng)絡適應于小樣本預測，但某些相對誤差相對較大，因此可用加權馬爾科夫鏈理論對結果進行進一步修正。

3.3 加權馬爾可夫鏈擬合及預測結果

(2)根據(jù)狀態(tài)區(qū)間確定各擬合值相對誤差對應的狀態(tài)如表2中“5”縱欄所示。

表2 灰色神經(jīng)網(wǎng)絡預測結果、相對誤差及狀態(tài)Tab.2 Grey neural network predicting, relative residuals and states

(3)不妨直接對37組相對誤差作馬氏性檢驗，根據(jù)步驟(5)和表2的資料可計算進一步狀態(tài)轉移頻數(shù)矩陣和一步轉移概率矩陣：

(4)根據(jù)步驟(4)計算各階自相關系數(shù)，如表3所示。

(5)根據(jù)步驟(5)計算相對誤差的K步狀態(tài)轉移概率矩陣p(K)：

表3 各階自相關系數(shù)和各種步長的馬爾可夫鏈權重Tab.3 The weights of Markov chain of various steps and various auto-relative coefficients

同理，將計算出的預測值加入原序列中重新構置馬爾可夫鏈，重復過程 (1)-(7)計算得2003年12月22日-2003年12月26日的預測值，相對誤差范圍為0.055%～1.625%，精度較高。本文又采用GM(1,1)模型進行對比預測，結果見表5和圖2。

表4 2003-12-21水平位移狀態(tài)預測Tab.4 Prediction of horizontal displacement state on 2003-12-21

表5 不同模型預測結果比較Tab.5 Comparison of the results of different prediction model

圖2 不同模型預測結果比較Fig.2 Comparison of the results of different prediction model

通過對比3種方法的預測結果，可知灰色神經(jīng)網(wǎng)絡-加權馬爾可夫鏈模型的預測誤差顯著降低，實測值與預報值更為接近，預測精度優(yōu)于其他兩種模型，為小樣本、波動性大的數(shù)據(jù)預測提供了新的方法參考。

4 結 語

本文針對蓄水初期階段大壩監(jiān)測數(shù)據(jù)有限且波動性大的特點，結合灰色理論適應于小樣本、BP神經(jīng)網(wǎng)絡高效準確擬合非線性函數(shù)、馬爾可夫鏈適應于波動性比較大數(shù)據(jù)的優(yōu)點，建立了灰色神經(jīng)網(wǎng)絡-加權馬爾可夫鏈模型。通過某碾壓混凝土壩4號壩段蓄水初期水平位移監(jiān)測資料進行預測，結果表明該模型對小樣本、數(shù)據(jù)序列波動大進行預測精度較高。由本文還可得到如下結論。

(1)灰色神經(jīng)網(wǎng)絡-加權馬爾可夫鏈模型不僅可對大壩變形進行短期高精度預報，隨著監(jiān)測數(shù)據(jù)的增加也可為中長期預報提供理論方法。

(2)在完成一步預測后，將預測信息重新加入原模型進行下一步預測，可對資料數(shù)據(jù)進行動態(tài)更新，跟蹤大壩變形的發(fā)展變化趨勢，也提高了預測精度。

□

[1] 吳中如.水工建筑物安全監(jiān)控理論及其應用[M].南京:河海大學出版社,1990.

[2] 楊 杰,吳中如,顧沖時，等.大壩變形監(jiān)測的BP網(wǎng)絡模型與預報研究[J].西安理工大學學報,2001,17(1):25-29.

[3] 何鮮峰,顧沖時,谷艷昌，等.分形-混沌混合預測模型在大壩安全監(jiān)測中的應用[J].武漢大學學報(工學版),2008,41(1):45-49.

[4] 閏 濱,周 晶.基于遺傳單純形神經(jīng)網(wǎng)絡的大壩變形監(jiān)控模型[J].水力發(fā)電學報,2007,26(4):110-114.

[5] 何政翔,張慧莉.模糊聚類方法應用于大壩變形監(jiān)測資料分析[J].水力發(fā)電,2013,(11):59-61.

[6] 徐 偉,何金平.基于多尺度小波分析的大壩變形自回歸預測模型[J].武漢大學學報(工學版),2012,45(3):286-289.

[7] 黃世秀,洪天求,高 飛，等.基于小波消噪及BP神經(jīng)網(wǎng)絡的大壩變形分析[J].人民長江,2011,42(9):90-93.

[8] 杜 川,梁秀娟,王中凱，等.改進灰色-馬爾科夫模型在年降水量預測中的應用研究[J].節(jié)水灌溉,2014,(6):32-36.

[9] 楊皓翔,梁 川.基于加權灰色一馬爾可夫鏈模型的城市需水預測[J].長江科學院院報,2015,32(7):15-21.

[10] 薛鵬松,馮民權,邢肖鵬.基于馬爾科夫鏈改進灰色神經(jīng)網(wǎng)絡的水質預測模型[J].武漢大學學報(工學版),2012,45(3) 319-324.

[11] 夏樂天,朱元牲.馬爾可夫鏈預測方法的統(tǒng)計試驗研究[J].水利學報,2007,10(增):372-378.

[12] 趙 玲,許宏科.基于灰色加權馬爾可夫SCGM(1,1)的交通事故預測[J].計算機工程與應用,2012,48(31):11-15.

[13] 夏樂天.馬爾可夫鏈預測方法及其在水文序列中的應用研究[D].南京：河海大學,2005.

[14] 顧沖時,吳中如.碾壓混凝土壩安全診斷與預警的理論和方法[M].南京:河海大學出版社，2012.