王玉振

(河南水利與環(huán)境職業(yè)學(xué)院，河南鄭州450008)

及時獲取大壩變形的演變規(guī)律并作出準(zhǔn)確預(yù)報，對大壩安全具有重要意義[1]。由于大壩變形受溫度、水壓和時間效應(yīng)等諸多因素影響，呈現(xiàn)出高度的非線性和非平穩(wěn)性，難以建立較為準(zhǔn)確的變形量與影響因素之間復(fù)雜的關(guān)系模型，因此，預(yù)測效果不佳[2]。目前，繼BP(Back Propagation，BP)神經(jīng)網(wǎng)絡(luò)之后的小波神經(jīng)網(wǎng)絡(luò)(Wavelet Neural Ntworks，WNN)模型被廣泛用于非線性、非平穩(wěn)性時間序列預(yù)測，其結(jié)合了小波變換和人工神經(jīng)網(wǎng)絡(luò)各自的特點，具備極強(qiáng)的自學(xué)習(xí)能力和小波變換局部化等優(yōu)點[3]，能夠有效解決局部極小值問題。但在實際應(yīng)用中發(fā)現(xiàn)，其存在小波分解無法實現(xiàn)的自適應(yīng)多分辨率分析問題[4-5]。而經(jīng)驗?zāi)B(tài)分解法[6-7](Empirical Mode Decomposition，EMD)在信號處理和分析中具有自適應(yīng)的良好特性，能夠有效彌補(bǔ)小波神經(jīng)網(wǎng)絡(luò)自適應(yīng)多分辨率分析的不足。本文考慮影響因素對大壩變形的影響，并采用組合模型的思想，將經(jīng)驗?zāi)B(tài)分解法和小波神經(jīng)網(wǎng)絡(luò)相結(jié)合，充分發(fā)揮兩者各自的優(yōu)勢，將有望對大壩變形中隱含的特征信息進(jìn)行多層次、多角度挖掘和分析，從而提高預(yù)測精度。

1 基于EMD的小波神經(jīng)網(wǎng)絡(luò)預(yù)測模型

1.1 經(jīng)驗?zāi)B(tài)分解原理

EMD是N.E.Huang等人提出的一種信號分解技術(shù)，無須預(yù)先設(shè)定基函數(shù)，就可以自適應(yīng)地將非平穩(wěn)、非線性時間序列中的的波動和趨勢分離開，形成若干個固有模態(tài)分量IMF和1個余項，其實質(zhì)是平穩(wěn)化數(shù)據(jù)處理。對于任意一變形時間序列S(t)，對其進(jìn)行經(jīng)驗?zāi)B(tài)分解，步驟參照文獻(xiàn)[8]。因此，原始序列可以表示為

(1)

式中，S(t)為原始數(shù)據(jù)監(jiān)測信號；t=1，2，3，…，m為時間采樣序列，其中，m為采樣總次數(shù)；IMFi為分解所得的第i個IMF分量；i=1，2，3，…，n為IMF分量的次序，是按頻率由高到低依次排列；B為單調(diào)殘余函數(shù)(EMD分解后的趨勢分量)。

1.2 游程判定法重構(gòu)原理

1.3 預(yù)測模型及流程

大壩變形的影響因素分為水位、溫度、時效三部分，設(shè)大壩在不同時間段的變形序列為{Y(x)=(tx，Hx，Tx)，x=1，2，…，n}，其中，t為時效量；H為上游水位影響因子；T為溫度影響因子。經(jīng)EMD對大壩變形原始數(shù)據(jù)Y分解得到n個IMF和一個剩余量B，然后通過游程判定法對各分量進(jìn)行重構(gòu)得到高頻項F1、中頻項F2和低頻項F3，設(shè)對應(yīng)的預(yù)測值分別為E1、E2和E3，模型預(yù)測流程如圖1所示，則最終的預(yù)測值

G=E1+E2+E3

(2)

圖1 模型預(yù)測流程

1.4 模型精度評定

為了評價各模型預(yù)測性能的優(yōu)劣，本文采用均方根誤差(Root Mean Square Error，RMSE)和平均絕對誤差(Mean Absolute Error，MAE)作為模型評價指標(biāo)。

平均絕對值誤差(MAE)

(3)

均方根誤差(RMSE)

(4)

2 算例

2.1 實驗數(shù)據(jù)獲取及分析

本文以豐滿大壩30號壩段1985年1月4日～1987年6月22日的壩頂水平位移及同期庫水位、溫度和時效的實測資料為例[10]，共140期觀測數(shù)據(jù)，建立模型進(jìn)行分析，如圖2所示。

圖2 大壩變形水平位移序列

由圖2可見，該大壩變形較劇烈，無明顯規(guī)律可循，若直接建立模型進(jìn)行預(yù)測分析，很難得到可靠的結(jié)果。本文先采用EMD對其進(jìn)行預(yù)處理，分解結(jié)果見圖3。

