劉 飛 熊賢亮

(1.安徽省旌德縣安全生產(chǎn)監(jiān)督管理局；2.中鋼集團(tuán)馬鞍山礦山研究院有限公司)

某螢石礦開(kāi)采近20 a，井下形成采空區(qū)較多，原巖應(yīng)力平衡被破壞，圍巖出現(xiàn)了變形、移動(dòng)。隨著繼續(xù)向深部開(kāi)采，圍巖變形、移動(dòng)和破壞將會(huì)進(jìn)一步加劇。因此，必須對(duì)地下開(kāi)采引起的上覆巖層的移動(dòng)規(guī)律及其對(duì)地表移動(dòng)的影響和采場(chǎng)地壓分布規(guī)律進(jìn)行研究[1]。

在特殊條件下，當(dāng)不適宜把主要開(kāi)拓巷道布置在巖石移動(dòng)范圍之外時(shí)，或者對(duì)已投產(chǎn)的礦井在井筒附近發(fā)現(xiàn)新礦體、礦體向下延伸使得井筒落入巖石移動(dòng)范圍之內(nèi)時(shí)，為了保護(hù)井筒及其建筑物，需要設(shè)置保安礦柱[2]。為實(shí)現(xiàn)該礦山的安全生產(chǎn)及可持續(xù)發(fā)展，在回采過(guò)程中，如何科學(xué)合理地確定礦山地下開(kāi)采的巖體移動(dòng)角，是一項(xiàng)亟待解決的問(wèn)題。

1 工程概況

某螢石礦采用平硐-盲豎井開(kāi)拓方式，生產(chǎn)規(guī)模為6萬(wàn)t/a，采用淺孔留礦采礦法和中深孔留礦采礦法，目前已經(jīng)開(kāi)采多年，井下形成較多空區(qū)，而且在生產(chǎn)探礦的過(guò)程中發(fā)現(xiàn)了新的礦體，如果按照傳統(tǒng)的移動(dòng)角圈定開(kāi)采移動(dòng)范圍，為保證盲豎井的安全，必須在Ⅰ#礦體+256，+216 m 及以下中段預(yù)留保安礦柱，保安礦柱礦量達(dá)6.78萬(wàn)t，資源浪費(fèi)嚴(yán)重，按照目前原礦石價(jià)格計(jì)算，直接經(jīng)濟(jì)損失將達(dá)到8100萬(wàn)元，因此，礦山必須對(duì)地表位移變化及井下地壓活動(dòng)規(guī)律進(jìn)行研究，在保證安全的基礎(chǔ)上最大限度地回收資源。

礦山采掘工程均位于花崗閃長(zhǎng)巖侵入體內(nèi)，巖體結(jié)構(gòu)均勻致密，裂隙不發(fā)育，呈完整狀，強(qiáng)度高，但局部有斷層，斷層中有多處可見(jiàn)的正角礫型礦石結(jié)構(gòu)，由于螢石礦本身性脆，受構(gòu)造影響，裂隙十分發(fā)育，礦石相當(dāng)破碎，其穩(wěn)固性較差。根據(jù)該螢石礦巖體實(shí)際的賦存情況，并結(jié)合礦山已有的巖石力學(xué)數(shù)據(jù)，確定該螢石礦的巖石物理力學(xué)指標(biāo)，見(jiàn)表1。

表1 花崗閃長(zhǎng)巖物理力學(xué)指標(biāo)

注：巖石為風(fēng)干狀態(tài)。

2 構(gòu)建筑物的損壞程度與地表變形的關(guān)系

地表構(gòu)建筑物的損壞是地表變形傳遞給構(gòu)筑物基礎(chǔ)而引起的。在井下回采過(guò)程中，地表可能產(chǎn)生下沉、傾斜、曲率、水平移動(dòng)和水平變形，構(gòu)建筑物則隨之產(chǎn)生相似的變形。但由于構(gòu)建筑物具有一定的承受能力，因此,地表變形與構(gòu)建筑物變形存在不一致性，兩者之間的關(guān)系與構(gòu)建筑物基礎(chǔ)的材質(zhì)、規(guī)格以及力學(xué)性質(zhì)有關(guān)。如果構(gòu)建筑物的變形值超過(guò)了規(guī)范要求的允許變形值，構(gòu)建筑物將遭受破壞。

