何榮興，張 晶，宋德林，劉 歡

(1.東北大學(xué)資源與土木工程學(xué)院，遼寧 沈陽(yáng) 110819；2.內(nèi)蒙古科技大學(xué)礦業(yè)與煤炭學(xué)院，內(nèi)蒙古 包頭 014010)

地表巖移是地下采礦引起的普遍現(xiàn)象，對(duì)礦山生產(chǎn)安全影響較大，尤其是傾斜厚礦體的開(kāi)采所引起的地表巖移，具有變形量大、陷落速度快和巖移陷落范圍大等特點(diǎn)，嚴(yán)重制約礦山安全生產(chǎn)，需要采取有效措施予以控制[1]。

金屬礦山地質(zhì)條件復(fù)雜多變，地下開(kāi)采引起的巖層移動(dòng)和地表變形受采礦方法、埋深、采高等諸多因素影響，規(guī)律性較差。很多學(xué)者對(duì)金屬礦山的地表巖移問(wèn)題進(jìn)行了大量的研究，如朱家橋[2]總結(jié)了程潮鐵礦東區(qū)采礦不同時(shí)期巖體冒落及塌陷坑的發(fā)展過(guò)程，認(rèn)為程潮鐵礦地表變形破壞形式以傾倒式拉裂破壞為主，巖體水平位移大于沉降位移；席人雙等[3]對(duì)程潮鐵礦東區(qū)進(jìn)行地表巖移的GPS監(jiān)測(cè)，積累近三年的變形數(shù)據(jù)，從而準(zhǔn)確地反映了的地表變形以向采空區(qū)的階梯變形為主等變形規(guī)律；蔡美峰等[4]建立金屬礦山開(kāi)采沉陷預(yù)計(jì)系統(tǒng)MSDAS-GIS，用于模擬預(yù)測(cè)北洺河鐵礦井下連續(xù)回采的開(kāi)采沉陷，對(duì)-50 m水平開(kāi)采后地表沉降進(jìn)行現(xiàn)場(chǎng)監(jiān)測(cè)，預(yù)測(cè)結(jié)果和監(jiān)測(cè)結(jié)果具有較好的相似性；趙淑霞[5]以地理信息系統(tǒng)(GIS)技術(shù)為基礎(chǔ)，開(kāi)發(fā)建立礦區(qū)地表動(dòng)態(tài)三維地質(zhì)立體模型，實(shí)現(xiàn)對(duì)段村-雷溝礦區(qū)地面地質(zhì)、地形、地面沉降監(jiān)測(cè)，研究了上覆巖層變形移動(dòng)規(guī)律。這些研究都是建立在地表巖移數(shù)據(jù)監(jiān)測(cè)和深部巖體地壓活動(dòng)觀測(cè)的基礎(chǔ)上，耗時(shí)較長(zhǎng)。

近年來(lái)，數(shù)值模擬軟件的更新?lián)Q代使得模擬結(jié)果更加貼近工程實(shí)際，且建模多樣性及獨(dú)特的算法使得數(shù)值模擬更適用于工程現(xiàn)象的規(guī)律研究。眾多學(xué)者[6-10]選擇數(shù)值模擬來(lái)研究地表巖移的規(guī)律和范圍，并取得很好的科研成果。本文根據(jù)弓長(zhǎng)嶺鐵礦東南區(qū)上含鐵帶、下含鐵帶的開(kāi)采現(xiàn)狀，利用FLAC3D軟件模擬對(duì)比分析了不考慮崩落散體作用和考慮散體作用兩種工況條件下的采動(dòng)巖移范圍和充填廢石散體對(duì)巖移的控制作用，提出了剝離上含鐵帶露天采場(chǎng)的廢石來(lái)充填下含鐵帶地采塌陷坑的巖移控制措施，并確定了下含鐵帶的極限開(kāi)采深度。

1 工程概況

