周冬冬

(上?？睖y設計研究院有限公司，上海 200434)

基坑監(jiān)測數(shù)據(jù)缺失的動態(tài)數(shù)學模型-參數(shù)插值恢復方法

周冬冬

(上海勘測設計研究院有限公司，上海 200434)

本文提出了一種恢復因短暫監(jiān)測中斷而缺失的基坑側壁變形數(shù)據(jù)的方法，稱之為動態(tài)數(shù)學模型-參數(shù)插值法?；觽缺谧冃问芏喾N因素影響，但支護結構類型和約束情況是主要因素，因此，本文主要研究側壁無約束和頂部有約束兩種情況下的深層水平位移隨深度變化規(guī)律，提出了分別采用對數(shù)型和正態(tài)分布型動態(tài)曲線模型，并對模型參數(shù)進行擬合插值，最終求解獲得缺失數(shù)據(jù)期的監(jiān)測數(shù)據(jù)。與常用的平均法恢復數(shù)據(jù)對比后發(fā)現(xiàn)，測斜孔上部2/3深度部位的深層水平位移，動態(tài)數(shù)學模型-參數(shù)插值法計算結果更接近實際監(jiān)測值，但測斜孔下部1/3深度部位的深層水平位移，平均法計算結果更接近實際監(jiān)測值。

基坑；深層水平位移；動態(tài)數(shù)學模型；參數(shù)插值；數(shù)據(jù)恢復

1 概 述

基坑深層水平位移是支護結構與巖土體相互作用的綜合結果，能反映基坑邊坡的變形與破壞特征，是判斷基坑穩(wěn)定[1-4]與否的重要依據(jù)。通常采用布設測斜管、使用測斜儀進行監(jiān)測[5-9]，獲取連續(xù)準確的基坑深層水平位移數(shù)據(jù)。但在實際的監(jiān)測過程中，常常由于某些因素導致監(jiān)測工作短期中斷，如測斜儀故障、惡劣天氣和現(xiàn)場施工干擾等，使監(jiān)測數(shù)據(jù)存在短期缺失，從而影響對測點水平位移隨時間變化趨勢的分析，甚至對及時發(fā)現(xiàn)滑動異常、確定滑動面位置及分析滑動面發(fā)展趨勢等方面帶來不利影響。為此，有必要提出一種合理的方法恢復監(jiān)測中斷期的深層水平位移數(shù)據(jù)。

國內外對恢復缺失數(shù)據(jù)有很多方法[10-13]，主要分為兩類：基于統(tǒng)計的方法和基于數(shù)據(jù)挖掘[14-15]的方法。前者基于概率分布的統(tǒng)計假設，后者則需要巨大的數(shù)據(jù)集，數(shù)據(jù)處理較為復雜，花費時間較長。每種方法都有其特定的適用性，對不同的數(shù)據(jù)缺失情況，需要采用不同的缺失數(shù)據(jù)恢復方法。

本文提出了基坑側壁深層水平位移動態(tài)數(shù)學模型-參數(shù)插值方法，利用該方法恢復的缺失期水平位移更符合實際。

2 基坑側壁深層水平位移垂向變化規(guī)律

基坑側壁深層水平位移在垂向上的變化規(guī)律因支護方式不同而存在明顯差異，其變化規(guī)律大致可以分為兩類，見圖1、圖2。從測斜管頂部到底部，圖1的水平位移逐漸減小，最大水平位移出現(xiàn)在頂部；圖2的水平位移先增大后減小，最大水平位移出現(xiàn)在測斜管中上部。本文針對這兩種支護方式下的水平位移監(jiān)測曲線進行缺失數(shù)據(jù)的恢復。

圖1 頂部無約束的深層水平位移曲線

圖2 頂部有約束的深層水平位移曲線

3 頂部無約束基坑深層水平位移缺失數(shù)據(jù)的恢復計算

3.1 數(shù)學模型

將最大水平位移位于測斜管頂部時的深層水平位移數(shù)據(jù)監(jiān)測輸入到Excel表中，繪制該曲線并添加其趨勢線，根據(jù)趨勢線的相關系數(shù)大小，判斷將對數(shù)函數(shù)作為側壁深層水平位移的數(shù)學模型見下式：

