吳會(huì)閣王 麗，葛建龍

(1.石家莊經(jīng)濟(jì)學(xué)院 工程學(xué)院，河北 石家莊 050031;2.河北省水資源可持續(xù)利用與開(kāi)發(fā)重點(diǎn)實(shí)驗(yàn)室，河北 石家莊050031;3.北京中鐵建筑工程設(shè)計(jì)院 天津分院，天津300011;4.天津振津工程集團(tuán)有限公司，天津300222)

0 前言

在越來(lái)越多的高層建筑和高重構(gòu)筑物設(shè)施中，群樁基礎(chǔ)以其承載力高、沉降小、生命力強(qiáng)等特點(diǎn)得到了廣泛的應(yīng)用。巖土工程參數(shù)的正確選取是群樁基礎(chǔ)可靠性的重要保證，在各參數(shù)中土的物理性質(zhì)和力學(xué)性質(zhì)的空間變異性最大，描述這種空間變異性較好的方法是文獻(xiàn)［1］提出的土工剖面的隨機(jī)場(chǎng)理論，而反映土工參數(shù)空間變異特性的重要指標(biāo)是土層的相關(guān)范圍。近年來(lái)，國(guó)內(nèi)外一些學(xué)者對(duì)土工參數(shù)的相關(guān)范圍進(jìn)行了一系列研究，文獻(xiàn)［2］用一維平穩(wěn)高斯隨機(jī)場(chǎng)理論，研究了用回歸模擬方法計(jì)算土工參數(shù)的相關(guān)范圍;文獻(xiàn)［3］用改進(jìn)的空間遞推平均法討論了土工參數(shù)相關(guān)范圍的確定;文獻(xiàn)［4］用非線性曲線擬合法計(jì)算了土性的相關(guān)范圍;文獻(xiàn)［5］用相關(guān)函數(shù)法討論了相關(guān)范圍確定中相關(guān)函數(shù)的選取及擬合區(qū)間;文獻(xiàn)［6］對(duì)目前確定土工參數(shù)相關(guān)范圍的幾種方法進(jìn)行了評(píng)述，指出空間遞推平均法雖計(jì)算方便，但平穩(wěn)點(diǎn)難以確定，其計(jì)算結(jié)果受人為因素影響較大，曲線擬合法算出的相關(guān)距離值偏小。

在對(duì)群樁基礎(chǔ)的可靠度分析時(shí)，本文擬采用相關(guān)函數(shù)法和變異函數(shù)法計(jì)算土工參數(shù)的相關(guān)范圍。

1 土工參數(shù)的空間變異性

在天然土層中，土性指標(biāo)在各點(diǎn)是確定的，根本不存在隨機(jī)性的問(wèn)題，“唯一”不確定的問(wèn)題是無(wú)法知道土工參數(shù)在土層中每一點(diǎn)的取值［7］。

土工參數(shù)的空間變異由以下3部分組成:

其中，E［X(Z)］為土工參數(shù)的趨勢(shì)分量，是非隨機(jī)的變量，可由最小二乘法回歸;X'(Z)為零均值，方差為σ2的正態(tài)分布隨機(jī)變量，{X'(Z)}為一平穩(wěn)隨機(jī)場(chǎng)，代表土性固有變異性;ε(Z)為測(cè)量誤差。

對(duì)于平穩(wěn)的隨機(jī)場(chǎng)，其統(tǒng)計(jì)特征不隨位置的不同而不同，其均值、方差以及自相關(guān)函數(shù)都可以通過(guò)樣本相應(yīng)的平均(隨機(jī)過(guò)程的時(shí)間平均)得到。對(duì)于有連續(xù)變化的X'(Z)而言，由Shannon抽樣定理，這時(shí)要求采樣能再現(xiàn)原來(lái)連續(xù)曲線中含有的全部信息，從這一角度出發(fā)，取樣間距越小、采樣點(diǎn)越多越好;但取樣間距太小，會(huì)增加點(diǎn)均值的統(tǒng)計(jì)誤差、低估點(diǎn)的變異性。當(dāng)取樣間距小于相關(guān)范圍時(shí)，由于土性自相關(guān)性的影響，用傳統(tǒng)的獨(dú)立、同分布序列統(tǒng)計(jì)方法計(jì)算“點(diǎn)”均值、“點(diǎn)”方差以及自協(xié)方差函數(shù)的無(wú)偏估計(jì)時(shí)，要加以修正。

2 兩種算法的基本原理

2.1 相關(guān)函數(shù)法

2.1.1 相關(guān)函數(shù)理論［8－10］

對(duì)于一均質(zhì)土層某一點(diǎn)t的土性記作x(t)，則x(t)的概率特性可表示為:

