順層巖質(zhì)高邊坡剪切滑動破壞的敏感性因素分析*

劉婷婷 張 軍

(武漢理工大學(xué)道路橋梁與結(jié)構(gòu)工程重點實驗室1) 武漢 430070)(武漢理工大學(xué)土木工程與建筑學(xué)院2) 武漢 430070) (中交第二公路勘察設(shè)計研究院有限公司3) 武漢 430056)

?

順層巖質(zhì)高邊坡剪切滑動破壞的敏感性因素分析*

劉婷婷1,2)張 軍3)

(武漢理工大學(xué)道路橋梁與結(jié)構(gòu)工程重點實驗室1)武漢 430070)(武漢理工大學(xué)土木工程與建筑學(xué)院2)武漢 430070) (中交第二公路勘察設(shè)計研究院有限公司3)武漢 430056)

依托山西永源花崗巖礦山邊坡，通過現(xiàn)場地質(zhì)調(diào)查、室內(nèi)試驗和極限平衡理論，研究了順層巖質(zhì)高邊坡沿結(jié)構(gòu)面組合發(fā)生滑動破壞的情況，基于簡化Bishop法，探討了潛在滑體的幾何特征(滑體高度、滑體寬度以及滑動面傾角)對邊坡穩(wěn)定性的影響，結(jié)果表明，敏感度大小順序為：滑動面傾角>滑體寬度>滑體高度，即滑動面傾角對邊坡穩(wěn)定性起主要控制作用，滑體寬度和滑體高度對邊坡穩(wěn)定性的影響相對較小.

順層邊坡；簡化Bishop法；滑坡；敏感度

0 引 言

順層巖質(zhì)邊坡是一種特殊類型的邊坡，在露天礦山、公路、鐵路等的建設(shè)中經(jīng)常遇到，并造成了重大損失[1].順層邊坡巖體往往沿著某個最不利的剪切滑動面向邊坡面方向滑動造成邊坡的失穩(wěn)破壞.分析順層巖質(zhì)邊坡的穩(wěn)定性問題，目前工程界主要應(yīng)用剛體極限平衡[2-3]和數(shù)值模擬[4-5]2種方法.在極限平衡方面，其中參數(shù)、值的選取是比較關(guān)鍵的因素，龔文惠等[6-7]通過算例計算分析了邊坡安全系數(shù)對、敏感程度.嚴(yán)明等[8]通過對一個順向坡實例的詳細(xì)調(diào)查發(fā)現(xiàn)，順向坡的滑移-彎曲破壞，在強烈彎曲-隆起之后，整體失穩(wěn)之前還存在一個中間狀態(tài)：碎裂-散體化，在對其形成條件進行初步分析的基礎(chǔ)上，提出了“四階段”的滑移-彎曲破壞地質(zhì)力學(xué)模式.梁專明等[9]以長沙商貿(mào)旅游職業(yè)技術(shù)學(xué)院工程邊坡為研究案例，探討了軟質(zhì)巖順層邊坡的非典型潰屈型滑坡機制.朱晗迓等[10]考慮水壓力和地震等多種荷載的共同作用，對層狀巖體的潰屈破壞模式進行了研究.Ebrahim等[11]利用離散元程序?qū)aroun順層滑坡進行反分析，發(fā)現(xiàn)結(jié)構(gòu)面參數(shù)對該邊坡的穩(wěn)定性起控制作用.針對軟硬巖互層順傾邊坡，夏開宗等[12]建立了順層軟硬巖互層邊坡穩(wěn)定性的分析模型，分析了層面強度、巖層厚度、巖層傾角以及水力作用等各項因素對邊坡穩(wěn)定性的影響.

上述學(xué)者主要討論了巖質(zhì)邊坡對參數(shù)c,φ值的敏感程度，鮮有學(xué)者討論潛在滑體的幾何特征對巖質(zhì)高邊坡穩(wěn)定性的影響.文中以山西永源花崗巖礦山邊坡為例，運用加拿大RocScience公司開發(fā)的SLIDE軟件分析了自重情況下潛在滑體的高度、厚度、滑動面傾角對邊坡穩(wěn)定性的影響，提出了運用SLIDE模擬結(jié)構(gòu)面的方法，為邊坡工程設(shè)計和施工提供指導(dǎo).