弓長(zhǎng)嶺鐵礦東南區(qū)礦體由多條平行礦帶組成，如圖1所示，礦體傾角為65°～85°，平均厚度為15 m，屬于急傾斜中厚礦體。上含鐵帶(圖1)為赤鐵礦，應(yīng)用露天采礦的方法進(jìn)行開(kāi)采，露天開(kāi)采臺(tái)階為+448～+388 m。下含鐵帶(圖1)為磁鐵礦，采用無(wú)底柱崩落法進(jìn)行開(kāi)采，開(kāi)采水平為192 m，已經(jīng)冒透地表，形成地表塌陷坑。因開(kāi)采規(guī)劃不合理，上含鐵帶、下含鐵帶的高差約300 m，最小水平距離約150 m，下含鐵帶的超深開(kāi)采使得近露天采場(chǎng)一側(cè)的圍巖發(fā)生巖移及塌陷，嚴(yán)重威脅露天采場(chǎng)的安全生產(chǎn)。該礦區(qū)采礦生產(chǎn)設(shè)計(jì)中，地采極限開(kāi)采深度為152 m，極限開(kāi)采深度以下礦體待露天開(kāi)采閉坑后再行開(kāi)采。此生產(chǎn)設(shè)計(jì)將導(dǎo)致該礦區(qū)下含鐵帶152 m水平開(kāi)采結(jié)束后產(chǎn)能接替失衡，影響礦山生產(chǎn)效益。

圖1 弓長(zhǎng)嶺鐵礦東南區(qū)開(kāi)采現(xiàn)狀圖Fig.1 Current mining situation of southeast area of Gongchangling Iron Mine

為保障露天采場(chǎng)和地下采場(chǎng)持續(xù)生產(chǎn)以滿足礦區(qū)產(chǎn)能需求，提出剝離上含鐵帶露天采場(chǎng)的廢石來(lái)充填下含鐵帶地采塌陷坑的技術(shù)方案，以期控制塌陷坑上盤圍巖移動(dòng)。因此，有必要探討充填廢石散體作用下的采動(dòng)巖移范圍和不影響上含鐵帶露天開(kāi)采的地采極限采深。

2 計(jì)算模型和計(jì)算方案

2.1 數(shù)值模型建立

目前弓長(zhǎng)嶺鐵礦東南區(qū)下含鐵帶地采塌陷坑體積較小，對(duì)數(shù)值模擬結(jié)果影響可以忽略，為此直接將礦體上部邊界延伸至地表。利用ANSYS軟件建立長(zhǎng)×寬×高為400 m×100 m×350 m的簡(jiǎn)化數(shù)值模型，其中礦體水平厚度為15 m，礦體傾角為68°，礦體埋深方向取350～0 m，沿礦體走向長(zhǎng)度為100 m，沿礦體傾向取400 m，劃分網(wǎng)格導(dǎo)入FLAC3D軟件數(shù)值模型，如圖2所示，計(jì)算域邊界采取位移約束。

圖2 下含鐵帶巖移模擬的FLAC3D模型Fig.2 FLAC3D model of rock movement on the iron-bearing belt mining

模型分五步開(kāi)挖用以模擬下含鐵帶回采，第一步開(kāi)挖至230 m水平模擬目前的塌陷坑，之后每次開(kāi)挖高度取40 m，至70 m水平時(shí)停止開(kāi)挖。位移監(jiān)測(cè)點(diǎn)布置于地表，按網(wǎng)格劃分的節(jié)點(diǎn)布置，各監(jiān)測(cè)點(diǎn)距z軸的距離依次為12.75 m、26.02 m、41.48 m、50.93 m、67.38 m、75.02 m、85.57 m、90.08 m、97.60 m、105.11 m、112.63 m、120.15 m、127.63 m、135.15 m、142.66 m、150.18 m、157.68 m、165.21 m、172.76 m、180.39 m、188.24 m、196.78 m、207.55 m、219.12 m、230.32 m、245.70 m、255.00 m。每步開(kāi)挖后均進(jìn)行數(shù)值求解并記錄各監(jiān)測(cè)點(diǎn)的水平位移和沉降位移，用以分析不同采深條件下的巖移規(guī)律及范圍。模擬下含鐵帶回采至70 m水平后逐步模擬充填，每次充填高度為40 m，直至地表。每次充填后進(jìn)行數(shù)值求解并記錄各監(jiān)測(cè)點(diǎn)的水平位移和沉降位移，用以分析充填散體作用下采動(dòng)巖移的規(guī)律和范圍。

2.2 巖體力學(xué)參數(shù)估算

