王 平 張開(kāi)新 徐 釗

(武漢科技大學(xué)資源與環(huán)境工程學(xué)院，湖北 武漢430081)

隨著經(jīng)濟(jì)的快速發(fā)展，全球?qū)ΦV產(chǎn)資源的需求量日益高漲，礦產(chǎn)資源被大量地開(kāi)采，與此同時(shí)礦山開(kāi)采所帶來(lái)的地表沉陷問(wèn)題逐漸地顯現(xiàn)出來(lái)，對(duì)人民的生命財(cái)產(chǎn)安全帶來(lái)了巨大的挑戰(zhàn)［1-2］。武鋼集團(tuán)程潮鐵礦經(jīng)過(guò)多年回采，目前，采礦作業(yè)場(chǎng)地和開(kāi)拓工程已經(jīng)進(jìn)入深部開(kāi)采，隨著礦石價(jià)格的不斷高漲，前期開(kāi)采遺留在－75 ～－150 m 水平的掛幫礦就成為了新的開(kāi)采對(duì)象，為了實(shí)現(xiàn)安全高效回采工作，有必要對(duì)掛幫礦開(kāi)采所引起的地表和覆巖移動(dòng)變形情況進(jìn)行研究，以確保掛幫礦安全回采和井下生產(chǎn)的安全進(jìn)行，為礦山取得良好的經(jīng)濟(jì)效益和社會(huì)效益。本研究以程潮鐵礦東區(qū)掛幫礦安全高效回采為研究背景，通過(guò)實(shí)驗(yàn)室相似模擬試驗(yàn)，模擬掛幫礦充填開(kāi)采的全過(guò)程，然后分析掛幫礦的開(kāi)采對(duì)地表變形和攔洪壩的影響以及覆巖移動(dòng)變形情況，以此指導(dǎo)礦山對(duì)掛幫礦的安全回采工作。

1 掛幫礦概況

武鋼集團(tuán)程潮鐵礦采用無(wú)底柱分段崩落采礦法。經(jīng)過(guò)多年回采，在礦區(qū)地表范圍內(nèi)形成了大量的塌陷坑，在其礦區(qū)東區(qū)－75 ～－150 m 水平由于歷史原因殘留部分零星小礦體未回采，由于該部分礦體存在于塌陷坑的邊緣并和攔洪壩臨近，因此我們將這部分零星礦體稱(chēng)之為“掛幫礦”。該零星礦體總共分為3 部分，將其命名為掛幫礦V1號(hào)礦體、掛幫礦V2號(hào)礦體和掛幫礦V3號(hào)礦體，礦體均為急傾斜礦體，沿走向方向長(zhǎng)度約為150 m，厚度約為20 m，礦體規(guī)模不大，其中掛幫礦V3號(hào)礦體位于地表攔洪壩水庫(kù)的下方。根據(jù)相關(guān)研究資料［3］，地表攔洪壩水庫(kù)已經(jīng)進(jìn)入地表變形區(qū)域，且在水庫(kù)的東西兩岸已經(jīng)產(chǎn)生了地表裂縫，礦山采取了注漿填堵的方式進(jìn)行了治理。因此有必要對(duì)掛幫礦開(kāi)采的安全性進(jìn)行分析，確保地表攔洪壩區(qū)域的安全和礦體安全高效回采的進(jìn)行。

2 掛幫礦開(kāi)采的物理相似模型試驗(yàn)

2.1 物理相似試驗(yàn)?zāi)Ｐ驮O(shè)計(jì)

(1)模型特點(diǎn)及選擇［4］。根據(jù)實(shí)驗(yàn)室相似材料模型的特征，結(jié)合程潮鐵礦東區(qū)掛幫礦分布范圍和工程地質(zhì)特點(diǎn)，綜合比較之后選用相似材料平面應(yīng)力模型進(jìn)行試驗(yàn)研究。

