露天轉(zhuǎn)地下開采誘發(fā)高邊坡滑移機(jī)制研究

孫世國 郭煒晨 劉文波 郭 珮 董彥飛

(北方工業(yè)大學(xué)土木工程學(xué)院，北京 100144)

·安全與環(huán)?！?/p>

露天轉(zhuǎn)地下開采誘發(fā)高邊坡滑移機(jī)制研究

孫世國 郭煒晨 劉文波 郭 珮 董彥飛

(北方工業(yè)大學(xué)土木工程學(xué)院，北京 100144)

露天轉(zhuǎn)地下開采過程中邊坡穩(wěn)定問題是礦山安全生產(chǎn)的關(guān)鍵技術(shù)和核心難題之一。以工程實(shí)例為基礎(chǔ)，應(yīng)用離散元數(shù)值模擬，選取合適的力學(xué)參數(shù)，運(yùn)用3DEC、Matlab等軟件系統(tǒng)，分析了地下開采對高大邊坡穩(wěn)定性的影響機(jī)制、滑移特點(diǎn)和后續(xù)變形發(fā)展趨勢。研究結(jié)果表明：受露天與地下開采兩者的復(fù)合采動作用的影響，滑移變形產(chǎn)生了疊加作用，地表的下沉曲線形成多個大小不同的下沉盆地，導(dǎo)致不同區(qū)域邊坡巖體變形破壞程度也不同。根據(jù)數(shù)值模擬結(jié)果可知：地下采區(qū)的不同開采階段誘發(fā)的移動角大小不同，開采厚度不同，上覆巖體的破壞程度也不同。所以，露天轉(zhuǎn)地下開采中，邊坡穩(wěn)定性主要受到采空區(qū)處理方式和開采厚度等多重因素的影響；在邊坡穩(wěn)定性評價時，需要結(jié)合具體開采情況進(jìn)行分析和評價。

高陡邊坡 露天轉(zhuǎn)地下 離散元 3DEC 邊坡穩(wěn)定

邊坡穩(wěn)定性是我國巖土工程研究中的一個熱門課題。我國有許多采用露天轉(zhuǎn)地下開采的礦區(qū)[1]，在露天轉(zhuǎn)地下開采過程中，邊坡會前后受到2次采動(露天開采+地下開采)影響，且第2次擾動影響會誘發(fā)邊坡產(chǎn)生滑移變形或破壞。邊坡受到2次開采影響的綜合疊加作用，稱為復(fù)合開采效應(yīng)[2]。而在復(fù)合開采效應(yīng)影響下，對邊坡進(jìn)行穩(wěn)定性分析時，應(yīng)充分考慮2次采動所產(chǎn)生的綜合作用[3]。由于在露天轉(zhuǎn)地下開采過程中，地下采空區(qū)會對上覆邊坡穩(wěn)定性造成影響[4-6]，因此需要控制采空區(qū)對其影響程度。本研究運(yùn)用更適合節(jié)理巖土應(yīng)力分析的離散單元法[7]，以紫金山金銅礦為工程依托，模擬采空區(qū)采用不同處理方法時上覆巖體的變形破壞，并總結(jié)露天轉(zhuǎn)地下開采時的邊坡滑移變形規(guī)律，對邊坡穩(wěn)定性評價進(jìn)行綜合評判，對實(shí)際工程具有指導(dǎo)意義。

1 工程概況

紫金山金銅礦是我國的一個特大型金銅礦[8]。設(shè)計紫金山金銅礦露天采場上口面積約1 450 m×1 500 m，下口(采坑底)面積約80 m×90 m，終了邊坡東幫最高標(biāo)高為+1 012 m，采場底標(biāo)高+100 m，最大高差達(dá)912 m。紫金山銅礦所處區(qū)域構(gòu)造活動十分強(qiáng)烈，以北西向和北東向?yàn)橹?，是本區(qū)的控巖控礦構(gòu)造。北西向斷裂總體走向300°～329°，傾向南西或北東，傾角40°～80°，是本區(qū)重要的控巖控礦構(gòu)造。北東向構(gòu)造總體呈北東走向的“S”形，向南西傾狀，軸面傾向北西，由震旦系組成其核部，兩翼主要為上泥盆統(tǒng)和石炭系。兩側(cè)相伴有官莊—珊湖、廟前—舊縣復(fù)向斜，與褶皺相平行的北東走向斷裂十分發(fā)育，紫金山礦田即位于該背斜核部的西南端。斷層主要為北東向斷裂，走向40°～50°，傾向北西或南東，屬基底斷裂[9-11]。

