Asaoka法預(yù)測軟基最終沉降的改進算法

鄧雄文

（廣州市地質(zhì)調(diào)查院，廣東廣州510440）

Asaoka法預(yù)測軟基最終沉降的改進算法

鄧雄文*

（廣州市地質(zhì)調(diào)查院，廣東廣州510440）

Asaoka法計算軟土地基最終沉降量的圖解過程中，通常對實測沉降數(shù)據(jù)點系列（sj-1、sj）采用直線擬合，通過最佳擬合直線與45°線相交來求取最終沉降量。針對數(shù)據(jù)點較分散時可能出現(xiàn)較大預(yù)測誤差的缺陷，采用非線性方法擬合數(shù)據(jù)點并結(jié)合工程實例驗證了改進算法的精度，具有一定參考和使用價值。

沉降預(yù)測；Asaoka法；改進算法

1　概述

軟土地基沉降預(yù)測可以有效指導(dǎo)施工、預(yù)防沉降地質(zhì)災(zāi)害的發(fā)生，具有重大的經(jīng)濟效益和社會效益。從軟土基地的實測沉降數(shù)據(jù)中尋找沉降和時間的擬合關(guān)系是預(yù)測最終沉降量的一種有效方法，常用的有雙曲線法、指數(shù)曲線法、泊松曲線法和Asaoka法等。本文主要基于Asaoka法計算過程中的不足，將實測沉降數(shù)據(jù)點系列（sj-1、sj）的擬合由線性改進為非線性，由一元一次改進為一元多次。

2　Asaoka法基本原理

Mikasa一維固結(jié)方程在一維條件下以體積應(yīng)變表示的固結(jié)方程如下：

式中：εv——與時間t和排水距離z有關(guān)的豎向應(yīng)變；

cv——固結(jié)系數(shù)。

Asaoka將式（1）近似的以級數(shù)形式的微分方程表示為：

式中：s——固結(jié)沉降量（包括瞬時沉降、主固結(jié)沉降、次固結(jié)沉降）；

α1、α2、…、αn——固結(jié)系數(shù)；

b——取決于土層邊界條件的常數(shù)。

Asaoka法的目標就是利用實測沉降資料求出這些常數(shù)，并利用這些常數(shù)求出最終沉降量和固結(jié)度。式（2）可用遞推關(guān)系表示如下：

式中：β0——沉降值；

βi——無維度常量。

上式中當(dāng)n=1時，解得：

式中：s0、s∞——初始和最終沉降量。

最終沉降量s∞的求解多采用圖解法，如圖1所示。將沉降監(jiān)測時間分為相等的時間段Δt，在實測沉降曲線中讀出t1、t2…對應(yīng)的沉降量s1、s2…，在sj-1-sj坐標系中繪出一系列數(shù)據(jù)點（）sj-1，sj，取最佳擬合直線與坐標系45°線交點處為預(yù)測最終沉降量s∞。

3　Asaoka法的改進

圖解過程中，實測數(shù)據(jù)點系列（sj-1、sj）在sj-1-sj坐標系上可能是較為分散的，有些情況下甚至本身就是非線性的，基于這種情況，本文采用非線性擬合的方法通過一個n次多項式的形式擬合數(shù)據(jù)點，取擬合曲線和45°線最合理交點為最終沉降量。其基本原理如圖2所示。

對于m組數(shù)據(jù)（xi、yi）（i=1，2，3…，m）設(shè)擬合的曲線多項式為：

當(dāng)n=1時即為常用的線性擬合。

圖1　Asaoka法示意圖

圖2　Asaoka法改進算法原理圖

可以解得k0、k1…kn的值，代入式（5）聯(lián)立y=x可得一多項式：令F（）

x=0，因為sj≥sj-1恒成立，結(jié)合曲線走勢及其凹凸性可以判斷當(dāng)n為偶數(shù)時方程F（）x=0可能無解，即擬合的曲線與45°線無交點。

當(dāng)n為奇數(shù)時方程F（）x=0的解可視為最終預(yù)測沉降量。

4　工程應(yīng)用

某小區(qū)工程建在分布有4層全新統(tǒng)海陸交互相淤泥、淤泥質(zhì)土的軟土地基上，從2011年1月18日起進行了約2年的沉降監(jiān)測。共設(shè)置了5個沉降監(jiān)測點，其平均沉降量隨時間的變化關(guān)系如圖3所示。

圖3　沉降—時間曲線

選取60d為時間間隔的8組實測數(shù)據(jù)點（即Δt=60，m=8），以常用的Asaoka法及其改進方法分別擬合數(shù)據(jù)點，作sj-1-sj曲線關(guān)系如圖4和圖5所示。

圖4　線性擬合sj-1-sj示意圖

經(jīng)過計算，線性擬合及各曲線擬合情況如表1所示。

圖5　非線性擬合sj-1-sj示意圖

表1　數(shù)據(jù)點擬合多項式

當(dāng)n為偶數(shù)時，方程無解，曲線與45°線無交點。線性擬合計算的最終沉降量是三次曲線擬合計算的最終沉降量的2.61倍，是五次曲線擬合計算的最終沉降量的2.41倍。事實上該小區(qū)經(jīng)過近8年的沉降固結(jié)，已經(jīng)進入較勻速的下沉過程且沉降速度極小，最大沉降速率與最小沉降速率之差僅為0.019mm/d，不均勻沉降很小。根據(jù)監(jiān)測數(shù)據(jù)綜合分析認為曲線擬合計算結(jié)果更加符合實際，預(yù)測精度大于常用的直線擬合計算方法，且n值越大，精度越高。據(jù)此，該小區(qū)最終預(yù)測沉降量為29.72mm。

表2　最終沉降量預(yù)測值

5　結(jié)論

［1］王志亮，黃景忠，李永池.沉降預(yù)測中的Asaoka法應(yīng)用研究［J］.巖土力學(xué)，2006，31（7）：2025-2028.

［2］潘林有，謝新宇.用曲線擬合的方法預(yù)測軟土地基沉降［J］.巖土力學(xué)，2007，26（7）：1053-1058.

［3］李凡，吳志昂，孫四平.軟土地基路堤最終沉降量推算方法研究［J］.合肥工業(yè)大學(xué)學(xué)報，2003，26（2）：277-280.

TU196

A

1004-5716（2016）02-0033-03

2015-02-05

鄧雄文（1969-），男（漢族），廣東五華人，高級工程師，現(xiàn)從事測繪技術(shù)工作。