辛保泉 譚欽文 謝羽佳 王永強 萬 露

(西南科技大學(xué)環(huán)境與資源學(xué)院，四川 綿陽 621010)

巖性劣化對露天礦山邊坡穩(wěn)定性的響應(yīng)規(guī)律

辛保泉 譚欽文 謝羽佳 王永強 萬 露

(西南科技大學(xué)環(huán)境與資源學(xué)院，四川 綿陽 621010)

針對巖性條件劣化導(dǎo)致邊坡穩(wěn)定性下降的問題,以某露天礦山邊坡為研究對象，運用SIR-20型地質(zhì)雷達探測礦山滑坡體內(nèi)部缺陷,并結(jié)合Slide軟件模擬分析了不同巖性條件下邊坡穩(wěn)定性的變化規(guī)律。結(jié)果表明：軟弱夾層、降雨入滲和爆破動荷載等內(nèi)外因共同作用造成邊坡巖性劣化，進而導(dǎo)致邊坡后期出現(xiàn)滑移失穩(wěn)；黏聚力c值降低是巖性劣化過程中的主要力學(xué)表現(xiàn)，安全系數(shù)與c值呈部分拋物線關(guān)系，且c值降為50 kPa是邊坡特殊工況的失穩(wěn)臨界值。在綜合分析礦山運行特征和巖性劣化條件基礎(chǔ)上，提出通過減小上部第四系臺階坡面角和最終邊坡角提高其永久穩(wěn)定性的建議，安全系數(shù)經(jīng)驗證符合要求。研究結(jié)果可為礦山后期安全開采和邊坡控制參數(shù)取值提供重要依據(jù)。

露天礦山 邊坡穩(wěn)定性 巖性劣化 數(shù)值模擬 黏聚力 安全系數(shù)

露天開采以其獨有的優(yōu)勢在我國礦山開采中占有較大比重，其中露采鐵礦占鐵礦總開采量的85%左右[1]。而臺階邊坡作為露天礦山重大的安全性工程，有近40%存在不同程度的穩(wěn)定性問題[2]，邊坡穩(wěn)定性研究已成為露天礦安全生產(chǎn)的關(guān)鍵技術(shù)問題。

邊坡巖體的變形、滑動和破壞是一個動態(tài)發(fā)展過程，科學(xué)合理的設(shè)計是保證露天邊坡安全穩(wěn)定的基礎(chǔ)。目前，露天礦邊坡設(shè)計中的穩(wěn)定性分析方法有很多[3]，但幾乎都是基于礦山初期巖性條件進行的[4]，缺乏對礦山不同工況和后期巖性參數(shù)動態(tài)變化的考慮，導(dǎo)致實際邊坡按照設(shè)計尺寸開采時，隨著巖性條件劣化，邊坡穩(wěn)定性迅速下降進而出現(xiàn)失穩(wěn)滑移?；诖耍狙芯恳阅陈短斓V開采邊坡為研究對象，采用地質(zhì)雷達技術(shù)探測礦山臺階邊坡內(nèi)部地質(zhì)缺陷[5]，分析推演巖性劣化原因，并結(jié)合有限元分析軟件Slide對計算開采邊坡穩(wěn)定性，探明巖性條件劣化對邊坡穩(wěn)定性的動態(tài)響應(yīng)規(guī)律，經(jīng)多次校核后提出符合該露天礦山開采邊坡穩(wěn)定性的建議尺寸。

1 某露天礦邊坡概況

1.1 滑坡分布及裂隙發(fā)育

根據(jù)現(xiàn)場調(diào)查和測繪，該采場邊坡屬典型破裂結(jié)構(gòu)巖質(zhì)邊坡，且存在具有貫穿趨勢的拉裂隙，呈圈椅狀，滑坡體縱向長度約151 m，橫向?qū)挾?34～252 m。當(dāng)前滑坡體已發(fā)生明顯坍塌及位移，邊坡全貌及滑坡分布如圖1、圖2所示。

圖1 邊坡全貌

圖2 采場邊坡滑坡分布

該開采臺階邊坡總體呈現(xiàn)2大類型：