s=alnh+c

(1)

考慮到圖1中的深層水平位移隨著深度的增加而變小，是一個單調函數(shù)。式(1)是一個關于深度h的單調函數(shù)，但為了使數(shù)學模型更符合實際，在式(1)中增加了同是單調函數(shù)的一次函數(shù)bh項，故

s=alnh+bch+c

(2)

式中s——水平位移，mm；

h——深度，m；

a、b、c——待定參數(shù)。

3.2 計算方法

a、b、c三個參數(shù)是變化的，每次監(jiān)測都有這三個參數(shù)。將缺失期前幾天的監(jiān)測數(shù)據(jù)以(h，s)形式輸入數(shù)學軟件Maple[16]中，并以式(2)的形式對每一天的數(shù)據(jù)進行擬合，得出每天的三個參數(shù)(ai,bi,ci)(i=1,2…，n)。

以第一天的參數(shù)為a1、b1、c1，第i天的參數(shù)為ai、bi、ci。將參數(shù)a1、a2、a3…以(1,a1)、(2,a2)、(3,a3)…的形式輸入Excel中，利用Excel畫趨勢線并能給出公式的功能，可以得出參數(shù)a關于i的函數(shù)：

a=f(i)

(3)

利用式(3)就能計算出缺失數(shù)據(jù)期的參數(shù)aq，同理也能得出缺失期的參數(shù)bq、cq，將這三個參數(shù)代入式(2)中，就能得出缺失期的深層水平位移關于深度h的函數(shù)：

s=aqlnh+bqh+cq

(4)

根據(jù)式(4)，即可計算出監(jiān)測數(shù)據(jù)缺失期在不同深度上的水平位移。

3.3 工程實例

上海某1號基坑開挖深度為3.5～6.5m，基坑側壁使用高壓旋噴樁進行支護，屬于頂部無約束基坑?；又ёo深度范圍內以淤泥質粉質黏土為主，在墻后土體中埋設測斜管cx1，埋設深度為14m，該管從2003年6月12日開始連續(xù)幾天水平位移變化量均比較大，最大水平位移均出現(xiàn)在測斜管的頂部。假設6月15日由于測斜儀故障導致未對測斜管進行監(jiān)測，在6月16日恢復了測量，但在6月15日對基坑側壁已經采用了竹釘加固處理。實際的監(jiān)測數(shù)據(jù)見表1。

表1 cx1深層水平位移監(jiān)測數(shù)據(jù) 單位：mm

從表1的監(jiān)測數(shù)據(jù)可以發(fā)現(xiàn)，6月16日測斜管不是向坑內而是向坑外變形，如果不知道6月15日的監(jiān)測數(shù)據(jù)，則很容易誤以為測斜管一直在向坑內變形，這就導致很難判斷竹釘是否起到了阻止基坑過大變形的作用，所以獲得較為準確的6月15日的監(jiān)測數(shù)據(jù)就顯得很重要。

將6月11—14日的監(jiān)測數(shù)據(jù)與深度分別輸入數(shù)學軟件Maple中，以式(1)進行擬合，并將每天的實測值與利用擬合公式計算出的值進行對比，以最大絕對值偏差(實測值與擬合值之差絕對值的最大值)來反映擬合的效果。擬合計算所得的參數(shù)a、b、c值見表2。

表2 cx1 6月11—14日監(jiān)測數(shù)據(jù)擬合參數(shù)值

將所得的a、b、c的值分別輸入Excel中，利用Excel的趨勢線功能，分別求出a、b、c與i的函數(shù)關系：

令i=5，即可求出缺失的6月15日深層水平位移參數(shù)aq=-7.4799、bq=0.0031、cq=20.3851,再利用式(4)，即可求出6月15日測斜管在不同深度上的水平位移，將其與實測值以及利用平均值法(取缺失期前后兩天監(jiān)測值的平均值)所得的平均值進行對比，見表3。

表3 擬合計算值與實測值及平均值對比 單位：mm

