高峰礦深部多空區(qū)下充填采礦地壓分布特征*

馬生徽 王文杰 鄧金燦 王越嶺 馬雄忠

(1.武漢科技大學(xué)資源與環(huán)境工程學(xué)院;2.廣西高峰礦業(yè)有限責(zé)任公司)

隨著金屬礦山開采深度不斷延伸，礦山進(jìn)入深部開采面臨著“三高一擾”的復(fù)雜環(huán)境:深部地壓活動(dòng)干擾愈加頻繁，導(dǎo)致深部采場(chǎng)與巷道的維護(hù)也變得困難，造成高損失、高貧化和高安全事故的技術(shù)難題。高峰礦在100號(hào)礦體延深開采以及105號(hào)礦體的開采中，由于開采深度增加以及原有民采及探礦過(guò)程中的亂采亂挖，礦巖的完整性受到很大破壞，出現(xiàn)許多不規(guī)則空區(qū)和大面積的破碎礦巖體，給礦床正常開發(fā)利用帶來(lái)了極大的安全隱患［1］。雖已對(duì)空區(qū)進(jìn)行了充填處理，但隨著開采深度的增加，地壓顯現(xiàn)將日漸顯著，其上部采空區(qū)充填體的穩(wěn)定性將受到影響，且直接危害下部開采的技術(shù)經(jīng)濟(jì)指標(biāo)和作業(yè)安全。因此有必要結(jié)合高峰礦現(xiàn)有地質(zhì)資料，利用ANSYS－Workbench建立高峰礦開采數(shù)值模型，導(dǎo)入FLAC3D計(jì)算分析，研究高峰礦深部多空區(qū)下充填采礦的地壓分布特征［2］，并針對(duì)高峰礦生產(chǎn)管理中存在的問(wèn)題，提出合理地壓控制措施，避免地質(zhì)災(zāi)害的發(fā)生，從而降低礦山安全事故和提高經(jīng)濟(jì)效益。

高峰大廠礦田西礦帶生物礁灰?guī)r中的100號(hào)礦體及深部延伸部分的垂向深度大于1 000 m，－79 m標(biāo)高以上稱為100號(hào)礦體，－79 m以下為100號(hào)礦體的延伸部分稱為105號(hào)礦體，礦體從標(biāo)高+690至－50 m，傾角由陡變緩，礦體由大變小;從標(biāo)高－50 m以下，傾角由緩變陡，礦體由小變大，礦體平面投影呈馬蹄狀。105號(hào)礦體總體走向近南北，總體傾向于東至南東，傾角陡，礦體在深部往南至南西方向側(cè)伏，礦體走向長(zhǎng)100～300 m，水平厚度10～70 m，礦體由致密塊狀礦石組成，與圍巖界線清楚，普氏系數(shù)f為8～10，礦巖構(gòu)造上既有褶皺和一系列斷層的切割，也有花崗巖脈的穿插，也形成一定范圍的破碎巖體。從采空區(qū)的分布來(lái)看，民采盜采形成的采空區(qū)主要分布在－166 m水平以上，且頂板暴露面積在2 000 m2以上，是影響安全開采的主因。－166 m水平以下的采空區(qū)規(guī)模較小，頂板暴露面積在2 000 m2以下，且采空區(qū)內(nèi)都留有一定的礦柱，采空區(qū)相互獨(dú)立處于相對(duì)穩(wěn)定狀態(tài)，加之高峰礦近幾年大量投入對(duì)采空區(qū)進(jìn)行充填處理，使采空區(qū)周圍地壓活動(dòng)趨于穩(wěn)定。

1 深部開采數(shù)值模擬

1.1 計(jì)算模擬的建立

由于礦體從地表淺部向下延伸超過(guò)1 000 m，其走向、傾向、傾角、各水平截面都是逐漸變化的，為世界罕見(jiàn)的巨型囊狀－透鏡狀特富礦體，復(fù)雜的礦體形態(tài)難以在FLAC3D中直接建模，因此利用ANSYSWorkbench建立高峰礦開采模型，劃分網(wǎng)格后，通過(guò)轉(zhuǎn)換程序?qū)隖LAC3D有限差分程序中進(jìn)行計(jì)算分析［3-4］。劃分網(wǎng)格時(shí)選用適應(yīng)性較強(qiáng)的四面體網(wǎng)格，模型如圖1所示，共380 653個(gè)單元，518 890個(gè)節(jié)點(diǎn)，模型主要分為圍巖、礦體和采空區(qū)。

圖1 導(dǎo)入FLAC3D后的計(jì)算模型

