FLAC3D數(shù)值模擬露天轉(zhuǎn)地下邊坡及采場穩(wěn)定性研究

李啟航，李小雙，2，3，耿加波，羅 浪

(1.江西理工大學(xué) 資源與環(huán)境工程學(xué)院，江西 贛州 341000；2.中鋼集團馬鞍山礦山研究總院股份有限公司，安徽 馬鞍山 243000；3.金屬礦山安全與健康國家重點實驗室，安徽 馬鞍山 243000)

根據(jù)前人運用數(shù)值模擬方法在礦山巖體力學(xué)中的研究和探索[1-7]，本文擬利用FLAC3D數(shù)值模擬分析方法，對緩傾斜薄至中厚磷礦床在單一地下開采及露天轉(zhuǎn)地下開采兩種不同開挖方式下坡高對采場頂板、圍巖的穩(wěn)定性及礦壓活動規(guī)律的影響。對未來的磷礦山地下開采工作具有一定的技術(shù)參考價值。

1 FLAC3D軟件簡介

FLAC3D是在FLAC基礎(chǔ)上研發(fā)的能夠進行土質(zhì)、巖石和其他材料的三維結(jié)構(gòu)受力特性模擬和塑性流動分析的軟件[8]。FLAC3D具有比較完善的后處理功能，對于數(shù)據(jù)和結(jié)果的處理比較容易實現(xiàn)。近年來，國內(nèi)外諸多科研工作者利用FLAC3D軟件對采礦過程中圍巖及邊坡穩(wěn)定性、采場結(jié)構(gòu)參數(shù)優(yōu)化、巷道與工作面的時空關(guān)系等一系列復(fù)雜力學(xué)問題展開了大量的研究，并取得了許多重要的研究成果[9]。

2 模型的構(gòu)建及計算方案

2.1 基本假設(shè)

本文的數(shù)值模擬主要對應(yīng)力與變形進行分析研究。

1)巖性的假設(shè)

假設(shè)巖石有各向同性、均質(zhì)的特性，應(yīng)力與應(yīng)變符合摩爾-庫倫準則。

2)水平構(gòu)造應(yīng)力的假設(shè)

開挖區(qū)域礦體構(gòu)造簡單，假設(shè)水平構(gòu)造應(yīng)力不予考慮。

3)礦房結(jié)構(gòu)的簡化

將數(shù)值模型中的巷道、天井、斜井、聯(lián)絡(luò)巷以及溜井等礦房結(jié)構(gòu)簡化為實體。

2.2 模型幾何尺寸

以云南晉寧周邊典型的緩傾斜薄至中厚磷礦床為研究背景，以層狀巖質(zhì)邊坡與地下采場圍巖及其上覆巖體組成的復(fù)合采動系統(tǒng)為研究對象。兩組邊坡模型尺寸分別為：200 m邊坡模型尺寸水平(X)方向長度為600 m，走向(Y)方向長度為300 m，垂直(Z)方向為350 m，整個模型共包含544 833個節(jié)點和526 550個單元；300 m邊坡模型尺寸水平(X)方向長度為600 m，走向(Y)方向長度為300 m，垂直(Z)方向為450 m，整個模型共包含623 934個節(jié)點和603 550個單元。在模型計算時簡化巖體應(yīng)力場，不考慮其他構(gòu)造應(yīng)力的作用，只考慮自身重力影響，因此只對模型施加-9.8 m/s2的反向應(yīng)力加速度。兩種方案的單一地下開采及露天轉(zhuǎn)地下開采三維模型如圖1所示。

圖1 不同邊坡高度下單一地下開采與露天轉(zhuǎn)地下開采三維模型圖Fig.1 3D model diagram of single underground mining and open-pit conversion to underground mining under different slope heights

2.3 模型計算本構(gòu)關(guān)系

本次數(shù)值模擬采用彈塑性本構(gòu)模型，F(xiàn)LAC3D巖石材料破壞準則采用摩爾-庫倫屈服準則來描述巖體的強度特征。

2.4 模型巖體參數(shù)

根據(jù)實地調(diào)研和相關(guān)巖石力學(xué)試驗研究結(jié)果，數(shù)值模擬計算所采用的磷礦巖體力學(xué)參數(shù)如表1所示。

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

