孔慶梅,喬海山,費(fèi)新鋒,武志剛,葛佳豪

(1.國(guó)家電投集團(tuán)青海黃河電力技術(shù)有限責(zé)任公司,西寧 810008；2.西安理工大學(xué)水利水電學(xué)院,西安 710048)

0 前 言

水庫(kù)大壩在運(yùn)行過(guò)程中不僅受各種動(dòng)、靜循環(huán)荷載和地震、滑坡等災(zāi)害的作用，還要受混凝土自身材料性質(zhì)改變和物理、化學(xué)作用導(dǎo)致的侵蝕與腐蝕，由此導(dǎo)致壩體在運(yùn)行期不同階段安全性態(tài)發(fā)生改變[1]。大壩失事后果嚴(yán)重，及時(shí)準(zhǔn)確了解大壩安全性態(tài)，保證其穩(wěn)定運(yùn)行十分重要[2]。大壩變形能夠客觀地判斷出工程是否處于安全性態(tài)，對(duì)大壩變形過(guò)程中的實(shí)測(cè)資料進(jìn)行分析，以此來(lái)建立誤差較小的安全監(jiān)控模型是檢驗(yàn)工程是否安全的主要手段之一[3-5]。工程應(yīng)用中針對(duì)具有復(fù)雜非線性特點(diǎn)的大壩變形監(jiān)測(cè)數(shù)據(jù)，已經(jīng)發(fā)展出了傳統(tǒng)數(shù)學(xué)統(tǒng)計(jì)模型、確定性模型、組合模型及混合模型用于大壩的安全監(jiān)控[6]。其中確定性模型和統(tǒng)計(jì)模型應(yīng)用較廣，但是由于壩體自身材料特性和運(yùn)行環(huán)境均具有一定的不確定性，導(dǎo)致模型比較復(fù)雜；傳統(tǒng)統(tǒng)計(jì)模型原理清楚且構(gòu)造簡(jiǎn)單，但預(yù)測(cè)效果一般[6]。

近年來(lái)，相關(guān)學(xué)者在大壩變形分析與變形預(yù)測(cè)方面作了開(kāi)展了深入的研究。胡德秀[7]提出了基于穩(wěn)健估計(jì)極限學(xué)習(xí)機(jī)ELM的大壩變形安全監(jiān)控模型，并通過(guò)實(shí)例驗(yàn)證了模型的抗粗差性和對(duì)大壩安全監(jiān)測(cè)非線性的適用性。對(duì)于環(huán)境影響因子之間的相關(guān)性問(wèn)題，陳詩(shī)怡[6]采用Copula函數(shù)非線性相關(guān)檢驗(yàn)確定最優(yōu)因子集，結(jié)合隨機(jī)森林算法建立大壩安全監(jiān)控模型。宋洋[8]利用自適應(yīng)噪聲完全集合經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)將原始變形數(shù)據(jù)分解為頻率不同、復(fù)雜度差異明顯的一組固有模態(tài)函數(shù)，并在此基礎(chǔ)上使用排列熵方法進(jìn)行合并重組，對(duì)若干重組序列分別構(gòu)建LSTM模型進(jìn)行預(yù)測(cè)，將預(yù)測(cè)結(jié)果相加得到最終變形預(yù)測(cè)值。馬廣臣[9]利用動(dòng)態(tài)權(quán)重粒子群算法對(duì)自適應(yīng)模糊神經(jīng)網(wǎng)絡(luò)中模糊層的適應(yīng)度值進(jìn)行參數(shù)尋優(yōu)，形成可以尋找最優(yōu)適應(yīng)度值的自適應(yīng)模糊神經(jīng)網(wǎng)絡(luò)。李其峰[5]結(jié)合貝葉斯優(yōu)化算法優(yōu)化門(mén)控制循環(huán)單元(GRU)超參數(shù)，建立大壩安全監(jiān)控模型。王麗蓉[10]等融合動(dòng)態(tài)時(shí)間規(guī)整算法和卷積長(zhǎng)短時(shí)記憶網(wǎng)絡(luò)，構(gòu)建了大壩安全監(jiān)測(cè)整體信息聯(lián)合監(jiān)控模型，該模型具有提取測(cè)點(diǎn)最新相關(guān)性信息和處理時(shí)序數(shù)據(jù)的能力。