1.2 力學(xué)參數(shù)的選取

數(shù)值模擬中主要考慮圍巖、礦體、充填體都為彈塑性材料，因此采用摩爾－庫(kù)侖屈服準(zhǔn)則，計(jì)算中把100號(hào)礦體和105號(hào)礦體的巖體參數(shù)區(qū)分開來(lái)，100號(hào)巖體的相關(guān)參數(shù)由高峰礦提供，105號(hào)礦體的相關(guān)參數(shù)在巖石力學(xué)實(shí)驗(yàn)的基礎(chǔ)上，通過(guò)巖體的RMR值與巖體參數(shù)之間關(guān)系轉(zhuǎn)化所得，礦巖物理參數(shù)如表1所示。計(jì)算中對(duì)于模型開挖后的空區(qū)選用FLAC3D內(nèi)置的空模型(Null)，對(duì)于開挖前以及開挖后的非空區(qū)部分都采用摩爾－庫(kù)侖塑性模型。

1.3 邊界條件與計(jì)算方案

根據(jù)高峰礦歷次地應(yīng)力測(cè)量的結(jié)果，進(jìn)行回歸分析后，得出地應(yīng)力與深度H之間的關(guān)系如式(1)所示。

表1 數(shù)值計(jì)算中礦巖的物理參數(shù)

按照快速應(yīng)力邊界法施加初始應(yīng)力［5］，由于地表起伏，取地表的平均標(biāo)高為模型上表面標(biāo)高，采用快速應(yīng)力邊界法施加初始條件并計(jì)算平衡后，進(jìn)行后續(xù)的開挖計(jì)算時(shí)要對(duì)模型施加位移邊界條件，即固定模型底部Z向位移，模型左右表面X向位移，模型前后表面Y向位移。

為了便于在FLAC3D里計(jì)算分析，對(duì)高峰礦開采模型進(jìn)行了適當(dāng)簡(jiǎn)化，忽略節(jié)理和裂隙等不連續(xù)面的影響，也不考慮巖體地下水、地震和爆破震動(dòng)等因素。計(jì)算中采用了各中段分步回采的計(jì)算方式，即開挖充填+690～+540 m、+540～+400 m、+400～+250 m、+250～+100 m、+100～－50 m和開挖充填所有空區(qū)，總共6步回采順序，每個(gè)階段的計(jì)算均按順序在前一階段計(jì)算基礎(chǔ)上連續(xù)進(jìn)行，從而客觀地反映了礦體逐步開采過(guò)程中，巖體應(yīng)力累積、變形延續(xù)與破壞逐漸發(fā)展的力學(xué)進(jìn)程，最后得出高峰礦在多空區(qū)群回采充填后的力學(xué)效應(yīng)，用以分析深部地壓分布特征。考慮礦山實(shí)際生產(chǎn)情況，數(shù)值模擬結(jié)果主要分析－151 m和－200 m 2個(gè)主要回采中段水平的地壓分布狀態(tài)以及典型多空區(qū)群開挖充填對(duì)地壓分布的影響。

2 深部地壓分布特征分析

2.1 －151 m水平應(yīng)力分布狀態(tài)

圖2(a)所示為多空區(qū)群充填前－151 m水平最大主應(yīng)力云圖(應(yīng)力符號(hào):壓應(yīng)力為負(fù)，拉應(yīng)力為正)，應(yīng)力值在－15.0～－59.8 MPa，圖中空白區(qū)域?yàn)?7號(hào)空區(qū)，承壓區(qū)主要分布在采空區(qū)四周、礦體北部以及脈外巷區(qū)域，19號(hào)采空區(qū)底板恰好位于礦體北部應(yīng)力集中區(qū)，可見(jiàn)承壓區(qū)形成的主因是空區(qū)開挖后應(yīng)力重新分布而在采空區(qū)四周出現(xiàn)次生應(yīng)力集中，承壓區(qū)中的礦巖因過(guò)度應(yīng)力集中易發(fā)生破壞，這與現(xiàn)場(chǎng)巖石條件相吻合。卸壓區(qū)主要分布在礦體上盤中部與下盤沿礦體走向區(qū)域，特別是在礦體下盤出現(xiàn)1條明顯的卸壓帶。圖2(b)為空區(qū)充填后－151 m水平最大主應(yīng)力云圖，應(yīng)力集中區(qū)保持不變而集中值明顯下降，其他區(qū)域應(yīng)力分布狀態(tài)基本與充填前一致，說(shuō)明充填空區(qū)對(duì)局部的應(yīng)力集中有所改善，但對(duì)整體的應(yīng)力分布狀態(tài)影響不大。

