汪 朝 周星宇 張文革

(1.西安建筑科技大學(xué)材料與礦資學(xué)院；2.四方金礦有限責(zé)任公司)

?

·采礦工程·

無底柱分段崩落法低貧化放礦數(shù)值模擬*

汪 朝1周星宇1張文革2

(1.西安建筑科技大學(xué)材料與礦資學(xué)院；2.四方金礦有限責(zé)任公司)

分析了無底柱分段崩落法現(xiàn)行截止品位放礦的缺點，提出了低貧化放礦管理方法。與前者相比，該方法的特點是以一個回采區(qū)域為統(tǒng)計目標(biāo)，上部回采時有意識的預(yù)留礦石，最后在最低水平集中出礦，從而減少礦石貧化過程。該方法在確?？偦厥章什唤档偷臈l件下，可以顯著減少礦石的貧化率。以某金礦為例，擬定了8種低貧化放礦方案，通過離散元數(shù)值模擬，以回貧差作為優(yōu)劣指標(biāo)，對各方案進行了比較，最終確定了上面兩個分段采用20%的礦石貧化率放礦，底下分段采用0.5 g/t的截止品位出礦時，回采指標(biāo)最優(yōu)。

低貧化放礦 PFC 數(shù)值模擬 回貧差

無底柱分段崩落法在開采礦床價值不高的厚大礦體時被廣泛使用。該采礦方法的優(yōu)點是生產(chǎn)能力大，設(shè)備機械化程度高，回采作業(yè)安全性好。缺點是在覆巖下放礦，損失率與貧化率較大[1]。因此，如何降低礦山回采過程中的兩率，成為該方法重點研究的問題。降低損失貧化的方法可以歸納為兩類：一類是優(yōu)化采場結(jié)構(gòu)參數(shù)，根據(jù)放礦理論，為了獲得最佳的出礦指標(biāo)，盡量使采場結(jié)構(gòu)參數(shù)組合(即分段高度、進路間距、崩礦步距)與放出體的形態(tài)相符合，在橢球體放礦理論中，提出了高分段和大間距的結(jié)構(gòu)布置方式[2]；另一類為實施放礦管理，通過控制出礦量與出礦過程來控制礦巖接觸面的形態(tài)，達到降低貧化損失的目的[3]。在實際生產(chǎn)中，這兩種方法一般會同時采用，本文擬對放礦管理中的低貧化放礦方案進行研究。

1 崩落法采礦礦巖移動過程與截止品位出礦

1.1 礦巖移動過程

