充填體下水平礦柱回采方案研究

廖偉成 劉 波 張耀平,3 徐佑民

(1.江西理工大學資源與環(huán)境工程學院，江西 贛州 341000；2.安徽龍橋礦業(yè)有限公司，安徽 廬江 231551；3.江西理工大學應用科學學院，江西 贛州 341000)

充填體下水平礦柱回采方案研究

廖偉成1劉 波2張耀平1,3徐佑民1

(1.江西理工大學資源與環(huán)境工程學院，江西 贛州 341000；2.安徽龍橋礦業(yè)有限公司，安徽 廬江 231551；3.江西理工大學應用科學學院，江西 贛州 341000)

安慶銅礦在生產(chǎn)過程中為確保高中段回采安全，在-385～-400 m之間留有約15 m厚的水平礦柱。隨著開采深度加大，預留礦柱量越來越多，為最大限度地回采礦石，現(xiàn)對水平礦柱進行回采。根據(jù)工程實際情況，對水平礦柱的回采方法提出上向中深孔嗣后充填，下向中深孔嗣后充填，上向進路分層充填3種方法，3種方法在回采效率和安全性上進行對比，著重考慮回采過程安全性，建議采用上向進路分層充填法進行回采。同時，針對上向充填法的回采順序提出中間向兩翼，左端向右端，右端向左端，兩翼向中間等4個方案，耦合surpac與FLAC3D構(gòu)建礦區(qū)數(shù)值模擬分析模型，對該4個方案在-385 m中段進行開挖模擬運算。分析結(jié)果得出：兩翼向中間回采方案垂直位移量最小、最大拉應力較小、塑性區(qū)域最小，為最優(yōu)回采方案。

水平礦柱 回采方案 數(shù)值模擬

為了維護采場的穩(wěn)定性及控制地壓一般需要留下各式各樣的礦柱。水平礦柱是礦柱中較典型的一種，一般為上中段的底板和下中段的頂板所構(gòu)成，起到承載采場及空區(qū)的作用，確保采場及空區(qū)的穩(wěn)定性[1-2]。隨著礦體開采的不斷持續(xù)，可開采礦體越來越少，礦柱量越來越多，若能安全有效地開采礦柱，能有效地提高礦山的經(jīng)濟效益，延長礦山的服務年限。

1 工程概況

安慶銅礦是一座大型矽卡巖銅鐵共生型地采礦山，礦體主要特點：埋藏較深、沿走向較短、儲量大[3]。1#礦體是礦區(qū)當前的主采礦體，埋藏范圍為-176 m至-640 m，屬于厚大急傾斜礦體。沿礦體走向方向，將礦體劃分為礦房、礦柱2部分，兩者寬度均是15 m，采場階段高度60 m～120 m，采礦方法采用VCR法，炮孔直徑達165 mm，先礦房后礦柱兩步驟、隔三采一的方式回采[4]。以安慶銅礦1#礦體預留的水平礦柱為研究背景，對其回采方法進行了論證，并利用數(shù)值模擬對回采方案進行研究。

2 回采方法選擇

基于水平礦柱處于充填體中現(xiàn)狀，礦柱被充填體包裹，其穩(wěn)定性較好。根據(jù)水平礦柱賦存的現(xiàn)狀，初步選擇以下3種回采方法：上向中深孔嗣后充填、下向深孔嗣后充填及上向分層充填[6-8]。

(1)上向中深孔嗣后充填。采場的寬度設定為15 m，長度即為水平礦柱賦存的厚度，垂直礦體走向方向布置。上向中深孔孔徑60 mm，底距1.8 m，排距1.4 m。采場充填采用尾砂膠結(jié)充填。其鑿巖巷道布置在采場中間位置。聯(lián)絡道布置在相鄰采場內(nèi)，出礦進路布置在出礦聯(lián)絡道內(nèi)。

(2)下向中深孔嗣后充填。采場的布置形式、參數(shù)及充填材質(zhì)和上向中深孔嗣后充填的方法一致。布置單側(cè)聯(lián)絡道在相鄰采場內(nèi)。鑿巖巷道布置在水平礦柱上方的充填體之中，在鑿巖巷道內(nèi)打下向孔，炮孔直徑100 mm，底距2.2 m，排距2.0 m。

(3)上向進路分層充填。采場的布置和前2種方法相同，將1個采場劃分為3個進路，單個進路寬度為5 m，回采過程中進路的順序先中間后兩側(cè)，分層高度為5 m，采用淺孔落礦方式，采場回采后采用尾砂膠結(jié)充填。

