溫軍鎖，程文文，繆國衛

(1.中國黃金集團公司，北京 100011；2.長春黃金研究院，吉林 長春 130012)

金廠峪金礦基于FLAC3D的殘礦回收與空區處理采礦方法優化

溫軍鎖1，程文文2，繆國衛2

(1.中國黃金集團公司，北京 100011；2.長春黃金研究院，吉林 長春 130012)

金廠峪金礦近年來井下空區數量不斷增大，空區條件愈加復雜，可能隨時出現地表塌陷等危害，亟需開展廢石尾砂膠結充填技術進行空區處理及殘礦回收研究工作。根據類似礦山經驗，提出4種殘礦回收方案，采用FLAC3D分別對4種方案進行數值模擬計算與分析，結果表明：第1種方案的豎向應力、位移沉降、塑性區分布以及最大主應力等指標優于其他方案，并能保證生產過程中的整體安全穩定性，推薦該方案在金廠峪金礦大規模推廣應用。

充填；FLAC3D；采礦方法

河北省遷西縣金廠峪礦業有限責任公司下屬的桑家峪金礦，位于河北省遷西縣城城北，公路里程30km。礦區位于西灣子村西0.5km，屬金廠峪鎮管轄。礦區東側有豐董公路通過，交通方便。金廠峪金礦桑家峪礦區為該礦新收購的礦山，之前進行過民采，已開采中段有110m、80m、50m中段。由于民采期間未進行合理規劃，即見礦就采，至此，在此三中段形成了約4萬m3的采空區，所形成的高大空區與地表接近，位于空區上方的山谷為一季節性河流，可能隨時出現河水灌入井下、地表塌陷等危害。因此，需進行空區的處理，并運用相應技術進行殘礦回收。

根據礦山的具體情況，決定采用廢石尾砂膠結充填技術進行空區處理及殘礦回收工作[1-2]，以達到處理空區確保不發生災害事故并安全回收殘礦的目的，提出幾種初步方案設計，并經數值模擬分析回采方案過程，確定最終回采設計方案。

1 數值模擬方法概述

巖體力學數值模擬相關技術的發展為礦塊回采工藝科學決策提供了條件，本研究以巖土工程專用分析軟件FLAC3D為輔助，在保證生產安全的前提下優化礦塊回采順序及工藝，模擬計算了在不同的回采工藝條件下，礦塊開采后礦房頂板、礦柱的應力、變形狀態及塑性區分布，根據計算結果確定安全合理回采工藝，提高礦山安全生產水平和資源的有效利用率[3-4]。

2 礦塊回采工藝優化數值模擬

2.1 模型的建立

本次模擬范圍從+30m到+130m共100m，設計地表標高為+200m，上部巖體對模型作用力通過施加豎向應力實現。根據設計選用采礦方法進行建模以及開挖分析等工作；最終選取建模參數為：長100m，寬50m，高100m，包含了+50m及+80m等中段巖體。礦體按照厚度15m，傾角70°，走向長度80m進行建模。整個模型分為兩個組塊，為簡化計算，將礦體外其余部分設為圍巖，其中對礦體網格進行了加密處理，以便得到更加準確的分析結果。模型共31920個單元體，34231個節點。

2.2 本構模型的選擇

本研究采用莫爾-庫侖(Mohr-Coulomb)準則。

2.3 模型力學參數

由于實驗室方法是對所取巖芯進行測試獲取巖體的力學屬性，其值較原位試驗獲取的巖體的力學屬性往往要大得多，因此在實踐中通常是將實驗室測試獲取的巖體的力學屬性進行必要的折減后作為模擬需要的物理力學參數[5]。

本次模擬主要巖性為斜長角閃巖，根據金廠峪巖石力學試驗總結報告選取模型巖石力學參數見表1。

3 模型開挖計算方案

3.1 開挖方案建立