2.1 構(gòu)建筑物破壞等級(jí)評(píng)判標(biāo)準(zhǔn)

我國(guó)目前對(duì)構(gòu)建筑物保護(hù)等級(jí)的劃分尚無(wú)統(tǒng)一規(guī)定，且非煤礦山在這方面研究工作較少。國(guó)內(nèi)有關(guān)的幾個(gè)煤炭礦務(wù)局均各自劃分標(biāo)準(zhǔn)，但評(píng)判標(biāo)準(zhǔn)差別不大，本次研究參照唐山煤研所根據(jù)棗莊和峰峰等礦區(qū)的實(shí)測(cè)資料總結(jié)的構(gòu)建筑物保護(hù)等級(jí)與地表變形值的關(guān)系[3],見(jiàn)表2。

表2 建筑物的保護(hù)等級(jí)及允許變形值

2.2 位移與變形的關(guān)系

目前我國(guó)對(duì)構(gòu)建筑物破壞等級(jí)標(biāo)準(zhǔn)[4]采用水平變形和傾斜變形來(lái)評(píng)判，而采取數(shù)值模擬計(jì)算的地表位移值應(yīng)進(jìn)行轉(zhuǎn)換。地表傾斜及水平變形計(jì)算圖見(jiàn)圖1。

圖1 地表傾斜、水平變形計(jì)算圖

地表傾斜變形i為地表下沉盤地沿某一方向的坡度值，即

(1)

式中,iAB為傾斜變形，mm/m；WA為A點(diǎn)沉降值,mm；WB為B點(diǎn)沉降值,mm；lAB為A、B兩點(diǎn)水平距離,m。

地表水平變形ε為移動(dòng)盆地內(nèi)一線段兩端點(diǎn)的水平移動(dòng)差與此線段長(zhǎng)度之比，即

(2)

式中，ε為地表水平變形，mm/m；UC為C點(diǎn)水平位移值，mm；UD為D點(diǎn)水平位移值，mm；lCD為C、D兩點(diǎn)水平距離，m。

3 巖層移動(dòng)角數(shù)值模擬

3.1 數(shù)值模型建立

傳統(tǒng)的剛體極限平衡計(jì)算分析適用于金屬材料，但對(duì)于巖體這類有塑性變形的材料來(lái)講，太過(guò)于理想化，不能體現(xiàn)巖體變形過(guò)程中應(yīng)力與應(yīng)變的影響，為了更好地分析不同因素對(duì)巖層移動(dòng)角的影響程度，本研究利用FLAC3D對(duì)螢石礦剖面建立數(shù)值模型[5]。該螢石礦礦體走向長(zhǎng)140余m，延深至+216 m 中段，平均厚5 m。為了減少模型邊界對(duì)計(jì)算結(jié)果的影響，模型邊界距開(kāi)挖邊界選取3倍左右的采場(chǎng)大小。模型高430 m左右，寬230 m，長(zhǎng)200 m，幾何計(jì)算模型見(jiàn)圖2。

由于地表第四系厚度較小，相對(duì)于模型可忽略不計(jì)，故本次模型未考慮地表第四系的影響。模型按不同的礦巖性質(zhì)及開(kāi)挖順序進(jìn)行分組：圍巖、開(kāi)采礦體、暫不開(kāi)采礦體，共包括90119個(gè)單元，20675

圖2 幾何計(jì)算模型

個(gè)節(jié)點(diǎn),三維計(jì)算模型見(jiàn)圖3。

圖3 三維計(jì)算模型

3.2 開(kāi)采方案模擬計(jì)算

只考慮自重應(yīng)力，模型Z方向的初始應(yīng)力場(chǎng)見(jiàn)圖4?？芍?，數(shù)值模擬計(jì)算的初始地應(yīng)力場(chǎng)分布與自重地應(yīng)力場(chǎng)很接近，僅在礦巖接觸區(qū)域(網(wǎng)格劃分不均勻造成)與端部存在地應(yīng)力等值線起伏和應(yīng)力集中現(xiàn)象，對(duì)模型計(jì)算域結(jié)果影響不大，可以滿足計(jì)算的要求。

圖4 垂直方向初始應(yīng)力云圖

