杜建廣 沈 尤 常旭輝

(中國電建集團北京勘測設計研究院有限公司, 北京 100024)

0 引言

近年來,我國高層建筑物數(shù)量快速增長,在影響高層建筑物安全建設的眾多因素中,基坑的沉降最為嚴重。建筑物基坑沉降不僅會造成基坑的變形與坍塌,還會對基坑保護區(qū)內其他建構筑物的安全造成影響。為了保障建筑物施工過程中基坑本身及周圍建構筑物的安全,基坑的沉降監(jiān)測尤為重要,能夠對基坑是否發(fā)生變形提供數(shù)據(jù)支撐,為了準確評估建筑物基坑的變形趨勢,控制變形量,基于基坑當前監(jiān)測數(shù)據(jù)對未來一段時間內的變形量進行預測可以有效對基坑的安全性進行評估。

目前常用的變形預測模型有時間序列分析模型、反向傳播(back propagation,BP)神經(jīng)網(wǎng)絡模型、灰色系統(tǒng)分析模型等[1]。相比于其他預測模型,灰色系統(tǒng)GM(1.1)模型的優(yōu)勢在于基于規(guī)律性與樣本量都較小的原始數(shù)據(jù)就可進行變形預測,有著較強的規(guī)律性搜索能力,解決了“小樣本、貧信息”的變形預測問題[2]。BP神經(jīng)網(wǎng)絡模型是一種智能化預測模型,通過自適應學習與計算可以對非線性信號進行較為準確的預測[3]。結合上述兩種預測模型的優(yōu)勢,不少學者進行了研究,文獻[4]有效結合灰色理論有BP神經(jīng)網(wǎng)絡模型,將組合模型應用于實際變形監(jiān)測項目中,得到的預測精度較高;文獻[5]將遺傳算法添加至灰色神經(jīng)網(wǎng)絡中,實現(xiàn)全局最優(yōu)求取的效果。然而,上述算法沒有考慮到非平穩(wěn)時間序列中的趨勢項,直接將灰色預測殘差作為BP神經(jīng)網(wǎng)絡輸入項會導致模型預測精度的受損。

本文在灰色神經(jīng)網(wǎng)絡模型的基礎上,提出了一種改進型組合模型。該組合模型通過經(jīng)驗模態(tài)分解(empirical mode decomposition,EMD)方法實現(xiàn)信號的自適應分解,得到不同頻段分量,對于穩(wěn)定性較高的低頻分量使用差分灰色神經(jīng)網(wǎng)絡模型進行預測,提高了灰色預測殘差的穩(wěn)定性,相比于灰色神經(jīng)網(wǎng)絡模型的預測精度更高;針對穩(wěn)定性較差的高頻分量,使用自相關性強的AR時間序列預測模型進行預測。將本文提出的組合模型應用到某在建建筑物基坑沉降監(jiān)測數(shù)據(jù)預測中,并將預測結果與傳統(tǒng)預測模型預測結果進行對比。

1 EMD分解與重構

EMD方法對信號進行處理的思想是:任何信號都是由一些頻率不同的本征模態(tài)函數(shù)(intrinsic mode function,IMF)組成[6],IMF分量滿足條件為：①上下包絡線均值為零,也就是上下包絡線關于時間軸對稱;②過零點個數(shù)與極值點個數(shù)相差在一個以內。信號的時頻特征通常由相鄰極值點的間隔進行反映，如果信號不存在極值點,可對信號進行多次求導獲得極值點。

使用EMD方法對高層建筑物基坑變形數(shù)據(jù)進行預處理,獲取不同頻率的IMF分量,得到理想分辨率下的低頻分量與高頻分量。針對高頻分量,由于其具有較強的自相關性,采用自回歸模型(autoregressive model,AR)進行預測;針對低頻分量,采用差分灰色神經(jīng)網(wǎng)絡模型進行預測。

2 差分灰色神經(jīng)網(wǎng)絡模型