采用上向水平分層充填法進行采礦，每隔體回采順序可從一側向另一側進行，每個礦塊作為單獨采場，每個分層為獨立充填單元，回采時采用隔一采一。

由于+80m上部有空區賦存，當殘礦回收至空+80m中段時存在著一定的安全隱患，為了避免在區下進行回采作業，因此考慮在+80m中段底板預先用尾砂膠結充填1m厚，既能提高下方殘礦回收作業的安全性，也在一定程度上對兩幫圍巖起到一定的支撐作用。對整體回收方案考慮四種方案進行分析和對比，見表2。

表1 巖石物理力學參數

表2 模型開挖方案

由于上向水平分層充填采礦工藝較為復雜，為了簡化開挖過程及便于分析對比方案，模擬過程對回采步驟進行必要的簡化、合并處理。

3.2 開挖過程分析模擬

殘礦回收模擬前要先形成+50m、+80m中段空區，在+50m中段充填至+65m后，將+80m底板尾砂膠結充填1m，見圖1。

模型對開挖過程按照采礦工藝過程分析相關豎向應力、位移沉降以及塑性區分布等參數，總結其變化規律，探討礦塊結構參數對圍巖的破壞機理及規律[6-9]。

圖1 +50m空區充填后有、無底板膠結示意圖

方案1開挖過程模擬，見圖2。方案1在+80m中段有1m后的尾砂膠結，回采順序為一側向另一側回采，分析豎向應力云圖可知，隨著回采的進行，殘礦巖體內豎向應力變化較小，這是因為殘留礦體上下空間無巖體賦存，其主要受兩幫圍巖的水平應力作用。而分層充填采礦方法的特點也使采場周邊圍巖得到有效支撐，主要變化集中在采區兩幫圍巖內，其豎向應力有所增大，但變化幅度較小。

3.3 開挖方案對比分析

上文對選取的開挖方案過程進行了分析模擬，通過分析豎向應力場、沉降位移變化，對礦巖的開采引起周圍巖層擾動規律有了具體的認識，下面通過上述參數對各個方案進行橫向比較分析，選取回采順序最優的方案。

由圖3可知，觀察各個方案回采結束后豎向應力對比云圖可知，各方案都在周邊圍巖形成應力集中，在+80m中段空區上部有應力貫通區域，整體應力分布大體一致。而無1m尾砂膠結的方案3、4周邊圍巖豎向應力較方案1、2要小，從而使圍巖內拉應力較方案1、2增大，從巖體整體穩定性來看是相對不利的，因為巖石抗拉強度相對抗壓強度要低很多。

這也說明預先施工的1m厚膠結尾砂對兩幫圍巖起到了支撐作用，從作業安全性及巖體穩定性兩方面考慮都是很有必要的。經過比較可知，四個方案中，方案1、2優于方案3、4，而方案1在前述分析中是優于方案2，可視為最優方案。

分析各個方案回采結束后豎向位移云圖見圖4，方案3、4周邊圍巖沉降幅度較方案1、2要大，由于方案3、4上部無膠結頂板，導致采場充填體內發生位移、移動的充填體范圍增大，綜合分析，可視方案1為最優方案。

分析各方案回采結束后最大主應力對比云圖見圖5，回采結束后在周邊圍巖及上部空區頂部均有拉應力出現，整體分布規律大體一致，且大部分區域拉應力在5MPa以下，低于巖體抗拉強度，整個采場可視為是安全穩定的。雖然在+50m采場端部出現較大拉應力集中，但是由于其范圍較小及處于端部位置，不會影響到整個采場區域的安全及穩定性。對比各方案主應力云圖可以看到，方案1其應力集中程度最小，整體安全穩定性最好，可視為最優方案。

圖3 各方案回采結束后豎向應力對比云圖

圖4 各方案回采結束后豎向位移對比云圖

4 結論

