基于強(qiáng)度折減法的土石混合體邊坡長(zhǎng)期穩(wěn)定性研究*

劉康琦 劉紅巖 祁小博

(①中國(guó)地質(zhì)大學(xué)(北京)工程技術(shù)學(xué)院， 北京 100083， 中國(guó)) (②自然資源部深部地質(zhì)鉆探技術(shù)重點(diǎn)實(shí)驗(yàn)室， 北京 100083， 中國(guó)) (③中國(guó)地質(zhì)環(huán)境監(jiān)測(cè)院， 北京 100081， 中國(guó))

0 引 言

土石混合體邊坡是第四紀(jì)以來(lái)形成的、由土和塊石組成的一種極度不均勻的巖土體邊坡。不同于一般的土質(zhì)邊坡或巖質(zhì)邊坡，由于塊石與土的共同存在，使土石混合體邊坡具有極大的不均勻性和不連續(xù)性。目前針對(duì)土石混合體邊坡的穩(wěn)定性已有較多研究，油新華等(2002)通過(guò)對(duì)土石混合體進(jìn)行野外水平推剪試驗(yàn)，得出土石混合體的變形特點(diǎn)和相關(guān)的抗剪強(qiáng)度參數(shù)，并發(fā)現(xiàn)土石混合體在剪切破壞時(shí)會(huì)繞過(guò)大的石塊。徐文杰等(2009)將土石混合體定義為第四紀(jì)以來(lái)形成的，由具有一定工程尺度、強(qiáng)度較高的塊石、細(xì)粒土體及孔隙構(gòu)成且具有一定含石量的極端不均勻松散巖土介質(zhì)系統(tǒng)，并通過(guò)可視粒徑、土石閾值、“土”與“石”的強(qiáng)度特征及含石量4個(gè)參數(shù)對(duì)土石混合體的概念進(jìn)行了補(bǔ)充。Zhou et al.(2019)利用室內(nèi)循環(huán)凍融試驗(yàn)和PFC3D數(shù)值模擬研究了凍融循環(huán)對(duì)土石混合體強(qiáng)度的影響。徐文杰(2009)基于對(duì)周家灣滑坡的大量地質(zhì)調(diào)查研究，對(duì)土石混合體滑坡體進(jìn)行了工程地質(zhì)分區(qū)，并建立了三維地質(zhì)模型，運(yùn)用有限元強(qiáng)度折減法分析了其穩(wěn)定性。楊忠平等(2017)基于PFC2D接觸黏結(jié)模型的離散元數(shù)值模型，研究了含石量變化對(duì)土石混合體剪切特性的影響。龔健等(2017)通過(guò)對(duì)不同含石量的土石混合體進(jìn)行靜力超載試驗(yàn)，得出土石混合體的力學(xué)行為與含石量的關(guān)系，另基于PIV分析方法，研究了剪切破壞時(shí)剪切帶的發(fā)展規(guī)律。張森等(2016)提出基于顆粒流區(qū)域填充技術(shù)的土石混合體有限元隨機(jī)生成技術(shù)，并利用強(qiáng)度折減法研究了含石量對(duì)邊坡穩(wěn)定性的影響規(guī)律。Sun et al.(2014)通過(guò)對(duì)冷谷水電站馬河堆積體進(jìn)行工程類(lèi)比、顆粒分析試驗(yàn)、直剪試驗(yàn)、現(xiàn)場(chǎng)大型推剪試驗(yàn)、數(shù)值模擬試驗(yàn)、參數(shù)反演等研究，對(duì)比分析了土石混合體的抗剪強(qiáng)度，并對(duì)馬河堆積體的穩(wěn)定性進(jìn)行了分析。Fan et al.(2015)使用離散元數(shù)值模擬方法，研究了土石混合體在循環(huán)荷載作用下的力學(xué)特性，結(jié)果表明峰值應(yīng)力會(huì)隨著動(dòng)應(yīng)力幅值和頻率的增大而增大，但會(huì)隨著塊石所占比例的增大而減小，同時(shí)彈性模量會(huì)隨著動(dòng)應(yīng)力幅值和頻率的增大而減小。鄭博寧等(2019)提出了基于CT掃描技術(shù)和PFC三維顆粒流程序的土石混合體三維建模技術(shù)，并通過(guò)與三軸剪切試驗(yàn)過(guò)程中礫石的空間運(yùn)動(dòng)規(guī)律進(jìn)行對(duì)比分析，驗(yàn)證了其建模方法的可行性。Wang et al.(2019)利用離心模型試驗(yàn)通過(guò)考慮不同體積塊石含量，研究了土石混合體邊坡在地表荷載作用下的破壞行為。Lu et al.(2018)使用離散元法研究了含石量、坡高比和塊石大小對(duì)土石混合體邊坡穩(wěn)定性的影響。Yang et al.(2019)利用修正數(shù)值流形法研究了連續(xù)開(kāi)挖條件下土石混合體邊坡的穩(wěn)定性。

