李 文

(1.湖南理大交通科技發(fā)展有限公司， 湖南 長沙 410006; 2.湖南交通職業(yè)技術(shù)學(xué)院， 湖南 長沙 410006)

0 引言

目前，公路路基沉降的預(yù)測主要是根據(jù)積累的實(shí)測數(shù)據(jù)，通過合適的模型進(jìn)行擬合分析，實(shí)現(xiàn)沉降預(yù)測，其預(yù)測效果主要取決于實(shí)測數(shù)據(jù)的準(zhǔn)確性和模型選取的合理性。路基沉降變化在不同施工階段內(nèi)有明顯差異，因此，在構(gòu)建預(yù)測模型時可以考慮對多個單一預(yù)測模型進(jìn)行組合，通過合理的權(quán)重分配來綜合各單項(xiàng)模型優(yōu)勢，取長補(bǔ)短[1-4]。本文根據(jù)項(xiàng)目路基沉降曲線的變化特點(diǎn)，研究不同施工階段內(nèi)單一預(yù)測模型的擬合效果和誤差，對模型精度和適用性進(jìn)行分析，以期找到一種科學(xué)、準(zhǔn)確的預(yù)測方法。

1 工程概況

湖北省某在建高速公路，其路線地處長江中下游地區(qū)，所經(jīng)區(qū)域以農(nóng)田為主，局部溝塘較多，部分路段為軟土地基。經(jīng)勘查分析，原地基存在軟弱夾層，不同路段地基土性質(zhì)差異性較大，主要不良地質(zhì)情況為軟土，該類土以淤泥、淤泥質(zhì)土等為主，具有含水量高、壓縮性高、強(qiáng)度低等特點(diǎn)，由于地基在附加應(yīng)力作用下發(fā)生排水固結(jié)，容易產(chǎn)生較大的沉降變形或差異沉降。

本項(xiàng)目軟基采取了分段、分措施處治，為了使沉降數(shù)據(jù)具有代表性，更好地反映本項(xiàng)目路基沉降的規(guī)律和特性，分別選取了3個典型填方斷面的沉降實(shí)測數(shù)據(jù)進(jìn)行分析，典型斷面基本情況如表1所示。

表1 典型斷面基本信息典型斷面樁號填土高度/m處治方案1K12+4023.72砂墊層加超載預(yù)壓2K12+4764.58水泥粉噴樁3K12+5054.83粉噴樁加等載預(yù)壓

軟土路基在填筑過程中，隨著填土高度增加，荷載不斷增大，地基在附加應(yīng)力作用下會發(fā)生排水固結(jié)，使整個路基出現(xiàn)持續(xù)沉降變形。其中，中樁位置沉降量要大于路肩處，因此可選取斷面中樁位置處的沉降數(shù)據(jù)對沉降過程和規(guī)律進(jìn)行分析，典型斷面觀測點(diǎn)的實(shí)測累計(jì)沉降量與時間變化曲線如圖1所示。

圖1 實(shí)測累計(jì)沉降-時間變化曲線

由圖1沉降變化曲線可知，不同處治措施的3個點(diǎn)位處的沉降變化趨勢基本相同。其累計(jì)沉降變化過程主要分為3個階段，分別為填筑期、恒載期和穩(wěn)定期。路基填筑期為0～120 d，該階段路基沉降速度較快，變化明顯；恒載期為路基填筑完成后的3個月左右(120～220 d)，沉降速率逐漸降低；穩(wěn)定期為220～380 d，此時沉降速率較小，沉降曲線也開始收斂，最終在380 d后達(dá)到基本穩(wěn)定狀態(tài)。

2 單一預(yù)測模型效果分析

結(jié)合本項(xiàng)目沉降曲線特點(diǎn)，參考同類型項(xiàng)目經(jīng)驗(yàn)，選擇項(xiàng)目所在地區(qū)工程實(shí)踐中較為常用的2種單一模型，即灰色GM(1，1)模型和Verhulst模型，分別就不同施工階段的沉降曲線進(jìn)行擬合，對預(yù)測模型的特點(diǎn)和適用性進(jìn)行分析。

2.1 單一預(yù)測模型

2.1.1GM(1，1)模型

GM(1，1)其預(yù)測值可以認(rèn)為是預(yù)測曲線下最優(yōu)曲線的延伸，該模型具有無限增長的性質(zhì)，模型表達(dá)式如下：

k=2，3，…，n

(1)

2.1.2灰色Verhulst模型

