孫 健, 陳小偉, 宮長亮, 安 晶, 任鵬召, 張 磊, 張少杰

(1.鞍鋼集團礦業(yè)有限公司大孤山分公司, 遼寧 鞍山 114046; 2.中國恩菲工程技術(shù)有限公司, 北京 100038;3.北京科技大學(xué), 北京 100083)

1 前言

目前,我國淺部礦產(chǎn)資源逐漸減少,大部分露天礦山相繼進入地下開采階段[1-2]。露天轉(zhuǎn)地下開采引起地壓活動規(guī)律改變,使開采盤區(qū)圍巖力學(xué)條惡化,威脅礦山企業(yè)的生產(chǎn)安全[3-5]。為有效控制盤區(qū)圍巖穩(wěn)定,設(shè)計中通常留設(shè)承載礦柱,相關(guān)研究主要涉及礦柱承載機理及厚度優(yōu)化、礦柱尺寸設(shè)計、盤區(qū)采場跨度、隔離礦柱回采方案優(yōu)化、礦柱承載與破壞特征等方面內(nèi)容。王文軍等[6]研究了深部開采中隔離庫礦柱的承載機理,采用理論分析和數(shù)值模擬的方法,確定原巖隔離礦柱為25 m時滿足安全要求;賴祖豪等[7]基于量綱分析理論,選取5因素水平建立境界礦柱厚度計算模型,采用計算模型得到的境界礦柱厚度與實際值接近,該模型具有較強的實用性,可為露天轉(zhuǎn)地下境界礦柱厚度的確定提供理論依據(jù);盧萍[8]等采用FLAC3D數(shù)值模擬方法,進行了不同厚度隔離礦柱的數(shù)值模擬,確定礦山隔離礦柱厚度20 m為最優(yōu)值;范文濤等[9]基根據(jù)不同礦柱長寬比模型計算,發(fā)現(xiàn)礦柱的長寬比值越大,承載能力降低,且脆性破壞特征明顯;史秀志等[10]建立了礦柱簡化力學(xué)模型,利用能量突變分析法,導(dǎo)出盤區(qū)礦柱的失穩(wěn)判據(jù),據(jù)此得到盤區(qū)礦柱的安全寬度為31 m;張少杰等[11]研究了深井超大面積開采盤區(qū)條形礦柱的應(yīng)力和塑性區(qū)分布特征,揭示了盤區(qū)礦柱的承載作用;羅斌玉等[12]建立了壓剪荷載作用下礦柱力學(xué)模型,分析不同應(yīng)力路徑對礦柱穩(wěn)定性的影響;隔離礦柱回采方案與采場結(jié)構(gòu)參數(shù)優(yōu)化[13-14],主要根據(jù)盤區(qū)回采順序和盤區(qū)隔離礦柱采場結(jié)構(gòu)參數(shù)進行圍巖變形和應(yīng)力分布特征對比研究,選出最優(yōu)方案;朱萬成等[15]認為礦柱的承載及失穩(wěn)破壞與采空區(qū)穩(wěn)定性緊密相關(guān),歸納總結(jié)單礦柱穩(wěn)定性研究方法和多礦柱載荷傳遞行為及失穩(wěn)破壞機制。采區(qū)通常為多礦柱布置,礦柱除承受上覆巖層荷載外,多礦柱之間也存在傳遞荷載,荷載形式極為復(fù)雜,相關(guān)研究相對較少,導(dǎo)致工程實踐中礦柱尺寸設(shè)計缺少科學(xué)依據(jù)。

本文對某鐵礦露天轉(zhuǎn)地下開采盤區(qū)礦柱進行了數(shù)值模擬分析,在保證地下礦體開采量最大化的基礎(chǔ)上,研究不同礦柱寬度對圍巖穩(wěn)定性的影響,研究成果可為露天轉(zhuǎn)地下大規(guī)模開采條件下的礦柱設(shè)計提供依據(jù)。

2 工程概況

2.1 工程地質(zhì)