1)通過對開挖方案建立模型，建立4個開挖方案，方案之間進行縱向分析并做橫向比較，分析選取的豎向應力、位移沉降、塑性區分布以及最大主應力等參數，最終得到不同回采順序的采場穩定性對比如下：方案1>方案2>方案3>方案4。最終推薦的礦塊結構參數為方案1，也就是在+80m中段底板預先施工1m厚膠結尾砂并由一側回采。一側回采不僅在采場安全及穩定性上優于中間向兩側回采，也便于施工組織，合理利用采準工程，加快工程進度。

2)分析各方案開挖過程的力學參數可以看到，方案之間的力學參數指標差異較小，方案的安全及穩定性都是可以接受的；但是方案3、4由于上部無預先施工膠結頂板，導致回采后期要在空區下作業，存在安全隱患，因此不建議采用。

圖5 各方案回采結束后最大主應力對比云圖

3)針對在分析模擬過程中出現的采場充填體會出現塑性破壞以及沉降等情況，建議礦山在生產采礦過程中，做好施工組織及規劃，在采場回采結束之前做好充填作業所需相關準備工作，回采結束后盡快進行充填，縮短回采周期，避免采場長時間的暴露；而在回采礦房時，在充填體附近區域作業時考慮調整爆破參數或減小一次爆破礦量，以此來減少爆破震動對充填體造成的影響。

[1] 古德生，李夕兵.現代金屬礦床開采科學技術 [M]．北京：冶金工業出版社，2007．

[2] 歐陽化兵．河臺金礦礦柱回采方法的改進[J]．黃金，2002，23(2)：22-25.

[3] 李培現，譚志祥，王磊，等.FLAC在老采空區地基穩定性評價中的應用研究[J].煤礦安全，2009，40(10):11-14.

[4] 黨建印，鄭建明，孫豁然，等.FLAC在東坪金礦空區處理中的試驗研究[J].有色礦冶，2001，17(6)：5-8.

[5] 楊濤波.充填體下水平礦柱回采優化及地壓控制研究[D].贛州：江西理工大學，2012.

[6] 劉杰，萬文，趙延林，等.緩傾斜中厚礦體采空區穩定性數值模擬分析[J].礦業工程研究，2011，26(4):44-48.

[7] 饒運章，楊慶明，吳國興，等.基于FLAC(3D)程序的馬坑鐵礦采空區穩定性分析[J].四川有色金屬，2011 (1)：9-13.

[8] 孫國權，李娟，胡杏保.基于FLAC(3D)程序的采空區穩定性分析[J].金屬礦山，2007，368(2)：29-32.

[9] 趙奎.殘留礦柱回采與空區穩定性的有限元模擬研究[J].采礦技術，2002(2)：51-53.

The optimization of mining methods by FLAC3Din stooping safely and pillar extraction of Jinchangyu Gold Mine

WEN Jun-suo1，CHENG Wen-wen2，MIAO Guo-wei2

(1.China National Gold Group Corporation，Beijing 100011，China;2.Changchun Gold Research Institute，Changchun 130012，China)

With the continuous increase of the number of underground mine cavity in Jinchangyu Gold Mine in recent years，the cavity condition become more and more complicated.It may cause the surface subsidence at any time.So it is very urgent to carry out cemented waste rock backfills technology for cavity treatment and the stoping of remnant ores research.According to the similar mine experience，it proposed 4 remnant ores recovery schemes，and use FLAC3D to do the numerical simulation and analysis for the 4 schemes.The results showed that the first scheme of vertical stress，displacement and settlement and distribution of plastic zone and the maximum principal stress index is better than the other schemes.It tan guarantee the the overall safety and stability of production process.Therefore the first scheme is recommend for large-scale application in Jinchangyu Gold mine.

filling;FLAC3D;mining methods

2014-12-10

溫軍鎖(1960-)，男，任職于中國黃金集團公司，主要從事采礦技術及生產管理工作。E-mail：cuisong0502@163.com。

TD82

A

1004-4051(2015)03-0111-05