基于非平穩(wěn)隨機場的水泥土堤壩邊坡可靠度分析

楊亮 劉輝 章榮軍 鄭俊杰

摘 要：采用水泥固化疏浚淤泥作為堤壩填料往往會使堤壩土體參數(shù)具有較高的變異性;同時，在水泥土堤壩填筑過程中，通常設(shè)置水平排水板加速土體排水固結(jié)，這使得主固結(jié)完成后水泥土強度隨排水板位置呈現(xiàn)一定空間變化趨勢，通常的平穩(wěn)或準(zhǔn)平穩(wěn)隨機場模型很難表征這一趨勢。為此，提出一種不排水抗剪強度垂直排水板方向衰減的非平穩(wěn)隨機場模型，并以水泥土堤壩邊坡為例，結(jié)合強度折減法和蒙特卡洛模擬，對比提出的非平穩(wěn)隨機場模型和傳統(tǒng)平穩(wěn)隨機場模型、深度趨勢非平穩(wěn)隨機場模型對邊坡安全系數(shù)和滑動面分布的影響，討論模型中趨勢參數(shù)對可靠度的影響。結(jié)果表明：平穩(wěn)隨機場模型及深度趨勢非平穩(wěn)隨機場模型的計算結(jié)果偏于危險，提出的非平穩(wěn)隨機場模型既能考慮土體不排水抗剪強度沿深度方向的趨勢項與隨機波動項，又能考慮其沿垂直排水板方向的衰減趨勢，與現(xiàn)有研究符合較好。

關(guān)鍵詞：非平穩(wěn)隨機場;水泥土邊坡;可靠度分析;不排水抗剪強度;堤壩

中圖分類號：TU457 文獻標(biāo)志碼：A 文章編號：2096-6717（2022）05-0078-09

收稿日期：2021-07-05

基金項目：國家自然科學(xué)基金（51978303）

作者簡介：楊亮（1998- ），男，主要從事巖土工程研究，E-mail：yangliang0815@hust.edu.cn。

章榮軍（通信作者），男，教授，博士生導(dǎo)師，E-mail：ce_zhangrj@hust.edu.cn。

Received：2021-07-05

Foundation items：National Natural Science Foundation of China （No. 51978303）

Author brief：YANG Liang （1998- ）， main research interests： geotechnical engineering， E-mail： yangliang0815@hust.edu.cn.

ZHANG Rongjun （corresponding author）， professor， doctorial supervisor， E-mail： ce_zhangrj@hust.edu.cn.

Reliability analysis of embankment slope in cement-mixed soil with non-stationary random field

YANG Liang， LIU Hui， ZHANG Rongjun， ZHENG Junjie

（Institute of Geotechnical and Underground Engineering， Huazhong University of Science and Technology， Wuhan 400074， P. R. China）

Abstract：The soil properties of embankment slope filled with cement solidified dredged sludge commonly feature high variability. Meanwhile， in construction process of cement-mixed embankment， horizontal drainage plate is generally used to accelerate the consolidation of cement-mixed soil， which makes the soil strength show a spatial distribution trend at the vertical direction （perpendicular to the drainage plate） after the primary consolidation finished. The traditional stationary or quasi-stationary random field model cannot well characterize this trend. This paper proposes a non-stationary model to characterize the attenuation trend of soil strength along the normal direction of drainage plate. The strength reduction method and Monte Carlo simulation are combined to evaluate the reliability of slope， and the proposed model is compared with two existing models （i.e.， traditional stationary random field model and depth trend non-stationary random model）. The influence of different random filed models on safety factor and slip surface are investigated， and the impact of trend parameters in the proposed model on slope reliability are discussed. The results show that the stationary random field model and depth trend non-stationary random field model may overestimate the reliability of embankment slope， and the proposed non-stationary random field model can consider both the trend term and random fluctuation term in vertical direction and the attenuation trend in normal direction of drainage plate， which is in good agreement with the existing research.

Keywords：non-stationary random field; cement-mixed soil slope; reliability analysis; undrained shear strength; embankment

