司家營鐵礦急傾斜厚大礦體大規(guī)模空場嗣后充填開采數(shù)值模擬*

翟會超 胡巍巍

(河北鋼鐵集團礦業(yè)公司)

司家營南區(qū)作為國內(nèi)設(shè)計年產(chǎn)1 500萬t級的急傾斜厚大地下大水貧鐵礦，礦區(qū)總長10 km，地質(zhì)總儲量約14.5億t，平均品位30.9%，礦體傾角55°，厚度150～300 m，礦體上部賦存100～120 m厚的第四系含水層，礦區(qū)內(nèi)斷裂構(gòu)造較發(fā)育，是影響礦坑充水的潛在危害，地表又有良田、村莊、河流，給地下鐵礦石安全開采提出更高的要求。采用空場嗣后充填法進行礦石資源回采，一二步礦房依次采用膠結(jié)材料和全尾砂材料進行嗣后充填。大規(guī)模地下開采活動勢必造成地壓規(guī)律變化，進而影響礦體圍巖、礦房礦柱以及充填體自穩(wěn)性，對礦山安全開采構(gòu)成了嚴重威脅[1-2]。為確保礦井安全開采，有必要對礦井大規(guī)模開采的穩(wěn)定性和安全性進行研究。

1 UDEC模型構(gòu)建與參數(shù)取值

針對司家營南區(qū)礦區(qū)斷裂構(gòu)造的發(fā)育特點及巖體力學(xué)特征，本研究采用UDEC軟件對大規(guī)模開采影響范圍內(nèi)的非連續(xù)介質(zhì)進行模擬計算，分析地下資源開采擾動帶來的地壓變化，進而分析規(guī)模化開采對司家營南區(qū)礦體及圍巖穩(wěn)定性的影響。通過對井下礦體進行大結(jié)構(gòu)參數(shù)設(shè)計，將礦塊分為一、二步礦房間隔開采，采用膠結(jié)及全尾砂充填采空區(qū)，礦塊內(nèi)礦房沿垂直礦體走向布置，礦房長度為50～52 m，寬度為20 m，階段高度為50 m，采用平行下向深孔爆破崩礦[3-5]。圖1為UDEC模型為平面二維模型[6-8]，設(shè)計模型X方向長度為0～1 200 m，Y方向長度為0～815 m，模型攻涉及了第四系松散巖、風(fēng)化帶、礦體、圍巖等巖性和7組優(yōu)勢結(jié)構(gòu)弱面，其中，巖石材料物理力學(xué)參數(shù)和結(jié)構(gòu)面力學(xué)參數(shù)分別見表1、表2，膠結(jié)充填材料和全尾砂充填材料的物理力學(xué)參數(shù)見表3。

圖1 節(jié)理模型

表1 巖石材料力學(xué)參數(shù)

表2 節(jié)理面力學(xué)參數(shù)

表3 井下充填體物理力學(xué)參數(shù)

2 數(shù)值計算結(jié)果

根據(jù)分步礦房空場嗣后充填以及膠結(jié)材料、全尾砂材料的充填體特點，重點分析開采活動引起的大范圍地應(yīng)力重分布以及地壓變化引起的位移變化特征。

2.1 應(yīng)力分析

利用UDEC軟件對非連續(xù)介質(zhì)模型中的一步礦房進行開采模擬，分析因一步礦房回采引起的應(yīng)力重分布，結(jié)果見圖2、圖3。

圖2 水平應(yīng)力云圖

圖3 垂直應(yīng)力云圖

分析圖2可知：由于一步礦房回采，使得礦體圍巖應(yīng)力場乃至地表以下礦巖應(yīng)力場發(fā)生了重新分布，形成了一個由二步礦房礦柱支撐的連跨礦房群拱形承載結(jié)構(gòu)作用場，該拱結(jié)構(gòu)存在3個應(yīng)力集中區(qū)域，分別位于回采礦塊(礦體)上盤邊界、礦塊(礦體)下盤邊界、風(fēng)化帶巖層，壓應(yīng)力集中值達到6 MPa，受拱形承載機理影響，回采礦塊垂直上部的第四系地表區(qū)域均顯示出壓應(yīng)力場的分布作用，未出現(xiàn)拉應(yīng)力區(qū)，對于確保回采礦塊上覆巖層穩(wěn)定性有利。由3個應(yīng)力集中區(qū)位置分析可知，該拱結(jié)構(gòu)即為不可逆三鉸拱結(jié)構(gòu)，跨度(礦體水平厚度)L與拱失f近似相等，由于二步回采礦房礦柱的支撐作用使得回采一步礦房直接頂板未出現(xiàn)拉應(yīng)力區(qū)，僅表現(xiàn)出弱壓應(yīng)力，對直接頂板有利。

由圖3可知：最大垂值壓應(yīng)力(20 MPa)作用于礦體上下盤邊界處，易于形成上下盤承壓結(jié)構(gòu)破損，進而導(dǎo)致該不可逆三鉸拱結(jié)構(gòu)發(fā)生破壞，影響頂板及圍巖穩(wěn)定性，甚至引發(fā)冒落，上下盤邊界為危險部位。