(2)試驗(yàn)剖面的選擇。通過(guò)現(xiàn)場(chǎng)考察和相關(guān)資料的查閱，并結(jié)合礦山相關(guān)工程技術(shù)人員的參考意見(jiàn)，研究決定選取的典型試驗(yàn)剖面如圖1 所示。

圖1 典型剖面圖Fig.1 Typical profile

(3)相似常數(shù)及模型尺寸的確定［5］。正確確定模型的相似比是物理相似模擬試驗(yàn)的重要一步，直接影響到5 個(gè)方面:①澆筑模型的工作量和經(jīng)濟(jì)性;②尋找合適的模型材料的難易程度;③實(shí)驗(yàn)室條件，包括模型試驗(yàn)設(shè)備、操作技術(shù)和測(cè)試技術(shù);④模擬范圍;⑤試驗(yàn)精度。為此，經(jīng)過(guò)與礦山技術(shù)人員商議，結(jié)合多方面因素綜合考慮比較分析，最終得到相似比如下:幾何相似比CL=1∶400;應(yīng)力相似比Cσ=1∶520;容重相似比Cγ=1∶1.3;時(shí)間相似比Ct=1∶20。本次相似模擬試驗(yàn)架全部采用型鋼制作，結(jié)合典型剖面圖最終確定試驗(yàn)的模型尺寸為長(zhǎng)180 cm，寬20 cm，高120 cm。

(4)地應(yīng)力模擬。由于礦體埋藏深度較淺，模型直接模擬至地表，因此不考慮水平構(gòu)造應(yīng)力的影響，只模擬自重應(yīng)力場(chǎng)。

(5)采礦過(guò)程模擬。本次相似模擬試驗(yàn)主要采用平底結(jié)構(gòu)的分段空?qǐng)鏊煤蟪涮畈傻V法向上進(jìn)行回采礦體，由于掛幫礦礦體主要賦存范圍為－ 75 ～－150 m水平，因此結(jié)合礦山的設(shè)計(jì)要求和礦山技術(shù)人員的參考意見(jiàn)，確定模擬開(kāi)采具體開(kāi)采順序是掛幫礦V1號(hào)礦體→掛幫礦V2號(hào)礦體→掛幫礦V3號(hào)礦體，從－150 m 水平高程起采，每個(gè)礦體分4 個(gè)分段從下往上依次進(jìn)行回采，開(kāi)采完1 個(gè)分段后，進(jìn)行該分段的充填工作，充填完成以后再進(jìn)行下一分段的開(kāi)采，依次類(lèi)推，直到最后一個(gè)分段的開(kāi)挖充填完成。

2.2 相似材料的研制

2.2.1 礦體及圍巖物理力學(xué)參數(shù)

巖體力學(xué)參數(shù)的選取較為復(fù)雜，目前獲取巖體力學(xué)參數(shù)的方法主要有現(xiàn)場(chǎng)或?qū)嶒?yàn)室試驗(yàn)、工程類(lèi)比法和經(jīng)驗(yàn)系數(shù)折減法等。本研究中，礦巖體的物理力學(xué)參數(shù)主要通過(guò)以下2 個(gè)方面獲取［6］:①現(xiàn)場(chǎng)和實(shí)驗(yàn)室試驗(yàn)，結(jié)合礦山提供的有關(guān)地質(zhì)資料，并對(duì)所研究區(qū)域進(jìn)行現(xiàn)場(chǎng)地質(zhì)調(diào)查，并取樣進(jìn)行了實(shí)驗(yàn)室?guī)r石力學(xué)試驗(yàn);②經(jīng)驗(yàn)系數(shù)折減法，該法主要是結(jié)合實(shí)驗(yàn)室實(shí)驗(yàn)并參照中科院武漢巖土所有關(guān)研究經(jīng)驗(yàn)將試驗(yàn)結(jié)果通過(guò)折減系數(shù)對(duì)其進(jìn)行弱化處理。最終得到實(shí)驗(yàn)?zāi)M所需的礦巖體的物理力學(xué)參數(shù)，見(jiàn)表1 所示。