礦區(qū)斷裂構(gòu)造對鐵礦床的影響最大,斷裂基本沿鐵礦帶近于平行展布,走向280°～310°,由數(shù)條平行或相互交匯復(fù)合的斷層組成,總體構(gòu)造性質(zhì)為南盤向東走滑的以走滑為主的斷裂構(gòu)造。巖石風(fēng)化程度一般,發(fā)育深度30～80 m,鉆孔揭露風(fēng)化裂隙發(fā)育最大深度為62.37 m。礦區(qū)內(nèi)斷裂構(gòu)造較為發(fā)育,個別位置受斷裂構(gòu)造、層間裂隙或巖脈侵入、穿插致使巖性頻繁變化的影響,使巖石質(zhì)量和巖體完整性受到一定的破壞,破碎位置多位于綠泥石英片巖和石英巖中,使其巖石質(zhì)量及巖體完整性均較差。露天采場斷層多與邊幫呈大角度相交,斷層傾角較陡,或與邊幫坡面反傾,邊幫巖體受斷層及次生結(jié)構(gòu)面的影響較大。

礦體及圍巖完整性多為中等-完整,巖體質(zhì)量中等-極好,局部構(gòu)造帶及巖脈侵入接觸帶附近巖體較破碎,完整性和質(zhì)量較差。礦區(qū)工程地質(zhì)條件復(fù)雜類型為中等類型。

2.2 開采現(xiàn)狀

某鐵礦經(jīng)過多年露天開采后,現(xiàn)形成東西長1 700 m、南北寬1 520 m、垂直深達400 m橢圓形露天坑,其地形地貌景觀是整個礦坑呈倒金字塔形?，F(xiàn)露天采場最低標高為-354 m,露天采場上口地面標高約70～110 m;礦體長約1 200 m,走向310°～315°,傾角60°～75°,平均厚度128 m,設(shè)計坑底標高-414 m。根據(jù)礦山產(chǎn)量安排,露天開采預(yù)計2024年結(jié)束,后續(xù)需對礦山轉(zhuǎn)入地下開采進行規(guī)劃、建設(shè)。

根據(jù)礦體的開采技術(shù)條件,采礦方法主要采用大直徑深孔空場嗣后充填法和分段空場嗣后充填法。將礦體劃分為盤區(qū),以盤區(qū)為回采單元組織生產(chǎn)。采場分為礦房和礦柱,進行兩步驟回采,礦房、礦柱寬度均定為18 m,采場長為90 m,中段高度60 m,采用隔三采一的方式開采,一步礦房回采結(jié)束后,再對二步礦柱進行回采。

3 數(shù)值模擬

3.1 航拍模型

采用測繪專業(yè)級設(shè)備無人機(經(jīng)緯M300_RTK)在礦區(qū)進行了三維實景復(fù)制,航測范圍周長約15 km,面積約16 km2,精細測量的重點區(qū)域主要有:礦坑、排巖場和尾礦庫。

無人機傾斜攝影技術(shù)的工作流程分為像控點布設(shè)、測線規(guī)劃、數(shù)據(jù)采集、數(shù)據(jù)處理、成果表達等流程。經(jīng)過無人機航拍得到的影像數(shù)據(jù),通過影像數(shù)據(jù)中的角點、邊緣等灰度值差異較大的點等一系列特征點進行匹配,匹配完成后對數(shù)據(jù)進行矯正,以保證得到數(shù)據(jù)的精度。對擁有同一特征點的三張及以上照片進行三角運算,獲得特征點的空間數(shù)據(jù),并對空間數(shù)據(jù)進行誤差運算,矯正照片的位置信息和姿態(tài)信息。矯正完成后對照片進行逐像素匹配運算,進而生成稠密點云,將點云中的三個點連接生成三角網(wǎng),形成初步模型,然后對模型進行平滑處理得到最終模型,將模型對應(yīng)的照片的紋理信息映射到模型上,生成三維實景模型,具體如圖1所示。

圖1 三維實景模型

3.2 FLAC3D模型

1)FLAC3D軟件介紹

FLAC3D是用于工程計算的大型三維顯式拉格朗日有限差分程序,能夠模擬計算巖體材料受力狀態(tài)下的各種力學(xué)行為。本次研究過程中,根據(jù)材料力學(xué)特征,圍巖和礦體均采用摩爾庫侖屈服準則,即

(1)

式中:fs——材料發(fā)生破壞的臨界值,MPa;

σ1、σ3——最大和最小主應(yīng)力,MPa;

c——黏結(jié)力,MPa;

φ——摩擦角,MPa。