對(duì)弓長(zhǎng)嶺鐵礦東南區(qū)進(jìn)行地質(zhì)調(diào)查和點(diǎn)荷載試驗(yàn)，獲得的巖體基本質(zhì)量指標(biāo)見(jiàn)表1。為獲取數(shù)值模擬所需的巖體力學(xué)參數(shù)，本文基于表1中的參數(shù)，結(jié)合Hoek-Brown強(qiáng)度準(zhǔn)則[11-12]和Mohr-Coulomb準(zhǔn)則[13]來(lái)估算數(shù)值模擬所需的巖體力學(xué)參數(shù)。

Hoek-Brown巖體破壞準(zhǔn)則表達(dá)式見(jiàn)式(1)。

σ1=σ3+σc(mσ3/σc+s)0.5

(1)

式中：σ1、σ3為巖體破壞時(shí)的最大主應(yīng)力、最小主應(yīng)力；m為Hoek-Brown常數(shù)，估算見(jiàn)式(2)；s為無(wú)量綱系數(shù)，估算見(jiàn)式(3)[14]。

(2)

(3)

式中：mi可根據(jù)巖體的Hoek-Brown常數(shù)表來(lái)選取[11-12]；RMR根據(jù)巖石的強(qiáng)度、RQD值、節(jié)理間距以及地下水影響的綜合評(píng)價(jià)巖體分類指標(biāo)，對(duì)應(yīng)巖石力學(xué)分類表(RMR)[15-16]獲得每項(xiàng)指標(biāo)評(píng)分求和得出。

當(dāng)σ1=0時(shí)，由式(1)可知，巖體的單軸抗拉強(qiáng)度σmt可表示為式(4)。

(4)

巖體彈性模量Em與地質(zhì)強(qiáng)度指標(biāo)GSI的關(guān)系可表示為式(5)[12]。

(5)

地質(zhì)強(qiáng)度指標(biāo)GSI取值可以用表1中的RMR值代替[14]，則巖體彈性模量Em可求。

如圖2所示，下含鐵帶巖移數(shù)值模型為一般地下開(kāi)挖問(wèn)題，可采用Mohr-Coulomb模型作為本構(gòu)模型，此時(shí)數(shù)值模擬計(jì)算還需明確巖體單軸抗壓強(qiáng)度σmc、內(nèi)摩擦角φ和內(nèi)凝聚力c等力學(xué)參數(shù)。

用主應(yīng)力表述的Mohr-Coulomb準(zhǔn)則可表示為式(6)[13]。

(6)

單軸抗壓條件下(σ3=0)，可知巖體單軸抗壓強(qiáng)度σmc可表示為式(7)。

(7)

令k=(1+sinφ)/(1-sinφ),則式(6)可表示為式(8)。

σ1=σmc+kσ3

(8)

當(dāng)0<σ3<0.25σc時(shí)，巖體的三軸實(shí)驗(yàn)結(jié)果同時(shí)符合Hoek-Brown強(qiáng)度準(zhǔn)則與Mohr-Coulomb準(zhǔn)則[12]。將通過(guò)式(1)計(jì)算獲得的σ1、σ3數(shù)據(jù)用式(8)回歸，可得到參數(shù)k及巖體單軸抗壓強(qiáng)度σmc，即可獲得巖體內(nèi)摩擦角φ和內(nèi)凝聚力c。聯(lián)立式(7)與參數(shù)k的表達(dá)式，由三角函數(shù)關(guān)系可知，內(nèi)凝聚力c可表示為式(9)。

(9)

將表1中數(shù)據(jù)代入各式進(jìn)行計(jì)算，最終獲得的礦巖力學(xué)參數(shù)見(jiàn)表2，塌陷坑回填廢石簡(jiǎn)化為強(qiáng)度較低的巖體[17]，其力學(xué)參數(shù)可按弓長(zhǎng)嶺鐵礦東南礦區(qū)現(xiàn)場(chǎng)調(diào)查資料獲取。

表1 巖體基本質(zhì)量指標(biāo)及計(jì)算參數(shù)Table 1 Index and calculated parameters for basic quality of rock mass

表2 Hoek-Brown準(zhǔn)則估算的礦巖力學(xué)參數(shù)Table 2 Mechanical parameters of ore and rock estimated by Hoek-Brown criterion

3 塌陷坑上盤巖體的巖移規(guī)律分析

