宿曉萍，潘明遠(yuǎn)，王亭亭

（長(zhǎng)春工程學(xué)院土木工程學(xué)院，長(zhǎng)春130012）

0 前言

滑坡存在的地質(zhì)條件復(fù)雜多樣，幾何形態(tài)也多種多樣，是最常見的自然災(zāi)害之一。即使是邊坡局部具有不穩(wěn)定性，同樣可能導(dǎo)致邊坡整體失穩(wěn)，甚而造成巨大的經(jīng)濟(jì)損失和人員傷亡?；禄瑒?dòng)通常是由于各因素共同作用而導(dǎo)致應(yīng)力重分布，一般用安全系數(shù)來(lái)衡量。在公路開挖形成的邊坡以及自然邊坡的穩(wěn)定性計(jì)算中，通常是通過巖質(zhì)邊坡穩(wěn)定性分析對(duì)邊坡平衡條件進(jìn)行分析［1］。然而，巖質(zhì)邊坡內(nèi)的滑動(dòng)通常被認(rèn)為是一種多因子影響的結(jié)果。早期研究認(rèn)為巖質(zhì)邊坡的穩(wěn)定性主要取決于巖體中的不連續(xù)面，而非巖石本身的強(qiáng)度。隨后有研究表明邊坡幾何形態(tài)與邊坡的滑動(dòng)方向是使邊坡破壞最大的因素，且邊坡的破壞往往發(fā)生在不連續(xù)面產(chǎn)狀的方向上［2］。因此，巖體的力學(xué)性質(zhì)和物理特性可以通過交叉的不連續(xù)節(jié)理組決定。尤其是大地構(gòu)造，如節(jié)理斷裂，隨巖體的強(qiáng)度特性和穩(wěn)定性問題的研究具有重要意義。

1 研究區(qū)的區(qū)域地質(zhì)

基于地貌的研究，研究區(qū)位于山地北邊的一部分，地質(zhì)研究表明，這一地區(qū)由一組鎂鐵質(zhì)、超鎂鐵的變質(zhì)巖和第四紀(jì)覆蓋層組成。根據(jù)研究區(qū)巖石的研究，巖石單元由變質(zhì)的玄武巖和變質(zhì)的輝長(zhǎng)巖以及變質(zhì)的沉積巖（葉片狀的千枚巖）組成，并經(jīng)過了高度的風(fēng)化。關(guān)于該地區(qū)的構(gòu)造方面，明顯存在逆斷層和褶皺的影響。該區(qū)域的斷層整體是沿著北西—南東方向延伸，方向?yàn)楸北睎|—南南西。

巖質(zhì)邊坡滑動(dòng)的危險(xiǎn)評(píng)估是驗(yàn)證此區(qū)域內(nèi)住宅樓穩(wěn)定性的關(guān)鍵問題。為此，本文選擇了該區(qū)域一個(gè)復(fù)雜的邊坡進(jìn)行研究，進(jìn)行邊坡穩(wěn)定性參數(shù)的評(píng)價(jià)與各項(xiàng)參數(shù)的研究。這一復(fù)雜的邊坡問題是研究區(qū)內(nèi)最大的問題之一，本文還將利用概率分析法和靈敏度分析法計(jì)算對(duì)其產(chǎn)生影響的其他參數(shù)（例如邊坡的產(chǎn)狀和高度、不連續(xù)面的產(chǎn)狀）以及一些地質(zhì)力學(xué)性質(zhì)（例如內(nèi)聚力、內(nèi)摩擦角和孔隙水壓力）。

2 研究方法

本研究是通過執(zhí)行標(biāo)準(zhǔn)的地質(zhì)技術(shù)來(lái)實(shí)現(xiàn)的。將數(shù)據(jù)收集和分析分為4個(gè)步驟：野外調(diào)查，室內(nèi)試驗(yàn)，穩(wěn)定性分析和靈敏度分析。

2.1 野外調(diào)查

在邊坡穩(wěn)定分析中的重要一步是實(shí)地調(diào)查。因?yàn)椴贿B續(xù)面在巖質(zhì)邊坡不穩(wěn)定研究中起著非常重要的作用，根據(jù)國(guó)際力學(xué)巖石學(xué)會(huì)的力學(xué)標(biāo)準(zhǔn)（ISRM1978）對(duì)巖質(zhì)邊坡的特征進(jìn)行了研究。首先，沿著溝內(nèi)110不連續(xù)面通過測(cè)線法進(jìn)行線性的統(tǒng)計(jì)，同時(shí)統(tǒng)計(jì)不連續(xù)面的產(chǎn)狀、張開度、填充情況、粗糙度以及節(jié)理間距等。其次，識(shí)別、確定地質(zhì)巖石類型和地質(zhì)結(jié)構(gòu)，溝傾向傾角以及地下水條件。最后，采集大量樣本，為室內(nèi)的實(shí)驗(yàn)做準(zhǔn)備。

