王雪妮， 韓國鋒

(1. 楊凌職業(yè)技術(shù)學(xué)院， 陜西 咸陽 712100； 2. 陜西鐵路工程職業(yè)技術(shù)學(xué)院， 陜西 渭南 714000)

趨勢項分離預(yù)測模型及重標(biāo)度極差分析在深基坑變形預(yù)測中的應(yīng)用研究

王雪妮1， 韓國鋒2

(1. 楊凌職業(yè)技術(shù)學(xué)院， 陜西 咸陽 712100； 2. 陜西鐵路工程職業(yè)技術(shù)學(xué)院， 陜西 渭南 714000)

為解決基坑變形預(yù)測精度低的問題，利用小波去噪和卡爾曼濾波對基坑變形序列進(jìn)行去噪處理，分離趨勢項及誤差項，并利用支持向量機(jī)和BP神經(jīng)網(wǎng)絡(luò)分別對趨勢項和誤差項進(jìn)行預(yù)測，以掌握基坑的變形規(guī)律及發(fā)展趨勢； 同時，采用重標(biāo)度極差分析(R/S分析)對基坑的變形趨勢進(jìn)行判斷，以驗證變形預(yù)測的可靠性。根據(jù)實例檢驗，得出小波去噪的去噪效果較好，且預(yù)測結(jié)果的相對誤差均值為1.03%，方差值為0.083，預(yù)測精度較高；基坑的變形序列與速率序列均具有持續(xù)增長的趨勢特征，與變形預(yù)測結(jié)果一致，驗證了預(yù)測思路的有效性。

深基坑； 變形預(yù)測； 去噪分析； 支持向量機(jī)； R/S分析； 趨勢判斷

0 引言

隨著近年來城市建設(shè)的快速發(fā)展，深基坑的應(yīng)用越來越廣泛，但由于基坑施工對地質(zhì)環(huán)境的改變，使得基坑開挖以后會出現(xiàn)不同程度的變形，而變形會對基坑自身、周邊建筑及關(guān)系帶來影響，甚至誘發(fā)安全問題，如廣州海珠區(qū)在2005年發(fā)生了一起深基坑坍塌的安全事故，杭州市在2008年也發(fā)生了一起深基坑的坍塌事故； 因此，對基坑的變形進(jìn)行監(jiān)測及研究，且根據(jù)監(jiān)測數(shù)據(jù)進(jìn)行分析及預(yù)測，已成為判斷基坑穩(wěn)定性及避免周邊建筑被破壞的重要手段[1-3]。

近些年，相關(guān)學(xué)者將支持向量機(jī)和神經(jīng)網(wǎng)絡(luò)引入到基坑的變形預(yù)測中。例如： 渠孟飛等[4]、林楠等[5]利用支持向量機(jī)對深基坑的變形進(jìn)行了預(yù)測研究，經(jīng)實例檢驗，得出該模型在深基坑中的預(yù)測精度較高，能有效判斷基坑的發(fā)展趨勢，達(dá)到對基坑安全監(jiān)控的目的；張世豪等[2]、曾暉等[6]利用BP神經(jīng)網(wǎng)絡(luò)對基坑的變形進(jìn)行預(yù)測，得到預(yù)測結(jié)果與實測值吻合度較好，且具有精度與效率兼顧的特點。然而基坑變形數(shù)據(jù)難免會受到環(huán)境因素的影響，使得監(jiān)測值與真實變形值存在一定差異，小波去噪和卡爾曼濾波方法能從監(jiān)測數(shù)據(jù)中提取出真實的變形數(shù)據(jù)[1,7]，具有較好的適用性。

上述研究雖然取得了一定的成果，但多是將趨勢分析及變形預(yù)測獨立進(jìn)行研究，缺少兩者的綜合應(yīng)用； 因此，采用小波去噪和卡爾曼濾波對基坑變形數(shù)據(jù)進(jìn)行去噪處理，將最優(yōu)去噪結(jié)果作為趨勢項和誤差項的分離依據(jù)，并利用支持向量機(jī)和BP神經(jīng)網(wǎng)絡(luò)分別對兩序列進(jìn)行預(yù)測，由兩者的預(yù)測結(jié)果綜合得到預(yù)測值。同時，為驗證預(yù)測結(jié)果的準(zhǔn)確性，再將R/S分析方法引入到基坑的變形趨勢判斷中[8]，以期達(dá)到對基坑變形的綜合變形趨勢判斷。