表1 礦體及圍巖物理力學(xué)參數(shù)Table 1 Physical and mechanical parameters of ore body and surrounding rock

2.2.2 相似材料的研制

相似實(shí)驗(yàn)配比過(guò)程中，分別將模擬巖石的抗壓強(qiáng)度和容重轉(zhuǎn)換為相似小試塊的抗壓強(qiáng)度和容重，再進(jìn)行試塊配比試驗(yàn)，其中本次實(shí)驗(yàn)參數(shù)的確定主要是以單軸抗壓強(qiáng)度和容重作為主要試驗(yàn)相似配比的標(biāo)準(zhǔn)，為了減小實(shí)驗(yàn)誤差，對(duì)于同一配比號(hào)，做3 個(gè)試件，取3 個(gè)試件容重和單軸抗壓強(qiáng)度的平均值分別作為這組試塊的容重和抗壓強(qiáng)度。經(jīng)過(guò)多種不同組分和配比的相似材料物理力學(xué)參數(shù)的測(cè)定，最終選取:①重晶石粉、細(xì)砂、石膏、水按67∶55∶4∶20 配比作為花崗巖的相似材料;②重晶石粉、細(xì)砂、石膏、水按50∶63∶2∶17 配比作為矽卡巖的相似材料;③重晶石粉、細(xì)砂、石膏、水按42∶68∶2∶15 配比作為閃長(zhǎng)巖的相似材料;④重晶石粉、細(xì)砂、石膏、水按35∶47∶1.5∶13 配比作為大理巖的相似材料;⑤重晶石粉、細(xì)砂、石膏、水按42∶68∶2∶15 配比作為亞砂土的相似材料;⑥重晶石粉、鐵礦粉、細(xì)砂、石膏、水按21∶78∶40∶8∶16 配比作為磁鐵礦的相似材料。

2.3 監(jiān)測(cè)手段及監(jiān)測(cè)內(nèi)容

為了分析地表位移隨工作面的推進(jìn)相應(yīng)沉降變化情況，結(jié)合實(shí)驗(yàn)室現(xiàn)有試驗(yàn)條件，在地表共放置5個(gè)百分表，其中地表攔洪壩由于是本次試驗(yàn)研究的重點(diǎn)觀察部位，所以該處設(shè)置3 個(gè)百分表，從左到右依次編號(hào)為1、2、3、4、5 號(hào)百分表。同時(shí)為有效反映隨著工作面的推進(jìn)礦體上覆巖層的應(yīng)力應(yīng)變變化，根據(jù)所研究區(qū)域的研究重點(diǎn)，盡可能地捕捉地下開(kāi)挖引起的地表變形特點(diǎn)，因此在地表攔洪壩的下部采用10 cm×15 cm 的網(wǎng)格式布置，布置3 排，且在掛幫礦V1、V2號(hào)礦體上方均勻布置7 個(gè)監(jiān)測(cè)點(diǎn)，具體監(jiān)測(cè)點(diǎn)編號(hào)布置見(jiàn)圖2 所示。在上述編號(hào)點(diǎn)位置用502 膠水作為黏接劑相應(yīng)地黏貼上電阻應(yīng)變片(規(guī)格BF －120 － 40AA)，將電阻應(yīng)變片通過(guò)導(dǎo)線(xiàn)連接上BX2011C1 型力學(xué)綜合參數(shù)測(cè)試儀，同時(shí)連接上計(jì)算機(jī)進(jìn)行相關(guān)調(diào)試，為正式測(cè)試做準(zhǔn)備。

圖2 模型監(jiān)測(cè)點(diǎn)布置Fig.2 Layout of monitoring point in the model

2.4 試驗(yàn)基本步驟