目前，基于智能算法的大壩安全監(jiān)控模型已經(jīng)在大壩安全性態(tài)分析中發(fā)揮了重要作用，但很多方法還是存在收斂速度慢和穩(wěn)定性差等問(wèn)題。

隨著計(jì)算科學(xué)技術(shù)的不斷發(fā)展和工程應(yīng)用，很多工程問(wèn)題都可以使用群智能優(yōu)化方法解決[11]。根據(jù)具體自然現(xiàn)象衍生出了如粒子群優(yōu)化算法[12]、人工蜂群算法[13]、細(xì)菌覓食算法[14]和螢火蟲(chóng)算法[15]等群智能優(yōu)化算法，并得到了較廣泛的工程應(yīng)用。但隨著科學(xué)問(wèn)題、工程問(wèn)題越來(lái)越復(fù)雜，很多群智能優(yōu)化算法還是存在計(jì)算時(shí)間過(guò)長(zhǎng)及易陷入局部最優(yōu)等的不足。在簡(jiǎn)化了麻雀?jìng)€(gè)體的覓食行為、反捕食行為的基礎(chǔ)上，薛建凱[16]等于2020年提出了麻雀搜索算法(Sparrow Search Algorithm, SSA)，其具有求解精度高、收斂快、穩(wěn)定性好等優(yōu)點(diǎn)[17]。

本文將具有全局優(yōu)化能力的麻雀優(yōu)化算法SSA引入深度學(xué)習(xí)中的長(zhǎng)短期記憶神經(jīng)網(wǎng)絡(luò)LSTM中，智能優(yōu)化神經(jīng)網(wǎng)絡(luò)超參數(shù)，以減少由于人工設(shè)定參數(shù)而造成過(guò)擬合的問(wèn)題，提高構(gòu)建模型的準(zhǔn)確度和收斂速度。在對(duì)比大壩安全監(jiān)控模型的同時(shí)，對(duì)受大壩異常情況影響的不同測(cè)點(diǎn)進(jìn)行對(duì)比分析。

1 麻雀優(yōu)化算法

SSA算法[16]預(yù)設(shè)發(fā)現(xiàn)者與加入者這兩種麻雀。算法運(yùn)行時(shí)，當(dāng)發(fā)現(xiàn)者發(fā)現(xiàn)食物充足的區(qū)域，加入者就會(huì)受到召喚并向該區(qū)域移動(dòng)，因此麻雀所在區(qū)域的食物是否充足決定了麻雀?jìng)€(gè)體能量的高低。另外，為了避免陷入局部最優(yōu)，在種群中按設(shè)定比例隨機(jī)選取部分麻雀作為預(yù)警者。當(dāng)?shù)?jì)算滿足一定條件時(shí)，預(yù)警者將會(huì)選擇相應(yīng)的更新方式進(jìn)行額外移動(dòng)。

對(duì)于發(fā)現(xiàn)者：

(1)

對(duì)于加入者：

(2)

對(duì)于預(yù)警者：

(3)

2 長(zhǎng)短期記憶神經(jīng)網(wǎng)絡(luò)LSTM

1997年，由Hochreiter[18]等改進(jìn)了循環(huán)神經(jīng)網(wǎng)絡(luò)(Recurrent Neural Network，RNN)中的隱藏層，提出了長(zhǎng)短期記憶神經(jīng)網(wǎng)絡(luò)(Long Short-Term Memory，LSTM)，且在RNN每個(gè)記憶神經(jīng)元上增加了一個(gè)記憶單元。同時(shí)其3個(gè)門(mén)結(jié)構(gòu)(即輸入門(mén)、遺忘門(mén)和輸出門(mén))選擇性記憶梯度下降時(shí)誤差函數(shù)反饋的修正參數(shù)，以保留有價(jià)值信息或丟棄過(guò)時(shí)信息。LSTM通過(guò)引入可控自循環(huán)有效避免了長(zhǎng)時(shí)間序列學(xué)習(xí)過(guò)程中經(jīng)常遇到的梯度消失或梯度爆炸問(wèn)題，適用于處理時(shí)序延遲和間隔冗長(zhǎng)的任務(wù)[19]。其結(jié)構(gòu)如圖1所示[5]。