1 基本原理

1.1 支持向量機(jī)

支持向量機(jī)是一種基于統(tǒng)計學(xué)開發(fā)的數(shù)據(jù)挖掘方法，其基本預(yù)測過程是采用最小化的結(jié)構(gòu)風(fēng)險準(zhǔn)則，通過非線性映射函數(shù)將輸入信息向高維空間進(jìn)行映射，建立起輸入信息與輸出信息間的關(guān)系，能有效解決非線性的預(yù)測問題，具有運算速度快、精度高等優(yōu)點。

在支持向量機(jī)的應(yīng)用過程中，首先引入超平面的概念，它指的是輸入信息與輸出信息在數(shù)量關(guān)系上的回歸線。支持向量機(jī)一元線性回歸的超平面可表示為

y=b+wx[9]。

(1)

式中： y為輸出信息； b為截距； w為斜率； x為輸入信息。

二元線性回歸的超平面則可表示為

y=b+w1x1+w2x2。

(2)

進(jìn)一步可將n個輸入變量的超平面表示為

(3)

式(3)中的W為斜率向量，可表示為

WT={w1,w2,…,wn}。

(4)

同時，由于在n維空間中可能找不到較為合適的超平面，可利用核函數(shù)將樣本間接地映射到高維空間，而核函數(shù)有若干形式，如多項式、線性、徑向基及Sigmoid。結(jié)合實例特點，將核函數(shù)選定為Sigmoid函數(shù)。在支持向量機(jī)的應(yīng)用過程中，利用損失函數(shù)對平面參數(shù)進(jìn)行估計，采用的損失函數(shù)為ε(不敏感損失函數(shù))，且在輸入變量較多的情況下，ε(不敏感損失函數(shù))會呈現(xiàn)“管道”現(xiàn)象。“管道”內(nèi)的點誤差可忽略，外側(cè)的點即為支持向量。

可將支持向量機(jī)的目標(biāo)函數(shù)定義為

(5)

在式(5)的應(yīng)用過程中，需滿足的約束條件如下：

(b+WTXi)-yi≤ε+ξi(i=1,2,…,m)；

(6)

(7)

在支持向量機(jī)的預(yù)測過程中，采用滾動預(yù)測的思想，以提高預(yù)測結(jié)果的準(zhǔn)確性[9]。滾動預(yù)測思想認(rèn)為預(yù)測節(jié)點xi的變形值與前n個節(jié)點相關(guān)，因此將序列{x0,x1,…,xn-1}作為輸入，以xn為輸出，以此類推進(jìn)行訓(xùn)練； 而在驗證或預(yù)測樣本中，將預(yù)測值替換為對應(yīng)時間節(jié)點的值，且結(jié)合實例特點，將n值確定為10。該方法能將有限樣本進(jìn)行較好的推廣預(yù)測，是一種較好的非線性識別方法。

1.2 BP神經(jīng)網(wǎng)絡(luò)

BP神經(jīng)網(wǎng)絡(luò)是一種多層前饋神經(jīng)網(wǎng)絡(luò)，具有較強(qiáng)的非線性映射能力，結(jié)構(gòu)圖如圖1所示。

圖1 BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)圖

在BP神經(jīng)網(wǎng)絡(luò)的參數(shù)設(shè)定過程中，選定train函數(shù)為訓(xùn)練函數(shù)，學(xué)習(xí)率設(shè)定為0.2，節(jié)點的連接權(quán)值區(qū)間設(shè)定為(0,1)，期望精度為0.01，最大迭代次數(shù)為 1 000次。同時，隱層節(jié)點數(shù)對預(yù)測精度及收斂速度具有較大的影響，結(jié)合文獻(xiàn)[10]的研究，采用試算法確定最優(yōu)隱層節(jié)點數(shù)，試算區(qū)間為5～12。