灰色Verhulst模型曲線呈“S”形，可反映路基沉降由前期加速變形到最終趨于穩(wěn)定的趨勢。Verhulst模型離散響應(yīng)方程如下[6]：

(k=1，2，…，n)

(2)

式中:a、b為擬合參數(shù)，S(1)為原始沉降觀測數(shù)據(jù)，S(0)為S(1)經(jīng)一次累減生成新數(shù)列，S(k)為第k次沉降數(shù)據(jù)。

2.2 實(shí)測數(shù)據(jù)等時距處理

由于現(xiàn)場實(shí)測沉降數(shù)據(jù)的時間間隔不一致，需要對原始數(shù)據(jù)進(jìn)行等時距處理，本文采用拉格朗日插值方法，分別計(jì)算各個斷面等時距各天的沉降值[5]。

非等時距沉降增量時間序列為：

(3)

(4)

(5)

t=n時：

(6)

當(dāng)t=2，3，…，n-1時，利用拉格朗日插值函數(shù)分段線形插值計(jì)算如下：

(7)

因此，等時距沉降數(shù)據(jù)時間序列：

(8)

按上述時距處理方法，選取沉降變化較大且特征較為明顯的3#斷面進(jìn)行分析，以20 d為單位對原始數(shù)據(jù)等時距處理，等時距沉降數(shù)據(jù)如表2。

表2 3#斷面等時距處理后的沉降值序號時間/d累計(jì)沉降/mm序號時間/d累計(jì)沉降/mm12031.4611220108.7024056.0212240109.7136069.7413260110.8348080.3314280111.61510088.0315300111.96612094.1816320112.20714098.6717340112.708160102.2718360113.349180105.2719380114.0010200107.50

2.3 單一預(yù)測模型階段預(yù)測效果

分別選取不同施工階段內(nèi)，累積實(shí)測數(shù)據(jù)序列中靠前的數(shù)據(jù)為樣本計(jì)算模型參數(shù)，并預(yù)留1組數(shù)據(jù)與模型預(yù)測值進(jìn)行比對對單個預(yù)測模型的適用性和效果進(jìn)行分析。

對填筑期6組沉降數(shù)據(jù)進(jìn)行回歸分析(見表3),模型擬合相對誤差對比見圖2。

表3 填筑期單一預(yù)測模型分析結(jié)果序號天數(shù)/d實(shí)測值/mmGM(1,1)模型Verhulst模型擬合值/mm相對誤差/%擬合值/mm相對誤差/%12031.4631.46—31.460 0—24056.0258.594.5951.966 7-7.2436069.7467.54-3.1569.972 50.3348080.3377.85-3.0981.060 10.91510088.0389.731.9386.469 1-1.77612094.18103.43(預(yù)測值)9.8289.767 2(預(yù)測值)-4.69

圖2 填筑期預(yù)測模型擬合相對誤差

由表3可以看出，路基填筑階段采用GM(1,1)模型和Verhulst模型進(jìn)行擬合，二者擬合值的相對誤差均保持在10%以內(nèi)，整體預(yù)測效果較好。GM(1，1)模型擬合值變化速率較大，主要表現(xiàn)為正向偏大，而Verhulst模型擬合值整體偏小且變化較平穩(wěn)。根據(jù)圖2，對模型的相對誤差進(jìn)行對比發(fā)現(xiàn)，2種模型的相對誤差正負(fù)走向相反，存在誤差的互補(bǔ)可能性。

路基持續(xù)填筑至設(shè)計(jì)高度后停止填土，直至鋪筑路面前，此時荷載處于相對穩(wěn)定，進(jìn)入恒載狀態(tài)，對該階段內(nèi)的12組沉降數(shù)據(jù)進(jìn)行回歸分析(見表4)，模型擬合相對誤差對比見圖3。

由表4和圖3可以看出，進(jìn)入恒載期后，隨著實(shí)測數(shù)據(jù)的增多，GM(1，1)和Verhulst模型的平均相對誤差分別為6.70%和4.51%，擬合效果有一定的提升。對單點(diǎn)預(yù)測值(12號)進(jìn)行分析可知，GM(1,1)模型的預(yù)測值為123.43 mm，相對誤差為12.51%，Verhulst模型預(yù)測值為106.51 mm，相對誤差為-2.92%，從誤差絕對值上看，GM(1,1)模型預(yù)測值明顯偏大，Verhulst模型更接近實(shí)測值，但預(yù)測值較小，結(jié)果偏保守。

