考慮數(shù)據(jù)異常及新舊程度影響有界性的地基沉降預測方法

曹文貴 譚濤

摘 ? 要：在現(xiàn)有沉降組合模型預測方法研究基礎上，通過考慮實測沉降數(shù)據(jù)新舊程度及其影響的有界性特點，提出了度量實測沉降數(shù)據(jù)新舊程度的新鮮度函數(shù)分析模型，它能同時反映實測點與預測點之間的時間距離和實測誤差以及實測數(shù)據(jù)新舊程度影響有界性對沉降預測的影響. 基于統(tǒng)計學理論提出了排除異常實測數(shù)據(jù)引起沉降預測不合理的實測數(shù)據(jù)樣本處理方法. 利用上述模型與方法，提出了可反映實測沉降數(shù)據(jù)異常和新舊程度及其影響程度有界性的地基或路基沉降組合預測新方法. 通過工程實例計算以及本文與現(xiàn)有同類方法預測與實測結果的比較分析，表明了本文地基或路基預測模型與方法的合理性與優(yōu)越性.

關鍵詞：地基或路基;沉降預測;組合模型;新鮮度函數(shù);實測數(shù)據(jù)異常;實測數(shù)據(jù)新舊程度影響有界性

中圖分類號：TU470 ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標志碼：A

A Prediction Method of Foundation Settlement Considering Anomaly

and Newness-oldness Degree Influence Boundedness of Measured Data

CAO Wengui？覮，TAN Tao

（Institute of Geotechnical Engineering，Hunan University，Changsha 410082，China）

Abstract：Based on the research of prediction methods for foundation or subgrade settlement by existing combination model，firstly，an analytical model of freshness function is proposed to measure the newness-oldness degree of measured settlement data，considering the features of newness-oldness degree and its influence boundedness of measured settlement data. The proposed method can also reflect the effect of time distance between the measured and predicted point，measured error and boundedness of newness-oldness as well as its influence degree of the measured data on the settlement prediction. Secondly，based on statistical theory，a disposal method of measured data sample is developed which can exclude the unreasonable settlement prediction caused by some abnormal measured data. Then，using the above-mentioned model and method，a new method for combined prediction of foundation or subgrade settlement is put forward，which can reflect the influence of abnormality of measured settlement data，their oldness and newness degree and boundedness of its influence degree on settlement prediction of foundation or subgrade. Finally，it shows the rationality and superiority of the proposed model and method by calculation of actual engineering examples and comparison and analysis among the prediction results by the proposed and the existing similar methods as well as the measured results.

Key words：foundation or subgrade;settlement prediction;combination model;freshness function;measured data anomaly; newness-oldness degree influence boundedness of measure data

眾所周知，實測沉降數(shù)據(jù)是地基或路基沉降預測的基礎，地基沉降預測的合理性除了取決于預測模型或方法與實測數(shù)據(jù)量的大小之外，還取決于實測數(shù)據(jù)的可靠程度及其新舊程度即實測數(shù)據(jù)點與預測點之間的時間距離，而且，實測數(shù)據(jù)新舊程度的影響并非無止境即存在有界性，也就是說，雖然實測與預測點之間的時間距離越小時，該實測點的數(shù)據(jù)應占有更高的權重，但是，這種影響并非隨時間距離減小而無限增大，或者說，隨時間距離增大而無限減小，因此，在進行地基或路基沉降合理預測時須對實測數(shù)據(jù)的可靠性與新舊程度影響的有界性予以充分重視，這正是本文研究的出發(fā)點.