將BP神經(jīng)網(wǎng)絡(luò)的預(yù)測流程[11]分述如下：

1)根據(jù)網(wǎng)絡(luò)結(jié)構(gòu)的設(shè)定，對各節(jié)點的權(quán)值進(jìn)行初始化賦值。

2)根據(jù)訓(xùn)練樣本進(jìn)行訓(xùn)練，得到輸出值Tk′。

3)若輸出值滿足期望要求，則停止訓(xùn)練；若不滿足期望要求，則采用梯度法進(jìn)行權(quán)值修正，且將輸出誤差的函數(shù)表示為

(8)

權(quán)值修正函數(shù)表示為

wij′=wij+Δwij=wij+ηδqi。

(9)

式(8)—(9)中：Tk為輸出目標(biāo)值；wij為原有權(quán)值；wij′為修正權(quán)值；η為學(xué)習(xí)因子；δqi為誤差修正系數(shù)。

4)在規(guī)定最大迭代次數(shù)內(nèi)，若訓(xùn)練誤差達(dá)到收斂誤差，則停止迭代；若未達(dá)到收斂誤差，則輸出最后預(yù)測，并對模型參數(shù)進(jìn)行重新設(shè)定，再重復(fù)進(jìn)行訓(xùn)練，直到誤差達(dá)到收斂誤差。

同時，本文BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)采用的是遞推模型結(jié)構(gòu)，是將xi節(jié)點的前n個節(jié)點的變形值作為輸入樣本，且n值設(shè)定為10。該模型具有滾動預(yù)測的特點，推廣性較強(qiáng)。

1.3 R/S分析

R/S分析首先由Hurst將其應(yīng)用于尼羅河的水文分析中，可對時間序列的長期相關(guān)性及分形特征進(jìn)行判斷，進(jìn)而確定時間序列的趨勢項。

若將隧道的時間序列表示為{Ni}(i=1,2,…,N)，將其劃分為A個長為n的子序列，設(shè)定每個子序列為Ia(a=1,2,…,A)，且子序列中的每個元素可表示為Nk,a，則子序列Ia的平均值可表示為

(10)

式中ea為子序列的均值。

根據(jù)子序列的基本元素和均值，計算各子序列的累計離差，即

(11)

進(jìn)一步計算得到每個子序列的極差

(12)

同時，每個子序列的標(biāo)準(zhǔn)差為

(13)

經(jīng)過前述對極差和標(biāo)準(zhǔn)差的求解，可得到各子序列的重標(biāo)度極差為

(14)

對于特定n值時，可將其重標(biāo)度極差的均值表示為

(15)

當(dāng)n值不同時，可得到多個(n，(R/S)n)散點，且根據(jù)Mandelbrot的研究成果，lgn和lg (R/S)n之間存在線性關(guān)系如下。

lg (R/S)n=lgC+H·lgn。

(16)

式中：H為Hurst值；C為常數(shù)。

根據(jù)Matlab的cftool工具箱對若干(n，(R/S)n)節(jié)點進(jìn)行最小二乘擬合，以得到Hurst指數(shù)，進(jìn)而判斷隧道的變形趨勢。

1.4 預(yù)測步驟

預(yù)測模型的預(yù)測步驟分述如下：

1)利用小波去噪和卡爾曼濾波對基坑變形的原始序列進(jìn)行去噪處理，并探討小波分解層數(shù)及小波函數(shù)對去噪效果的影響；同時，以信噪比和均方根誤差為基礎(chǔ)指標(biāo)，以其歸一化值累加作為綜合評價指標(biāo)，對比篩選出最優(yōu)去噪結(jié)果，將其作為趨勢項和誤差項的分離依據(jù)。

2)利用支持向量機(jī)和BP神經(jīng)網(wǎng)絡(luò)分別對趨勢項和誤差項進(jìn)行預(yù)測，并將兩者的結(jié)果進(jìn)行疊加，即得到基坑變形預(yù)測的綜合結(jié)果；將預(yù)測結(jié)果與文獻(xiàn)[12]的預(yù)測結(jié)果進(jìn)行對比，以驗證預(yù)測思路的有效性。

