深井硬巖礦山采場結(jié)構(gòu)參數(shù)優(yōu)化的數(shù)值模擬研究

許振華,喬 蘭

(1.北京科技大學(xué)，北京 100083, 2.中國瑞林工程技術(shù)有限公司，江西 南昌 330031；)

眾所周知，采礦方法是地下開采的核心[1]。對于深井硬巖礦床的開采而言，采礦方法和回采順序不僅要克服深井帶來的諸多不利影響，而且要簡化生產(chǎn)組織，適應(yīng)大型無軌設(shè)備強(qiáng)化高效開采的要求。但大設(shè)備就要求與之相適應(yīng)的大斷面，這無疑給采場的穩(wěn)定性帶來較大挑戰(zhàn)，因此采場單元尺寸大小及其比例關(guān)系不僅是保證安全回采的必要條件，同時對分步驟回采的綜合技術(shù)經(jīng)濟(jì)效果也具有重要影響。此外，影響采礦和地下工程穩(wěn)定性的因素很多，特別是各種影響因素的相互耦合，使得采場即使在相同的原巖應(yīng)力場及巖體強(qiáng)度參數(shù)條件下，不同的采場結(jié)構(gòu)參數(shù)，其圍巖及礦柱的應(yīng)力分布狀態(tài)、位移狀態(tài)也表現(xiàn)出不同的特性，這直接關(guān)系到采場自身的穩(wěn)定性。因此，如何在保證生產(chǎn)安全的前提下，合理選取采場結(jié)構(gòu)參數(shù)成為了是現(xiàn)場工程技術(shù)人員和設(shè)計工程師普遍關(guān)注的焦點。

1 有限元數(shù)值模擬

1.1 巖體力學(xué)參數(shù)的確定

數(shù)值計算結(jié)果的可靠程度在一定程度上取決于巖體宏觀力學(xué)參數(shù)選取的合理性，某銅礦巖體力學(xué)基礎(chǔ)研究初步結(jié)果，并依據(jù)巖塊試樣的物理力學(xué)參數(shù)和巖體分級指標(biāo)RMR 值，按照Hoek-Brown 提出的節(jié)理巖體經(jīng)驗破壞準(zhǔn)則對巖體宏觀力學(xué)參數(shù)進(jìn)行工程處理。經(jīng)過處理的巖體力學(xué)參數(shù)匯總于表1。

表1 經(jīng)工程處理的巖體力學(xué)參數(shù)匯總

1.2 力學(xué)模型的建立

本次計算采用日本軟腦株式會社開發(fā)的有限元結(jié)構(gòu)應(yīng)力分析系統(tǒng)3D-σ。為了簡化研究結(jié)構(gòu)巖體,本設(shè)計選取模型尺寸為300m×300m。地應(yīng)力場按自重應(yīng)力場計算，最不利開采部位距地表的平均深度約800m。模型底部采用垂直方向約束，其余側(cè)面均采用水平方向約束。模型未開挖前，共有單元76440個，節(jié)點總數(shù)為324153個，滿足一般數(shù)值分析的要求[2-4]。

2 采場跨度優(yōu)化計算結(jié)果及分析

某銅礦擬采用階段大直徑深孔嗣后充填法開采，即將礦體劃分為礦房和礦柱，第1步驟回采礦房，用膠結(jié)充填形成人工礦柱；第2步驟回采礦柱，用全尾砂充填。在回采順序上擯棄了傳統(tǒng)的“隔三采一”，采用更為靈活先進(jìn)的“隔一采一”二步回采順序。那么，當(dāng)回采一步采場時，由于礦房充填不及時，隨著礦房暴露面積的增大，采場內(nèi)圍巖頂、底板及礦柱所承受的應(yīng)力也逐漸增大。當(dāng)一步回采結(jié)束時，若其頂板仍然處于穩(wěn)定狀態(tài)，則可認(rèn)為在回采過程中采場處于穩(wěn)定狀態(tài)[5]。因此，研究不同跨度下采場內(nèi)圍巖頂、底板及礦柱所承受的應(yīng)力狀態(tài)及其變化規(guī)律具有重要意義。

作為方案優(yōu)化研究來說，最主要的是選擇有關(guān)的特征量進(jìn)行比較。由于巖體材料對拉伸作用的抵抗能力極差，這里主要選取最大拉應(yīng)力作為方案比較的依據(jù)。其他部分指標(biāo)，如最大壓應(yīng)力、最大剪應(yīng)力和位移值作為方案比較的輔助參量。本文對采場高度為80m的硬巖深井采場穩(wěn)定性進(jìn)行分析研究，分別進(jìn)行了采場跨度為15m、20m、25m、30m工況下的應(yīng)力傳播規(guī)律、變形機(jī)制及安全率進(jìn)行分析，分析計算結(jié)果見表2。

表2 采場高度80m時五種不同采場跨度尺寸的有關(guān)計算結(jié)果

2.1 應(yīng)力分析