與其他采礦方法相比，崩落法最大的特點是放礦過程中礦巖接觸面在不斷變化。根據(jù)放礦橢球體理論，無邊界限制的條件下，單漏斗放出體形態(tài)為一近似橢球體，放出后在礦巖堆中產(chǎn)生移動的部分稱為松動橢球體，在松動范圍內(nèi)各水平礦巖接觸面呈凹漏斗狀，稱為放出漏斗。假設(shè)放出體高度為Hf，放出漏斗水平面高度為h。放礦初始，礦巖接觸面都在放出橢球體上方(h>Hf)，為移動漏斗；隨著礦石的放出，放出漏斗逐漸移動,與放出體頂部平齊(h=Hf)，為降落漏斗；礦石繼續(xù)放出，漏斗水平將小于放出體頂部(h

1.2 現(xiàn)行截止品位出礦及其缺點

放礦時，若當(dāng)次礦石品位下降到一定數(shù)值時便停止出礦，此時的品位稱為截止品位。為了更多的放出礦石，現(xiàn)行放礦截止品位都取得比較低，該放礦方式存在以下問題[6-7]：

(1)采用盈虧平衡法確定出礦界線，片面強調(diào)單個步距的礦石回收率，勢必造成廢石混入過多，貧化率高；由于無底柱分段崩落法本中段損失的礦石可以在下中段回收，因此該放礦方式不是最優(yōu)。

(2)出礦時每個分段礦巖都會有大量廢石混合，貧化次數(shù)隨著分段數(shù)的增加而增加。

(3)截止品位放礦是以每個放礦步距單元內(nèi)的礦石被充分回收為核心，因此要求在放礦過程中必須隨時掌握采場的品位變化，以指導(dǎo)采場出礦。但由于目前國內(nèi)尚未有一種高效、及時并能適合井下特殊條件的品位分析儀表能滿足該要求，造成礦山生產(chǎn)時無法進行有效的相關(guān)控制。

綜上所述，由于截止品位放礦本身的工藝特點要求，放礦過程中的礦巖混合程度高，總體礦石貧化率大，且在采場生產(chǎn)時受到檢測儀表的限制而較難控制。在實際生產(chǎn)中為了降低礦石的貧化，建議實施低貧化放礦。

2 低貧化放礦原理

上部分段出礦時，截止品位較高，放礦截止時當(dāng)次貧化率比較低，允許適當(dāng)?shù)V石殘留在采場內(nèi)，殘留的礦石既是減少礦巖石的混合，同時獲得較高品位的采出礦石，并可以作為覆蓋層使用。最后一個分段放礦時，再采用現(xiàn)行截止品位出礦，盡可能多的放出礦石。與現(xiàn)行截止品位出礦相比，低貧化放礦的特點是[7-9]：①放礦截止品位較高，當(dāng)次貧化率較小，前期放出礦石品位較高；②有意識的預(yù)留部分礦石在采場里面，減少礦巖混合程度，綜合貧化率較??；③上分段預(yù)留的礦石在下分段被回收，增加了下分段的純礦石量；④根據(jù)噸品位曲線，采用控制放出礦量實施低貧化放礦管理，簡化了管理流程；⑤上分層部分礦石暫時殘留在空區(qū)內(nèi)，下分層可以在尚未形成足夠頂板巖石冒落的“覆蓋層”下繼續(xù)進行開采，即在“礦石+巖石”組成的覆蓋層下繼續(xù)進行下分層的開采，給上部巖石頂板的冒落提供了充分的時間；⑥在最終采出指標(biāo)上，該工藝在保持回收率基本相當(dāng)?shù)那闆r下，可以顯著降低礦石貧化率，提高產(chǎn)品質(zhì)量，減少無效費用，提高礦山經(jīng)營效果。

在生產(chǎn)實際中，低貧化放礦可根據(jù)礦山的具體情況，各分段采用不同的品位降低程度，組合出最優(yōu)的低貧化放礦方案。當(dāng)?shù)V山資金允許時可以采用大幅度降低貧化而提高采出品位方案，當(dāng)?shù)V山資金比較緊張則可以采用適當(dāng)降低貧化而穩(wěn)步提高采出品位方案，以避免礦山在應(yīng)用初期需要投入大量的資金進行三級礦量儲備。

3 低貧化放礦PFC數(shù)值模擬

3.1 礦山概況及低貧化方案確定

某金礦礦體厚大，品位低，圍巖中也含礦石，故采用無底柱分段崩落法開采，年產(chǎn)量100萬t，礦石平均品位1.2 g/t，廢石平均品位0.3 g/t。礦山現(xiàn)采用0.5 g/t的截止品位出礦，礦石損失貧化較大，經(jīng)常出現(xiàn)礦石被廢石包裹出不來的現(xiàn)象。本課題以該礦1 190 m中段為研究對象，上分段盡可能少出或不出廢石，最下面分段按現(xiàn)行截止品位出礦，研究最優(yōu)的低貧化放礦方案。根據(jù)礦山實際，擬定了8種放礦方案(見表1)，通過模擬分析，選定最優(yōu)方案。

放礦數(shù)值模擬采用PFC2D軟件，通過編程建立放礦離散模型，利用其自帶的fish語言，對放礦過程進行伺服控制，記錄放礦過程中的損失與貧化信息。

本次模擬程序的工作流程為[10]：①建立墻體并設(shè)置墻體參數(shù)，包括采場邊界、放礦進路；②生成廢石顆粒與中段礦石顆粒，并設(shè)置顆粒參數(shù)；③分段放礦模擬，根據(jù)各放礦方案設(shè)置的放礦截止品位進行放礦，統(tǒng)計放出的廢石與礦石體積；④放礦結(jié)果的輸出與分析。