目前，地基或路基沉降預測方法研究主要存在兩條基本途徑.其一就是單項模型預測方法，例如，雙曲線模型、指數(shù)曲線模型、Asaoka法、Logistic模型、Weibull模型、MMF模型等[1-7]，其基本思路是依據(jù)地基或路基沉降發(fā)展變化規(guī)律，利用已有實測沉降數(shù)據(jù)或曲線，采用單一函數(shù)模型進行擬合以獲得沉降與時間之間經(jīng)驗關系，進而實現(xiàn)對地基或路基沉降的預測.該方法的突出優(yōu)越性表現(xiàn)在分析計算簡單，但其缺陷也是明顯的，例如，預測精度或效果往往不盡如人意，而且，普遍適應性差.究其根本原因，由于地基或路基沉降發(fā)展變化的影響因素非常多且沉降機理復雜，難以找到具有普遍適用性的單項模型描述地基或路基沉降發(fā)展變化規(guī)律.其二就是組合模型預測方法，例如，組合預測法、最優(yōu)組合預測模型、非線性組合預測、變權重組合預測法等[8-15]，該方法的基本思路是考慮地基或路基沉降發(fā)展變化規(guī)律的復雜多樣性，利用可描述不同沉降發(fā)展變化規(guī)律的多個單項模型，并依據(jù)各單項模型預測誤差確定權重，實現(xiàn)對地基或路基沉降預測.該方法的突出優(yōu)越性表現(xiàn)在其沉降預測精度或效果明顯優(yōu)于單項模型預測方法，其原因在于組合預測模型能反映地基或路基沉降發(fā)展變化規(guī)律和沉降力學機理的多樣性，彌補了采用單項模型描述地基或路基沉降發(fā)展變化規(guī)律的片面性.因此，組合模型預測已成為目前地基或路基沉降的主流預測方法.

盡管如此，上述沉降組合模型預測方法并未考慮實測數(shù)據(jù)新舊程度對地基或路基沉降預測效果的影響，因此，曹文貴等[16-17]首次引進實測沉降數(shù)據(jù)新鮮度的概念，并采用諸如f（t） = t的新鮮度函數(shù)描述實測沉降數(shù)據(jù)的新舊程度，提出了考慮實測數(shù)據(jù)新舊程度影響的地基或路基沉降預測新方法，這在一定程度上改善了沉降預測效果，但仍然存在明顯不足.一方面，由于其采用的新鮮度函數(shù)是無界的，故無法反映實測數(shù)據(jù)新舊程度影響的客觀有界性，也不能反映實測誤差對新鮮度函數(shù)構建的影響，另一方面，沒有考慮實測數(shù)據(jù)異常的影響.由于人為因素和沉降監(jiān)測設備與方法可能引起的沉降監(jiān)測系統(tǒng)誤差，使不同實測沉降數(shù)據(jù)的可靠程度差別很大，明顯不合理的實測數(shù)據(jù)參與沉降預測必然嚴重影響預測精度或效果，甚至引起沉降預測的謬誤，故在進行地基或路基沉降預測之前，必須首先剔除這些不合理的異常數(shù)據(jù).

綜合上面所述，現(xiàn)有地基或路基沉降預測方法仍然存在較明顯的不足與局限性，為此，本文將在地基或路基沉降組合模型預測方法研究基礎上，重點考慮實測數(shù)據(jù)異常和新舊程度影響有界性對沉降組合預測的影響，對地基或路基沉降組合模型預測方法進行研究，以期完善地基或路基沉降預測的理論與方法.

1 ? 沉降組合預測新方法

設已獲得地基或路基實測沉降數(shù)據(jù)樣本S（ti），它表示ti時刻對應的某點實測沉降（其中，i為該實測點序號，i = 1，2，…，N1;N1為實測點總數(shù)）.如果選取某單項預測模型j，可表示為

因此，由式（3）可得到N2個方程，將其聯(lián)立求解即可確定出各單項預測模型j的全部模型參數(shù)ajk，從而實現(xiàn)采用單項預測模型對t時刻沉降進行預測.

考慮到沉降實測與預測點之間的時間距離不同對地基或路基沉降預測的影響不同以及不同單項預測模型對時間距離影響的依賴性不同，如果對應于單項預測模型j選取不同的反映時間距離即實測數(shù)據(jù)新舊程度影響的新鮮度函數(shù)fj（t），則依據(jù)式（2）可獲得反映實測數(shù)據(jù)新舊程度影響的地基或路基沉降預測的沉降誤差平方和新函數(shù)Q′j，可表示為