2.5 礦山開采方式

根據(jù)礦山實際情況，選取境界礦柱厚度24m，由于整個地下開采礦體傾向長度為348 m，垂直開采深度為84.2 m。因此設(shè)置3個開采階段，每個開采階段包含7個礦房和6個礦柱，單一開采階段沿礦層傾向長度為106 m，垂向深度為25 m，階段與階段之間設(shè)置15 m的階段礦柱。同時選取每步開挖的礦房10 m，礦柱6 m，且礦柱為連續(xù)型礦壁。另外，考慮到模型的邊界效應(yīng)，沿模型開挖方向(模型Y方向)兩側(cè)各留50 m邊界影響區(qū)域。開采順序為分階段開采，先沿礦層傾向逐一開采每個階段的礦房，待礦房開采完畢后再沿傾向?qū)ΦV柱進行逐一回采，回采完畢后即進入下一階段礦體的開挖。采場俯視圖如圖2所示。

圖2 采場俯視圖Fig.2 Top view of stope

2.6 數(shù)值模擬計算方案

本文展開了兩種模擬方案：方案1)先進行單一的露天開挖，待露天開挖結(jié)束后再轉(zhuǎn)入地下開采；方案2)在同樣地質(zhì)地貌條件下，直接進行地下開采。利用FLAC3D數(shù)值模擬軟件分析計算得到不同開采階段的應(yīng)力分布情況及數(shù)值大小來揭示露天終了邊坡在坡角一定的條件下，不同邊坡高度對地下采場的影響以及邊坡與地下開采耦合作用下覆巖采動響應(yīng)特征。

3 計算結(jié)果及分析

3.1 露天邊坡安全系數(shù)

為了確定在地下開采前邊坡是否安全穩(wěn)定需要對礦山露天終了邊坡的安全系數(shù)進行求解[10]。求解結(jié)果如表2所示。

表2 邊坡安全系數(shù)

邊坡穩(wěn)定性安全系數(shù)是指沿指定滑裂面的抗滑力與滑動力的比值，比值大于1時坡體穩(wěn)定；等于1時坡體處于極限平衡狀態(tài)；小于1時坡體發(fā)生破壞。由表2可知，兩種不同坡高的邊坡均整體穩(wěn)定，但隨著邊坡高度增加，邊坡發(fā)生滑移破壞的可能性增大。在不同開采階段下，根據(jù)FLAC3D數(shù)值模擬最大垂直應(yīng)力的結(jié)果可知：在本文模擬實驗環(huán)境中不考慮滑坡風(fēng)險。

3.2 坡高200 m模型分析

由圖3可知，第一階段礦體開挖完畢后，由于受到階段礦體的開挖卸荷，兩種開采方案的地下采場最大垂直應(yīng)力均出現(xiàn)在采空區(qū)前后兩端，此時采區(qū)兩側(cè)及露天轉(zhuǎn)地下開采形成的露天終了邊坡坡腳處形成了局部的應(yīng)力集中區(qū)域。第二階段礦體開挖完畢后，單一地下開采最大垂直應(yīng)力僅出現(xiàn)在兩個采空區(qū)的階段礦柱上。第三階段礦體開挖完畢后，其應(yīng)力相較第二階段再一次釋放，應(yīng)力降低區(qū)域的范圍再次擴大。在不同開采階段下，根據(jù)FLAC3D數(shù)值模擬最大垂直應(yīng)力的結(jié)果可知：采空區(qū)兩側(cè)第一階段到第三階段單一地下開采最大垂直應(yīng)力由20.9 MPa增加到28.8 MPa再增加到40.0 MPa；露天轉(zhuǎn)地下開采最大垂直應(yīng)力由19.4 MPa增加到26.5 MPa再增加到39.3 MPa；采空區(qū)頂板第一階段到第三階段單一地下開采最大垂直應(yīng)力由4.0 MPa增加到4.3 MPa再增加到4.7 MPa；露天轉(zhuǎn)地下開采最大垂直應(yīng)力由4.2 MPa減少到3.9 MPa再增加到4.6 MPa；通過最大垂直應(yīng)力數(shù)值變化特點，三個不同開采階段下單一地下開采最大垂直應(yīng)力比露天轉(zhuǎn)地下開采都要大。從第一階段到第三階段開挖，隨著礦體開挖埋深的增加，采空區(qū)頂板巖體向下沉降，導(dǎo)致間柱頂部剪切塑性區(qū)范圍擴大，但未到達邊坡臨空面，此時地下采場及整個巖體相對穩(wěn)定。