為了處置大量的河湖疏浚淤泥，一種理想的做法是采用水泥固化高含水率泥漿，再結(jié)合真空預(yù)壓、堆載預(yù)壓等多種方式將水泥土作為工程（如堤壩）填料，但由于疏浚淤泥雜質(zhì)眾多、施工過程中攪拌不勻、養(yǎng)護條件不一等多種因素，疏浚淤泥固化處理后形成的水泥土堤壩邊坡材料參數(shù)分布存在較大變異性。

在固化淤泥堤壩填筑過程中，常采用真空預(yù)壓加速土體固結(jié)，需要在堤壩土體內(nèi)鋪設(shè)多層水平排水板，針對這種條件下堤壩土體內(nèi)抗剪強度的分布規(guī)律，學(xué)者們已開展相關(guān)試驗研究。周洋等通過試驗發(fā)現(xiàn)，僅鋪設(shè)水平排水板時，垂直于水平排水板方向，距離越遠(yuǎn)土體含水率越高，并且從土體含水率分布情況來看，僅鋪設(shè)水平排水板的效果要優(yōu)于僅鋪設(shè)豎向排水板的效果;張文彬等發(fā)現(xiàn)，隨著離排水板距離的增加，土體的十字板抗剪強度逐漸減小;高會強發(fā)現(xiàn)，鋪設(shè)水平排水板時，加固后土體表層和靠近排水板的土體強度較高，而遠(yuǎn)離排水板的土體強度很低，鋪設(shè)豎向排水板時，隨著深度的增加，十字板抗剪強度呈線性增加趨勢。根據(jù)已有的試驗結(jié)果可知，在堤壩內(nèi)鋪設(shè)水平排水板會使土體的強度特性分布受到影響，且一般的規(guī)律是：在一定范圍內(nèi)，距離排水板越遠(yuǎn)，含水率越高，土體強度越低。

為了描述土體參數(shù)的空間分布變異性并進行可靠度分析，常采用平穩(wěn)隨機場模型。所謂平穩(wěn)隨機場即認(rèn)為土體參數(shù)的均值、標(biāo)準(zhǔn)差等統(tǒng)計特征沿埋深保持不變，參數(shù)之間的空間相關(guān)性只與兩點之間的相對距離相關(guān)而與其絕對距離無關(guān)。這種模型很難表征排水板固結(jié)帶來的土體強度空間變化趨勢。相對于平穩(wěn)隨機場，也有學(xué)者提出了非平穩(wěn)隨機場分析模型，在非平穩(wěn)隨機場的表征過程中，土體參數(shù)離散變量一般可視為趨勢項、殘余項與波動分量的多項式，以此表征土體參數(shù)在水平方向上的波動性和埋深方向上的趨勢性，許多學(xué)者展開了相關(guān)研究，提出了非平穩(wěn)隨機場的表征方法與趨勢項及殘余項的處理方法。蔣水華等綜合分析了兩種非平穩(wěn)隨機場模型的特征：1）采取去趨勢分析方法，將土體參數(shù)模擬為某一平穩(wěn)隨機場，再基于該平穩(wěn)隨機場，進一步考慮土體參數(shù)隨埋深的變化趨勢分量的影響，從而實現(xiàn)對土體參數(shù)二維非平穩(wěn)隨機場的表征，該方法所得非平穩(wěn)隨機場的均值與標(biāo)準(zhǔn)差均隨埋深而變化，但其變異系數(shù)保持不變;2）取趨勢項為趨勢分量與土體重度和埋深的乘積，殘余項為地表處土體參數(shù)值，繼而將趨勢分量參數(shù)模擬為服從某一分布的平穩(wěn)隨機場，忽略殘余項的不確定性，得到二維非平穩(wěn)隨機場，該方法所得土體參數(shù)的均值、標(biāo)準(zhǔn)差與變異系數(shù)均隨埋深而變化;在此基礎(chǔ)上，Li等提出了一種新的二維非平穩(wěn)隨機場表征方法，驗證了所提方法的可行性，系統(tǒng)比較了所提方法與現(xiàn)有二維非平穩(wěn)隨機場與隨機場的區(qū)別，該方法的優(yōu)勢在于能夠同時考慮趨勢分量與波動分量的不確定性。豆紅強等通過去趨勢分析方法，將飽和滲透系數(shù)均值視為隨埋深遞減的線性函數(shù)，殘余項為地表處平均飽和滲透系數(shù)，趨勢分量為飽和滲透系數(shù)隨埋深的變化梯度，建立了表征土體飽和滲透系數(shù)變異性的一維非平穩(wěn)隨機場模型，并探討了土體飽和滲透系數(shù)的空間變異性對邊坡可靠度的影響。Li等通過不排水抗剪強度、有效豎向應(yīng)力和黏土超固結(jié)比（OCR）之間的經(jīng)驗關(guān)系，研究了不排水剪切強度隨埋深的變化趨勢，在此基礎(chǔ)上得到了趨勢項隨埋深變化的二維非平穩(wěn)隨機場。Shu等取殘余項為泥水分界線處的不排水剪切強度，取趨勢項為埋深的線性函數(shù)，且趨勢項與殘余項滿足對數(shù)正態(tài)分布，建立了不排水剪切強度的二維非平穩(wěn)隨機場，討論了自相關(guān)距離對深?；A(chǔ)平均承載力與破壞機理的影響;Wu等取趨勢項為地表處非零且隨埋深線性變化的不排水剪切強度均值，保持變異系數(shù)不變，通過高斯平穩(wěn)隨機場的轉(zhuǎn)換與映射，得到了不排水剪切強度的二維非平穩(wěn)隨機場表征，研究了變異系數(shù)、分布函數(shù)、垂直自相關(guān)距離等因素對土體承載力的影響。

