曹建立， 譚寶會(huì)， 劉 洋， 任鳳玉

(東北大學(xué) 資源與土木工程學(xué)院， 遼寧 沈陽(yáng) 110819)

由于歷史原因，在我國(guó)礦業(yè)發(fā)展過(guò)程中不少地下礦山在淺層形成了大量的民采空區(qū)，這些民采空區(qū)數(shù)量大、規(guī)模不一、分布零散，破壞了礦體的完整性和穩(wěn)定性，且隨時(shí)可能冒落，嚴(yán)重影響了礦山生產(chǎn).因此開展多空區(qū)礦體冒落機(jī)理及回采工藝研究對(duì)于提升多空區(qū)礦體采礦效率及安全性有重要意義.在機(jī)理研究方面，Wang等在考慮巖體流變特性的基礎(chǔ)上，利用新建立的黏彈性模型，對(duì)采空區(qū)頂板巖層的蠕變破壞進(jìn)行了分析，并利用該分析方法，對(duì)邢臺(tái)石膏礦頂板巖層的破壞時(shí)間進(jìn)行了估算，所得結(jié)果與實(shí)測(cè)結(jié)果基本一致[1]；Sun等對(duì)弓長(zhǎng)嶺三礦區(qū)多采空區(qū)穩(wěn)定性進(jìn)行了分析評(píng)價(jià)，并研究了采礦過(guò)程中采空區(qū)與非采空區(qū)地層之間的相互作用關(guān)系[2]；鄭磊等采用FLAC3D軟件對(duì)金屬礦多空區(qū)的崩落放頂過(guò)程進(jìn)行了數(shù)值模擬分析，模擬結(jié)果反映了采空區(qū)范圍的地壓活動(dòng)規(guī)律及礦柱穩(wěn)定性變化情況[3].在采礦工藝方面，趙興明采用RAMAC/GPR探地雷達(dá)和瑞雷波探測(cè)技術(shù)，對(duì)復(fù)雜采空區(qū)進(jìn)行了綜合探測(cè)研究，并以此為基礎(chǔ)提出多空區(qū)下分段鑿巖連續(xù)崩礦回采方案[4]；李興尚等針對(duì)民采多空區(qū)礦體，提出了采用全尾砂上向分層充填開采技術(shù)回采殘礦[5].

限于民采空區(qū)的復(fù)雜性，常規(guī)的采礦方法很難真正實(shí)現(xiàn)多空區(qū)礦體的安全高效回采.誘導(dǎo)冒落法是近年來(lái)發(fā)展起來(lái)的一項(xiàng)新型改進(jìn)的采礦方法，在破碎礦體高效回采及采礦區(qū)安全高效處理方面取得了顯著成效.本文針對(duì)大嶺西山菱鎂礦存在著大量民采空區(qū)，提出了利用誘導(dǎo)冒落法回采多空區(qū)礦體的技術(shù)方案，采用3DEC軟件研究了多空區(qū)礦體的誘導(dǎo)冒落機(jī)理，并模擬了4種不同拉底順序的誘導(dǎo)冒落方案，確定出了最優(yōu)拉底方案.

1 多空區(qū)礦體誘導(dǎo)冒落采礦方案

1.1 工程背景

大嶺西山菱鎂礦遭受多年無(wú)序開采，在+137 m水平以上分布有大量采空區(qū)，采空區(qū)層疊交錯(cuò)(見圖1).傳統(tǒng)采礦方法的回采安全問題難以解決，需要從礦山多空區(qū)礦體的特點(diǎn)、礦巖結(jié)構(gòu)及穩(wěn)定性特征出發(fā)，研究出適合該礦山的安全高效采礦方案.

1.2 多空區(qū)礦體誘導(dǎo)冒落回采方案

誘導(dǎo)冒落采礦法是由東北大學(xué)任鳳玉提出的一種新型改進(jìn)的地下采礦方法，該采礦方法具有產(chǎn)能大、成本低、結(jié)構(gòu)簡(jiǎn)單、應(yīng)用靈活等優(yōu)點(diǎn)[6-7].目前該方法已經(jīng)成功在國(guó)內(nèi)多個(gè)礦山的礦石回采及采空區(qū)處理中得到應(yīng)用.