至此已建立出考慮實測數(shù)據(jù)新舊程度影響的地基或路基沉降的組合預測模型即式（5）. 值得注意，利用上述模型預測地基或路基沉降能否反映實測數(shù)據(jù)新舊程度及其影響程度有界性以及實測數(shù)據(jù)異常的影響，還須探討實測數(shù)據(jù)新鮮度函數(shù)分析模型的構建方法及實測沉降數(shù)據(jù)異常的處理方法.

2 ? 實測數(shù)據(jù)新鮮度函數(shù)模型構建方法

曹文貴等[16-17]提出了不同單項預測模型采用相同新鮮度函數(shù)f（t）=t來描述實測數(shù)據(jù)新舊程度對沉降預測的影響，其雖然改善了地基或路基沉降預測效果，但明顯存在不足，一方面，對于不同單項預測模型，采用相同新鮮度函數(shù)模型不能反映不同單項預測模型對實測數(shù)據(jù)新舊程度影響的差異性，第二方面，不能反映實測數(shù)據(jù)新舊程度影響的有界性;另一方面，不能反映不同時刻實測沉降與其客觀沉降之間的差異對新鮮度函數(shù)模型的影響，因此，有必要重新探討實測數(shù)據(jù)新鮮度函數(shù)分析模型的構建方法.為此，須首先掌握實測數(shù)據(jù)新鮮度函數(shù)分析模型具有的特性，因此，下面將進一步對此進行探討.

1）反映時間距離即預測時刻t與實測時刻tj之差的影響，而且，時間距離越大，則該實測沉降數(shù)據(jù)對沉降預測的重要性越低，反之亦然，也就是說，時間距離對地基或路基沉降預測的影響呈現(xiàn)單調變化規(guī)律;

2）反映tj時刻實測沉降誤差的影響，而且，誤差越大，該實測點沉降數(shù)據(jù)對沉降預測的重要性越大，反之亦然，也就是說，實測沉降誤差對地基或路基沉降預測的影響也呈現(xiàn)單調變化規(guī)律;

3）不同時刻實測沉降數(shù)據(jù)對地基或路基沉降預測的影響程度是不同的，也就是說，不同實測點實測沉降數(shù)據(jù)的新鮮度是不同的;

4）實測沉降數(shù)據(jù)新鮮度對單項預測模型具有依賴性，也就是說，對于不同單項預測模型，相同時刻實測數(shù)據(jù)新鮮度是不同的;

5）實測沉降數(shù)據(jù)新舊程度對地基或路基沉降預測的影響并非無止境，也就是說，實測數(shù)據(jù)新鮮度函數(shù)分析模型是一個有界函數(shù).

于是，考慮上述實測數(shù)據(jù)新鮮度函數(shù)分析模型的基本特性，經(jīng)深入研究，本文構建出了新型實測數(shù)據(jù)新鮮度函數(shù)模型，可表示為

上述實測數(shù)據(jù)新鮮度函數(shù)分析模型具有如下

特點：

1）由于Eji是根據(jù)單項模型j計算出的ti時刻實測點的沉降誤差，故反映了新鮮度對單項模型的依賴性，也反映了不同測點沉降誤差對新鮮度影響的差異性，同時還反映了新鮮度隨沉降誤差增大而減小的內(nèi)在規(guī)律.

2）由于t和ti分別為預測點和實測點時間，故（ti - t）/t體現(xiàn)了實測點與預測點之間的時間長度大小即實測點沉降數(shù)據(jù)的新舊程度，而且，也由于（ti - t）/t必然為負值，e必然隨ti增大而增大，故也反映了越新的實測點沉降數(shù)據(jù)對沉降預測的影響越大的合理客觀現(xiàn)象.

3）由于0 ≤ ti < t，因此，e必然是一個有界函數(shù)，因此，fj（ti）也必然是一個有界函數(shù)，故本文建立的新鮮度函數(shù)反映了實測沉降數(shù)據(jù)新舊程度對沉降預測影響的有界性.