然而目前對(duì)土石混合體邊坡穩(wěn)定性的研究，基本都忽略了土體的蠕變效應(yīng)。研究表明(陳衛(wèi)兵等， 2008； 蔣海飛等， 2013； 陳廷君等， 2019)，土體蠕變會(huì)降低坡體的穩(wěn)定性系數(shù)，對(duì)邊坡穩(wěn)定性造成不利影響。本文首先利用數(shù)字圖像處理技術(shù)對(duì)我國(guó)某水電站庫(kù)區(qū)的一個(gè)土石混合體邊坡進(jìn)行處理，最大化的獲得土石混合體邊坡的細(xì)觀(guān)模型，然后利用FLAC3D程序中強(qiáng)度折減法對(duì)其進(jìn)行穩(wěn)定性分析，并考慮了土體的蠕變性質(zhì)，對(duì)比分析了考慮土體蠕變的土石混合體邊坡與不考慮土體蠕變的土石混合體邊坡的變形特征。然后通過(guò)對(duì)采用不同土體蠕變參數(shù)的土石混合體邊坡進(jìn)行穩(wěn)定性分析，討論了土體蠕變參數(shù)對(duì)土石混合體邊坡穩(wěn)定性的影響。

1 土石混合體邊坡模型的建立

1.1 數(shù)字圖像處理

數(shù)字圖像處理(Digital Image Processing)被最早用來(lái)改善圖像的質(zhì)量，而隨著計(jì)算機(jī)科學(xué)的不斷發(fā)展，其開(kāi)始廣泛應(yīng)用于航空航天、醫(yī)學(xué)、軍事、工程等領(lǐng)域。近年來(lái)，數(shù)字圖像處理技術(shù)也越來(lái)越多地被應(yīng)用到巖土工程界(Chen et al.， 2004； Xu et al.， 2008; 徐文杰等， 2008； 廖秋林等， 2010； 王學(xué)濱， 2018)，它將圖像經(jīng)過(guò)數(shù)學(xué)變換后轉(zhuǎn)換為數(shù)字信號(hào)并利用計(jì)算機(jī)對(duì)其進(jìn)行分析處理，可將圖像定義為一個(gè)非零的和有限的二維函數(shù)f(x，y)，其中(x，y)表示空間或平面坐標(biāo)，f表示圖像在該點(diǎn)處的灰度值。在對(duì)土石混合體邊坡進(jìn)行處理時(shí)，可根據(jù)土與塊石灰度值的不同，通過(guò)設(shè)定閾值，將土與石塊分別開(kāi)來(lái)，這樣利于提取土石混合體邊坡中土與塊石的分布信息。

在進(jìn)行圖像處理時(shí)依托于OpenCV強(qiáng)大的圖像處理功能及Python語(yǔ)言的簡(jiǎn)潔高效，OpenCV是由C++語(yǔ)言編寫(xiě)的一種計(jì)算機(jī)視覺(jué)庫(kù)(Bradski， 2000)，具有輕量、高效且開(kāi)源的優(yōu)點(diǎn)，并提供了多種語(yǔ)言的接口。

