康晨陽，劉益平，葛海明

(江蘇省電力設(shè)計(jì)院有限公司，江蘇 南京 211102)

1 概述

工程勘察中，土樣從地下取出要經(jīng)過鉆探取樣、運(yùn)輸、貯藏、試樣制備等步驟，然后通過儀器測量其物理力學(xué)性質(zhì)參數(shù)。顯然，每一環(huán)節(jié)對(duì)土樣都有不可避免的影響，而對(duì)于靈敏度高，結(jié)構(gòu)性強(qiáng)的軟土，為更加經(jīng)濟(jì)、安全反映其原位土的特性，研究擾動(dòng)狀況對(duì)了解和認(rèn)識(shí)軟土的工程特性具有重要意義。

國內(nèi)外學(xué)者研究了取樣后土樣擾動(dòng)度的定義和評(píng)價(jià)方法，概括起來有殘余有效應(yīng)力法、殘余孔隙水壓力法、不排水變形模量對(duì)比法、固結(jié)屈服應(yīng)力法、體積壓縮法等。其中，Z. Hong 和K.Onitsuka在雙對(duì)數(shù)坐標(biāo)中的ln(1+e)－lg p 土體壓縮曲線呈雙直線特性的基礎(chǔ)上，提出了修正的體積壓縮法，該方法得到廣泛應(yīng)用。

傳統(tǒng)的修正體積壓縮法通過作圖實(shí)現(xiàn)擾動(dòng)度的評(píng)價(jià)，對(duì)于壓縮曲線斜率和屈服點(diǎn)的確定具有一定的經(jīng)驗(yàn)性和偶然性。本文基于修正體積壓縮法的原理，采用門限分位數(shù)自回歸的數(shù)學(xué)模型，自動(dòng)捕捉屈服點(diǎn)，并進(jìn)行參數(shù)檢驗(yàn)，通過統(tǒng)計(jì)分析軟件R語言進(jìn)行編程計(jì)算，對(duì)土樣進(jìn)行擾動(dòng)度評(píng)價(jià)，分析擾動(dòng)度和屈服應(yīng)力及十字板抗剪強(qiáng)度之間的關(guān)系，進(jìn)而討論了取樣擾動(dòng)的恢復(fù)處理方法。

2 擾動(dòng)度的評(píng)價(jià)方法

R.Butter fi eld通過對(duì)大量的固結(jié)試驗(yàn)數(shù)據(jù)分析，提出在雙對(duì)數(shù)坐標(biāo)ln(1+e)－lgp下，土體的壓縮曲線在屈服前階段(p＜p'y)和屈服后階段(p＞p'y)均為兩條直線。取樣擾動(dòng)影響這兩條直線的斜率。

Z. Hong 和K. Onitsuka根據(jù)Butter－ fi eld的研究結(jié)果，改進(jìn)了傳統(tǒng)的體積壓縮法，得到擾動(dòng)度的定義如下：

式中：CCLB，CCLR分別為擾動(dòng)樣和重塑樣屈服前l(fā)n(1+e)－lgp坐標(biāo)中壓縮曲線的斜率(見圖1)。

圖1 擾動(dòng)度的定義

如圖1所示，CCLA為擾動(dòng)樣屈服后壓縮曲線的斜率。定義p'0為土體上覆有效應(yīng)力，p'y為屈服應(yīng)力。當(dāng)SD=0%時(shí)，土樣未擾動(dòng)；當(dāng)SD=100%時(shí)土樣完全擾動(dòng)；當(dāng)SD為0%～100%，壓縮曲線在原位與完全重塑土的壓縮曲線之間。

Z. Hong 和K. Onitsuka通過研究發(fā)現(xiàn)，CCLR為液限的函數(shù)，即

3 門限分位數(shù)自回歸模型