由此可以看出，本文實測數(shù)據(jù)新鮮度函數(shù)分析模型明顯較曹文貴等[16-17]提出的僅能反映時間距離影響的新鮮度函數(shù)分析模型更具合理性.

3 ? 實測沉降數(shù)據(jù)異常的處理方法

在進行地基或路基沉降監(jiān)測時，由于人為因素或監(jiān)測設備操作與測量方法不當，難以避免實測沉降數(shù)據(jù)的異?；虿徽鎸崳绻蛔魈幚碇苯訉⑷繉崪y數(shù)據(jù)參與地基或路基沉降預測，必然會導致地基或路基沉降預測結果的不合理甚至謬誤.因此，在預測出某時刻t地基或路基沉降之前須剔除這些異常實測沉降數(shù)據(jù).為此，下面將依據(jù)統(tǒng)計學理論探討實測沉降數(shù)據(jù)異常的處理方法.

于是，依據(jù)統(tǒng)計學理論的“3σ”原則[19]，當E″（tk） ∈[μ - 3σ，μ + 3σ]時，tk時刻實測沉降數(shù)據(jù)合理，應參與t時刻沉降組合預測，否則應將其剔除即不參與沉降預測，也就是說，t時刻沉降預測的實測沉降數(shù)據(jù)樣本大小應為N1 - 1.

值得注意，由于σ和μ與參與t時刻沉降組合預測的實測沉降數(shù)據(jù)樣本直接相關，因此，必須采用迭代分析方法進行反復迭代計算，直至完全滿足上述條件為止，此時依據(jù)異常實測沉降數(shù)據(jù)處理后的實測沉降數(shù)據(jù)樣本進行沉降預測，方可獲得最終地基或路基沉降的合理預測結果.另外還須注意，利用“3σ”原則剔除異常數(shù)據(jù)時，需要有足夠大的數(shù)據(jù)樣本，這也是合理沉降預測的基本前提條件，否則，無法進行沉降預測.

4 ? 工程實例分析與驗證

上述已建立出考慮實測數(shù)據(jù)異常與新舊程度影響有界性的地基沉降預測新方法，但其可行性與合理性仍需驗證.為此，本文將采用兩個工程的實測沉降資料進行分析.

工程實例一為杭浦高速公路K90+769～ K111+419段[20]，工程實例二為深圳濱海大道K1+800測試段沉降板T11沉降實測結果[21]，其實測沉降數(shù)據(jù)如表1所示. 由于要驗證本文方法的合理性，只能由理論預測與實測結果進行比較分析，因此，須將兩個工程的實測數(shù)據(jù)分為兩部分，第一部分作為沉降預測的原始數(shù)據(jù)樣本，第二部分作為預測與實測的比較數(shù)據(jù)，因此，對于工程實例一和二，分別取前11和15組數(shù)據(jù)作為原始數(shù)據(jù)，即對應N1分別為11和15.

為了進行沉降預測，首先要選定若干單項預測模型，于是，通過對兩個工程實測沉降曲線變化規(guī)律進行分析，選定Logistic模型[4]、Weibull模型[5]、MMF模型[6]3個單項沉降預測模型進行組合模型預測（即N2 = 3），其中，Logistic模型[4]可表示為

5 ? 結 ? 論

本文在現(xiàn)有組合模型沉降預測方法研究基礎上，重點考慮實測沉降數(shù)據(jù)異常與新舊程度及其影響程度有界性對地基或路基沉降預測精度和效果的影響，進一步深入探討了地基或路基沉降預測的模型與方法. 由此可得如下結論：

1）構建了新型沉降實測數(shù)據(jù)的新鮮度函數(shù)分析模型，它不僅可反映實測沉降數(shù)據(jù)的時間距離與誤差對單項預測模型的依賴性的影響，還可體現(xiàn)不同實測點數(shù)據(jù)的新鮮度函數(shù)不同以及新鮮度的影響程度有界性的特點.

