翟文龍 周漢民

(1.北京礦冶研究總院；2.金屬礦山智能開采技術(shù)北京市重點實驗室)

簡化Bishop法和余推力法在高臺階排土場穩(wěn)定性計算中的應(yīng)用

翟文龍1，2周漢民1，2

(1.北京礦冶研究總院；2.金屬礦山智能開采技術(shù)北京市重點實驗室)

邊坡穩(wěn)定性計算的方法有很多種，但適合排土場的尚未有統(tǒng)一的方法。介紹了簡化Bishop法和余推力法的基本理論，分析了2種計算方法的不同假設(shè)條件，以某高臺階排土場為例，選取了5個關(guān)鍵剖面分別采用簡化Bishop法和余推力法進(jìn)行了正常工況和地震工況穩(wěn)定性計算分析。結(jié)果顯示，正常工況下，計算安全系數(shù)值均大于1.4；地震作用下，安全系數(shù)均產(chǎn)生了相同趨勢的變化，安全儲備有限，需加強(qiáng)排土場防震和抗震措施。同種工況下同一剖面2種方法的計算結(jié)果相近，且滑弧的位置也幾近重合。計算分析結(jié)果可為排土場后期生產(chǎn)管理提供一定的參考。

排土場 余推力法 簡化Bishop法 穩(wěn)定性分析

排土場也稱廢石場，是指礦山采礦排棄物集中堆放的場所[1]。我國露天開采占用土地面積為露天采場本身面積的2～11倍，其中排土場面積占40%～60%[2]。我國排土場出現(xiàn)過很多事故，造成了巨大的生命和財產(chǎn)損失，其主要原因是我國排土場規(guī)范出臺較晚，且現(xiàn)行標(biāo)準(zhǔn)只針對有色金屬礦山排土場，對于排土場規(guī)模遠(yuǎn)大于有色金屬礦山的黑色金屬礦山，還沒有相應(yīng)的規(guī)范[3]。有色金屬礦山排土場規(guī)范中規(guī)定對土質(zhì)邊坡和較大規(guī)模的巖質(zhì)邊坡，采用瑞典法進(jìn)行計算。產(chǎn)生平面滑動的，直接利用滑面上的抗滑力除以下滑力；對于折線滑面，采用傳遞系數(shù)法。且規(guī)范說明中指出，計算公式參考現(xiàn)行國家標(biāo)準(zhǔn)規(guī)范[4]?！督饘俜墙饘俚V山排土場安全生產(chǎn)規(guī)則》(AQ 2005—2005)中規(guī)定排土場正常級需滿足用余推力法計算的安全系數(shù)大于1.15[5]。當(dāng)滑動面是圓弧時，余推力法和簡化Bishop法的計算安全系數(shù)是非常接近的，而且兩者所搜索到的臨界滑弧位置幾乎重合[6]。為進(jìn)一步論證2種方法在排土場中的適用性，本文以某黑色金屬礦山高臺階排土場為例，從簡化Bishop法和余推力法基本原理出發(fā),分析2種計算方法的不同假設(shè)條件，對排土場邊坡穩(wěn)定性進(jìn)行計算，對比分析計算結(jié)果，為類似排土場穩(wěn)定性計算方法選擇提供借鑒。

1 簡化Bishop法的基本理論

簡化Bishop法是國際公認(rèn)的一種比較精確的計算方法[6-7]，假定條塊間力的合力為法向水平力，不考慮條間切向力，計算分析簡圖見圖1。

圖1 簡化Bishop法條塊受力圖

經(jīng)推導(dǎo)，滑面的穩(wěn)定系數(shù)可表示為

(1)

式中，c為條塊底面黏聚力，kN；li為第i條塊底面長度，cm；Wi為第i條塊的自重力，kN；ui為滲透水壓力，kN；αi為第i條塊底面傾角，(°)；QAi為地震等產(chǎn)生的水平力，kN；φ為條塊底面摩擦角，(°)；ei、R分別為第i條塊重心到滑弧圓心的垂直距離和滑弧的半徑，cm。