基于上述垂直排水板方向強度衰減的一般規(guī)律及非平穩(wěn)隨機場相關(guān)理論，筆者提出一種水泥土堤壩內(nèi)垂直排水板方向的不排水抗剪強度隨距離衰減的非平穩(wěn)隨機場模型，該模型既能表征土體抗剪強度沿埋深方向增加的趨勢，又能考慮垂直排水板方向的強度衰減。在此基礎(chǔ)上，給出土體參數(shù)三維非平穩(wěn)隨機場模擬計算方法流程，并分析了土體參數(shù)非平穩(wěn)分布特征對邊坡可靠度的影響。

1 不排水抗剪強度非平穩(wěn)隨機場模型

1.1 深度趨勢隨機場

土體參數(shù)的空間變異性一般由趨勢項和隨機波動項聯(lián)合表征。以無側(cè)限抗壓強度隨機場為例，試驗值q沿深度方向的不平穩(wěn)特征可表示為

q（h）=t（h）+ω（1）

式中：h為土體埋深;t（h）為與埋深有關(guān)的趨勢項，一般取相應(yīng)深度下土體參數(shù)的均值;ω為隨機波動項，用來描述土體空間內(nèi)土體參數(shù)大小無趨勢的隨機波動特征，其均值和標(biāo)準(zhǔn)差與埋深無關(guān)。通常，土體的無側(cè)限抗壓強度隨深度增加而增加。簡單起見，用線性函數(shù)描述土體無側(cè)限抗壓強度隨深度的變化關(guān)系，即

q=q+b·σ′=q+γ·h·b（2）

式中：q為地表土體的無側(cè)限抗壓強度（地表UCT試驗值）;b為無側(cè)限抗壓強度隨埋深增加的速率;σ′=γh，為豎向有效應(yīng)力，γ為土體重度。參考文獻[10]的做法，采用去趨勢分析方法，首先用一均值為μ、標(biāo)準(zhǔn)差為σ的對數(shù)正態(tài)平穩(wěn)隨機場模擬q，再在此基礎(chǔ)上疊加一個沿埋深線性增加的趨勢分量，最終得到三維深度趨勢非平穩(wěn)隨機場

q=q（x，y，z）+γ·h·b（3）

1.2 垂直排水板不排水抗剪強度衰減非平穩(wěn)隨機場模型 基于已有的試驗研究，可假設(shè)土體在排水固結(jié)后，抗剪強度與其離排水板的距離線性相關(guān)，此時土體隨機場可表示為

q（x，y，z）

=q+γ·h·b-d（x，z）·b·exp[ω（x，y，z）]（4）

式中：q、b的含義同前所述，只是此處q為均值，為μ、標(biāo)準(zhǔn)差為σ的隨機變量;d（x，z）為土體空間內(nèi)一點與排水板的垂直距離（排水板沿y方向鋪設(shè)）;b為垂直排水板方向的趨勢項;exp[ω（x，y，z）]為隨機波動項。