圖 1 實(shí)際的土石混合體邊坡圖像Fig. 1 Image of soil-rock mixture slope

圖 1為我國(guó)某水電站庫(kù)區(qū)的一個(gè)土石混合體邊坡的斷面圖像。對(duì)其具體處理時(shí)過(guò)程如下：

(1)新建.py文件。與需要處理的土石混合體邊坡圖像放置在同一文件夾下。

(2)邊坡圖像二值化處理。依次對(duì)圖像進(jìn)行均值濾波器去噪處理→圖像灰度化處理，將原來(lái)的RGB彩色圖像轉(zhuǎn)化為灰度圖像→灰度圖像高斯濾波器去噪處理，以消除灰度圖內(nèi)的高斯噪聲→自適應(yīng)閾值二值化處理，可在圖像局部范圍內(nèi)自適應(yīng)一個(gè)閾值對(duì)圖像進(jìn)行二值化處理，將圖像上石塊的灰度值設(shè)置為0，呈黑色，土體的灰度值設(shè)置為255，呈白色，使用自適應(yīng)閾值二值化可最大化消除光照等因素對(duì)圖像造成的影響→腐蝕算法去噪處理。

(3)細(xì)節(jié)調(diào)整。使用Photoshop軟件對(duì)二值化后的邊坡圖像內(nèi)不理想的區(qū)域進(jìn)行手工處理，使二值化圖像最大化接近真實(shí)土石分布情況。

(4)光滑及矢量化處理。由于生成的二值化圖像內(nèi)石塊邊界粗糙會(huì)對(duì)數(shù)值模擬造成困難甚至?xí)箶?shù)值模擬失敗，因此有必要對(duì)其進(jìn)行邊界光滑處理。另數(shù)字圖像并不能直接用于數(shù)值模擬，因此需要對(duì)其進(jìn)行矢量化處理，生成AutoCAD可接受的DXF文件。圖 2即為對(duì)圖 1進(jìn)行矢量化處理后的結(jié)果，土石比約為14︰1，最大粒徑為1.37im。

圖 2 矢量化土石混合體邊坡模型Fig. 2 Vector model of soil-rock mixture slope

1.2 建立模型

在進(jìn)行數(shù)值計(jì)算時(shí)塊石距離邊坡邊界過(guò)近會(huì)產(chǎn)生邊界效應(yīng)(徐文杰等， 2008)，為減少由于邊界效應(yīng)產(chǎn)生的影響，在邊坡模型的右側(cè)及左側(cè)坡腳處使用均質(zhì)土體進(jìn)行了延拓，在AutoCAD中面域化后生成DXF文件導(dǎo)入ANSYS軟件進(jìn)行網(wǎng)格劃分，然后轉(zhuǎn)化為可供FLAC3D軟件接受的文件格式。最終在FLAC3D中生成模型(圖 3)，共生成節(jié)點(diǎn)10i906個(gè)，劃分網(wǎng)格單元5455個(gè)，其中土體單元4736個(gè)，塊石單元719個(gè)。

圖 3 FLAC3D中生成的土石混合體邊坡模型圖Fig. 3 Model of soil-rock mixture slope in FLAC3D

2 考慮蠕變性質(zhì)的強(qiáng)度折減法

2.1 強(qiáng)度折減法

由于傳統(tǒng)的靜力平衡法需要預(yù)先指定滑動(dòng)面，不適于內(nèi)部塊石分布復(fù)雜的土石混合體邊坡，因此采用強(qiáng)度折減法對(duì)邊坡進(jìn)行穩(wěn)定性分析(Roth et al.， 1999)。強(qiáng)度折減法中將安全系數(shù)定義為邊坡剛好達(dá)到臨界失穩(wěn)狀態(tài)時(shí)，對(duì)巖土體的抗剪強(qiáng)度的折減程度。以莫爾-庫(kù)侖強(qiáng)度屈服準(zhǔn)則為例，采用強(qiáng)度折減法對(duì)邊坡進(jìn)行穩(wěn)定性分析時(shí)，通過(guò)不斷的對(duì)邊坡的黏聚力c和內(nèi)摩擦角φ除以折減系數(shù)進(jìn)行試算，直到邊坡達(dá)到失穩(wěn)狀態(tài)，此時(shí)的折減系數(shù)即是該邊坡的安全系數(shù)，其折減過(guò)程如下：