表1 擬定放礦方案

生成放礦模型時，根據(jù)設(shè)定的孔隙率，在模型邊界內(nèi)生成全部顆粒，然后通過PFC中的擠壓排斥法將顆粒分開，為了避免在局部區(qū)域造成大面積的顆粒重疊，需要通過初始循環(huán)計算，使顆粒動能為零，達到初始應(yīng)力環(huán)境和條件，盡量與礦山的實際情況相同[11]。本次PFC模擬中放礦動力主要來源于礦石和廢石顆粒的自重，放礦阻力位顆粒間的摩擦力，兩種放礦方式的截止條件可以根據(jù)放出的礦石顆粒和廢石顆粒的比例進行控制。

3.2 模型建立

選取某一放礦工作面的垂直剖面建立放礦模型，如圖1所示。該模型共5個分段，寬121 m，高83.5 m，廢石覆蓋層厚度33.4 m。采礦結(jié)構(gòu)參數(shù)(高×寬)為16.7 m×22 m，進路尺寸(寬×高)為3 m×3 m，邊孔角為55°。模型的邊界條件：通過限制水平移動實現(xiàn)側(cè)面的固定要求，限制垂直移動實現(xiàn)底面的固定要求。

圖1 放礦計算模型(單位：m)

由于模型尺寸較大，崩落礦巖的平均直徑為30 cm，如果計算時將整個區(qū)域顆粒一次性全部生成，大約有112 328個球體，放礦時計算量太大、速度慢。因此，本研究采用分段建模、分段放礦。即先生成廢石與第1分段礦石，然后執(zhí)行放礦程序；待1分段放礦完畢后，再生成第2分段顆粒，執(zhí)行放礦；最后生成第3分段顆粒，執(zhí)行放礦直至結(jié)束。以方案1計算為例，各分段放礦PFC計算模型見圖2，放礦模型信息見表2。

圖2 方案1各分段放礦PFC模型

3.3 模型參數(shù)選取

(1)顆粒粒徑。根據(jù)現(xiàn)場測量的礦巖粒徑分布，本次模擬中，廢石粒徑取0.3～0.5 m，礦石粒徑取0.2～0.3 m，均服從均勻分布。

(2)顆粒密度。礦石密度取 2.75×103kg/m3，廢石密度取2.75×103kg/m3。

(3)顆粒間黏結(jié)強度。根據(jù)現(xiàn)場實際情況，爆破后的礦石和廢石都是松散狀態(tài)的堆積體，礦石顆粒和廢石顆粒間的連接關(guān)系已經(jīng)被破壞，但是仍舊存在較小的摩擦力。為了將該因素加入到數(shù)值模擬中，取礦巖顆粒間、墻體的法向剛度為108N/m，切向剛度為108N/m，墻體與顆粒間的摩擦系數(shù)為0.5，礦石顆粒的摩擦系數(shù)為0.1；廢石顆粒的摩擦系數(shù)為0.3，計算參數(shù)見表3。

表2 放礦模型信息

注：H為分段高度

表3 顆粒參數(shù)

3.4 計算過程控制

在本次PFC2D模擬放礦過程中，主要是依靠礦石和巖石顆粒的自重進行放礦。放礦過程和截止條件利用PFC2D自帶內(nèi)嵌的fish語言進行編寫，再根據(jù)編寫的fish程序指導(dǎo)計算機進行放礦[12]。模擬放礦過程的程序如下：①刪除分段各進路頂板，模擬爆破開挖；②運行30 000步，模擬爆破后礦巖運動過程；③刪除進路中的礦巖，模擬出礦過程，同時統(tǒng)計刪除的礦石與廢石顆粒個數(shù)、體積，并做累計；④計算本次出礦的廢石體積與礦石體積之比，與放礦截止值對比，如果大于放礦截止品位，則此進路放礦結(jié)束,若小于放礦截止值，回到第2步;⑤分段所有進路放礦結(jié)束則進入下分段放礦模擬。