經(jīng)過(guò)現(xiàn)場(chǎng)調(diào)查及取樣測(cè)試后得出大嶺西山菱鎂礦穩(wěn)定性級(jí)別為中等穩(wěn)定，圍巖穩(wěn)定性級(jí)別為中等穩(wěn)定至不穩(wěn)定，比較適合誘導(dǎo)冒落方案的開展.采用誘導(dǎo)冒落法開采多空區(qū)礦體時(shí)，在首采分段通過(guò)有計(jì)劃的中深孔爆破形成一個(gè)貫通的大規(guī)模采空區(qū)，誘導(dǎo)頂板發(fā)生冒落(見圖2)；在首采分段下部通常至少布置1～2個(gè)采礦分層，以確保冒落礦石能夠被充分回收；位于最底部的分段通常是采場(chǎng)內(nèi)礦石最后的回收機(jī)會(huì)，因此為減少礦石損失可以通過(guò)加密進(jìn)路的方式來(lái)提升礦石的回采率.利用誘導(dǎo)冒落法回采多空區(qū)礦體，不僅可以大幅度提升采礦率，而且可有效避免空區(qū)群組失穩(wěn)對(duì)采礦作業(yè)帶來(lái)的危害，提高了采礦作業(yè)安全性.

2 多空區(qū)礦體誘導(dǎo)冒落機(jī)理

3DEC作為一款基于離散單元法的基本理論來(lái)描述離散介質(zhì)力學(xué)行為的計(jì)算分析程序，在模擬巖體的開裂失穩(wěn)以及冒落等方面具有獨(dú)到的優(yōu)勢(shì)，目前已經(jīng)被廣泛應(yīng)用于研究邊坡失穩(wěn)、頂板冒落、巷道失穩(wěn)等工程問題[8-10].基于所研究問題的特性，本文采用3DEC軟件對(duì)多空區(qū)下誘導(dǎo)冒落機(jī)理展開研究.

2.1 模型的建立

本次所建模型尺寸為長(zhǎng)×高×寬=2 400 mm×1 100 mm×200 mm，模型中共包含了17個(gè)分布在不同高度的采空區(qū)，為便于建模，空區(qū)形態(tài)均簡(jiǎn)化為矩形，空區(qū)高度在5～10 mm之間，長(zhǎng)度在10～30 mm之間.模型從下至上共劃分了9個(gè)巖層，模型底部第1巖層設(shè)定為拉底水平.模型具有6個(gè)自由面，對(duì)應(yīng)的6個(gè)邊界條件為：模型左右邊界約束在X=0和X=2 400 mm,限制X方向的速度；前后邊界約束在Y=0和Y=200 mm,限制Y方向速度；底部邊界約束在Z=0為固定邊界，頂部邊界在Z=1 100 mm施加均布荷載，其值為0.145 MPa，應(yīng)力相似常數(shù)aσ=10，模擬埋深為50 m，模型達(dá)到應(yīng)力平衡后進(jìn)行開挖.由于空區(qū)埋深較淺，本次數(shù)值模擬只考慮上覆巖層的垂直應(yīng)力，對(duì)水平構(gòu)造應(yīng)力不考慮.同時(shí)為了便于后續(xù)分析，將模型中的空區(qū)按照上下重疊交錯(cuò)的位置關(guān)系，以及空區(qū)群與模型中軸線的空間位置關(guān)系，共劃分出了4個(gè)空區(qū)群，如圖3所示.

2.2 力學(xué)模型和參數(shù)的選取

數(shù)值模型中巖體的力學(xué)模型采用莫爾-庫(kù)倫模型，節(jié)理本構(gòu)模型采用庫(kù)侖-滑移模型，數(shù)值模型包含了節(jié)理彈性剛度、摩擦特性、剪脹特性等功能，符合巖體的基本特征.數(shù)值模型中巖體及節(jié)理的力學(xué)參數(shù)如表1和表2所示.

表1 巖體力學(xué)參數(shù)

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

2.3 多空區(qū)礦體冒落機(jī)理分析