對(duì)模型進(jìn)行數(shù)值模擬，將開(kāi)采礦體部分進(jìn)行開(kāi)挖(圖5)，垂直方向的應(yīng)力分布云圖見(jiàn)圖6?？芍?，在垂直方向上，主要的應(yīng)力是由巖體的自重產(chǎn)生，由于開(kāi)采擾動(dòng)的影響，破壞原有應(yīng)力平衡，致使應(yīng)力重新分布，達(dá)到最終平衡，采空區(qū)周邊縱向應(yīng)力為2～6 MPa。

模型共布置14個(gè)監(jiān)測(cè)點(diǎn)，各監(jiān)測(cè)單元水平開(kāi)挖后間距為10 m。監(jiān)測(cè)點(diǎn)布置見(jiàn)圖7。

由于監(jiān)測(cè)點(diǎn)較多，僅列出具有代表性的10#、11#監(jiān)測(cè)點(diǎn)豎直及水平位移曲線，見(jiàn)圖8、圖9。各個(gè)監(jiān)測(cè)點(diǎn)水平及豎直位移統(tǒng)計(jì)見(jiàn)表3。

4 結(jié) 論

根據(jù)礦山巖石力學(xué)性質(zhì)，通過(guò)數(shù)值模擬的方法建立模型，可以判斷出礦體開(kāi)采已經(jīng)造成沉降，影響到了地表建構(gòu)筑物。數(shù)值模擬計(jì)算結(jié)果表明， 1#～10#監(jiān)測(cè)點(diǎn)豎直及水平位移均小于Ⅰ級(jí)保護(hù)類型的允許變形值，其中10#監(jiān)測(cè)點(diǎn)豎向位

圖5 開(kāi)挖計(jì)算模型

圖6 開(kāi)挖后垂直方向應(yīng)力云圖

圖7 開(kāi)采位移監(jiān)測(cè)點(diǎn)布置

圖8 10#監(jiān)測(cè)點(diǎn)位移變化情況

圖9 11#監(jiān)測(cè)點(diǎn)位移變化情況

移為-43.8 mm，水平位移為-15.8 mm，傾斜變形為2.18 mm/m，水平變形為0.96 mm/m，在允許變形范圍內(nèi)；11#監(jiān)測(cè)點(diǎn)豎向位移為-79.5 mm，水平位移為-34.2 mm，傾斜變形為3.57 mm/m，水平變形為1.84 mm/m，均大于允許變形值。結(jié)合我國(guó)建構(gòu)筑物破壞等級(jí)的標(biāo)準(zhǔn)對(duì)受開(kāi)采影響的建構(gòu)筑物進(jìn)行評(píng)判，模型中10#監(jiān)測(cè)點(diǎn)對(duì)應(yīng)的開(kāi)采移動(dòng)角為77°，即礦山最大允許的開(kāi)采移動(dòng)角應(yīng)小于77°。礦山上、下盤圍巖為同一巖性，故確定礦山上、下盤圍巖開(kāi)采最大移動(dòng)角均為77°。在保證盲豎井安全的前提下，為進(jìn)一步確保礦山生產(chǎn)安全，該螢石礦開(kāi)采移動(dòng)角確定為75°。

表3 各監(jiān)測(cè)點(diǎn)水平及豎直位移統(tǒng)計(jì)

注：移動(dòng)角數(shù)據(jù)由圖7測(cè)量所得。

[1] 高延法,張慶松.礦山巖體力學(xué)[M].徐州：中國(guó)礦業(yè)大學(xué)出版社,2000.

[2] 解世俊.金屬礦床地下開(kāi)采(修訂版)[M].北京：冶金工業(yè)出版社,1979.

[3] 河南理工大學(xué).峰峰集團(tuán)有限公司萬(wàn)年礦村莊下采煤采動(dòng)損害判定與研究[R].焦作：河南理工大學(xué)，2004.

[4] 杜維吾,劉寶琛.金屬礦山地表移動(dòng)與變形規(guī)律[C]//第五屆全國(guó)礦山巖體力學(xué)學(xué)術(shù)會(huì)議論文集.北京：中國(guó)有色金屬學(xué)會(huì)采礦學(xué)術(shù)委員會(huì)巖體力學(xué)專業(yè)委員會(huì)，1992：35-38.

[5] 麻鳳海,王永嘉.地層沉降控制的可變形離散單元模擬[J].巖石力學(xué)與工程學(xué)報(bào)，1999(2):176-179.