那么，土層中任意兩點(diǎn)t，t+h之間的相關(guān)函數(shù)可用下式表示:

式中Γ2(h)為方差折減系數(shù)。

常用的理論相關(guān)函數(shù)有三角型、指數(shù)型、指數(shù)余弦型等，其對(duì)應(yīng)的函數(shù)形式和相關(guān)范圍(或變程)計(jì)算公式見(jiàn)表1。

表1 理論相關(guān)函數(shù)模型

2.1.2 相關(guān)范圍的確定

用相關(guān)函數(shù)確定相關(guān)范圍，就是用一定類型的理論相關(guān)函數(shù)模型去擬合實(shí)測(cè)相關(guān)函數(shù)曲線，求出相應(yīng)參數(shù)，代入對(duì)應(yīng)的相關(guān)范圍計(jì)算公式求得結(jié)果。

相關(guān)函數(shù)法求解相關(guān)范圍的步驟:

(1)原始數(shù)據(jù)的處理:將試驗(yàn)測(cè)得的原始數(shù)據(jù)X0i標(biāo)準(zhǔn)化后作為樣本值Xi，即:

式中，μ，σ分別為原始數(shù)據(jù)的均值和標(biāo)準(zhǔn)差。處理后，樣本均值為0，方差為1。

(2)取間距τ1=△Z，△Z為取樣間距，計(jì)算相關(guān)函數(shù)值:

(3)分別取間距 τ2=2τ1，τ3=3τ1，…，τm=mτ1，…(m ＜ n)，分別計(jì)算出此時(shí)的相關(guān)函數(shù)值:

(4)此時(shí)，所得到的相關(guān)函數(shù)為實(shí)測(cè)的相關(guān)函數(shù)，可做出實(shí)測(cè)相關(guān)函數(shù)曲線，根據(jù)實(shí)測(cè)相關(guān)函數(shù)曲線選擇合適的理論相關(guān)函數(shù)模型，對(duì)實(shí)測(cè)曲線進(jìn)行擬合，把擬合出的相應(yīng)參數(shù)代入計(jì)算公式即可。

2.2 變異函數(shù)法

2.2.1 變異函數(shù)理論

土的空隙度、相對(duì)密度、塑性指標(biāo)、滲透系數(shù)、壓縮模量、抗剪(壓)強(qiáng)度以及某一特定持力土層的厚度都可看作是區(qū)域化變量。區(qū)域化變量一般可反映巖土參數(shù)的尺寸效應(yīng)，不同程度的連續(xù)性和不同種類的異向性及空間變化的可遷性等特征。

對(duì)于一個(gè)代表巖土參數(shù)的區(qū)域化變量Z(x)，變異函數(shù)可定義為向量h相隔的Z(x)、Z(x+h)兩個(gè)區(qū)域化變量間增量平方的數(shù)學(xué)期望的一半，或兩區(qū)域化變量增量的方差的一半(也稱半變異函數(shù))，本文仍稱之為變異函數(shù):

由上式可知:變異函數(shù)是點(diǎn)x和向量h兩者的函數(shù)。由于在任意點(diǎn)上不可能獲得一系列樣本值，因此要求Z(x)滿足本征假設(shè)條件，即可認(rèn)為［Z(x)－Z(x+h)］只依賴于h，而不依賴于x，則可用實(shí)際觀測(cè)值計(jì)算變異函數(shù)的估計(jì)量，即:

式中，N(h)是以向量h相隔的試驗(yàn)數(shù)據(jù)對(duì)Z(xi)和Z(xi+h)的數(shù)目。

常用的理論變異函數(shù)模型有球狀模型、指數(shù)模型、高斯(Gauss)模型等，其中最常用的是球狀模型。各種模型的函數(shù)形式及其對(duì)應(yīng)的相關(guān)范圍(或變程)值見(jiàn)表2和圖1。

表2 理論半變異函數(shù)類型

圖1 理論變異函數(shù)模型

2.2.2 相關(guān)范圍的確定

與相關(guān)函數(shù)法相似，變異函數(shù)法是用理論變異函數(shù)去擬合實(shí)測(cè)變異函數(shù)曲線，從而得出相應(yīng)的參數(shù)，代入公式求得相關(guān)范圍的值。

變異函數(shù)法求解相關(guān)范圍的步驟:

(1)試驗(yàn)測(cè)得的原始數(shù)據(jù)作為樣本值{X(t)}，這時(shí)，樣本均值為 μ，方差為 σ2。

(2)△Z為取樣間距，取間距τ1=△Z，計(jì)算這時(shí)的變異函數(shù)值:

其中N1表示相距τ1的樣本點(diǎn)對(duì)數(shù)。

