排土場散體物料抗剪強度的安全系數(shù)反演分析

張 默 汪 斌 周玉新 朱君星

(1.中鋼集團(tuán)馬鞍山礦山研究院有限公司，安徽 馬鞍山 243000;2.金屬礦山安全與健康國家重點實驗室，安徽 馬鞍山 243000；3.華唯金屬礦產(chǎn)資源高效循環(huán)利用國家工程研究中心有限公司，安徽 馬鞍山243000)

排土場散體物料抗剪強度的安全系數(shù)反演分析

張 默1,2,3汪 斌1,2,3周玉新1,2朱君星1,2

(1.中鋼集團(tuán)馬鞍山礦山研究院有限公司，安徽 馬鞍山 243000;2.金屬礦山安全與健康國家重點實驗室，安徽 馬鞍山 243000；3.華唯金屬礦產(chǎn)資源高效循環(huán)利用國家工程研究中心有限公司，安徽 馬鞍山243000)

在初步判斷排土場邊坡所處的穩(wěn)定狀態(tài)的基礎(chǔ)上設(shè)定一個安全系數(shù)，利用安全系數(shù)反演計算排土場散體物料的抗剪強度參數(shù)。以常用的Bishop圓弧法和余推力折線法作為排土場邊坡穩(wěn)定性計算的基本方法，構(gòu)建非線性方程組，迭代求解，同時反演出排土場散體物料的抗剪強度參數(shù)。針對西南地區(qū)某個大型露天礦排土場邊坡進(jìn)行了實際反演分析,得出：含第四系表土較多的細(xì)粒散體巖石，反分析黏聚力c比三軸試驗小0.026 MPa，而摩擦角φ則大3.94°；以粗粒土為主的堅硬輝長巖，反分析黏聚力c比三軸試驗大0.021 MPa，而摩擦角φ則小4.06°。雖然2種散體巖土反演分析與試驗結(jié)果都有差異，但是通過安全系數(shù)反演分析的結(jié)果，更能從整體上把握排土場散體物料的力學(xué)強度參數(shù)。該方法不像室內(nèi)試驗和現(xiàn)場試驗?zāi)菢邮芘磐辽Ⅲw物料不均勻和離散性的影響，同時獲得的抗剪強度參數(shù)可以作為后期動態(tài)排土規(guī)劃穩(wěn)定性分析的依據(jù)。

排土場 安全系數(shù) 反演分析 抗剪強度 排土規(guī)劃

排土場是礦山開采集中排棄剝離的廢石土，又稱廢石場。它是人為堆積的散體邊坡，處于欠固結(jié)狀態(tài)，物理力學(xué)參數(shù)很低，若沒有完善的排土工藝控制排土場有效受土措施，極易導(dǎo)致邊坡失穩(wěn)，致使排土場邊坡發(fā)生滑坡破壞；或因排土場不善的安全管理，在雨季期間沒有疏暢的截排水措施，致使排土場形成一個飽水體而發(fā)生泥石流現(xiàn)象[1]。因此而發(fā)生的排土場安全問題將牽涉一系列人身財產(chǎn)和環(huán)境、地質(zhì)災(zāi)害等問題[2-3]。

排土場的安全穩(wěn)定性主要通過散體物料的力學(xué)參數(shù)和有效的排土工藝參數(shù)(臺階高度、安全平臺寬度和排土速度等)確定，而力學(xué)參數(shù)除自身巖性決定外，還主要受排土工藝影響。目前，排土場安全穩(wěn)定性分析的力學(xué)參數(shù)可以通過室內(nèi)粗粒土三軸試驗和野外現(xiàn)場直剪試驗確定。室內(nèi)試驗可進(jìn)行各種巖性和粒度組成的對比試驗，比較靈活；野外現(xiàn)場試驗較室內(nèi)試驗更接近于實際情況，但受排土不均勻和各種現(xiàn)場因素的影響[4]，試驗數(shù)量有限，并且試驗困難和費用高。排土場散體物料在臺階上自然滾落堆積，呈現(xiàn)出明顯的下部大塊石含量多，上部細(xì)粒土含量多，致使散體物料的宏觀力學(xué)特性表現(xiàn)強烈的不均勻性和離散性[5]。雖然在穩(wěn)定性分析時，通常進(jìn)行室內(nèi)試驗，通過粒徑分部組成，重新配比粒徑級配代替原型粒徑級配的力學(xué)試驗[6]，在同一個臺階上部和下部采用不同的抗剪強度取值參數(shù)，但也很難從整體上準(zhǔn)確把握排土場散體物料的力學(xué)參數(shù)。根據(jù)邊坡穩(wěn)定性評價[7]，對于處在暫且穩(wěn)定的邊坡(邊坡處于欠穩(wěn)定狀態(tài))，安全系數(shù)可取[1.00，1.05]，進(jìn)行力學(xué)強度參數(shù)反分析[8-9]。本研究根據(jù)排土場所處的穩(wěn)定狀態(tài)，首先確定個安全系數(shù)，然后根據(jù)安全系數(shù)反分析強度參數(shù)，進(jìn)而可以從整體上較準(zhǔn)確地掌握排土場散體物料的抗剪強度參數(shù)。