在布置好觀測(cè)點(diǎn)及設(shè)備以后，記錄初始讀數(shù)，即可開(kāi)始進(jìn)行開(kāi)挖與沖填。本次試驗(yàn)?zāi)M的開(kāi)挖充填順序依據(jù)圖1 中3 部分礦體依次分為3 步進(jìn)行回采，開(kāi)挖先后順序?yàn)閽鞄偷VV1號(hào)礦體→掛幫礦V2號(hào)礦體→掛幫礦V3號(hào)礦體，每個(gè)礦體再按照礦山設(shè)計(jì)參數(shù)分為4 個(gè)分段進(jìn)行開(kāi)采，從下往上，依次分為－130～－150 m，－110 ～－130 m，－90 ～－110 m，－75 ～－90 m 計(jì)4 個(gè)分段。整個(gè)試驗(yàn)將礦體分為12 個(gè)步驟進(jìn)行開(kāi)采，第1 ～4 步依次從下向上開(kāi)采V1號(hào)礦體，第5 ～8 步依次從下向上開(kāi)采V2號(hào)礦體，第9 ～12 步依次從下向上開(kāi)采V3號(hào)礦體，每一次開(kāi)挖和充填完成以后相應(yīng)地記錄百分表和力學(xué)綜合參數(shù)測(cè)試儀的數(shù)據(jù)，為后期地表位移和覆巖應(yīng)力應(yīng)變分析做好準(zhǔn)備。

3 試驗(yàn)結(jié)果與分析

3.1 掛幫礦開(kāi)采對(duì)地表變形的影響分析

根據(jù)在模型地表布置的5 個(gè)百分表測(cè)點(diǎn)的觀測(cè)結(jié)果，繪制出地表沉降位移隨開(kāi)采步數(shù)的變化曲線(xiàn)如圖3 所示。

圖3 地表沉降位移與開(kāi)挖步數(shù)變化曲線(xiàn)Fig.3 Variation curve of surface subsidence displacement with excavation steps

由圖3 可以看出，在開(kāi)挖至第4 步即V1號(hào)礦體－150 ～－75 m 段開(kāi)挖完畢之后，地表除了礦體上方的百分表5 測(cè)點(diǎn)下沉量明顯以外，其他各測(cè)點(diǎn)的下沉位移量很小，其中百分表5 最大測(cè)量值為－0.21 mm(相當(dāng)于工程原型的8.4 cm)。在開(kāi)挖至第8 步即V2號(hào)礦體－150 ～－75 m 段開(kāi)挖完畢之后，地表布置的百分表4 沉降量有明顯增加，結(jié)合圖2 充分說(shuō)明了礦體地下開(kāi)采所造成的地表沉陷區(qū)域主要位于采空區(qū)正上方，并以采空區(qū)為中心向外發(fā)展。當(dāng)?shù)?2 步即V3號(hào)礦體－150 ～－75 m 段開(kāi)挖完畢之后，發(fā)現(xiàn)隨著開(kāi)挖步次的不斷推進(jìn)，地表各監(jiān)測(cè)點(diǎn)沉降位移都有增加，其中地表最大沉降位移發(fā)生在5 號(hào)百分表測(cè)點(diǎn)上，最大沉降位移為－0.35 mm(相當(dāng)于工程原型的14 cm)。由于百分表5 位于V1號(hào)礦體上方，從其位移變化數(shù)據(jù)可以看出，隨著開(kāi)挖步次的增加，其沉降位移在緩慢地增大，說(shuō)明地表沉降明顯滯后于采礦，地表沉降的時(shí)間效應(yīng)非常明顯。

3.2 掛幫礦開(kāi)采對(duì)攔洪壩的影響分析