3)利用重標(biāo)度極差分析(R/S分析)對基坑變形的原始序列進(jìn)行變形趨勢判斷，并將判斷結(jié)果與預(yù)測結(jié)果進(jìn)行對比分析，相互驗證各自的準(zhǔn)確性。

預(yù)測步驟集合了多種方法的優(yōu)點，將其進(jìn)行結(jié)構(gòu)組合，以增強(qiáng)預(yù)測模型的穩(wěn)定性，同時利用R/S分析法對預(yù)測結(jié)果進(jìn)行驗證，進(jìn)一步核實了預(yù)測模型的有效性，為基坑變形預(yù)測提供了一種新的預(yù)測思路。另外，由于小波去噪[7]、卡爾曼濾波[13]及R/S分析[10]已在相關(guān)文獻(xiàn)中進(jìn)行了詳細(xì)的原理介紹，這里不再贅述。

2 實例分析

2.1 工程概況

某基坑[12]開挖深度較大，最大深度達(dá)32 m，其所處的地貌多以山間沖積洼地和剝蝕殘丘為主，地形起伏變化較小，相對高差為6.75 m。上覆土層主要是第四系覆蓋層，以人工填土、沖洪積層和殘積土層為主，下覆巖層為侵入巖，巖性為花崗巖?；拥叵滤皇艿匦渭敖涤甑挠绊戄^大，埋深在1.9～6.5 m，以孔隙潛水和基巖微承壓水為主。同時，基坑所處區(qū)域的地質(zhì)問題較多，如軟土震陷、砂土液化等，在施工擾動的作用下，易使土層結(jié)構(gòu)遭到破壞，進(jìn)而導(dǎo)致土體壓縮變形，影響基坑的穩(wěn)定，也易誘發(fā)鄰近建筑物的變形破壞。因此，基坑開挖過程中的施工監(jiān)測就顯得格外重要，根據(jù)監(jiān)測成果進(jìn)行預(yù)測分析，能有效實現(xiàn)信息反饋，以達(dá)到指導(dǎo)施工的目的。預(yù)測以Z70監(jiān)測點的沉降數(shù)據(jù)為例，共計110個周期，每周期為1周，其變形曲線如圖2所示。

圖2 基坑沉降變形曲線

由圖2可得出： 基坑的沉降變形具有明顯的階段性特征，表現(xiàn)為“緩—陡—緩—陡”的趨勢，最大變形量達(dá)504.2 mm； 其變化原因與施工階段相關(guān)，且累計變形量大，說明變形預(yù)測研究的必要性強(qiáng)。

再對基坑的變形速率進(jìn)行統(tǒng)計作圖，見圖3，得出基坑的沉降速率波動起伏較大，規(guī)律性較差，說明基坑變形受施工及環(huán)境因素的擾動較大，隨機(jī)性和波動性較強(qiáng)，這也從側(cè)面反映出基坑變形預(yù)測的困難，體現(xiàn)出對基坑變形綜合研究的價值?；映两底冃嗡俾实淖畲笾岛妥钚≈捣謩e為19.9 mm/周期和0.1 mm/周期，平均變形速率為4.58 mm/周期，說明基坑沉降變形波動較大。

圖3 基坑沉降速率變形曲線

2.2 去噪分析

采用小波去噪和卡爾曼濾波對基坑沉降的原始變形序列進(jìn)行去噪處理。由文獻(xiàn)[7]的研究，得到硬閾值方法和啟發(fā)式閾值標(biāo)準(zhǔn)的去噪效果較好，因此，僅對小波分解層數(shù)及小波函數(shù)的去噪效果進(jìn)行研究。其中，小波函數(shù)包含db小波系和sym小波系，分解層數(shù)分別為6、8、10、12、14層，小波系的去噪結(jié)果如表1和表2所示。

表1 db小波系去噪結(jié)果

在db小波系的去噪結(jié)果中，小波函數(shù)的階次和分解層數(shù)對去噪效果均有較大的影響。其中： 隨著分解層數(shù)的增加，去噪效果總體表現(xiàn)為先變優(yōu)后變差的特點，以12層的去噪效果為最優(yōu)，且去噪效果的穩(wěn)定性變化差異不大；在小波階次對去噪效果的影響方面，總體表現(xiàn)為隨著階次的增加，去噪效果相對更優(yōu)，以第6階次的去噪效果為最優(yōu)。同時，在db小波系的所有去噪結(jié)果中，以db7小波函數(shù)在10層小波分解時的去噪效果為最優(yōu)。