3.1 塌陷坑上盤巖體位移特征和破壞模式

由圖3可知，在地表靠近空區(qū)位移等值線為逆傾向微凸曲線，等值線由近空區(qū)到近地表傾角減小；在地表遠(yuǎn)離空區(qū)的位移等值線近似直線或折線，采深越大傾角變化越大，靠近地表等值線傾角越大(近似垂直)、越密集，表明隨著采深的增加，靠近空區(qū)的上盤巖體豎向位移的增加越大，遠(yuǎn)離空區(qū)的近地表巖體主要是以橫向位移為主，初期地表塌陷坑上盤巖體發(fā)生直線型滑移，隨著采深的增加，形成折線型破壞滑移，并在近地表處會(huì)形成近似豎直的裂縫，即隨著采深的增加，在沒(méi)有充填散體作用下，塌陷坑邊壁巖體會(huì)向坑內(nèi)片落，塌陷范圍逐漸擴(kuò)大，為上盤漸進(jìn)的破壞模式[18]。

圖3 上盤巖體位移等值線圖Fig.3 Displacement contour of the rock of hanging wall

3.2 水平變形和傾斜度分析

各監(jiān)測(cè)點(diǎn)的水平位移(圖4(a))變化，按橫坐標(biāo)112 m為分界線可劃分2個(gè)區(qū)段：0～112 m間的監(jiān)測(cè)點(diǎn)位移曲線上升，稱為上升區(qū)，112～260 m間的監(jiān)測(cè)點(diǎn)位移近乎一條水平直線，近似收斂于某一固定值，稱為收斂區(qū)。采深在280 m時(shí)，即開(kāi)采至70 m水平，位移急劇增加，水平位移最大值達(dá)778 mm。對(duì)應(yīng)各采深條件下的水平變形最大值分別為0.12 mm、0.27 mm、0.56 mm、1.13 mm、24.80 mm。隨著采深的增加，水平變形的最大值呈指數(shù)增加，當(dāng)采深達(dá)280 m，水平變形最大值驟增。借鑒煤礦中對(duì)地表移動(dòng)變形指標(biāo)[19](傾斜3 mm/m，水平變形2 mm/m)和陷落變形指標(biāo)(傾斜10 mm/m，水平變形6 mm/m)，由采深為280 m時(shí)的水平位移曲線可見(jiàn)，第3監(jiān)測(cè)點(diǎn)和第4監(jiān)測(cè)點(diǎn)間的水平變形為3.1 mm/m，超過(guò)了危險(xiǎn)移動(dòng)變形臨界值，第5監(jiān)測(cè)點(diǎn)和第4監(jiān)測(cè)點(diǎn)間的水平變形為6.9 mm/m，超過(guò)了陷落的臨界值，因此，可以第4監(jiān)測(cè)點(diǎn)作為陷落范圍的邊界，其距離開(kāi)采礦體205 m(模型橫坐標(biāo)為50 m)，已經(jīng)進(jìn)入到上含鐵帶的露天采場(chǎng)范圍。綜上，在沒(méi)有考慮充填散體作用條件下，要確保露天采場(chǎng)不受巖移影響，下含鐵帶的極限采深宜為240 m，即回采至110 m水平。

監(jiān)測(cè)點(diǎn)下沉曲線(圖4(b))也可以近似分為2個(gè)區(qū)段，第一個(gè)區(qū)段近似一條直線，而且下沉值都很小，稱為穩(wěn)定區(qū)，第二個(gè)區(qū)段為一條下傾的曲線，采深越大，下傾程度越大，稱為下沉區(qū)。當(dāng)采深280 m(即開(kāi)采至70 m水平)時(shí)，下沉最大值達(dá)-668 mm；各采深條件下的傾斜度最大值分別為0.17 mm/m、0.24 mm/m、0.37 mm/m、0.60 mm/m、4.29 mm/m。隨著采深的增加，傾斜度最大值是指數(shù)增加的過(guò)程，當(dāng)采深達(dá)280 m即開(kāi)采至70 m水平時(shí)，傾斜最大值增加到4.3 mm/m，還未達(dá)到塌陷臨界變形值。對(duì)比開(kāi)采過(guò)程中水平位移變化曲線可知，在形成塌陷后的急傾斜中厚礦體后續(xù)開(kāi)采過(guò)程中，地表巖移破壞以水平移動(dòng)為主。

