張耀斌　楊　春

(1.西北有色地質勘查局七一七總隊；2.湖南省深部金屬礦開發(fā)與災害控制重點實驗室)

?

大尺寸采場充填體穩(wěn)定性數(shù)值分析

張耀斌1楊春2

(1.西北有色地質勘查局七一七總隊；2.湖南省深部金屬礦開發(fā)與災害控制重點實驗室)

摘要運用數(shù)值模擬軟件對某大型地下金屬礦山回采設計方案進行模擬，從最大主應力、最小主應力、塑性區(qū)、垂直位移4個方面,分析了處于最大暴露面積下的充填體穩(wěn)定性。結果表明：采用50 m高度回采時，充填體可能出現(xiàn)片幫；充填體與頂板相接位置形成了5 m左右的貫通塑性變形層，充填礦柱喪失了支撐能力，該回采方案存在較大安全隱患，需進行調整。

關鍵詞充填體礦柱頂板穩(wěn)定性

為擴大產能，實現(xiàn)采場高效開采，國內某地下金屬礦山委托相關科研單位進行了采場結構參數(shù)優(yōu)化、充填工藝等課題研究，提出了礦塊間隔布置的高效空場嗣后充填采礦法。由于回采過程中充填體最大暴露高度達50 m，有必要對其安全性進行評價，以便實現(xiàn)礦體的安全高效開采[1-3]。

目前數(shù)值模擬分析軟件在巖土工程領域應用廣泛，眾多學者利用數(shù)值模擬軟件進行采場穩(wěn)定性分析。白銀[4]采用理論和數(shù)值模擬相結合的方法，進行采礦方法安全性的數(shù)值模擬分析和采場結構參數(shù)的確定，最終選擇技術可行、經濟合理的采礦方法；王新民[5]對9種典型采場結構模型的模擬結果進行了比較分析,確定了采場最優(yōu)結構參數(shù),并對當前采礦方案提出修正和調整建議，以便高效生產。本文通過數(shù)值模擬軟件對礦體回采及充填過程進行模擬，從最大主應力、最小主應力、塑性區(qū)、垂直位移4個方面對采場穩(wěn)定性進行分析，研究礦塊間隔布置的高效空場嗣后充填采礦法在礦山的適用性。

1數(shù)值模擬的相關理論

在采礦、水工、地下工程中，常用的數(shù)值模擬方法包括有限差分法、有限單元法、邊界元法、半解析法、離散元法和無界元法等。有限差分程序是一個專門模擬地下巖體開挖工程應力及變形的三維分析軟件，對于模擬地下礦山開采、支護、隧道開挖，尤其是多步驟開挖有著很好的適用性。通過數(shù)值計算結果，可以得出不同開采步驟下的力學、位移變化特性，對采礦工程的穩(wěn)定性分析以及采場優(yōu)化研究有重要作用。

有限元分析軟件在模型構建及網格劃分方面較為方便，因此本次數(shù)值模擬在有限元分析軟件中建立計算模型并劃分網格，編組后導入有限差分軟件賦予材料參數(shù)進行計算，節(jié)省了有限差分軟件中的建模成本。

2數(shù)值模型構建

2.1基本假定

為便于計算，對采場結構特別是對回采支撐的充填體框架在模型中需作進一步處理，對礦體介質性質及計算模型作必要的假設：①將礦體視為理想的彈塑性材料，選用摩爾-庫侖模型，進行非線性分析；②忽略斷層、節(jié)理裂隙等不連續(xù)面對采場穩(wěn)定性的影響；③忽略地震波、爆炸沖擊波、地下水等因素對采場穩(wěn)定性的影響；④忽略地表起伏，根據圣維南原理，開挖影響范圍之外的力學作用轉化為相應荷載或約束。

2.2數(shù)值模型構建