(1)局部臺階坍塌，最終臺階結(jié)構(gòu)受損。在采場左側(cè)中部臺階由于巖層產(chǎn)狀、降雨和風(fēng)化等因素影響，已出現(xiàn)較大面積垮塌，嚴(yán)重影響終了邊坡局部結(jié)構(gòu)完整性和穩(wěn)定性，威脅礦山開采和運輸安全。

(2)中上部整體滑移。在采場右側(cè)中上部出現(xiàn)“圈椅狀”拉裂縫，且其走勢沿典型巖性分界面和節(jié)理裂隙呈現(xiàn)出多臺階貫通趨勢。

1.2 滑坡發(fā)展趨勢及危害性分析

滑坡現(xiàn)狀為欠穩(wěn)定—不穩(wěn)定狀態(tài)，局部存在緩慢移動變形。暴雨作用使滑坡體后緣裂隙急速發(fā)展，若坡面暴雨積水進入場地后緣，將給場地的礦山開采帶來潛在致命的威脅。該滑坡體后部現(xiàn)已形成一級明顯滑坡，在強降雨、爆破震動等因素影響下，下部陡坎處也已形成1個潛在剪出滑移面。整個滑坡體呈現(xiàn)整體滑坡勢態(tài)。

此外，采礦活動對礦體上的滑坡體有較大的影響，如遇大開挖、大爆破、雨水浸潤，或是地震等誘發(fā)因素極易導(dǎo)致坡體滑動[6]。

2 滑坡體雷達探測結(jié)果分析

地質(zhì)雷達可探測、揭示臺階邊坡內(nèi)部地質(zhì)構(gòu)造和潛在缺陷。本研究運用SIR-20型地質(zhì)雷達檢測礦山滑坡體及裂隙情況[7]，探明該區(qū)域出現(xiàn)臺階嚴(yán)重垮塌及滑坡的主要原因。分別在采區(qū)滑坡頂部、中部和底部布置3條測線進行探測，并選取裂隙較為典型的測線Ⅰ、Ⅱ(見圖2中測線Ⅰ、測線Ⅱ)分析，結(jié)果如圖3、圖4所示。圖中橫軸為水平位置，縱軸表示探測深度。

圖3 測線Ⅰ掃描圖像

圖4 測線Ⅱ掃描圖像

圖3中，地面12 m以下雷達掃描圖信號急劇衰減至缺失，說明雷達電磁波衰減、吸收強烈，表示土體含水率較高；在12 m以上，圖像呈現(xiàn)反射異常，說明反射區(qū)發(fā)生了振幅反向，由于電磁波從介電常數(shù)小的介質(zhì)進入介電常數(shù)大的介質(zhì)造成。分析礦山滑坡體觀測情況可知，所測測線表面裂縫區(qū)域以下存在深約10 m裂隙或斷層，即圖3中標(biāo)有同相軸錯斷區(qū)域所示，說明局部土體松散或存在貫穿裂隙。

圖4中，測線區(qū)域深度9 m以下的雷達電磁波衰減、吸收較強烈，表明該區(qū)域土體松散，含水率較高；在距測線起始點水平位置為98～100 m的地下，存在雷達云圖同相軸缺失現(xiàn)象，表明該區(qū)域存在裂隙、土體不密實或局部空洞。這對滑坡體的穩(wěn)定十分不利。

按照同樣方法解析各測線不同里程的地層、巖性構(gòu)造，結(jié)果見表1。

表1 采區(qū)雷達探測結(jié)果匯總

根據(jù)雷達探測結(jié)果可知該區(qū)域裂隙發(fā)育，深度在1.5～9 m。采區(qū)地面以上3 m基本為松散土體，地下約10 m巖土體出現(xiàn)疏松。分析上述雷達探測結(jié)果，可將巖性劣化的原因歸為：

(1)關(guān)鍵內(nèi)因。坡面內(nèi)部存在軟弱夾層，滑坡物質(zhì)的互層、破碎給滑坡的滑移提供了物質(zhì)基礎(chǔ)。

(2)外在誘因。礦山開采改變了坡面坡面形態(tài)和巖土支撐，導(dǎo)致應(yīng)力重新分布[8]；降雨及爆破震動等外因致使滑帶土飽和，降低了抗滑能力，為滑坡體的滑動提供了外部條件。