圖1 LSTM原理及結(jié)構(gòu)

LSTM神經(jīng)網(wǎng)絡(luò)的計(jì)算公式如下：

(4)

公式(4)中:it、ft、ct、ot分別表示輸入門(mén)、遺忘門(mén)、更新后的細(xì)胞狀態(tài)及輸出門(mén)；xt為輸入信息；ht為獲取的輸出信息；W為權(quán)重系數(shù)；bi、bf、bc、bo為偏置量;*表示哈達(dá)瑪積符號(hào)；下標(biāo)t表示累計(jì)時(shí)長(zhǎng)，d。

輸入門(mén)用來(lái)更新信息；遺忘門(mén)決定從上一個(gè)狀態(tài)中丟棄信息或保留信息；輸出門(mén)控制計(jì)算后的狀態(tài)并輸出。LSTM通過(guò)門(mén)結(jié)構(gòu)使之前的輸入信息能夠繼續(xù)向后傳播而不消失，從而具有較好的長(zhǎng)、短時(shí)記憶能力，在預(yù)測(cè)中又有較好的擬合效果，能很好地進(jìn)行變形時(shí)間序列的預(yù)測(cè)分析[8]。

3 混凝土拱壩變形預(yù)測(cè)模型構(gòu)建

建立基于麻雀搜索算法(SSA)的長(zhǎng)短期記憶神經(jīng)網(wǎng)絡(luò)(LSTM)混凝土拱壩變形預(yù)測(cè)模型，其主要有4個(gè)步驟:① 數(shù)據(jù)預(yù)處理；② 模型超參數(shù)進(jìn)行優(yōu)化；③ 模型訓(xùn)練與預(yù)測(cè)；④ 模型性能評(píng)價(jià)。

因?yàn)椴煌沫h(huán)境量如水位、溫度、時(shí)效和效應(yīng)量如位移、應(yīng)力等監(jiān)測(cè)項(xiàng)目具有不同的量綱，也會(huì)有不同的數(shù)量級(jí)[5]。為便于分析，淡化量綱和數(shù)量級(jí)對(duì)預(yù)測(cè)值的影響，對(duì)原始時(shí)間序列進(jìn)行處理。本文通過(guò)公式(5)所示的方法對(duì)數(shù)據(jù)進(jìn)行歸一化處理，將原始時(shí)間序列縮放到[0，1]區(qū)間。

(5)

公式(5)中：Xstd為歸一化后的時(shí)間序列；Xmin為原始時(shí)間序列的最小值；Xmax為原始時(shí)間序列的最大值。

然后選取具有代表性的壩頂測(cè)點(diǎn)徑向位移數(shù)據(jù)進(jìn)行訓(xùn)練。為提高模型的精度，將LSTM嵌入到麻雀優(yōu)化算法中，對(duì)其超參數(shù)進(jìn)行優(yōu)化。搜索模型中最大迭代次數(shù)、初始學(xué)習(xí)率、隱藏單元個(gè)數(shù)共3個(gè)參數(shù)的最優(yōu)化值，本文設(shè)定的搜索范圍分別為[100,500]、[0.0001,0.001]和[10,400]，miniBatchSize設(shè)置為128。將尋找到的最優(yōu)超參數(shù)組合設(shè)定給LSTM，向其輸入訓(xùn)練集的實(shí)測(cè)數(shù)據(jù)得到預(yù)測(cè)數(shù)據(jù)，對(duì)預(yù)測(cè)結(jié)果進(jìn)行反歸一化處理，得到預(yù)測(cè)結(jié)果。

為驗(yàn)證本文所提出模型的準(zhǔn)確性和適用性，將SSA-LSTM模型的預(yù)測(cè)結(jié)果與BP神經(jīng)網(wǎng)絡(luò)、極限學(xué)習(xí)機(jī)(ELM)和支持向量機(jī)(SVM)等傳統(tǒng)監(jiān)控模型的預(yù)測(cè)結(jié)果進(jìn)行對(duì)比分析。并通過(guò)其他測(cè)點(diǎn)的實(shí)測(cè)數(shù)據(jù)驗(yàn)證該方法的適用性。綜上，基于SSA-LSTM的混凝土拱壩變形預(yù)測(cè)流程見(jiàn)圖2。