由式(1)可以看出，ln(1+e)－lgp坐標(biāo)中壓縮曲線的斜率直接影響擾動(dòng)度的評(píng)價(jià)。在實(shí)際操作中，考慮到室內(nèi)土工試驗(yàn)結(jié)果的離散性，通過作圖的方法獲得曲線斜率和屈服點(diǎn)，具有一定的經(jīng)驗(yàn)性和偶然性，進(jìn)而影響擾動(dòng)度的評(píng)價(jià)結(jié)果。為更好的實(shí)現(xiàn)雙對(duì)數(shù)壓縮曲線的兩階段線性擬合，自動(dòng)捕捉屈服點(diǎn)，并檢驗(yàn)擬合誤差，引入金融經(jīng)濟(jì)系統(tǒng)預(yù)測中常用的非線性門限自回歸模型。

門限自回歸模型(TAR)最早由經(jīng)濟(jì)學(xué)家Tong(1978，1983)提出，針對(duì)金融經(jīng)濟(jì)系統(tǒng)預(yù)測中，內(nèi)部變量具有的非線性和異質(zhì)性等特征，通過門限變量來控制分段線性機(jī)制，進(jìn)而能夠較好地預(yù)測金融時(shí)間序列處于不同階段時(shí)，所呈現(xiàn)的非對(duì)稱特征。

本文利用兩階段門限自回歸模型(TQAR)，建模過程包括模型表示、參數(shù)估計(jì)、模型定階及診斷檢驗(yàn)。

3.1 模型表示

假設(shè)時(shí)間序列{yt}為一維響應(yīng)變量，{xt=(1，yt-1，yt-2，…yt-p，)}為p+1維向量組成的解釋變量，{qt}為門限變量。數(shù)學(xué)模型如下：

3.2 參數(shù)估計(jì)

TQAR模型的參數(shù)估計(jì)可通過優(yōu)化下式得到：

式中：T為樣本量；p為滯后期；St(θ(τ)，γ(τ))為目標(biāo)函數(shù)；ρτ(u)為非對(duì)稱損失函數(shù)，滿足

式(3)的求解過程可通過兩步法實(shí)現(xiàn)。

3.3 模型定階

通過AIC準(zhǔn)則來確定最優(yōu)滯后階數(shù)p，定義如下：

最優(yōu)滯后階數(shù)p*的選擇標(biāo)準(zhǔn)為：使用盡可能簡單的結(jié)構(gòu)，最大程度地?cái)M合樣本數(shù)據(jù)，即一般通過網(wǎng)格搜索方法，獲得p*的取值。

3.4 診斷檢驗(yàn)

TQAR模型的診斷檢驗(yàn)，主要包括兩方面，第一，門限效應(yīng)存在性檢驗(yàn)；第二，回歸方程及回歸系數(shù)顯著性檢驗(yàn)。本文采用Koenker等構(gòu)造的似然比檢驗(yàn)，通過標(biāo)準(zhǔn)差來反映模型中參數(shù)的精度。

4 工程實(shí)例

4.1 工程背景

某工程位于連云港徐圩開發(fā)區(qū)，地貌單元為海積平原。地層主要由第四系全新統(tǒng)海相沉積成因的黏土、淤泥、粉質(zhì)黏土、粉質(zhì)黏土夾粉土和粉砂等組成。場地20 m以內(nèi)各土層的物理力學(xué)參數(shù)見表1。

表1 土層物理力學(xué)參數(shù)

4.2 擾動(dòng)度評(píng)價(jià)

層②淤泥為連云港典型深厚軟土，現(xiàn)場鉆探通過靜壓敞口式薄壁取樣器獲得該層不同深度處原狀土樣共17件，進(jìn)行室內(nèi)高壓固結(jié)試驗(yàn)，最大壓力為1600 kPa。

利用統(tǒng)計(jì)分析軟件R語言分析計(jì)算，實(shí)現(xiàn)門限自回歸模型的數(shù)值擬合，以1號(hào)土樣為例，擬合曲線見圖1。