探測結(jié)果表明，地質(zhì)構(gòu)造、降雨入滲和爆破震動等因素的共同影響造成邊坡巖性劣化嚴(yán)重，定性地解釋了該露天礦現(xiàn)狀邊坡穩(wěn)定性較差的原因。

3 邊坡穩(wěn)定性數(shù)值模擬分析

3.1 計算剖面

選取垮塌區(qū)圈椅狀滑坡體中線(圖2中A-A剖面)為典型剖面進行分析，并與其同期進行的地質(zhì)雷達結(jié)果互相補充和印證。

現(xiàn)狀邊坡尺寸(即原設(shè)計尺寸)如圖5所示：臺階高10 m；安全平臺寬4 m；清掃運輸平臺寬10 m；每隔2個安全平臺設(shè)置1個清掃平臺；最終臺階坡面角上部第四系45°～50°，下部60°，最終邊坡角42°。構(gòu)成本開采剖面的礦巖主要為片理化蝕變基性火山巖，其物理力學(xué)性質(zhì)如表2所示。

圖5 設(shè)計尺寸邊坡剖面簡圖

表2 礦巖物理力學(xué)性質(zhì)

注：RQD值表示礦巖工程質(zhì)量系數(shù)。

從表2知，近6 a邊坡巖體破碎現(xiàn)象突出，吸水性提高，質(zhì)量下降明顯，巖性條件劣化嚴(yán)重。其中，表征礦巖抗剪強度的黏聚力c值降幅最大(如圖6所示)，自2009年至今已降低了75%，平均25 kPa/a，目前已降為50 kPa并趨于穩(wěn)定。

圖6 黏聚力-時間變化曲線

3.2 計算方法

該露天礦山巖體松散且破壞面復(fù)雜，所處區(qū)域降雨充沛，宜選用圓弧形滑動面并根據(jù)簡化Bishop法計算分析邊坡穩(wěn)定性[9]，同時用簡化Janbu法校核。簡化Bishop法是假設(shè)每個條塊側(cè)面上的作用力是水平方向的，即條塊之間無摩擦，作用在第i條塊上的力見圖7所示。邊坡穩(wěn)定程度采用邊坡安全系數(shù)表示，即邊坡安全系數(shù)F定義為抗滑力矩和滑動力矩之比或抗剪強度與剪應(yīng)力之比[10]。

圖7 簡化Bishop法穩(wěn)定系數(shù)計算

根據(jù)垂直方向的分力為零，并考慮孔隙壓力比ru，通過求解方程得出穩(wěn)定性系數(shù)計算公式為

(1)

其中，c為黏聚力，kPa；φ為內(nèi)摩擦角，(°)；γ為土的重度,kN/m3;bi為條塊寬度，m；R為滑弧半徑，m；Qi為地震慣性力(水平方向)，kN；Wi為條塊的巖土重力，kN；θi為條塊滑面的傾角，(°)。

3.3 計算結(jié)果

根據(jù)邊坡設(shè)計尺寸和實測巖體力學(xué)參數(shù)，采用(Slide軟件的)隨機搜索方法，分析和確定圓弧的最危險滑動面，并計算安全系數(shù)。計算時，礦巖密度ρ取2.70 t/m3，滲透系數(shù)取2×10-6cm/s，設(shè)計基本地震加速度值為0.10g，計算所用浸潤線位置由現(xiàn)場勘查和地質(zhì)雷達探測結(jié)果所得，平均埋深為7.169 m，坡面下部最小埋深1.912 m。計算結(jié)果見表3。

表3 現(xiàn)狀邊坡穩(wěn)定性計算結(jié)果

由計算結(jié)果可知，在實際黏聚力c值降至50 kPa時，若繼續(xù)按原設(shè)計尺寸放坡采礦，正常工況(干燥無雨且施工放坡狀態(tài)良好)邊坡較為穩(wěn)定；暴雨影響下由簡化Bishop法計算出的安全系數(shù)較小，簡化Janbu法計算的安全系數(shù)僅為1.097，已接近臨界值；特殊工況(暴雨、地震共同作用)下邊坡穩(wěn)定性計算結(jié)果:簡化Bishop法的安全系數(shù)為1.004，而Janbu法僅為0.883，已不滿足穩(wěn)定性要求，極易出現(xiàn)失穩(wěn)垮塌情況。