采場(chǎng)應(yīng)力平衡后，在最底部分層由中央向兩側(cè)進(jìn)行退采拉底.礦塊在拉底破壞擾動(dòng)作用下應(yīng)力重新調(diào)整，由于礦體受水平層面及空區(qū)間礦柱強(qiáng)度的限制，調(diào)整后的應(yīng)力使得礦體沿層面破壞或者超過(guò)空區(qū)間礦柱的強(qiáng)度而發(fā)生空區(qū)失穩(wěn).礦體破壞后又形成新的應(yīng)力平衡，在不斷的應(yīng)力調(diào)整、礦體破壞、空區(qū)失穩(wěn)的循環(huán)中，最終達(dá)到應(yīng)力平衡，使整個(gè)采場(chǎng)中的巖體趨于穩(wěn)定.圖4給出了模型拉底后空區(qū)間礦柱發(fā)生裂紋擴(kuò)展的現(xiàn)象以及空區(qū)失穩(wěn)后的礦體冒落現(xiàn)象.

數(shù)值模擬研究證明空區(qū)群的存在嚴(yán)重影響了礦體的完整性和冒落規(guī)律.空區(qū)群弱面的存在使得礦體往往沿空區(qū)弱面發(fā)生冒落.為了更好地解釋多空區(qū)礦體冒落過(guò)程的應(yīng)力演化及礦柱破壞的先后順序，對(duì)礦體整個(gè)冒落過(guò)程產(chǎn)生的裂紋進(jìn)行了編號(hào)，見圖5，圖中紅色線條代表裂紋，數(shù)字代表裂紋產(chǎn)生的先后序號(hào).

在采場(chǎng)拉底形成后，裂紋發(fā)展及冒落過(guò)程如下：① 在右1空區(qū)群弱面位置，裂紋1和裂紋2開始萌生、擴(kuò)展；② 在裂紋1和裂紋2擴(kuò)展過(guò)程中，拉底擾動(dòng)產(chǎn)生的應(yīng)力場(chǎng)在不斷地調(diào)整，重疊交錯(cuò)的左1空區(qū)群和右2空區(qū)群中的裂紋3、裂紋4、裂紋5、裂紋6開始萌生、擴(kuò)展；③ 隨后左1空區(qū)群中的裂紋7和裂紋8在應(yīng)力調(diào)整過(guò)程中也開始萌生；④ 待裂紋1～8貫通后，底部開始形成一個(gè)整體，最后在底部第5層巖層層面張開，從而產(chǎn)生了底部第2層到第5層礦體冒落；⑤ 底部第2層到第5層礦體冒落后，底部第6層的礦體受到瞬間卸荷作用，礦體受到強(qiáng)烈的擾動(dòng)作用，在應(yīng)力調(diào)整的過(guò)程中，以及第7層巖體層面張開下，最后第6層和第7層巖體開始冒落；⑥ 在第6層巖體垮落后，第7層巖體和第8層巖體受到強(qiáng)烈的瞬間卸荷作用，在垮落后的兩幫位置，巖體產(chǎn)生應(yīng)力集中，內(nèi)部開始萌生出新生裂紋9和裂紋10，然后擴(kuò)展貫通，該巖層層面在巖層巖體的重力作用下張開，最終垮落，見圖4b.

數(shù)值模擬結(jié)果表明，空區(qū)群的存在打破了無(wú)空區(qū)礦體的應(yīng)力均勻分布的狀態(tài)，拉底后空區(qū)間的礦柱上產(chǎn)生了較大的應(yīng)力集中，礦柱尺寸越小、周圍空區(qū)越多，則應(yīng)力集中程度越大，在集中應(yīng)力作用下礦柱很容易產(chǎn)生裂隙而破壞.根據(jù)產(chǎn)生裂隙的兩空區(qū)相對(duì)位置不同，裂隙產(chǎn)生機(jī)理分為剪切破壞和剪切-拉伸破壞，如圖6所示，當(dāng)上下層空區(qū)邊界位于(或接近)同一豎直平面內(nèi)時(shí)，礦柱受到剪切破壞產(chǎn)生裂隙；當(dāng)上下層空區(qū)部分重疊且邊界拉開一定距離時(shí)，礦柱受到拉伸-剪切破壞產(chǎn)生裂隙.當(dāng)?shù)V柱間產(chǎn)生的裂隙與礦體層面或破壞面貫通后，礦體便發(fā)生垮塌冒落.

3 多空區(qū)礦體誘導(dǎo)冒落拉底方案優(yōu)選

