胡 超 丁成功 金長宇 東龍賓 陳五九

（1.東北大學(xué)深部金屬礦山安全開采教育部重點實驗室，遼寧 沈陽 110819；2.中鋼集團馬鞍山礦山研究總院股份有限公司，安徽 馬鞍山 243000；3.中冶北方工程技術(shù)有限公司，遼寧 大連 116600；4.安徽馬鋼張莊礦業(yè)有限責(zé)任公司，安徽 六安 237484）

隨著國家對礦山開采造成的地表塌陷、土地占用等環(huán)境問題越來越重視，空場嗣后充填采礦法逐漸成為金屬礦床地下開采的重要選擇［1-4］。為了增加礦石產(chǎn)量降低礦石回采成本，空場嗣后充填法的階段高度呈現(xiàn)出不斷上升的趨勢，同時采場穩(wěn)定作為地下礦床開采的前提條件，又在一定程度上限制了采場階段高度的提升。采場穩(wěn)定性涉及的影響因素眾多，使得采場穩(wěn)定性分析存在極大困難。目前國內(nèi)外學(xué)者主要采用相似材料試驗、數(shù)值模擬、理論分析、現(xiàn)場監(jiān)測等［5-9］方法，對地下采場礦柱圍巖的應(yīng)力變形及破壞規(guī)律進行研究。杜翠鳳等［10］采用現(xiàn)場鉆孔監(jiān)測與相似材料試驗相結(jié)合的方法，得到了無底柱分段崩落法在開采過程中頂板圍巖的變形規(guī)律；宋衛(wèi)東等［11］、徐文彬等［12］采用數(shù)值模擬方法研究了階段空場嗣后充填采場圍巖的破壞規(guī)律和采空區(qū)破壞機理；吳迪等［13］、李剛等［14］采用相似材料試驗方法研究了階段空場嗣后充填采場頂板和間柱的變形規(guī)律；徐偉蘭等［15］采用數(shù)值模擬方法對大直徑深孔階段空場嗣后充填采場穩(wěn)定性進行了研究；吳亞斌等［16］、吳振坤等［17］通過數(shù)值模擬研究了礦房回采充填順序?qū)Σ蓤龇€(wěn)定性的影響；張海波等［18］、王龍等［19］采用理論分析和數(shù)值模擬相結(jié)合的方法，研究了空場嗣后充填采場在不同礦柱間距情況下礦巖應(yīng)力和位移變化規(guī)律。以上學(xué)者采用相似材料試驗、數(shù)值模擬、理論分析、現(xiàn)場監(jiān)測等方法對不同礦山進行了研究，實踐證明以上方法均能為礦山安全生產(chǎn)提供理論指導(dǎo)。目前，對于急傾斜厚大礦體大尺寸階段空場嗣后充填開采使用數(shù)值模擬方法較多，在采場穩(wěn)定性方面很少使用相似材料試驗這種直觀的研究方法。礦山開采環(huán)境各異，任何方法在面對不同礦山時均可能有不足之處，理論分析簡單易行，但對于復(fù)雜地質(zhì)因素的討論過于理想化，理論模型過度簡化會導(dǎo)致較大誤差；現(xiàn)場監(jiān)測具有真實可靠的優(yōu)點，但監(jiān)測工作量大、周期長、成本高；數(shù)值模擬具有豐富的本構(gòu)模型，能夠求解各種復(fù)雜地質(zhì)條件中的工程問題，在工程穩(wěn)定性的判斷界線上不夠清楚；相似材料試驗將實際工程縮小到模型中，對模型按照實際工程方案進行模擬，能夠直觀反映采場應(yīng)力、變形規(guī)律［20］，但相似材料試驗的定量描述不夠精確。因此，本研究采用相似材料試驗和數(shù)值模擬相結(jié)合的方法，對急傾斜厚大礦體礦房階段高度提升后新中段采場內(nèi)礦柱和圍巖的應(yīng)力、變形變化規(guī)律進行研究，通過模型直觀判斷采場破壞情況和穩(wěn)定性，結(jié)合應(yīng)力監(jiān)測和變形監(jiān)測，研究頂板變形規(guī)律、礦柱破壞位置和破壞模式，采用數(shù)值模擬方法分析礦柱寬度和礦房回采充填順序?qū)Σ蓤鰬?yīng)力及變形規(guī)律的影響。