(4)此時(shí)，所得到的變異函數(shù)為實(shí)測(cè)變異函數(shù)，可做出實(shí)測(cè)變異函數(shù)曲線，根據(jù)實(shí)測(cè)變異函數(shù)曲線形式選擇合適的理論變異函數(shù)模型對(duì)實(shí)測(cè)曲線進(jìn)行擬合，求出參數(shù)后，便得到相關(guān)范圍(或變程)值。

3 工程計(jì)算實(shí)例

圖2為實(shí)測(cè)的天津地區(qū)某臨近場(chǎng)地的靜力觸探曲線，fs為雙橋探頭測(cè)得的側(cè)摩阻，kPa;qc為作用于探頭上的阻力，MPa;h為相關(guān)范圍，m。

圖2 試樁靜力觸探曲線

基于指數(shù)模型，對(duì)各孔土層的相關(guān)范圍分別采用相關(guān)函數(shù)法和變異函數(shù)法進(jìn)行了計(jì)算。以J-5孔為例，具體計(jì)算結(jié)果如表3所示，考慮到計(jì)算結(jié)果的用途，本文只計(jì)算到樁端持力層，且計(jì)算結(jié)果為各土層側(cè)阻的相關(guān)范圍，對(duì)于端阻的相關(guān)范圍，只考慮了樁端持力層。由表3可見(jiàn):兩種方法計(jì)算出的絕大多數(shù)土層的相關(guān)范圍很接近，說(shuō)明在對(duì)群樁基礎(chǔ)的可靠性分析時(shí)，利用相關(guān)函數(shù)法和變異函數(shù)法均具有一定的可靠性，計(jì)算結(jié)果精度較高;少數(shù)結(jié)果有一定差距，為了降低計(jì)算結(jié)果的誤差，推薦采用兩種算法的平均結(jié)果。其余各孔的計(jì)算結(jié)果如表4所示，該結(jié)果為兩種計(jì)算方法的平均值。

表3 J-5孔土層的相關(guān)范圍

表4 各孔土層的相關(guān)范圍

4 結(jié)論

(1)在鉆孔灌注群樁基礎(chǔ)的可靠度分析中，隨機(jī)場(chǎng)理論是考慮土工參數(shù)空間變異性的重要方法。

(2)采用相關(guān)函數(shù)法和變異函數(shù)法對(duì)群樁基礎(chǔ)中的相關(guān)范圍進(jìn)行研究，通過(guò)比較發(fā)現(xiàn)兩種算法都具有較高可靠性。

(3)采用相關(guān)函數(shù)法和變異函數(shù)法計(jì)算結(jié)果的平均值，作為各土層的相關(guān)范圍，可進(jìn)一步減小如人為因素等原因引起的誤差。

［1］Vanmarcke E H.Probabilistic Modeling of Soil Profiles［J］.Journal of the Geotechnical Engineering Division，ASCE，1977，103(11):1227－1246.

［2］傅旭東.土工參數(shù)相關(guān)范圍及相關(guān)距離的計(jì)算方法［J］.西南交通大學(xué)學(xué)報(bào)，1996，31(5):510－515.

［3］張梅，丁繼輝，宋向東.對(duì)求解土性相關(guān)距離的空間遞推平均法的分析與改進(jìn)［J］.河北農(nóng)業(yè)大學(xué)學(xué)報(bào)，1999，22(1):85－89.

［4］丁繼輝，宇云飛，張慶宏，等.土性相關(guān)距離計(jì)算的非線性曲線擬合法［J］.勘察科學(xué)技術(shù)，2000(3):13－16.

［5］程強(qiáng)，羅書(shū)學(xué)，高新強(qiáng).相關(guān)函數(shù)法計(jì)算相關(guān)距離的分析探討［J］.巖石力學(xué)，2000，21(3):281－283.

［6］劉春原.基于GIS系統(tǒng)的巖土參數(shù)隨機(jī)場(chǎng)特性研究［D］.天津:天津大學(xué)，2003.

［7］姚若軍.土性參數(shù)空間變異和相關(guān)情況下的土坡可靠度分析［D］.武漢:武漢大學(xué)，2005.

［8］郭懷志，彭大鵬.材料性能的隨機(jī)場(chǎng)特性參數(shù)的檢定方法［J］.巖土工程學(xué)報(bào)，1994，16(3):21－29.

［9］宋向東，任建華.土工參數(shù)隨機(jī)場(chǎng)模型中相關(guān)函數(shù)的選擇［J］.數(shù)理統(tǒng)計(jì)與管理，2007，26(5):773－775.

［10］尹利華，王曉謀，張留俊.天津軟土土性指標(biāo)概率分布統(tǒng)計(jì)分析［J］.巖石力學(xué)，2010，31(2):462－468.