鄔 金 李元輝 司呈斌 徐 帥

(1.深部金屬礦山安全開(kāi)采教育部重點(diǎn)實(shí)驗(yàn)室，遼寧 沈陽(yáng) 110819；2.東北大學(xué)資源與土木工程學(xué)院，遼寧 沈陽(yáng)110819；3.遼寧廣播電視大學(xué)人文學(xué)院，遼寧 沈陽(yáng) 110161；4.山東黃金礦業(yè)(萊州)有限公司焦家金礦，山東 煙臺(tái) 261441)

深埋厚大礦體采場(chǎng)結(jié)構(gòu)參數(shù)優(yōu)化

鄔 金1,2,3李元輝1,2司呈斌4徐 帥1,2

(1.深部金屬礦山安全開(kāi)采教育部重點(diǎn)實(shí)驗(yàn)室，遼寧 沈陽(yáng) 110819；2.東北大學(xué)資源與土木工程學(xué)院，遼寧 沈陽(yáng)110819；3.遼寧廣播電視大學(xué)人文學(xué)院，遼寧 沈陽(yáng) 110161；4.山東黃金礦業(yè)(萊州)有限公司焦家金礦，山東 煙臺(tái) 261441)

合理的采場(chǎng)結(jié)構(gòu)參數(shù)可使采場(chǎng)處于有利的力學(xué)狀態(tài)，使圍巖的應(yīng)力、應(yīng)變分布趨于均勻化，在保證開(kāi)采系統(tǒng)穩(wěn)定和生產(chǎn)安全的前提下，減少支護(hù)工作量，提高采礦強(qiáng)度和生產(chǎn)效率。在深入分析思山嶺鐵礦地質(zhì)概況與采礦方法的基礎(chǔ)上，對(duì)影響礦房回采穩(wěn)定性的礦房高度、礦房寬度、采場(chǎng)長(zhǎng)度、礦柱寬度、礦柱充填方式等5個(gè)關(guān)鍵因素進(jìn)行2水平正交設(shè)計(jì)，獲得8種試驗(yàn)方案。運(yùn)用大型巖土軟件FLAC3D對(duì)盤(pán)區(qū)內(nèi)不同方案的采場(chǎng)結(jié)構(gòu)參數(shù)進(jìn)行數(shù)值模擬研究，分析其在不同結(jié)構(gòu)參數(shù)下應(yīng)力、位移、塑性區(qū)等特征，初步得出采場(chǎng)處于最有利力學(xué)狀態(tài)時(shí)的結(jié)構(gòu)參數(shù)方案(采場(chǎng)高60 m、采場(chǎng)長(zhǎng)60 m，礦房寬18 m、礦柱寬20 m的參數(shù)方案)。計(jì)算結(jié)果表明：回采過(guò)程中，采場(chǎng)長(zhǎng)度對(duì)頂板應(yīng)力和頂板位移的影響最大，采場(chǎng)越長(zhǎng)，應(yīng)力值越大，且壓應(yīng)力主要在盤(pán)區(qū)間柱集中，頂?shù)装逄幊霈F(xiàn)拉應(yīng)力集中。分析結(jié)果可為盤(pán)區(qū)礦房礦柱的安全高效回采提供技術(shù)支持。

正交設(shè)計(jì) 深部厚大礦體 階段礦房嗣后充填法 結(jié)構(gòu)參數(shù)優(yōu)化 數(shù)值模擬

礦產(chǎn)資源是經(jīng)濟(jì)發(fā)展的基礎(chǔ)，國(guó)民經(jīng)濟(jì)高速發(fā)展與淺部資源的逐漸枯竭之間的矛盾，促使深埋厚大礦體開(kāi)采成為未來(lái)礦業(yè)發(fā)展的重要方向[1-2]。深埋厚大礦體處于較高的應(yīng)力環(huán)境下，為強(qiáng)化開(kāi)采而進(jìn)行的深孔爆破等采礦行為勢(shì)必導(dǎo)致高應(yīng)力場(chǎng)的重新分布，產(chǎn)生新的應(yīng)力集中和降低區(qū)，導(dǎo)致巷道和采場(chǎng)周?chē)鷰r體變形、破壞，極易誘發(fā)巖爆、大體積塌方等工程災(zāi)害，危及人員設(shè)備安全。地壓活動(dòng)誘發(fā)的工程災(zāi)害與采礦方法、采場(chǎng)結(jié)構(gòu)參數(shù)、開(kāi)采強(qiáng)度、回采順序等緊密相關(guān)，尤其是采場(chǎng)結(jié)構(gòu)參數(shù)，對(duì)圍巖二次應(yīng)力的分布影響較大。合理的采場(chǎng)結(jié)構(gòu)參數(shù)可使采場(chǎng)處于有利的力學(xué)狀態(tài)，使圍巖的應(yīng)力、應(yīng)變分布趨于均勻化，保證開(kāi)采系統(tǒng)穩(wěn)定和生產(chǎn)安全的前提下，減少支護(hù)工作量，提高采礦強(qiáng)度和生產(chǎn)效率，降低礦石損失與貧化，實(shí)現(xiàn)礦山安全、高效、經(jīng)濟(jì)開(kāi)采。