作為對灰色神經(jīng)網(wǎng)絡的一種改進,差分灰色神經(jīng)網(wǎng)絡模型實現(xiàn)變形數(shù)據(jù)的預測與分析的原理為:對經(jīng)灰色預測得到殘差進行差分化處理,將差分化處理得到的平穩(wěn)序列輸入至BP神經(jīng)網(wǎng)絡中,利用神經(jīng)網(wǎng)絡的自適應學習預測變形數(shù)據(jù)。

2.1 GM(1.1)模型

假設有一時間序列X(0)={x(0)(1),x(0)(2),…,x(0)(n)},x(0)(n)表示第n期觀測值。通過GM(1.1)模型實現(xiàn)時間序列預測的步驟為[7]:

(1)通過對原始觀測數(shù)據(jù)進行累加:X(1)={x(1)(1),x(1)(2),…,x(1)(n)}。

(2)緊鄰均值生成序列:z(1)(k)=1/2(x(1)(k)+x(1)(k-1)),k=2,3,…,n,得到新數(shù)列z(1)={z(1)(2),z(1)(3),…,z(1)(n)}。

(3)通過對新序列構建白化方程,得到式(1)。

(1)

式中,b表示灰色作用量;a表示發(fā)展系數(shù),a可通過最小二乘原理求解得到,其中a=[ab]T=(BB)-1BTY,B與Y可表示為式(2)。

(2)

求解白化方程,得到式(3)。

(3)

通過累減對原序列進行還原,得到:ˉx(0)(k+1)=ˉx(1)(k+1)-x(1)(k),同時得到預測值x(0)(n)。

2.2 灰色預測的殘差差分化

將灰色模型預測殘差作為BP神經(jīng)網(wǎng)絡的輸入數(shù)據(jù)可能會造成網(wǎng)絡無法收斂。因此,需要對GM(1.1)模型預測得到的殘差進行差分處理,避免因殘差序列的不平穩(wěn)造成的模型收斂速度慢以及預測精度降低。

將GM(1.1)模型預測得到的殘差序列記為e(0),對e(0)進行差分處理,得到一階差分項e(1)，見式(4)。

(4)

將平穩(wěn)化殘差序列e(1)(t)輸入至BP神經(jīng)網(wǎng)絡中進行訓練。

2.3 BP神經(jīng)網(wǎng)絡模型

BP是一種多層前饋神經(jīng)網(wǎng)絡,信息向前傳播、誤差逆向傳遞。使用BP神經(jīng)網(wǎng)絡進行數(shù)據(jù)預測不需要事先將輸入與輸出的映射關系進行確認,而是通過純機器學習的方法,在算法中輸入集,算法自動對輸入集學習、訓練,從而得到某種特定的規(guī)則[8]。

對網(wǎng)絡進行訓練的具體步驟為[9]:

(1)網(wǎng)絡初始化。根據(jù)輸入序列X與輸出序列Y對輸入層節(jié)點n、隱含層節(jié)點l與輸出層節(jié)點m進行確定,初始化輸入層與隱含層神經(jīng)元間的連接權值wij以及隱含層與輸出層神經(jīng)元間的連接權值wjk,同時還要初始化的參數(shù)有神經(jīng)元激勵函數(shù)、給定學習速率lr以及隱含層閾值與輸出層閾值a1、a2。

(2)隱含層輸出計算。通過隱含層閾值a1與輸入層、隱含層神經(jīng)元間的連接權值wij對隱含層輸出進行確定，見式(5)。

(5)

式中,Hj表示隱含層輸出;f表示隱含層激勵函數(shù),本文選擇sigmoid作為隱含層激勵函數(shù)，見式(6)。

(6)

(3)輸出層輸出計算。通過輸出層閾值a2、隱含層與輸出層神經(jīng)元間的連接權值wjk以及隱含層輸出Hj對預測輸出進行確定，見式(7)。

(7)

(4)誤差計算。通過期望輸出與預測輸出計算得到預測誤差ek，見式(8)。

(8)

式中,Y表示期望輸出值。

(5)權值更新。使用預測誤差ek對各層之間的權值wij與wjk進行反向更新。

(6)閾值更新。使用預測誤差ek對節(jié)點閾值a1、a2進行反向更新,見式(9)。

