某石膏礦采空區(qū)坍塌影響數(shù)值模擬研究

張勝光 ，卿自強(qiáng) ，戴軍 ，徐超 ，王淇萱

(1.長(zhǎng)沙礦山研究院有限責(zé)任公司， 湖南 長(zhǎng)沙 410012；2.金屬礦山安全技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室， 湖南 長(zhǎng)沙 410012)

0 引言

房柱法是地下石膏礦產(chǎn)資源的主要采礦方法，房柱法開(kāi)采形成了大面積的懸頂采空區(qū)群[1]。大采空區(qū)隨時(shí)間推移在風(fēng)化、流變等因素作用下會(huì)使采空區(qū)產(chǎn)生大規(guī)模坍塌[2-3]，造成巨大損失，因大面積采空區(qū)失穩(wěn)而導(dǎo)致的地壓災(zāi)害，其后果往往是極其嚴(yán)重的，會(huì)造成重大財(cái)產(chǎn)損失或人員傷亡。2005年11月6日，河北省邢臺(tái)縣尚汪莊石膏礦區(qū)發(fā)生特別重大坍塌事故，造成33人死亡，38人受傷，井下4人失蹤，直接經(jīng)濟(jì)損失774萬(wàn)元[4]。2015年12月25日，山東省臨沂市平邑縣萬(wàn)莊石膏礦區(qū)發(fā)生采空區(qū)坍塌事故，造成1人死亡，13人失蹤，直接經(jīng)濟(jì)損失4133.9萬(wàn)元[5]。

目前，國(guó)內(nèi)一些學(xué)者對(duì)采空區(qū)穩(wěn)定性數(shù)值模擬開(kāi)展了相關(guān)研究工作。郝勇浙等[6]為了解內(nèi)蒙古某金礦豎井的穩(wěn)定情況，采用3DMine-Rhino-FLAC3D耦合建模方法對(duì)內(nèi)蒙古某金礦地表、礦體和采空區(qū)進(jìn)行建模，對(duì)豎井周?chē)膽?yīng)力、位移、塑性區(qū)等方面進(jìn)行了分析，并提出了治理措施。楊勇等[7]以紅嶺多金屬礦復(fù)雜采空區(qū)為例，通過(guò)數(shù)值模擬及K.B.魯別涅依他公式法分別研究了空區(qū)間柱、頂板和上下盤(pán)圍巖的穩(wěn)定性，并與現(xiàn)場(chǎng)實(shí)際情況做對(duì)比分析，得到復(fù)雜采空區(qū)的穩(wěn)定性情況以及潛在破壞區(qū)域。容宇[8]以某銅礦采空區(qū)為例，采用 ANSYS軟件建立采空區(qū)、地層、斷層和巷道單元模型，再利用FLAC3D軟件進(jìn)行計(jì)算，分析了采場(chǎng)開(kāi)挖后-120 m中段及各縱剖面的沉降和彈塑性區(qū)域分布，分析結(jié)果可為后續(xù)采空區(qū)周邊巷道開(kāi)挖提供指導(dǎo)。

本研究以某石膏礦為例，利用二維數(shù)值模擬有限元軟件Phase2，建立了相關(guān)地下采場(chǎng)分析數(shù)值模型，對(duì)-150 m水平上部采空區(qū)坍塌后巖層移動(dòng)和應(yīng)力重新分布趨勢(shì)進(jìn)行了模擬計(jì)算，計(jì)算結(jié)果可為礦山企業(yè)下一步的采空區(qū)治理方案提供參考依據(jù)。

1 工程概況

某石膏礦礦體傾向西南，淺部礦體傾角為25°～40°，深部礦體傾角為 19°～40°，礦體平均傾角為25°～30°，總的趨勢(shì)為由東向西逐漸變陡，淺部及斷層附近常出現(xiàn)次級(jí)褶皺，礦區(qū)范圍內(nèi)由上至下共有Ⅰ、Ⅱ、Ⅲ、Ⅳ4層礦體，主要開(kāi)采Ⅱ、Ⅲ號(hào)礦體。礦山采用斜井開(kāi)拓，建有主副兩條斜井，采用房柱法開(kāi)采，前進(jìn)式回采，留設(shè)條形礦柱支撐頂板。為了減少矸石提升、排棄量，利用井下矸石對(duì)采空區(qū)進(jìn)行部分充填。

-50 m、-83 m、-116 m、-150 m中段已經(jīng)回采結(jié)束并部分區(qū)域垮塌，-183 m中段已開(kāi)采完畢并部分密閉，-216 m、-248 m中段已回采Ⅱ、Ⅲ號(hào)膏層部分礦體，Ⅱ、Ⅲ號(hào)礦體采空區(qū)部分重疊。大量采空區(qū)垮塌致使礦山開(kāi)采條件惡化，會(huì)對(duì)地表、下部中段及井筒產(chǎn)生擾動(dòng)影響，給礦山生產(chǎn)和安全管理帶來(lái)嚴(yán)重影響。因此，對(duì)采空區(qū)坍塌后巖層移動(dòng)和應(yīng)力重分布趨勢(shì)進(jìn)行研究分析是非常必要的，對(duì)保障礦山安全生產(chǎn)具有重要意義。

2 數(shù)值模擬

2.1 數(shù)值計(jì)算模型建立

