黃 敏

(1.長沙礦山研究院有限責(zé)任公司, 湖南長沙 410012；2.金屬礦山安全技術(shù)國家重點(diǎn)實(shí)驗(yàn)室, 湖南長沙 410012)

地下開采對山體邊坡影響的數(shù)值模擬分析

黃 敏1，2

(1.長沙礦山研究院有限責(zé)任公司, 湖南長沙 410012；2.金屬礦山安全技術(shù)國家重點(diǎn)實(shí)驗(yàn)室, 湖南長沙 410012)

采用FLAC3D模擬在重力作用下礦體開挖對山體邊坡的影響，研究表明，對山體邊坡影響較大的區(qū)域位于17?！?0?？碧骄€，對礦體影響較大的區(qū)域位于17?！?8?？碧骄€，這與礦山生產(chǎn)現(xiàn)狀基本吻合，表明在上部采場開采過程中應(yīng)注意山體邊坡的穩(wěn)定性。

地下開采；邊坡穩(wěn)定性；數(shù)值模擬

0 工程概況

國內(nèi)某銅礦礦體為多層、緩傾斜-傾斜、中厚-厚礦體,沿走向和傾向上礦體的產(chǎn)狀和厚度變化很大,且厚大部分多靠近地表。針對礦體的不同產(chǎn)狀,目前采礦方法主要為底盤漏斗空場法、房柱法和全面法等。該銅礦礦體為硫化礦,圍巖以大理巖為主。按照設(shè)計(jì)要求分為上、下2個采場進(jìn)行開采,3150標(biāo)高以上為上部采場,3150標(biāo)高以下為下部采場,采場勘探線剖面圖如圖1所示。經(jīng)過多年的開采,山體邊坡穩(wěn)定性受到破壞,曾出現(xiàn)過較大的滑坡災(zāi)害。為真實(shí)地反映礦山地下開采對山體邊坡所造成的影響,運(yùn)用三維數(shù)值模擬軟件進(jìn)行模擬分析。

通過建立符合實(shí)際的三維模型,并采用FLAC3D數(shù)值模擬軟件對開采過程進(jìn)行數(shù)值模擬研究,模擬主要目的是分析礦體回采過程中山體邊坡應(yīng)力、位移分布演化及塑性區(qū)的發(fā)展過程。

圖1 采場勘探線剖面

1 數(shù)值模型

在現(xiàn)場調(diào)查及提供的平、剖面圖的基礎(chǔ)上,建立數(shù)值計(jì)算模型,由于巖土工程問題的復(fù)雜性,數(shù)值模擬計(jì)算不可能考慮到所有因素,需要作一定簡化,根據(jù)圣維南原理以及相關(guān)理論分析可知:模型范圍的選擇只需要重點(diǎn)考慮工程地質(zhì)體區(qū)域輪廓線尺寸的2.5～3倍,開挖后產(chǎn)生的應(yīng)力和位移對此范圍之外的影響甚微,可忽略不計(jì)。最終建立好的FLAC3D單元體模型見圖2,模型共劃分不274590節(jié)點(diǎn), 207840個單元,x方向長為600m；y方向?yàn)殚L為750 m；z方向?yàn)?00 m,底部標(biāo)高為3000 m,頂部標(biāo)高為3400 m,建模范圍為0?？碧骄€至27?？碧骄€,模擬在重力作用及礦體開挖過程中對山體邊坡的影響。

圖2 FLAC3D單元體模型

根據(jù)實(shí)際開采情況,整個模擬計(jì)算過程分為3 步:原始應(yīng)力平衡階段,該步驟無任何開挖充填過程,只施加自重應(yīng)力場,先開挖下部采場,最后開挖上部采場。

計(jì)算區(qū)域邊界采取位移約束,即模型底部所有節(jié)點(diǎn)采用x、y、z 3個方向約束,模型x方向的兩端采用x方向約束,模型y方向的兩端采用y方向約束。模型坡面為自由邊界。約束后的模型見圖1。數(shù)值計(jì)算只考慮重力作用,選取重力加速度為9.8 m/s2,方向垂直向下,選用莫爾-庫侖(Mohr-Coulomb)模型,并賦材料力學(xué)參數(shù)(見表1)。

表1 材料力學(xué)參數(shù)