目前國(guó)內(nèi)外礦山生產(chǎn)過(guò)程中優(yōu)化采場(chǎng)結(jié)構(gòu)參數(shù)的主要方法有：工程類(lèi)比法[4]、半經(jīng)驗(yàn)半理論方法、遺傳算法[5]、數(shù)值模擬法[6-7]、神經(jīng)網(wǎng)絡(luò)優(yōu)化等方法[8-9]。其中，隨著計(jì)算機(jī)的發(fā)展和計(jì)算技術(shù)的提高，數(shù)值分析方法現(xiàn)已成為解決地下工程問(wèn)題的有力工具，在采礦工程中得到廣泛應(yīng)用，目前應(yīng)用比較多的數(shù)值計(jì)算軟件是基于有限元思想的ANSYS和基于有限差分思想的FLAC3D。如鄂西緩薄高磷赤鐵礦的采場(chǎng)結(jié)構(gòu)參數(shù)優(yōu)化、三鑫金銅桃花嘴礦區(qū)厚大礦體采場(chǎng)結(jié)構(gòu)參數(shù)優(yōu)化研究、盤(pán)區(qū)隔離礦柱采場(chǎng)結(jié)構(gòu)參數(shù)數(shù)值優(yōu)化、新城金礦深部采場(chǎng)結(jié)構(gòu)參數(shù)優(yōu)化研究等[8-14]。在國(guó)外，采場(chǎng)結(jié)構(gòu)優(yōu)化問(wèn)題早就引起來(lái)人們的廣泛關(guān)注。在解決采場(chǎng)結(jié)構(gòu)參數(shù)問(wèn)題上，主要采用經(jīng)驗(yàn)法和計(jì)算機(jī)模擬法來(lái)進(jìn)行優(yōu)化，進(jìn)一步通過(guò)礦山現(xiàn)場(chǎng)試驗(yàn)對(duì)結(jié)果進(jìn)行檢驗(yàn)。

1 思山嶺鐵礦采礦方法

1.1 采礦方法選擇

思山嶺鐵礦埋深404～1 934 m，東西長(zhǎng)1 500 m，南北寬960 m，垂直深度1 580 m，礦體平均厚246.84 m，傾角約70°，屬隱伏盲礦床。區(qū)內(nèi)巖性較簡(jiǎn)單，巖石較完整，礦體呈厚層狀產(chǎn)出。礦體賦存于赤鐵石英巖與磁鐵石英巖中，頂板較完整，穩(wěn)定性較好。依據(jù)礦山設(shè)計(jì)，一期開(kāi)采深度為-960 m中段以下磁鐵礦，埋深1 400 m，產(chǎn)能1 500萬(wàn)t/a，采用大直徑深孔階段礦房嗣后充填法進(jìn)行開(kāi)采，采礦方法如圖1所示。

圖1 采礦方法示意

1.2 礦塊和采場(chǎng)布置

設(shè)計(jì)將礦體劃分為盤(pán)區(qū)，盤(pán)區(qū)是回采生產(chǎn)基本單元。盤(pán)區(qū)長(zhǎng)230 m，寬200 m。盤(pán)區(qū)中部垂直走向留20 m寬的間柱，間柱內(nèi)布置穿脈巷道，在盤(pán)區(qū)端部留10 m間柱，間柱盡量布置在夾石或者赤鐵礦中。盤(pán)區(qū)間沿礦體走向布置整體條形礦柱，其中布置盤(pán)區(qū)沿脈巷道和溜井。每個(gè)盤(pán)區(qū)內(nèi)劃分礦房和礦柱。

