王炳文，黎 林，張書金，武茂超

(1.中國(guó)礦業(yè)大學(xué)(北京)能源與礦業(yè)學(xué)院，北京 100083；2.山東黃金礦業(yè)(玲瓏)有限公司，山東 招遠(yuǎn) 265419)

上向進(jìn)路充填法是金屬礦山內(nèi)常用的一種采礦方法，因礦石回收率高、貧化率低等優(yōu)點(diǎn)被廣泛應(yīng)用[1]。同時(shí)，許多礦山由于采場(chǎng)的跨度過于保守，也存在著生產(chǎn)效率低、出礦成本較高等缺點(diǎn)[2]。因此，上向進(jìn)路充填法擁有一個(gè)合理的采場(chǎng)結(jié)構(gòu)參數(shù)是十分必要的。

在采場(chǎng)結(jié)構(gòu)參數(shù)優(yōu)化方面，許多學(xué)者進(jìn)行了相關(guān)的研究。歐任澤等[3]基于3D-σ有限元法對(duì)礦山可能采用的四種采場(chǎng)結(jié)構(gòu)參數(shù)方案在不同埋藏深度時(shí)的穩(wěn)定性進(jìn)行了數(shù)值模擬研究，獲得了不同埋藏深度時(shí)最佳采場(chǎng)結(jié)構(gòu)參數(shù)；杜宇翔等[4]運(yùn)用成本分析模型對(duì)程潮鐵礦崩落法采場(chǎng)就應(yīng)優(yōu)化，得到了開采經(jīng)濟(jì)效益最優(yōu)的結(jié)構(gòu)參數(shù)；程文文等[5]應(yīng)用多因素層次分析法對(duì)金山金礦淺孔房柱采礦法采場(chǎng)結(jié)構(gòu)參進(jìn)行了數(shù)優(yōu)化，解決了采場(chǎng)生產(chǎn)能力較低、貧損指標(biāo)不理想等問題；陳順滿等[6]基于響應(yīng)面法對(duì)破碎圍巖條件下的采場(chǎng)結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化，實(shí)現(xiàn)了白牛廠礦破碎圍巖條件下安全、高效回采；葛文杰等[7]、毛貴林[8]利用MIDAS-GTS得到了采場(chǎng)結(jié)構(gòu)參數(shù)的最優(yōu)方案；易善華等[9]基于ANSYS和模糊評(píng)價(jià)對(duì)采場(chǎng)結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化，得到了礦房和礦柱寬度的最優(yōu)組合方案；李寧等[10]運(yùn)用正交試驗(yàn)原理對(duì)某厚大鐵礦采場(chǎng)進(jìn)行參數(shù)優(yōu)化，實(shí)現(xiàn)了高效經(jīng)濟(jì)開采。

玲瓏金礦東風(fēng)礦區(qū)采用上向進(jìn)路膠結(jié)充填采礦法，生效率低，工人勞動(dòng)強(qiáng)度大。針對(duì)玲瓏金礦的現(xiàn)場(chǎng)情況，采用FLAC3D軟件分別模擬采場(chǎng)跨度為4 m、5 m、6 m、7 m、8 m、9 m時(shí)的頂板應(yīng)力位移分布，并結(jié)合實(shí)際情況，得出最合理的采場(chǎng)跨度參數(shù)。

1 工程背景

東風(fēng)礦區(qū)內(nèi)的礦脈主要表現(xiàn)為沿?cái)嗔哑扑閹С涮罱淮纬傻奈g變帶。 主礦脈為171號(hào)脈，分布于160線至56線間，礦脈長(zhǎng)2 600 m，寬300 m左右，走向NE60°～70°，傾向SE，傾角35°～45°。礦脈中賦存有1711號(hào)主礦體，呈大脈狀，嚴(yán)格受斷裂構(gòu)造控制。其上盤次級(jí)裂隙發(fā)育，巖石也較破碎，由碎裂狀花崗巖、絹英巖化碎裂巖等組成；下盤巖石相對(duì)完整，由絹英巖及絹英巖質(zhì)碎裂巖、花崗質(zhì)碎裂巖組成。局部具有明顯的裂面，在主裂面發(fā)育有斷層泥。區(qū)內(nèi)地層簡(jiǎn)單，主要為第四系陸相沉積，沿山間谷地分布。區(qū)內(nèi)構(gòu)造為脆性斷裂構(gòu)造,受多期次的區(qū)域性構(gòu)造作用，礦區(qū)內(nèi)斷裂構(gòu)造發(fā)育。北東向斷裂是主要的斷裂構(gòu)造，控制著含金礦脈的分布，北西向斷裂構(gòu)造次之。按地下水類型及富水特征，劃分為第四系松散巖類孔隙水、基巖風(fēng)化裂隙水、構(gòu)造裂隙水和隔水層(體)。