當fs>0時,材料將發(fā)生剪切破壞。在通常應(yīng)力狀態(tài)下,巖體的抗拉強度很低,因此可根據(jù)抗拉強度準則(σ3≥σt(抗拉強度,MPa))判斷巖體是否產(chǎn)生拉破壞。

2)模型建立

在礦區(qū)工程地質(zhì)調(diào)查、無人機實景復(fù)制的基礎(chǔ)上,綜合礦區(qū)的主要地質(zhì)構(gòu)造和工程巖體分組,采用大型有限差分軟件FLAC3D,構(gòu)建了該礦山的真三維數(shù)值模型,對露天轉(zhuǎn)地下開采全過程進行模擬計算。模型尺寸為3 800 m×3 800 m×900 m,具體如圖2所示。

圖2 三維數(shù)值模型

數(shù)值模型中的應(yīng)力邊界條件按照地應(yīng)力測量得到的地應(yīng)力分布規(guī)律(式(2)～式(4))施加,模型側(cè)面限制水平移動,底部限制垂直移動。模型中礦體與礦柱之間的位置關(guān)系如圖3所示。

圖3 模型垂直Y軸切面示意圖

根據(jù)礦山水壓致裂地應(yīng)力測量研究的數(shù)據(jù)顯示:礦區(qū)應(yīng)力場以水平應(yīng)力為主導(dǎo),最大水平主應(yīng)力與垂直應(yīng)力的平均比值為2.42,最小水平主應(yīng)力與垂直應(yīng)力的平均比值為1.49。平均最大主應(yīng)力方向為N62.6°E。礦區(qū)測點的三個主應(yīng)力隨深度變化的線性回歸方程為

SH=0.051 22H+5.158 21

(2)

Sh=0.031 04H+3.342 73

(3)

SV=0.028 19H

(4)

式中:SH——最大水平主應(yīng)力,MPa;

Sh——最小水平主應(yīng)力,MPa;

SV——垂直應(yīng)力,MPa

H——垂直深度,m。

3.3 力學(xué)參數(shù)

本研究計算所需的參數(shù)是在室內(nèi)巖石力學(xué)試驗的基礎(chǔ)上,并通過Hoek-Brown強度準則折減之后獲得,能夠真實反映現(xiàn)場巖體力學(xué)性質(zhì),具體參數(shù)取值見表1。

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

3.4 設(shè)計方案與模擬步驟

1)開采順序

設(shè)計對-324～-674 m之間的礦體進行回采,以-494 m水平為界分為上部采區(qū)和下部采區(qū),上部采區(qū)礦體分為-380 m、-430 m、-490 m三個無軌開采中段,下部采區(qū)礦體分為-550 m、-610 m、-670 m三個無軌開采中段。上部和下部兩個采區(qū)同時生產(chǎn),各采區(qū)均按從下往上開采順序進行開采。首采中段為-490 m中段和-670 m中段。

2)模型方案

盤區(qū)沿走向布置,長度120 m,中段高度60 m,盤區(qū)內(nèi)垂直走向劃分采場,5個采場組成一個盤區(qū)。為了得到最優(yōu)的盤區(qū)礦柱設(shè)計尺寸,將盤區(qū)礦柱的寬度設(shè)計為20 m、30 m、40 m以及34 m與26 m交替布置的四種情況,分別進行模擬計算。在模型中布置垂直Y軸的剖面A-A′和垂直Z軸的剖面B-B′,具體如圖4所示。擬分析的A-A′剖面布置在穿過深部礦體正中央的位置;擬分析的B-B′剖面布置在-494 m水平。

圖4 模型剖面布置

3)模擬步驟

采礦工程的力學(xué)特點是巖體力學(xué)行為與開采歷史和開采過程有關(guān),結(jié)合礦山地下開采規(guī)劃方案,確定數(shù)值模擬的計算過程。

(1)計算在給定邊界力學(xué)與位移邊界條件下模型的初始狀態(tài)。

(2)模擬露天坑的形成過程,模擬掛幫礦的開采過程。

(3)在盤區(qū)礦柱分別為20 m、30 m、40 m以及34 m和26 m間隔布置四種方案下,模擬地下開采過程,分析研究地下開采過程中預(yù)留礦柱及周邊巖體的穩(wěn)定性特征。

4 礦柱穩(wěn)定性分析

4.1 數(shù)值模擬結(jié)果