圖4 無(wú)充填散體作用下監(jiān)測(cè)點(diǎn)位移變化曲線圖Fig.4 Displacement curves of monitoring point without the action of granular rock

3.3 充填散體對(duì)巖移作用分析

在采深為280 m的模型條件下，逐漸增加散體高度進(jìn)行模擬。圖5為不同監(jiān)測(cè)點(diǎn)位移隨散體充填高度的變化曲線，隨著散體充填高度的增加，地表水平位移和下沉值都隨之減小，當(dāng)散體充填高度從40 m增加至280 m過(guò)程中，對(duì)應(yīng)最大水平位移值分別為473.3 mm、210.9 mm、168.8 mm、136.5 mm、106.2 mm、79.7 mm、55.3 mm，對(duì)應(yīng)最大下沉值分別為-399.3 mm、-211.5 mm、-178.5 mm、-154.9 mm、-132.7 mm、-111.7 mm、-89.9 mm，表明散體對(duì)控制地表巖移具有顯著作用，隨散體充填高度繼續(xù)增加，地表位移減小幅度逐漸降低。

圖5 不同充填散體高度下監(jiān)測(cè)點(diǎn)位移變化曲線圖Fig.5 Displacement curves of monitoring point at different heights of granular rock

圖6為地表水平變形和傾斜變形的極大值隨散體充填量的變化曲線。由圖6可知，充填散體控制傾斜變形的作用要優(yōu)于控制水平變形，當(dāng)充填散體高度超過(guò)80 m后，繼續(xù)充填散體對(duì)控制傾斜和水平變形的作用不明顯，表明利用廢石充填對(duì)控制圍巖塌陷具有臨界值[20]。當(dāng)充填散體量為160 m時(shí)，監(jiān)測(cè)點(diǎn)最大傾斜變形為0.64 mm/m，最大水平變形為1.77 mm/m，也就是說(shuō)，下含鐵帶開(kāi)采到70 m水平，要保證露天采場(chǎng)不受影響，至少要保證下含鐵帶塌陷坑中充填有160 m高的廢石散體。當(dāng)散體充填量為280 m時(shí)，與無(wú)充填條件的位移相比分別下降92.9%、86.5%，再次說(shuō)明充填散體對(duì)控制水平移動(dòng)的作用要優(yōu)于控制地表下沉。地表水平位移和下沉最大值仍高達(dá)55.3 mm、-89.9 mm，表明充填散體對(duì)控制巖體大變形具有重要作用，但對(duì)地表巖體連續(xù)小變形移動(dòng)的控制作用不明顯。這是因?yàn)椋涮钌Ⅲw顆粒之間存在較多的孔隙而導(dǎo)致散體具有可壓縮性，隨著巖體移動(dòng)，充填散體被壓密，當(dāng)充填散體具有足夠的承壓能力后，巖體進(jìn)一步移動(dòng)將受到抑制，從而起到控制巖體大變形或者塌陷的作用[21-22]。

圖6 不同充填散體高度下的監(jiān)測(cè)點(diǎn)移動(dòng)變形曲線圖Fig.6 Movement deformation curves of monitoring points at different heights of granular rock

4 現(xiàn)場(chǎng)監(jiān)測(cè)及應(yīng)用效果

根據(jù)數(shù)值模擬結(jié)果，在不考慮散體作用時(shí)，下含鐵帶開(kāi)采巖移不影響露天采場(chǎng)的極限開(kāi)采水平為110 m水平，考慮充填散體作用時(shí)，可開(kāi)采至70 m水平，但至少要充填160 m以上的散體。綜合考慮，為了確保上含鐵帶露天采場(chǎng)的安全，確定下含鐵帶的極限開(kāi)采水平為110 m，并且需要及時(shí)向下含鐵帶塌陷坑充填廢石散體。

為觀察方案實(shí)施過(guò)程中弓長(zhǎng)嶺東南礦區(qū)地表裂縫的新生位置和發(fā)展范圍，在25#勘探線上設(shè)置A1、A2兩個(gè)位移監(jiān)測(cè)點(diǎn)(圖7)，監(jiān)測(cè)點(diǎn)分別距地采塌陷坑30 m和66 m，監(jiān)測(cè)周期為4個(gè)月。