如圖1所示，P為三維空間中截取的某一x-z平面內(nèi)一點，因排水板（鋪設(shè)）方向平行于y軸方向，故點P垂直排水板距離d（指距離最近的排水板）只與點P的x、z坐標(biāo)有關(guān)，而與x-z平面所處的位置（y坐標(biāo)）無關(guān)，表示為d（x，z）。

2 算例分析

2.1 算例描述

算例為一高陡水泥土堤壩邊坡，斜坡高度6 m，坡度1∶0.5，堤壩頂面寬度為18 m，堤壩頂面作用20 kPa的交通荷載（假設(shè)為均勻分布）。為了提高計算效率，取半結(jié)構(gòu)進行建模分析，模型斷面圖和三維網(wǎng)格圖如圖2所示。

2.2 材料參數(shù)取值

參照Zhang等、Lee等的研究，認(rèn)為水泥土堤壩邊坡無側(cè)限抗壓強度服從對數(shù)正態(tài)分布，固化淤泥形成的水泥土模量E與無側(cè)限抗壓強度q呈正比關(guān)系，E/q=135～140?？紤]到固化淤泥形成的水泥土材料滲透系數(shù)很低，參考Liu等的做法，假定主固結(jié)完成后水泥土堤壩為完全不排水條件，采用內(nèi)摩擦角為0的Mohr-Coulomb模型來模擬固化淤泥形成的水泥土，取不排水抗剪強度c=q/2、E/q=140。相應(yīng)地，堤壩地基土體為老黏土，相關(guān)材料參數(shù)見表1。

為了探究垂直排水板衰減非平穩(wěn)隨機場模型的可行性，采用3種隨機場模型模擬水泥土堤壩部分對比計算結(jié)果，下部老黏土采用一般的平穩(wěn)隨機場模型。

模型1為平穩(wěn)隨機場模型，采用地表處的UCT試驗均值μ=86 kPa，不考慮其趨勢變化，UCT試驗值空間隨機場采用對數(shù)正態(tài)分布，變異系數(shù)COV=0.3。采用高斯型自相關(guān)函數(shù)，三維空間內(nèi)水平、豎向相關(guān)距離分別取為38、3.8 m。

模型2為深度趨勢隨機場模型，三維空間隨機場生成公式為

q（x，y，z）=q（x，y，z）+γ·（6-z）·b（5）

式中：q（x，y，z）為對數(shù)正態(tài)分布的平穩(wěn)隨機場，離散取值同上。土體趨勢分量參數(shù)b隨深度變化，參照蔣水華的統(tǒng)計結(jié)果，取b=0.2，為定值。不排水抗剪強度的均值和標(biāo)準(zhǔn)差分別為

μ（h）=μ+γhb

σ（h）=σ（6）

模型3為垂直排水板衰減非平穩(wěn)隨機場模型，在模型計算時，認(rèn)為排水板附近土體的抗剪強度不斷減小。同時，采用靜力計算，不考慮排水帶來的滲透固結(jié)作用，也不考慮排水板自身對土體強度分布的影響。此時，三維空間隨機場生成模型為

q（x，y，z）=q+γ·（6-z）·b-

d（x，z）·b·exp[w（x，y，z）]（7）

式中：q為對數(shù)正態(tài)分布的隨機變量，μ=86 kPa，變異系數(shù)COV=0.3;趨勢參數(shù)b=0.2，為定值;b為因真空預(yù)壓排水固結(jié)導(dǎo)致的q隨距排水板距離衰減趨勢項（排水板處最大），理論上與施工條件有關(guān)，參考b的取值結(jié)果，考慮其隨機性，取b為對數(shù)正態(tài)分布的隨機變量，均值μ=2.0，變異系數(shù)COV=0.3;exp[w（x，y，z）]為隨機波動項，w（x，y，z）取為均值μ=0、標(biāo)準(zhǔn)差σ=0.3的正態(tài)分布隨機場，水平、豎向相關(guān)距離分別取38、3.8 m。不排水抗剪強度的均值和標(biāo)準(zhǔn)差分別為

μ（h）=μ+γhb+dμ

σ（h）=σ+dexp（σ）

