劉 賀，張弘強，劉 斌

1.吉林大學交通學院，長春 130022

2.吉林農業(yè)工程職業(yè)技術學院，吉林 四平 136001

3.遼寧城建設計院有限公司，遼寧 撫順 113008

0 引言

隨著地下空間的不斷開發(fā)，深基坑工程呈現(xiàn)出規(guī)模不斷擴大、深度不斷加深的發(fā)展趨勢，因此，確?；邮┕ぐ踩蔀樯罨庸こ痰氖滓蝿?。深基坑工程的變形不僅涉及土的力學特性，還受到土與結構相互作用、地下水位和環(huán)境溫度等因素的影響，開挖過程其變形具有明顯的時間和空間效應。所以，深基坑安全施工是一個動態(tài)、復雜的系統(tǒng)工程。深基坑變形預測是確保深基坑施工安全的一項重要措施，在施工方案調整和施工風險規(guī)避等方面具有重要的實用價值［1］。

針對深基坑變形預測方法，國內外專家學者展開了相應的研究，將數(shù)值模擬方法應用到基坑變形預測中［2－4］。該方法的基本思想是建立基坑有限元分析模型，在有限元分析模型中輸入計算參數(shù)，對基坑的變形進行計算。由于深基坑工程是一個動態(tài)非線性過程，加上內部和外部參數(shù)的不確定性，導致有限元模型的輸入?yún)?shù)難以準確計算。

由于數(shù)值模擬方法難以對基坑變形進行準確預測，系統(tǒng)分析方法在基坑變形預測中得以廣泛應用。常用的基坑變形預測系統(tǒng)分析方法主要包括時間序列分析［5］、灰色系統(tǒng)理論［6］、模擬退火算法［7］、遺傳算法［8－9］以及人工神經(jīng)網(wǎng)絡理論［10－11］等。 其中人工神經(jīng)網(wǎng)絡理論具有較強的自適應和自學習能力，對于高度復雜的非線性問題具有良好的模擬能力，能較好地模擬基坑變形的發(fā)展規(guī)律。針對復雜多因素影響系統(tǒng)，人工神經(jīng)網(wǎng)絡能夠克服其他系統(tǒng)分析方法的不足，建立準確的預測模型，實現(xiàn)基坑變形預測。但神經(jīng)網(wǎng)絡算法存在對初始權值和閥值較為敏感、容易陷入局部極小值和收斂速度慢等缺點。利用粒子群優(yōu)化算法對其進行優(yōu)化，能夠改善神經(jīng)網(wǎng)絡性能，提高深基坑變形預測的準確性［12］。

粒子群優(yōu)化算法（particle swarm optimization，PSO）是由J.Kennedy和R.C.Eberhart［13］根據(jù)鳥群覓食過程中的遷徙和聚集現(xiàn)象提出的一種基于群體智能原理的非線性優(yōu)化算法。它采用種群的全局搜索策略和簡單位移－速度模式在解空間搜索最優(yōu)解；通過粒子的不斷更新，完成解空間的動態(tài)全面搜索［14］；采用粒子群優(yōu)化算法改進神經(jīng)網(wǎng)絡，通過粒子迭代求出全局優(yōu)化初始權值和閾值，應用于神經(jīng)網(wǎng)絡，可以提高神經(jīng)網(wǎng)絡的性能［15］。

筆者采用粒子群算法對神經(jīng)網(wǎng)絡模型中的閥值和初始權值進行優(yōu)化，建立了基坑變形預測模型，提出了深基坑變形預測方法。將前期實測變形數(shù)據(jù)作為神經(jīng)網(wǎng)絡的輸入?yún)?shù)，對神經(jīng)網(wǎng)絡進行訓練，訓練的神經(jīng)網(wǎng)絡能夠對基坑后續(xù)變形進行準確預測。

1 基坑變形預測方法的基本理論及原理

1.1 神經(jīng)網(wǎng)絡模型的設計

三層BP網(wǎng)絡能夠模擬復雜非線性的輸入輸出關系［16］。筆者采用三層BP神經(jīng)網(wǎng)絡作為基坑變形預測模型進行基坑變形預測。采用的BP神經(jīng)網(wǎng)絡結構由輸入層、隱含層和輸出層組成。

1.1.1 輸入層設計

基坑的變形受開挖形式、開挖深度、地下水位以及土壓力等眾多因素的影響。隨著基坑的開挖，變形是一個連續(xù)動態(tài)過程，具有極強的時間效應。把影響因素直接作為神經(jīng)網(wǎng)絡輸入變量，難以對基坑變形進行準確預測。綜合考慮各因素對變形的影響，把前期實測變形值作為神經(jīng)網(wǎng)絡模型輸入變量，預測基坑后續(xù)變形值。

