基于卷積神經(jīng)網(wǎng)絡(luò)的拋石防波堤沉降預(yù)測(cè)

陳少青，陳文強(qiáng)，程 林，包勝利

（1.天津市北洋水運(yùn)水利勘察設(shè)計(jì)研究院有限公司，天津 300460；2.天津理工大學(xué) 管理學(xué)院，天津 300384；3.中交天航港灣建設(shè)工程有限公司，天津 300450）

拋石防波堤是通過(guò)拋填粗細(xì)石料的方法修建的一種斜坡式防波堤，能夠起到阻斷波浪入侵、圍護(hù)港池、維持水面平穩(wěn)的作用，因其對(duì)地基承載力要求不高且結(jié)構(gòu)簡(jiǎn)單、施工方便，比較廣泛的應(yīng)用于港口建設(shè)中。防波堤沉降量的控制是影響其防護(hù)功能正常發(fā)揮的重要因素，目前國(guó)內(nèi)外學(xué)者對(duì)防波堤沉降的研究主要集中在地基固結(jié)變形和堤身沉降監(jiān)測(cè)兩方面。崔廣強(qiáng)等［1］根據(jù)實(shí)際觀測(cè)資料，運(yùn)用灰色Ver?hulst模型預(yù)測(cè)了堤身各時(shí)刻的沉降量；葉劍紅等［2］基于有限元法和時(shí)間積分法對(duì)大型復(fù)合防波堤固結(jié)過(guò)程進(jìn)行了仿真模擬；王偉等［3］根據(jù)拋石防波堤拋填施工和GPS監(jiān)測(cè)技術(shù)特點(diǎn)，提出了一種實(shí)時(shí)監(jiān)測(cè)防波堤沉降變形的新方法；黃其歡等［4］提出了一種基于PSI技術(shù)的防波堤沉降監(jiān)測(cè)方法。可見(jiàn)，關(guān)于防波堤沉降的預(yù)測(cè)方法較為豐富，現(xiàn)階段主要采用的是根據(jù)沉降分布規(guī)律提出的經(jīng)驗(yàn)公式法；但由于防坡堤施工工序較為復(fù)雜，沉降變形具有多變性，經(jīng)驗(yàn)公式法在沉降參數(shù)的動(dòng)態(tài)調(diào)整上有很大的局限性。

本文依托天津港大沽口港區(qū)東、北防波堤潛堤段工程，通過(guò)分析拋石防波堤施工過(guò)程中影響沉降的因素，提出拋石防波堤施工沉降的深度學(xué)習(xí)預(yù)測(cè)模型；同時(shí)利用現(xiàn)場(chǎng)監(jiān)測(cè)數(shù)據(jù)對(duì)深度學(xué)習(xí)模型進(jìn)行訓(xùn)練，有效提升拋石防波堤沉降預(yù)測(cè)的準(zhǔn)確性。

1 研究方法

為提升預(yù)測(cè)的準(zhǔn)確性，本文基于卷積神經(jīng)網(wǎng)絡(luò)(Convolutional Neural Networks,CNN)構(gòu)建深度網(wǎng)絡(luò)模型。卷積神經(jīng)網(wǎng)絡(luò)包含卷積層、池化層、權(quán)重及偏置[5]，將n組歸一化后的參數(shù)代入模型中，通過(guò)卷積層和池化層的雙重濾網(wǎng)過(guò)濾系統(tǒng)，實(shí)現(xiàn)對(duì)n個(gè)參數(shù)的訓(xùn)練，通過(guò)不斷優(yōu)化和調(diào)整原有池化層、權(quán)重及偏置參數(shù)，判斷代價(jià)函數(shù)（平方誤差）CT的值是否無(wú)限趨近或等于0。若代價(jià)函數(shù)CT無(wú)限趨近或等于0，則輸出的結(jié)果為最優(yōu)，完成拋石防波堤沉降的預(yù)測(cè)。

1.1 卷積操作

由式（1）得到的卷積值還需要加上一個(gè)不依賴i、j變化的數(shù)bFk。

根據(jù)式（2）得到加權(quán)輸入矩陣，作為卷積層一個(gè)子層。將Sigmoid函數(shù)作為激活函數(shù)，卷積層輸出矩陣為

1.2 池化操作

卷積神經(jīng)網(wǎng)絡(luò)中有用于壓縮卷積層信息的池化層，通過(guò)最大池化法，將卷積層中2個(gè)×2個(gè)神經(jīng)單元壓縮為1個(gè)神經(jīng)單元，這些壓縮后的神經(jīng)單元集合就形成了池化層，池化層輸入和輸出分別為

計(jì)算輸出層輸入和輸出

式中：——輸出層第n個(gè)神經(jīng)單元給池化層神經(jīng)單元的輸出（k=1，2，3；i=1，2；j=1，2）分配的權(quán)重；

——輸出層第n個(gè)神經(jīng)單元的偏置。

