易露，張華

(大冶有色金屬有限責任公司銅山口銅礦， 湖北 大冶市 435122)

0 引言

銅山口銅礦是大冶有色金屬有限責任公司下屬的一座年產(chǎn)量上百萬噸的大型露天銅礦山，隨著多年開采，露天坑越來越深，露天邊坡也越來越高，面臨的安全風險也隨之增加[1]，有必要對其穩(wěn)定性進行評價[2]。當前，對露天礦山邊坡穩(wěn)定性的研究較多，鄧文建等[3]使用數(shù)值模擬的方法研究了降雨對邊坡穩(wěn)定性的影響；李學鋒等[4]使用Rockfall數(shù)值軟件分析了小型露天采石場高陡邊坡滾石的危險性；張云鵬等[5]通過監(jiān)測觀測點爆破后的位移情況，使用極限平衡法分析了邊坡在爆破影響下的安全性；吳朝松等[6]研究了研山鐵礦東幫邊坡在地下水長期作用下的穩(wěn)定性。他們的研究往往是針對邊坡某種工況下的安全性，因此不夠全面。本文設(shè)定了包括自重、地下水、爆破和地震的4種組合工況，采用極限平衡法計算了銅山口銅礦露天邊坡穩(wěn)定性，并運用FLAC3D數(shù)值模擬對邊坡進行建模，驗證了數(shù)值計算的結(jié)果，證實銅山口銅礦露天坑邊坡是穩(wěn)定的。

1 工程概況

銅山口銅礦礦區(qū)面積約2.88 km2，東西方向，長度為1008 m，埋藏深度從-196 m至地表，礦體上表層覆蓋較薄，露天開采是最佳開采方法。經(jīng)過多年開采，已經(jīng)形成了一個南高北低的巨大露天坑，最大邊坡高度達180 m，地面以下深度達54 m。邊坡臺階的高度為12 m，坡面角為65°～70°。邊坡巖體主要為碳酸鹽巖和巖漿巖，石炭、二疊、三疊系碳酸鹽巖層是區(qū)內(nèi)主要含水層，以巖溶裂隙水及裂隙巖溶水為主。根據(jù)露天坑實際情況，將其邊坡按照方位分為東（N）、南（S）、西（W）、北（E）4個區(qū)，如圖1所示，各分區(qū)的邊坡中N區(qū)露采邊坡較緩；W區(qū)、E區(qū)露采邊坡較陡，邊坡高度約150～180 m。

圖1 露天坑邊坡現(xiàn)狀

2 邊坡穩(wěn)定性分析

銅山口銅礦露天邊坡穩(wěn)定性極限平衡分析采用的是垂直裂隙加圓弧型滑動模式，采用Slide 6.0軟件自動計算拆分和搜索潛在滑面位置并計算潛在滑動面的安全系數(shù)。主要分為3個步驟：

（1）根據(jù)邊坡實際情況，選擇邊坡的分析剖面；

（2）確定邊坡坡面的工況條件和對應(yīng)的安全系數(shù)；

（3）計算邊坡坡面的安全系數(shù)，分析邊坡坡面的穩(wěn)定性。

2.1 選擇邊坡分析剖面

建立剖面是建立地質(zhì)模型進行邊坡穩(wěn)定性研究的基礎(chǔ)，銅山口露天邊坡穩(wěn)定性研究分析剖面的布置是按照工程地質(zhì)分區(qū)進行。根據(jù)露天邊坡的實際長度和條件，在N、E、S、W 4個不同邊坡區(qū)域共布置剖面10條，其中：N區(qū)3條；E區(qū)2條；S區(qū)2條；W區(qū)3條，如圖2中所示的N1～N3、E1～E2、S1～S2、W1～W3線剖面。部分剖面的斷面如圖3所示，剖面的邊坡高度為43～218 m，剖面線的巖性依據(jù)銅山口銅礦所做的地質(zhì)詳查資料。

圖2 分析剖面布置

圖3 邊坡部分剖面

2.2 邊坡的工況

對10個剖面進行表1中所示的4種載荷組合下的受力情況分析，根據(jù)工程經(jīng)驗[7]確定安全系數(shù)的許用值[K]，在計算出的安全系數(shù)K＞許用值[K]時邊坡穩(wěn)定；1＜K＜[K]時邊坡基本穩(wěn)定，但是預(yù)留的安全系數(shù)不足，在發(fā)生極端天氣時可能存在其他風險導(dǎo)致邊坡失穩(wěn)；在K＜1時，認為邊坡處于不穩(wěn)定狀態(tài)，會隨時失穩(wěn)。