在上述分析的基礎(chǔ)上，對一步礦房膠結(jié)充填、二步礦房開采，以及二步礦房回采后全尾砂充填的主應(yīng)力分布情況進行了分析，結(jié)果見圖4、圖5。

圖4 二步礦房充填前主應(yīng)力分布

分析圖4、圖5可知：一步礦房膠結(jié)充填后，對二步礦房礦柱進行回采后出現(xiàn)了礦塊及圍巖全范圍壓應(yīng)力場區(qū)，上下盤處最大壓應(yīng)力高達40～50 MPa，一步膠結(jié)充填體承受的最大壓應(yīng)力為15 MPa，并在上覆頂板存在一顯著拱形支撐結(jié)構(gòu)，拱失f近似等于礦塊的水平寬度L(開采礦體厚度)，表明一步膠結(jié)充填體正與頂板、上下盤圍巖自適性拱形支撐結(jié)構(gòu)協(xié)同維穩(wěn)，須對二步回采礦房及時進行全尾砂快速充填。

圖5 二步礦房充填后主應(yīng)力分布

對比分析二步充填前的應(yīng)力場特點可知，最大主壓應(yīng)力在上下盤處有所降低(35～45 MPa，降幅為10%～12.5%)，最小主壓應(yīng)力在上下盤處則有所增加(10～12 MPa)，對于巖石材料而言，由于二步礦房尾砂充填體的作用使地壓作用得到緩解。但在第四系和風(fēng)化帶交接處衍生出1處高壓應(yīng)力集中區(qū)，以及在風(fēng)化帶內(nèi)形成了長帶形主拉應(yīng)力區(qū)。

2.2 位移分析

分析圖6可知：一步礦房回采后，采場頂板垂直上覆巖層移動最為嚴重，礦塊中心兩側(cè)各80 m，共160 m范圍內(nèi)的地表沉降達到了0.2 m以上，最大值為0.28 m。應(yīng)對該范圍內(nèi)的地表建構(gòu)筑物進行長期觀測，并設(shè)置地表導(dǎo)線測點。

圖6 地表沉降觀測

一步礦房膠結(jié)充填后，礦房頂板位移分布特征如圖7所示。該圖揭示了在膠結(jié)充填體支撐作用的下，由于大規(guī)模礦房群采活動帶來的地應(yīng)力重分布，使得礦塊頂板出現(xiàn)了0.5～2.0 m的不均勻正態(tài)分布沉降。故需提前在礦塊礦房回采前對采場頂板進行支護處理，以防礦房鑿巖硐室、頂板遭受破壞，并確保采空區(qū)充填接頂質(zhì)量。

根據(jù)井下開采引發(fā)的地壓分布特點和位移變化特征，針對礦體上下盤交接處的高應(yīng)力集中區(qū)、采場頂板及地表沉降問題，應(yīng)在礦體上下盤采用中空注漿錨桿聯(lián)合網(wǎng)噴技術(shù)進行支護，在加強承載能力的同時防止礦體與圍巖交接處突發(fā)涌水，運用中空注漿錨索聯(lián)合網(wǎng)噴技術(shù)進行礦塊頂板加固，并充分利用外部加固和礦巖及充填體協(xié)同維穩(wěn)技術(shù)共同確保礦井安全開采。此外，還需對礦井開采影響范圍內(nèi)的地表及地表附著物進行長期跟蹤監(jiān)測。

圖7 垂直位移

3 結(jié) 論

(1)針對司家營大規(guī)模地下鐵礦開采，將優(yōu)勢結(jié)構(gòu)弱面轉(zhuǎn)化為二維平面產(chǎn)狀數(shù)據(jù)信息構(gòu)建的UDEC離散元數(shù)值模型，對于研究群采礦房圍巖地壓分布規(guī)律及其穩(wěn)定性可靠性較強。

(2)礦房群回采過程中，在采場頂板巖層中存在一顯著拱形承壓結(jié)構(gòu)，其拱失近似等于回采水平的總跨度，表明充填體與礦房礦柱、采場頂板、上下盤圍巖實現(xiàn)了協(xié)同維穩(wěn)。

(3)最終高應(yīng)力集中在上下盤處，為40～50 MPa，一步膠結(jié)充填體承壓近15 MPa，且頂板沉降量達0.5～2.0 m，需在礦房回采前對頂板及上下盤交接處進行聯(lián)合支護，以避免礦體圍巖遭受破壞。

(4)二步礦房回采充填后，采場周圍地壓集中得到緩解，但在第四系、風(fēng)化帶交接處及風(fēng)化帶內(nèi)部衍生出1處高壓應(yīng)力集中區(qū)和長帶形主拉應(yīng)力區(qū)，故應(yīng)對地下水和斷層裂隙帶進行監(jiān)測，防范突發(fā)涌水。

(5)針對該礦井下大規(guī)模開采，礦體上部地表與礦房頂板易發(fā)生不均勻沉降的問題，需對采場頂板和地表進行長期觀測。