從表3可以發(fā)現(xiàn)，在測斜管上部，利用擬合公式所得的計算值更加接近于實測值，而在測斜管下部，平均值則更加接近實測值，這是因為測斜管下部位移變化相對較小，有時甚至幾乎沒有變化。盡管測斜管下半部分所得的擬合值與實測值相比誤差較大，但是對于底部水平位移很小的測斜孔，其上半部分水平位移的變化情況無疑更為重要。因此，為了更加準確地恢復缺失數(shù)據(jù)，取擬合計算值作為測斜管上部2/3部分的缺失值，而測斜管下部1/3部分則取平均值作為缺失值，最終所得結果見表3。

4 頂部有約束基坑深層水平位移缺失數(shù)據(jù)的恢復計算

4.1 數(shù)學模型

最大水平位移在頂部以下時深層水平位移的變化規(guī)律見圖2，根據(jù)其中間大兩頭小的特征，考慮其數(shù)學模型為正態(tài)分布的概率密度函數(shù)模型：

(8)

對式(8)進行化簡，將其合并為

(9)

事實上，圖2中的深層水平位移變化規(guī)律并非一個完整的正態(tài)分布模型，其對稱軸位于測斜管中部以上的部分。因此，在計算擬合值的過程中，對深度在0.5m以上部分的深層水平位移擬合值予以舍棄。

4.2 求解方法

與3.2的計算方法一致。

4.3 工程實例

同樣是上海某基坑，開挖深度為6.8m，基坑側壁采用水泥攪拌樁形成的連續(xù)墻進行支護，攪拌樁頂部有圈梁壓頂，墻后土體中埋設測斜管cx2，埋設深度為20m，其監(jiān)測數(shù)據(jù)見表4，假設7月24日由于天氣原因沒有測量，在7月25日，對基坑壁進行了土釘支護處理，且進行了底板澆筑。

表4 cx2深層水平位移監(jiān)測數(shù)據(jù) 單位：mm

將7月20—23日的深度值與位移變形值以(h,s)的形式輸入Maple中，并以式(7)的形式進行正態(tài)擬合，Maple將以式(8)的形式顯示擬合的結果。將實測值與利用擬合公式所得的計算值對比，并以最大絕對值偏差來評估擬合的效果。擬合計算所得參數(shù)值見表5。

表5 cx2測斜孔7月20日至7月23日監(jiān)測數(shù)據(jù)擬合參數(shù)

將所得的a、b、c的值分別輸入Excel中，分別求出a、b、c與天數(shù)i的函數(shù)關系：

a=0.042469i+0.432774

(10)

b=0.000253i-0.018567

(11)

c=-0.057613i+2.532177

(12)

令i=5，即可求出缺失期7月24日的深層水平位移參數(shù)aq=0.645119、bq=-0.017304、cq=2.244110，再將其代入式(9)，即可求出7月24日測斜管在不同深度上的水平位移，將其與實測值以及利用平均值法所得的平均值進行對比，見表6。

表6 擬合計算值與實測值及平均值對比 單位：mm

續(xù)表

從表6可以發(fā)現(xiàn)，在測斜管的中上部，擬合值與實測值更加接近，而在下部，則是平均值更加符合實際情況。因此，在測斜管上部2/3部分，取擬合值作為最終值；在下部1/3部位，以平均值作為最終值。

5 結 語

根據(jù)基坑不同的支護結構類型和變形規(guī)律分別建立基坑側壁深層水平位移關于深度的動態(tài)數(shù)學模型，通過模型參數(shù)的擬合插值，計算基坑監(jiān)測短暫中斷期缺失的監(jiān)測數(shù)據(jù)，這種方法是可行的。本文針對基坑側壁深層水平位移的兩種典型變形規(guī)律提出了兩種不同的動態(tài)數(shù)學模型-參數(shù)插值方法，擬合計算的結果與常用的平均法計算結果對比顯示：測斜孔上部2/3深度部位的深層水平位移，本擬合計算方法所得到的結果更接近實際監(jiān)測值，但測斜孔下部1/3深度部位的深層水平位移，平均法計算結果更接近實際監(jiān)測值。

[1] 石林珂,賀為民,孫懿斐,等.深基坑穩(wěn)定分析中的等弧長條分法及可視化軟件設計[J].巖石力學與工程學報,2002,21(10):1568-1572.