圖2 數(shù)值擬合結(jié)果

根據(jù)門限自回歸模型的數(shù)值擬合結(jié)果，試樣擾動(dòng)度SD、不同擾動(dòng)程度時(shí)土體的屈服應(yīng)力見表2。

表2 土樣擾動(dòng)度評(píng)價(jià)結(jié)果

由以上分析結(jié)果可以看出，采用門限自回歸模型計(jì)算所得評(píng)價(jià)指標(biāo)的標(biāo)準(zhǔn)差均小于1%，精度高，數(shù)據(jù)擬合結(jié)果較好。

根據(jù)修正的體積壓縮法的評(píng)價(jià)結(jié)果，通過薄壁取土器獲得的軟土原狀試樣擾動(dòng)度SD在13.1%～31.4%之間，平均擾動(dòng)度為22.0%，說明連云港軟土具有較強(qiáng)的結(jié)構(gòu)性，經(jīng)取樣、運(yùn)輸、制樣過程后，土樣均受到了不同程度的擾動(dòng)。

4.3 擾動(dòng)與屈服應(yīng)力的關(guān)系

圖3 擾動(dòng)度與屈服應(yīng)力關(guān)系

鄧永峰在總結(jié)M. Nagaraj，S. G. Chung，S. Leroueil等人的研究成果上，建立了不同擾動(dòng)度時(shí)土體的屈服應(yīng)力和擾動(dòng)度SD之間的關(guān)系如下式：

圖3表明：土樣擾動(dòng)度與屈服應(yīng)力的關(guān)系與鄧永峰的研究結(jié)果基本一致；屈服應(yīng)力p'y和擾動(dòng)度SD之間呈負(fù)冪指數(shù)關(guān)系，土體受擾動(dòng)程度越大，土體的屈服應(yīng)力衰減的程度越大。

4.4 擾動(dòng)與原位強(qiáng)度參數(shù)的關(guān)系

根據(jù) S. Leroueil 等的研究結(jié)果，軟土實(shí)測的十字板抗剪強(qiáng)度cu與前期屈服應(yīng)力比值為0.2～0.4，且為塑性指數(shù)的函數(shù)，關(guān)系式為：

圖4 十字板抗剪強(qiáng)度和屈服應(yīng)力隨深度變化

如圖5建立擾動(dòng)度SD與的關(guān)系，可見，軟土的十字板抗剪強(qiáng)度cu與屈服應(yīng)力的比值約為0.25～0.55，較S. Leroueil研究結(jié)果偏大，這是由于土樣擾動(dòng)后屈服應(yīng)力有所折減。

圖5 擾動(dòng)度與十字板抗剪強(qiáng)度及屈服應(yīng)力關(guān)系

4.5 取樣擾動(dòng)的恢復(fù)方法

采用K0固結(jié)預(yù)處理法對(duì)比試樣處理前后的無側(cè)限抗壓強(qiáng)度，并同實(shí)測的十字板抗剪強(qiáng)度比較分析，見圖6。

圖6 不同試驗(yàn)方法的抗剪強(qiáng)度

從圖6看出，未預(yù)處理時(shí)試樣強(qiáng)度僅為天然強(qiáng)度的47.3%，固結(jié)預(yù)處理后試樣強(qiáng)度達(dá)天然強(qiáng)度的79.2%。說明經(jīng)K0固結(jié)預(yù)處理，試樣強(qiáng)度恢復(fù)明顯，試驗(yàn)成果更接近天然強(qiáng)度，但仍有20.1%未恢復(fù)。因此，固結(jié)預(yù)處理措施雖然能消除部分取樣擾動(dòng)的影響，但不能全部恢復(fù)，抗剪強(qiáng)度指標(biāo)的選擇須結(jié)合室內(nèi)試驗(yàn)和原位測試綜合確定。

5 結(jié)論