表2 sym小波系去噪結(jié)果

在sym小波系的去噪結(jié)果中，分解層數(shù)對去噪效果的影響規(guī)律與db小波系相同，只是最優(yōu)分解層數(shù)為10層，而最優(yōu)小波階次則為第6階次。同時，sym小波系的最優(yōu)去噪結(jié)果是sym6小波函數(shù)在10層分解時的去噪結(jié)果。

為進(jìn)一步對比2個小波系的去噪結(jié)果，再對2個小波系的去噪結(jié)果進(jìn)行對比分析，如圖4和圖5所示。

圖4 分解層數(shù)的影響對比

圖5 小波階次的影響對比

由圖4可知： 6層、8層和10層的去噪效果均是sym小波系相對更優(yōu)，且以6層的去噪效果差距最大；12層和14層的去噪效果以db小波系相對更優(yōu)。說明sym小波系更適用于分解層數(shù)較少的情況，而db小波系更適用于分解層數(shù)較多的情況。

由圖5可知，2個小波系均是隨小波階次的增加，去噪效果相對更優(yōu)，且隨著階次的增加，sym小波系的去噪效果要優(yōu)于db小波系的去噪效果，但差異均不大。說明小波階次越高，去噪效果越好，但影響不大。

對比2個小波系的穩(wěn)定性，得出分解層數(shù)去噪效果的穩(wěn)定性以db小波系相對略優(yōu)，而小波階次去噪效果的穩(wěn)定性則以sym小波系相對更優(yōu)。最后，綜合得到在小波去噪的結(jié)果中，以sym6小波函數(shù)在10層分解時的去噪效果為最優(yōu)。

卡爾曼濾波的去噪結(jié)果如表3所示。由表3可知，自適應(yīng)卡爾曼濾波的評價指標(biāo)值為1.582，而半?yún)?shù)卡爾曼濾波的評價指標(biāo)值為1.621，相對更優(yōu)。

表3 卡爾曼濾波的去噪結(jié)果

對比小波去噪和卡爾曼濾波去噪的結(jié)果，得出不同去噪方法所能達(dá)到的效果也具有差異，就最優(yōu)結(jié)果而言，以小波去噪的效果為最優(yōu)； 因此，以sym6小波函數(shù)在10層分解時的結(jié)果作為趨勢項和誤差項序列的分離依據(jù)。

2.3 變形預(yù)測

采用支持向量機(jī)對趨勢項序列進(jìn)行預(yù)測，預(yù)測結(jié)果如表4所示。根據(jù)趨勢項的預(yù)測結(jié)果，相對誤差值均小于2%，其中，最大相對誤差值為1.4%，最小相對誤差值為0.8%，且相對誤差平均值為1.07%，方差值為0.057 9，得出趨勢項的預(yù)測精度及穩(wěn)定性均較好，驗證了支持向量機(jī)在趨勢項預(yù)測中的有效性。

表4 趨勢項預(yù)測結(jié)果

采用BP神經(jīng)網(wǎng)絡(luò)對趨勢項序列進(jìn)行預(yù)測，采用試算法對最優(yōu)隱層節(jié)點進(jìn)行試算，結(jié)果如表5所示。通過試算得出隱層節(jié)點數(shù)對預(yù)測效果的影響是存在的，說明采用試算法確定最優(yōu)節(jié)點的方法是可行的。隨著隱層節(jié)點數(shù)的增加，誤差平方和表現(xiàn)為先減小后增大的規(guī)律，節(jié)點數(shù)為10時的誤差平方和最小，值為1.713 7 mm2，且后期增加的幅度小于前期減小的幅度，說明最優(yōu)節(jié)點數(shù)前的節(jié)點對預(yù)測精度的敏感性要大于最優(yōu)節(jié)點數(shù)后的節(jié)點對預(yù)測精度的敏感性。綜合試算結(jié)果，隱層節(jié)點確定為10，將該節(jié)點對應(yīng)的預(yù)測結(jié)果作為誤差項的預(yù)測結(jié)果。