為研究盤區(qū)礦柱寬度對其穩(wěn)定性的影響,根據(jù)數(shù)值模擬結(jié)果,得到三種不同盤區(qū)礦柱寬度(20 m、30 m和40 m)方案下地下開采后Y軸的剖面A-A′和水平剖面B-B′的塑性區(qū)分布及最大主應(yīng)力分布情況,具體如圖5至圖8所示。

圖5 地下開采后A-A′剖面的塑性區(qū)分布

圖6 地下開采后B-B′剖面的塑性區(qū)分布

圖7 地下開采后A-A′剖面最大主應(yīng)力分布

圖8 地下開采后B-B′剖面最大主應(yīng)力分布

數(shù)值模擬結(jié)果顯示,當盤區(qū)礦柱寬為20 m時,地下開采結(jié)束后,礦柱處于高應(yīng)力狀態(tài),最大主應(yīng)力值約為23 MPa,處于塑性破壞狀態(tài);當盤區(qū)礦柱寬為30 m時,地下開采結(jié)束后,礦柱處于高應(yīng)力狀態(tài),但未進入塑性狀態(tài),盤區(qū)礦柱的最大主應(yīng)力值約為19 MPa,礦柱較好的保護了開采區(qū)域的穩(wěn)定性;當盤區(qū)礦柱寬為40 m時,地下開采結(jié)束后,礦柱處出現(xiàn)應(yīng)力集中,但未進入塑性狀態(tài),盤區(qū)礦柱的最大主應(yīng)力值約為16 MPa,礦柱保護了開采區(qū)域的穩(wěn)定性,但經(jīng)濟性相對較差。

綜合考慮三種盤區(qū)礦柱寬度方案下的模擬結(jié)果,認為盤區(qū)礦柱寬度為30 m的方案優(yōu)于盤區(qū)礦柱寬度為20 m和40 m的方案。

4.2 優(yōu)化設(shè)計

為了最優(yōu)化盤區(qū)礦柱的尺寸,盤區(qū)礦柱寬度設(shè)置為34 m與26 m交替變化方案下地下開采后的塑性區(qū)分布及最大主應(yīng)力分布情況如圖9至圖12 所示。

圖9 地下開采后A-A′剖面的塑性區(qū)分布

圖10 地下開采后B-B′剖面的塑性區(qū)分布

圖11 地下開采后A-A′剖面最大主應(yīng)力分布

圖12 地下開采后B-B′剖面最大主應(yīng)力分布

當盤區(qū)礦柱寬為34 m與26 m交替變化時,地下開采結(jié)束后,礦柱處于高應(yīng)力狀態(tài),但未進入塑性狀態(tài),盤區(qū)礦柱的最大主應(yīng)力值約為18 MPa,礦柱較好的保護了開采區(qū)域的穩(wěn)定性。相較于盤區(qū)礦柱均設(shè)置為30 m,礦柱的最大主應(yīng)力值基本不變。

綜合考慮上述幾種盤區(qū)礦柱寬度方案下的模擬結(jié)果,認為盤區(qū)礦柱寬度為34 m與26 m交替的方案最優(yōu)。

5 結(jié)論

本文采用數(shù)值模擬方法對某露天礦山盤區(qū)礦柱的塑性分布和應(yīng)力場特征進行了分析,得出以下結(jié)論。

(1)地下開采盤區(qū)礦柱具有良好的承載作用,采用數(shù)值模擬方法,能夠反映多礦柱開采工況時上覆巖層及礦柱自身的穩(wěn)定情況,是地下開采礦柱設(shè)計的有利工具。

(2)盤區(qū)礦柱寬較小時,礦柱處于相對較高應(yīng)力和塑性破壞狀態(tài);隨著礦柱寬度的增大,塑性破壞分布區(qū)域減少,最大主應(yīng)力值降低,礦柱較好的保護了開采區(qū)域的穩(wěn)定性;從經(jīng)濟性方面考慮,認為盤區(qū)礦柱寬度為30 m的方案優(yōu)于盤區(qū)礦柱寬度為20 m和40 m的方案。

(3)采用盤區(qū)礦柱寬度34 m與26 m交替變化方案,礦柱處于高應(yīng)力狀態(tài),但未進入塑性狀態(tài),礦柱能夠較好保護開采區(qū)域的穩(wěn)定性,為最優(yōu)方案。