1 工程概況及計算參數(shù)

1.1 工程概況

山西永源石材有限公司渾源縣中莊鋪鄉(xiāng)人崖山花崗石礦(以下簡稱永源礦)1993年6月建成投產(chǎn)，批準(zhǔn)最高開采高程為2 070 m，最低開采標(biāo)高1 920 m，則會形成最高約150 m高的采場邊坡.現(xiàn)在處于生成后期，邊坡穩(wěn)定問題已經(jīng)成為制約露天礦生產(chǎn)的決定性因素.

根據(jù)現(xiàn)場地質(zhì)情況，選取西邊坡中露頭較好的位置進行了地表調(diào)查，各個調(diào)查點(區(qū)域)的位置見圖1；西邊坡由礦體(輝綠巖)和圍巖(片麻巖)兩部分組成，1～9號調(diào)查點出露的基巖大都是片麻巖，而10～14號調(diào)查點出露的基巖大都是輝綠巖.

圖1 地質(zhì)調(diào)查點分布圖

圖2 邊坡地質(zhì)模型

基于現(xiàn)場地質(zhì)調(diào)查，建立了邊坡的地質(zhì)模型，見圖2所示，并對該邊坡可能的主要破壞模式進行了初步分析，歸納為以下4種.

1) 潰曲破壞 潰屈破壞大多發(fā)育在插入式順層邊坡中，坡面沒有順層臨空面，邊坡高度較高，在重力或其他荷載作用下，邊坡坡腳應(yīng)力集中，導(dǎo)致坡腳巖層發(fā)生揉屈，進而發(fā)生潰屈破壞.永源礦采場西邊坡為中傾反向邊坡，發(fā)育有與坡向相同的陡傾結(jié)構(gòu)面，且此陡傾結(jié)構(gòu)面貫通性較好，采場邊坡高度較高，有可能在重力荷載作用下發(fā)生潰屈破壞.

2) 傾倒破壞 傾倒破壞大多發(fā)生在陡傾順層或陡傾反向巖質(zhì)邊坡中.永源礦采場西邊坡坡面上部分地方順坡向陡傾節(jié)理密集發(fā)育，容易發(fā)生傾倒破壞，見圖3.

圖3 邊坡局部坡面節(jié)理密集發(fā)育

3) 崩塌、墜落破壞 邊坡中產(chǎn)狀為340°∠82°的節(jié)理很發(fā)育，該組結(jié)構(gòu)面與坡面走向垂直.坡面上由該組結(jié)構(gòu)面與傾向坡內(nèi)的陡立結(jié)構(gòu)面(70°∠60°)、坡面切割出來的小型塊體容易發(fā)生崩塌或墜落.

4) 剪切滑移破壞 剪切破壞是巖質(zhì)邊坡常見的破壞模式，一般是沿著順層結(jié)構(gòu)面發(fā)生.對于永源礦西邊坡，一組產(chǎn)狀為70°∠60°陡傾順坡節(jié)理控制邊坡的穩(wěn)定性，節(jié)理平行排列，構(gòu)成滑移體的后緣邊界，而邊坡體上存在一組產(chǎn)狀為59°∠15°的結(jié)構(gòu)面，該組結(jié)構(gòu)面構(gòu)成了滑移控制面，即在邊坡中容易發(fā)生形如塊體ACDE(及ABIG)及AJF(及AIG)塊體的滑移型破壞.

上述4種潛在破壞模式，前3種破壞模式主要發(fā)生在邊坡表層或局部，潛在危害較小，而對于剪切滑移破壞，破壞深度和規(guī)模都很大，潛在危害嚴(yán)重，故需要對這種模式進行分析.

1.2 結(jié)構(gòu)面參數(shù)