表4 恒載期單一預(yù)測模型分析結(jié)果序號天數(shù)/d實(shí)測值/mmGM(1,1)模型Verhulst模型擬合值/mm相對誤差/%擬合值/mm相對誤差/%12031.4631.46—31.46—24056.0269.5224.1046.60-16.8236069.7473.635.5862.92-9.7848080.3377.98-2.9377.58-3.42510088.0382.59-6.1888.720.78612094.1887.47-7.1296.182.12714098.6792.63-6.12100.742.108160102.2798.11-4.07103.381.099180105.27103.90-1.30104.87-0.3810200107.50110.042.36105.69-1.6811220108.70116.547.21106.13-2.3612240109.71123.43(預(yù)測值)12.51106.51(預(yù)測值)-2.92

圖3 恒載期預(yù)測模型擬合相對誤差

在施工完成后路基進(jìn)入穩(wěn)定期并逐漸達(dá)到基本穩(wěn)定狀態(tài)，此時沉降速率逐漸趨近為0，累積沉降曲線尾端開始收斂，對路基填筑至穩(wěn)定期內(nèi)的全過程的18組沉降數(shù)據(jù)進(jìn)行回歸分析(見表5)，模型擬合相對誤差對比見圖4。

由表5和圖4分析可知，以該斷面路基沉降全過程數(shù)據(jù)為樣本進(jìn)行預(yù)測時，2種模型的早期擬合數(shù)據(jù)變化速率較快，具有較明顯上升趨勢，中后期變化速率相對降低，但GM(1,1)模型仍保持持續(xù)增長，Verhulst模型則逐漸收斂。從后3組預(yù)測值可以看出，GM(1,1)模型預(yù)測值增長速度仍然較大，而Verhulst模型增速則趨近0，與現(xiàn)場實(shí)際情況相符，整體上Verhulst模型擬合效果較好，更適用于中長期預(yù)測。

表5 沉降全過程單一預(yù)測模型分析結(jié)果序號天數(shù)/d實(shí)測值/mmGM(1,1)模型Verhulst模型擬合值/mm相對誤差/%擬合值/mm相對誤差/%12031.4631.46—31.46—24056.0279.3641.6645.09-19.5136069.7481.5716.9660.15-13.7548080.3383.844.3774.52-7.23510088.0386.18-2.1086.46-1.78612094.1888.58-5.9595.281.17714098.6791.04-7.73101.242.608160102.2793.58-8.50105.042.719180105.2796.18-8.63107.371.9910200107.5098.86-8.04108.761.1711220108.70101.61-6.52109.580.8112240109.71104.44-4.80110.060.3213260110.83107.35-3.14110.34-0.4414280111.61110.33-1.15110.50-0.9915300111.96113.411.30110.60-1.2116320112.20116.563.89110.65-1.3817340112.70119.816.31110.68-1.7918360113.34123.148.65110.70-2.3319380114.00126.571 2(預(yù)測值)11.03110.717 4(預(yù)測值)-2.8820400-130.094 8(預(yù)測值)110.720 9(預(yù)測值)21420-133.716 4(預(yù)測值)110.722 9(預(yù)測值)

圖4 沉降全過程預(yù)測模型擬合相對誤差

綜上所述，在不同路基施工階段中沉降表現(xiàn)出一定的階段性，GM(1，1)模型能夠較好的反映路基填筑期沉降變化，Verhulst模型則更適合中長期預(yù)測。在不同的沉降階段，2種預(yù)測模型相對誤差變化趨勢不同，GM(1，1)模型誤差較大，且主要為正向增加，Verhulst模型整體誤差較小，擬合值也相對偏小，二者適用范圍且具有一定的互補(bǔ)性。

3 熵權(quán)法組合預(yù)測模型

3.1 組合模型構(gòu)建

通過對2種模型擬合值效果和誤差進(jìn)行分析，考慮預(yù)測擬合值與實(shí)測值之間的殘差及相對誤差之間具有不確定性特點(diǎn)，采用熵權(quán)法來構(gòu)建組合模型，提高預(yù)測精度和模型的適用范圍。熵權(quán)法的基本原理是用熵來表征系統(tǒng)無序程度，如果指標(biāo)的信息熵越小，該指標(biāo)提供的信息量越大且權(quán)重越大。熵權(quán)法中對于某預(yù)測問題，選用的第i種單項(xiàng)預(yù)測模型在第t時刻得到的預(yù)測值與實(shí)測值之間的相對誤差可表達(dá)為：

(9)

式中：i=(1，2，3,…,m)；t=(1，2，3,…,n)；xit為第i種單項(xiàng)預(yù)測方法在第t時刻得到的預(yù)測值。