圖3 單一地下開采(圖左)與露天轉(zhuǎn)地下開采(圖右)各階段開挖后垂直應(yīng)力分布特征Fig.3 The vertical stress distribution features after each stage of excavation in single underground mining(picture left)and open-pit to underground mining(picture right)

3.3 坡高300 m模型分析

由圖4可知，第一階段礦體開挖完后，單一地下開采最大垂直應(yīng)力分布在采空區(qū)前后兩端，由于接近露天邊坡的臨空面，露天轉(zhuǎn)地下開采的最大垂直應(yīng)力只出現(xiàn)在采空場的前端，第二階段礦體開挖完畢后，單一地下開采最大垂直應(yīng)力僅出現(xiàn)在兩個采空區(qū)的階段礦柱右下側(cè)，露天轉(zhuǎn)地下開采最大垂直應(yīng)力出現(xiàn)在二階段采場前后兩端。第三階段礦體開挖完畢后，單一地下開采最大垂直應(yīng)力出現(xiàn)在兩個間柱右下側(cè)部位，露天轉(zhuǎn)地下開采最大垂直應(yīng)力僅出現(xiàn)在第二個階段礦柱的右下側(cè)。在不同開采階段下，根據(jù)FLAC3D數(shù)值模擬最大垂直應(yīng)力的結(jié)果可知：采空區(qū)兩側(cè)第一階段到第三階段單一地下開采最大垂直應(yīng)力由24.8 MPa增加到35.7 MPa再增加到44.5 MPa；露天轉(zhuǎn)地下開采最大垂直應(yīng)力由20.3 MPa增加到26.5 MPa再增加到31.8 MPa；采空區(qū)頂板第一階段到第三階段單一地下開采最大垂直應(yīng)力由3.9 MPa增加到4.0 MPa再增加到4.7 MPa；露天轉(zhuǎn)地下開采最大垂直應(yīng)力由4.1 MPa減少到3.8 MPa再增加到4.8 MPa。

圖4 單一地下開采(圖左)與露天轉(zhuǎn)地下開采(圖右)各階段開挖后垂直應(yīng)力分布特征Fig.4 The vertical stress distribution features after each stage of excavation in single underground mining(picture left)and open-pit to underground mining(picture right)

隨著工作面的開挖推進，兩種開采方案下的最大垂直應(yīng)力也跟著動態(tài)前移，當(dāng)露天開采完成后，露天邊坡坡腳附近出現(xiàn)應(yīng)力集中區(qū)域，單一地下開采在三個不同開采階段最大垂直應(yīng)力都比露天轉(zhuǎn)地下開采大。從第一階段到第三階段采場塑性區(qū)已經(jīng)發(fā)育充分，貫通至邊坡臨空面，整個地下采場及邊坡體已處于失穩(wěn)狀態(tài)，極易造成整個采場及邊坡失穩(wěn)垮塌。

4 結(jié)論

當(dāng)?shù)叵碌V體開挖不斷推進時，采場頂板最大垂直應(yīng)力點也沿著礦層傾向而動態(tài)前移。單一地下開采的頂板垂直應(yīng)力相較于露天轉(zhuǎn)地下開采更大。露天轉(zhuǎn)地下開采的終了邊坡越高，邊坡對采場覆層的影響效應(yīng)越明顯。露天轉(zhuǎn)地下開采，坡高為200 m時，巖體比較穩(wěn)定，坡高為300 m時，第三階段采空區(qū)覆巖體已處于失穩(wěn)狀態(tài)極易發(fā)生邊坡垮塌現(xiàn)象。因此在相似環(huán)境下的磷礦山開采工作中，建議控制露天邊坡高度在200 m以內(nèi)。