陳銀萍 金愛(ài)兵 王 凱 馬 松 姚鵬飛 吳佐漢

(1.北京科技大學(xué)土木與環(huán)境工程學(xué)院，北京 100083；2.四川省安全科學(xué)技術(shù)研究院,四川 成都 610016)

基于SURPAC-FLAC3D露天礦邊坡設(shè)計(jì)及穩(wěn)定性分析

陳銀萍1金愛(ài)兵1王 凱1馬 松2姚鵬飛1吳佐漢1

(1.北京科技大學(xué)土木與環(huán)境工程學(xué)院，北京 100083；2.四川省安全科學(xué)技術(shù)研究院,四川 成都 610016)

以金鑫銅鉬礦為工程背景，利用SURPAC和FLAC3D軟件建立礦山實(shí)體計(jì)算模型，同時(shí)獲得礦山剝采比和采出礦石量等經(jīng)濟(jì)數(shù)據(jù)。在不降低邊坡安全系數(shù)的前提下，通過(guò)適當(dāng)提高單臺(tái)階角度、高度和減小清掃平臺(tái)寬度3種途徑來(lái)增加露天礦最終邊坡角度，使得最終邊坡角提高2°，剝采比從3.66降低到3.52，減少巖石剝離量約1350萬(wàn)t，有效增加礦山企業(yè)經(jīng)濟(jì)效益。并在此基礎(chǔ)上，通過(guò)FLAC3D對(duì)加陡后的露天礦山邊坡進(jìn)行穩(wěn)定性分析，分析加陡之后的各邊坡水平位移、塑性區(qū)破壞分布情況，并對(duì)加陡后較為危險(xiǎn)的北邊幫進(jìn)行跟蹤點(diǎn)監(jiān)測(cè)。分析結(jié)果表明，各邊坡除小范圍局部不穩(wěn)定外，整體處于穩(wěn)定狀態(tài)，可以通過(guò)局部加固來(lái)達(dá)到礦山安全開(kāi)采的要求。

邊坡設(shè)計(jì) SURPAC FLAC3D生產(chǎn)剝采比 邊坡穩(wěn)定性

露天礦山開(kāi)采過(guò)程中，適當(dāng)增加邊坡角，可以有效降低平均剝采比。自上世紀(jì)80年代后期以來(lái)，國(guó)內(nèi)外部分礦山在進(jìn)行加陡最終邊坡角的研究[1]，以達(dá)到提高礦山企業(yè)經(jīng)濟(jì)效益的目的。據(jù)測(cè)算，一座中等規(guī)模的露天礦山，若最終邊坡角提高1°，即可減少巖石剝離量約1 000萬(wàn)m3[2]。在相同的穩(wěn)定條件下，加大露天礦邊坡角，可以減少剝離巖石[3]。尤其是對(duì)大型露天礦山，提高邊坡角是充分回收資源、降低生產(chǎn)成本、增加開(kāi)采效益的手段之一[4]。但加大邊坡角將帶來(lái)的邊坡穩(wěn)定性問(wèn)題。

1 工程概況

金鑫銅鉬礦是中鋼集團(tuán)投資興建的大型露天礦山。礦山生產(chǎn)能力300萬(wàn)t/a，設(shè)計(jì)服務(wù)年限20 a以上。礦體為隱狀礦體，賦存于火山穹窿構(gòu)造的北西部的白堊系粗面質(zhì)角礫凝灰?guī)r中，礦體呈不規(guī)則的長(zhǎng)扁豆?fàn)?，平面呈略向北突出的弧線，走向自西向東變化為40°～74°～105°。礦體主要成分為細(xì)粒鉀長(zhǎng)花崗巖、粗粒蝕變輝綠巖、閃長(zhǎng)巖、細(xì)粒蝕變輝綠巖、石英斑巖及花崗斑巖。

采場(chǎng)封閉圈標(biāo)高為696 m水平，696 m標(biāo)高以上為山坡露天礦，696 m水平以下為深凹露天礦。礦山原設(shè)計(jì)確定的平均剝采比為3.66，為了提高企業(yè)盈利能力，加陡露天礦邊坡從而降低剝采比是有效方法之一。

2 加陡方案

金鑫銅鉬礦露天開(kāi)采原設(shè)計(jì)方案：臺(tái)階高度為14 m，安全平臺(tái)寬為4 m，清掃平臺(tái)寬為10 m，單臺(tái)階坡面角為65°，運(yùn)輸平臺(tái)寬度為15 m，最小工作平臺(tái)寬度為45 m。初始境界平面圖見(jiàn)圖1。作露天采場(chǎng)A-A′剖面圖與I-I′剖面圖，如圖2、3所示。將露天礦的邊幫簡(jiǎn)化為圖中四個(gè)邊幫進(jìn)行分析，依次命名為東邊幫(A-A′右)、南邊幫(I-I′右)、西邊幫(A-A′左)、北邊幫(I-I′左)。