回采分兩步進(jìn)行，一步驟回采礦房，二步驟回采礦柱。由下而上梯段式分段側(cè)向崩礦，爆下的礦石用25 t電動(dòng)鏟運(yùn)機(jī)在采場(chǎng)底部出礦。待礦房礦石全部出完后，進(jìn)行膠結(jié)或尾砂充填，待兩面或三面礦房采完充填好并養(yǎng)護(hù)2個(gè)月后，開(kāi)采礦柱，礦柱的回采方式與礦房相同，回采結(jié)束后進(jìn)行充填。

2 優(yōu)化方案及數(shù)值模擬

2.1 影響礦房回采穩(wěn)定性因素的分析

正交實(shí)驗(yàn)設(shè)計(jì)是實(shí)驗(yàn)設(shè)計(jì)中一種優(yōu)良的方法，它是指用正交表合理地安排實(shí)驗(yàn)，利用數(shù)理統(tǒng)計(jì)的原理科學(xué)地分析實(shí)驗(yàn)結(jié)果，是一種處理多因素實(shí)驗(yàn)的科學(xué)方法。影響礦房回采穩(wěn)定性的因素很多，根據(jù)礦山實(shí)際，本次選取礦房的高度、礦房寬度、采場(chǎng)長(zhǎng)度、礦柱寬度、礦柱充填方式等5方面的因素?？紤]以上5個(gè)因素，各因素所選取的水平見(jiàn)表1。

表1 影響因素選取及其水平

注：A為礦房的高度，m；B為礦房寬度，m；C為礦柱寬度，m；D為采場(chǎng)長(zhǎng)度，m；E為礦柱充填方式。

2.2 正交試驗(yàn)?zāi)M方案

選取了影響礦房穩(wěn)定性的5因素2水平，采用正交分析來(lái)安排模擬試驗(yàn)，選用L8(27)比較合適。需要做8次模擬試驗(yàn)，試驗(yàn)方案如表2所示。

表2 正交分析的所有模擬試驗(yàn)組號(hào)

2.3 數(shù)值模擬研究

2.3.1 三維有限元模型的建立

為確定合理的采場(chǎng)結(jié)構(gòu)參數(shù)，根據(jù)正交試驗(yàn)原理設(shè)計(jì)了如表2所示的8種方案，并基于AutoCAD建立采場(chǎng)三維模型，采用ANSYS與FLAC3D耦合的方式進(jìn)行模型的數(shù)值計(jì)算。思山嶺鐵礦礦巖均屬于脆性硬巖，因此本構(gòu)模型采用彈塑性模型，塑性屈服條件為摩爾-庫(kù)倫準(zhǔn)則。巖石力學(xué)參數(shù)通過(guò)試驗(yàn)取得，如表3所示，計(jì)算模型如圖2所示。

2.3.2 模擬結(jié)果分析

本次模擬是在盤(pán)區(qū)礦房礦柱均已開(kāi)采充填之后，分析其頂板的應(yīng)力狀態(tài)及位移狀態(tài)。各方案正交試驗(yàn)的結(jié)果(頂板最大主應(yīng)力、最小主應(yīng)力比較及頂板位移等)如表4所示，對(duì)正交試驗(yàn)的結(jié)果進(jìn)行極差分析如表5所示。

表3 思山嶺鐵礦主要巖體力學(xué)參數(shù)

圖2 計(jì)算分析模型

表4 正交試驗(yàn)結(jié)果

注：σmax、σmin分別為頂板最大主應(yīng)力和頂板最小主應(yīng)力；ε為頂板中點(diǎn)的位移。

從表4中各個(gè)方案的最大主應(yīng)力、最小主應(yīng)力與垂向位移可以看出，方案二與其余方案相比，具有較明顯的優(yōu)勢(shì)，能夠保證采場(chǎng)圍巖的穩(wěn)定性。該方案的結(jié)構(gòu)參數(shù)為：采場(chǎng)高60m、采場(chǎng)長(zhǎng)度60m，礦房寬18m、礦柱寬20m。

從頂板最大主應(yīng)力云圖(從略)可以看出礦房開(kāi)挖完畢后，應(yīng)力向采場(chǎng)周?chē)小Ｓ杀?可以看出，方案5的最大主應(yīng)力值最大，達(dá)到92.53MPa。