2.2 室內(nèi)試驗(yàn)

基于莫爾—庫(kù)倫準(zhǔn)則，不連續(xù)面的剪切強(qiáng)度是研究巖質(zhì)邊坡穩(wěn)定性最重要的參數(shù)之一，該參數(shù)可以通過內(nèi)聚力和內(nèi)摩擦角來(lái)計(jì)算得出。

其過程簡(jiǎn)要如下：計(jì)算巖土體參數(shù)，試驗(yàn)樣本為2個(gè)鉆井挖掘出的巖芯樣本。樣本是從同一個(gè)地區(qū)中被挖掘而出的，樣本在105℃下烘24h。試驗(yàn)的巖心樣本的直徑為54～63mm，高度為63～150mm。樣品在采集過程中利用塑料袋包裝保存以保證其天然含水率。此外，一些樣本是通過地質(zhì)錘從溝渠內(nèi)巖石的露頭部分中未風(fēng)化的部分中取得，這樣可以評(píng)價(jià)巖石的類型和巖石參數(shù)。

根據(jù)規(guī)程ASTM D2938（2002），本文選取沒有節(jié)理的樣本進(jìn)行無(wú)側(cè)限抗壓強(qiáng)度實(shí)驗(yàn)，實(shí)驗(yàn)是在保持了天然含水率條件下進(jìn)行的。同樣，根據(jù)規(guī)程ASTM D5607（1995）選取一個(gè)不連續(xù)面的樣本進(jìn)行直剪試驗(yàn)，用來(lái)確定剪切強(qiáng)度。樣本分為2部分：一部分是在飽和條件下，即置于水中48h；另一部分是干燥條件下，即在105℃下烘24h。樣品的剪切強(qiáng)度是在不同的正壓力條件下重復(fù)進(jìn)行剪切強(qiáng)度測(cè)定來(lái)衡量的。此外，測(cè)定巖石樣本的強(qiáng)度參數(shù)，包括內(nèi)聚力、內(nèi)摩擦角、比重、無(wú)側(cè)限抗壓強(qiáng)度等。本例節(jié)理的特征見表1。

2.3 穩(wěn)定性分析

通常，動(dòng)態(tài)分析被認(rèn)為是評(píng)價(jià)潛在的不穩(wěn)定塊體和山坡上的幾何塊體的一種普遍方法。動(dòng)態(tài)分析是在沒有外力作用在滑動(dòng)面上時(shí)，處理邊坡傾向傾角和不連續(xù)面之間的關(guān)系，從未計(jì)算出失效破壞的可能性。

一個(gè)不穩(wěn)定塊體運(yùn)動(dòng)的可能性可以通過運(yùn)動(dòng)分析進(jìn)行評(píng)估。因此，本文采用了這種不穩(wěn)定性分析機(jī)制。為了進(jìn)行運(yùn)動(dòng)學(xué)分析，使用Wyllie and Mah（2004）的方法，運(yùn)用Dip軟件求得平面破壞、楔形體破壞以及傾倒破壞的可能性。利用Wyllie and Mah的方法對(duì)安全系數(shù)進(jìn)行確定性分析，即采用為每個(gè)參數(shù)的常量和結(jié)果代表指定的值安全系數(shù)（FS）。具體用到的參數(shù)包括：不連續(xù)面產(chǎn)狀、巖石的密度、邊坡的產(chǎn)狀、內(nèi)聚力、內(nèi)摩擦角和地下水條件。地震力作用在分析中沒有考慮。此外，確定性分析、概率分析以及敏感度分析都是基于莫爾—庫(kù)倫準(zhǔn)則。