1 工程概況

礦區(qū)位于霍邱縣城西北31 km，地處周集鎮(zhèn)和馮井鎮(zhèn)之間，北距淮河10 km、阜陽45 km。礦區(qū)東側(cè)為G105國道，西側(cè)有濟廣高速和阜陽—六安鐵路，向南與312國道、寧西鐵路相連，可至六安、合肥、南京等地，交通十分方便。礦床類型為大型鞍山式沉積變質(zhì)鐵礦床，礦體埋藏在146～196 m厚的第四系覆蓋層之下，屬隱伏礦床。第四系之下為基巖古風(fēng)化帶，厚度為20～70 m，一般為40 m左右，接近古風(fēng)化帶礦石為氧化礦，氧化礦帶厚度為28～61 m。氧化帶以下為原生磁鐵礦體，主要為石英磁鐵礦，礦體較堅硬、完整。礦體頂板主要為黑云片巖，巖石堅硬，整體性好，底板主要為角閃斜長片麻巖、黑云斜長片麻巖，頂?shù)装宸€(wěn)固性中等。礦體走向長2 300 m，分布在2～24號勘探線之間，最大埋深近700 m，傾角30°～70°，自西向東由I、II、III號礦體組成。

目前-450 m中段的回采工作已經(jīng)接近尾聲，需要對新中段進行開拓，擬將-390 m中段的階段礦房高度從上中段的60 m提升至90 m，階段礦房高度在原來的基礎(chǔ)上增加了50%，采場應(yīng)力變形規(guī)律及穩(wěn)定性有待研究。

2 相似材料試驗

2.1 相似理論

相似材料實質(zhì)是將工程按照一定比例縮小制作成模型，并有相似理論作為基礎(chǔ)，模型與原型之間相同量綱參數(shù)具有相同相似比，無量綱參數(shù)相似比為1，保證了模型和原型之間的相似性。幾何相似比CL、容重相似比Cγ、摩擦系數(shù)相似比Cf、內(nèi)摩擦角相似比Cφ、泊松比相似比Cμ、應(yīng)變相似比Cε、彈性模量相似比CE、應(yīng)力相似比Cσ之間滿足如下關(guān)系：

試驗選取3個盤區(qū)礦房設(shè)計相似材料試驗，盤區(qū)礦房長100 m，礦柱寬18 m，盤區(qū)內(nèi)有3個寬度為15 m的一步膠結(jié)充填礦房和2個寬度為18.5 m的二步全尾砂充填礦房，模型長度對應(yīng)工程長為480 m；豎直方向上為標(biāo)高在-460～-240 m之間的區(qū)域，包含-450 m中段正在回采階段高度為60 m的采區(qū)和-390 m擬采中段階段高度為90 m的采區(qū)，總高220 m。試驗臺尺寸為2 400 mm×300 mm×1 100 mm，根據(jù)相似原理，確定幾何相似常數(shù)CL=1∶200，容重相似常數(shù)Cγ=1∶2，應(yīng)力相似常數(shù)Cσ=1∶400。

2.2 相似材料配比

試驗為二維平面相似材料試驗，整個模型材料主要為礦體、膠結(jié)充填體和尾砂充填體3種材料，選用EPDM橡膠顆粒作為尾砂充填體的相似材料，珍珠棉為膠結(jié)充填體的相似材料，河沙、重晶石粉和石蠟為礦體的相似材料成分。根據(jù)相似原理，通過試驗，測試不同配比條件下材料的力學(xué)參數(shù)，得到的礦體相似材料配比見表1。

2.3 模型制備

按照河沙∶重晶石粉∶石蠟=8∶2∶1的配比進行試驗?zāi)Ｐ投呀?，制備過程如圖1所示。將河沙和重晶石粉加入攪拌機中進行混合預(yù)熱，預(yù)熱完成后倒入融化的石蠟攪拌均勻，將制備好的材料裝填到試驗臺上，在設(shè)計的應(yīng)力監(jiān)測位置預(yù)埋應(yīng)力傳感器，并壓實模型。模型養(yǎng)護完成后拆除擋板，在模型上確定礦體界線，并在模型上粘貼非編碼點，在試驗臺上設(shè)置編碼點。