2）基于統(tǒng)計學理論，提出了實測沉降數(shù)據(jù)樣本的數(shù)據(jù)處理方法，避免了異常實測沉降數(shù)據(jù)導致地基或路基沉降預測精度下降甚至謬誤的不合理現(xiàn)象.

3）在上述模型與方法研究中，引進組合預測理論，建立了新型地基或路基沉降預測新方法，并通過工程實例計算以及本文與現(xiàn)有同類方法預測和實測結果的比較分析，表明了本文模型與方法的合理性與優(yōu)越性.

4）值得注意，本文模型與方法的沉降預測精度和效果與所選擇單項預測模型的種類和數(shù)量直接相關，在實際工程應用中，應予以高度重視.

參考文獻

[1] ? ?TAN S A. Hyperbolic method for settlements in clays with vertical drains[J]. Canadian Geotechnical Journal，1994，31（31）：125—131.

[2] ? ?李濤，張儀萍，張土喬. 軟土路基沉降的優(yōu)性組合預測[J]. 巖石力學與工程學報，2005，24（18）：3282—3286.

LI T，ZHANG Y P，ZHANG T Q. Superior combination forecasting settlement of soft clay embankments[J]. Chinese Journal of Rock Mechanics and Engineering，2005，24（18）：3282—3286. （In Chinese）

[3] ? ASAOKA A. Observational procedure of settlement prediction[J]. Soils and Foundations，1978，18（4）：30—34.

[4] ? ?徐洪鐘，施斌，李雪紅. 全過程沉降量預測的Logistic生長模型及其適用性研究[J]. 巖土力學，2005，26（3）：387—391.

XU H Z，SHI B，LI X H. Logistic growth model and its applicability for predicting settlement during the whole process[J]. Rock and Soil Mechanics，2005，26（3）：387—391. （In Chinese）

[5] ? ?MURALEEDHARAN G，NAIR N U，KURUP P G. Application of Weibull model for redefined significant wave height distributions[J]. Proceedings of the Indian Academy of Sciences-Earth and Planetary Sciences，1999，108（3）：149—153.

[6] ? ?廖衛(wèi)紅，王軍保. MMF模型在地基沉降預測中的應用研究[J]. 地下空間與工程學報，2011，7（4）：807—811.

LIAO W H，WANG J B. Study on application of MMF model to prediction of foundation settlement[J]. Chinese Journal of Underground Space and Engineering，2011，7（4）：807—811. （In Chinese）

[7] ? ?路志英，任一墨，葛路琨. 基于樣條估計分位數(shù)回歸的光伏功率回歸模型[J]. 湖南大學學報（自然科學版），2017，44（10）：96—103.

LU Z Y，REN Y M，GE L K. Photovoltaic power regression model based on spline estimation and quantile regression[J]. Journal of Hunan University （Natural Sciences），2017，44（10）：96—103. （In Chinese）

[8] ? ?冷伍明，楊奇，聶如松，等. 高速鐵路橋梁樁基工后沉降組合預測研究[J]. 巖土力學，2011，32（11）：3341—3348.

LENG W M，YANG Q，NIE R S，et al. Study of post-construction settlement combination forecast method of high-speed railway bridge pile foundation[J]. Rock and Soil Mechanics，2011，32（11）：3341—3348. （In Chinese）

[9] ? ?馮震，王娜，王連俊. 沉降組合預測模型研究及其應用[J]. 北京交通大學學報，2008，32（1）：29—31.

FENG Z，WANG N，WANG L J. Research of settlement combination forecasting model and its application[J]. Journal of Beijing Jiaotong University，2008，32（1）：29—31. （In Chinese）

[10] ?BATES J M，GRANGER C W J. The combination of forecasts [J]. Operational Research Quarterly，1969，20（4）：451—468.

[11] ?李猛，謝正文，付艷. 最優(yōu)沉降組合預測模型及應用[J]. 西部探礦工程，2008，20（8）：216—219.

LI M，XIE Z W，F(xiàn)U Y. Optimization settlement combined forecasting model and its application[J]. West-China Exploration Engineering，2008，20（8）：216—219. （In Chinese）