1 滑坡反分析計算

排土場邊坡強度參數(shù)反演分析時，首先需要確定滑體穩(wěn)定狀態(tài)及滑動面位置等參數(shù)，排土場滑坡反演分析的滑面位置及穩(wěn)定狀態(tài)考慮如下。

(1)由于排土場初期散體沒有顯著的峰值強度，邊坡處于欠穩(wěn)定狀態(tài)，安全系數(shù)Fs儲備較低，本次取Fs=1.05進(jìn)行反演分析。

(2)由于滑體部分滑面未知，故按正分析方法，假設(shè)排土場散體有足夠的力學(xué)強度，求出臨界位置作為滑坡的滑面位置。

排土場滑坡模型包含有單一介質(zhì)散體巖土的圓弧破壞及散體巖體與地基土層組合破壞兩種形式[3]，對于不同類型的滑坡反分析是，應(yīng)使用適合不同類型破壞模式的穩(wěn)定性計算公式。

1.1 單一介質(zhì)圓弧破壞

單一介質(zhì)圓弧破壞采用Bishop穩(wěn)定性計算公式進(jìn)行分析：

(1)

式中，c、φ為黏結(jié)力和內(nèi)摩擦角；Wi為條塊質(zhì)量；bi為條塊寬度；Fs為安全系數(shù)，取1.05。

用2個滑坡剖面的參數(shù)可以組成一個非線性方程組，通過計算求解即可求得c、φ值。

1.2 組合滑動面破壞

滑面穿過2種不同介質(zhì)的組合滑面破壞共有4個參數(shù)。其反分析采用適合計算任意滑面安全的余推力法進(jìn)行分析[10]。

余推力反分析公式：

(2)

式中，Δα=αi-1-αi；Di為第i條塊的推力。

(2)式的定解必須有4個滑坡剖面的參數(shù)，對每一剖面可建立一個包含c、φ值的非線性方程，最后聯(lián)立成一個非線性方程組。

1.3 非線性方程組的求解

非線性方程組采用解線性代數(shù)非線性組的下降法，具體步驟如下(以2個參數(shù)Bishop法為例)。

目標(biāo)函數(shù)：

(3)

令上兩式分別為f1(c、φ)和f2(c、φ)，

定義函數(shù)

).

(4)

(1)取1組初值c0、φ0(不能為零)。

(2)設(shè)計算到m步；求解Fm(c、φ)。

(3)若Fm<ε(給定的任意小常數(shù)，例如10-5)則cm、φm即為所求值，否則計算c、φ的偏導(dǎo)數(shù)

(5)

式中，Δc=ε·cm，Δφ=ε·φm。

(4)計算

(6)

式中，

由此得到新的c、φ值，重復(fù)(2)、(3)、(4)步的計算。

2 案例分析