在Mohr-Coulomb準(zhǔn)則中，c、φ為常量，其他符號意義同上。通過反復(fù)迭代，直至前后相鄰兩次算出的穩(wěn)定系數(shù)非常接近(或滿足預(yù)設(shè)精度要求)為止，此時的Fs即為特定滑面的穩(wěn)定系數(shù)。

2 余推力法的基本理論

余推力法又稱平衡推力傳遞法[7]，在假設(shè)第一條塊合力方向與上一條塊底面平行的條件下，求得極限平衡狀態(tài)下的安全系數(shù)，是工程中一種實用方法。計算分析簡圖見圖2。

圖2 余推力法條塊受力圖

根據(jù)力系分析，條塊界面受力為

Pi=Wisinαi+KcWicosαi+ψiPi-1-[Cili

+(Wicosαi-Ui-KcWisinαi)tanψi]/Fs，

(3)

ψi=cos(αi-1-αi)-sin(αi-1-αi)tanφi/Fs，

(4)

式中，Pi為第i條塊作用力，kN；Kc為地震作用力系數(shù);Pi-1為第i-1條塊作用力，kN；Ci為第i條塊底面黏聚力，kN；Ui為第i條塊水壓力，kN；αi-1為第i-1條塊底面傾角，(°)；φi為第i條塊底面摩擦角，(°);其他符號意義同上。

在計算過程中，若出現(xiàn)Pi<0，對下一條塊取Pi=0，因為模型假定條塊間不承受相互的拉力作用。

無外力作用時，P0=PN=0，因而計算時先給不同的Fs值，自上而下計算Pi，使PN=0的Fs為所求的安全系數(shù)。2種計算方法的區(qū)別見表1。

3 工程實例

3.1 排土場概況

某排土場原始地形為中低山區(qū)溝谷地貌，溝谷呈Y形分布，近于北東—南西走向；溝谷斷面呈V形，山體呈渾圓狀，坡度一般為10°～30°，山體大部分為殘坡積成因的碎石土及粉質(zhì)黏土，局部為遷西群三屯營組片麻巖出露。

經(jīng)鉆探揭露及排土場周邊工程地質(zhì)調(diào)查，該排土場及周圍100 m范圍內(nèi)，不存在斷裂構(gòu)造、滑坡、崩塌、泥石流、塌陷、液化土層等不良地質(zhì)作用。

表1 2種計算方法主要區(qū)別

該排土場基底坡度相對較緩，為平地型排土場，若排土場邊坡發(fā)生滑坡破壞,多為沿排土場內(nèi)部發(fā)生圓弧形滑動。排土場堆排物料共計1.9億m3，堆置總高度為250 m,屬于一等排土場。

3.2 計算剖面選擇

該排土場共分5個臺階堆排，臺階標(biāo)高分別為160,195,240,285,330 m，各排土平臺自下而上命名為第一～第五臺階；其中，第一臺階單臺階高35～80 m，第二臺階高35 m，第三～第五臺階高45 m，平臺寬30 m；臺階坡面角為自然安息角37.5°。

根據(jù)該排土場現(xiàn)場實際情況，針對危險區(qū)域(下游存在構(gòu)筑物)的分布共選擇5個剖面，建立排土場穩(wěn)定性分析概化模型，見圖3。

圖3 排土場穩(wěn)定性計算典型剖面

3.3 計算參數(shù)的選取

根據(jù)地質(zhì)調(diào)查及鉆探，該排土場地層主要由人工填土排土料(Qr)、第四系坡積土(Qel+dl)、黑云混合片麻巖(Ars)組成，人工填土區(qū)域下覆蓋的第四系殘坡積土較薄。具體排土料及地基土層物理力學(xué)參數(shù)見表2。

表2 巖體力學(xué)參數(shù)

根據(jù)《建筑抗震設(shè)計規(guī)范》(GB 50011—2010)，該場區(qū)地震設(shè)防烈度為7度，取地震峰值加速度0.15 g。