巖質(zhì)邊坡穩(wěn)定性分析的主要問題是確定巖土體參數(shù)的不確定性，這與邊坡的天然物質(zhì)組成相關(guān)。不連續(xù)面的產(chǎn)狀方向變化非常大，甚至巖體內(nèi)部的節(jié)理組產(chǎn)狀變化也很大。同時(shí)，沿著不連續(xù)面，巖體抗剪強(qiáng)度會(huì)隨著節(jié)理的充填情況和粗糙程度的不同而不同。巖體剪切特性和外力條件，如水位和地震力等都隨時(shí)間而變化。因此，選擇一個(gè)恰當(dāng)?shù)挠行?shù)值在利用確定性方法評(píng)價(jià)邊坡穩(wěn)定性的過程中是非常關(guān)鍵的。因此，本文旨在確定FS的概率。概率分析法可用來(lái)克服利用前述的確定性方法解決不了的問題。

概率分析是一個(gè)系統(tǒng)化的方法，通過建模確定不確定性和變化的評(píng)價(jià)參數(shù)。在這個(gè)分析中，F(xiàn)S表示為一個(gè)隨機(jī)變量，可以取代失效的概率，從而評(píng)價(jià)邊坡穩(wěn)定值。概率分析包括2個(gè)步驟：首先，確定巖土參數(shù)的分析基本統(tǒng)計(jì)元素（平均值和標(biāo)準(zhǔn)偏差）和概率密度函數(shù)，用來(lái)預(yù)測(cè)巖土參數(shù)的隨機(jī)數(shù)據(jù)；其次，在第一步結(jié)果的基礎(chǔ)上進(jìn)行不穩(wěn)定的概率分析［3］。蒙特卡洛模擬方法是概率分析的基礎(chǔ)方法，在使用統(tǒng)計(jì)分布的基礎(chǔ)上確定控制滑動(dòng)面穩(wěn)定性的參數(shù)。因此，F(xiàn)S是通過數(shù)據(jù)分布中的隨機(jī)數(shù)值所確定的。

對(duì)于每個(gè)參數(shù)，本文使用不同的值，利用概率統(tǒng)計(jì)方法來(lái)確定安全系數(shù)，因?yàn)樵摲ń鉀Q了滑動(dòng)分析高度不確定性的問題。這就需要用以下參數(shù)來(lái)分析：邊坡高度、傾角和傾向、斷層傾角和傾向、內(nèi)部凝聚力及內(nèi)摩擦角。在分析過程中，密度和邊坡節(jié)理組的產(chǎn)狀定位為某一常數(shù)，而其他數(shù)據(jù)所服從的正態(tài)分布情況通過野外調(diào)查、實(shí)驗(yàn)室數(shù)據(jù)和基本統(tǒng)計(jì)參數(shù)來(lái)獲得。

考慮到在一些不連續(xù)面上沒有物質(zhì)的填充，同時(shí)為了明確在概率分析和敏感性分析法中無(wú)內(nèi)聚力的情況，本文統(tǒng)一假設(shè)最小的內(nèi)聚力值為零。不穩(wěn)定概率在輸入基本參數(shù)值后通過概率分析法求得。

2.4 靈敏度分析

如前所述，通過利用參數(shù)的范圍以及FS概率分布可以預(yù)測(cè)邊坡失效的可能性大小，而在靈敏度分析中，每一個(gè)參數(shù)的最大值和最小取值都是用來(lái)確定滑坡的臨界條件。不同參數(shù)之間的關(guān)系以及它們對(duì)FS的作用將進(jìn)行評(píng)估和討論。由于只進(jìn)行概率分析法并不能確定各參數(shù)之間的關(guān)系，因此，靈敏度分析是非常重要的環(huán)節(jié)。為解決這一問題，本文進(jìn)行了靈敏度分析，共3個(gè)步驟［4－5］：

1）在所有的參數(shù)之間，比如邊坡高度、邊坡和節(jié)理組的產(chǎn)狀、內(nèi)聚力及其他所有的參數(shù)中，只選取一個(gè)參數(shù)為變量，其他參數(shù)都選做常量。這樣在所有的統(tǒng)計(jì)因素都確定的情況下，F(xiàn)S可以通過這些因素而確定。這樣，分別更換不同的參數(shù)為變量，多次進(jìn)行這一運(yùn)算過程。

2）為研究節(jié)理組產(chǎn)狀、邊坡高度、邊坡產(chǎn)狀對(duì)FS的影響，通過保持內(nèi)摩擦角不變、改變內(nèi)聚力的值來(lái)進(jìn)行計(jì)算。每次分析都是其中的1個(gè)值在不斷變化，另外2個(gè)選為常數(shù)。

3）類似的為研究節(jié)理組產(chǎn)狀、邊坡高度、邊坡產(chǎn)狀對(duì)FS的影響，通過保持改變內(nèi)聚力的值不變，改變內(nèi)摩擦角的值來(lái)進(jìn)行計(jì)算。