[12] ?趙明華，龍照，鄒新軍. 基于Logistic曲線和Gompertz曲線的最優(yōu)組合沉降預測模型及其應用[J]. 公路交通科技，2007，24（12）：1—4.

ZHAO M H，LONG Z，ZOU X J. Optimal combinational predictive model of settlement based on Logistic curve and Gompertz curve and its application[J]. Journal of Highway and Transportation Research and Development，2007，24（12）：1—4. （In Chinese）

[13] ?歐陽煉. GM（1，1）-Logistic路基沉降組合預測模型研究與應用[J]. 鐵道科學與工程學報，2010，7（4）：56—60.

OU Y L. Research and application of GM（1，1） and Logistic combination forecasting model for subgrade settlement[J]. Journal of Railway Science and Engineering，2010，7（4）：56—60. （In Chinese）

[14] ?吳清海，李慧芳. 變權重組合模型在沉降預測中的應用[J]. 測繪科學技術學報，2009，26（2）：118—124.

WU Q H，LI H F. Application of the changeable weight combination model in building settlement predicting[J]. Journal of Geomatics Science and Technology，2009，26（2）：118—124. （In Chinese）

[15] ?趙明華，劉煜，曹文貴. 軟土路基沉降變權重組合S型曲線預測方法研究[J]. 巖土力學，2005，26（9）：1443—1447.

ZHAO M H，LIU Y，CAO W G. Study on variable-weight combination forecasting method of S-type curves for soft clay embankment settlement[J]. Rock and Soil Mechanics，2005，26（9）：1443—1447. （In Chinese）

[16] ?曹文貴，印鵬，賀敏，等. 考慮實測數(shù)據(jù)新舊程度的工后沉降單項模型預測方法[J]. 水文地質工程地質，2015，42（6）：65—70.

CAO W G，YIN P，HE M，et al. A prediction method for post-construction settlement of a single model with the consideration of the new or old degree of the measured data[J]. Hydrogeology & Engineering Geology，2015，42（6）：65—70. （In Chinese）

[17] ?曹文貴，印鵬，賀敏，等. 基于數(shù)據(jù)新舊程度和預測取值區(qū)間調整的沉降組合預測方法[J]. 巖土力學，2017，38（2）：534—540.

CAO W G，YIN P，HE M，et al. A combination method for predicting settlement based on new or old degree of data and adjustment of value interval of prediction[J]. Rock and Soil Mechanics，2017，38（2）：534—540. （In Chinese）

[18] ?卿啟湘，王永和，李光耀. 標準擊實試驗的一種數(shù)據(jù)處理方法[J]. 湖南大學學報（自然科學版），2004，31（2）：48—51.

QING Q X，WANG Y H，LI G Y. A numerical method for analysis of the result to indoor standard compaction test[J]. Journal of Hunan University （Natural Sciences），2004，31（2）：48—51. （In Chinese）

[19] ?周建普，李獻民，王永和. 粘土邊坡可靠性分析方法研究[J]. 湖南大學學報（自然科學版），2002，29（5）：92—97.

ZHOU J P，LI X M，WANG Y H. Research for method of reliability analysis for clayey soil slop[J]. Journal of Hunan University （Natural Sciences），2002，29（5）：92—97. （In Chinese）

[20] ?顧建武，何良德，李林青. 基于速率-沉降特征分析工后沉降預測方法的適用性[J]. 交通科技，2012（2）：85—88.

GU J W，HE L D，LI L Q. Applicability of post-construction settlement prediction model based on velocity-settlement relationship[J]. Transportation Science & Technology，2012（2）：85—88. （In Chinese）

[21] ?上海市政工程設計研究院. 深圳市濱海大道市政工程軟基處理分析報告[R]. 上海：上海市政工程設計研究院，1998：16—38.

Shanghai Municipal Engineering Design & Research Institute. The analysis report for soft ground treatment of costal road of Shenzhen municipal engineering [R]. Shanghai：Shanghai Municipal Engineering Design & Research Institute，1998：16—38. （In Chinese）