圖1 初始境界平面圖Fig.1 The plan of initial boundary

圖2 露天采礦A-A′剖面圖Fig.2 Section A-A′ of open-pit mining

圖3 露天采礦I-I′剖面圖Fig.3 Section I-I′ of open-pit mining

在不降低邊坡安全系數(shù)的前提下，根據(jù)礦山邊坡的實(shí)際參數(shù)，選擇具體加陡方案：?jiǎn)闻_(tái)階坡面角從65°提高到66°，臺(tái)階高度由14 m提高到15 m，安全平臺(tái)寬度不變，將清掃平臺(tái)寬度由10 m減小到8 m，運(yùn)輸平臺(tái)寬度15 m保持不變。加陡之后，東邊幫的最終邊坡角由50°提高到52°，南邊幫的最終邊坡角由43°提高到45°，西邊幫的最終邊坡角由45°提高到47°，北邊幫的最終邊坡角由49°提高到51°。東、北邊幫由于不存在運(yùn)輸平臺(tái)，所以最終邊坡角較南、西邊幫陡。

3 加陡前后采剝礦巖對(duì)比分析

根據(jù)加陡前后的邊坡參數(shù)，利用SURPAC軟件露天礦設(shè)計(jì)功能，建立加陡后的露天境界DTM見(jiàn)圖4。

圖4 加陡后露天境界DTMFig.4 The DTM of open-pit limit after increasing the overall slope angle

根據(jù)加陡前后的2種境界塊體模型，利用SURPAC軟件塊體模型的報(bào)告系統(tǒng)功能，分別計(jì)算出露天礦加陡前后的剝離巖石量、開(kāi)采礦石量和鉬金屬量如圖5～圖7和表1所示。

圖5 加陡前后剝離巖石對(duì)比Fig.5 Comparison of the amount of extracted rocks before and after increasing the overall slope angle■—加陡前；●—加陡后

圖6 加陡前后開(kāi)采礦石量對(duì)比Fig.6 Comparison of the amount of extracted ore before and after increasing the overall slope angle■—加陡前；●—加陡后

由圖5～圖7和表1可知，邊坡加陡之后，剝離巖石的總量減少了約1 350.62萬(wàn)t，剝采比從3.66減小到了3.52。

圖7 加陡前后開(kāi)采鉬金屬量對(duì)比Fig.7 Comparison of the amount of the extracted molybdenum before and after increasing the overall■—加陡前；◆—加陡后表1 加陡前后采剝礦巖量對(duì)比Table 1 Comparison of the amount of the extracted rock before and after increasing the overall slope angle

方 案剝離巖石總量/萬(wàn)t鉬礦石總量/萬(wàn)t剝采比加陡前15699.556835.813.66加陡后14348.936715.083.52

4 數(shù)值模擬

4.1 數(shù)值模擬模型

采用SURPAC與FLAC3D相結(jié)合的方式建立礦山計(jì)算模型，SURPAC提供了六面體單元形狀，根據(jù)地質(zhì)體的特征、計(jì)算精度要求以及單元形狀的空間展布特點(diǎn)，改變六面體單元的大小[5]。設(shè)置邊界條件，固定模型邊界所有點(diǎn)y方向速度，固定模型底部邊界的x、z方向速度、固定模型兩端的x方向速度。在開(kāi)挖模擬前將整個(gè)系統(tǒng)各個(gè)方向的位移和速度設(shè)置為0，即將初始應(yīng)力生成的變形去除，使得露天礦邊坡變形完全由露天開(kāi)挖引起。設(shè)置地質(zhì)體模型的力學(xué)參數(shù)得到加陡之后的計(jì)算模型見(jiàn)圖8，其中中間部分為開(kāi)挖部分，周邊部分表示最終留下的露天邊坡部分。

圖8 SURPAC-FLAC3D計(jì)算模型Fig.8 Computational model of SURPAC-FLAC3D

4.2 位移計(jì)算結(jié)果分析

當(dāng)開(kāi)挖到最終狀態(tài)時(shí)得到各邊幫在x方向的位移云圖如圖9所示。

x方向位移是指向坑內(nèi)的位移。從位移云圖可知：最大位移主要發(fā)生在北邊幫底部的局部幾個(gè)連續(xù)臺(tái)階處，為7.76 cm，其他3個(gè)邊幫指向坑內(nèi)的位移都較小，東邊幫最大位移為2.61 cm，南邊幫最大位移為6.01 cm，西邊幫最大位移為2.66 cm，均發(fā)生在邊幫底部的連續(xù)幾個(gè)臺(tái)階處，從圖中可以看出加陡之后各邊幫指向露天礦坑內(nèi)的水平位移都較小，平均位移在2～3 cm，從底部到頂部呈帶狀分布。