根據(jù)現(xiàn)場(chǎng)調(diào)查所知，現(xiàn)有開采方法采場(chǎng)跨度為3.3 m，由于采場(chǎng)結(jié)構(gòu)參數(shù)過于保守，導(dǎo)致了采場(chǎng)生產(chǎn)能力低、生產(chǎn)成本高等問題，且采場(chǎng)斷面過小，無(wú)法引入大型機(jī)械化鑿巖鉆孔設(shè)備，降低了工作效率。因此，需要對(duì)現(xiàn)場(chǎng)結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化，提高礦山的生產(chǎn)效率。

2 數(shù)值模擬分析

利用FLAC3D軟件進(jìn)行數(shù)值模擬，分析不同跨度采場(chǎng)的頂板應(yīng)力和位移分布，對(duì)玲瓏金礦東風(fēng)礦區(qū)的采場(chǎng)結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化

2.1 模型建立

計(jì)算模型的幾何范圍，垂直礦體走向400 m，礦體走向500 m，豎直方向-570 m到-770 m，高度為200 m。模擬區(qū)域礦體傾角平均40°，礦體厚度取40 m。模型尺寸為500 m×400 m×200 m，如圖1所示。其中，坐標(biāo)軸Z方向?yàn)樨Q直方向，坐標(biāo)軸Y方向沿礦體走向，坐標(biāo)軸X方向垂直礦體走向。

結(jié)合實(shí)際情況，邊界條件確定如下：充填采礦法的采動(dòng)影響范圍有限，在采場(chǎng)水平四周和較遠(yuǎn)處的巖體移動(dòng)值將很小，計(jì)算域邊界采用位移約束，模型的4個(gè)側(cè)面X方向和Z方向位移約束為零，模型底面Y方向位移約束為零，上部平面承受21.6 MPa的面載荷。采場(chǎng)埋置較深，礦區(qū)主應(yīng)力屬于水平構(gòu)造應(yīng)力主導(dǎo)的高應(yīng)力場(chǎng)，因此，計(jì)算域初始條件考慮原巖構(gòu)造應(yīng)力場(chǎng)和自重應(yīng)力場(chǎng)的共同作用。

2.2 模擬方案

模型計(jì)算采用莫爾-庫(kù)倫強(qiáng)度準(zhǔn)則，對(duì)已開挖巖體賦予空殼模型。根據(jù)玲瓏金礦地質(zhì)調(diào)查結(jié)果及室內(nèi)巖石力學(xué)參數(shù)試驗(yàn)，礦巖體力學(xué)參數(shù)見表1。根據(jù)現(xiàn)場(chǎng)調(diào)查可知，礦區(qū)內(nèi)采場(chǎng)高度統(tǒng)一為3.2 m，因此本次模擬只對(duì)采場(chǎng)跨度進(jìn)行優(yōu)化，共分為6個(gè)方案，各方案參數(shù)見表2。

圖1 礦體FLAC3D模型Fig.1 FLAC3D model of ore deposit

表1 模型計(jì)算力學(xué)參數(shù)Table 1 The mechanical parameters of model calculation

表2 各模擬方案采場(chǎng)結(jié)構(gòu)參數(shù)Table 2 Stope structure parameters of each simulation scheme

2.3 各方案模擬結(jié)果分析

根據(jù)所模擬的結(jié)果采場(chǎng)對(duì)頂板應(yīng)力分布和位移分布進(jìn)行分析，鑒于篇幅，僅給出方案5(跨度8 m)和方案6(跨度9 m)的頂板應(yīng)力位移分布圖以及各方案匯總的曲線圖。

2.3.1 應(yīng)力分析

從模擬結(jié)果可以看出，隨著采場(chǎng)跨度的增大，頂板拉應(yīng)力也越來(lái)越大，當(dāng)超過礦體的抗拉強(qiáng)度時(shí)，頂板就會(huì)發(fā)生破壞。圖2和圖3分別為方案5(跨度8 m)和方案6(跨度9 m)的應(yīng)力分布圖，當(dāng)采場(chǎng)跨度為8 m時(shí)，開挖后，采場(chǎng)頂板迅速出現(xiàn)應(yīng)力集中，最大拉應(yīng)力為2.40 MPa，十分接近礦體的最大抗拉強(qiáng)度；當(dāng)采場(chǎng)跨度為9 m時(shí)，頂板的最大拉應(yīng)力為2.70 MPa，已經(jīng)超過礦體的最大抗拉強(qiáng)度，頂板發(fā)生破壞。