通過上述計算和分析可以認為，巖石黏聚力c值顯著降低是巖性劣化的主要表現(xiàn)，也是造成露天礦山開采邊坡穩(wěn)定性急劇下降的關(guān)鍵內(nèi)因?；诖?，通過改變c值(其他參數(shù)不變)，進一步計算分析實際邊坡在巖性劣化過程中的穩(wěn)定性變化情況，探明巖性劣化對邊坡穩(wěn)定性的響應(yīng)規(guī)律。計算結(jié)果如圖8所示。

圖8 邊坡穩(wěn)定性隨黏聚力減小的變化趨勢

從圖8可看出，邊坡安全系數(shù)隨黏聚力下降呈明顯減小趨勢。通過多次曲線擬合計算和對比得出，二者為多項式關(guān)系，呈部分拋物線趨勢，3個方程的決定系數(shù)R2達到了0.999 9，擬合結(jié)果極為可靠。當(dāng)c值降到50 kPa時，邊坡特殊工況時的安全系數(shù)將達到臨界值。由表2知，原設(shè)計邊坡的黏聚力取值不小于200 kPa，證明邊坡在巖性劣化前穩(wěn)定性較好，黏聚力隨巖性劣化逐漸減小，進而直接降低了邊坡的穩(wěn)定性，導(dǎo)致臺階邊坡失穩(wěn)滑移。

3.4 修改建議

考慮到該露天礦開采已接近中后期，邊坡平臺不宜再做改動的現(xiàn)實狀況，通過反復(fù)分析演算，建議將上部第四系臺階坡面角調(diào)整為40°～45°，最終邊坡角為40°，其他尺寸值不變。對修改尺寸后的臺階邊坡再次進行穩(wěn)定性驗算，結(jié)果如表4所示。

表4 修改尺寸后的安全系數(shù)驗算結(jié)果

由驗算結(jié)果可知，3種工況下2種方法的安全系數(shù)均滿足穩(wěn)定性要求，修改后的尺寸合理可行。

4 結(jié) 論

(1)巖性劣化的主要原因是基巖存在軟弱夾層和破損結(jié)構(gòu)面且裂隙發(fā)育，抗水蝕能力低，極端降雨入滲，以及爆破動荷載導(dǎo)致滑動體內(nèi)部應(yīng)力重新分布。

(2)露天礦邊坡穩(wěn)定性是是動態(tài)變化的，其安全系數(shù)與黏聚力c值呈部分拋物線關(guān)系。

(3)邊坡3種工況下的安全系數(shù)均較小，尤其是在暴雨和特殊工況時已臨近失穩(wěn)，調(diào)整最終坡面角至40°后可滿足邊坡穩(wěn)定性要求。

(4)建議在設(shè)計初期對邊坡穩(wěn)定性進行計算和分析時，須結(jié)合礦山實際情況，考慮“巖性劣化”對采場邊坡后期運行的動態(tài)影響，否則極可能造成邊坡實際穩(wěn)定性較原設(shè)計偏小，存在滑坡危險。

[1] 楊天鴻，張鋒春，于慶磊，等.露天礦高陡邊坡穩(wěn)定性研究現(xiàn)狀及發(fā)展趨勢[J].巖土力學(xué)，2011，32(5):1437-1444. Yang Tianhong，Zhang Fengchun，Yu Qinglei,et al.Research status and development trend of open pit slope stability[J].Rock and Soil Mechanics，2011，32(5):1437-1444.

[2] 王啟明.我國非煤露天礦山大中型邊坡安全現(xiàn)狀及對策[J].金屬礦山，2007(10):1-6. Wang Qiming.Coal open-pit mine safety situation and counter measures in medium and large slope and nonmetal mines in China[J].Metal Mine,2007(10):1-6.

[3] 欒婷婷，謝振華，張雪冬,等.露天礦山高陡邊坡穩(wěn)定性分析及滑坡預(yù)警技術(shù)[J].中國安全生產(chǎn)科學(xué)技術(shù)，2013， 9(4):11-16. Luan Tingting， Xie Zhenhua， Zhang Xuedong，et al.The high steep slope stability analysis and landslide early warning technology of open pit mines[J].China Safety Science and Technology，2013，9(4)：11-16.

