陳禮石，雷 濤

(1.化工部長沙設(shè)計研究院， 湖南 長沙 410116；2.中南大學 資源與安全工程學院， 湖南 長沙 410083)

地下臺階式礦柱安全性的數(shù)值模擬研究

陳禮石1，雷 濤2

(1.化工部長沙設(shè)計研究院， 湖南 長沙 410116；2.中南大學 資源與安全工程學院， 湖南 長沙 410083)

根據(jù)杉樹埡磷7盤區(qū)的工程實際，提出了可適應礦體賦存條件的臺階式礦柱方案，基于細觀損傷力學理論，利用RFPA軟件對臺階式礦柱的安全性進行了數(shù)值模擬研究。結(jié)果表明，臺階式礦柱可以保證礦房的安全，礦房頂板圍巖的最大拉應力遠小于其抗拉強度；隨著開挖的進行，礦柱中的壓應力增大，損傷也隨之增加，但這些破壞僅發(fā)生在局部，對礦柱整體的安全性影響不大；臺階式礦柱具有常規(guī)礦柱相同的功能，且能提高礦石的回收率，具有一定的實用價值。

臺階式礦柱；細觀損傷；數(shù)值模擬；安全性

0 引 言

在地下采礦中，留下一定規(guī)格和數(shù)量的礦柱以保證采礦空間的安全是一種常用的方法[1]。因此，礦柱穩(wěn)定性的研究歷來就是采礦工作者們關(guān)注的熱點，并對比開展了卓有成效的研究。同時，合理設(shè)計礦柱，在保證安全的情況下盡量提高礦石的回收率也是設(shè)計者們面臨的重要課題，并在實際工程中取得了一定的效果。

但是，在這些研究中，礦柱的組成巖性和形式通常都是比較均勻的，而對橫截面形狀或巖性變化較大的礦柱的研究較為缺乏。為此，本文以實際工程為例，引入細觀有限元分析軟件RFPA，對巖性和橫截面形狀均存在差異的地下臺階式礦柱的穩(wěn)定性進行了數(shù)值模擬研究，并通過與傳統(tǒng)礦柱模擬結(jié)果的對比，指出這種礦柱的適用條件。

1 臺階礦柱概況

1.1 工程背景

杉樹埡磷礦位于湖北省宜昌市夷陵區(qū)，是一座設(shè)計年產(chǎn)量150萬t的大型地下磷礦山。目前主要的開采礦體為緩傾斜礦體，傾角一般為3°～10°，平均6°。根據(jù)礦體賦存條件及開采技術(shù)條件，主要采用無軌鏟運機出礦的盤區(qū)房柱法進行采礦。

但是，礦山7盤區(qū)的礦體厚度變化較大，并且礦體分布沒有明顯的規(guī)律性，頂板較破碎，穩(wěn)定性較差，且在厚度較大的地方，礦石品位很不均勻，表現(xiàn)出明顯的上層為富礦(中磷層)下層為貧礦(下貧礦)的特征。為了盡可能地回收資源，要求中磷層和下貧礦分開回采，實現(xiàn)貧礦與富礦分采分運，并提出了臺階式礦柱的方案。

1.2 臺階礦柱的形式

7盤區(qū)礦體中，中磷層的平均厚度約為6 m，下貧礦的平均厚度約為5 m，典型的臺階式礦柱如圖1所示。

圖1 臺階式礦柱的形式

由圖1可知，臺階式礦柱主要是為了多采富礦，提高資源回收率，故礦柱呈現(xiàn)為下部大上部小的形式。在圖1中，上部礦柱的截面尺寸為5×5 m2，下部礦柱的截面尺寸為6×6 m2。

可見，與截面為6×6 m2的均勻礦柱相比，圖1所示的臺階式礦柱可多回收66 m3的高品位礦石。

2 數(shù)值模型構(gòu)建

對保證生產(chǎn)活動的安全進行而言，礦柱自身的穩(wěn)定性具有非常重要的作用，在此對臺階式礦柱的穩(wěn)定性及相應采場的安全性進行數(shù)值模擬研究。根據(jù)所采用的采礦方法和臺階式礦柱的布置特點，數(shù)值模擬采用幾何模型如圖2所示。

圖2 數(shù)值計算幾何模型

圖2中，考慮到礦體傾角較小，為計算簡便，假定礦體的傾角為0°。同時為了更好地研究礦柱在生產(chǎn)過程的安全性，對圖2中的礦房按3步驟開挖，即先開挖中間礦房，再開挖左側(cè)礦房，最后開挖右側(cè)礦房。

根據(jù)礦體的賦存情況，取礦體埋深為400 m，即通過在數(shù)值模型的上部添加4×2.7 MPa的均布荷載反應豎直方向地應力的情況。同時在模型底部邊界上添垂直方向的位移約束，在兩個側(cè)面邊界均添加垂直該面的位移約束。

模擬采用RFPA軟件進行，該軟件基于細觀損傷力學，可直觀地模擬巖體材料中裂隙的產(chǎn)生和演化規(guī)律[2]。