表5 隱層節(jié)點試算結(jié)果

誤差項序列的預(yù)測結(jié)果如表6所示。在誤差項的預(yù)測過程中，預(yù)測誤差絕對值的最大值和最小值分別為0.7 mm和0.4 mm，且相對誤差絕對值在6%～12%，得出誤差項預(yù)測結(jié)果的特點總體表現(xiàn)為： 預(yù)測的絕對誤差值小，但相對誤差值大，與誤差序列的隨機(jī)性強(qiáng)相關(guān)。

表6 誤差項序列預(yù)測結(jié)果

對比趨勢項和誤差項的預(yù)測結(jié)果，在對應(yīng)監(jiān)測周期的相對誤差值，均是誤差項序列的值較大，且差異明顯。說明趨勢項序列的預(yù)測精度要明顯高于誤差項序列，這與誤差項序列含有較多的誤差因素有關(guān)，也從側(cè)面反映出通過去噪處理分離趨勢項和誤差項的方法是可行的，達(dá)到了預(yù)期的效果和目的。

通過對趨勢項和誤差項的預(yù)測，將兩者的結(jié)果進(jìn)行疊加，得到基坑變形的預(yù)測結(jié)果如表7所示。得出預(yù)測結(jié)果相對誤差的均值為1.03%，方差值為0.083，分離序列的預(yù)測精度及穩(wěn)定性均有所提高。同時，為對比預(yù)測結(jié)果的有效性，將文獻(xiàn)[12]中的預(yù)測結(jié)果一并進(jìn)行統(tǒng)計，得到其相對誤差均值為2.98%，最大相對誤差達(dá)-7.27%，方差值為7.368 4，相對本文預(yù)測結(jié)果，預(yù)測精度或穩(wěn)定性均較差，驗證了本文預(yù)測模型具有相對更好的優(yōu)越性。另外，對最后3個周期進(jìn)行遞推預(yù)測，得出基坑在后3個周期的變形仍呈增長趨勢，最大變形預(yù)測值已達(dá)516.01 mm。

表7 基坑變形預(yù)測結(jié)果

2.4 變形趨勢判斷

為進(jìn)一步判斷基坑的變形趨勢及驗證前文預(yù)測結(jié)果的可信度，采用R/S分析對基坑變形的變形序列和速率序列進(jìn)行重標(biāo)度極差分析，分析數(shù)據(jù)來源于基坑變形的實測數(shù)據(jù)，即前110周期的變形數(shù)據(jù)。經(jīng)過計算得到分析結(jié)果如表8所示。

表8 R/S分析結(jié)果統(tǒng)計

由表8可知： 變形序列的Hurst指數(shù)值為0.961，速率序列的Hurst指數(shù)值為0.585，均大于0.5。說明2個序列均具有長期的記憶性和持續(xù)性，是持續(xù)增強(qiáng)序列，基坑的變形將會繼續(xù)增加，這與基坑的變形預(yù)測結(jié)果一致，驗證了前期預(yù)測模型的準(zhǔn)確性； 但變形序列較速率序列的Hurst指數(shù)值要大，且速率序列的Hurst指數(shù)值接近0.5，說明變形序列的趨勢性較速率序列強(qiáng)，而速率序列相對具有更強(qiáng)的隨機(jī)性。2個序列的擬合度均較為接近1，誤差平方和也較小，說明重標(biāo)度極差分析結(jié)果的可信度較高，其中速率序列的擬合度略高于變形序列的擬合度。2個序列的相關(guān)性指標(biāo)CM均大于0，得出2個序列均為正相關(guān)，且變形序列的相關(guān)性程度明顯高于速率序列。

綜合對比變形預(yù)測及趨勢判斷的結(jié)果，得出基坑的變形具有持續(xù)增長的趨勢，需要加強(qiáng)現(xiàn)場監(jiān)測，并采取必要的控制變形措施。

3 結(jié)論與討論