圖2 SSA-LSTM的混凝土拱壩變形預(yù)測(cè)流程

4 工程實(shí)例

將建立的LSTM模型應(yīng)用在拉西瓦拱壩垂線監(jiān)測(cè)的徑向位移成果中，以檢驗(yàn)?zāi)Ｐ偷挠行院汪敯粜浴?/p>

4.1 工程概況

拉西瓦水電站擋水建筑物為雙曲薄拱壩，最大壩高250 m，建基面高程2 210.00 m，壩頂高程2 460.00 m，壩頂中心線弧長(zhǎng)466.63 m；兩岸拱座采用半徑向布置,如圖3。該拱壩的水位在2011年已經(jīng)保持較高水位運(yùn)行，2011—2013年間又逐漸上升，但上升幅度較小，2013年后則基本穩(wěn)定。因此作為一個(gè)空間殼體結(jié)構(gòu)，拉西瓦拱壩壩頂?shù)膹较蛭灰票憩F(xiàn)出了明顯的季節(jié)性波動(dòng)。

圖3 拉西瓦拱壩

拉西瓦拱壩垂線監(jiān)測(cè)布置見(jiàn)圖4。本文以壩頂PL2-1、PL4-1和PL5-1共3個(gè)測(cè)點(diǎn)的徑向位移實(shí)測(cè)數(shù)據(jù)為例，檢驗(yàn)SSA-LSTM混凝土拱壩變形預(yù)測(cè)方法的精度和適用性，并與傳統(tǒng)監(jiān)控模型進(jìn)行性能對(duì)比。PL2-1、PL4-1和PL5-1測(cè)點(diǎn)的徑向位移實(shí)測(cè)數(shù)據(jù)見(jiàn)圖5。

圖4 拉西瓦拱壩垂線監(jiān)測(cè)布置 單位：m

圖5 壩頂垂線PL2-1、PL4-1和PL5-1測(cè)點(diǎn)實(shí)測(cè)徑向絕對(duì)位移過(guò)程線

如圖5，由于2013年壩體下游面左岸壩段的保溫板失效，導(dǎo)致左岸壩段PL5-1和拱冠梁PL4-1測(cè)點(diǎn)在2013年底冬季、2014年夏季徑向位移變幅明顯增大，與失效前和修復(fù)后的徑向位移周期性規(guī)律明顯不同，而右岸壩段的PL2-1測(cè)點(diǎn)受影響較小，因此將這3個(gè)典型測(cè)點(diǎn)進(jìn)行對(duì)比，以分析異常位移對(duì)監(jiān)控模型的影響。

4.2 模型結(jié)構(gòu)和參數(shù)

混凝土拱壩在服役期內(nèi)主要受水壓力、溫度、時(shí)效等因素作用發(fā)生變形[20]。因此大壩的徑向位移δ為：

δ=δH+δT+δθ

(6)

公式(6)中：δH、δT、δθ分別為水壓、溫度和時(shí)效因子，具體公式如下：

(1) 水壓因子

(7)

(2) 溫度因子

(8)

(3) 時(shí)效因子

δθ=c1θ+c2ln(1+θ)

(9)

公式(9)中：時(shí)效因子主要包括線性函數(shù)θ和對(duì)數(shù)函數(shù)ln(1+θ)，以模擬拱壩運(yùn)行初期的趨勢(shì)性變形。

以上述8項(xiàng)水壓因子、4項(xiàng)溫度因子和2項(xiàng)時(shí)效因子作為模型輸入，壩頂PL2-1、PL4-1和PL5-1測(cè)點(diǎn)徑向位移的實(shí)測(cè)數(shù)據(jù)作為輸出，分別應(yīng)用SSA-LSTM模型和傳統(tǒng)的BP神經(jīng)網(wǎng)絡(luò)、極限學(xué)習(xí)機(jī)(Extreme Learning Machine，ELM)、支持向量機(jī)(Support Vector Machine，SVM)進(jìn)行訓(xùn)練和預(yù)測(cè)，以進(jìn)行對(duì)比研究。采用GA算法對(duì)SVM模型優(yōu)化，得到最優(yōu)懲戒因子系數(shù)為128，核函數(shù)參數(shù)為0.2；ELM采用默認(rèn)傳遞函數(shù)Sigmoid，隱藏神經(jīng)元數(shù)為15。數(shù)據(jù)樣本為2011—2021年間的共417個(gè)數(shù)據(jù)(一周一測(cè))，各模型均指定85%數(shù)據(jù)用于訓(xùn)練，15%的數(shù)據(jù)用于預(yù)測(cè)和驗(yàn)證。