結(jié)構(gòu)面是具有一定方向、延展較大而厚度較小的二維面狀地質(zhì)界面.它在巖體中的變化非常復(fù)雜.結(jié)構(gòu)面的存在，使巖體顯示構(gòu)造上的不連續(xù)性和不均質(zhì)性，巖體力學(xué)性質(zhì)與結(jié)構(gòu)面的特性密切相關(guān).結(jié)構(gòu)面最重要的力學(xué)性質(zhì)之一是抗剪強度.從結(jié)構(gòu)面的變形分析可以看出，結(jié)構(gòu)面在剪切過程中的力學(xué)機制比較復(fù)雜，構(gòu)成結(jié)構(gòu)面抗剪強度的因素是多方面的，大量試驗結(jié)果表明，結(jié)構(gòu)面抗剪強度一般可以用庫倫準(zhǔn)則表述

(1)

式中：c，φ分別為結(jié)構(gòu)面上的粘聚力和內(nèi)摩擦角；σn為作用在結(jié)構(gòu)面上的法向應(yīng)力.

式中：σ為剪切面上的正應(yīng)力，MPa；τ為剪切面上的切應(yīng)力，MPa；S為剪切面積，cm2；Pn為法向荷載，kg；Pτ為切向荷載，kg.

試驗后結(jié)構(gòu)面見圖4，根據(jù)試驗后破壞巖樣可以看出，試件的最終破壞按試件結(jié)構(gòu)面破壞.圖5為結(jié)構(gòu)面剪切試驗剪切應(yīng)力-剪切變形曲線.

圖4 結(jié)構(gòu)面試樣

圖5 結(jié)構(gòu)面剪切應(yīng)力-剪切變形曲線

以結(jié)構(gòu)面的直剪試驗結(jié)果和和已有的巖石力學(xué)參數(shù)為基礎(chǔ)，綜合以上巖體基本質(zhì)量指標(biāo)BQ計算、Hoek-Brown經(jīng)驗方程估算和折減系數(shù)法等3種取值方法，給出了建議采用的結(jié)構(gòu)面計算參數(shù)，見表1.

表1 邊坡計算參數(shù)

2 計算方法與模型

2.1 計算方法

采用簡化Bishop法計算巖質(zhì)邊坡折減滑動破壞，簡化Bishop法是英國著名土力學(xué)專家Bishop[13]提出的邊坡穩(wěn)定計算方法，被認(rèn)為是計算圓弧滑動面安全系數(shù)最好的方法,張香渝等[14]把簡化Bishop法推廣到非圓弧滑動面安全系數(shù)的求解上，對于非圓弧巖質(zhì)邊坡極限平衡分析，簡化Bishop法分析步驟如下：(1)假定只存在水平條間力,豎向條間力為零；(2)各土條滿足豎向力平衡條件；(3)滑體滿足總體力矩平衡條件.巖塊受力分析示意圖見圖6.

圖6 簡化Bishop條分法受力分析圖

可建立如下平衡條件.由巖塊i豎向力平衡可知

(2)

由巖塊i底面切向力平衡可知

(3)

潛在滑體整體力矩平衡條件

根據(jù)式(1)～(3)，可得

(5)

(6)

(7)

2.2 敏感度計算

敏感度分析是指通過計算一個自變量發(fā)生變化時所引起的因變量的改變率從而分析該自變量對因變量的貢獻.邊坡安全系數(shù)對第i個影響因素的敏感度計算公式為

(8)

2.3 計算模型

沿結(jié)構(gòu)面組合的剪切滑動破壞是巖質(zhì)高邊坡失穩(wěn)的一種主要形式，為了研究潛在滑體的幾何特征(高度h、厚度w以及滑動面傾角δ)對巖質(zhì)高邊坡穩(wěn)定性的影響，文中以山西永源花崗石礦邊坡為例，探討了自重情況下邊坡安全系數(shù)對潛在滑體的幾何特征的敏感度，計算中h=15，30，45，60，75，90m，w=20，40，60，80，100m；δ=0°，15°，30°，共計算了6×5×3=90種情況.計算模型見圖7，計算參數(shù)見表1.