{（μ+σ）[exp（σ）-1]+σ}（8）

2.3 水泥土堤壩邊坡可靠度計算分析流程

采用Karhunen-Loeve級數(shù)展開法離散隨機場，每一工況生成500次隨機場，結(jié)合蒙特卡洛模擬和強度折減法分析堤壩邊坡失效概率，主要分析計算流程如圖3所示，其中，強度折減法將邊坡安全系數(shù)定義為使邊坡剛好達到臨界破壞狀態(tài)時對其強度的折減程度。采用Mohr-Coulomb破壞準(zhǔn)則時，其折減的強度參數(shù)是黏聚力c和內(nèi)摩擦角φ（對于本例的不排水條件，僅折減黏聚力c），一般通過二分法迭代計算，通過不斷擴大（或縮?。┱蹨p值，直至折減系數(shù)的上下限滿足精度要求，將安全系數(shù)F定義為初始黏聚力c和破壞時的黏聚力c之比，即F=c/c。

隨后統(tǒng)計一個工況內(nèi)500組隨機場模型計算結(jié)果，可得到安全系數(shù)的分布，并借此計算可靠指標(biāo)、失效概率。

若安全系數(shù)的分布呈正態(tài)分布，則可靠指標(biāo)

β=μ-1σ（9）

當(dāng)安全系數(shù)呈對數(shù)正態(tài)分布時

β=lnμ1+Vln（1+V）（10）

式中：β為可靠指標(biāo);μ為500組安全系數(shù)均值;σ為500組安全系數(shù)標(biāo)準(zhǔn)差，V=σ /μ。利用可靠指標(biāo)β與失效概率P之間的數(shù)學(xué)關(guān)系，可以用式（11）計算出可靠指標(biāo)。

P=1-Φ（β）（11）

3 3種隨機場模型對比

3.1 3種隨機場模型實現(xiàn)情況對比

圖4為3個模型堤壩土部分同一豎線上（參考圖5（a）、（b），y=0平面上、x=3.6 m這條線）不同c的統(tǒng)計值（c=q/2）?？梢钥闯觯瑢τ谀Ｐ?，c大小隨深度并無明顯增加趨勢，但沿深度呈現(xiàn)出隨機波動趨勢，且波動幅度較大。模型2呈現(xiàn)出c均值隨深度增加的趨勢，能夠在一定程度上表征土體的非平穩(wěn)特征。

模型3則表現(xiàn)出了沿深度增加和周期性波動的趨勢，波動幅度明顯小于模型1和模型2。這種周期性波動的原因是：土體不排水抗剪強度c沿深度方向有增加的趨勢，但由于排水板加速局部的排水固結(jié)，使得局部強度最高，局部強度最大值即為排水板鋪設(shè)位置，遠(yuǎn)離排水板則強度減小。這一效果與前述試驗規(guī)律基本一致，且隨著深度的增加，波動幅度增大，說明模型3的c均值沿深度方向增加。

圖5為模型3三維網(wǎng)格不排水抗剪強度c分布結(jié)果，采用局部平均方法得到，圖5（b）、（c）分別為三維模型的兩方向上的斷面，排水板鋪設(shè)方向平行于邊坡橫截面，其橫截面分布在三維邊坡內(nèi)的鋪設(shè)點為圖5（b）中的白點，可以看出，排水板四周土體的強度最高，隨著離排水板距離的增加，土體強度逐漸降低，也基本符合前文所述研究規(guī)律。在平行于排水板方向，邊坡橫截面上，不排水抗剪強度c分布則近似于深度趨勢隨機場，在深度方向上呈增大的趨勢，波動性亦很強。

3.2 3種隨機場模型計算結(jié)果對比

對上述3種典型工況分別進行500組蒙特卡洛模擬，利用強度折減法得出500組安全系數(shù)，模型3的安全系數(shù)分布情況如圖6所示，對分布情況進行K-S檢驗發(fā)現(xiàn)，采用對數(shù)正態(tài)分布擬合時，其p=1，說明對數(shù)正態(tài)分布符合較好。經(jīng)檢驗，其他工況安全系數(shù)分布同樣符合對數(shù)正態(tài)分布。