巖體的破壞采用Mohr-Coulomb強度準則來判定，初始巖體力學參數(shù)如表1所示。

表1 初始巖體力學參數(shù)

3 結(jié)果分析

3.1 采場頂板最大拉應力分析

圖3 采場頂板最大拉應力變化曲線

在地下采礦中，采場頂板會發(fā)生由受壓向受拉轉(zhuǎn)變，其破壞形式一般也以受拉破壞為主。圖3即為在開挖過程中，采場頂板拉應力的變化曲線。

由圖3可知，采場頂板中拉應力隨著開采范圍的增加而不斷增大，但最大值也僅在3.5 MPa左右，遠小于上盤圍巖的抗拉強度，說明在臺階式礦柱的支撐作用下，采場頂板圍巖不會發(fā)生較大規(guī)模的塌落。

3.2 礦柱中的最大壓應力

在礦房的形成過程中，礦柱的整體破壞主要由壓應力來控制，礦柱中的最大壓應力變化情況如圖4所示。

圖4 礦柱中的最大拉應力變化曲線

由圖4可知，礦柱的最大壓應力隨開采的進行而逐漸增大，說明在礦體開采過程中，壓力不斷向礦柱轉(zhuǎn)移和集中，礦柱對采場的穩(wěn)定性具有明顯的支撐作用。同時值得注意的是，雖然仍保持增長的態(tài)勢，但第3時步時礦柱內(nèi)最大壓應力的增長率有減小的趨勢。

在完成3個開挖時步后，礦柱內(nèi)的壓應力達到最大值，最大壓應力達到23.7 MPa，但仍遠小于礦體的抗壓強度，說明礦柱發(fā)生整體破壞的可能性很小，是滿足采礦生產(chǎn)安全需求的。

3.3 礦柱的細觀損傷特征

如前所述，RFPA軟件是基于細觀損傷力學的原理開發(fā)出來的，因此在本模擬中可以通過裂隙的產(chǎn)生和發(fā)展直觀地展示礦柱的損傷破壞情況。礦柱隨開挖時步的細觀損傷變化特征如圖5所示。

圖5展示了第1到第3開挖時步中每個礦房形成過程中礦柱中損傷變化情況。由圖5可知，隨著開挖的進行，礦柱及其礦房附近圍巖中的裂隙越來越多，但這些裂隙還是主要出現(xiàn)在最后形成的礦柱中，這說明礦柱承受了最主要的地應力轉(zhuǎn)移，導致在礦柱中出現(xiàn)了比較明顯的應力集中。

值得注意的是，兩個礦柱的破壞情況并不完全一致，左側(cè)礦柱中裂隙數(shù)量和范圍明顯比右側(cè)礦柱要大，這說明礦柱的損傷特征與礦房的開挖順序有一定的關(guān)聯(lián)，先形成的礦柱將比后形成的礦柱承擔更大的壓力。但是，從圖5(c)中可知，盡管兩個礦柱中均有損傷破壞的情況產(chǎn)生，但這些破壞僅發(fā)生在局部范圍，且在臺階式礦柱的上部破壞情況明顯比下部要小，這說明礦柱是安全的，臺階式礦柱可以同均勻礦柱一樣，足以保證采礦生產(chǎn)的安全進行。

4 結(jié) 論

根據(jù)杉樹埡磷礦7盤區(qū)礦體的實際情況，提出了留臺階式礦柱的盤區(qū)房柱法進行采礦，并利用RFPA軟件對礦柱的安全性進行了數(shù)值模擬研究，結(jié)論如下：

(1) 模擬結(jié)果顯示，臺階式礦柱可以保證采場的安全，礦房頂板的最大拉應力及礦柱內(nèi)部的最大壓應力均遠小于巖體本身的強度；

(2) 隨著開挖時步的推進，礦柱中的最大壓應力和損傷情況均增加，表明礦柱發(fā)揮了支撐作用，承受了轉(zhuǎn)移的地應力，且礦柱的破壞情況與礦房的開挖順勢存在一定的關(guān)聯(lián)；

圖5 礦柱的細觀損傷變化特征

(3) 臺階式礦柱具有同常規(guī)礦柱同樣的功能，對一些礦體品位出現(xiàn)明顯分層現(xiàn)象的礦體具有較強適用性，在設(shè)計合理的前提下，可以在滿足采礦生產(chǎn)安全的同時提高礦石的回收率。

[1]古德生，李夕兵.現(xiàn)代金屬礦床開采科學技術(shù)[M].北京：冶金工業(yè)出版社，2006.

[2]梁正召，唐春安，張永彬，等.巖石直接拉伸破壞過程及其分形特征的三維數(shù)值模擬研究[J].巖石力學與工程學報，2008，27(7)：1402-1410.

2014-07-11)

陳禮石(1979-)，男，湖南臨武人，工程師，主要從事礦山設(shè)計工作，Email：43612618@qq.com。