該礦山一期開采礦體規(guī)模最大，形態(tài)及產狀較穩(wěn)定，平均水平厚度207 m(包括夾石)，傾角40°～50°，屬極厚大傾斜礦體。主要采礦方法為空場嗣后充填采礦法，在初步設計報告中推薦使用礦塊留間柱布置回采方案：盤區(qū)沿礦體走向布置，長度120 m，寬為礦體厚，高100 m(階段高度)，盤區(qū)間柱15 m。在盤區(qū)范圍內沿礦體走向和垂直方向劃分回采礦塊，礦塊長軸方向與礦體走向一致。一步開采采用膠結充填，二步以尾砂充填為主。一步礦塊連續(xù)布置，分步回采，礦塊長60 m，寬20 m，高100 m。二步礦塊亦連續(xù)布置，礦塊中間用間柱分割，礦塊長56 m，寬20 m，高100 m，間柱長20 m，寬8 m。一步采和二步采礦塊間隔排列如圖1所示。

圖1　礦塊布置

為了節(jié)省計算成本，選取采場5個礦塊作為研究對象。根據圣維南原理，開挖擾動影響范圍約為開挖尺寸的3～5倍，為簡化模型，不考慮間柱，以節(jié)省計算成本，垂直走向取5個礦塊，單個礦塊長60 m，寬20 m，高100 m，完整模型尺寸700 m×420 m×700 m，坐標原點位于模型幾何中心，礦塊模型如圖2所示。

圖2　礦塊有限元模型

一步驟1#、3#、5#礦塊按10 m分層高度依次上向回采，每一分層回采完成后采用膠結充填。待1#、3#、5#礦塊采完最后一個分層并充填接頂后，2#、4#礦塊按50 m分層高度分兩次回采，回采后形成的采空區(qū)采用尾砂非膠結充填。2#、4#礦體回采過程中1#、3#、5#膠結充填體最大暴露高度為50 m，此時膠結充填體處于最危險時刻，礦塊分層開挖模型如圖3所示。

圖3　礦塊有限元模型

2.3礦巖物理力學參數(shù)

2.3.1礦石力學參數(shù)

礦石力學參數(shù)通過室內試驗獲得，如表1所示。

表1　礦石物理力學參數(shù)

2.3.2非膠結充填材料力學參數(shù)

非膠結充填體材料采用尾砂，根據相關試驗參數(shù)進行類比[6]，參數(shù)見表2。

表2　非膠結充填材料力學參數(shù)

2.3.3膠結充填材料力學參數(shù)

根據科研單位試驗，研發(fā)了高強度膠結充填材料，試驗數(shù)據如表3所示。視充填體為均一材料，將試驗參數(shù)取均值作為計算參數(shù)。

2.4本構模型

構造模擬圍巖的數(shù)值模型需定義巖體性質的數(shù)學模型，即本構模型。由于計算研究范圍涉及的巖體(灰?guī)r、砂巖)、礦體、廢石膠結充填體均屬彈塑性材料，適應于莫爾-庫侖破壞準則，其力學模型為：

表3　全尾砂試驗參數(shù)

(1)

式中，σ1、σ3分別為最大和最小主應力；c、φ分別為材料粘聚力和內摩擦角；fs為破壞判斷系數(shù)，當fs≥0時，材料處于塑性流動狀態(tài),當fs≤0時，材料處于彈性變形階段。

3采場穩(wěn)定性分析

3.1一步膠結充填完成

礦塊分層回采后頂板失去支撐，產生較大垂直位移，拉應力也集中在頂板中央，壓應力則集中在空區(qū)側幫。礦體為彈塑性材料，當變形超過一定范圍即形成塑性區(qū)，每一分層開采過程中均有塑性區(qū)產生。在1#、3#、5#礦塊逐層膠結充填接頂后，充填體對采空區(qū)頂板側幫具有一定的支撐作用，頂板與側幫應力極值減小。因采空區(qū)塑性變形不可逆，故充填接頂后頂?shù)装寮皞葞臀恢玫呢Q直位移、塑性區(qū)范圍較未充填時無明顯變化。

如圖4所示，1#、3#、5#礦塊10個分層膠結充填完畢，采場無暴露空間，20 m×60 m×100 m膠結充填體與2#、4#待采礦柱支撐頂板。最大拉應力0.46 MPa，最大壓應力24.5 MPa，相比采場頂板有暴露空間時的應力場，充填體承擔了較大應力；10個分層充填單元周圍垂直位移12～20 cm，采場塑性變形主要集中在充填體與礦體接觸面附近。從第一分層的開挖與充填到第十分層開挖充填完成，采場的應力場趨于較為安全狀態(tài)，各項指標均未超過礦體強度極限。