西南地區(qū)某大型露天鐵礦，設(shè)計年產(chǎn)礦石量630萬t/d，年采剝總量3 043萬t。排土場設(shè)計容積36 474萬m3，采用覆蓋式多臺階排土。自投產(chǎn)受土以來，由于不良的安全排土計劃，到處亂排、集中過快棄土，發(fā)生了數(shù)次滑坡、電機(jī)車箱傾覆和路基懸空等事故，造成了重大的經(jīng)濟(jì)損失，并威脅到礦山安全生產(chǎn)。該排土場因排棄散體物料中含有大量第四系表土和強風(fēng)化巖石，同時由于臺階較高而多次在發(fā)生排土場Ⅰ線發(fā)生內(nèi)部滑坡，在排土場II線沿軟弱 “昔格達(dá)”地層發(fā)生過底鼓而引起的滑坡。目前排土場多數(shù)臺階后緣出現(xiàn)了不同程度的張拉裂縫，已經(jīng)被停止受土，改善處理，可以初步判定邊坡處于欠穩(wěn)定狀態(tài)。在詳細(xì)實測了部分滑坡的滑體資料的基礎(chǔ)上，選取貼合實際的滑面，即分別選取排土場Ⅰ線和II線進(jìn)行散體物料內(nèi)部滑坡和沿軟弱地基的滑坡反分析計算。排土場平面見圖1所示。

圖1 排土場平面

分別按照上述的計算方法，對排土場Ⅰ線和II線進(jìn)行計算，計算剖面和滑弧位置見圖2所示。

通過上述的排土場滑面強度參數(shù)反演分析結(jié)果列于表1，排土場散體三軸試驗結(jié)果見表2。

由滑坡條件，排土場內(nèi)部滑坡反演分析強度為含有大量第四系表土的散體強度，與散體三軸試驗以細(xì)粒為主的試驗條件相近；沿地基軟弱層滑坡反演分析的散體強度為排棄以輝長巖為主的物料強度，與散體三軸試驗以粗粒為主的試驗條件相近：從2種方法所得結(jié)果看，含第四系表土較多的細(xì)粒散體巖石，反分析黏聚力c比三軸試驗小0.026 MPa，而摩擦角φ則大3.94°；以粗粒土為主的堅硬輝長巖反分析黏聚力c比三軸試驗大0.021 MPa，而摩擦角φ則小4.06°。雖然兩種散體巖土反演分析與試驗結(jié)果都有差異，但總的來看，兩者強度相差不大。所以通過安全數(shù)反演分析，更能從整體上把握排土場散體物料的力學(xué)強度參數(shù)。

圖2 計算剖面滑弧

表1 排土場滑坡反演分析計算結(jié)果

表2 排土場散體三軸試驗結(jié)果

3 結(jié) 論

(1)首先通過潛在的滑動面，根據(jù)排土場邊坡的穩(wěn)定狀態(tài)，設(shè)定已知的穩(wěn)定性安全系數(shù)，然后運用常規(guī)的條分法，構(gòu)建非線性方程組，迭代求解，同時求出排土場散體物料的抗剪強度參數(shù)。

(2)運用安全系數(shù)反演分析強度參數(shù)，可以從整體上掌握排土場的抗剪強度，不像室內(nèi)試驗和現(xiàn)場試驗?zāi)菢邮芘磐辽Ⅲw物料不均勻和離散性的影響。同時獲得的抗剪強度參數(shù)可以作為后期動態(tài)排土規(guī)劃穩(wěn)定性分析的依據(jù)。

[1] 中華人民共和國建設(shè)部.GB 50021—2001 巖土工程勘察規(guī)范[S].北京:中國建筑工業(yè)出版社,2009. Ministry of Construction P.R.China.GB50021-2001 Code for Investigation of Geotechnical Engineering[S].Beijing:China Building Industry Press，2009.

[2] 鄭穎人,陳祖煜,王恭先,等.邊坡與滑坡工程治理[M].2版.北京:人民交通出版社，2010. Zheng Yingren，Chen Zuyu，Wang Gongxian，et al.Engineering Treatment of Slope & Landslide[M].2nd Edition.Beijing:China Communications Press,2010.

[3] 李廣信,張丙印,于玉貞.巖土力學(xué)[M].2版.北京:清華大學(xué)出版社,2013. Li Guangxin，Zhang Bingyin，Yu Yuzhen.Soil Mechanics[M].2nd Edition.Beijing:Tsinghua University Press,2012.

[4] 張傳慶,馮夏庭,周 輝,等.隧洞圍巖收斂損失位移的求取方法及應(yīng)用[J].巖土力學(xué),2009，30(4):997-1003. Zhang Chuanqing，F(xiàn)eng Xiating，Zhou Hui，et al.Method of obtaining loss convergence displacement and its application to tunnel engineering[J].Rock and Soil Mechanics,2009,30(4):997-1003.