(a)東邊幫在x方向位移云圖

(b)南邊幫在x方向位移云圖

(c)西邊幫在x方向位移云圖

(d)北邊幫在x方向位移云圖 圖9 開(kāi)挖至最終狀態(tài)時(shí)各邊幫x方向位移云圖Fig.9 Contours of x-displacement on four slopes of final state

4.3 塑性區(qū)破壞分析

露天邊坡失穩(wěn)破壞可以看作是塑性區(qū)逐漸發(fā)展、擴(kuò)大直達(dá)貫通而進(jìn)入完全塑流狀態(tài)、無(wú)法繼續(xù)承受載荷的過(guò)程。若發(fā)生塑性變形的區(qū)域相互貫通，則邊坡有可能發(fā)生整體失穩(wěn)[6]。利用FLAC3D數(shù)值模擬得到開(kāi)挖到最終狀態(tài)時(shí)各邊幫的塑性區(qū)破壞分布圖，如圖10所示。

從塑性區(qū)破壞分布圖中可以看出4個(gè)邊幫的坡體內(nèi)部分均未發(fā)生明顯塑性變形。在開(kāi)挖計(jì)算過(guò)程中，4個(gè)邊坡的坡面都曾進(jìn)入過(guò)屈服狀態(tài)，其中坡體底部坡腳主要以拉伸破壞為主，坡面主要以剪切破壞和拉伸破壞為主。但邊坡只在極小范圍內(nèi)發(fā)生破壞，其中，東邊幫的坡腳處有約75 m2處于剪切破壞，邊坡坡面大約有20 m2處于拉伸破壞狀態(tài)。南邊幫的坡腳處大約有20 m2處于剪切破壞狀態(tài)狀態(tài)，邊坡坡面有2處面積分別約為15 m2和5 m2處于拉伸破壞狀態(tài)，但不貫通。西邊幫的坡腳處有約20 m2處于剪切破壞狀態(tài)，邊坡坡面有約60 m2處于拉伸破壞狀態(tài)。北邊幫只有邊坡坡面有1處約15 m2處于拉伸破壞狀態(tài)。但是由于這些塑性變形的區(qū)域沒(méi)有相互貫通，所以表明邊坡不會(huì)發(fā)生整體失穩(wěn)。只需要對(duì)各個(gè)邊坡坡面的幾個(gè)局部發(fā)生破壞的區(qū)域進(jìn)行加固即可。

(a)東邊幫塑性區(qū)破壞分布圖

(b)南邊幫塑性區(qū)破壞分布圖

(c)西邊幫塑性區(qū)破壞分布圖

(d)北邊幫塑性區(qū)破壞分布圖圖10 開(kāi)挖至最終狀態(tài)時(shí)各邊幫塑性區(qū)分布Fig.10 Distributions of plastic zone on four slopes of final state

為了更好地觀察開(kāi)挖對(duì)露天礦邊幫穩(wěn)定性的影響情況，須針對(duì)加陡后較為危險(xiǎn)的北邊幫進(jìn)行跟蹤點(diǎn)監(jiān)測(cè)。從北邊幫頂部、中部和底部依次設(shè)置了5個(gè)監(jiān)測(cè)點(diǎn)，編號(hào)分別為N1、N2、N3、N4、N5。提取這5個(gè)點(diǎn)在開(kāi)挖過(guò)程中的水平位移的變化量，得到圖11所示的監(jiān)測(cè)點(diǎn)的位移曲線圖。

由圖11可以看出，隨著開(kāi)挖步數(shù)的增加，監(jiān)測(cè)點(diǎn)指向坑內(nèi)的位移漸漸增加，但最終趨于穩(wěn)定。其中，位于邊坡頂部和中上部的監(jiān)測(cè)點(diǎn)N1、N2的水平位移偏向較小，位于邊坡中部的監(jiān)測(cè)點(diǎn)N3的位移變化不明顯，位于邊坡底部的2個(gè)監(jiān)測(cè)點(diǎn)的位移變化較大，最終穩(wěn)定在了4.5 cm左右。

圖11 北邊幫監(jiān)測(cè)點(diǎn)的位移曲線Fig.11 The displacement curve of monitoring sites on the northern slope■—N1；●—N2；▲—N3；▼—N4；?—N5

5 結(jié) 論

(1)利用SURPAC和FLAC3D軟件可以方便地建立符合礦山實(shí)際的數(shù)值計(jì)算模型，同時(shí)能獲得礦山剝采比、采出礦石量等經(jīng)濟(jì)數(shù)據(jù)。