通過(guò)圖3 百分表記錄的數(shù)據(jù)可以看出，隨著掛幫礦的不斷開(kāi)采，地表的影響范圍在逐漸地增大。其中對(duì)攔洪壩區(qū)域的監(jiān)測(cè)主要是利用百分表1、2、3 的實(shí)時(shí)監(jiān)測(cè)，從圖3 中變化曲線(xiàn)趨勢(shì)還可以知道，當(dāng)掛幫礦V1和V2號(hào)礦體開(kāi)采完畢時(shí)，攔洪壩區(qū)域地表沉降較小，其最大沉降值發(fā)生在第8 步開(kāi)采結(jié)束后，此時(shí)最大沉降值為－0.05 mm(相當(dāng)于工程原型的2.0 cm)。當(dāng)繼續(xù)開(kāi)挖掛幫礦V3號(hào)礦體時(shí)，百分表1、2、3的數(shù)據(jù)都有突變，這表明掛幫礦V3號(hào)礦體的開(kāi)采對(duì)地表攔洪壩的沉降產(chǎn)生了較大影響，隨著V3號(hào)礦體全部開(kāi)采完畢，地表攔洪壩范圍內(nèi)最大沉降值達(dá)到－0.26 mm(相當(dāng)于工程原型的10.4 cm)，對(duì)地表攔洪壩水庫(kù)的安全運(yùn)行構(gòu)成了威脅。

3.3 上覆巖層應(yīng)變變化特征分析

為了監(jiān)測(cè)地下掛幫礦開(kāi)采對(duì)上覆巖層擾動(dòng)的應(yīng)變變化規(guī)律，在巖層的關(guān)鍵位置設(shè)置了應(yīng)變片進(jìn)行監(jiān)測(cè)，利用應(yīng)變片可以監(jiān)測(cè)巖層變化的微小應(yīng)變。隨著礦體的開(kāi)挖不斷推進(jìn)，各個(gè)測(cè)點(diǎn)的相對(duì)位置逐漸在變化，各測(cè)點(diǎn)的壓力變化不同，測(cè)得的應(yīng)變值也不同，利用力學(xué)綜合參數(shù)測(cè)試儀采集應(yīng)變數(shù)據(jù)，各監(jiān)測(cè)點(diǎn)隨開(kāi)挖步數(shù)推進(jìn)的應(yīng)變變化規(guī)律如圖4 所示。

圖4 各監(jiān)測(cè)點(diǎn)應(yīng)變變化Fig.4 Strain figure of monitoring points

由圖4(a)可以觀察到當(dāng)掛幫礦從第1 步開(kāi)采至第8 步即V1和V2號(hào)礦體開(kāi)采完畢時(shí)，攔洪壩區(qū)域以下的監(jiān)測(cè)點(diǎn)應(yīng)變值變化很小，幾乎不變，結(jié)合監(jiān)測(cè)點(diǎn)布置圖可知，由于掛幫礦V1和V2號(hào)礦體遠(yuǎn)離攔洪壩，因此開(kāi)采掛幫礦V1和V2號(hào)礦體對(duì)攔洪壩區(qū)域的影響較小。從圖4(a)中拐點(diǎn)出現(xiàn)的開(kāi)挖步數(shù)可知，當(dāng)開(kāi)采至第9 步即掛幫礦V3號(hào)礦體開(kāi)始回采時(shí)，攔洪壩區(qū)域以下的監(jiān)測(cè)點(diǎn)應(yīng)變突然變大，說(shuō)明掛幫礦V3號(hào)礦體的開(kāi)采對(duì)攔洪壩產(chǎn)生了顯著影響，當(dāng)開(kāi)挖至第12 步時(shí)，位于礦體正上方的最近2 排監(jiān)測(cè)點(diǎn)的應(yīng)變值從負(fù)值變化為正值，即表示其應(yīng)力變化由壓應(yīng)力突然轉(zhuǎn)化為拉應(yīng)力，拉應(yīng)力的出現(xiàn)對(duì)開(kāi)采安全產(chǎn)生了影響，容易造成冒頂片幫事故的發(fā)生，當(dāng)拉應(yīng)力達(dá)到巖石的極限抗拉強(qiáng)度時(shí)，圍巖會(huì)發(fā)生破壞，可能造成攔洪壩區(qū)域的節(jié)理裂隙的不斷發(fā)育，最終連通攔洪壩水庫(kù)可能造成淹井事故的發(fā)生。