考慮各個(gè)因素的極差能夠反應(yīng)試驗(yàn)指標(biāo)的變動(dòng)幅度，進(jìn)而說(shuō)明，極差越大，該因素對(duì)試驗(yàn)指標(biāo)的影響越大，并由此可以判斷出因素的主次順序。從表5最大主應(yīng)力的極差分析結(jié)果可以看出，影響礦體回采頂板穩(wěn)定性的5個(gè)因素的極差由大到小分別為采場(chǎng)長(zhǎng)度、礦柱寬度、充填方式、礦房寬度、礦房高，即采場(chǎng)長(zhǎng)度對(duì)采場(chǎng)回采頂板穩(wěn)定性影響最為明顯，而礦房高度影響最小。由于礦房的高度對(duì)頂板的暴露面積沒(méi)有作用，所以礦房高度的影響必然最小。而在影響暴露面積的幾個(gè)因素中，礦柱的寬度的影響明顯大于礦房的寬度。這可能是因?yàn)榈V房開(kāi)挖時(shí)是在原巖中開(kāi)挖，而開(kāi)挖礦柱時(shí)兩側(cè)被充填體所替代，所以造成礦柱的寬度對(duì)頂板的影響明顯大于礦房的寬度。因而在設(shè)計(jì)采場(chǎng)時(shí)應(yīng)該盡量少增加礦柱的寬度。

表5 正交試驗(yàn)極差分析結(jié)果

礦房回采后，礦房出現(xiàn)了拉應(yīng)力集中，且主要集中在盤(pán)區(qū)礦房的頂?shù)装逄?，礦柱充填體也出現(xiàn)了拉應(yīng)力。從表4中可以看出方案3的拉應(yīng)力值最大，達(dá)到4.32 MPa，但是沒(méi)有達(dá)到礦體的抗拉強(qiáng)度，所以采場(chǎng)處于安全范圍。但由于巖石抗壓不抗拉的特性，采場(chǎng)的拉應(yīng)力對(duì)采場(chǎng)的穩(wěn)定性有很大的影響。

從表5的最小主應(yīng)力的極差分析結(jié)果可以看出，影響采場(chǎng)回采拉應(yīng)力大小的5個(gè)因素的極差值由大到小分別為采場(chǎng)長(zhǎng)度、礦房的高度、礦柱充填方式、礦房寬度、礦柱寬度。對(duì)采場(chǎng)拉應(yīng)力影響最大的因素是采場(chǎng)的長(zhǎng)度，其次是礦房的高度，礦房的寬度和礦柱的寬度影響均較小。從頂板位移云圖(從略)可以看出，盤(pán)區(qū)中部采場(chǎng)的位移明顯大于盤(pán)區(qū)端部采場(chǎng)的位移，并且盤(pán)區(qū)間柱位置也處在較大位移區(qū)域。由表4可以看出，方案5的頂板位移最大，達(dá)到了19.49 cm，方案2最小。在圍巖位移狀態(tài)中，垂直位移比水平位移重要，因?yàn)檫^(guò)量的垂直位移將導(dǎo)致頂板冒落和地表沉陷。因此，從此時(shí)位移狀態(tài)看，方案2為最好，方案5最差。

從表5的頂板垂直位移的極差分析結(jié)果可以看出，影響礦體回采頂板位移的5個(gè)因素的極差值由大到小分別為采場(chǎng)長(zhǎng)度、礦房高度、礦柱寬度、礦房寬度、礦柱充填方式，即采場(chǎng)長(zhǎng)度對(duì)采場(chǎng)回采頂板位移影響最為明顯，而第二步是否采用膠結(jié)充填影響最小。礦柱寬度的影響明顯要大于礦房的寬度。礦房的高度對(duì)采場(chǎng)頂板位移的影響明顯要大于對(duì)應(yīng)力的影響。

3 結(jié) 論

(1)基于正交實(shí)驗(yàn)建立了5因素2水平8種不同采場(chǎng)結(jié)構(gòu)參數(shù)模型，開(kāi)展數(shù)值模擬分析，并針對(duì)數(shù)值模擬的結(jié)果分別從最大主應(yīng)力、位移、拉應(yīng)力等角度進(jìn)行了比較分析。初步得出了采場(chǎng)高60 m、采場(chǎng)長(zhǎng)度60 m，礦房寬18 m、礦柱寬20 m的參數(shù)方案為思山嶺鐵礦大直徑深孔嗣后膠結(jié)充填法較優(yōu)參數(shù)。

(2)正交試驗(yàn)分析表明對(duì)頂板拉應(yīng)力影響最大的因素是采場(chǎng)長(zhǎng)度，其次是礦房高度、礦房的寬度和礦柱的寬度影響均較小。對(duì)于頂板壓應(yīng)力分布，礦柱寬度的影響明顯大于礦房的寬度。因此，在采場(chǎng)結(jié)構(gòu)參數(shù)選取過(guò)程中應(yīng)該選擇較小的礦房長(zhǎng)度和礦柱寬度。同時(shí)為了達(dá)到產(chǎn)量要求可以選擇較大的礦房寬度和高度。