2.4 試驗監(jiān)測方案

試驗采用光學(xué)三維攝影測量系統(tǒng)和靜態(tài)應(yīng)力采集系統(tǒng)兩種監(jiān)測方法，對模型在回采充填過程中頂板圍巖變形和頂板、礦柱應(yīng)力進行監(jiān)測。在模型上布置了變形監(jiān)測區(qū)域和4個應(yīng)力監(jiān)測點，如圖2所示。1、2步回采充填為1號礦房，3、4步回采充填為2號礦房，5、6步回采充填為3號礦房，7、8步回采充填為4號礦房，9、10步回采充填為5號礦房，11、12步回采充填為6號礦房。試驗中盤區(qū)礦房回采充填順序與現(xiàn)場施工方案一致，即1～6號礦房依次回采充填。

2.5 充填回采試驗

因試驗臺尺寸限制，模型頂部對應(yīng)實際工程中埋深H為230 m的位置，礦區(qū)以自重應(yīng)力場為主，物理模型無法模擬的上部巖層，使用試驗臺上方加載系統(tǒng)進行加壓等效上覆巖層產(chǎn)生的壓應(yīng)力，模型兩側(cè)采用位移限制條件，磁鐵礦平均密度ρ為3 598.57 kg/m3，應(yīng)力相似比C為400，試驗機加載壓頭面積S為0.18 m2，通過公式（3）、（4）計算可得到每個加載壓頭的等效壓力。

式中，σ′為模型加載應(yīng)力，MPa；F為每個加載壓頭的加載力，N。

試驗開始前對模型加載，所有壓頭加載至3.725 kN保持穩(wěn)定，使用位移監(jiān)測系統(tǒng)記錄原始場，開啟應(yīng)力采集儀采集應(yīng)力數(shù)據(jù)，然后對模型進行回采充填試驗，加載保持穩(wěn)定直至試驗結(jié)束后卸載。

圖3（a）和圖3（b）分別為一步膠結(jié)充填和二步全尾砂充填過程，礦房內(nèi)3個一步膠結(jié)充填礦房同時開挖，形成的空區(qū)一次充填完成，膠結(jié)充填完成后對礦房內(nèi)2個二步尾砂充填礦房進行回采，開挖結(jié)束形成空區(qū)后，使用有機玻璃板封堵模型兩側(cè)，然后進行散體充填。

按照圖3所示的兩步驟回采充填方案依次對1～6號礦房進行了回采充填，試驗過程持續(xù)采集應(yīng)力監(jiān)測點的數(shù)據(jù)，并在各中段回采充填結(jié)束后對位移場進行了記錄。-450 m中段回采充填完成和-390 m中段回采充填完成后的試驗結(jié)果，如圖4所示。分析可知：-450 m中段回采充填完成后盤區(qū)礦柱完整無明顯破壞現(xiàn)象，采場穩(wěn)定性良好。-390 m中段在5號礦房的回采充填過程中，左礦柱在-370 m水平出現(xiàn)開裂破壞，在6號礦房的回采充填過程中，右礦柱在-380 m水平出現(xiàn)開裂破壞，左右礦柱均在礦柱兩側(cè)礦房被回采之后發(fā)生破壞，回采充填完成后礦柱破壞情況如圖5所示。

3 試驗結(jié)果分析

3.1 變形分析

通過各中段回采充填完成后的監(jiān)測位移場與原始場進行對比分析，可以得到-450 m中段和-390 m中段回采充填完成后相對于試驗前的變形情況，如圖6所示。分析得到：中間盤區(qū)礦房頂板沉降量大于兩側(cè)盤區(qū)礦房頂板沉降量。-450 m中段回采充填完成后，2號礦房頂板最大沉降量約0.332 mm，1、3號礦房頂板最大沉降量約0.291 mm，圍巖礦柱變形連續(xù)，最大變形出現(xiàn)在中間礦房頂板處，采場穩(wěn)定性良好。-390 m中段回采充填完成后，5號礦房頂板最大沉降量為1.491 mm，4、6號礦房頂板最大沉降量約為1.193 mm，礦柱中部出現(xiàn)變形不連續(xù)現(xiàn)象，采場基本保持穩(wěn)定。