2 模擬計(jì)算與分析

通過模擬計(jì)算可知,在初始應(yīng)力平衡狀態(tài)下,模型中僅有自重應(yīng)力場作用,主應(yīng)力從上往下依次增大,主應(yīng)力線呈水平狀態(tài),最大主應(yīng)力-28.5 MPa,最小主應(yīng)力0.653 MPa,整個模型位移垂直向下,位移從上部到下部依次減小,最大位移從上部1.2 m到模型底部的1.5 cm。從位移矢量圖來看(見圖3),山體邊坡坡面位移豎直向下,可知山體邊坡處于穩(wěn)定狀態(tài)下。

圖3 原始狀態(tài)下13?？碧骄€上位移矢量

由圖4可以看出,下部采場開采以后,應(yīng)力變化不大。山體邊坡和礦體均發(fā)生了位移,最大位移8 cm,山體邊坡最大位移位于16?！?0?？碧骄€,礦體最大位移位于17?！?8?？碧骄€。下部采場礦體開挖后,塑性區(qū)和原始狀態(tài)相比,塑性區(qū)面積擴(kuò)大,山體邊坡和礦體均受到開挖擾動,其穩(wěn)定性均有所降低。從位移矢量圖上看(見圖5),位移方向指向凌空面,說明山體受開挖擾動,出現(xiàn)了滑坡趨勢,山體穩(wěn)定性降低。

由圖6可以看出,上部采場礦體開采以后,山體邊坡和礦體位移繼續(xù)增大,最大位移分別達(dá)到18 cm和14 cm,山體最大位移位于17＃～20?？碧骄€,礦體最大位移位于17＃～18?？碧骄€,與實(shí)際情況相符合。上部采場礦體開挖后,塑性區(qū)發(fā)生了變化,塑性區(qū)面積擴(kuò)大,山體和礦體均受到開挖擾動,其穩(wěn)定性均有所降低。從位移矢量圖上看(見圖6),位移方向指向凌空面,加之位移較大,山體出現(xiàn)了滑坡,這與實(shí)際開采情況基本一致。

圖4 下部采場開挖后山體最大位移分布

圖5 下部采場開挖后13?？碧骄€位移矢量

圖6 上部采場開挖后山體最大位移分布

3 結(jié) 論

采用FLAC3D數(shù)值模擬軟件對礦段地下開采現(xiàn)狀進(jìn)行了再現(xiàn)模擬研究,從應(yīng)力、位移和塑性區(qū)綜合分析來看,開采至今,對山體影響較大的區(qū)域位于17?！?0＃勘探線,對礦體影響較大的區(qū)域位于17?！?8?？碧骄€,這與礦山生產(chǎn)現(xiàn)狀基本吻合。因此,在上部采場開采過程中應(yīng)注意山體邊坡的穩(wěn)定性。

[1]鄧紅衛(wèi),朱和玲,周科平,等.基于FLAC3D數(shù)值模擬的前后處理優(yōu)化研究[J].礦業(yè)研究與開發(fā),2007,28(2):60-62.

[2]彭文斌.FLAC3D實(shí)用教程[M].北京:機(jī)械工業(yè)出版社,2008.

[3]陳育民,徐鼎平.FLAC/FLAC3D基礎(chǔ)與工程實(shí)例[M].北京:中國水利水電社出版社,2009.

[4]封凡忠,馬云龍.FLAC3D技術(shù)在邊坡穩(wěn)定性分析上的應(yīng)用[J].山西建筑,2010,36(26):97-98.

[5]李 博,楊志強(qiáng),高 謙.深凹露天轉(zhuǎn)地下開采對高陡邊坡穩(wěn)定性影響數(shù)值分析[J].礦業(yè)研究與開發(fā),2015,35(2):65-68.

[6]劉 波,韓彥輝.FLAC原理、實(shí)例及應(yīng)用指南[M].北京:人民交通出版社,2005.

[7]周偉永,饒運(yùn)章,汪 弘,等.基于FLAC3D的采場穩(wěn)固性數(shù)值模擬研究[J].礦業(yè)研究與開發(fā),2014,34(2):13-17.

2014-02-17)

黃 敏(1985-)，男，湖南長沙人，碩士，助理工程師，主要從事采礦技術(shù)及工程，Email:huangmin923923＠163. com。