葉樹學，陶 磊，胡京濤

(長沙有色冶金設(shè)計研究院有限公司， 湖南 長沙 410019)

采場跨度優(yōu)化數(shù)值模擬研究

葉樹學，陶 磊，胡京濤

(長沙有色冶金設(shè)計研究院有限公司， 湖南 長沙 410019)

采場跨度是影響采場結(jié)構(gòu)穩(wěn)定性的重要因素。為確定采場回采的最優(yōu)跨度，以礦山實際地質(zhì)條件為基礎(chǔ)，借助FLAC3D數(shù)值計算軟件，分別對采場跨度為6，8，10 m和12 m 4種情況進行數(shù)值模擬，并從最大主應力、垂直方向沉降及塑性區(qū)3方面對不同跨度情況進行對比分析，最終確定采場回采最優(yōu)跨度為8 m。

采礦跨度；數(shù)值模擬；采場沉降

1 礦山概況

某鋁土礦山某采場下分段位于北區(qū)1150分段和1160分段之間，采場標高確定巖體埋藏深度為220 m，礦體沿走向長50 m，沿傾向為70 m，礦體平均真厚度為4.43 m，傾角8°左右。礦塊呈似層狀單斜產(chǎn)出，總體傾向北西，傾向為330°～340°。采場周邊地質(zhì)構(gòu)造整體比較簡單，南部為F32斷層，北部為F19斷層的紅花寨礦塊。采場屬沉積型鋁土礦床，含礦巖系為下石炭統(tǒng)大塘組(C1d)。含礦巖系厚度具明顯二分性，上段為含鋁巖系，下段為含鐵巖系。上段含鋁巖系(C1d2)由粘土巖、黃鐵礦、鋁土巖、鋁土礦等組成，為鋁土礦礦體的產(chǎn)出層位，具工業(yè)價值的鋁土礦礦體產(chǎn)于含鋁巖系中。

為準確獲得鋁土礦床的物理力學性質(zhì)，在現(xiàn)場在礦體、上下盤圍巖等具有代表性的位置及非爆破開挖方式產(chǎn)生的大塊上進行取樣，以減少爆破損傷對巖石力學性質(zhì)的影響。并在實驗室進行一系列巖石力學物理力學試驗，確定礦床相關(guān)巖體的巖石物理力學參數(shù)(見表1)。

表1 試驗所得巖石物理力學參數(shù)

2 FLAC3D數(shù)值模擬方法

在采礦、水利和地下工程中，經(jīng)常會遇到巖石力學問題，這些問題包括確定工程開挖前后的應力和位移及其在開挖過程中力學參數(shù)的變化，以便為設(shè)計和施工提供可靠的依據(jù)。已有的巖石力學理論解只能解決圓形或橢圓形硐室等簡單工程形狀的問題，對于幾何形狀復雜的硐室或采場只能用數(shù)值方法來近似求解。

FLAC3D數(shù)值模擬軟件是一款連續(xù)介質(zhì)力學分析軟件，它的基本分析原理是源于流體力學的拉格朗日差分法。包含11種材料本構(gòu)摸型、5種計算模式和多種單元形式及邊界條件，可以模擬巖土體的多種施工工況，具有諸多優(yōu)點：

(1) 采用混合離散法，所有的應力、應變等變量均在網(wǎng)格劃分后的四面體上進行，而六面體單元上的變量可通過加權(quán)平均計算得到，采用此法來模擬土體材料的塑性破壞和塑性流動；

(2) 采用有限差分法，無論是非線性本構(gòu)關(guān)系還是線性本構(gòu)關(guān)系，其實質(zhì)上并無差別，只需根據(jù)已知的應變增量就可以很方便地求出應力增量和不平衡力，無需形成剛度矩陣，也無需推導大變形本構(gòu)模型；

(3) 在求解過程中通過增加阻尼力來使系統(tǒng)運動衰減至平衡狀態(tài)，采用此法在處理不穩(wěn)定問題時不會遇到數(shù)值上的障礙。采用FLAC3D軟件對巖體進行摸擬的計算流程見圖1。