將3種回采方法的優(yōu)缺點對比后可得：上向及下向中深孔嗣后充填法均采用中孔落礦，出礦量大，回采效率高；但鑿巖巷道布置在充填體中，作業(yè)安全性較差，且采場空間較大，頂板的穩(wěn)定性難以確保，其充填體有可能產(chǎn)生垮落對礦柱回采造成不利影響。上向分層充填法采切工程量小，采場結(jié)構(gòu)參數(shù)小，頂板及周邊圍巖容易控制，支護比較容易，損失及貧化率低。缺點是回采效率低 。

著重考慮回采過程的安全性，結(jié)合三者優(yōu)缺點，建議采用上向分層充填法回采該礦柱。

3 數(shù)值模型建立及賦值

3.1 模型建立

采用surpac建立塊體模型，給塊體模型各單元賦以坐標等屬性數(shù)據(jù)，將塊體模型以數(shù)據(jù)文件的形式導出，利用access對導出的數(shù)據(jù)進行編輯轉(zhuǎn)化成FLAC3D識別的數(shù)據(jù)格式，再通過FLAC3D中的impgrid命令調(diào)入數(shù)據(jù)建立FLAC3D數(shù)值模型[9]。將模型向外圍發(fā)散原有的3～5倍，得整體的數(shù)值模型，如圖1。

圖1 幾何模型及沿礦體走向縱剖面圖

3.2 模型參數(shù)賦值

根據(jù)前期地質(zhì)調(diào)查可知該礦區(qū)最大主應力方向由北東逐漸到東西的趨勢。礦區(qū)地應力場回歸方程如下：

σ1=0.019 8H+11.033,

σ2=0.010 1H+5.356,

σ3=0.024 5H+0.298,

式中，σ1，σ2，σ3分別為最大、最小水平主應力及垂直方向主應力的值，MPa；H為礦體賦存深度，m。

采取左右兩端邊界約束水平方向的位移，底部約束z方向的位移設邊界限制。初始應力采取在計算區(qū)域內(nèi)直接賦值，垂直方向施加σ3，水平面的兩個邊界應力根據(jù)最大主應力的方向分別施加作用力。

數(shù)值模擬涉及材料參數(shù)主要分成閃長巖、大理巖、礦體、充填體5類。充填體分為膠結(jié)充填體及尾砂充填體2種。數(shù)值模擬的全部材料力學參數(shù)參見表1。

表1 各種材料類型的力學參數(shù)

4 回采方案數(shù)值模擬

4.1 回采方案的模擬

水平礦柱上、下盤圍巖巖性不同，賦存環(huán)境復雜，再加上礦體的回采活動影響，使得不同部位應力顯現(xiàn)不同，因此，本次數(shù)值模擬將整個水平礦柱作為研究對象。對水平礦柱回采順序提出4個回采方案，綜合數(shù)值模擬分析結(jié)果，從而選擇最優(yōu)方案。

方案一：由中間向兩翼回采，即從水平礦柱中間向兩翼的推進回采。

方案二：由左端向右端回采，即從水平礦柱左端向右端推進回采。

方案三：由右端向左端回采，即從水平礦柱右端向左端推進回采。

方案四：由兩翼向中間回采，即從水平礦柱兩翼分別向中間推進回采。

4.2 模擬結(jié)果分析

本次對水平礦柱回采提出4個方案，并分別進行模擬，受限于篇幅，每個方案模擬的結(jié)果只分析-385 m中段的位移、最小主應力及塑性區(qū)范圍大小。在此基礎上對水平礦柱回采之后的穩(wěn)定性作出評價，從而得出最優(yōu)回采方案。

(1)位移分析。按照4個方案的回采順序開挖，依次將水平礦柱全部開挖完之后，選取水平礦柱頂部-385 m水平垂直位移等色圖對水平礦柱回采完之后的位移進行分析，如圖2～圖5。由圖可知方案一最大垂直下沉位移值量大約為1.79 cm，最大上升位移為0.7 cm，大部分區(qū)域的位移為下沉其值為0.75～1 cm；方案二下沉移量最大值大約為1.71 cm，上升位移最大值為0.65 cm，大部分變形為下沉其位移量為0.75～1 cm；方案三下沉位移量的最大值1.76 cm，上升位移最大值約為0.69 cm，大部分區(qū)域的位移為下沉其值0.75～1 cm；方案四的下沉位移最大值約為1.62 cm，上升位移最大值為0.66 cm，位移區(qū)域以下沉為主其值主要在0.50～0.75 cm之間。從4個方案的-385 m水平的整體垂直位移分析，方案四即從兩端向中間回采的垂直位移值最小，優(yōu)于其他3個方案。