1.1.2 隱含層設計

確定隱含層的節(jié)點個數(shù)是隱含層設計的主要任務。隱含層節(jié)點個數(shù)過多，會降低神經(jīng)網(wǎng)絡算法的泛化能力，增加樣本訓練時間；節(jié)點個數(shù)過少，則會降低神經(jīng)網(wǎng)絡算法獲取訓練樣本有效信息的能力。筆者采用公式（1）對隱含層節(jié)點數(shù)進行估算，確定隱含層由10個節(jié)點組成，即

式中：a為隱含層節(jié)點數(shù)；b為輸入層節(jié)點數(shù)；l為輸出層節(jié)點數(shù)。

1.1.3 輸出層設計

Pureline函數(shù)是一個線性函數(shù)，函數(shù)形式如圖1所示。筆者采用Pureline函數(shù)作為輸出層的傳遞函數(shù)，能夠使輸出層對輸出值具有較好的敏感性。

圖1 Pureline函數(shù)曲線Fig.1 Pureline function curve

1.2 粒子群優(yōu)化算法

粒子群優(yōu)化算法［14］的基本思想為M個無質量、無體積的粒子構成一個粒子群，在D維搜索空間內以一定速度飛行，尋找全局最優(yōu)解，飛行位置和速度分別由位置向量和速度向量確定。粒子在每一時刻跟蹤2個極值運動，一個是粒子群的全局最優(yōu)位置，另一個是粒子自身最優(yōu)位置。

粒子群優(yōu)化算法通過迭代搜尋每個粒子的當前最優(yōu)解，采用適應度函數(shù)評價解的優(yōu)劣程度，比較粒子的適應度函數(shù)值，確定全局最優(yōu)解。

1.3 粒子群優(yōu)化神經(jīng)網(wǎng)絡算法

粒子群優(yōu)化神經(jīng)網(wǎng)絡算法中，粒子的個體極值點和全局最優(yōu)極值點以及神經(jīng)網(wǎng)絡模型中權值和閾值的優(yōu)化終止條件都是基于粒子的適應度來確定的。粒子的適應度函數(shù)為

式中：Yi，j為第i個樣本的第j個理想輸出值；yi，j為第i個樣本的第j個實際輸出值；n為樣本數(shù)；m＝1，2，…，M，M為粒子種群規(guī)模，即粒子數(shù)。

將粒子群優(yōu)化算法的誤差要求作為神經(jīng)網(wǎng)絡算法閾值和權值的優(yōu)化終止條件，誤差采用下式進行計算：

粒子群優(yōu)化神經(jīng)網(wǎng)絡算法的流程如下：

1）確定神經(jīng)網(wǎng)絡結構。

2）初始化神經(jīng)網(wǎng)絡權值和閾值。

3）建立粒子群與權值和閾值的關系，對粒子群進行初始化。

4）將基坑前期實測變形值作為輸入值，進行神經(jīng)網(wǎng)絡訓練，通過式（4）、（5）計算適應度，對適應度進行比較和評價，按式（2）、（3）更新每個粒子的位置和速度。

5）通過式（6）計算粒子群優(yōu)化算法的誤差。如果誤差滿足要求，則終止優(yōu)化，輸出權值和閾值優(yōu)化結果；若誤差不滿足要求，則繼續(xù)優(yōu)化，重復步驟4），直到誤差滿足要求。

6）在神經(jīng)網(wǎng)絡模型中，采用優(yōu)化后的權值和閾值進行樣本訓練，形成基坑變形預測模型。

7）根據(jù)已有的實測數(shù)據(jù)，利用神經(jīng)網(wǎng)絡模型預測測點即將發(fā)生的變形。

2 實際工程應用

2.1 工程概況

選取長春市火車站綜合交通換乘中心北廣場為實體工程?；訓|西方向長550m，南北方向寬150 m；深度為7.45～14.00m，采用明挖施工形式。由于基坑東側較淺，采用兩級放坡開挖，網(wǎng)噴混凝土護面對土質邊坡進行防護。其余3面開挖深度較深，采用1級放坡開挖，坡率為1∶1。同時采用鉆孔灌注樁和網(wǎng)噴混凝土護面進行防護，圍護樁徑1.0m，中心間距1.2m?？辈旖Y果表明，擬建場地地層沉積具有一定的規(guī)律性，場地地層主要由三部分組成：地表一般分布有道路結構層和人工堆積雜填土層，厚度為1.00～4.20m；中間部分為第四系沖擊黏性土和沖洪積砂土，由粉質黏土、黏土和粗砂組成，其中粉質黏土的厚度為5.1～21.8m，黏土厚度為5.20～10.80m，粗砂厚度為0.70～4.00m；最下層為白堊紀泥巖。