假設(shè)輸出層神經(jīng)單元的3個(gè)輸出對(duì)應(yīng)的正解分別為t1、t2、t3

將全體學(xué)習(xí)數(shù)據(jù)的平方誤差綜合就是代價(jià)函數(shù)CT，通過(guò)求解CT最小達(dá)到最優(yōu)模型解的權(quán)重及偏差，將最優(yōu)模型運(yùn)用至測(cè)試樣本，從而達(dá)到沉降的預(yù)測(cè)效果。見(jiàn)圖1。

圖1 卷積和池化過(guò)程

2 工程應(yīng)用

2.1 工程概況

天津港大沽口港區(qū)東、北防波堤潛堤段工程位于大沽口港區(qū)大沽沙航道口門附近，東、北防波堤均為永久性海堤，設(shè)計(jì)使用年限為50 a，結(jié)構(gòu)安全等級(jí)為二級(jí)，其中東防波堤潛堤段全長(zhǎng)1 500 m，北防波堤潛堤段全長(zhǎng)2 350 m。堤身結(jié)構(gòu)為復(fù)合布+編織布大型充填袋裝砂棱體；內(nèi)外坡均為1∶1.75；護(hù)面結(jié)構(gòu)厚度1.1 m（堤頭段1.6 m），為規(guī)格300~500 kg拋石墊層+6 t扭王字塊人工護(hù)面塊體；反濾結(jié)構(gòu)為600 g/m2復(fù)合土工布+厚度0.7 m規(guī)格10~40 kg的拋石墊層；護(hù)底結(jié)構(gòu)為砂肋軟體排+厚度1.0 m規(guī)格60~100 k g的拋石；地基鋪設(shè)砂被墊層，打設(shè)塑料排水板，間距1.1 m（局部1.0 m），底標(biāo)高-16 m。

2.2 監(jiān)測(cè)數(shù)據(jù)

為了充分利用現(xiàn)有監(jiān)測(cè)數(shù)據(jù)，預(yù)測(cè)防波堤施工期間沉降量，實(shí)時(shí)掌握防波堤堤身沉降情況，以動(dòng)態(tài)指導(dǎo)和控制施工，施工期間進(jìn)行了堤身沉降、深層土體水平位移和孔隙水壓力觀測(cè)。將堤身高度、孔隙水壓力、深層水平位移作為卷積神經(jīng)網(wǎng)絡(luò)的輸入層參數(shù)，將拋石防波堤沉降量作為卷積神經(jīng)網(wǎng)絡(luò)的輸出層參數(shù)。

對(duì)各樣本中的堤身高度、孔隙水壓力、深層水平位移指標(biāo)參數(shù)進(jìn)行歸一化處理。

式中：xij——第i個(gè)參數(shù)的第j個(gè)參數(shù)；

xmax——參數(shù)的最大值；

xmin——參數(shù)的最小值。

卷積神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)樣本及參數(shù)歸一化結(jié)果見(jiàn)表1。

表1 學(xué)習(xí)樣本和歸一化樣本

2.3 卷積神經(jīng)網(wǎng)絡(luò)（CNN）預(yù)測(cè)結(jié)果

將1～13組樣本作為學(xué)習(xí)樣本，基于式（1）～（8）進(jìn)行訓(xùn)練學(xué)習(xí)，求解平方誤差CT并進(jìn)行優(yōu)化權(quán)重和偏置，使CT最小，最終得到最小平方誤差CT=0.002，將最優(yōu)過(guò)濾器、權(quán)重和偏置代入到模型中，作為訓(xùn)練好的卷積神經(jīng)網(wǎng)絡(luò)。

將14～17組樣本作為檢驗(yàn)樣本，代入已經(jīng)訓(xùn)練好的卷積神經(jīng)網(wǎng)絡(luò)進(jìn)行檢驗(yàn)，輸出的沉降量結(jié)果分別為362、368、374、928 mm，與原始數(shù)據(jù)進(jìn)行比較，發(fā)現(xiàn)與實(shí)際結(jié)果基本一致，最大誤差為1.3%，說(shuō)明本文建立的卷積神經(jīng)網(wǎng)絡(luò)方法滿足要求。堤施工期間的沉降量，預(yù)測(cè)模型根據(jù)輸入?yún)?shù)實(shí)測(cè)數(shù)據(jù)能動(dòng)態(tài)指導(dǎo)和控制施工進(jìn)度。□■

3 結(jié)論

1）卷積神經(jīng)網(wǎng)絡(luò)具有降采樣的特點(diǎn)，預(yù)測(cè)精度較高。本文將卷積神經(jīng)網(wǎng)絡(luò)應(yīng)用在拋石防波堤施工沉降預(yù)測(cè)中，在防坡堤施工沉降預(yù)測(cè)中具有較高的可靠性。

2）卷積神經(jīng)網(wǎng)絡(luò)模型選擇拋石防波堤堤身高度、孔隙水壓力和深層水平位移3個(gè)輸入?yún)?shù)來(lái)估算防波