圖7 算例邊坡幾何形狀

基于SLIDE的特點，文中提出了一種用SLIDE計算潛在滑體沿結(jié)構(gòu)面破壞的安全系數(shù)的方法，解決了SLIDE無法計算巖體沿弱面發(fā)生破壞的問題.具體思路是：在對巖體參數(shù)賦值時，巖石的重度不變，但是巖石強度參數(shù)c，φ賦的是結(jié)構(gòu)面的強度參數(shù)值.

3 結(jié)果與討論

圖8～10為不同δ下的安全系數(shù)與滑體高度關(guān)系的曲線.由圖8～10可知，當(dāng)w和δ一定時，隨著h的增大，安全系數(shù)Fos非線性減小，但是變化率即dFos/dh隨h的增大而不斷減小，說明滑體高度較小時，安全系數(shù)隨之增大而顯著減小，但當(dāng)h增大到一定程度后，安全系數(shù)隨滑體高度的增大變化不明顯.

圖8 安全系數(shù)與滑體高度關(guān)系曲線(δ=0°)

圖9 安全系數(shù)與滑體高度關(guān)系曲線(δ=15°)

圖10 安全系數(shù)與滑體高度關(guān)系曲線(δ=30°)

圖11 安全系數(shù)對滑體高度的敏感度與滑體寬度關(guān)系曲線

由圖8～10可知，當(dāng)δ≤15°時，隨滑體寬度的增大，邊坡安全系數(shù)不斷增大，但是當(dāng)δ>15°，隨滑體寬度的增大Fos不斷減小.上述結(jié)論說明滑體寬度對邊坡安全系數(shù)的影響受滑動面傾角影響較大，δ不同，F(xiàn)os-w關(guān)系曲線走勢可能不同.

圖12 安全系數(shù)對滑體寬度的敏感度與滑體寬度關(guān)系曲線

由圖8～10可知，在w和h一定的情況下，滑動面傾角越大，邊坡安全系數(shù)越小，例如當(dāng)δ=0°時，F(xiàn)os最大達(dá)25.58，相同情況下，δ=30°時，F(xiàn)os僅為1.97，兩者竟相差13倍，說明滑動面傾角對邊坡安全系數(shù)的影響非常顯著.

圖13 安全系數(shù)對滑動面傾角的敏感度與滑體高度關(guān)系曲線

4 結(jié) 論

1) 潛在滑體的幾何特征中，敏感度大小順序為：滑動面傾角>滑體寬度>滑體高度，說明滑動面傾角對邊坡穩(wěn)定性起主要控制作用，滑體寬度和滑體高度對邊坡穩(wěn)定性的影響相對較小.

2) 邊坡穩(wěn)定性隨滑體高度的增大非線性減小，并且邊坡安全系數(shù)對滑體高度的敏感度受滑體寬度的影響不大，但是滑動面傾角對其影響較明顯，滑動面傾角越大，安全系數(shù)對滑體高度的敏感度越小，反之越大.

3) 滑體寬度對邊坡穩(wěn)定性的影響受滑動面傾角影響較大，當(dāng)滑動面傾角≤15°時，滑體寬度越大，邊坡穩(wěn)定性越高，但是當(dāng)滑動面傾角>15°時，滑體寬度越大，邊坡穩(wěn)定性越低.

4) 滑動面傾角對邊坡穩(wěn)定性的影響最大，相同情況下，當(dāng)滑動面傾角由0°增大到30°時，邊坡安全系數(shù)由25.58驟降到1.97.安全系數(shù)對滑動面傾角的敏感度隨滑體高度的增大而非線性減小，而隨滑體寬度增加而增大.

[1]李維光，張繼春.爆破振動作用下順層巖質(zhì)邊坡穩(wěn)定性分析[J].爆炸與沖擊，2007，27(5)：426-430.

[2]孟剛，余志雄.順層巖質(zhì)邊坡穩(wěn)定性影響因素分析[J].巖土力學(xué)，2005，26(1)：52-56.