模型1、模型2、模型3分別對應(yīng)的可靠指標(biāo)為1.34、2.79、1.30。從式（5）可見，模型2深度趨勢模型的各項均為正，強度不存在衰減，理論上是“最安全”的工況，結(jié)果也驗證了這一點。

參考李劍等的做法，選取邊坡達到極限狀態(tài)時所有垂線上最大剪應(yīng)變增量的位置，通過最小二乘法擬合圓弧滑動面，3種典型工況的滑動面分布結(jié)果如圖7所示。3種典型工況的邊坡滑動面都主要發(fā)生在水泥土堤壩土體內(nèi)，但模型1、模型2的離散性更強，滑動面頂端與坡頂距離不一，滑動體體積也差距較大，模型3的滑動面則主要為從坡底貫穿坡頂?shù)幕瑒用?，分布集中?/p>

為了探究地表UCT實驗值q的變異系數(shù)對模型結(jié)果的影響，分別對3種模型加算COV=0.2、0.4的工況，其最終結(jié)果如圖8所示，隨著變異系數(shù)的增加，3種模型計算得到的安全系數(shù)標(biāo)準(zhǔn)差增大，而安全系數(shù)均值變化不明顯，導(dǎo)致可靠指標(biāo)均減小。3種工況的安全系數(shù)均值μ分別為1.39、1.70、1.55。若僅考慮安全系數(shù)均值μ，3種工況均能滿足規(guī)范要求，但可靠度僅模型2（深度趨勢）能滿足，并且模型3的安全系數(shù)離散程度高，變異性強，采用確定性分析或者平穩(wěn)隨機場模型偏于危險，討論模型3隨機分布情況就顯得更為重要。

3.3 趨勢參數(shù)b對安全系數(shù)分布的影響

為了探究模型3中趨勢參數(shù)b對安全系數(shù)分布的影響，分別單獨取隨機場離散公式（7）中b均值為1、1.5、2.5，b變異系數(shù)為0.1、0.2、0.4，其他參數(shù)保持不變，共6種工況，將其計算結(jié)果分別匯總于表2、表3，并與基礎(chǔ)工況（b均值為2、變異系數(shù)為0.3）對比。從表中可以看出，無論是提高b的均值還是變異系數(shù)，最終都表現(xiàn)出可靠指標(biāo)降低，失效概率增大的趨勢。其中，b均值的增加帶來式（7）中最終強度的減小，安全系數(shù)均值μ減小，標(biāo)準(zhǔn)差σ增大，變異性增強;b變異系數(shù)增大，安全系數(shù)均值μ基本未改變，但安全系數(shù)標(biāo)準(zhǔn)差σ增大，變異性增強。

4 結(jié)論

結(jié)合強度折減法、蒙特卡洛模擬方法，比較了3種隨機場模型對水泥土堤壩邊坡的失效風(fēng)險和滑動面位置的影響，得到以下結(jié)論：

1）對于水泥土邊坡的可靠度分析問題，采用平穩(wěn)隨機場或深度趨勢隨機場可能會高估邊坡的可靠度，使計算結(jié)果偏于危險;即便是安全系數(shù)均值滿足規(guī)范要求，但其標(biāo)準(zhǔn)差過大，最終的可靠指標(biāo)并不能滿足工程要求。

2）提出的不排水抗剪強度垂直排水板衰減非平穩(wěn)模型既能考慮沿深度方向的趨勢項，又能考慮垂直排水板方向的衰減趨勢，綜合表現(xiàn)為沿深度方向的近似周期性波動，與現(xiàn)有研究符合較好。在該模型中，趨勢參數(shù)b對隨機場模型的影響較大，當(dāng)其值和標(biāo)準(zhǔn)差增加，最終都將降低計算模型的可靠指標(biāo)，失效概率增大。

3）提出的不排水抗剪強度垂直排水板衰減非平穩(wěn)模型的邊坡破壞形式較為集中，大部分破壞面集中在堤壩土土體內(nèi)，且多為從坡底貫穿坡頂?shù)幕瑒用妗?/p>

參考文獻：

[1] GRIFFITHS D V， HUANG J S， FENTON G A. Influence of spatial variability on slope reliability using 2-D random fields [J]. Journal of Geotechnical and Geoenvironmental Engineering， 2009， 135（10）： 1367-1378.