4.3 模型應(yīng)用及評(píng)價(jià)

4.3.1評(píng)價(jià)指標(biāo)

為了檢驗(yàn)預(yù)測(cè)模型的性能，本文選取了4個(gè)參數(shù)作為評(píng)價(jià)指標(biāo)來(lái)分析模型的預(yù)測(cè)精度，分別為平均絕對(duì)誤差(MAE)、絕對(duì)誤差百分比(MAPE)、均方根誤差(RMSE)和決定系數(shù)(R2)。其計(jì)算公式如下：

(10)

(11)

(12)

(13)

4.3.2模型訓(xùn)練及性能評(píng)價(jià)

使用麻雀搜索算法對(duì)LSTM模型參數(shù)的尋優(yōu)結(jié)果為：最大迭代次數(shù)Max_Epochs=500、初始學(xué)習(xí)率Initial_Learn_Rate=0.005、隱藏單元個(gè)數(shù)num_hidden_units=256。SVM核為徑向基函數(shù)；ELM采用默認(rèn)傳遞函數(shù)SIGMOID，隱藏神經(jīng)元數(shù)為15。不同模型在PL4-1測(cè)點(diǎn)的徑向位移的預(yù)測(cè)結(jié)果如圖6所示，3個(gè)測(cè)點(diǎn)采用不同模型的預(yù)測(cè)性能評(píng)價(jià)指標(biāo)見(jiàn)表1。

圖6 多個(gè)模型在PL4-1測(cè)點(diǎn)的預(yù)測(cè)結(jié)果與實(shí)測(cè)值

表1 PL2-1、PL4-1和PL5-1實(shí)測(cè)徑向位移多個(gè)模型預(yù)測(cè)性能指標(biāo)

由圖6可知，PL2-1、PL4-1和PL5-1測(cè)點(diǎn)的徑向位移均具有明顯的周期性變化，同時(shí)還有較明顯的趨勢(shì)性變化。其中PL5-1受2013年保溫板失效的影響，其徑向位移在2014年變幅明顯變大，PL4-1和PL2-1受影響較小但仍有一定變化。

由圖6可知4種模型均能較好預(yù)測(cè)PL4-1測(cè)點(diǎn)的徑向位移規(guī)律，預(yù)測(cè)結(jié)果基本與實(shí)測(cè)值呈現(xiàn)相同的周期性規(guī)律。其中，BP神經(jīng)網(wǎng)絡(luò)的預(yù)測(cè)結(jié)果較不穩(wěn)定，波動(dòng)較大；SVM模型預(yù)測(cè)結(jié)果的周期性變幅變明顯小于實(shí)測(cè)值；ELM模型預(yù)測(cè)效果較好，但周期性峰值、谷值對(duì)應(yīng)的時(shí)間點(diǎn)明顯晚于實(shí)測(cè)值。SSA-LSTM模型的預(yù)測(cè)效果明顯好于其他模型，較好地體現(xiàn)了該測(cè)點(diǎn)的周期性變化，變幅預(yù)測(cè)也更為準(zhǔn)確，在幾個(gè)對(duì)比方法中表現(xiàn)最優(yōu)。

從預(yù)測(cè)值的評(píng)價(jià)指標(biāo)來(lái)看，使用SSA-LSTM模型的預(yù)測(cè)評(píng)價(jià)指標(biāo)均比其他模型要好而且更加突出：SSA-LSTM模型的平均絕對(duì)誤差、絕對(duì)誤差百分比、均方根誤差和復(fù)相關(guān)系數(shù)分別為0.56 mm、8.59%、0.78 mm和0.98，與其他對(duì)比模型相比明顯更優(yōu)。從預(yù)測(cè)值評(píng)價(jià)指標(biāo)復(fù)相關(guān)系數(shù)R2來(lái)看，相較于BP神經(jīng)網(wǎng)絡(luò)、SVM和ELM方法的預(yù)測(cè)值，SSA-LSTM模型預(yù)測(cè)值的復(fù)相關(guān)系數(shù)分別提高了0.22、0.06和0.03，這表明SSA-LSTM模型所估的偏差范圍較小，相關(guān)性更強(qiáng)，擬合的精準(zhǔn)度也較高；從平均絕對(duì)誤差、絕對(duì)誤差百分比和均方根誤差來(lái)看，SSA-LSTM預(yù)測(cè)的位移誤差相對(duì)其他方法均更低，說(shuō)明該方法的預(yù)測(cè)值與實(shí)測(cè)值更加接近，預(yù)測(cè)效果更好。