(9)

(7)迭代是否結束判斷。若迭代結束,將預測值輸出,若沒有結束,重復上述步驟。

3 AR自回歸模型

在構建自回歸模型前首先計算時間序列的自相關系數(shù)。假設有一時間序列X,該時間序列的相關系數(shù)為r1,r2,…,rn。t時刻觀測值為xt,t+k時刻觀測值為xt+k,二者之間的自相關系數(shù)rk表示為[10]式(10)。

(10)

式中,ˉx表示時間序列X均值。

AR模型能夠得到廣泛使用的原因在于在建模時不需要涉及白噪聲序列解算,只要對線性方程組進行求解即可，見式(11)。

(11)

式中,φ1,φ2,…,φp表示自回歸系數(shù);p表示階數(shù);εt表示與xt不相關的白噪聲序列。

4 基于EMD的組合預測模型

綜合上述算法的優(yōu)勢,本文提出了一種基于EMD的差分灰色神經(jīng)網(wǎng)絡-AR模型,該組合預測模型的具體算法為:

(1)將建筑物基坑沉降監(jiān)測數(shù)據(jù)進行EMD自適應分解,得到一系列不同頻率的IMF分量,其中c1表示低頻分量,d1、d2、d3表示高頻分量。

(2)對于低頻分量c1,使用本文提出的差分灰色神經(jīng)網(wǎng)絡進行預測,得到擬合值與預測值y。

(3)對于高頻分量d1、d2,使用AR模型進行預測,分別得到擬合值與預測值y1、y2。

(4)對于高頻分量d3,將其視為噪聲項予以剔除。

(5)綜上可以,建筑物基坑變形監(jiān)測數(shù)據(jù)X可表示為X=c1+d1+d2+d3,c1、d1、d2的預測值分別為y,y1、y2,將各自預測值進行累加即可得到建筑物基坑沉降預測值Y=y+y1+y2。

圖1為本文提出算法的技術路線。

圖1 本文提出算法的技術路線

5 實例分析

5.1 工程概況

本文以北京市朝陽區(qū)某在建居民樓項目的基坑沉降監(jiān)測數(shù)據(jù)為例進行算法的驗證,該基坑深度約6 m。為保障項目在建設過程中的安全,經(jīng)實地踏勘,按照三級基坑控制的思想布設7個監(jiān)測點,分別為J1～J7。使用二等水準按閉合路線,結合成本與效率的考慮對7個監(jiān)測點實施復測。

5.2 變形預測與分析

選取該基坑監(jiān)測點J4點連續(xù)40期累計沉降監(jiān)測數(shù)據(jù)為例進行變形預測,其中前30期監(jiān)測數(shù)據(jù)用于建模,后10期監(jiān)測數(shù)據(jù)作為預測數(shù)據(jù)期望值,對模型預測的效果進行檢驗。

使用EMD分解前30期監(jiān)測數(shù)據(jù)序列。使用差分灰色神經(jīng)網(wǎng)絡對低頻分量c1進行預測,其中參數(shù)設置為:輸入層、隱含層、輸出層神經(jīng)元個數(shù)分別為4、9、1,分別選擇tan sig、purelin作為隱含層與輸出層的神經(jīng)元傳遞函數(shù)。選擇最大最小值法對訓練樣本進行歸一化處理,以避免訓練結果受數(shù)據(jù)本身影響。

本文通過計算排列熵(permutation entropy,PE)區(qū)分高頻分量與低頻分量,排列熵是一種對時間序列隨機性程度進行檢測的方法,信號的復雜程度越高,排列熵值就會越大；反之,排列熵值越小。多尺度排列熵(multi-scale permutation entropy,MPE)是對排列熵方法的一種改進,它能夠更加準確地對信號的隨機性進行分析。多尺度排列熵的基本思想是對時間序列的多尺度粗?；?將多尺度閾值熵的閾值設置為0.6,也就是說,多尺度排列熵的值大于0.6的分量為高頻分量,多尺度排列熵的值小于0.6的分量為低頻分量。