[3]李安洪，周德培，馮君.順層巖質(zhì)路塹邊坡破壞模式及設(shè)計對策[J].巖石力學(xué)與工程學(xué)報，2010，28(1)：2915-921.

[4]盧敦華，曲艷偉，何忠明，等.結(jié)構(gòu)面強度參數(shù)對層狀邊坡穩(wěn)定性影響的三維分析[J].中南大學(xué)學(xué)報(自然科學(xué)版)，2011，42(1)：147-151.

[5]倪衛(wèi)達(dá)，唐輝明，劉曉，等.考慮結(jié)構(gòu)面震動劣化的巖質(zhì)邊坡動力穩(wěn)定分析[J].巖石力學(xué)與工程學(xué)報，2013，32(3)：492-500.

[6]姚雪貴.面板堆石壩結(jié)構(gòu)面參數(shù)的敏感性分析[J].山西建筑，2010，36(35)：356-357.

[7]陳佳偉，齊魯，翟開源，等.巖質(zhì)高邊坡結(jié)構(gòu)面參數(shù)的敏感性分析[J].四川水力發(fā)電，2011，30(6)：122-123.

[8]嚴(yán)明，陳劍平，黃潤秋，等.巖質(zhì)邊坡滑移-彎曲破壞中間狀態(tài)的工程地質(zhì)分析[J].水利水電技術(shù)，2004，36(11):41-44.

[9]梁專明，龔道民，賀友鋒，等.軟質(zhì)巖順層邊坡非典型潰屈型滑坡機制初探[J].自然災(zāi)害學(xué)報，2010，19(9):31-39.

[10]朱晗迓，馬美玲，尚岳全，等.順傾向?qū)訝顜r質(zhì)邊坡潰屈破壞分析[J].浙江大學(xué)學(xué)報，2004，38(9):1144-1149.

[11]EBRAHIM F S, SAEED H. Slope stability assessment using both empirical and numerical methods: a case study[J]. Environ Geol,2015,76:13-25.

[12]夏開宗，陳從新，魯祖德，等.軟硬巖互層邊坡穩(wěn)定性的敏感性因素分析[J].武漢理工大學(xué)學(xué)報(交通科學(xué)與工程版)，2013，37(4)：729-733.

[13]BISHOP A W. The use of the slip circle in the stability analysis of earth slopes[J]. Geotechnique，1955，5(1)：7-17.

[14]張香渝，鄭穎人.簡化Bishop法的擴展及其在非圓弧滑面中的應(yīng)用[J].巖土力學(xué)，2004，25(6)：927-930.Analysis of Sensitivity Factors to Sliding-shear Failure of High Dipping Rock Slope

LIU Tingting1,2)ZHANG Jun3)

(KeyLaboratoryofRoadwayBridgeandStructureEngineering,WuhanUniversityofTechnology,Wuhan430070,China)1)(SchoolofCivilEngineeringandArchitecture,WuhanUniversityofTechnology,Wuhan430070,China)2)(CCCCSecondHighwayConsultantsCo.Ltd,Wuhan430070,China)3)

Selecting Yongyuan granite mine slope as the site for analysis, the study on the mechanism of sliding-shear failure along the combination of steep and gently inclined joints is conducted through field investigation, laboratory tests and limit equilibrium theory. Based on simplified bishop method, the influence of potential landslide geometric characteristics(landslide height, landslide width and sliding surface inclination) to slope stability is explored for sliding-shear failure of high rock slope that along the combinations of joints and various fracture zones. The sensitivity level of each geometric characteristic is established. The orders of sensitivity from large to small are the sliding surface inclination, landslide width and land slide height. It means that sliding surface inclination performs a control role on the slope stability and the influence of landslide height and landslide width is relatively weak.

dipping rock slope; simplified bishop method; landslide; sensitivity

2016-08-24

*中國博士后面上基金資助(2015M582289)

TU457 doi:10.3963/j.issn.2095-3844.2016.05.007

劉婷婷(1986- )：女，博士后，主要研究領(lǐng)域為應(yīng)力波傳播和邊坡穩(wěn)定性分析