4 放礦結(jié)果及分析

4.1 放礦模擬結(jié)果

本課題對每種方案的放礦過程進行了模擬，并對每種方案的放礦結(jié)果進行了統(tǒng)計與計算(篇幅有限，這里不予陳列)。放礦完畢后礦巖界線如圖3所示(以方案1為例)。

圖3 方案1各分段放礦完畢狀態(tài)

4.2 結(jié)果分析

崩落法放礦過程中，每一次放礦都希望礦石的回收率越大越好，礦石的貧化率越小越好。但在實際過程中，一味追求高回收率往往會導(dǎo)致礦石貧化率增大，為了減小礦石的貧化率，則要以犧牲回收率為代價，二者不能兼顧。因此，為綜合反映兩個指標(biāo)對放礦效果的影響，一般采用回收率與貧化率之差作為綜合指標(biāo)，即回貧差[13-14]。若用D表示回貧差，N表示礦石回收率，P表示貧化率，則回貧差可表示為：D=N-P。

通過比較8種放礦方案的統(tǒng)計結(jié)果，分別繪出礦石回收率、貧化率和回貧差在3個放礦分段的技術(shù)指標(biāo)，見圖4～圖7。

從圖中可以得出：

(1)在不同放礦方案條件下，相同分段的礦石回收率有差別；相同方案條件下，各分段的礦石回收率差別較大，整體趨勢來看，第1分段的回收率較小，到第2分段礦石的回收率有所增大，最后一個分段放礦時，礦石的回收率很大，甚至超過了本分段的爆破礦量。這是由于各分段采用了不同的放礦截止品位，前兩個分段截止品位較高，最后一個分段采用入選品位作為截止放礦品位。

圖4 不同方案分段礦石回收率

圖5 不同方案分段礦石貧化率

圖6 不同方案分段礦石回貧差

圖7 不同方案礦石綜合回貧差

(2)在不同放礦方案條件下，相同分段的礦石貧化率有差別；相同方案條件下，各分段的礦石貧化差別較大，整體趨勢來看，第1分段放礦的礦石貧化率較小，第2分段放礦的礦石貧化率次之，第3分段放礦的礦石貧化率最大。這是由于1、2分段采用了較大的放礦截止品位，控制了放礦時礦石的混入率，而第3分段放礦時采用入選品位，放礦截止品位較低，增大了礦石的混入。

(3)在不同放礦方案條件下，相同分段的礦石回貧差有一定的波動；相同方案條件下，各分段的回貧差差別較大。整體趨勢來看，第1分段的回貧差明顯比2、3分段的小，2、3分段的回貧差值相似，均在0.9左右。

(4)從方案的綜合回貧差曲線可以看出，方案1和方案8的回貧差明顯比其他方案要大，方案7的回貧差最小，其次是方案3、方案5、方案6、方案2、方案4。

5 結(jié) 論

(1)低貧化放礦法的礦石回收率要低于現(xiàn)行截止品位放礦法，但是隨著放礦分段的下移，兩者間的差距在逐漸減小，這是由于低貧化放礦法在首個分段需要部分礦石與巖石混合，形成礦石貧化層，在較高的品位時即停止放礦，而殘留的礦石將在下一個分段放出。

(2)低貧化放礦的廢石混入率要低于現(xiàn)行截止品位放礦法，而且隨著分段的下移，低貧化放礦的廢石混入率還將持續(xù)降低，而現(xiàn)行截止品位放礦法隨著分段的下移廢石混入率卻在不斷增加。

(3)低貧化放礦的綜合回貧差比現(xiàn)行截止品位放礦法要好，不同低貧化放礦方案條件下回貧差值也不一樣。本次研究中，當(dāng)上分段的截止品位越高時，回貧差越大，例如方案1、方案8分別為貧化20%，無貧化時，綜合回貧差分別為0.798、0.791。

綜上所述，低貧化放礦管理方式在確?；厥章实幕A(chǔ)上，對降低礦石的貧化率有著積極的作用，可以推廣應(yīng)用。