表1 4種工況下的許用安全系數(shù)[K]

2.3 剖面穩(wěn)定性分析

對于劃分的 10個邊坡剖面，利用軟件自動對邊坡進行分層，選擇極限平衡法對各個分層的安全系數(shù)進行計算，圖4為是W1剖面的計算結(jié)果，圖中工況 1，2，3，4條件下，表層邊坡的安全系數(shù)依次是1.691，1.639，1.580和1.127。分別計算10個剖面的安全系數(shù)，見表2。

圖4 W1剖面設(shè)計邊坡受力4工況下安全系數(shù)穩(wěn)定性分析

表2中10個剖面在4種工況下計算出來的極限安全系數(shù)都大于表1中設(shè)定的許用安全系數(shù)[K]，表明現(xiàn)狀邊坡整體是穩(wěn)定的。此外，現(xiàn)狀條件下，生產(chǎn)爆破對邊坡穩(wěn)定性的影響約為 5%，天然地震對邊坡穩(wěn)定性的影響約10%，即工況2，3較工況1的安全系數(shù)分別降低約5%～10%；生產(chǎn)爆破對邊坡穩(wěn)定性的影響約 5%，天然地震對邊坡穩(wěn)定性的影響約9%，即工況2，3較工況1的安全系數(shù)分別降低約5%～9%；降雨對邊坡穩(wěn)定性的影響為降低安全系數(shù)約3%～15%。

表2 邊坡各剖面的安全系數(shù)計算結(jié)果

3 數(shù)值模擬

3.1 構(gòu)造模型

將繪制的露天邊坡結(jié)構(gòu)圖利用kubrix進行工程結(jié)構(gòu)的網(wǎng)格劃分，并導(dǎo)入FLAC3D中生成銅山口礦露天邊坡數(shù)值分析模型[8]，采用中空單元模型和Mohr-Coulomb的2種本構(gòu)模型建立邊坡模型，礦區(qū)的東西方向共1008 m，模型范圍x=887～1895 m；礦區(qū)的南北方向共1225 m，模型范圍y=878～2103 m；模型范圍z=-196 m至地表。模型劃分單元2 791 008個，節(jié)點388 284個，模型部分剖面如圖5所示。圖5中不同的顏色表示不同的巖性，假定每一種巖土體都為各向同性的連續(xù)介質(zhì)。

圖5 整體模型和巖性

3.2 分析結(jié)果

圖6為整體模型的應(yīng)力和應(yīng)變圖，從圖6可知，模型的應(yīng)力大小受地形的影響，隨深度的增加而增大，越接近坡面主應(yīng)力值越小，其方向從坡體深部的垂直方向逐漸轉(zhuǎn)為淺部的順坡方向；受到巖性條件影響，在巖層交接部位、應(yīng)力的擾動狀態(tài)較大；在坡面附近、地表和巖層淺部存在一定的拉應(yīng)力，其中東部高邊坡面附近存在范圍較大的連續(xù)拉應(yīng)力區(qū)，最大拉應(yīng)力值約為0.48 MPa。模型邊坡巖體的位移主要是向采坑臨空面移動，西北坑底附近巖體向上隆起，最大隆起高度有 10.2 cm；東部由于地勢比較高，東部山坡最高標高為180 m，預(yù)留臺階比較窄，存在兩三個臺階并段，東部位移最大達到0.25 m。坡面巖體附近的位移值不大，最大值只有幾厘米，表明邊坡穩(wěn)定，無安全風險。

圖6 模型應(yīng)力和位移

4 結(jié)論

（1）采用極限平衡法分析評價銅山口銅礦開采現(xiàn)狀邊坡的10個剖面在4種工況下的穩(wěn)定性。結(jié)果表明：計算的邊坡穩(wěn)定性安全系數(shù)都大于規(guī)定的許用安全系數(shù) 1.2，銅山口銅礦的現(xiàn)狀邊坡是穩(wěn)定的。

（2）采用極限平衡法和FLAC3D數(shù)值模擬方法相結(jié)合來分析露天礦邊坡穩(wěn)定性，將理論和模擬相結(jié)合是一種合理可行的方法，為露天礦山提供更為可靠的邊坡穩(wěn)定性邊坡評價及參數(shù)設(shè)計方案。

（3）邊坡巖體的強度較高，整體邊坡的高度為中邊坡，分析安全系數(shù)較高，邊坡穩(wěn)定好。