圖2 －151 m水平最大主應(yīng)力分布

2.2 －200 m水平應(yīng)力分布狀態(tài)

圖3(a)為空區(qū)充填前－200 m水平最大主應(yīng)力云圖，應(yīng)力值在－36.0～－52.6 MPa，承壓區(qū)主要在礦體北部和南部花崗斑巖附近。在垂直礦體走向上，應(yīng)力分布狀態(tài)表現(xiàn)為應(yīng)力值先減小，后逐漸增大，最后趨于該水平平均水平。由于受多空區(qū)群影響較小，在－200 m中段水平承壓區(qū)都處于礦體上盤，卸壓區(qū)位于礦體下盤，在礦體下盤形成1條明顯的卸壓帶。圖3(b)為空區(qū)充填后－200 m水平最大主應(yīng)力云圖，與空區(qū)開挖后基本趨于一致。在南部花崗斑巖附近的承壓區(qū)內(nèi)井巷工程較為密集，水倉(cāng)、變電站等也布置在此，出現(xiàn)的局部應(yīng)力集中造成巖體較破碎，支護(hù)量大、成本高，巷道變形劇烈，應(yīng)盡量加強(qiáng)該區(qū)域圍巖的監(jiān)測(cè)以及巷道的支護(hù)。

圖3 －200 m水平最大主應(yīng)力分布

2.3 采空區(qū)群對(duì)地壓分布影響

由于深部空區(qū)群眾多，這里選取各采空區(qū)中相鄰空間區(qū)間頂?shù)装搴穸刃『捅┞睹娣e大的空區(qū)群作為研究對(duì)象［6-8］。圖4為2號(hào)、4號(hào)、15號(hào)和23號(hào)共4個(gè)采空區(qū)的空間位置圖，最上面的2號(hào)空區(qū)位于105號(hào)礦體的－79～－103 m標(biāo)高，頂板暴露面積為3 632 m2，4號(hào)空區(qū)在－106～－110 m標(biāo)高，并與15號(hào)、23號(hào)空區(qū)部分重疊。

圖4 典型多空區(qū)群相對(duì)位置

取典型多空區(qū)群的縱剖面進(jìn)行分析。圖5(a)為計(jì)算截面上充填前的采空區(qū)群最大主應(yīng)力分布狀態(tài)，可以看出在采空區(qū)群的頂?shù)装鍏^(qū)域?yàn)樾秹簠^(qū)，易出現(xiàn)拉應(yīng)力集中，拉應(yīng)力的存在與擴(kuò)展對(duì)采空區(qū)穩(wěn)定性構(gòu)成較大威脅。在采空區(qū)兩側(cè)形成承壓區(qū)，并逐漸向外延伸，應(yīng)力最大值為－69.9 MPa，巖體也容易發(fā)生壓縮破壞，承壓區(qū)主要在15號(hào)空區(qū)與23號(hào)空區(qū)兩側(cè)。圖5(b)為充填后采空區(qū)群最大應(yīng)力分布狀態(tài)，可以發(fā)現(xiàn)卸載區(qū)和承載區(qū)的位置未發(fā)生變化，但是承載區(qū)的應(yīng)力集中值有所下降，說(shuō)明充填后有利于巖體中的能量向充填體轉(zhuǎn)移，有效減少次生應(yīng)力集中。

圖5 各空區(qū)的最大主應(yīng)力分布

從采空區(qū)的地質(zhì)調(diào)查來(lái)看，前期民采形成的空區(qū)群底板標(biāo)高最低為－180 m水平，且暴露面積大的空區(qū)群在－166 m水平以上，結(jié)合－151 m水平和－200 m水平的地壓分布特征，可以發(fā)現(xiàn)多空區(qū)群的存在對(duì)－151 m水平的地壓分布影響較大，導(dǎo)致局部形成了明顯的承壓區(qū)和卸壓區(qū)，上部多空區(qū)群對(duì)－200 m水平地壓分布影響較小。

2.4 深部地壓控制措施建議