2 計算模型與參數(shù)

計算模型如圖1所示。本次數(shù)值模擬所采用的力學(xué)參數(shù)取自巖石力學(xué)實(shí)驗(yàn)結(jié)果，并充分考慮了以往紫金山金銅礦地質(zhì)勘探中所記錄的力學(xué)參數(shù)，具體力學(xué)參數(shù)如表1所示。

圖1 紫金山金銅礦三維離散元計算模型

表1 紫金山金銅礦巖體力學(xué)參數(shù)

3 空區(qū)處理對高邊坡穩(wěn)定性的影響

邊坡穩(wěn)定性分析采用3DEC軟件進(jìn)行數(shù)值模擬,如圖1所示。將地下采區(qū)開挖分3個階段:分別從+30 m開采至0 m階段、-50 m階段和-100 m階段；地下開采模擬計算按采空區(qū)充填和不充填2種開采方式進(jìn)行：第1種是采空區(qū)按崩落處理，第2種采空區(qū)采用尾礦充填。選取5-5剖面(見圖1)，其地下開采步驟如圖2所示。

圖2 地下開采步驟

根據(jù)2種開采方法下的數(shù)值模擬結(jié)果，提取三維模型表面節(jié)點(diǎn)和不同開采階段的邊坡水平位移和豎直位移數(shù)據(jù)導(dǎo)入Matlab，可得出紫金山金銅礦+30 m開采至-100 m階段三維下沉盆地位移色譜圖，如圖3所示；在每次開挖過后，沿各剖面提取邊坡巖體坡頂至坡腳的關(guān)鍵點(diǎn)的豎直位移值，可畫出邊坡巖體在地下開采過程中的下沉曲線對比圖，如圖4所示。

圖3 +30 m開采至-100 m階段三維下沉盆地位移色譜圖

圖4 沉降曲線變化規(guī)律

由圖4可以看出，崩落法采空區(qū)對邊坡巖體下沉的影響遠(yuǎn)大于露天開采的影響。當(dāng)?shù)叵麻_采由+30 m開挖至-100 m階段時，采空區(qū)上覆巖體最大下沉值為-7.68 m，其中坡頂最大下沉值為-1.95 m；地表的變形破壞比較劇烈，采場地表最大沉降出現(xiàn)于采空區(qū)上方，對露天采場邊坡的穩(wěn)定性有很大影響，需要采取相應(yīng)的治理措施。而采用充填法處理采空區(qū)時，從填充開采整體三維下沉位移云圖和水平位移云圖(圖3(b))可以看出，采用填充開采時最大下沉區(qū)域出現(xiàn)在東幫，由于采用填充法所以地表下沉區(qū)域和下沉量顯著變小，且最大下沉位置向坡頂移動。根據(jù)填充后的邊坡體頂部各監(jiān)測點(diǎn)位移曲線變化(圖4(b))可以看出，填充采空區(qū)開采時地表的最大下沉位移為-0.34 m，遠(yuǎn)小于崩落采空區(qū)開采產(chǎn)生的沉降，因此，充填開采方法對邊坡穩(wěn)定性影響相對較小。

4 位移場變化分析