3.4 計算結(jié)果及分析

采用2種計算方法分別對該排土場5個關(guān)鍵剖面進(jìn)行正常工況和地震工況穩(wěn)定性計算，計算結(jié)果見表3。

表3 計算結(jié)果統(tǒng)計

該排土場邊坡均為沿排土場內(nèi)部的圓弧形滑動，符合該排土場易發(fā)生淺層圓弧形滑動的規(guī)律。5個關(guān)鍵剖面通過簡化Bishop法計算的正常工況下安全系數(shù)范圍在1.4～1.464，地震工況下安全系數(shù)范圍在1.184～1.235；通過余推力法計算的正常工況下安全系數(shù)范圍在1.403～1.469，地震工況下安全系數(shù)范圍在1.189～1.237。2種計算方法對5個剖面計算的安全系數(shù)變化規(guī)律是相同的，地震作用下，安全系數(shù)均有一定幅度的降低。簡化Bishop法計算的地震工況下相比正常工況安全系數(shù)降低幅度在0.216～0.229；余推力法計算的地震工況下相比正常工況安全系數(shù)降低幅度在0.214～0.232。

4 結(jié) 語

(1)穩(wěn)定性計算過程中，簡化Bishop法和余推力法主要區(qū)別在于對條塊間相互作用力作了不同的處理，前者假定了條塊間作用力的方向，后者假定了條塊間水平作用力的位置。

(2)余推力法適用于任意滑動面，簡化Bishop法是計算圓弧形破壞常用的計算方法。針對該黑色金屬礦山高臺階排土場5個關(guān)鍵剖面正常工況和地震工況分別采用2種方法計算的穩(wěn)定性安全系數(shù)結(jié)果相近，安全系數(shù)相差均在小數(shù)點后第三位。當(dāng)確定排土場邊坡滑動面是圓弧滑動時，可選擇簡化Bishop法與余推力法同時計算復(fù)核，二者不但計算安全系數(shù)相近，且臨界滑弧幾近重合。

(3)2種方法的計算結(jié)果均能較好與該排土場現(xiàn)場實際情況吻合。排土場處于穩(wěn)定狀態(tài)，地震工況下安全儲備有限，建議加強(qiáng)排土場防震、抗震措施及監(jiān)測工作，嚴(yán)禁盜挖坡角，防止破壞排土場整體穩(wěn)定性。

[1] 翟文龍，周漢民,陳 斌，等.降雨對某山坡型高臺階排土場的影響[J].有色金屬：礦山部分，2013，65(4)：87-90.

[2] 李根福.土地復(fù)墾知識[M].北京：冶金工業(yè)出版社，1991.

[3] 李小春，任 偉，王少泉，等.論金屬礦山排土場設(shè)計規(guī)范中邊坡極限平衡計算方法的選取[J].巖石力學(xué)與工程學(xué)報，2011，30(s2)；4136-4142.

[4] 中國有色金屬工業(yè)協(xié)會.GB 50421—2007 有色金屬礦山排土場設(shè)計規(guī)范[S].北京：中國計劃出版社，2007.

[5] 中鋼集團(tuán)馬鞍山礦山研究院有限公司.AQ 2005—2005 金屬非金屬礦山排土場安全生產(chǎn)規(guī)則[S].北京：中國計劃出版社，2005.

[6] 蔡美峰，何滿朝，劉東燕.巖石力學(xué)與工程[M].北京：北京科學(xué)出版社，2004.

[7] 朱大勇，錢七虎，周早生，等.基于余推力法的邊坡臨界滑動場[J].巖石力學(xué)與工程學(xué)報，1999，18(6)：667-670.

2015-03-13)

*國家科技支撐計劃(編號：NO.2012BAB08B06，NO.2013BAE04B01)。

翟文龍(1984—)，男，工程師，碩士，100160 北京市大興區(qū)北興路(東段)22號。