從圖4中的頂板最大拉應(yīng)力曲線可以看出，頂板最大拉應(yīng)力的總體變化趨勢(shì)是隨著采場(chǎng)跨度的增大而增加，在跨度超過6 m后增長(zhǎng)較為明顯，當(dāng)跨度超過8 m后，頂板發(fā)生破壞。

圖2 方案5(跨度8 m)應(yīng)力分布Fig.2 Scheme 5 (span 8 m) stress distribution

圖3 方案6(跨度9 m)應(yīng)力分布Fig.3 Scheme 6 (span 9 m) stress distribution

圖4 頂板最大拉應(yīng)力和最大下沉量變化曲線Fig.4 Curve of maximum tensile stress and maximum subsidence of roof

2.3.2 位移分析

頂板的最大位移量也是評(píng)價(jià)采場(chǎng)穩(wěn)定性的一個(gè)重要指標(biāo)，因?yàn)?，通過觀察頂板的位移量也可以了解頂板此時(shí)的穩(wěn)定狀態(tài)。

圖5和圖6分別為方案5(跨度8 m)和方案6(跨度9 m)的位移分布圖。結(jié)合圖4中的頂板最大下沉量曲線可以看出，在進(jìn)行開挖之后，頂板會(huì)出現(xiàn)不同程度的下移，總體趨勢(shì)是隨著采場(chǎng)跨度的擴(kuò)大下沉量不斷增加；從圖4中還可以看出，跨度為4 m到6 m時(shí)下沉量增加較快，6 m到7 m時(shí)趨于平緩，7 m之后又快速增加。

圖5 方案5(跨度8 m)位移分布Fig.5 Scheme 5 (span 8 m) displacement distribution

圖6 方案6(跨度9 m)位移分布圖Fig.6 Scheme 6 (span 9 m) displacement distribution

2.3.3 分析結(jié)論

由模擬結(jié)果可以看出，目前技術(shù)條件下采場(chǎng)跨度8 m為極限跨度，但由于數(shù)值模擬的簡(jiǎn)化性以及下場(chǎng)施工條件的復(fù)雜性，建議生產(chǎn)時(shí)采場(chǎng)跨度不得超過6 m。

3 工程應(yīng)用

3.1 現(xiàn)場(chǎng)試驗(yàn)情況

東風(fēng)礦區(qū)上向水平分層大斷面進(jìn)路膠結(jié)充填采礦工業(yè)試驗(yàn)采場(chǎng)位于-260水平三分段，136勘探線附近，礦體平均厚度18 m，走向30°左右，傾向NE，傾角31°左右，平均品味1.2 g/t。根據(jù)試驗(yàn)地點(diǎn)采場(chǎng)巖體穩(wěn)定性，將原來(lái)設(shè)計(jì)的采場(chǎng)跨度3.3 m更改為6 m，進(jìn)路沿著走向布置，爆破方式為聚能光面爆破，從現(xiàn)場(chǎng)施工情況來(lái)看，采場(chǎng)壁面光滑整潔，部分破碎區(qū)域需要支護(hù)，試驗(yàn)效果良好。

3.2 技術(shù)經(jīng)濟(jì)對(duì)比

采場(chǎng)跨度變更為6 m后，對(duì)試驗(yàn)點(diǎn)的技術(shù)經(jīng)濟(jì)指標(biāo)進(jìn)行了統(tǒng)計(jì)，并與之前跨度為3.3 m時(shí)進(jìn)行了對(duì)比，具體數(shù)據(jù)見表3。

表3 試驗(yàn)前后技術(shù)經(jīng)濟(jì)指標(biāo)對(duì)照Table 3 Comparison of technical and economic index before and after the test

從表3可以看出，采場(chǎng)跨度擴(kuò)大到6 m后，損失率和貧化率明顯降低，出礦量增加，生產(chǎn)成本減小，實(shí)現(xiàn)礦山的高效經(jīng)濟(jì)開采，具有很大的實(shí)用價(jià)值。

4 結(jié) 論

1) FLAC3D軟件模擬結(jié)果表明，在該礦山地質(zhì)條件下，8 m為采場(chǎng)的極限跨度，超過8 m頂板就會(huì)發(fā)生破壞，鑒于模擬的簡(jiǎn)化性，建議現(xiàn)場(chǎng)采場(chǎng)跨度不得超過6 m。

2) 現(xiàn)場(chǎng)工業(yè)試驗(yàn)表明，采場(chǎng)跨變更為6 m后，損失率和貧化率明顯降低，出礦效率提高，效果良好；斷面擴(kuò)大后引進(jìn)的大型鑿巖鉆孔設(shè)備，減輕了工人的勞動(dòng)強(qiáng)度，實(shí)現(xiàn)了礦山的安全、高效、經(jīng)濟(jì)開采。