計算經(jīng)EMD分解沉降監(jiān)測數(shù)據(jù)得到的IMF分量的MPE值,區(qū)分出高頻IMF分量與低頻IMF分量,各分量的MPE值如表1所示。

表1 各IMF分量多尺度排列熵

從表1中可以看到,前3個IMF分量的MPE值均大于0.6,故可將前3個IMF分量定為高頻分量。通過頻譜計算得知,第1個IMF分量主頻位于0.25～0.5 Hz頻段,振幅小于1 mm,呈現(xiàn)噪聲特征,可將前第1個IMF分量認定為高頻噪聲分量,予以剔除;第2、第3個IMF分量的主頻位于0.1～0.25 Hz頻段,可認為此頻段反映了建筑物基坑的沉降特性,將第2、第3個IMF分量認定為高頻分量d1、d2;將剩余分量認定為低頻分量c1。

使用AR預測模型對高頻分量d1、d2進行建模,使用PACF定階段法確定d1、d2分量均為2階。

使用差分灰色神經(jīng)網(wǎng)絡對低頻分量進行預測得到結果如圖2所示;使用AR模型對高頻分量進行預測得到結果如圖3所示,同時給出了預測結果誤差百分比。

對上述各分量預測結果進行重構,最終得到建筑物基坑沉降擬合值與預測值。將GM(1.1)模型、灰色神經(jīng)網(wǎng)絡模型與本文提出預測模型預測結果進行對比,如表2所示為利用3種預測模型對建筑物基坑J4點第31～40期沉降預測結果。

圖2 低頻分量預測結果

圖3 高頻分量預測結果

表2 3種模型預測結果單位:mm

通過表2計算可知,利用GM(1.1)模型進行建筑物基坑沉降監(jiān)測數(shù)據(jù)預測結果的殘差平均絕對值為0.794 mm;利用差分灰色神經(jīng)網(wǎng)絡進行建筑物基坑沉降監(jiān)測數(shù)據(jù)預測結果的殘差平均絕對值為0.547 mm;利用本文提出的基于EMD的組合預測模型進行建筑物基坑沉降監(jiān)測數(shù)據(jù)預測結果的殘差平均絕對值為0.518 mm。相比于前兩種預測模型,本文提出的預測模型預測效果更好,驗證了本文算法的可靠性與優(yōu)越性。傳統(tǒng)變形預測模型僅考慮單因子變化,不適用于具有非線性、非穩(wěn)定性特征的建筑物基坑變形預測的情況。充分考慮建筑物基坑非線性、非穩(wěn)定性的變形特征,提出了基于EMD的組合預測模型并應用到建筑物基坑變形監(jiān)測中,通過差分灰色神經(jīng)網(wǎng)絡模型消除了由于灰色預測殘差不平穩(wěn)而導致的預測誤差,同時發(fā)揮了AR時間序列模型自相關強的優(yōu)勢[12-15]。

6 結束語

針對傳統(tǒng)時間序列預測模型在針對非平穩(wěn)序列預測中的劣勢,本文結合多種預測模型優(yōu)勢,提出了一種基于EMD方法的改進預測模型。該預測模型有效消除了灰色預測殘差不穩(wěn)定的問題,通過EMD方法對基坑沉降監(jiān)測數(shù)據(jù)進行分解,得到不同頻段的IMF分量。針對低頻段穩(wěn)定性更好的序列,使用差分神經(jīng)網(wǎng)絡模型進行預測;針對高頻段不穩(wěn)定的序列,使用AR模型進行預測,充分利用了AR模型自相關強的優(yōu)勢。將不同頻段分量經(jīng)預測模型預測得到的序列進行重構得到最終預測結果。

將灰色GM(1.1)模型、差分灰色神經(jīng)網(wǎng)絡模型與本文提出組合模型應用于北京市朝陽區(qū)某在建小區(qū)建筑物基坑沉降監(jiān)測數(shù)據(jù)預測中,通過對比3種模型預測結果,得到改進后的模型預測精度更高,能夠會建筑物基坑沉降的變形趨勢做出更加準確的預報。