朱軍桃,沙 穎,李 濤
(1.桂林理工大學(xué) 土木與建筑工程學(xué)院,廣西 桂林 541004;2.中交四航工程研究院有限公司,廣東 廣州 510288)
沉降監(jiān)測(cè)貫穿著工程建筑物的設(shè)計(jì)期、施工期和運(yùn)營期整個(gè)過程,無論是公路工程、建筑工程還是水利工程,地基沉降安全通常被視為工程成敗的關(guān)鍵.在目前的沉降預(yù)測(cè)中,通常選取累積沉降樣本作為訓(xùn)練樣本,但究竟是采用累積沉降樣本還是間隔沉降樣本才能使預(yù)測(cè)結(jié)果相對(duì)較好,尚無說明.本文利用同一沉降樣本中的累積樣本和間隔樣本,用小波分析、神經(jīng)網(wǎng)絡(luò)和小波神經(jīng)網(wǎng)絡(luò)模型對(duì)這些樣本分別進(jìn)行了分析預(yù)測(cè),得出了相應(yīng)結(jié)論.
含有噪聲的變形監(jiān)測(cè)數(shù)據(jù)列可建模如下[1]:
s(n)=f(n)+e(i),i=1,2…n
(1)
其中,s(n)是監(jiān)測(cè)數(shù)據(jù),f(n)是變形數(shù)據(jù),它包含實(shí)際變形量和確定性噪聲(如系統(tǒng)誤差),e(i)是隨機(jī)噪聲,e(i)~N(0,σ2).在小波分析中,變形信號(hào)表現(xiàn)為低頻信號(hào)或是一些比較平穩(wěn)的信號(hào),而噪聲信號(hào)主要被包含在小波分解的高頻層中,所以以門限閾值形式對(duì)小波系數(shù)進(jìn)行處理,減少噪聲部分的值,再通過重構(gòu)恢復(fù)信號(hào),達(dá)到消噪的目的.以累積沉降樣本為例,具體的流程見圖1.
圖1 累積沉降量小波分析示意圖Fig.1 Analysis sketch map of accumulated settlement based on wavelet model
yn=(ym,ym+1,…yn-1),
(2)
在(2)式中,yn表示第n期沉降累積或間隔量,經(jīng)過分析研究和大量試算,樣本構(gòu)造如下:
y13=f(y1,y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12),
y14=f(y2,y3,y4,y5,y6,y7,y8,y9,y10,y11,y12,y13)…
具體步驟如下:
(1) 輸入層元素取12,分別為累積或間隔沉降量yn-1,yn-2,yn-3,yn-4,yn-5,yn-6,yn-7,yn-8,yn-9,yn-10,yn-11,yn-12;
(2)隱含層節(jié)點(diǎn)數(shù)取20[4];
(3)輸出層元素取1,為累積/間隔沉降量yn.
小波與神經(jīng)網(wǎng)絡(luò)的結(jié)合通常有2種方式[5]:一種是輔助式結(jié)合,另一種是嵌套式結(jié)合.
首先,用小波分析將樣本數(shù)據(jù)進(jìn)行預(yù)處理,然后將預(yù)處理過的數(shù)據(jù)作為神經(jīng)網(wǎng)絡(luò)的輸入向量,再用神經(jīng)網(wǎng)絡(luò)進(jìn)行處理,得出最終結(jié)果,具體流程如圖2所示.
圖 2 小波神經(jīng)網(wǎng)絡(luò)運(yùn)行圖Fig.2 Operation map of wavelet neural network
小波與神經(jīng)網(wǎng)絡(luò)的嵌套式結(jié)合,用小波函數(shù)直接代替隱含層函數(shù),基于多分辨率分析理論的小波神經(jīng)網(wǎng)絡(luò)本結(jié)構(gòu)如圖3所示,其中,xi(i=1,2,…,n) 為輸入樣本,ψk(k=1,2,…,l)為小波基函數(shù),fj(j=1,2,…,m)為網(wǎng)絡(luò)的輸出,用ωki表示隱含層第k個(gè)神經(jīng)元與輸入層第i個(gè)神經(jīng)元間的連接權(quán)值,ωjk表示輸出層第j個(gè)神經(jīng)元與隱含層第k個(gè)神經(jīng)元間的連接權(quán)值,根據(jù)所選取的小波基函數(shù)的連續(xù)性的不同,可以將該模型分為連續(xù)參數(shù)的小波神經(jīng)網(wǎng)絡(luò)和基于小波框架的小波神經(jīng)網(wǎng)絡(luò).
(1)連續(xù)參數(shù)的小波神經(jīng)網(wǎng)絡(luò)
(3)
含層節(jié)點(diǎn)激勵(lì)函數(shù)一般用Morlet和Mexican Hat等連續(xù)小波.
圖3 隱含層為小波函數(shù)的神經(jīng)網(wǎng)絡(luò)Fig.3 Neural network with wavelet function in hidden layer
(2)離散小波神經(jīng)網(wǎng)絡(luò)
(4)
隱含層節(jié)點(diǎn)的激勵(lì)函數(shù)一般采用離散小波或小波框架.
由小波函數(shù)直接代替隱含層函數(shù)構(gòu)成的神經(jīng)網(wǎng)絡(luò)模型,其訓(xùn)練與學(xué)習(xí)可以采取與傳統(tǒng)的神經(jīng)網(wǎng)絡(luò)完全相同的方法進(jìn)行.
利用Matlab軟件的數(shù)值計(jì)算和仿真分析功能[6],編寫與調(diào)用3個(gè)模型程序,輸出各種模型的結(jié)果圖和數(shù)據(jù).
本文數(shù)據(jù)來自于某大廈沉降數(shù)據(jù)中的1#點(diǎn)的累積沉降量和間隔沉降量樣本,具體見表1.
表1 沉降觀測(cè)成果累積沉降和間隔沉降匯總表Tab.1 Observation results of accumulated settlement and interval settlement
表2 累積沉降樣本和間隔沉降樣本的小波去噪結(jié)果Tab.2 Processing of accumulated sample data and interval sample data based on wavelet model
從以上小波去噪的結(jié)果可以看出,17期的沉降去噪結(jié)果中,有12期的結(jié)果是間隔沉降樣本處理結(jié)果優(yōu)于累積沉降樣本的.
應(yīng)用累積沉降數(shù)據(jù)樣本數(shù)據(jù)計(jì)算,本模型本項(xiàng)目采用lr為0.1、mc為0.9(lr為學(xué)習(xí)速率,mc為動(dòng)量系數(shù))[7].這里對(duì)此方案進(jìn)行訓(xùn)練,使樣本集的訓(xùn)練精度目標(biāo)為0.001 mm,訓(xùn)練2 000次,經(jīng)神經(jīng)網(wǎng)絡(luò)訓(xùn)練后結(jié)果與原數(shù)據(jù)的對(duì)比見表3.
表3 基于神經(jīng)網(wǎng)絡(luò)的累積沉降樣本和間隔沉降樣本的數(shù)據(jù)處理結(jié)果Tab.3 Processing of accumulated sample data and interval sample data based on neural network
圖4 累積樣本的神經(jīng)網(wǎng)絡(luò)運(yùn)行圖 Fig.4 Operation map of neural network of accumulated sample
圖5 間隔樣本的神經(jīng)網(wǎng)絡(luò)運(yùn)行圖Fig.5 Operation map of neural network of interval sample
從以上應(yīng)用神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)的沉降結(jié)果中可以看到,5組預(yù)測(cè)值中有4組的應(yīng)用間隔沉降樣本的結(jié)果優(yōu)于應(yīng)用累積沉降樣本.
下面將經(jīng)過小波分析處理后的沉降數(shù)據(jù)代入神經(jīng)網(wǎng)絡(luò)模型中,神經(jīng)網(wǎng)絡(luò)參數(shù)采用lr為0.1、mc為0.9.首先是累積沉降樣本的訓(xùn)練,使樣本集的訓(xùn)練精度目標(biāo)為0.01 mm,訓(xùn)練2 000次,經(jīng)神經(jīng)網(wǎng)絡(luò)訓(xùn)練后結(jié)果與原數(shù)據(jù)的對(duì)比如表4所示.
表4 基于小波神經(jīng)網(wǎng)絡(luò)的累積沉降樣本和間隔沉降樣本的數(shù)據(jù)處理結(jié)果Tab.4 Processing of accumulated sample data and interval sample data based on wavelet neural network
圖6 小波神經(jīng)網(wǎng)絡(luò)的累積樣本的神經(jīng)網(wǎng)絡(luò)運(yùn)行Fig.6 Operation map of wavelet neural network of accumulated sample
圖7 小波神經(jīng)網(wǎng)絡(luò)的間隔樣本的神經(jīng)網(wǎng)絡(luò)運(yùn)行圖Fig.7 Operation map of wavelet neural network of interval sample
從以上2種樣本的小波神經(jīng)網(wǎng)絡(luò)的預(yù)測(cè)結(jié)果可以看出,5期的預(yù)測(cè)結(jié)果都是間隔沉降樣本的預(yù)測(cè)結(jié)果優(yōu)于累積沉降樣本.
綜合以上幾種模型分析的結(jié)果,采用間隔沉降樣本預(yù)測(cè)在多數(shù)情況下可以取得較好的預(yù)測(cè)效果.本文演示了一個(gè)沉降實(shí)例的預(yù)測(cè)結(jié)果比較,作者曾對(duì)多個(gè)沉降實(shí)例進(jìn)行試驗(yàn),綜合結(jié)果顯示,對(duì)于處于快速下沉狀態(tài)或者是由于其他原因的影響下沉加劇的監(jiān)測(cè)點(diǎn)來說,采用間隔沉降樣本的預(yù)測(cè)結(jié)果基本上都優(yōu)于累積沉降樣本的預(yù)測(cè)結(jié)果.通過本文還可以看出,采用小波神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)沉降相對(duì)于其他2種模型,輸出結(jié)果的誤差精度更易收斂,輸出結(jié)果的整體性也更好.
參考文獻(xiàn):
[1] 黃聲亨,尹暉,蔣征.變形監(jiān)測(cè)數(shù)據(jù)處理[M].武漢:武漢大學(xué)出版社,2002.
[2] 張立明.人工神經(jīng)網(wǎng)絡(luò)的模型及應(yīng)用[M].上海:復(fù)旦大學(xué)出版社,1999.
[3] 岳榮花.小波神經(jīng)網(wǎng)絡(luò)在沉降預(yù)測(cè)中的應(yīng)用研究[D].南京:河海大學(xué),2007.
[4] 周紅曉,蔡俊,任德官.一種優(yōu)化多層前饋神經(jīng)網(wǎng)絡(luò)中隱節(jié)點(diǎn)數(shù)的算法[J].浙江師范大學(xué)學(xué)報(bào):自然科學(xué)版,2002,25(3):268-271.
[5] 任俊玲,郭軍.基于小波理論的神經(jīng)網(wǎng)絡(luò)模型構(gòu)造[J].電腦開發(fā)與應(yīng)用,2004(8):19-21.
[6] 周開利,康耀紅.神經(jīng)網(wǎng)絡(luò)模型及其Matlab仿真程序設(shè)計(jì)[M].北京:清華大學(xué)出版社,2005:7.
[7] 張正祿,汪宏晨.滑坡變形分析與預(yù)報(bào)的新方法[J].武漢大學(xué)學(xué)報(bào):信息科學(xué)版,2009, 34(12):1387-1389.