根據(jù)既定的總體回采順序，某銅礦井下采場直接頂板為石英閃長玢巖，其極限抗拉強(qiáng)度為2.60MPa左右，以此值作為判別標(biāo)準(zhǔn)。通過線性插值，當(dāng)主應(yīng)力σ1最大拉應(yīng)力為2.6MPa時，對應(yīng)的采場跨度為25m左右；而當(dāng)正應(yīng)力(σy)的最大拉應(yīng)力為2.6MPa時，對應(yīng)的采場跨度介于25～30m之間?？梢?，二者對應(yīng)的采場跨度值均為25m多，不到30m，但其拉應(yīng)力值基本都達(dá)到巖體極限抗拉強(qiáng)度的99%以上，甚至超出巖體極限抗拉強(qiáng)度。因此綜合考慮到現(xiàn)場巖體穩(wěn)定狀況的差異，建議將采場高度為80m時的采場跨度值定為20m。

同時，由圖1～8可知，主應(yīng)力σ1最大拉應(yīng)力隨著采場跨度的增大，呈現(xiàn)線性增長關(guān)系，且主要分布在采場兩幫的中部以及礦柱的腰部。但當(dāng)采場跨度增大至25m時，主應(yīng)力σ1最大拉應(yīng)力急劇下降，這是由于采場兩幫已超出巖體本身的塑性破壞區(qū)域，發(fā)生了應(yīng)力釋放，從而使得拉應(yīng)力減少。σy方向上最大拉應(yīng)力主要分布在采場頂板和底板上，在數(shù)值上隨著采場跨度增大而增加，但在15～25m之間有一個平緩的過渡，證明在此跨度范圍內(nèi)，應(yīng)力相對平衡，調(diào)整幅度較小，采場相對穩(wěn)定，但一旦跨度達(dá)到30m，σy方向上最大拉應(yīng)力就急劇增大，這是巖體進(jìn)入了屈服區(qū)的典型特征；最大剪應(yīng)力均分布在采場的四個邊角，同時在數(shù)值上也表現(xiàn)出與σy方向上最大拉應(yīng)力相似的特征；σ3最大壓應(yīng)力主要集中在采場四個邊角，在數(shù)值上達(dá)到43.79MPa，為巖體極限抗壓強(qiáng)度的73.5%，因此，在生產(chǎn)中應(yīng)該盡可能避免這種直角造成的應(yīng)力集中，采用圓弧過渡。

圖1 σ1最大拉應(yīng)力-跨度曲線

圖2 采場跨度為20m時σ1最大拉應(yīng)力云圖

圖3 σy最大拉應(yīng)力-跨度曲線

圖4 采場跨度為20m時σy最大拉應(yīng)力云圖

圖5 最大剪應(yīng)力-跨度曲線

圖6 采場跨度為20m時最大剪應(yīng)力云圖

圖7 σ3最大壓應(yīng)力-跨度曲線

圖8 采場跨度為20m時σ1最大拉應(yīng)力云圖

2.2 位移分析

由圖9～10可知，采場開挖后，位移均指向采空區(qū)，即頂板下沉，底板上移；隨著開采區(qū)域的擴(kuò)大，位移值呈逐漸增大趨勢。

圖9 最大位移量-跨度曲線

圖10 采場跨度為20m時位移云圖

2.3 安全率分析

采用安全率指標(biāo)S 定量表示巖體材料接近破壞的程度，S>1 表示位于屈服面內(nèi)部，S≤1表示位于屈服面上部或屈服面外部。因此，S 值愈小說明破壞程度愈厲害。由圖11可知，當(dāng)采場跨度達(dá)到20m時，采場頂板和兩幫，以及礦柱的安全率均大于1，證明此時的采場是穩(wěn)定的。但從數(shù)值上看，采場頂板及兩幫局部區(qū)域安全率均為1.1，說明此時礦柱及采場頂板均處于屈服狀態(tài)，但已接近臨近狀態(tài)。

圖11 采場跨度為20m時的安全率曲線

3 結(jié)論

綜合以上巖石力學(xué)數(shù)值分析結(jié)果，某深井硬巖銅礦基于充填采礦“隔一采一”的回采順序情況下，當(dāng)采場高度為80m時，采場跨度定為20m，無論從應(yīng)力分布、位移量，還是安全率上均為合理，既保證了采場的穩(wěn)定，又最大程度上發(fā)揮了采場的生產(chǎn)能力，有利于大型高效采礦設(shè)備的應(yīng)用。

[1] 童光煦.高等硬巖采礦學(xué)[M].北京：冶金工業(yè)出版社,1995:143-148.

[2] B.H.G.布雷迪,E.T.布朗.地下采礦巖石力學(xué)[M].馮樹仁等,譯.北京:煤炭工業(yè)出版,1990.

[3] 王清來，許振華.復(fù)雜采空區(qū)條件下殘礦回收與采區(qū)穩(wěn)定性的有限元數(shù)值模擬研究[J].金屬礦山，2010,409(7)：37-40.

[4] 沈慧明，許振華, 朱利平,等.殘礦回采順序優(yōu)化與復(fù)雜采空區(qū)穩(wěn)定性的有限元模擬研究[J].中國礦業(yè),2011,20(1):78-81,103.

[5] 陳何,程國江.垂直分條充填采礦法采場參數(shù)的數(shù)值分析[J].礦冶,2001,3(10):1-5,10.