(1)

(2)

式中，ctrial、φtrial為土體折減后的黏聚力、內(nèi)摩擦角；Ftrial為折減系數(shù)；c、φ為巖土體折減前的黏聚力、內(nèi)摩擦角。

2.2 蠕變本構(gòu)模型的選取

在對(duì)蠕變邊坡進(jìn)行穩(wěn)定性分析時(shí)，選取合適的黏彈塑性蠕變本構(gòu)模型及參數(shù)可充分反映蠕變對(duì)邊坡變形及穩(wěn)定性的影響，直接關(guān)系到分析結(jié)果的合理性與正確性。一般來(lái)說(shuō)，在選取蠕變模型時(shí)，應(yīng)先進(jìn)行室內(nèi)試驗(yàn)，在試驗(yàn)數(shù)據(jù)的基礎(chǔ)上利用蠕變理論選取合適的本構(gòu)模型，然后利用回歸分析和最小二乘法得出模型的參數(shù)。本文選取FLAC3D中Burgers黏彈性模型與Mohr-coulomb模型串聯(lián)而成的Cvisc模型，其構(gòu)成如下(圖4)：

圖 4 FLAC3D中Cvisc模型示意圖Fig. 4 Illustrations of Cvisc rheological model in FLAC3D

在FLAC3D中，Cvisc模型由Kelvin體、Maxwell體和一個(gè)塑性體構(gòu)成，其中，σ表示作用在巖土體上的應(yīng)力；σt表示巖土體的屈服強(qiáng)度；EM、EK分別表示彈性模量和黏彈性模量；ηM、ηK分別表示Maxwell黏性參數(shù)和Kelvin黏性參數(shù)。在巖土體的蠕變變形中，若其具有穩(wěn)定蠕變或加速蠕變性質(zhì)，則即使不對(duì)巖土體的黏聚力和內(nèi)摩擦角進(jìn)行折減，邊坡也會(huì)一直變形直至破壞，這樣則無(wú)法體現(xiàn)邊坡的變形失穩(wěn)與巖土體強(qiáng)度之間的關(guān)系。因此需要將巖土體設(shè)定為只具有衰減蠕變性質(zhì)，在FLAC3D中，只需將Maxwell黏性系數(shù)ηM設(shè)置為無(wú)限大即可，此時(shí)，Cvisc模型則退化為廣義Kelvin黏彈塑性模型。

三維形式下模型應(yīng)變?cè)隽科坑蒑axwell體、Kelvin體和塑性體3部分組成。分別為：

(3)

(4)

(5)

(6)

2.3 邊坡失穩(wěn)判據(jù)的選取

在使用強(qiáng)度折減法對(duì)邊坡進(jìn)行穩(wěn)定性分析時(shí)，需要選取合適的失穩(wěn)判據(jù)來(lái)判斷邊坡是否達(dá)到臨界失穩(wěn)狀態(tài)。傳統(tǒng)的強(qiáng)度折減法中，一般有3種方法用來(lái)作為判據(jù)：

(1)數(shù)值計(jì)算不收斂。當(dāng)巖土體的強(qiáng)度參數(shù)折減到一定程度使邊坡達(dá)到失穩(wěn)破壞狀態(tài)時(shí)，F(xiàn)LAC3D程序無(wú)法得到一個(gè)確定的解來(lái)滿(mǎn)足靜力平衡和應(yīng)力-應(yīng)變關(guān)系。