(2)通過(guò)適當(dāng)加陡單臺(tái)階角度、提高臺(tái)階高度和縮短清掃平臺(tái)寬度來(lái)加陡金鑫銅鉬礦最終邊坡角度，使得生產(chǎn)剝采比從3.66減小至3.52。巖石開(kāi)挖量減小了約1 350萬(wàn)t。

(3)穩(wěn)定性數(shù)值模擬結(jié)果表明，加陡后露天礦各幫除局部區(qū)域外，整體處于穩(wěn)定狀態(tài)，通過(guò)局部加固的方法，完全可以滿足礦山安全開(kāi)采的要求。

[1] 李軍財(cái)，李角群，文 俊，等.大孤山鐵礦深凹露天礦邊坡工程[J].有色金屬，2004，56(3):86-89. Li Juncai，Li Jiaoqun，Wen Jun，et al.Rock slop engineering of Dagushan Open Pit Mine[J].Nonferrous Metals，2004，56(3):86-89.

[2] 姜立春，吳愛(ài)祥，李青松,等.基于陡幫開(kāi)采的新型安全清掃平臺(tái)穩(wěn)定性數(shù)值模擬分析[J].巖石力學(xué)與工程學(xué)報(bào)，2004，23(8):1263-1268. Jiang Lichun，Wu Aixiang，Li Qingsong，et al.Analysis and numerical simulation on stability of new safe-cleaning beach in mining with steep working slope[J].Chinese Journal of Rock Mechanics and Engineering，2004，23(8):1263-1268.

[3] 岳樹(shù)宇，張世雄，劉雁鷹,等.金堆城凹陷露天礦加陡深部邊坡角的空間原理研究[J].中國(guó)鉬業(yè)，2004，28(6):10-13. Yue Shuyu，Zhang Shixiong，Liu Yanying，et al.Study on spatial principle of steeping the slope angle in Jinduicheng Deep Open Pit[J].China Molybdenum Industry，2004，28(6):10-13.

[4] 楊天鴻，張鋒春，于慶磊,等.露天礦高陡邊坡穩(wěn)定性研究現(xiàn)狀及發(fā)展趨勢(shì)[J].巖土力學(xué)，2011，32(5):1439-1444. Yang Tianhong，Zhang Fengchun，Yu Qinglei，et al，Research situation of open-pit mining high and steep slope stability and its developing trend[J].Rock and Soil Mechanics，2011，32(5):1439-1444.

[5] 林 杭，曹 平，李江騰,等.基于SURPAC的FLAC3D三維模型自動(dòng)構(gòu)建[J].中國(guó)礦業(yè)大學(xué)學(xué)報(bào)，2008，37(3):339-342. Lin Hang，Cao Ping，Li Jiangteng，et al，Automatic generation of FLAC3Dmodel based on SURPAC[J].Journal of China University of Mining & Technology，2008，37(3):339-342.

[6] 林 杭，曹 平，李江騰,等.邊坡臨界失穩(wěn)狀態(tài)的判定標(biāo)準(zhǔn)[J].煤炭學(xué)報(bào)，2008，33(6):643-646. Lin Hang，Cao Ping，Li Jiangteng，et al.The standards for critical failure state of slope[J].Journal of China Coal Society，2008，33(6):643-646.

(責(zé)任編輯 石海林)

Slope Design and Stability Analysis of Open Pit Based on SURPAC and FLAC3D

Chen Yinping1Jin Aibing1Wang Kai1Ma Song1,2Yao Pengfei1Wu Zuohan1

(1.SchoolofCivilandEnvironmentEngineering，UniversityofScienceandTechnologyBeijing，100083Beijing，China；2.SichuanAcademyofSafetyScienceandTechnology，Chengdu610016，China)

With Jinxin copper molybdenum mine as engineering background，SURPAC and FLAC3Dare adopted to build the real calculation model of the mine，and the economic data such as the stripping ratio and the quantity of extracted ore are obtained simultaneously.On the premise of keeping safety，three ways of properly increasing the angle and height of single bench，and decreasing the width of cleaning bench can increase the final slope angle of open pit by 2°，drop the stripping ratio from 3.66 to 3.52 and lessen the amount of stripped rocks by 13.5 Mt.All these can efficiently improve the economic benefit of the mine enterprise.Based on this，the stability of the open pit slope after increasing the overall slope angle is analyzed，and its horizontal displacements and failures of the plastic zones are demonstrated.Real-time monitoring on the northern sidewall that is in a high risk is made.The analysis indicated that each new slope is basically in the stable condition and only small part of slope requires local ground support to guarantee the safety in mining.

Slope design，SURPAC，F(xiàn)LAC3D，The stripping ratio，Slope stability

2014-05-11

陳銀萍(1988—)，女，碩士研究生。

TD325，TD854.6

A

1001-1250(2014)-10-166-05