對(duì)比不同測(cè)點(diǎn)預(yù)測(cè)結(jié)果的評(píng)價(jià)指標(biāo)可以發(fā)現(xiàn)：采用不同方法對(duì)PL2-1和PL5-1兩個(gè)測(cè)點(diǎn)的預(yù)測(cè)效果整體上基本與PL4-1一致，不同方法的平均絕對(duì)誤差、絕對(duì)誤差百分比、均方根誤差和復(fù)相關(guān)系數(shù)4個(gè)評(píng)價(jià)指標(biāo)均表現(xiàn)出了同樣的趨勢(shì)，且SSA-LSTM模型均明顯優(yōu)于其他方法。但受大壩下游面保溫板失效導(dǎo)致的異常變形的影響，PL5-1測(cè)點(diǎn)的異常變幅在3個(gè)測(cè)點(diǎn)中最大，因此其預(yù)測(cè)效果相比PL4-1和PL2-1兩個(gè)測(cè)點(diǎn)較差，這也符合異常位移對(duì)預(yù)測(cè)模型的影響規(guī)律。

綜上所述，使用SSA-LSTM模型可以很好地預(yù)測(cè)混凝土拱壩非線性變形序列的優(yōu)勢(shì)。相比于其他模型，基于SSA-LSTM的混凝土變形模型可以將深度學(xué)習(xí)在非線性時(shí)間序列預(yù)測(cè)方面的優(yōu)點(diǎn)發(fā)揮出來(lái)，從而使預(yù)測(cè)結(jié)果更加準(zhǔn)確。

5 結(jié) 論

本文將麻雀搜索算法和深度學(xué)習(xí)中的長(zhǎng)短期記憶神經(jīng)網(wǎng)絡(luò)結(jié)合，用于混凝土拱壩徑向位移變形監(jiān)測(cè)分析中，構(gòu)建了基于SSA-LSTM模型的混凝土拱壩變形預(yù)測(cè)方法。通過(guò)實(shí)例驗(yàn)證了本模型是比較可靠且實(shí)用的，主要結(jié)論如下：

(1) 通過(guò)麻雀搜索算法優(yōu)化LSTM模型中的最大迭代次數(shù)、初始學(xué)習(xí)率、隱藏單元個(gè)數(shù)3個(gè)參數(shù)，避免了人工設(shè)置參數(shù)帶來(lái)的不確定性，避免模型陷入局部最優(yōu)，提高了模型的精度。

(2) 通過(guò)與BP神經(jīng)網(wǎng)絡(luò)、SVM和ELM模型進(jìn)行對(duì)比，充分證明了SSA-LSTM模型具有較強(qiáng)的預(yù)測(cè)能力，可用于大壩安全預(yù)警預(yù)報(bào)。

(3) 由于突發(fā)事件導(dǎo)致的大壩某測(cè)點(diǎn)發(fā)生異常位移，會(huì)導(dǎo)致大壩預(yù)測(cè)模型的精度出現(xiàn)一定程度的降低。

本文僅針對(duì)單個(gè)測(cè)點(diǎn)逐個(gè)進(jìn)行回歸預(yù)測(cè)分析，但大壩是一個(gè)高度復(fù)雜的系統(tǒng)，具有較強(qiáng)的不確定性。因此，在開(kāi)展多測(cè)點(diǎn)的變形預(yù)測(cè)研究時(shí)，有必要需要綜合考慮多個(gè)測(cè)點(diǎn)的變形協(xié)調(diào)關(guān)系，開(kāi)展多測(cè)點(diǎn)的變形預(yù)測(cè)研究。