對整個基坑開挖過程中的基坑變形進行了監(jiān)測，在基坑周圍共布設34個水平位移測點和32個地面沉降測點，測點布置如圖2所示（其中H表示水平位移測點，V表示地面沉降測點）。由于在基坑監(jiān)測過程中實際的測點數(shù)較多，為了驗證本文方法的有效性，只對8號水平位移測點H8和9號地面沉降測點V9的數(shù)據(jù)進行預測。這兩個位移監(jiān)測點位于基坑長度的中點，即深度最深的位置。

圖2 測點布置示意圖Fig.2 Layout of deformation measuring point

2.2 基坑變形預測模型參數(shù)選取

利用前15次的監(jiān)測數(shù)據(jù)來構造訓練樣本。指定測點的某次監(jiān)測數(shù)據(jù)為輸出，測點之前的6次監(jiān)測數(shù)據(jù)為輸入，則該測點的6次輸入數(shù)據(jù)和1次輸出數(shù)據(jù)構成了1組訓練樣本。利用前15次的監(jiān)測數(shù)據(jù)共構建了10組訓練樣本。

BP神經(jīng)網(wǎng)絡模型具有6個輸入變量和1個輸出變量，隱含層具有10個節(jié)點。根據(jù)連接權值和閾值的個數(shù)，確定粒子群的規(guī)模為20。

2.3 預測結果

利用粒子群優(yōu)化神經(jīng)網(wǎng)絡預測方法對8號水平位移測點和9號地面沉降點的變形數(shù)據(jù)進行了預測，并利用均方根誤差（RMSE）和平均百分比誤差（MAPE）對預測結果進行評價：

式中：ti為第i個實測值；oi為第i個預測值；n為數(shù)據(jù)數(shù)。

8號測點的水平位移預測結果如圖3所示，9號地面沉降點的預測結果如圖4所示。

圖3 8號水平位移測點的預測結果Fig.3 Prediction result of No.8horizontal displacement measuring point

圖4 9號地面沉降點的預測結果Fig.4 Prediction result of No.9land subsidence point

由預測結果可知：1）8號水平位移測點預測結果的RMSE＝3.78%，MAPE＝5.48%；9號地面沉降點預測結果的RMSE＝5.62%，MAPE＝3.23%。一般來說，在基坑監(jiān)測過程中，變形監(jiān)測的誤差小于8%即可，本文預測誤差均小于該誤差限值，這說明本文提出的方法可以較為準確地預測出深基坑開挖過程中的變形。2）8號水平位移測點和9號地面沉降點的變化趨勢不一致，這是由于2個測點的觀測時間不同導致的，根據(jù)施工進度，9號地面沉降點的布置及觀測時間要比8號水平位移測點早3個月左右。

3 結論與建議

1）針對神經(jīng)網(wǎng)絡算法在變形預測中存在的問題，采用粒子群優(yōu)化算法對神經(jīng)網(wǎng)絡模型參數(shù)進行優(yōu)化，從而形成了基于粒子群優(yōu)化神經(jīng)網(wǎng)絡算法的深基坑變形預測方法。該方法僅利用前期的實測變形數(shù)據(jù)來預測同一測點的后續(xù)變形值。實際應用表明，變形預測誤差均小于6%，滿足基坑變形監(jiān)測的要求，說明本方法在深基坑開挖的變形預測中可得到有效應用。

2）基坑變形是評價基坑開挖過程中安全性的一個重要參數(shù)，本方法可以預測出一個較為準確的理論參考值，在基坑監(jiān)測中可以把變形實測值與該理論參考值進行比較來評價基坑的安全性。本研究為同類工程的變形監(jiān)測和預測提供了一種不同的思路。

在部分監(jiān)測數(shù)據(jù)失真或誤差過大的情況下，該方法如何給出一個準確的預測值，還需要進行進一步的研究，以便提高該方法的實用性。

（References）：

［1］王寧，黃銘.開挖作用下的深基坑變形神經(jīng)網(wǎng)絡監(jiān)測模型［J］.上海交通大學學報，2009，43（6）：990－994.Wang Ning，Huang Ming.Deformation Monitoring Artificial Neural Network Model of Deep Foundation Pit Considering the Excavation Effect［J］.Journal of Shanghai Jiaotong University，2009，43（6）：990－994.

［2］陳興賢，駱祖江，安曉宇，等.深基坑降水三維變參數(shù)非穩(wěn)定滲流與地面沉降耦合模型［J］.吉林大學學報：地球科學版，2013，43（5）：1572－1578.Chen Xingxian，Luo Zujiang，An Xiaoyu，et al.Coupling Model of Groundwater Three Dimensional Variable－Parametric Non－Steady Seepage and Land－Subsidence［J］.Journal of Jilin University：Earth Science Edition，2013，43（5）：1572－1578.