圖2 方案一-385 m水平z方向位移

圖3 方案二-385 m水平z方向位移

圖4 方案三-385 m水平z方向位移

圖5 方案四-385 m水平z方向位移

(2)應力分析。圖6～圖9為4個方案回采結(jié)束后，水平礦柱頂板-385 m水平最小主應力分布圖。從圖中可以得出水平礦柱回采結(jié)束后，方案一的最小主應力最小值大約為-30.8 MPa，在水平礦柱中部上方及左下角形成明顯的應力集中，水平礦柱區(qū)域最小主應力值西部區(qū)域大部分為-14 MPa左右，東部區(qū)域大部分為-20 MPa左右，最小主應力值下降值較大，范圍主要集中在水平礦柱西部；方案二最小主應力最小值大約為-29.1 MPa，水平礦柱區(qū)域最小主應力西部區(qū)域大部分為-16 MPa左右，東部區(qū)域最小主應力在-20～-24 MPa之間；方案三最小主應力最小值最小值大約為-29.2 MPa，水平礦柱最小主應力西部區(qū)域大部分為-16 MPa左右，東部區(qū)域最小主應力在-20～-24 MPa之間；方案四最小主應力最小值約為-29.3 MPa，水平礦柱最小主應力東部區(qū)域大部分為-18.1 MPa，東部區(qū)域最小主應力在-20～-24 MPa之間。通過分析-385 m水平最小主應力可得知，方案一最小主應力最大，其他3個方案相當；方案一最小主應力值下降較其他3個方案稍微小一點；方案一在水平礦柱中部上方及左小角的應力集中也更明顯，在東部存在的應力集中區(qū)域較其他3個方案大，其他3個方案在東部的應力集中區(qū)域相當。整體分析可得，其他3個方案效果相當，優(yōu)于方案一。

圖6 方案一-385 m水平最小主應力

圖7 方案二-385 m水平最小主應力

圖8 方案三-385 m水平最小主應力

圖9 方案四-385 m水平最小主應力

(3)塑性區(qū)分析。圖10～圖13是水平礦柱4個方案回采之后，沿走向剖面圍巖整體塑性區(qū)域分布圖。從4個方案模擬結(jié)果的整體塑性區(qū)域來看，遠離水平礦柱區(qū)域的上下盤圍巖，4個方案的塑性區(qū)基本一致，這說明這些區(qū)域基本上沒有受到水平礦柱開采的影響。對比分析4張圖可以看出，水平礦柱回采結(jié)束后，方案一和方案二的水平礦柱塑性區(qū)域最大，且主要集中在兩端；方案三的塑性區(qū)域較方案一、二的要小，分布也主要集中在兩端；水平礦柱塑性區(qū)域最小的為方案四。通過分析塑性區(qū)域的分布情況，方案四要優(yōu)于其他3個方案。

圖10 方案一整體塑性區(qū)域分布

圖12 方案三整體塑性區(qū)域分布

圖13 方案四整體塑性區(qū)域分布

從垂直位移、最大主應力、最小主應力、塑性破壞區(qū)域大小4個方面對4種回采方案的數(shù)值模擬結(jié)果進行了綜合分析，從而得出最優(yōu)的回采方案。將綜合對比的結(jié)果列入表2。

表2 回采順序綜合分析結(jié)果

5 結(jié) 論

(1)結(jié)合現(xiàn)場實際情況，對水平礦柱回采3種方法進行比對，最終采用上向分層充填回采水平礦柱。

(2)用 surpac軟件建模簡便、快捷的特點，基于surpac的FLAC3D自動生成復雜礦體計算模型的方法，克服了FLAC3D前處理能力不足的缺陷。使得模型更接近工程實際。

(3)對水平礦柱的回采的4個方案進行數(shù)值模擬，對模擬結(jié)果進行詳細的分析論證，最后得出最優(yōu)的回采順序為兩翼向中間回采，確保回采過程的安全。

[1] 劉 波.復雜條件下水平礦柱穩(wěn)定性分析及其開采方案研究 [D].贛州：江西理工大學，2013. Liu Bo.Level Pillar Stability Analysis under Complicated Conditions and Mining Method Study[D].Ganzhou：Jiangxi University of Science and Technology,2013.