(2)關(guān)鍵點(diǎn)位移突變?？梢酝ㄟ^(guò)監(jiān)測(cè)若干關(guān)鍵點(diǎn)的位移來(lái)判定邊坡是否破壞，當(dāng)強(qiáng)度折減到一定程度使邊坡達(dá)到臨界破壞狀態(tài)時(shí)，位移會(huì)突然增大。

(3)具有貫通的塑性區(qū)。當(dāng)塑性區(qū)從坡頂一直貫通到坡腳時(shí)，可認(rèn)為邊坡已達(dá)到失穩(wěn)狀態(tài)。在FLAC3D中，可通過(guò)查看剪應(yīng)變速率云圖來(lái)觀(guān)察塑性區(qū)是否貫通。

以上3種方法即為傳統(tǒng)的強(qiáng)度折減法中邊坡是否失穩(wěn)的判據(jù)，而在考慮蠕變的土石混合體邊坡穩(wěn)定性分析中并不適用。一方面，由于塊石的存在，很可能使邊坡在破壞時(shí)不會(huì)出現(xiàn)明顯貫通的塑性區(qū)； 另一方面，在考慮巖土體的蠕變性質(zhì)時(shí)，很難以數(shù)值收斂作為邊坡是否失穩(wěn)破壞的判據(jù)。因此，可將關(guān)鍵點(diǎn)位移在經(jīng)過(guò)一定時(shí)間后是否穩(wěn)定及位移在強(qiáng)度折減到某個(gè)程度時(shí)是否突變作為土石混合體邊坡是否失穩(wěn)的判據(jù)(蔣海飛等， 2013)。

2.4 模型驗(yàn)證

為驗(yàn)證本文方法的正確性，使用自編蠕變強(qiáng)度折減法對(duì)陳衛(wèi)兵等(2008)的算例進(jìn)行對(duì)比。在對(duì)算例進(jìn)行試算時(shí)，取土體參數(shù)如下：重度γ=20ikN·m-3，黏聚力c=35 kPa，內(nèi)摩擦角φ=18°，剪脹角Ψ=18°，彈性模量和黏彈性模量分別為E=50 MPa、Ek=300 MPa，黏性系數(shù)ηk=5 GPa·d， 泊松比v=0.￣35，蠕變時(shí)間為1ia。通過(guò)不斷折減土體強(qiáng)度參數(shù)，最后求得在考慮蠕變效應(yīng)時(shí)該算例的安全系數(shù)為1.48，這與陳衛(wèi)兵等(2008)的結(jié)果一致，因此可以認(rèn)定本文所使用自編強(qiáng)度折減法是可靠的。圖 5為該算例與使用自編強(qiáng)度折減法取折減系數(shù)為1.48時(shí)的剪應(yīng)變速率云圖。

圖 5 折減系數(shù)為1.48時(shí)剪應(yīng)變速率云圖Fig. 5 Contours of shear strain rate at reduction factor=1.48a. 本文計(jì)算結(jié)果； b. 文獻(xiàn)計(jì)算結(jié)果(陳衛(wèi)兵等，2008)

3 土石混合體邊坡穩(wěn)定性分析

計(jì)算模型的邊界條件為：模型右側(cè)及左側(cè)邊界采用水平約束，邊坡底面采用水平、垂直全部約束。

3.1 考慮蠕變的土石混合體邊坡穩(wěn)定性分析

為作對(duì)比，首先研究不考慮蠕變效應(yīng)時(shí)的土石混合體邊坡穩(wěn)定性，土體和塊石均采用Mohr-Coulomb模型。基于蔣海飛等(2013)的研究，本文所選巖土體模型參數(shù)如表 1所示。

表 1 不考慮蠕變時(shí)模型參數(shù)表Table 1 Physico-mechanical parameters of S-RMS without considering rheology

圖 6 不考慮蠕變效應(yīng)時(shí)邊坡剪應(yīng)變速率云圖Fig. 6 Contours of shear strain rate without considering rheology

表 2 土體蠕變模型參數(shù)表Table 2 Physico-mechanical parameters of soil with considering rheology