3.2 應(yīng)力分析

使用應(yīng)力采集系統(tǒng)記錄了4個監(jiān)測點的應(yīng)力隨回采充填步驟變化的數(shù)據(jù)，得到如圖7所示的監(jiān)測點應(yīng)力隨時間的變化曲線。分析得到：1號監(jiān)測點的應(yīng)力在前8步回采充填過程中基本保持穩(wěn)定，在第9步回采充填過程中出現(xiàn)應(yīng)力跌落，導(dǎo)致此現(xiàn)象的原因是1號監(jiān)測點位于5號礦房的頂板，第9步回采充填發(fā)生在5號礦房內(nèi)，且位于監(jiān)測點正下方；2號監(jiān)測點的應(yīng)力在前10步回采充填過程中基本保持穩(wěn)定，在第11步回采充填過程中出現(xiàn)應(yīng)力跌落，導(dǎo)致此現(xiàn)象的原因是2號監(jiān)測點位于6號礦房的頂板，第11步回采充填發(fā)生在6號礦房內(nèi)，且位于監(jiān)測點正下方。1、2號頂板監(jiān)測點從試驗開始至結(jié)束應(yīng)力處于先穩(wěn)定再下降后保持穩(wěn)定的趨勢，在-450 m中段回采充填時1、2號監(jiān)測點處于壓力拱內(nèi)，隨著-390 m中段回采充填工作的進行應(yīng)力拱位置上移，1、2號監(jiān)測點位于應(yīng)力拱下方，開挖卸荷作用和應(yīng)力重分布導(dǎo)致了1、2號頂板監(jiān)測點的應(yīng)力出現(xiàn)下降。4號礦房內(nèi)的第7步回采充填使得左礦柱被揭露出來，導(dǎo)致了3號監(jiān)測點的應(yīng)力突躍；5號礦房內(nèi)的第9步回采充填使得右側(cè)礦柱被揭露出來，導(dǎo)致了4號監(jiān)測點的應(yīng)力突躍。3、4號礦柱監(jiān)測點從試驗開始至結(jié)束，礦柱中應(yīng)力處于持續(xù)增長趨勢，表明原巖礦柱作為采場主要支承結(jié)構(gòu)并未失穩(wěn)，一直發(fā)揮著支承作用。

將沉降量乘以幾何相似比CL，應(yīng)力監(jiān)測值乘以應(yīng)力相似比Cσ可得到實際值，表2為各礦房頂板在-450 m中段和-390 m中段回采充填完成后的實際變形量，表3為各監(jiān)測點在-450 m中段和-390 m中段回采充填完成后的實際應(yīng)力值。在新中段開采前后，中間礦房頂板最大沉降量從66 mm增長為298 mm，在原來的基礎(chǔ)上增長了3.5倍；兩側(cè)礦房頂板最大沉降量從49 mm增長為236 mm，在原來的基礎(chǔ)上增長了3.8倍。3、4號礦柱應(yīng)力從19.4 MPa增長到32.8 MPa，應(yīng)力水平在原來的基礎(chǔ)上提升了69%。新中段的回采導(dǎo)致了礦房頂板沉降量和礦柱內(nèi)應(yīng)力集中程度顯著增大，有必要采取相關(guān)措施和技術(shù)減少采場環(huán)境惡化帶來的危害。

4 數(shù)值模擬

4.1 模擬方案

通過相似材料試驗中的應(yīng)力和變形監(jiān)測數(shù)據(jù)發(fā)現(xiàn)，在-390 m新中段回采充填過程中，礦柱內(nèi)應(yīng)力集中程度明顯提升，并發(fā)生了破壞現(xiàn)象，礦房頂板是采場內(nèi)變形量最大的部位，但并沒有發(fā)生明顯破壞現(xiàn)象。因此實際回采中礦柱是最有可能發(fā)生破壞的部位，為提高礦柱承載能力，增加礦柱寬度是普遍做法。本研究通過數(shù)值模擬方法研究不同礦柱寬度（表4）和不同回采順序?qū)Σ蓤龇€(wěn)定性的影響。