3.1 拉底方案設(shè)計(jì)

由于在拉底過(guò)程中，各空區(qū)群所受的應(yīng)力狀態(tài)各不一樣，空區(qū)間礦柱裂紋貫通的時(shí)間也不一樣，因此不同的拉底順序?qū)?huì)導(dǎo)致不同的冒落時(shí)間和形態(tài).此處設(shè)定了4種不同的拉底方案進(jìn)行數(shù)值模擬計(jì)算，最后，通過(guò)比較冒落的先后順序、冒落過(guò)程的平穩(wěn)度來(lái)確定最優(yōu)拉底方案，原則上最先發(fā)生冒落且冒落過(guò)程平緩持續(xù)進(jìn)行者為最優(yōu)方案.4種拉底方案如圖7所示.其中，方案1以左1空區(qū)群為中心向左拉底；方案2以左1空區(qū)群為中心向兩側(cè)拉底；方案3以右2空區(qū)群為中心向左拉底；方案4以右2空區(qū)群為中心向兩側(cè)拉底.4個(gè)方案中均以每次擴(kuò)展0.2 m進(jìn)行拉底.

3.2 模擬結(jié)果及方案優(yōu)選

圖8給出了四種拉底方案分別在拉底長(zhǎng)度為0.4 m和0.8 m時(shí)的數(shù)值模擬計(jì)算結(jié)果.

由圖8a可知，當(dāng)拉底長(zhǎng)度為0.4 m時(shí)，第2層開始冒落；當(dāng)拉底長(zhǎng)度達(dá)到0.8 m時(shí)，冒落高度達(dá)到第8層礦體.冒落范圍涉及到左1空區(qū)群和左2空區(qū)群，且隨著冒落向上發(fā)展，冒落范圍有向右側(cè)右1空區(qū)群發(fā)展的趨勢(shì).主要原因是因?yàn)樽?空區(qū)群左側(cè)為完整礦體，未受空區(qū)破壞，而右側(cè)為多個(gè)空區(qū)群，礦體強(qiáng)度降低較多，拉底造成頂板能量聚集導(dǎo)致的礦體破壞通常會(huì)朝強(qiáng)度較低一側(cè)發(fā)展.

由圖8b可知，當(dāng)拉底長(zhǎng)度為0.4 m時(shí)，第3層巖層冒落；當(dāng)拉底長(zhǎng)度達(dá)到0.8 m時(shí)，底部第2層礦體開始發(fā)生冒落，同時(shí)帶動(dòng)其上大范圍冒落，冒落高度達(dá)到了第8層，拉底空間左側(cè)的冒落范圍基本與拉底范圍在同一豎直平面內(nèi)，而在右側(cè)則超出拉底范圍豎直面較多.主要原因是以左1空區(qū)群為中心同時(shí)向左右兩側(cè)拉底時(shí)，拉底范圍同時(shí)包含了左2空區(qū)群和右1空區(qū)群，且右1空區(qū)群向右側(cè)延伸較大，由此造成冒落范圍超出拉底范圍較大.

由圖8c計(jì)算結(jié)果可知，當(dāng)拉底長(zhǎng)度為0.4 m時(shí)巖層未發(fā)生冒落；當(dāng)拉底長(zhǎng)度為0.8 m時(shí)，突然發(fā)生較大規(guī)模冒落.主要原因是由于初始拉底位置上方是右2空區(qū)群的中心，上部載荷被轉(zhuǎn)移到空區(qū)兩側(cè)的礦柱上，拉底空間上方屬于卸壓區(qū)，此時(shí)施加在拉底空間上方完整巖層上的應(yīng)力較低，無(wú)法使巖層破壞；當(dāng)拉底長(zhǎng)度超過(guò)0.8 m時(shí)，拉底范圍超過(guò)了空區(qū)兩側(cè)礦柱，此時(shí)底部第2層巖層受到上部礦柱集中應(yīng)力的作用而發(fā)生破壞，一旦兩個(gè)礦柱失穩(wěn)，則其上部和中間的空區(qū)均失穩(wěn)破壞，發(fā)生較大規(guī)模的瞬間冒落，形成冒落危害.