圖 6為不考慮蠕變效應(yīng)時(shí)土石混合體邊坡強(qiáng)度折減法的分析結(jié)果，可以看到，由于塊石的存在，在邊坡內(nèi)部形成了多條滑帶，而且滑帶在靠近塊石的部位有明顯的“繞石”現(xiàn)象，由強(qiáng)度折減法得出的折減系數(shù)為1.85，即不考慮蠕變效應(yīng)時(shí)此土石混合體邊坡的安全系數(shù)為1.85。

在對(duì)考慮蠕變效應(yīng)的土石混合體邊坡進(jìn)行穩(wěn)定性分析時(shí)，在可能的滑帶上及坡頂坡腳處選擇關(guān)鍵點(diǎn)來(lái)監(jiān)測(cè)其水平位移情況(陳衛(wèi)兵等， 2008)。石塊采用Mohr-Coulomb模型，土體采用Cvisc黏彈塑性蠕變模型，蠕變時(shí)間為1ia。其Cvisc模型參數(shù)如表 2所示。

圖 7 考慮蠕變效應(yīng)時(shí)邊坡剪應(yīng)變速率云圖Fig. 7 Contours of shear strain rate with considering rheology

圖 8 折減系數(shù)為1.46時(shí)監(jiān)測(cè)點(diǎn)位移曲線(xiàn)Fig. 8 Horizontal displacement vs time of monitoring points at reduction factor=1.46

圖 9 折減系數(shù)為1.47時(shí)監(jiān)測(cè)點(diǎn)位移曲線(xiàn)Fig. 9 Horizontal displacement vs time of monitoring points at reduction factor=1.47

圖 10 關(guān)鍵點(diǎn)水平位移與折減系數(shù)的關(guān)系曲線(xiàn)Fig. 10 Horizontal displacements of monitoring points at different reduction factors with considering rheology

通過(guò)自編強(qiáng)度折減法來(lái)對(duì)考慮蠕變效應(yīng)的土石混合體邊坡進(jìn)行穩(wěn)定性分析，其中蠕變參數(shù)不變，折減系數(shù)從1.30增加到1.47。圖 7為折減系數(shù)為1.47時(shí)的剪應(yīng)變速率云圖，圖 8和9為折減系數(shù)分別是1.46和1.47時(shí)監(jiān)測(cè)點(diǎn)水平位移曲線(xiàn)，可以看到，當(dāng)折減系數(shù)為1.46時(shí)，各監(jiān)測(cè)點(diǎn)水平位移在經(jīng)過(guò)短暫的瞬時(shí)變形后，位移增加量逐漸減小，最終趨于穩(wěn)定，位移增加量近似為0。當(dāng)折減系數(shù)增加到1.47時(shí)，雖然監(jiān)測(cè)點(diǎn)3和監(jiān)測(cè)點(diǎn)4的位移最終趨于穩(wěn)定，但監(jiān)測(cè)點(diǎn)1和監(jiān)測(cè)點(diǎn)2在經(jīng)歷1ia的蠕變時(shí)間之后，位移仍未收斂。

各關(guān)鍵點(diǎn)水平位移與強(qiáng)度折減系數(shù)的關(guān)系曲線(xiàn)如圖 10所示，可以看出，折減系數(shù)由1.30增加到1.46時(shí)，關(guān)鍵點(diǎn)總位移并未發(fā)生太大變化，而當(dāng)折減系數(shù)增加到1.47后，可以明顯地看到位移曲線(xiàn)有一個(gè)突變，位移增大較為明顯，而此時(shí)其蠕變曲線(xiàn)也未到達(dá)穩(wěn)定狀態(tài)，因此可以認(rèn)定，考慮蠕變時(shí)該土石混合體邊坡的安全系數(shù)為1.47，比不考慮蠕變時(shí)該土石混合體邊坡的安全系數(shù)下降了20.5%。以監(jiān)測(cè)點(diǎn)2為例，取折減系數(shù)為1.46，不考慮土體蠕變效應(yīng)時(shí)，其水平位移僅為0.002im，而考慮土體蠕變效應(yīng)時(shí)，水平位移為0.32im，是不考慮蠕變效應(yīng)時(shí)的160倍，可見(jiàn)土體蠕變對(duì)邊坡變形影響極大； 在不考慮蠕變效應(yīng)的情況下取折減系數(shù)為1.84，監(jiān)測(cè)點(diǎn)2的水平位移為0.15im，因此若將0.15im作為預(yù)測(cè)邊坡失穩(wěn)的位移量則顯得過(guò)于保守。