由圖4(b)中監(jiān)測(cè)點(diǎn)應(yīng)變變化曲線(xiàn)可知，當(dāng)掛幫礦V1號(hào)礦體開(kāi)采過(guò)程中也即開(kāi)挖步次從第1 步至第4 步，均是測(cè)點(diǎn)2 －15 的應(yīng)變值最大，且2 －13、2 －14、2 －15 測(cè)點(diǎn)的應(yīng)變值在同步次監(jiān)測(cè)值里都是較大的，表明礦體開(kāi)挖主要影響區(qū)域是礦體開(kāi)采周?chē)鷶_動(dòng)區(qū)域，離采空區(qū)越遠(yuǎn)的監(jiān)測(cè)點(diǎn)，受到的擾動(dòng)越小，應(yīng)變值越小。當(dāng)掛幫礦V2號(hào)礦體回采過(guò)程中開(kāi)挖至第8步次時(shí)，測(cè)點(diǎn)2 －12、2 －13 應(yīng)變值陡然上升，說(shuō)明隨著礦體開(kāi)挖越接近地表，其頂板上部會(huì)出現(xiàn)急劇的應(yīng)力集中，在開(kāi)采過(guò)程中，一定要注意對(duì)頂板的及時(shí)監(jiān)測(cè)和處理。

4 掛幫礦開(kāi)采安全隱患及防治措施

通過(guò)對(duì)實(shí)驗(yàn)室相似模擬試驗(yàn)結(jié)果的分析，并綜合考慮礦山地質(zhì)、掛幫礦礦體特征、開(kāi)采技術(shù)條件等因素，結(jié)合現(xiàn)場(chǎng)實(shí)際勘察情況分析掛幫礦開(kāi)采存在的安全隱患主要包括:①水庫(kù)發(fā)生潰壩事故;②井下淹井事故;③地表突然沉陷災(zāi)害;④冒頂片幫事故等。根據(jù)掛幫礦開(kāi)采潛在的安全隱患，特制定掛幫礦開(kāi)采安全防治措施如下:①對(duì)掛幫礦采用充填法進(jìn)行開(kāi)采;②加強(qiáng)對(duì)水庫(kù)安全監(jiān)測(cè)和地表裂縫注漿封堵處理;③對(duì)掛幫礦采取試探開(kāi)采的原則并合理控制回采進(jìn)度;④礦山組織設(shè)計(jì)應(yīng)急防水避災(zāi)線(xiàn)路并定期進(jìn)行井下避災(zāi)演習(xí);⑤加強(qiáng)地表移動(dòng)、變形監(jiān)測(cè)工作，達(dá)到提前預(yù)防災(zāi)害的目的。

5 結(jié) 論

(1)通過(guò)實(shí)驗(yàn)室相似模擬試驗(yàn)結(jié)果可知，隨著掛幫礦礦體的回采，地表變形范圍不斷加大，采空區(qū)上方的圍巖產(chǎn)生變形，局部地方由壓應(yīng)力轉(zhuǎn)化為拉應(yīng)力，將對(duì)地表攔洪壩的安全運(yùn)行產(chǎn)生危害。

(2)當(dāng)掛幫礦V1和V2號(hào)礦體開(kāi)采時(shí)，攔洪壩區(qū)域地表沉降和應(yīng)變變化較小，其中掛幫礦V3號(hào)礦體的開(kāi)采會(huì)導(dǎo)致攔洪壩的變形以及攔洪壩與V3號(hào)礦體之間的圍巖變形急劇增加，使攔洪壩與V3號(hào)礦體之間的裂隙進(jìn)一步擴(kuò)展，危及程潮鐵礦主采區(qū)的安全，因此建議V3號(hào)礦體不予開(kāi)采，V1和V2號(hào)礦體應(yīng)該按照開(kāi)采技術(shù)規(guī)范要求安全回采。