[1] 劉興國,張志貴.無底柱分段崩落法不貧化放礦[J].東北大學(xué)學(xué)報:自然科學(xué)版,1998(5):448-450.

[2] 吳愛祥,武力聰.無底柱分段崩落法結(jié)構(gòu)參數(shù)研究[J].中南大學(xué)學(xué)報:自然科學(xué)版,2012,43(5):1846-1850.

[3] 胡杏保.低貧化放礦工藝現(xiàn)狀及應(yīng)用前景[J].金屬礦山,2002(1):16-18.

[4] 孫 浩,金愛兵,高永濤.多放礦口條件下崩落礦巖流動特性[J].工程科學(xué)學(xué)報,2015(10):1251-1258.

[5] 徐 帥,安 龍.急斜薄礦脈崩落礦巖散體流動規(guī)律研究[J].采礦與安全工程學(xué)報,2013(4):512-517.

[6] 胡建華，郭福鐘，羅先偉.緩傾斜中厚礦體崩落開采礦石流動規(guī)律仿真與放礦參數(shù)優(yōu)化[J].中南大學(xué)學(xué)報:自然科學(xué)版,2015(5):1772-1777.

[7] Alford C G.Computer simulation models for the gravity flow of ore in sublevel caving[D].Melbourne:University of Melbourne,1978:383.

[8] 胡杏保.硬巖無底柱無巖覆蓋層開采技術(shù)[J].金屬礦山,2013(12):1-4.

[9] 孫文勇，陳星明.復(fù)雜礦體條件礦石損失貧化原因及對策[J].遼寧工程技術(shù)大學(xué)學(xué)報:自然科學(xué)版,2013(10):1339-1342.

[10] Zhou C B，Yao Y K，Guo L W，et al.Numerical simulation of independent advance of ore breaking in the non-pillar sublevel caving method[J].Journal of China University of Mining and Technology，2007，17(2)：295-300.

[11] 安 龍，徐 帥.基于多方法聯(lián)合的崩落法崩礦步距優(yōu)化[J].巖石力學(xué)與工程學(xué)報,2013(4):754-759.

[12] David J.Computer simulation of the movement of ore and waste in an underground mining pillar[J].The Canadian Mining and Metallurgical Bulletin，1968(61)：854-859.

[13] 王培濤,楊天鴻.無底柱分段崩落法放礦規(guī)律的PFC2D模擬仿真[J].金屬礦山,2010(8):123-127.

[14] 王連慶,高 謙.自然崩落采礦法的顆粒流數(shù)值模擬[J].北京科技大學(xué)學(xué)報,2007(6):557-560.

Numerical Simulation of Low Dilution Ore-drawing of Non-bottom Pillar Sublevel Caving Method

Wang Zhao1Zhou Xingyu1Zhang Wenge2

(1.School of Materials and Mineral Resources,Xi'an University of Architecture and Technology; 2.Sifang Gold Mine Co., Ltd.)

The weakness of current ore drawing by cut-off grade in non-bottom pillar sublevel caving method is analyzed,the low dilution ore-drawing method is proposed.Compared with the former method,the characteristics of low dilution ore-drawing method is the ore in the upper sublevel is reserved consciously,the large number of ore is drawn in the bottom sublevel,so,ore dilution is reduced.The method of low dilution ore-drawing method can significantly reduce the ore dilution rate during ore-drawing,at the same time the recovery rate can be lowered obviously.Take a gold mine as an example,eight eight schemes of low dilution ore-drawing are proposed.Using the value of recovery rate minus the dilution rate for the indicator, the eight schemes are compared by conducting the discrete element numerical simulation.The results show that when the top two sublevels drawing by cut-off grade with 20% dilution,and the bottom sublevel drawing by cut-off grade of 0.5 g/t, the value of recovery rate minus the dilution rate is the highest indicator.

Low dilution ore-drawing,PFC,Numerical modeling,Difference of the recovery and dilution

*國家自然科學(xué)青年基金(編號：51404182)；陜西省大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練計劃項目(編號：1204)。

2016-07-10)

汪 朝(1987—)，男，碩士，710055 陜西省西安市碑林區(qū)雁塔路13號。