高峰礦深部地壓分布特征為:在沿礦體走向上，上盤形成1條承壓帶，下盤形成1條卸壓帶，受多空區(qū)的影響，在空區(qū)四周也形成次生應(yīng)力集中。而高峰礦大量采準(zhǔn)開拓工程布置在承壓帶上，造成采場(chǎng)與巷道頂板變形大，支護(hù)困難，地壓顯現(xiàn)頻繁，應(yīng)加強(qiáng)該區(qū)域采場(chǎng)與巷道的監(jiān)測(cè)，提高支護(hù)等級(jí)，井巷工程應(yīng)盡量布置在卸壓帶上，并采用卸壓、錨桿支護(hù)等多種方式維護(hù)頂板的穩(wěn)定［9-11］。原則上應(yīng)先回采承壓區(qū)內(nèi)的礦體，后回采卸壓區(qū)內(nèi)的礦體，采用前進(jìn)式回采順序，以控制地壓有序的釋放，防止中間礦帶始終處于應(yīng)力集中狀態(tài)而導(dǎo)致巷道與炮孔易發(fā)生變形與破壞。采場(chǎng)應(yīng)強(qiáng)采強(qiáng)出強(qiáng)充，減少頂板暴露面積與時(shí)間。因此，在合理回采順序的前提下，優(yōu)化井巷工程布置，加強(qiáng)目前承壓區(qū)內(nèi)井巷工程的維護(hù)是保證高峰礦深部安全高效回采的關(guān)鍵。

3 結(jié)論

數(shù)值模擬研究表明，深部空區(qū)群同時(shí)開挖充填后，對(duì)離多空區(qū)較近水平地壓分布影響大，對(duì)離多空區(qū)群較遠(yuǎn)的深部水平地壓分布影響較小，充填開采在高峰深部礦體上下盤產(chǎn)生明顯的承壓區(qū)與卸壓區(qū)，而高峰礦大量井巷工程布置在上盤承壓帶上，這成為導(dǎo)致高峰礦深部地壓活動(dòng)頻繁的主要因素。針對(duì)高峰礦深部地壓分布特征，采用合理的回采順序與采場(chǎng)結(jié)構(gòu)參數(shù)，優(yōu)化布置井巷工程等措施能有效控制地壓災(zāi)害的發(fā)生，使地壓活動(dòng)規(guī)律能為生產(chǎn)作業(yè)服務(wù)，保證高峰礦深部回采安全高效地進(jìn)行。

［1］ 顏榮貴，方建勤，曹 陽(yáng)，等.耦合于暗空?qǐng)龅牡V山地壓演化規(guī)律研究［J］.金屬礦山，2004(S.):53-59.

［2］ 李 輝，陳廣平，馬海濤，等.三維數(shù)值模擬在礦山地壓災(zāi)害評(píng)價(jià)中的應(yīng)用［J］.金屬礦山，2011(8):136-139.

［3］ 劉心庭，唐輝明.FLAC3D復(fù)雜網(wǎng)格模型的構(gòu)建及其工程應(yīng)用［J］.金屬礦山，2010(11):108-111.

［4］ 廖秋林，曾錢幫，劉 彤，等.基于ANSYS平臺(tái)復(fù)雜地質(zhì)體FLAC3D模型的自動(dòng)生成［J］.巖石力學(xué)與工程學(xué)報(bào)，2005，24(6):1010-1013.

［5］ 李仲奎，戴 榮，姜逸明，等.FLAC3D分析中的初始應(yīng)力場(chǎng)生成及在大型地下洞室群計(jì)算中的應(yīng)用［J］.巖石力學(xué)與工程學(xué)報(bào)，2002，21(S2):2387-2392.

［6］ 黃應(yīng)盟，鄧金燦，毛建華，等.高峰礦區(qū)100號(hào)礦體采空區(qū)治理與地壓監(jiān)測(cè)研究［J］.礦業(yè)研究與開發(fā)，2005，25(3):18-19.

［7］ 唐紹輝，黃英華，潘 懿，等.災(zāi)害性群空區(qū)深部開采采場(chǎng)結(jié)構(gòu)參數(shù)優(yōu)化研究［J］.中國(guó)礦業(yè)，2009，18(11):65-68.

［8］ 王湖鑫，陳 何，吳志安，等.赤峰國(guó)維礦多空區(qū)復(fù)雜條件下殘礦資源回采技術(shù)研究［J］.中國(guó)礦業(yè)，2011，20(3):69-70.

［9］ 王文杰.高應(yīng)力區(qū)巷道卸壓支護(hù)技術(shù)及應(yīng)用［J］.金屬礦山，2010(3):23-25.

［10］ 王文杰.高應(yīng)力區(qū)卸壓與錨桿支護(hù)對(duì)巷道穩(wěn)定性影響分析［J］.金屬礦山，2010(10):1-5.

［11］ 馬生徽，王文杰.基于FLAC3D的管縫錨桿與玻璃鋼錨桿支護(hù)效果對(duì)比［J］.金屬礦山，2012(6):9-12.