數(shù)值模擬計(jì)算方法不僅可以分析礦山采空區(qū)坍塌圍巖應(yīng)力、位移的變化趨勢(shì)，還能直觀地展現(xiàn)礦山回采對(duì)構(gòu)筑物的影響機(jī)理和規(guī)律。本次數(shù)值計(jì)算模擬了-150 m水平上部采空區(qū)坍塌后巖層移動(dòng)和應(yīng)力重新分布的趨勢(shì)，分析了礦山采空區(qū)坍塌對(duì)地表、下部中段及井筒的影響。

沿主井做垂直剖面，如圖1所示，此剖面與礦體、巖層走向垂直，兩側(cè)礦體長(zhǎng)度均超過(guò)200 m，符合平面應(yīng)變問(wèn)題前提條件，可采用二維數(shù)值模擬有限元軟件Phase2計(jì)算分析。

以圖1所示剖面為平面模型，建立數(shù)值計(jì)算網(wǎng)格圖，如圖2所示。

圖1 沿主井垂直礦體走向剖面

圖2 數(shù)值計(jì)算網(wǎng)絡(luò)模型

2.2 巖石力學(xué)參數(shù)的選取

根據(jù)數(shù)值模擬計(jì)算分析所需要的基本數(shù)據(jù)，在礦山開(kāi)采現(xiàn)場(chǎng)選取了礦石及圍巖巖樣，主要進(jìn)行了礦巖容重、單軸壓縮條件下礦、巖的力學(xué)特性、劈裂抗拉試驗(yàn)4項(xiàng)實(shí)驗(yàn)室試驗(yàn)。并根據(jù)巖體結(jié)構(gòu)特征（見(jiàn)圖3）、分級(jí)指標(biāo)及室內(nèi)巖石力學(xué)試驗(yàn)結(jié)果，運(yùn)用Roclab軟件計(jì)算分析得出了巖體的力學(xué)參數(shù)，計(jì)算結(jié)果見(jiàn)表1。

圖3 單軸壓縮破壞事件

表1 礦巖力學(xué)參數(shù)試驗(yàn)結(jié)果匯總

3 數(shù)值模擬計(jì)算結(jié)果分析

3.1 位移計(jì)算結(jié)果分析

-150 m水平上部采空區(qū)坍塌后，在-150 m水平上部采空區(qū)地表形成下沉盆地，巖移最大點(diǎn)位于采空區(qū)上方，產(chǎn)生的最大地表沉降為 0.3 m。采空區(qū)被圍巖塌落巖石充填，沿-150 m水平以上坍塌采場(chǎng)至地表巖層位移量逐漸減小，巖層變形量由空區(qū)頂板的3.6 m逐漸減小至地表的0.3 m。巖層移動(dòng)主要發(fā)生在坍塌采空區(qū)上部，下部中段采場(chǎng)圍巖位移量在8 mm以內(nèi)，影響微弱。主井及副井位于礦體底板，-150 m水平上部采空區(qū)坍塌產(chǎn)生的巖層移動(dòng)對(duì)主、副井無(wú)影響。-150 m水平以上采空區(qū)坍塌計(jì)算位移結(jié)果如圖4所示。

圖4 -150 m水平以上采空區(qū)坍塌計(jì)算位移結(jié)果

3.2 應(yīng)力計(jì)算結(jié)果分析

-150 m水平上部采空區(qū)正上方是覆巖的應(yīng)力釋放區(qū)，覆巖剪應(yīng)力、拉應(yīng)力破壞區(qū)均位于坍塌采空區(qū)頂板。上部水平采空區(qū)坍塌造成的-150 m以下礦體圍巖應(yīng)力變化為：最小主應(yīng)力變化范圍在±0.1 MPa；最大主應(yīng)力變化范圍在1.5 MPa以內(nèi)，在已回采結(jié)束的-183 m水平變化最大。因此，-150 m水平上部采空區(qū)坍塌只對(duì)-183 m水平產(chǎn)生較弱的影響，對(duì)目前-216 m、-248 m生產(chǎn)水平及下部-283 m開(kāi)拓水平無(wú)影響。主井、副井在-150 m標(biāo)高以下，-150 m以上采空區(qū)的坍塌不會(huì)對(duì)主井、副井造成影響。-150 m以上采空區(qū)坍塌計(jì)算應(yīng)力結(jié)果如圖5所示。

圖5 -150 m水平以上采空區(qū)坍塌應(yīng)力計(jì)算結(jié)果

4 結(jié)論

（1）石膏礦-150 m水平上部采空區(qū)坍塌后，巖層移動(dòng)主要發(fā)生在坍塌采空區(qū)上部，下部中段采場(chǎng)圍巖位移量在8 mm以內(nèi)，影響微弱，主井及副井位于礦體底板，-150 m水平上部采空區(qū)坍塌對(duì)主、副井無(wú)影響。

（2）石膏礦-150 m水平上部采空區(qū)正上方是覆巖的應(yīng)力釋放區(qū)，覆巖剪應(yīng)力、拉應(yīng)力破壞區(qū)均位于坍塌采空區(qū)頂板，-150 m水平上部采空區(qū)坍塌對(duì)-183 m水平產(chǎn)生較弱的影響，對(duì)目前-216 m、-248 m生產(chǎn)水平及下部-283 m開(kāi)拓水平無(wú)影響。

（3）石膏礦主井、副井在-150 m標(biāo)高以下，-150 m以上采空區(qū)的坍塌不會(huì)對(duì)主井、副井造成影響。

（4）礦山應(yīng)加強(qiáng)對(duì)-183 m水平采空區(qū)的監(jiān)測(cè)，以提前采取空區(qū)失穩(wěn)階治措施。