圖7 地表位移監(jiān)測(cè)點(diǎn)布置圖Fig.7 Layout of the surface displacement monitoring points

25#勘探線上的地表監(jiān)測(cè)點(diǎn)的水平位移和垂直位移變化曲線如圖8所示。由圖8可知，下含鐵帶繼續(xù)向下開(kāi)采時(shí)，A1監(jiān)測(cè)點(diǎn)和A2監(jiān)測(cè)點(diǎn)的水平位移和垂直位移均有較大增長(zhǎng)，通過(guò)及時(shí)向下含鐵帶塌陷坑充填廢石散體，使得水平位移和垂直位移的增長(zhǎng)率逐步減小，特別是2016年11月1日之后，位移的增大率均明顯降低，此時(shí)廢石散體充填高度為108 m，這表明塌陷坑內(nèi)的充填散體可抑制上盤圍巖的采動(dòng)巖移。

圖8 25#勘探線上地表監(jiān)測(cè)點(diǎn)位移變化曲線Fig.8 Displacement curve of surface monitoring point on the 25# exploration line

根據(jù)現(xiàn)場(chǎng)調(diào)查可知，因塌陷坑上盤圍巖移動(dòng)新生的地表裂縫主要集中在A1監(jiān)測(cè)點(diǎn)和近A1側(cè)的塌陷坑邊緣之間，且近塌陷坑上盤側(cè)存在較淺深度的張裂縫，其裂縫深度和張開(kāi)程度遠(yuǎn)小于充填前。在A1監(jiān)測(cè)點(diǎn)和近A2側(cè)的露天采場(chǎng)邊緣之間，未出現(xiàn)張裂縫和突然沉降的現(xiàn)象。這表明在充填廢石散體作用下，地表巖移未對(duì)露天采場(chǎng)造成影響，下含鐵帶可開(kāi)采至110 m水平。

截至2017年底，該礦區(qū)累積向塌陷坑安全排放廢石總量5 040萬(wàn)t。下含鐵帶開(kāi)采至110 m水平的過(guò)程中，上含鐵帶露天采場(chǎng)持續(xù)安全生產(chǎn)。相較原采礦生產(chǎn)設(shè)計(jì)中地采極限開(kāi)采水平為152 m而言，剝離上含鐵帶露天采場(chǎng)的廢石來(lái)充填下含鐵帶地采塌陷坑的技術(shù)方案延長(zhǎng)了地下開(kāi)采的年限和深度，緩解了生產(chǎn)銜接的難題，保障了露天地下的協(xié)同安全開(kāi)采。

5 結(jié) 論

1)隨著采深的增加，越靠近空區(qū)巖體的沉降位移的增加量越大，遠(yuǎn)離空區(qū)的巖體主要是以橫向位移為主，初期地表塌陷坑上盤巖體發(fā)生直線滑移，隨著采深的增加，形成折線型破壞滑移，并在地表會(huì)形成豎直的裂縫，符合上盤漸進(jìn)的破壞模式。

2)地表位移曲線存在分區(qū)的特征，隨著采深的增加，地表水平變形和傾斜變形值增加，水平變形增加量遠(yuǎn)大于傾斜變形增加量，表明在形成塌陷后的急傾斜中厚礦體后續(xù)開(kāi)采過(guò)程中，地表巖移破壞以水平變形為主。

3)在不考慮散體作用條件下，當(dāng)開(kāi)采到70 m水平時(shí)，露天采場(chǎng)內(nèi)將出現(xiàn)危險(xiǎn)移動(dòng)變形；開(kāi)采到70 m水平，保證塌陷坑內(nèi)填散體高度高于160 m時(shí)，可避免露天采場(chǎng)內(nèi)出現(xiàn)危險(xiǎn)移動(dòng)變形。

4)充填散體對(duì)控制地表巖移，尤其是水平移動(dòng)作用明顯。充填散體控制巖移作用具有臨界值，當(dāng)充填散體高度達(dá)臨界值之后，控制巖移作用減弱。散體的充填對(duì)于控制地表移動(dòng)范圍作用不明顯，也說(shuō)明充填散體的主要作用在于控制巖體的大變形。