[2] PHOON K K， KULHAWY F H. Characterization of geotechnical variability [J]. Canadian Geotechnical Journal， 1999， 36（4）： 612-624.

[3] 張明， 趙有明. 基于隨機場理論的吹填土土性參數(shù)空間變異性[J]. 中國港灣建設(shè)， 2012， 32（3）： 7-9， 21.

ZHANG M， ZHAO Y M. Spatial variability of property parameters of dredged fill based on random field theory [J]. China Harbour Engineering， 2012， 32（3）： 7-9， 21. （in Chinese）

[4] 周洋， 蒲訶夫， 李展毅， 等. 水平排水板-真空預(yù)壓聯(lián)合處理高含水率疏浚淤泥模型試驗研究[J]. 巖石力學(xué)與工程學(xué)報， 2019， 38（Sup1）： 3246-3251.

ZHOU Y， PU H F， LI Z Y， et al. Experimental investigations on treatment of dredged slurry by vacuum-assisted prefabricated horizontal drains [J]. Chinese Journal of Rock Mechanics and Engineering， 2019， 38（Sup1）： 3246-3251. （in Chinese）

[5] 張文彬， 王賢奔， 謝錦波， 等. 水平排水板真空預(yù)壓法處理吹填流泥試驗研究[J]. 水運工程， 2016（7）： 122-127.

ZHANG W B， WANG X B， XIE J B， et al. Experimental study on dredger fill flow mud improved by vacuum preloading method employing horizontal drains [J]. Port & Waterway Engineering， 2016（7）： 122-127. （in Chinese）

[6] 高會強. 真空預(yù)壓下塑料排水板有效影響范圍的研究分析[D]. 廣州： 華南理工大學(xué)， 2017.

GAO H Q. Study on the effective range of plastic vertical drain under vacuum preloading [D]. Guangzhou： South China University of Technology， 2017. （in Chinese）

[7] 曹宇春. 固結(jié)后軟黏土不排水抗剪強度簡化計算方法比較[J]. 土木工程學(xué)報， 2014， 47（10）： 107-116.

CAO Y C. A comparison of simplified calculation methods of undrained shear strength of soft clays after consolidation [J]. China Civil Engineering Journal， 2014， 47（10）： 107-116. （in Chinese）

[8] 蔣水華， 李典慶， 周創(chuàng)兵， 等. 考慮參數(shù)空間變異性的非飽和土坡可靠度分析[J]. 巖土力學(xué)， 2014， 35（9）： 2569-2578.

JIANG S H， LI D Q， ZHOU C B，et al. Reliability analysis of unsaturated slope considering spatial variability [J]. Rock and Soil Mechanics， 2014， 35（9）： 2569-2578. （in Chinese）

[9] 鄭俊杰， 喬雅晴， 章榮軍. 被動加固區(qū)參數(shù)變異性對軟土深基坑變形行為的影響[J]. 土木與環(huán)境工程學(xué)報（中英文）， 2019， 41（6）： 1-8.

ZHENG J J， QIAO Y Q， ZHANG R J. Effect of spatial variability of engineering properties on holistic performance of passive zone improvement for deep excavation in soft soil [J]. Journal of Civil and Environmental Engineering， 2019， 41（6）： 1-8. （in Chinese）

[10] 蔣水華， 曾紹慧， 楊建華， 等. 不排水抗剪強度非平穩(wěn)隨機場模擬及邊坡可靠度分析[J]. 巖土力學(xué)， 2018， 39（3）： 1071-1081.

JIANG S H， ZENG S H， YANG J H， et al. Slope reliability analysis by simulation of non-stationary random field of undrained shear strength [J]. Rock and Soil Mechanics， 2018， 39（3）： 1071-1081. （in Chinese）

[11] GRIFFITHS D V， HUANG J， FENTON G A. Probabilistic slope stability analysis using RFEM with non-stationary random fields [C]//Geotechnical Safety and Risk V， 2015： 704-709.

[12] 祁小輝， 李典慶， 周創(chuàng)兵， 等. 考慮不排水抗剪強度空間變異性的條形基礎(chǔ)極限承載力隨機分析[J]. 巖土工程學(xué)報， 2014， 36（6）： 1095-1105.