(3)礦山應(yīng)采取合理的安全防范措施保證掛幫礦的安全生產(chǎn)，建議采用充填法對(duì)礦體進(jìn)行回采，平時(shí)加強(qiáng)對(duì)攔洪壩水庫(kù)水位的監(jiān)測(cè)，及時(shí)對(duì)地裂縫進(jìn)行注漿封堵處理，并采用試探開(kāi)采的原則全面保證礦體的安全回采。

［1］ 余學(xué)義，張恩強(qiáng).開(kāi)采沉陷學(xué)［M］.2 版.北京:煤炭工業(yè)出版社，2010.

Yu Xueyi，Zhang Enqiang. Subsidence Science［M］. 2nd ed. Beijing:Coal Mining Industry Press，2010.

［2］ 李向陽(yáng)，李俊平，周創(chuàng)兵，等. 采空?qǐng)龈矌r變形數(shù)值模擬與相似模擬比較研究［J］.巖土力學(xué)，2005，26(12):1907-1912.

Li Xiangyang，Li Junping，Zhou Chuangbing，et al.Comparative study on numerical simulation and similarity simulation of overburden deformation in abandoned stope［J］. Rock and Soil Mechanics，2005，26(12):1907-1912.

［3］ 中國(guó)科學(xué)院武漢巖土力學(xué)研究所.程潮鐵礦塌陷區(qū)周邊地表變形監(jiān)測(cè)和預(yù)測(cè)［R］. 武漢:中國(guó)科學(xué)院武漢巖土力學(xué)研究所，2013.

Institute of Rock and Soil Mechanics，Chinese Academy of Sciences.Chengchao Iron Mine Subsidence Area Deformation Monitoring and Prediction［R］. Wuhan:Institute of Rock and Soil Mechanics，Chinese Academy of Sciences，2013.

［4］ 蔡嗣經(jīng)，陳清運(yùn)，明世祥.金山店鐵礦平行礦體地下開(kāi)采地表沉降物理模擬預(yù)測(cè)研究［J］.中國(guó)安全生產(chǎn)科學(xué)技術(shù)，2006，2(5):13-19.

Cai Sijing，Chen Qingyun，Ming Shixiang. Study on subsidence prediction of underground mining of parallel ore bodies in Jinshandian mine by physical modeling［J］.Journal of Safety Science and Technology，2006，2(5):13-19.

［5］ 馬 俊，許夢(mèng)國(guó)，王 平，等. 無(wú)底柱分段崩落法塌陷區(qū)回填體移動(dòng)規(guī)律研究［J］.金屬礦山，2011(8):37-40.

Ma Jun，Xu Mengguo，Wang Ping，et al.Study on the movement laws of subsided backfilling body with pillarless sublevel caving［J］.Metal Mine，2011(8):37-40.

［6］ 張小義.金山店鐵礦深部開(kāi)采地壓活動(dòng)規(guī)律研究及回采順序優(yōu)化［D］.武漢:武漢科技大學(xué)，2013.

Zhang Xiaoyi.The Study of the Rule of Underground Pressure Activity and Optimization of Stoping Sequence in the Deep Mining of Jinshandian Iron Mine［D］. Wuhan:Wuhan University of Science and Technology，2013.

［7］ 陳順滿(mǎn)，許夢(mèng)國(guó)，王 平，等. 某鐵礦掛幫礦開(kāi)采對(duì)地表攔洪壩的影響［J］. 金屬礦山，2014(2):150-153.

Chen Shunman，Xu Mengguo，Wang Ping，et al.Stability analysis of the effect of the hanging wall ore excavation on earth's surface floodcontrol dam in an iron mine［J］.Metal Mine，2014(2):150-153.