圖3 EMD分解結(jié)果

由圖3可知，EMD將變形序列自適應(yīng)分解為4個IMF和一個殘余分量B，使各分量變化曲線比原來變形曲線(圖2)更光滑和平穩(wěn)。EMD可以降低其非平穩(wěn)性，將相同和類似信息獨立出來，有利于模型的建立和分析。但經(jīng)EMD分解得到的分量個數(shù)多，會增加建模次數(shù)。因此，有必要采用游程判定法對各分量進(jìn)行波動程度分析。IMF1、IMF2、IMF3、IMF4、B的游程總數(shù)分別為66、15、7、3、2。

表1 各模型預(yù)測結(jié)果與實際值對比

圖4 重構(gòu)結(jié)果

由于游程總數(shù)反映出各分量的波動程度，所以將IMF1作為高頻分量，IMF2和IMF3疊加作為中頻分量，IMF4和余量B作為低頻分量。重構(gòu)結(jié)果如圖4所示。

經(jīng)重構(gòu)后，可以減少建模次數(shù)和簡化數(shù)據(jù)處理過程，提高預(yù)測效率。為驗證本文算法的可行性和有效性，建立4種方案進(jìn)行對比分析：①方案1。多元回歸模型。②方案2。BP神經(jīng)網(wǎng)絡(luò)模型。③方案3。WNN模型。④方案4?；贓MD的小波神經(jīng)網(wǎng)絡(luò)大壩變形預(yù)測模型。

2.2 預(yù)測結(jié)果對比分析

本文分別利用WNN訓(xùn)練和測試經(jīng)EMD分解重構(gòu)后的高、中和低頻分量，并疊加各個預(yù)測值即為最終大壩變形預(yù)測結(jié)果，各模型預(yù)測結(jié)果見表1。

由表1可知，本文模型求出的殘差較小，最大殘差為0.7 mm，最小只有0.04 mm，其余殘差都小于1 mm，而其他模型的殘差均有大于1 mm的情況，不能保證較好的全局預(yù)測精度。因此，本文模型能較好地表達(dá)大壩變形分量與其影響因素之間復(fù)雜的關(guān)系，預(yù)測值與實際值吻合最好，均優(yōu)于其他模型，預(yù)測精度較高。為了更直觀地探討各模型的預(yù)測效果，將各模型的預(yù)測值與大壩變形的真實值進(jìn)行對比分析，見圖5。

圖5 各模型預(yù)測結(jié)果與實際對比

由圖5知，方案1、2和3預(yù)測不穩(wěn)定，部分預(yù)測值與實際值相比偏差較大，尤其在變形幅度較大的第121到127期、129到135期，波動比較大。而方案4的預(yù)測值與實際值吻合最好，表明基于EMD的小波神經(jīng)網(wǎng)絡(luò)預(yù)測模型能較好地反映大壩變形規(guī)律，說明EMD能自適應(yīng)地分析和提取變形時間序列中隱含的特征信息，將相同或類似的信息獨立出來，使得模型的自身性得到充分發(fā)揮，有效地彌補(bǔ)了小波神經(jīng)網(wǎng)絡(luò)無法實現(xiàn)自適應(yīng)多分辨率分析的不足，取得了較好的效果。同時，采用游程判定法將相同或相似波動程度的信息重新組合，減少了建模次數(shù)，提高了預(yù)測效率。為了進(jìn)一步探討各模型的預(yù)測效果，各模型預(yù)測精度見表2。

由表2可看出，多元回歸、BP和WNN模型精度指標(biāo)較差，不能很好地表達(dá)大壩變形與其影響因素之間復(fù)雜的函數(shù)關(guān)系。而本文利用EMD對大壩變形原始數(shù)據(jù)進(jìn)行分解，能有效挖掘其內(nèi)部規(guī)律，建立影響因素與變形分量之間復(fù)雜的關(guān)系模型，預(yù)測精度均優(yōu)于WNN、BP和多元回歸模型，具有較好的泛化能力，可信度較高。這是由于EMD與WNN的結(jié)合，充分發(fā)揮了兩者各自的優(yōu)勢，提高了預(yù)測精度。綜上，說明有必要先對非線性、非平穩(wěn)時間序列進(jìn)行EMD分析，這為大壩變形預(yù)測研究提供了一種新思路。

表2 各模型精度對比

3 結(jié) 語

目前，WNN是時間序列預(yù)測方法中較為理想的一種預(yù)測方法，但是其具有在預(yù)測中無法實現(xiàn)自適應(yīng)多分辨率分析的不足，而利用EMD可以解決此問題。本文提出一種基于EMD的小波神經(jīng)網(wǎng)絡(luò)大壩變形預(yù)測模型，經(jīng)理論和算例分析，并與多元回歸、BP神經(jīng)網(wǎng)絡(luò)和WNN模型對比分析表明，該算法預(yù)測精度較高，可以用于大壩變形預(yù)測。