[2] 劉鑫,洪寶寧,韓尚宇.地鐵車站基坑穩(wěn)定風險值計算[J].中南大學學報(自然科學版),2012,43(6):2331-2337.

[3] 張玉成,楊光華,吳舒界,等.土釘支護結構變形與穩(wěn)定性關系探討[J].巖土力學,2014,35(1):238-247.

[4] S.L.Huang.Stability analysis of mine pit slopes in interior Alaska[J].International Journal of Rock Mechanics and Mining Sciences,2004,41:816-820.

[5] 高俊和,趙維炳,施建勇,等.基坑測斜技術研究[J].河海大學學報,1998,26(3):120-122.

[6] 蔡干序.基坑監(jiān)測工程測斜技術的探討[J].建筑科學,2009,25(11):99-102.

[7] 劉金龍,朱建群,王吉利,等.測斜儀在路基水平位移監(jiān)測中的若干問題探討[J].湖南科技大學學報(自然科學版),2011(9):159.

[8] Lucia Simeoni,Luigi Mongiovi.Inclinometer Monitoring of the Castelrotto Landslide in Italy[J].Journal of Geotechnical and Geoenvironmental Engineering,2007,133(6):653-666.

[9] Timothy D.Stark,Hangseok Choi.Slope inclinometer for landslides[J].2008,5(3):339-350.

[10] 潘立強,李建中,駱吉洲.傳感器網絡中一種基于時空相關性的缺失值估計算法[J].計算機學報,2010,33(1):1-11.

[11] 彭紅毅,朱思銘,蔣春福.數(shù)據(jù)挖掘中基于ICA的缺失數(shù)據(jù)值的估計[J].計算機科學,2005,32(12):203-205.

[12] 武森,馮小東,單志廣.基于不完備數(shù)據(jù)聚類的缺失數(shù)據(jù)填補方法[J].計算機學報,2012,35(8):1726-1738.

[13] Tseng S,Wang K,Lee C.A pre-processing method to deal with missing values by integrating clustering and regression technique[J].Applied Artificial Intelligence,2003,17(5-6):535-544.

[14] 李德仁,王樹良,李德毅,等.論空間數(shù)據(jù)挖掘和知識發(fā)現(xiàn)的理論與方法[J].武漢大學學報(信息科學版),2002,27(3):221-233.

[15] Vladimir Nedic,Slobodan Cvetanovic,DanijelaDespotovic,et al.Data mining with various optimization methods[J].Expert Systems with Applications,2014,41(8):3993-3999.

[16] 晏水平,黃樹紅.生成多體系統(tǒng)運動方程的MAPLE程序設計[J].華中科技大學(自然科學版),2002,30(5):62-64.

Dynamic Mathematical Model-parameter Interpolation Recovery Method for Monitoring Data Missing of Foundation Pits

ZHOU Dongdong

(SHANGHAIINVESTIGATION,DESIGN&RESEARCHINSITUTE,Shanghai200434,China)

In the paper, a method of restoring missing foundation pit side wall deformation data due to short monitoring interruption is proposed. It is called dynamic mathematical model-parameter interpolation method. The foundation pit side wall deformation is affected by many factors. However, supporting structure type and constraint condition are the main factors. Therefore, the deep horizontal displacement change law with depth under the two conditions of no constraint on the side wall and constraint at the top is mainly studied in the paper. It is proposed that logarithmic and normal distribution dynamic curve models are adopted respectively for fitting interpolation on model parameters. Finally, the monitoring data during the data missing period is solved and obtained. The method is compared with the average method through data recovery. It is discovered that the deep horizontal displacement in the deep part 1/3 above the inclined hole is measured. The calculation result of the dynamic mathematical model-parameter interpolation method is closer to the measured value. However—when the deep horizontal displacement in the deep part 1/3 above the inclined hole is measured. The calculation result of the average method is closer to the measured value.

foundation pit; deep horizontal displacement; dynamic mathematical model; parameter interpolation; data recovery

10.16617/j.cnki.11-5543/TK.2017.07.015

TV551.4

A

1673-8241(2017)07- 0052- 07