圖4　一步充填(剖面位置Y=0)

3.211分層回采

此步驟回采是整個研究過程的關鍵，2#、4#礦塊11分層開挖后充填體有了暴露空間，特別是中央3#充填體下部雙側暴露。根據模擬結果，11分層回采后3#充填體上部主要表現(xiàn)為拉伸應力，應力值為0～1 MPa，下部為壓縮狀態(tài)，壓應力值為0～2.5 MPa；12分層側幫與充填體相接觸，與完整巖石相比此接觸面上粘結力較小，因此頂板產生較大垂直位移，約25 cm；充填體塑性區(qū)主要產生在兩幫及頂板附近，從應力角度，充填體中有大部分區(qū)域應力值已經超過了其強度極限。見圖5。

圖5　11分層開挖(剖面位置Y=0)

3.312分層回采

12分層開挖后3#充填體再次形成50 m的暴露高度，后續(xù)作業(yè)均在暴露頂板下進行，有必要對頂板及充填體礦柱進行穩(wěn)定性分析。如圖6所示，頂板與3#充填體接觸面附近產生了較大垂直位移，約26 cm，并且頂板全部產生塑性變形，與3#充填體接觸面上主要表現(xiàn)為拉伸塑性變形,說明充填體承受了較大豎向壓力；充填體內最大拉應力為1.63 MPa，已經超過其極限強度，同時3#充填體側幫塑性區(qū)已經上下貫通。從側面截圖也可以看出，在12分層回采后塑性區(qū)變化較為明顯，與頂板接觸的地方形成了高度約5 m且內部貫通的塑性變形層。

圖6　12分層回采頂板

4結論

通過對模型進行多方位剖視，從最大主應力、最小主應力、垂直位移、塑性區(qū)4個方面對采場穩(wěn)定性進行綜合分析。結果表明：11分層開挖時，3#充填礦柱大部分區(qū)域應力已經超過其強度極限，12分層開挖后充填體外圍形成了貫通的塑性區(qū)，暴露面可能產生大范圍冒落，從而削弱其支撐能力。故此方案存在極大的安全隱患，需進行調整，如選用高強度的充填材料或降低階段高度等。

參考文獻

[1]郭愛國,張華興.我國充填采礦現(xiàn)狀及發(fā)展[J].礦山測量,2005(1):60-61.

[2]周旭.我國充填采礦技術應用及存在問題探討[J].世界有色金屬,2014(3):28-30.

[3]張海波,宋衛(wèi)東,許英霞.充填采礦技術應用發(fā)展及存在問題研究[J].黃金,2013,31(1):23-25.

[4]白銀，王星.緩傾斜中厚礦體回采數(shù)值模擬分析[J].現(xiàn)代礦業(yè)，2015(8)：33-35.

[5]王新民,王長軍,張欽禮,等.基于ANSYS程序下的采場穩(wěn)定性分析[J].金屬礦山,2008(8)：17-25.

[6]曾照凱,張義平,王永明.高階段采場充填體強度及穩(wěn)定性研究[J].金屬礦山,2010(1):31-34.

(收稿日期2015-10-20)

Numerical Simulation Analysis of Filling Body Stability of Large Scale Stope

Zhang Yaobin1Yang Chun2

(1. 717 Team,Northwest Nonferrous Geological Exploration Bureau; 2. Hunan Key Laboratory of Mineral Resources Exploitation and Hazard Control for Deep Metal Mines)

AbstractThe stoping scheme of a large underground metal mine is simulated by numerical simulation software, the filling body stability under maximum exposure area is analyzed from four aspects of the maximum principal stress, the minimum principal stress, the plastic zone and the vertical displacement. The simulation results show that the filling body might appear wall falling phenomenon with 50 m mining height;The 5 m plastic deformation layer is occurred in the following location of filling body and roof,the supporting capacity of filling pillars is lost.Therefore,big security hidden dangers of the stoping scheme is existed,it is necessary to improved it.

KeywordsFilling body, Mine pillar, Roof, Stability

張耀斌(1973—)，男，工程師，721004 陜西省寶雞市金臺區(qū)。