3 采場跨度數(shù)值模擬結(jié)果及分析

3.1 最大主應力分析

4種不同跨度的采場剖面最大主應力分布云圖見圖2。通過對比分析，可得如下結(jié)論：

圖1 FLAC3D計算流程

圖2 各方案最大主應力分布云圖

(1) 采場的開挖會引起周圍區(qū)域的應力重分布，采場的頂板和相鄰礦房位置呈現(xiàn)出一定的應力集中現(xiàn)象。

(2) 當跨度為6 m、8 m時，頂板最大主應力值分別為0.088 MPa、0.125 MPa，二者均處于拉應力狀態(tài)，且應力集中范圍較小，說明采場處于穩(wěn)固狀態(tài)；當跨度增加到10 m時，最大主應力值為0.405 MPa，呈現(xiàn)為拉應力狀態(tài)，處于護頂?shù)V壁極限抗拉強度值附近，但直接頂板部分巖體破碎不穩(wěn)固，在該拉應力作用下可能發(fā)生破壞，因此頂板圍巖存在一定程度的失穩(wěn)危險；當跨度增加到12 m時，最大主應力值為0.9 MPa，呈現(xiàn)為拉應力狀態(tài)，且其值遠大于抗拉強度，說明采場頂板處一定發(fā)生剪切破壞或拉伸破壞。

(3) 隨著采場跨度的增加，相鄰礦房部分處于較大壓應力范圍內(nèi)，但該值沒有超過巖體的單軸抗壓強度，不足以造成采場圍巖的屈服破壞。

3.2 垂直方向沉降分析與比較

4種不同跨度的采場剖面垂直方向位移分布云圖見圖3。通過對比分析，可得如下結(jié)論：

圖3 各方案垂直方向位移云圖

(1) 采場開挖之后，頂?shù)装宥汲霈F(xiàn)不同程度的位移變化，其中頂板主要呈現(xiàn)出下沉現(xiàn)象，而底板則呈現(xiàn)出底鼓現(xiàn)象。

(2) 當跨度為6，8，10，12 m時，頂板沉降量最大值分別為0.98 cm、1.14 mm、1.3 cm、1.47 cm，其值較為接近，因此可以認為在單純開挖一個礦房時，跨度的增加對頂板的沉降量沒有太大的影響；從垂直沉降的影響范圍來看，跨度從10 m增加到12 m后，影響范圍有了很大的變化，因此可以認為12 m的跨度將引起大范圍的沉降，對礦體開采不利。

(3) 因為該礦埋深較淺，地應力較小，底板的底鼓現(xiàn)象并不嚴重，在實際工作中該影響可忽略不計。

3.3 塑性區(qū)分析與比較

4種不同跨度采場的剖面塑性區(qū)分布云圖見圖4。通過對比分析，可得如下結(jié)論：

圖4 各方案塑性區(qū)分布云圖

(1) 采場開挖后，采場頂?shù)装寮皟蓭途霈F(xiàn)不同程度的塑性區(qū)，但隨著開挖跨度的增加，兩幫塑性區(qū)范圍幾乎沒有改變，而采場頂板的塑性區(qū)范圍則有所擴大，說明隨著采場跨度的增加，頂板容易發(fā)生剪切破壞，需要加強頂板的支護工作。

(2) 當采場跨度為6 m、8 m時，護頂?shù)V壁和直接頂板塑性區(qū)表現(xiàn)為shear-p、tension-p，即為彈性，小范圍區(qū)域內(nèi)出現(xiàn)剪切破壞區(qū)域(shear-n)，因此可以認為該范圍內(nèi)巖體在荷載作用下仍處于彈性變形階段，但開挖區(qū)域周邊圍巖經(jīng)歷過應力峰值，說明頂板仍處在較為穩(wěn)固的狀態(tài)，局部破碎區(qū)域需要加強支護；當采場跨度為10 m、12 m時，護頂?shù)V壁和直接頂板出現(xiàn)剪切破壞的區(qū)域(shear-n)，且直接頂板與間接頂板的接觸面上出現(xiàn)復合破壞區(qū)域(shear-n、tension-n)，此時采場已經(jīng)處于失穩(wěn)狀態(tài)。