[4] 宋志飛，辛?xí)云?某露天礦邊坡穩(wěn)定性分析與優(yōu)化設(shè)計[J].煤礦安全，2014，45(7):216-218. Song Zhifei，Xin Xiaoping.Stability analysis and optimization design of an open pit mine slope[J].Mine Safety，2014，45(7)：216-218.

[5] 王永強，曹 竹，譚欽文,等.露天礦滑坡體的探地雷達檢測技術(shù)[J].煤炭學(xué)報，2011，36(7):1093-1097. Wang Yongqiang，Cao Zhu，Tan Qinwen，et al.Detection technique of ground penetrating radar open-pit mine landslide body[J].Journal of China Coal Society，2011，36(7)：1093-1097.

[6] 王思文，張緒杰，宋玉金.露天礦邊坡穩(wěn)定性的安全因素分析及防治建議[J].露天采礦技術(shù)，2011(3):21-23. Wang Siwen，Zhang Xujie，Song Yujin.Stability safety factor analysis and prevention advice of open-pit mine slope[J].Opencast Mining Technology，2011(3):21-23.

[7] 李大心.探地雷達方法與應(yīng)用[M].北京:地質(zhì)出版社，1994：20-21. Li Daxin.Method and Application of Ground Penetrating Radar[M].Beijing:Geological Publishing House,1994：20-21.

[8] 盧坤林，朱大勇.坡面形態(tài)對邊坡穩(wěn)定性影響的理論與試驗研究[J].巖石力學(xué)與工程學(xué)報，2014，33(1):35-43. Lu Kunlin，Zhu Dayong.Theoretical and experimental research on the influence of slope shape on the stability of the slope[J].Chinese Journal of Rock Mechanics and Engineering， 2014，33(1):35-43.

[9] 朱大勇，鄧建輝，臺佳佳.簡化Bishop法嚴(yán)格性的論證[J].巖石力學(xué)與工程學(xué)報，2007，26(3):455-458. Zhu Dayong，Deng Jianhui，Tai Jiajia.Theoretical verification of rigorous nature of simplified Bishop method[J].Chinese Journal of Rock Mechanics and Engineering，2007，26(3)：455-458.

[10] 黃景華，陳朝暉，莫 玻,等.參數(shù)特性及分布形式對邊坡穩(wěn)定可靠性的影響分析[J].四川大學(xué)學(xué)報:工程科學(xué)版， 2014，46(3):23-30. Huang Jinghua，Chen Zhaohui，Mo Bo，et al.Analysis of the influence of characteristics and distribution of form parameters on reliability of the slope stability[J].Journal of Sichuan University：Engineering Science，2014，46(3)：23-30.

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

Response Law of Lithology Deterioration on Slope Stability of Open-pit Mines

Xin Baoquan Tan Qinwen Xie Yujia Wang Yongqiang Wan Lu

(SchoolofEnvironmentandResource，SouthwestUniversityofScienceandTechnology，Mianyang621010，China)

In view of the fact that lithology deterioration condition leads to the decline of slope stability，and with a slope of open-pit mine as the research object，the internal defects of landslide were detected by SIR-20 geological radar，and the change law of slope stability under the condition of different lithology was analyzed with Slide software.Results showed that：slope lithology deterioration is caused by the interaction of weak interlayer，the rainfall infiltration and the blasting dynamic load，which lead to the slope sliding and instability in the late;cvalue reduction is the main mechanical performance in the process of lithology deterioration; Safety factor andcvalue showed as a partial parabola relationship，and the reducedcvalue 50 kPa is the critical value of slope instability under special working conditions.According to a comprehensive analysis of mine operating characteristics and the lithology deterioration condition，a proposal to improve its permanent stability by decreasing the upper quaternary bench slope angle and final slope angle is put forward.It is proven that the safety factor could meet the requirements.The results provide an important basis for safety mining and determining slope control parameters in the late.

Open-pit mine，Slope stability，Lithology deterioration，Numerical simulation，Cohesion，Safety factor

2015-02-01

國家自然科學(xué)基金項目(編號：51404200)，四川省教育廳科研項目(編號：14zd1106)。

辛保泉(1990—)，男，碩士研究生。

TD804，TU457

A

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