取5-5剖面如圖1所示，采用不同開采方法時其位移場對比如圖5所示，可以看出采用崩落法處理采空區(qū)時，邊坡巖體的移動值是充填法的22.6倍，說明采區(qū)充填與否，上覆巖體的破壞程度差異非常大。開采厚度不同其影響也不同，厚度越大其影響也越大，開采一、二、三階段的移動影響角分別為51°、47°和43°，說明移動開采厚度大，其對地表的影響范圍也越大，從而對邊坡穩(wěn)定性的影響范圍也越大，從數(shù)值模擬結(jié)果來看，采空區(qū)影響范圍基本覆蓋整個邊坡巖體(圖5)。可以看出由+30 m開采至0 m階段時，部分邊坡巖體已經(jīng)開始失穩(wěn)變形，由+30 m開采至-100 m階段時，采空區(qū)上覆巖體發(fā)生大面積的失穩(wěn)變形破壞。而從采用充填法開采時的位移場可以看出，充填之后地下采空區(qū)對露天邊坡體的影響大幅度的降低；從模擬結(jié)果來看，從+30 m開采至-100 m階段時，邊坡巖體局部才產(chǎn)生破壞(圖5(b))，且破壞程度遠(yuǎn)小于崩落法開采，說明充填法是控制聯(lián)合開采中邊坡穩(wěn)定程度的最有效的手段。

圖5 +30 m至-100 m開采階段5-5剖面位移對比

5 結(jié) 論

(1)針對復(fù)雜工程地質(zhì)條件高陡邊坡，采用了節(jié)理裂隙的離散元分析模型，解決非連續(xù)介質(zhì)模擬問題，并在建立模型時充分考慮特殊工程地質(zhì)條件和接觸面參數(shù)，為建立符合現(xiàn)場實(shí)際的模型提供了理論基礎(chǔ)。

(2)在露天轉(zhuǎn)地下開采過程中，由于兩者的復(fù)合作用，地表的下沉曲線是2次開采所造成的滑移變形的疊加結(jié)果。所以，在下沉盆地中產(chǎn)生了多個大小不同的下沉盆地，不同區(qū)域邊坡巖體滑移變形與破壞程度也不同。

(3)通過數(shù)值模擬研究，高陡邊坡在地下有采空區(qū)存在的情況下，不同開采階段的移動角大小是不同的。說明開采厚度不同，其對上覆巖體的破壞程度是不一樣的，在邊坡穩(wěn)定性評價時需要結(jié)合具體情況進(jìn)行分析。另外，采空區(qū)的處理方法不同，其移動量差距非常大，對于重要區(qū)段的邊坡必須采取充填方法，以確保邊坡的穩(wěn)定性。

[1] 孫世國,蔡美峰，王思敬.露天轉(zhuǎn)地下開采邊坡巖體滑移機(jī)制的探討[J].巖石力學(xué)與工程學(xué)報,2000,19(1):126-129. Sun Shiguo，Cai Meifeng，Wang Sijing.Study of sliding mechanism for slope due to the excavation from open pit into underground mine[J].Chinese Journal of Rock Mechanics and Engineering，2000,19(1):126-129.

[2] 孫世國,王思敬.地下與露天復(fù)合采動影響下邊坡巖體穩(wěn)定性評價方法的研究[J].工程地質(zhì)學(xué)報,1998， 6(4):312-318. Sun Shiguo，Wang Sijing.A method for evaluation of slope rock-mass stability under the combined action pit and underground mining[J].Journal of Engineering Geology,1998，6(4):312-318.

[3] 孫世國,蔡美峰,王思敬.地下與露天復(fù)合采動效應(yīng)及邊坡變形機(jī)理[J].巖石力學(xué)與工程學(xué)報,1999,18(5):563-566. Sun Shiguo，Cai Meifeng，Wang Sijing.Effect of combined underground and open pit mining and deformation mechanism of slope[J].Chinese Journal of Rock Mechanics and Engineering,1999,18(5)：563-566.

[4] 李遠(yuǎn)耀,唐朝暉,陳仁全.廣西合山煤田淺埋煤層采空區(qū)塌陷機(jī)理數(shù)值分析[J].金屬礦山,2014(3):26-30. Li Yuanyao，Tang Zhaohui，Chen Renquan.Numerical analysis on the mechanism of mining collapse in shallow seam of Heshan Coal Field,Guangxi[J].Metal Mine,2014(3):26-30.