(3)在回采過(guò)程中，采場(chǎng)頂板出現(xiàn)了拉應(yīng)力集中區(qū)，雖沒(méi)達(dá)到巖體抗拉強(qiáng)度，但也是構(gòu)成采場(chǎng)不穩(wěn)定的重要因素，在采場(chǎng)設(shè)計(jì)過(guò)程中，應(yīng)適當(dāng)減小采場(chǎng)長(zhǎng)度，減少拉應(yīng)力的產(chǎn)生。

(4)在礦房開(kāi)挖過(guò)程中，壓應(yīng)力向盤(pán)區(qū)間柱集中，由于間柱中布置有穿脈巷道，應(yīng)對(duì)巷道采取必要的支護(hù)措施。

[1] 古德生,李夕兵,等.現(xiàn)代金屬礦床開(kāi)采科學(xué)技術(shù)[M].北京:冶金工業(yè)出版社,2006. Gu Desheng，Li Xibing，et al.Modern Mining Science and Technology for Metal Mineral Resources[M].Beijing：Metallurgical Industry Press，2006.

[2] 何滿潮,謝和平,彭蘇萍,等.深部開(kāi)采巖體力學(xué)研究[J].巖石力學(xué)與工程學(xué)報(bào),2005,24(16)：2803-2813. He Manchao，Xie Heping，Peng Suping，et al.Study on rock mechanics in deep mining engineering[J].Chinese Journal of Rock Mechanics and Engineering,2005,24(16)：2803-2813.

[3] 盧 萍.深部采場(chǎng)結(jié)構(gòu)參數(shù)及回采順序優(yōu)化研究[D].重慶：重慶大學(xué),2008. Lu Ping.Optimization of Deep stope Structure Parameters and Mining Sequence[D].Chongqing:Chongqing University，2008.

[4] 朱維申,李術(shù)才,程 峰.能量耗散模型在大型地下洞群施工順序優(yōu)化分析中的應(yīng)用[J].巖石力學(xué)與工程學(xué)報(bào),2001,23(3):333-336. Zhu Weishen，Li Shucai，Cheng Feng.Application of energy dissipation model to optimization of construction order for large underground caverns[J].Chinese Journal of Rock Mechanics and Engineering，2001,23(3):333-336.

[5] 劉泉聲,王崇革.巖石時(shí)-溫等效原理的理論與實(shí)驗(yàn)研究[J].巖石力學(xué)與工程學(xué)報(bào),2002,21(2):193-198. Liu Quansheng，Wang Chongge.The theoretical and experimental study on time-temperature equivalent principle for rock[J].Chinese Journal of Rock Mechanics and Engineering，2002,21(2):193-198.

[6] 周科平.采場(chǎng)結(jié)構(gòu)參數(shù)的遺傳優(yōu)化[J].礦業(yè)研究與開(kāi)發(fā)，2000(3):7-10. Zhou Keping.Genetic optimization of stope structure parameters[J].Mining Research and Development,2000(3)：7-10.

[7] 楊朝暉,劉浩吾.地下工程圍巖穩(wěn)定性分類(lèi)的人工神經(jīng)網(wǎng)絡(luò)模型[J].四川聯(lián)合大學(xué)學(xué)報(bào)：工程科學(xué)版,1993,3(4):66-72. Yang Zhaohui，Liu Haowu.Artificial neural network model for the stability classification of adjoining rock of underground construction[J].Journal of Sichuan Union University：Engineering Science,1993,3(4)：66-72.

[8] 朱伯芳.有限單元法原理與應(yīng)用[M].北京:水利電力出版社,1979. Zhu Bofang.The Finite Element Method Theory and Applications[M].Beijing：Hydraulic Publisher，1979.

[9] 葉海旺,常 劍,周 磊.基于FLAC3D的采場(chǎng)結(jié)構(gòu)參數(shù)優(yōu)化[J].金屬礦山,2010(12)：6-9. Ye Haiwang，Chang Jian，Zhou Lei.Parameters optimization of stope structure based on FLAC3D[J].Metal Mine,2010(12):6-9.