計算模型如圖8所示，模型中-450 m中段和-390 m中段各5個盤區(qū)礦房，計算采用了與現(xiàn)場一致的順序回采方案，即按照1～10號礦房依次回采充填，同時采用了間隔回采方案，即按照1、3、5、2、4、6、8、10、7、9號礦房的順序回采充填。

4.2 計算分析

對表4中的3種方案分別進行建模計算，以塑性區(qū)、變形和應(yīng)力為指標(biāo)，分析對比不同方案之間各指標(biāo)的變化趨勢，圖9為不同回采方案與不同回采充填順序的計算結(jié)果。

分析圖9可知：不同回采方案均在礦柱內(nèi)出現(xiàn)應(yīng)力集中現(xiàn)象，最大值出現(xiàn)在礦柱底部，隨著礦柱寬度增加礦柱內(nèi)最大主應(yīng)力逐漸減小，順序回采充填與間隔回采充填的最大主應(yīng)力基本一致；不同回采方案的最大拉應(yīng)力均出現(xiàn)在兩側(cè)礦房的頂角處，間隔回采充填比順序回采充填形成的最大拉應(yīng)力?。徊煌夭煞桨缸畲蟪两底冃纬霈F(xiàn)在8號礦房頂板中，兩側(cè)礦房頂板逐漸減小，呈現(xiàn)出對稱趨勢，間隔回采充填比順序回采充填的頂板沉降量??；不同回采充填方案的塑性區(qū)主要集中在礦柱中下部和礦房頂板中，隨著礦柱寬度的增加礦柱中的塑性區(qū)范圍逐漸減小，順序回采比間隔回采在礦柱中形成的塑性區(qū)更小。由于礦房內(nèi)礦體的回采充填，采場內(nèi)應(yīng)力平衡狀態(tài)被打破并重新分布，位移和塑性區(qū)隨著回采充填而不斷積累，最大沉降出現(xiàn)在中間礦房頂板，導(dǎo)致中間2、3號礦柱應(yīng)力集中程度和塑性區(qū)范圍大于兩側(cè)1、4號礦柱，礦房內(nèi)礦體的回采充填使得采區(qū)上方形成應(yīng)力拱，導(dǎo)致在兩側(cè)6、10號礦房邊墻頂部形成了剪切破壞區(qū)域。

各方案及不同順序回采情況下采場內(nèi)各項指標(biāo)數(shù)據(jù)統(tǒng)計結(jié)果如表5所示。分析可知：隨著礦柱寬度增加，采場內(nèi)頂板沉降量、最大主應(yīng)力、最小主應(yīng)力、塑性區(qū)范圍都呈現(xiàn)出下降趨勢。各方案采場沉降量順序回采比間隔回采大，但相差較小均不超過10 mm，采場最大拉應(yīng)力值順序回采大于間隔回采，順序回采與間隔回采在礦柱內(nèi)的應(yīng)力集中程度基本行相同。

5 結(jié) 論

通過相似材料試驗和數(shù)值模擬對采場穩(wěn)定性及采場應(yīng)力變形規(guī)律進行研究，得出了以下結(jié)論：

（1）采用18 m寬礦柱進行-390 m中段回采充填時，礦柱會在-380～-370 m水平發(fā)生剪切破壞，但礦柱并未完全失穩(wěn)且能夠發(fā)揮主要支承作用，采場基本穩(wěn)定。

（2）礦柱是階段空場嗣后充填采場中主要支撐結(jié)構(gòu)，在礦體回采過程中當(dāng)?shù)V柱被揭露時，礦柱內(nèi)應(yīng)力會明顯增大，隨著礦柱應(yīng)力逐步積累會出現(xiàn)剪切破壞現(xiàn)象。

（3）巖石質(zhì)量中等以上的急傾斜厚大礦體在進行階段空場嗣后充填開采時，采場礦柱內(nèi)會產(chǎn)生應(yīng)力集中，最大應(yīng)力出現(xiàn)在礦柱底部；采場最大變形出現(xiàn)在中間盤區(qū)礦房頂板，兩側(cè)盤區(qū)礦房頂板變形呈對稱性減??；盤區(qū)礦房間隔回采充填引起采場頂板變形比順序回采充填小，有利于控制采場變形，提高采場穩(wěn)定性。