[5] 王 鵬，周傳波.大冶鐵礦高陡露天邊坡F9滑體力學(xué)參數(shù)反分析[J].金屬礦山,2008(1):36-38. Wang Peng，Zhou Chuanbo.Mechanics parameter back analysis of F9 sliding mass of high steep open-pit slope in Daye Iron Mine[J].Metal Mine,2008(1):36-38.

[6] 高 瑋，鄭穎人.采用快速遺傳算法進(jìn)行巖土工程反分析[J].巖土工程學(xué)報,2001,23(1):120-122. Gao Wei，Zheng Yingren.Back analysis in geotechnical engineering based on fast-convergent genetic algorithm[J].Chinese Journal of Geotechnical Engineering，2001,23(1):120-122.

[7] 徐漢斌,王 軍.反算法中滑坡穩(wěn)定系數(shù)的取值問題[J].四川地質(zhì)學(xué)報,1999,19(1)68-70. Xu Hanbin，Wang Jun.Selection of landslide stability factor in inverse calculation[J].Acta Geologica Sichuan，1999，19(1):68-70.

[8] 徐 軍，鄭穎人.基于響應(yīng)面方法的圍巖參數(shù)隨機(jī)反分析[J].巖土力學(xué),2001,22(2):167-170. Xu Jun，Zheng Yingren.Random back analysis of field geotechnical parameter by response surface method[J].Rock and Soil Mechanics,2001,22(2):167-170.

[9] 魏正友，華 躍，吳世雄.基于反分析法滑坡穩(wěn)定性分析[J].現(xiàn)代礦業(yè)，2012(5):61-62. Wei Zhengyou，Hua Yue，Wu Shixiong.Anti-based analysis of landslide stability analysis[J].Modern Mining,2012(5):61-62.

[10] 高德軍,徐衛(wèi)亞,郭其達(dá).長江三峽大石板滑坡計算參數(shù)反分析[J].河海大學(xué)學(xué)報:自然科學(xué)版,2006，34(1):74-78. Gao Dejun，Xu Weiya，Guo Qida.Back analysis of calculation parameters for Dashiban landslide at Three Gorges Reservoir on Yangtze River[J].Journal of Hohai University:Natural Science,2006,34(1):74-78.

(責(zé)任編輯 徐志宏)

Inversion Analysis of Safety Factor on Shear Strength of Waste Dump Bulk Materials

Zhang Mo1,2,3Wang Bin1,2,3Zhou Yuxin1,2Zhu Junxing1,2

(1.SinosteelMaanshanInstituteofMiningResearchCo.，Ltd.，Maanshan243000，China;2.StateKeyLaboratoryofSafetyandHealthforMetalMines，Maanshan243000，China；3.HuaweiNationalEngineeringResearchCenterofHighEfficientCyclicandUtilizationofMetalMineralResourcesCo.，Ltd.，Maanshan243000，China)

On the basis of preliminary judging the steady state of dump slope，a safety factor is set，the shear strengths of waste dump bulk materials are inversed with the safety factor.With the use of the common Bishop method and Residual thrust method，the nonlinear equation group is built to make the iterative solution.Meanwhile，the shear strengths of waste dump bulk materials are calculated out.Aiming at the actual inversion analysis on dump slope of a large open pit mine in southwest China，it is concluded that:for fine-grained granular rock with more quaternary topsoil，its cohesioncby inversion analysis is smaller than triaxial test by 0.026 MPa，but friction angleφis bigger by 3.94°；For hard gabbro with coarse grained soil as the main，its cohesioncby inversion analysis is bigger than that at triaxial test by 0.021 MPa，but its friction angleφis smaller by 4.06°.Although the two kinds of bulk materials have different inversion analysis and experiment results，the inversion analysis results based on the safety factor can better grasp the shear strength parameters of waste dump bulk materials as a whole.This method is superior to the indoor test and field test which are affected by the unevenness and discreteness of bulk materials.At the same time，the shear strength parameters obtained can be used as the basis of dynamic stability analysis in the late dump planning.

Waste dump，Safety factor，Inversion analysis，Shear strength，Dump plan

2015-02-16

“十二五”國家科技支撐計劃項目(編號：2012BAK09B05)。

張 默(1986—)，男，工程師。

TD824

A

1001-1250(2015)-05-171-04