[10] 朱志彬,劉成平.厚大礦體回采順序和采場(chǎng)結(jié)構(gòu)參數(shù)優(yōu)化[J].中國(guó)礦山工程,2008,37(4):20-24. Zhu Zhibin，Liu Chengping.Stoping sequence and optimum of structural parameter of thick orebody[J].China Mine Engineering，2008,37(4):20-24.

[11] 羅周全,管佳林,馮富康,等.盤(pán)區(qū)隔離礦柱采場(chǎng)結(jié)構(gòu)參數(shù)數(shù)值優(yōu)化[J].采礦與安全工程學(xué)報(bào),2012，29(2):261-264. Luo Zhouquan，Guan Jialin，F(xiàn)eng Fukang，et al.Stope structural parameters of panel isolation pillar numerical optimization[J].Journal of Mining and Safety Engineering，2012，29(2):261-264.

[12] 喬 蘭.新城金礦深部采場(chǎng)結(jié)構(gòu)參數(shù)和開(kāi)采順序優(yōu)化研究[J].金屬礦山,2001(6):11-15. Qiao Lan.Study on the optimization of deep stope structure parameters and mining sequence of Xincheng Gold Mine[J].Metal Mine,2001(6):11-15.

[13] 關(guān) 凱,朱萬(wàn)成,張洪訓(xùn),等.新城金礦中段頂柱盤(pán)區(qū)回采順序優(yōu)化[J].金屬礦山,2013(9):44-49. Guan Kai，Zhu Wancheng，Zhang Hongxun，et al.Optimization of panel mining sequence for the middle-section top pillar of Xincheng Gold Mine[J].Metal Mine，2013(9):44-49.

[14] 楊繼勝,余偉健.復(fù)雜地質(zhì)條件下礦體開(kāi)采順序的數(shù)值分析[J].金屬礦山,2012(7):13-18. Yang Jisheng，Yu Weijian.Numerical analysis on mining sequence under complex geological conditions[J].Metal Mine,2012(7):13-18.

(責(zé)任編輯 徐志宏)

Optimization of Structure Parameters in Deep & Large-Sized Orebody

Wu Jin1,2,3Li Yuanhui1,2Si Chengbin4Xu Shuai1,2

(1.Key Laboratory of Ministry of Education for Safe Mining，Shenyang 110819，China；2.School of Resources & Civil Engineering，Northeastern University，Shenyang 110819，China；3.School of Humanity and Social Science，Liaoning Radio and TV University，Shenyang 110161，China；4.Jiaojia Gold Mine，Shandong Gold Mining Co.，Ltd.，Yantai 261441，China)

Rational stope structural parameters may lead to a favorable mechanic state for the stopes and the stress and strain distribution of the wall rocks may become more homogenized.Under this condition，support workload is reduced and the mining intensity as well as the production efficiency can be improved on the premise of ensuring the stable exploring system and safe production.Based on the geological situation and mining methods of Sishanling iron ore，two-level orthogonal experimental design is made to decide the key five factors affecting the room exploring stability，i.e.room height，room width，stope length，pillar width and pillar filling method.Then，eight experimental schemes are concluded.The finite-difference software FLAC3Dis adopted to make numerical simulation on the stope structural parameters of different scheme at panel，and to obtain the feature of the stress，displacement，and plastic zone under different structural parameters.It is preliminarily obtained that the optimal structure parameters which can stabilize the stope is the height of stope 60 m，the length of stope 60 m，the width of room 18m，and the width of pillar 20 m.The calculating results show that the stope length of the roof influences the stress and roof displacement maximum，i.e.the longer the stope，the greater the stress.The compressive stress mainly concentrates on the panel pillars，while the tensile stress concentrates on the roof and the bottom.The results may provide technical support for safe and efficient mining of rooms and pillars within panel．

Orthogonal experimental design，Deep & large-sized orebody，Stage room mining method with subsequent filling，Optimization of Structural Parameters，Numerical simulation

2014-09-09

國(guó)家自然科學(xué)基金項(xiàng)目(編號(hào)：51204031，51274055)，“十二五”國(guó)家科技支撐計(jì)劃項(xiàng)目(編號(hào)：2013BAB02B03)，教育部項(xiàng)目基本科研業(yè)務(wù)費(fèi)項(xiàng)目(編號(hào)：130401007)，教育部項(xiàng)目博士點(diǎn)基金博導(dǎo)類(lèi)項(xiàng)目(編號(hào)：20130042110010)。

鄔 金(1982—)，女，博士研究生。通訊作者 李元輝(1968—)男，教授，博士研究生導(dǎo)師。

TD853.36

A

1001-1250(2014)-11-011-05