1)通過小波去噪能很好地剔除誤差信息，分離趨勢項及誤差項，但不同的去噪?yún)?shù)對去噪效果的影響較大，通過試算法確定最優(yōu)去噪?yún)?shù)的可行性較好。在去噪分析過程中，sym小波系的去噪效果要優(yōu)于db小波系的去噪效果，而分解層數(shù)的去噪效果主要表現(xiàn)為隨分解層數(shù)的增加去噪效果先變優(yōu)后變差，存在最優(yōu)分解層數(shù)。

2)對比趨勢項和誤差項序列的預(yù)測結(jié)果，得出趨勢項的預(yù)測精度及穩(wěn)定性均明顯高于誤差項的預(yù)測精度及穩(wěn)定性，這是由于誤差序列含有更多的誤差因素所致，也從側(cè)面反映出去噪處理的有效性。

綜上所述，通過建立綜合預(yù)測模型進(jìn)行基坑的變形預(yù)測，并通過重標(biāo)度極差分析方法對變形趨勢進(jìn)行判斷，能進(jìn)一步驗證變形預(yù)測的可靠性，可為基坑的變形及穩(wěn)定性評價提供一種新的思路。

[1] 王寧偉, 朱豐, 韓旭, 等. 基于小波去噪的灰色馬爾柯夫基坑變形預(yù)測[J]. 水利與建筑工程學(xué)報, 2015(4): 215. WANG Ningwei, ZHU Feng, HAN Xu, et al. The deformation prediction of deep foundation pits based on a wavelet denoised grey Markov Model [J]. Journal of Water Resources and Architectural Engineering, 2015 (4): 215.

[2] 張世豪, 朱德華. 基于BP神經(jīng)網(wǎng)絡(luò)的深基坑變形預(yù)測[J]. 南京工程學(xué)院學(xué)報(自然科學(xué)版), 2016, 14(3): 17. ZHANG Shihao, ZHU Dehua. Deformation prediction of deep foundation pits based on BP neural network [J]. Journal of Nanjing Institute of Technology (Natural Science Edition), 2016, 14 (3): 17.

[3] 吳歡, 周春波, 秦昆. 基于改進(jìn)支持向量機(jī)的深基坑變形預(yù)測[J]. 地礦測繪, 2015, 31(4): 25. WU Huan, ZHOU Chunbo, QIN Kun. Deformation prediction of deep foundation pit based on improved support vector machine[J]. Surveying and Mapping of Geology and Mineral Resources, 2015, 31 (4): 25.

[4] 渠孟飛,謝強(qiáng),李朝陽,等. 基于支持向量機(jī)的膨脹土深基坑變形預(yù)測研究[J]. 施工技術(shù),2016,45(19): 30. QU Mengfei， XIE Qiang, LI Zhaoyang, et al. Prediction of slope displacement of deep foundation excavation based on support vector machine [J]. Construction Technology，2016, 45(19): 30.

[5] 林楠, 李偉東, 張文春, 等. 最小二乘支持向量機(jī)在深基坑變形預(yù)測中的應(yīng)用[J]. 遼寧工程技術(shù)大學(xué)學(xué)報, 2014(11): 1471. LIN Nan, LI Weidong, ZHANG Wenchun, et al. Application of least squares support vector machine to deformation prediction of deep foundation pit [J]. Journal of Liaoning Technical University, 2014 (11): 1471.

[6] 曾暉, 胡俊, 鮑俊安. 基于BP人工神經(jīng)網(wǎng)絡(luò)的基坑圍護(hù)結(jié)構(gòu)變形預(yù)測方法研究[J]. 鐵道建筑, 2011(1): 70. ZENG Hui, HU Jun, BAO Jun′an. Study of prediction method of deformation of foundation pit enclosure-structure based on BP (error back-propagation) artificial nerve network [J]. Railway Engineering, 2011 (1): 70.

[7] 楊哲峰, 羅林, 賈東彥, 等. 基于小波去噪的深基坑變形預(yù)測研究[J]. 人民長江, 2014(19): 41. YANG Zhefeng, LUO Lin, JIA Dongyan, et al. Prediction of deep foundation pit deformation based on wavelet denoising [J]. Yangtze River, 2014 (19): 41.