[2] 周熊超.礦柱回采穩(wěn)定性分析及地壓顯現(xiàn)規(guī)律研究[D].贛州：江西理工大學，2010. Zhou Xiongchao.Pillar Stope Stability Analysis and Pressure Patterns Appeared[D].Ganzhou：Jiangxi University of Science and Technology,2010.

[3] 董 誠，王連捷，楊小聰，等.安慶銅礦地應力測量[J].地質(zhì)力學學報，2001，7(2)：260-263. Dong Cheng，Wang Lianjie，Yang Xiaocong,et al.Stress measurement of Anqing Coper Mine[J].Journal of Geomechanics，2001，7(2)：260-263.

[4] 楊志強，吳福和，楊小聰，等.安慶銅礦不規(guī)則二步驟礦柱采場回采實踐[J].中國礦業(yè)，2011，20(2)：78-80. Yang Zhiqiang，Wu Fuhe，Yang Xiaocong,et al.Practice of irregular jamp exploitation in Anqing Copper Mine[J].China Mining Magazine，2011，20(2)：78-80.

[5] 張海軍，李 濤，衣淑鈺，等.特大型水平礦柱底柱資源回收的技術問題分析及對策[J].有色金屬：礦山部分，2011，63(2)：1-5. Zhang Haijun，Li Tao，Yi Shuyu,et al.Analysis and countermeasure of exploiting the horizontal pillar ore resources[J].Nonferrous Metals：Mining Section，2011，63(2)：1-5.

[6] 趙其禎，郭慧高，張海軍.特大型水平礦柱穩(wěn)定性數(shù)值模擬[J].有色金屬：礦山部分，2008，60(3)：28-31. Zhao Qizhen，Guo Huigao，Zhang Haijun.Numerical simulation on stability of extra large horizontal pillar[J].Nonferrous Metals：Mining Section，2008，60(3)：28-31.

[7] 史啟明.鳳凰山銅礦水平礦柱回采實踐[J].礦業(yè)快報，2001(9)：10-11. Shi Qiming.Mining practice of level pillar in Phoenix Mountain Copper Mine[J].Express Information of Mining Industry，2001(9)：10-11.

[8] 董凱程，吳春平，馮盼學.充填體下水平礦柱開采技術研究[J].中國礦業(yè)，2012，21(9)：103-105. Dong Kaicheng，Wu Chunping，F(xiàn)eng Panxue.Mining technology research on horizontal pillar under backfill[J].China Mining Magazine，2012，21(9)：103-105.

[9] 林 杭，曹 平，等.基于SURPAC的FLAC3D三維模型自動構(gòu)建[J].中國礦業(yè)大學學報，2008，37(3):339-342. Lin Hang，Cao Ping,et al.Automatic generation of FLAC3Dmodel based on SURPAC[J].Journal of China University of Mining and Technology，2008，37(3):339-342.

(責任編輯 石海林)

Study on Stoping Scheme of Horizontal Pillars under the Back-fill

Liao Weicheng1Liu Bo2Zhang Yaoping1,3Xu Youmin1

(1.SchoolofResourcesandEnvironmentEngineering，JiangxiUniversityofScienceandTechnology,Ganzhou341000，China；2.AnhuiLongqiaoMiningCo.，Ltd.，Lujiang231551，China；3.InstituteofAppliedScience，JiangxiUniversityofScienceandTechnology，Ganzhou341000，China)

A horizontal pillar with 15 m of thickness was kept between -385 and -400 m to ensure the stoping safety of high or medium section during production period in Anqing copper mine.With the mining deepening，more and more pillars are reserved.It's going to excavate these horizontal pillars in order to recover copper mine as much as possible.Three mining plans including upward medium length hole stopping with subsequent back-fill，downward medium length hole stopping with subsequent back-fill，and upward drift stopping method are proposed for the horizontal pillar mining according to the real working situation.These three plans are contrasted from its efficiency and safety，especially safety.The upward medium length hole stopping with subsequent back-fill is recommended because of its safety.Four kinds of excavation orders are given for this plan，showing as from the middle to the side，from the left to the right，from the right to the left and from the side to the middle.Then，mine numerical model is created based on FLAC3Dand Surpac software to calculate the four excavation processes at -385 m.The analysis shows that the vertical displacement，principle stress and plastic distribution are minimum by adopting the excavation order from the side to the middle.So it is the optimal scheme.

Horizontal pillar，Stopping scheme，Numerical simulation

2015-03-05

廖偉成(1990—)，男，碩士研究生。

TD853.391

A

1001-1250(2015)-05-061-05