(1)根據(jù)修正的體積壓縮法的原理，采用門限自回歸模型自動(dòng)捕捉屈服點(diǎn)，評(píng)價(jià)指標(biāo)的標(biāo)準(zhǔn)差均小于1%，精度高，很好的反映雙對(duì)數(shù)曲線的兩階段性，數(shù)值擬合結(jié)果較好。

(2)通過敞口式薄壁取土器獲得的連云港地區(qū)軟土試樣平均擾動(dòng)程度為22.0%，說明連云港軟土具有較強(qiáng)結(jié)構(gòu)性，經(jīng)取樣、運(yùn)輸、制樣過程后土樣均受到了不同程度的擾動(dòng)。

(3)建立土樣擾動(dòng)度與屈服應(yīng)力的關(guān)系，屈服應(yīng)力和擾動(dòng)度之間呈負(fù)冪指數(shù)關(guān)系，土體受擾動(dòng)程度越大，土體的屈服應(yīng)力衰減的程度越大，這與鄧永峰的研究結(jié)果基本一致。

(4)根據(jù)理論公式推算得到的十字板剪切強(qiáng)度較實(shí)測值平均減少37%，且當(dāng)深度越大時(shí)，實(shí)測值與理論值偏差越大，進(jìn)而說明通過“原狀樣”獲得的室內(nèi)強(qiáng)度指標(biāo)有一定程度的失真。

(5)通過K0固結(jié)預(yù)處理的方法，能很大程度上消除取樣擾動(dòng)的影響，試樣抗剪強(qiáng)度大幅提高，但不能完全恢復(fù)天然強(qiáng)度，抗剪強(qiáng)度指標(biāo)的選擇須結(jié)合室內(nèi)試驗(yàn)和原位測試綜合確定。

[1] 胡中雄．土力學(xué)與環(huán)境土工學(xué)[M]．上海： 同濟(jì)大學(xué)出版社，1997．

[2] Ladd C C，Lambe T W．The strength of undisturbed clay determined from undrained test[J] ASTMSTP，1963，(85)．

[3] Raymond G P．The effect of sampling on the undrained soil properties of a Leda clay[J]．Canadian Geotechnical Journal，1971，(8)．

[4] Nagaraj M，Chung S G．Analysis and assessment of sampling disturbance of soft sensitive clays[J]．Geotechnique，2003，53(7)．

[5] Hong Z，Onitsuka K．A method of correcting yield stress and compression index of Ariake clays for sample disturbance[J]．Soils and Foundations，1998，38(2)．

[6] 劉海濤，郭院成．軟土結(jié)構(gòu)性分析與擾動(dòng)分析[J]．鄭州大學(xué)學(xué)報(bào)(工學(xué)版)，2014．38(5)．

[7] BUTTERFIELD R．A natural compression law for soils(an advanceon e－lgp) [J]．Geotechnique，1979，29(4)．

[8] Tong H．Threshold models in non～linear time series analysis [M]．Berlin：Springer－Verlag，1983．

[9] 康寧，荊科．門限分位數(shù)自回歸模型的預(yù)測方法及應(yīng)用[J]．?dāng)?shù)量經(jīng)濟(jì)技術(shù)經(jīng)濟(jì)研究，2016，(3)．

[10] Koenker R．，Machado J．A．F．，Goodness of fi t and related inference processes for quantile regression [J]．Journal of the American statistical association，1999，94(448)．

[11] 鄧永峰，劉松玉．?dāng)_動(dòng)對(duì)軟土強(qiáng)度影響規(guī)律研究[J]．巖石力學(xué)與工程學(xué)報(bào)，2007，26(9)．

[12] LEROUEIL S，MAGNAN J P，TAVENAS F．Embankments on soft clays[M]．England： Eliis Horwood Limited，1990．

[13] 劉華清，等．取樣擾動(dòng)對(duì)土天然強(qiáng)度指標(biāo)的影響和處理方法[J]．巖土工程技術(shù)，2002， (3)．