[8] 龐建成, 隆然, 馬嬌. 基于改進(jìn)R/S分析方法的深基坑變形預(yù)報預(yù)測研究[J]. 施工技術(shù), 2015, 44(13): 56. PANG Jiancheng, LONG Ran, MA Jiao. Study of deformation prediction of deep foundation pit based on improved R/S analysis [J]. Construction Technology, 2015, 44 (13): 56.

[9] 師旭超, 梁醒培. 深基坑變形預(yù)測的進(jìn)化支持向量機(jī)模型[J]. 武漢理工大學(xué)學(xué)報, 2010(1): 158. SHI Xuchao, LIANG Xingpei. An evolutionary support vector machine model for deformation prediction of deep foundation pit [J]. Journal of Wuhan University of Technology, 2010 (1): 158.

[10] 羅林, 左昌群, 趙連, 等. 基于BP神經(jīng)網(wǎng)絡(luò)和R/S分析的隧道仰坡沉降變形預(yù)報預(yù)測[J]. 施工技術(shù), 2014(11): 80. LUO Lin, ZUO Changqun, ZHAO Lian, et al. Prediction of settlement deformation of tunnel on the basis of BP neural network and R/S analysis [J]. Construction Technology, 2014 (11): 80.

[11] 龍浩, 高睿, 孔德新,等. 基于BP神經(jīng)網(wǎng)絡(luò)-馬爾科夫鏈模型的隧道圍巖位移預(yù)測[J]. 長江科學(xué)院院報, 2013, 30(3): 40. LONG Hao, GAO Rui, KONG Dexin, et al. Forecast of tunnel′s surrounding rock displacement by BP neural network and Markov Chain [J]. Journal of Yangtze River Scientific Research Institute, 2013, 30 (3): 40.

[12] 黃傳勝. 地鐵深基坑開挖變形預(yù)測方法及工程應(yīng)用研究[D]. 長沙： 中南大學(xué), 2011. HUANG Chuansheng. Study of deformation prediction method and engineering application of deep foundation pit excavation [D]. Changsha： Central South University, 2011.

[13] 王成, 何美琳, 覃婕,等. 半?yún)?shù)Kalman濾波模型在GPS變形數(shù)據(jù)處理中的應(yīng)用[J]. 施工技術(shù), 2015(增刊2)： 831. WANG Cheng, HE Meilin, QIN Jie, et al. Semiparametric Kalman filtering model in the application of GPS data processing of deformation [J]. Construction Technology, 2015 (S2): 831.

Study of Application of Trend Term Separation Prediction Model and Rescaled Range (R/S) Analysis to Deformation Prediction of Deep Foundation Pit

WANG Xueni1, HAN Guofeng2

(1.YanglingVocational&TechnicalCollege,Xianyang712100,Shaanxi,China;2.ShaanxiRailwayInstitute,Weinan714000,Shaanxi,China)

The deformation prediction accuracy of foundation pit is low nowadays. The denoising is carried out for deformation sequence of foundation pit by wavelet denoising and Calman filter, the trend term and error term are separated, and the trend term and error term is predicted by support vector machine (SVM) and BP neural network respectively. Meanwhile, the deformation trend of foundation pit is predicted by rescaled range (R/S) analysis so as to verify the feasibility of deformation prediction results. The case study shows that: 1) For wavelet denoising method, the denoising effect is superior; the average relative error and variance of the prediction results is 1.03% and 0.083 respectively; and the prediction accuracy is much higher. 2) The deformation sequence and deformation velocity sequence of foundation pit are prone to increasing, which coincide with the prediction results and verify the effectiveness of the prediction idea.

deep foundation pit; deformation prediction; denoising analysis; support vector machine (SVM); rescaled range (R/S) analysis; trend judgment

2017-01-04;

2017-04-24

王雪妮(1984—)，女，陜西富平人，2011年畢業(yè)于西安建筑科技大學(xué)，安全技術(shù)及工程專業(yè)，碩士，講師，主要從事建筑施工及安全管理工作。E-mail: 568756588@qq.com。

10.3973/j.issn.1672-741X.2017.08.012

U 452.1+1

A

1672-741X(2017)08-0990-07