以上3個(gè)步驟中，用到了野外和室內(nèi)試驗(yàn)中參數(shù)的最小值、最大值、平均值和標(biāo)準(zhǔn)偏差。同時(shí)，在分析過程中，邊坡和節(jié)理組的傾角認(rèn)為是常數(shù)。隨著參數(shù)值的變化，F(xiàn)S對(duì)不同參數(shù)的靈敏度不同，根據(jù)靈敏度的大小從而確定影響穩(wěn)定性最大的參數(shù)。靈敏度分析結(jié)果如圖1所示。

表1 節(jié)理組特征統(tǒng)計(jì)

2.5 結(jié)果與討論

概率分析有效地確定不同破壞方式概率的大小，根據(jù)這一分析，雖然平面滑動(dòng)和楔形體破壞的可能性都很高，但基于現(xiàn)場(chǎng)調(diào)查，平面滑動(dòng)的可能性更大一些。因此，圖1中的靈敏度分析是基于平面滑動(dòng)的。靈敏度分析是為確定有效參數(shù)對(duì)FS的影響程度。

基于研究結(jié)果，坡角在降低到84°時(shí)，F(xiàn)S開始有大的變化，當(dāng)降低到82°時(shí)，F(xiàn)S達(dá)到2.11（如圖1（a））。FS的變化值和邊坡高度之間呈線性關(guān)系，高度從20m降低到10m過程中，F(xiàn)S從0.53增加到0.91（如圖1（b））。

根據(jù)圖1（c），節(jié)理組的傾角在81°時(shí)，F(xiàn)S突然顯著增大。在角度小于81°時(shí)，F(xiàn)S僅僅有很小的變化。同樣，可以看到內(nèi)聚力和FS之間的關(guān)系，在內(nèi)聚力為0時(shí)，F(xiàn)S為0.12（如圖1（d））。

內(nèi)聚力的變化對(duì)FS的影響很小，增加內(nèi)聚力只能得到FS較小的增加。根據(jù)圖1（e）可以看出內(nèi)摩擦角在可控范圍內(nèi)變化，對(duì)FS的影響不明顯。而地下水對(duì)FS有相反的影響，當(dāng)含水率達(dá)到27%時(shí)，F(xiàn)S突然降到0。

圖1 平面滑動(dòng)的靈敏度分析

3 結(jié)語(yǔ)

巖石和不連續(xù)面的巖土參數(shù)視為研究邊坡穩(wěn)定性的有效因素，但由于這些參數(shù)的不確定性，為其選擇合適的量化值在巖石邊坡的穩(wěn)定性分析中成為重要的問題。本文基于靈敏度分析評(píng)估了不同的參數(shù)對(duì)邊坡穩(wěn)定性的影響，結(jié)果清楚表明：滑坡破壞存在著潛在的楔形破壞和平面破壞2種可能。通過使用敏感性分析有效地分析參數(shù)的影響。靈敏度分析的最終結(jié)果表明：內(nèi)摩擦角對(duì)FS值沒有顯著影響，但改變破壞面的傾角FS表現(xiàn)出較大的變化，同時(shí)當(dāng)改變含水率時(shí)，F(xiàn)S表現(xiàn)出的變化也很大。因此，在本文案例研究中，減少邊坡傾角和控制水泄漏是提高邊坡穩(wěn)定最有效的方法。

［1］張勇慧，李紅旭，盛謙，等.基于模糊綜合評(píng)判的公路巖質(zhì)邊坡穩(wěn)定性分級(jí)研究［J］.巖土力學(xué)，2010，10：3151－3156.

［2］張社榮，彭敏瑞，董紹堯.巖質(zhì)邊坡穩(wěn)定性分析方法及工程應(yīng)用［J］.中國(guó)農(nóng)村水利水電，2007（4）：94－96.

［3］彭小云，高德彬，畢選生，等.高陡邊坡穩(wěn)定性的概率分析［J］.長(zhǎng)安大學(xué)學(xué)報(bào)（地球科學(xué)版），2003（3）：67－70.

［4］佴磊，徐燕，代樹林.邊坡工程［M］.北京：科學(xué)出版社，2000.

［5］唐輝，陳潔，錢會(huì).飽和指數(shù)在水－巖作用研究中的應(yīng)用及其靈敏度分析［J］.水資源與水工程學(xué)報(bào)，2012，23（6）：180－183.