QI X H， LI D Q， ZHOU C B， et al. Stochastic analysis of ultimate bearing capacity of strip footing considering spatial variability of undrained shear strength [J]. Chinese Journal of Geotechnical Engineering， 2014， 36（6）： 1095-1105. （in Chinese）

[13] LI D Q， QI X H， PHOON K K， et al. Effect of spatially variable shear strength parameters with linearly increasing mean trend on reliability of infinite slopes [J]. Structural Safety， 2014， 49： 45-55.

[14] 豆紅強， 王浩. 非平穩(wěn)隨機場下飽和滲透系數(shù)空間變異性的無限長邊坡穩(wěn)定概率分析[J]. 土木工程學(xué)報， 2017， 50（8）： 105-113， 128.

DOU H Q， WANG H. Probabilistic analysis of spatial variability of saturated hydraulic conductivity on infinite slope based on the non-stationary random field [J]. China Civil Engineering Journal， 2017， 50（8）： 105-113， 128. （in Chinese）

[15] LI D Q， QI X H， CAO Z J， et al. Reliability analysis of strip footing considering spatially variable undrained shear strength that linearly increases with depth [J]. Soils and Foundations， 2015， 55（4）： 866-880.

[16] SHU S， GAO Y F， WU Y X， et al. Bearing capacity and reliability analysis of spudcan foundations embedded at various depths based on the non-stationary random finite element method [J]. Applied Ocean Research， 2020， 100： 102182.

[17] WU Y X， ZHOU X H， GAO Y F， et al. Effect of soil variability on bearing capacity accounting for non-stationary characteristics of undrained shear strength [J]. Computers and Geotechnics， 2019， 110： 199-210.

[18] ZHANG R J， HASAN M S M S， ZHENG J J， et al. Effect of spatial variability of engineering properties on stability of a CSMC embankment [J]. Marine Georesources & Geotechnology， 2018， 36（1）： 91-99.

[19] LEE F H， LEE Y， CHEW S H， et al. Strength and modulus of marine clay-cement mixes [J]. Journal of Geotechnical and Geoenvironmental Engineering， 2005， 131（2）： 178-186.

[20] LIU Y， CHEN E J， QUEK S T， et al. Effect of spatial variation of strength and modulus on the lateral compression response of cement-admixed clay slab [J]. Géotechnique， 2015， 65（10）： 851-865.

[21] 瞿秀祥， 林杭， 王敏， 等. FLAC V邊坡穩(wěn)定剪切強度折減法的多線程并行算法[J]. 鐵道科學(xué)與工程學(xué)報， 2018， 15（7）： 1693-1699.

QU X X， LIN H， WANG M， et al. Study on the slope stability shear strength reduction method based on the multi-thread parallel computing technology in the FLAC V software [J]. Journal of Railway Science and Engineering， 2018， 15（7）： 1693-1699. （in Chinese）

[22] 李劍， 陳善雄， 余飛. 基于最大剪應(yīng)變增量的邊坡潛在滑動面搜索[J]. 巖土力學(xué)， 2013， 34（Sup1）： 371-378.

LI J， CHEN S X， YU F. A method for searching potential failure surface of slope based on maximum shear strain increment [J]. Rock and Soil Mechanics， 2013， 34（Sup1）： 371-378. （in Chinese）

[23] 水利水電工程邊坡設(shè)計規(guī)范： SL 386—2007 [S]. 北京： 中國水利水電出版社， 2007.

Design code for engineered slopes in water resources and hydropower projects： SL 386-2007 [S]. Beijing： China Water & Power Press， 2007. （in Chinese）

[24] 公路工程結(jié)構(gòu)可靠度設(shè)計統(tǒng)一標(biāo)準(zhǔn)： JTG 2120—2020 [S].北京： 人民交通出版社， 2020.

Unified standard for reliability design of highway engineering structures： JTG 2120-2020 [S]. Beijing： China Communications Press， 2020. （in Chinese）

[25] 公路路基設(shè)計規(guī)范： JTG D30—2015 [S]. 北京： 人民交通出版社， 2015.

Specifications for design of highway subgrades： JTG D30-2015 [S]. Beijing： China Communications Press， 2015. （in Chinese）

（編輯 胡玲）