［3］劉志凱，鄭毅，左廣州.基坑開挖與土釘支護的數(shù)值模擬［J］.吉林大學學報：地球科學版，2004，34（增刊1）：103－106.Liu Zhikai，Zheng Yi，Zuo Guangzhou.Numerical Simulation of the Foundation Digging and the Soil Nail Protection［J］.Journal of Jilin University：Earth Science Edition，2004，34（Sup.1）：103－106.

［4］曹國金，姜弘道.信息化施工技術在地下工程中的應用［J］.巖土力學，2002，23（6）：795－799.Cao Guojin， Jiang Hongdao.Application of Information Technology in Underground Engineering［J］.Rock and Soil Mechanics，2002，23（6）：795－799.

［5］袁金榮，趙福勇.基坑變形預測的時間序列分析［J］.土木工程學報，2001，34（6）：55－59.Yuan Jinrong，Zhao Fuyong.Time Series Analysis of Foundation Pit Deformation Prediction［J］.China Civil Engineering Journal，2001，34（6）：55－59.

［6］陳曉斌，張家生，安關峰.GM（1，1）與 GM（2，1）模型在基坑工程預測中的應用［J］.巖土工程學報，2006，28（增刊）：1401－1405.Chen Xiaobin， Zhang Jiasheng， An Guanfeng.Application of GM（1，1）Model and GM（2，1）Model in Prediction of Foundation Pit［J］.Chinese Journal of Geotechnical Engineering，2006，28（Sup.）：1401－1405.

［7］張發(fā)明，龐正江，陳祖爆，等.糯扎渡電站溢洪道邊坡穩(wěn)定分析及系統(tǒng)錨固方法［J］.巖土力學，2005，（增刊1）：255－259.Zhang Faming，Pang Zhengjiang，Chen Zubao，et al.Slope Stability Analysis of Nuozhadu Power Station Spillway and System Anchoring Method［J］.Rock and Soil Mechanics，2005，（Sup.1）：255－259.

［8］周海清，陳正漢.面向對象的深度搜索遺傳算法及其工程應用：I：算法與程序［J］.巖石力學與工程學報，2005，24（11）：1996－2002.Zhou Haiqing，Chen Zhenghan.Object－Oriented Depth Search Genetic Algorithm and Its Engineering Application：I：Algorithm and Program［J］.Chinese Journal of Rock Mechanics and Engineering，2005，24（11）：1996－2002.

［9］周海清，陳正漢.面向對象的深度搜索遺傳算法及其工程應用：II：算法與程序［J］.巖石力學與工程學報，2005，24（12）：2194－2206.Zhou Haiqing，Chen Zhenghan.Object－Oriented Depth Search Genetic Algorithm and Its Engineering Application：II：Algorithm and Program［J］.Chinese Journal of Rock Mechanics and Engineering，2005，24（12）：2194－2206.

［10］Mayoraz F，Vulliet L.Neural Networks for Slope Movement Prediction［J］.The International Journal of Geotechnics，2002（2）：153－174.

［11］賀志勇，鄭偉.基于BP神經(jīng)網(wǎng)絡的深基坑變形預測［J］.華南理工大學學報：自然科學版，2008，36（10）：92－96.He Zhiyong，Zheng Wei.Deformation Prediction of Deep Foundation Pit Based on BP Neural Network［J］.Journal of South China University of Technology：Natural Science Edition，2008，36（10）：92－96.

［12］周蕾.粒子群算法的改進及其在人工神經(jīng)網(wǎng)絡中的應用［D］.西安：西安電子科技大學，2010.Zhou Lei.Improved Particle Swarm Optimization Algorithm and the Application of Artificial Neural Network［D］.Xi’an：Xi’an University of Electronic Science and Technology，2010.

［13］Kennedy J，Eberhart R C.Particle Swarm Optimization［J］.Proceedings of the IEEE International Joint Conference on Neural Networks，1995，4：1942－1948.

［14］紀震，廖惠連，吳青華.粒子群算法及應用［M］.北京：科學出版社，2009.Ji Zhen，Liao Huilian，Wu Qinghua.Particle Swarm Optimization and Application［M］.Beijing：Science Press，2009.

［15］韓月嬌.粒子群算法的改進及其在BP神經(jīng)網(wǎng)絡中的應用［D］.南昌：南昌航空大學，2012.Han Yuejiao.Improvement on Particle Swarm Algorithm and Its Application in BP Neural Network［D］.Nanchang：Nanchang Aeronautical University，2012.

［16］Mehrjoo M，Khaji N，Moharrami H，et al.Damage Detection of Truss Bridge Joints Using Artificial Neural Networks［J］.Expert Systems with Applications，2008，35：1122－1131.