權(quán)值具體計(jì)算方法如下：

1) 計(jì)算第i種單項(xiàng)預(yù)測方法的預(yù)測值與實(shí)測值在t時刻的相對誤差比重

(10)

(11)

2) 根據(jù)熵值與其變異程度相反的原則，求相對誤差變異系數(shù)

di=1-hi

(12)

3) 再得單項(xiàng)預(yù)測模型的加權(quán)值

(13)

由式(13)計(jì)算出權(quán)系數(shù)后，組合預(yù)測模型表達(dá)式為：

(14)

單項(xiàng)預(yù)測誤差序列的熵值越大，其變異程度越小、權(quán)值越大。熵權(quán)法充分考慮了不同預(yù)測方法相對誤差的大小，對組合預(yù)測模型的權(quán)重進(jìn)行分配，使單一模型之間能夠相互補(bǔ)充，權(quán)值確定結(jié)果，見表6。

3.2 預(yù)測模型效果對比分析

本文選用的預(yù)測模型均屬于曲線回歸法，即基于實(shí)測數(shù)據(jù)采用回歸計(jì)算的方法反算模型參數(shù)，可采用相關(guān)系數(shù)和相對誤差對模型精度進(jìn)行評價，根據(jù)熵權(quán)法確定GM(1，1)模型和Verhulst模型的權(quán)重值，計(jì)算相關(guān)系數(shù)R和預(yù)測值的相對誤差結(jié)果見表7。

表6 權(quán)值計(jì)算結(jié)果預(yù)測模型權(quán)值分配填筑期填筑-恒載期全過程GM(1,1)模型0.3590.3960.403Verhulst模型0.6410.6040.597

表7 組合模型和單一模型效果評定預(yù)測模型填筑期填筑-恒載期全周期相關(guān)系數(shù)R相對誤差/%相關(guān)系數(shù)R相對誤差/%相關(guān)系數(shù)R相對誤差/%GM(1,1)模型0.970 5 9.820.901 2 10.300.837 8 11.03Verhulst模型0.987 9 -4.690.991 2 -2.920.990 8 -2.88組合模型0.999 7 1.060.994 2 2.320.993 9 2.02

由表7可知，組合模型在路基沉降各個時期中較單一模型而言具有良好的擬合效果，且預(yù)測值相對誤差保持在2%左右，相關(guān)系數(shù)R值保持在0.99以上，具有較高的預(yù)測精度。組合模型曲線變化由權(quán)系數(shù)分配較大的單一預(yù)測模型控制，另一種模型作為完善和補(bǔ)充，使組合模型曲線與實(shí)測結(jié)果更為接近，體現(xiàn)了組合預(yù)測模型的優(yōu)越性。此外，組合模型可根據(jù)實(shí)測數(shù)據(jù)的變化特點(diǎn)，對權(quán)重分配和調(diào)整，具備良好的適用性。

4 結(jié)論

本文根據(jù)路基沉降實(shí)測數(shù)據(jù)變化特點(diǎn)，按施工階段的不同，對灰色GM(1,1)和Verhulst 2種模型的擬合效果和適用性進(jìn)行了分析，并對預(yù)測值與實(shí)測值之間相對誤差進(jìn)行了研究，構(gòu)建了組合模型，主要結(jié)論如下：

1) 該項(xiàng)目路基累計(jì)沉降變化具有明顯的階段變化，GM(1，1)模型能夠較好的反映路基填筑期的沉降變化，而Verhulst模型則更適用于中長期預(yù)測。從模型整體的相對誤差進(jìn)行對比發(fā)現(xiàn)GM(1，1)模型相對誤差較大，Verhulst模型誤差較小，2種中預(yù)測模型相對誤差變化不同，具有一定的互補(bǔ)性。

2) 基于熵權(quán)法構(gòu)建的組合預(yù)測模型具有可行性和有效性，權(quán)值確定過程中考慮單一預(yù)測模型中預(yù)測值與實(shí)測值的相對誤差不確定性，綜合了單一模型的優(yōu)勢，預(yù)測效果有明顯提升。驗(yàn)證結(jié)果表明，組合模型在路基沉降的各階段中均有較好地?cái)M合效果，預(yù)測值相對誤差保持在2%左右。同時，組合模型可根據(jù)項(xiàng)目沉降曲線的變化特點(diǎn)，對權(quán)重進(jìn)行靈活調(diào)整，具備更好的適用性，可為同類型路基沉降預(yù)測提供參考。