[5] 許年春.巖體節(jié)理多參數(shù)反演的反射波法研究[D].重慶:重慶大學(xué),2007. Xu Nianchun.Research on Multi-parameter Inversion for Rock-mass Joint Using Reflected Wave[D].Chongqing:Chongqing University,2007.

[6] 袁海平，王繼倫，趙 奎，等.采空區(qū)形態(tài)對地表塌陷分布影響研究[J].金屬礦山，2011(11)：25-28. Yuan Haiping，Wang Jilun，Zhao Kui，et a1.Study on the goaf form influence on distribution of surface subsidence[J].Metal Mine，2011(11)：25-28.

[7] Itasca Consulting Group Inc.3DEC Three-dimensional Distinct Element Code Version 3.0 User's Manual[R].Itasca:Itasca Consulting Group Inc.,1988.

[8] 李長洪,張 永,蔡美峰,等.基于有限單元技術(shù)的紫金山金銅礦區(qū)地應(yīng)力場反演分析[J].金屬礦山，2008(7)：1-18. Li Changhong，Zhang Yong，Cai Meifeng，et al.FEM-based inversion analysis of in-situ stress field of Zijinshan Gold-copper Mining Area[J].Metal Mine,2008(7):1-18.

[9] 崔玉柱.連續(xù)與非連續(xù)介質(zhì)的數(shù)值模擬與拱壩地基系統(tǒng)安全分析[D].北京:清華大學(xué),2001. Cui Yuzhu.Numerical Modeling of Continuum and Discontinuous Media and Its Application to Safety Analysis of Arch Dam-foundation Systems[D].Beijing:Tsinghua University,2001.

[10] 劉獻(xiàn)華.紫金山金礦采空區(qū)處理技術(shù)研究[J].有色礦山,2002,31(1):21-43. Liu Xianhua.Technical research of goaf disposal in Zijinshan Gold Mine[J].Nonferrous Mines,2002,31(1):21-43.

[11] 饒運(yùn)章,柴 煒,黃奔文.ANSYS數(shù)值計算在紫金山金銅礦礦柱穩(wěn)定性分析中的應(yīng)[J].黃金,2008,29(7):22-26. Rao Yunzhang，Chai Wei，Huang Benwen.ANSYS numerical calculation of ore pillar stability in Zijinshan Gold-copper Mine[J].Gold,2008,29(7):22-26.

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

Study on High Slope Sliding Deformation Mechanism induced by Transiting from Opencast into Underground Mining

Sun Shiguo Guo Weichen Liu Wenbo Guo Pei Dong Yanfei

(SchoolofCivilEngineering，NorthChinaUniversityofTechnology，Beijing100144，China)

The slope stability during transition from opencast into underground mining is one of the key techniques and the core issues in safe mining production.Based on engineering project，the numerical simulation analysis of discrete element method is applied to select appropriate mechanical parameters.With the use of 3DEC，Matlab and other software systems，the slope stability，slip deformation characteristics and follow trends of high slopes at underground mining are analyzed.As the result，Due to the combined impact of both underground mining and open-pit mining，the surface subsidence curves show as multiple sink basins of different sizes under superposition，so the rock slope deformation extent at different regions are different.According to the results of numerical simulation，the mining-induced displacement angle at underground mining area are different at the different exploitation stage，indicating that with different mining thickness，the damage level of the overlying rock is different.So during the transition from opencast into underground mining，the slope stability is mainly impacted by multiple factors such as goaf treatment and different mining thickness.Furthermore，the slope stability evaluation requires a combination of the specific circumstances of the mining.

High and steep slope，Transition from open pit to underground mining，Discrete element method，3DEC，Slope stability

2015-03-02

國家自然科學(xué)基金項(xiàng)目(編號：41172250)，“十二五”國家科技支撐計劃項(xiàng)目(編號：2012BAK09B06)，北京市創(chuàng)新團(tuán)隊(duì)提升計劃項(xiàng)目(編號：IDHT20140501)。

孫世國(1961—)，男，教授，博士研究生導(dǎo)師。

TD853

A

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