圖 11 不同黏性系數(shù)時(shí)邊坡安全系數(shù)Fig. 11 FOS with different viscosity coefficient

3.2 蠕變參數(shù)對(duì)安全系數(shù)的影響

巖土體的黏性系數(shù)是反映其蠕變程度的一個(gè)重要參數(shù)，為探索黏性系數(shù)對(duì)土石混合體邊坡安全系數(shù)的影響關(guān)系，在其他力學(xué)參數(shù)不變的情況下，分別取Cvisc模型中的Kelvin黏性系數(shù)為1.19e3iPa·d、1.19e4iPa·d、1.19e5iPa·d和1.19e6iPa·d，并對(duì)邊坡進(jìn)行穩(wěn)定性分析。圖 11為取不同黏性系數(shù)時(shí)的邊坡安全系數(shù)關(guān)系圖，由圖可以看出，黏性系數(shù)對(duì)于邊坡的變形及穩(wěn)定影響較大，隨著黏性系數(shù)的增大，邊坡安全系數(shù)隨之減小，說(shuō)明黏性系數(shù)的增大對(duì)于邊坡的變形及穩(wěn)定具有不利影響。

4 結(jié) 論

土石混合體邊坡是自然界中廣泛存在的一種巖土體邊坡，以我國(guó)某水電站庫(kù)區(qū)的一個(gè)土石混合體邊坡為例，首先運(yùn)用數(shù)字圖像處理技術(shù)對(duì)該邊坡進(jìn)行處理，得到土石混合體邊坡的細(xì)觀(guān)物理模型，然后在FLAC3D中運(yùn)用強(qiáng)度折減法對(duì)該邊坡進(jìn)行了穩(wěn)定性分析，并考慮了土體的蠕變性質(zhì)，與未考慮蠕變性質(zhì)的土石混合體邊坡穩(wěn)定性分析結(jié)果進(jìn)行比較，得出以下結(jié)論：

(1)土石混合體邊坡不同于一般均質(zhì)土體或巖質(zhì)邊坡，在達(dá)到臨界失穩(wěn)狀態(tài)時(shí)，在邊坡內(nèi)部會(huì)出現(xiàn)多條滑帶，且滑帶具有明顯的“繞石”現(xiàn)象。

(2)不考慮土體蠕變性質(zhì)的土石混合體邊坡的安全系數(shù)為1.85，而考慮了土體蠕變性質(zhì)的土石混合體邊坡的安全系數(shù)為1.47，相比下降了20.5%； 在取折減系數(shù)同為1.46的情況下，考慮土體蠕變效應(yīng)的邊坡監(jiān)測(cè)點(diǎn)水平位移是不考慮蠕變效應(yīng)時(shí)的160倍之多。因此巖土體的蠕變性質(zhì)不利于邊坡的穩(wěn)定，在對(duì)巖土體邊坡進(jìn)行穩(wěn)定性分析時(shí)不應(yīng)忽略其蠕變性質(zhì)。

(3)研究了Kelvin黏性系數(shù)對(duì)邊坡進(jìn)行長(zhǎng)期穩(wěn)定性的影響，認(rèn)為邊坡安全系數(shù)隨著Kelvin黏性系數(shù)的增大而減小，說(shuō)明黏性系數(shù)對(duì)于巖土體邊坡的變形及穩(wěn)定具有不利影響。