由圖8d可知，當(dāng)拉底長(zhǎng)度為0.4 m時(shí)，未發(fā)生冒落，此時(shí)拉底空間上部存在一個(gè)梯形礦柱，該礦柱下部尺寸較大，應(yīng)力集中程度低，因此即使礦柱完全位于拉底空間上方時(shí)仍未發(fā)生冒落；當(dāng)拉底長(zhǎng)度達(dá)到0.8 m時(shí)，雖然拉底空間上部巖體發(fā)生大面積破壞，但依舊沒有發(fā)生冒落，可以推測(cè)，若繼續(xù)增大拉底面積，則極有可能突然發(fā)生大規(guī)模冒落.

對(duì)比4個(gè)方案可知，方案1和方案2冒落較早，方案3次之，方案4最難冒落.雖然方案1和方案2均能實(shí)現(xiàn)早冒落、持續(xù)冒落的目標(biāo)，但方案2所需的時(shí)間更短，考慮到拉底作業(yè)應(yīng)盡快撤出空區(qū)下方，盡可能減小冒落沖擊傷害的發(fā)生，因此選擇方案2作為最優(yōu)方案.數(shù)值模擬實(shí)驗(yàn)表明，對(duì)于多空區(qū)礦體而言，應(yīng)盡可能地將拉底工程布置在最易破壞的礦柱下方，并以其為中心向四周不斷擴(kuò)大拉底面積，從而實(shí)現(xiàn)上部多空區(qū)礦體的持續(xù)、平穩(wěn)冒落.

4 大嶺西山菱鎂礦誘導(dǎo)冒落拉底方案

根據(jù)數(shù)值模擬研究結(jié)果，將大嶺西山菱鎂礦拉底初始位置設(shè)定在+120 m水平4號(hào)進(jìn)路與2號(hào)聯(lián)巷的交叉點(diǎn)，如圖9所示.在此處拉切割槽，

向四周退采，形成階梯狀退采形態(tài)，彼此保持 10 m 左右的超前距離，當(dāng)各條進(jìn)路退采接近1號(hào)和3號(hào)聯(lián)巷時(shí)，與聯(lián)巷保持約5 m長(zhǎng)度即停止回采，以便于最邊界處的16號(hào)進(jìn)路運(yùn)輸暢通.此時(shí)，拉底空間足以使首采區(qū)域上部采空區(qū)冒落，由于多條進(jìn)路與冒落空區(qū)聯(lián)通，即使有冒落氣浪產(chǎn)生，也會(huì)被多條進(jìn)路分流，從而降低了氣浪沖擊的可能性.

5 結(jié) 論

1) 空區(qū)群的存在破壞了礦體的完整性和穩(wěn)定性，難以利用常規(guī)采礦方法實(shí)現(xiàn)多空區(qū)礦體的安全高效開采，利用誘導(dǎo)冒落法回采多空區(qū)礦體時(shí)，合理布置拉底位置和順序可利用地壓快捷地使多空區(qū)礦體發(fā)生自然冒落，再通過(guò)底部結(jié)構(gòu)回收礦石，從而實(shí)現(xiàn)多空區(qū)礦體的安全高效回采.

2)數(shù)值模擬實(shí)驗(yàn)結(jié)果表明，對(duì)多空區(qū)層狀礦體實(shí)施誘導(dǎo)冒落時(shí)，拉底后空區(qū)之間的礦柱上容易產(chǎn)生較大的應(yīng)力集中，而且礦柱尺寸越小、周圍空區(qū)越多，則應(yīng)力集中程度越大，在集中應(yīng)力作用下礦柱很容易產(chǎn)生破壞裂隙，隨著裂隙的產(chǎn)生，應(yīng)力也在不斷調(diào)整，直至裂隙與層狀礦體的層理面貫通時(shí)，礦體便會(huì)發(fā)生冒落.

3)通過(guò)對(duì)4種不同退采順序的拉底方案進(jìn)行數(shù)值模擬研究，結(jié)果表明對(duì)于多空區(qū)礦體實(shí)施誘導(dǎo)冒落時(shí)，應(yīng)盡可能地將拉底工程布置在最易破壞的礦柱下方，并以其為中心向四周不斷擴(kuò)大拉底面積，更有利于實(shí)現(xiàn)上部多空區(qū)礦體的持續(xù)、平穩(wěn)冒落.