4 結(jié) 論

FLAC3D軟件針對采場跨度分別為6，8，10 m和12 m 4種情況，分別進行數(shù)值模擬分析，從最大主應力方面考慮，當跨度增加到10 m時，礦壁圍巖存在一定程度的失穩(wěn)危險；從垂直方向沉降方面考慮，當跨度增加到12 m時，將引起大范圍沉降；從塑性區(qū)分布方面考慮，當跨度增加至10 m時，圍巖出現(xiàn)大范圍剪切破壞及部分復合破壞，采場易失穩(wěn)；綜合可得，采場極限跨度為8～10 m區(qū)間內(nèi)，為保障一定的安全邊際，確定采場跨度的最優(yōu)值為8 m。

[1]李啟月,劉 愷,李夕兵.基于協(xié)同回采的深部厚大礦體分段充填采礦法[J].工程科學學報,2016(11):1515-1521.

[2]牛 勇,李克鋼.房柱采礦法采場結(jié)構(gòu)參數(shù)優(yōu)化理論研究及應用[J].黃金科學技術(shù),2015(06):75-80.

[3]許振華,喬 蘭.深井硬巖礦山采場結(jié)構(gòu)參數(shù)優(yōu)化的數(shù)值模擬研究[J].中國礦業(yè),2014(03):84-87.

[4]劉 冬,金長宇,夏自鋒,等.焦家金礦破碎礦體采場進路的跨度優(yōu)化研究[J].金屬礦山,2014(02):18-21.

[5]董金奎,馮夏庭,張希巍,等.地下采場破碎巖體穩(wěn)定性評價與參數(shù)優(yōu)化[J].東北大學學報(自然科學版),2013(09):1322-1326.

[6]黃英華,焦文宇,唐紹輝,等.充填體下隔離中段采場結(jié)構(gòu)參數(shù)優(yōu)化[J].湖南科技大學學報(自然科學版),2012(04):11-17.

[7]石 飛.采場跨度優(yōu)化有限元模擬研究[J].有色冶金設(shè)計與研究,2012(05):6-9.

[8]申 延.寺莊礦區(qū)工程巖體穩(wěn)定性及采場合理跨度研究[D].長沙：中南大學,2012.

[9]彭 康,李夕兵,彭述權(quán),等.基于響應面法的海下框架式采場結(jié)構(gòu)優(yōu)化選擇[J].中南大學學報(自然科學版),2011(08):2417-2422.

[10]劉洋樹,李安平,王 剛,等.VCR法采場結(jié)構(gòu)參數(shù)優(yōu)化的相似模型實驗[J].有色礦冶,2011(02):10-15.

[11]葉海旺,常 劍,周 磊.基于FLAC(3D)的采場結(jié)構(gòu)參數(shù)優(yōu)化[J].金屬礦山,2010(12):6-8,26.

[12]李夕兵,劉志祥,彭 康,等.金屬礦濱海基巖開采巖石力學理論與實踐[J].巖石力學與工程學報,2010(10):1945-1953.

[13]劉培正,張傳信,胡永泉,等.厚大礦體分段空場嗣后充填采場結(jié)構(gòu)參數(shù)優(yōu)化研究[J].金屬礦山,2009(11):10-13,123.

[14]薛美英.超大跨度采場結(jié)構(gòu)參數(shù)優(yōu)化研究[J].礦業(yè)研究與開發(fā),2016,36(09):14-18.

[15]張海波,李示波,張 揚,等.金屬礦山嗣后充填采場頂板合理跨度參數(shù)研究及建議[J].金屬礦山,2014(06):21-24.

[16]翟會超,閆滿志.壁式削壁充填法采場頂板臨界跨度計算[J].巖土力學,2013(S1):337-339.

[17]姚高輝,吳愛祥,王洪江.白牛廠礦節(jié)理化巖體力學參數(shù)與采場安全跨度評價[J].礦業(yè)研究與開發(fā),2010,30(04):59-62.

2017-03-29)