劉 泉

(四川省冶金設計研究院，四川 成都610041)

邊坡是露天礦山的安全隱患［1］，在極端惡劣環(huán)境下(如大量降水、地震和爆破等)，露天礦山邊坡的軟弱結構面極易發(fā)生傾倒和滑動破壞［2］。邊坡穩(wěn)定性評價是當今研究熱點，然而，以往研究只針對某一惡劣條件(如暴雪［3］、降雨［4］、地震［5］等)，沒有考慮到多重惡劣條件同時發(fā)生的可能性。西藏高原某露天銅礦，海拔在4 600 m 以上，且礦區(qū)邊坡可能同時面對幾種極端惡劣環(huán)境的考驗，本研究針對這一情況，以礦區(qū)邊坡A－A 剖面為典型研究對象，采用ANSYS 軟件［6］建立礦區(qū)A －A 剖面的計算模型，在普通開挖條件下加載積雪和地震荷載，運用FLAC3D軟件［7］分析邊坡及臺階面的穩(wěn)定性。

1 工程概況

西藏高原某露天銅礦，海拔在4 600 m 以上，邊坡表面常年積雪，人工開挖邊坡最高達1 000 m 以上，為國內罕見，露采邊坡的物質組成以巖石為主，節(jié)理裂隙發(fā)育，規(guī)模雖不大，但它是控制臺階邊坡穩(wěn)定性的主要內在因素。斷層、層理和片理延伸較長，規(guī)模較大，是影響大邊坡穩(wěn)定性的重要因素。A －A 剖面(圖1)為礦區(qū)的典型地質剖面，地形起伏較大，第四系松散堆積物較發(fā)育，山坡表面密灌叢生或為腐殖覆蓋或為寒凍風化形成的倒石堆覆蓋，厚度一般在10 ～30 m 不等;A－A 剖面下部主要為凝灰?guī)r與閃長玢巖等，第四系碎石土與凝灰?guī)r之間發(fā)育有厚度5 m不等的全強風化層。

圖1 A－A 工程地質剖面(單位:m)Fig.1 A-A geological section

2 參數及模型建立

2.1 巖體力學參數選取

礦山邊坡長期受到極端惡劣條件的影響，其力學性質比中深部巖土體差，為了增加模擬的可靠性，根據工程經驗對實驗數據進行折算，折算系數取0.82［8］，折減結果如表1，其中，E 為彈性模量，μ 為泊松比，ρ 為密度，C 為抗剪強度，φ 為內摩擦角。

表1 巖體力學參數推薦值Table 1 Recommended values of rock mass mechanics parameters

2.2 計算模型

建立A－A 剖面地質模型，安全平臺6 m，清掃平臺12 m，接渣平臺30 m;巖質臺階坡面角70°，碎石土覆蓋層坡面角45°，一期最終邊坡角38.5°。采用ANSYS 軟件建立A－A 剖面邊坡計算模型，然后導入FLAC3D進行后處理，形成計算模型，如圖2 所示，建立了詳細的邊坡空間狀形態(tài)及地層結構。共劃分單元4 119 個、節(jié)點9 113 個。模型對于部分地區(qū)進行了網格細化。

圖2 A－A 剖面計算模型Fig.2 A-A profile calculation model

2.3 添加荷載

(1)根據相關邊坡規(guī)范，在考慮地震載荷的情況下，可以忽略爆破載荷的作用。礦區(qū)地震基本烈度為Ⅶ度，據相關地震規(guī)范，取地震加速度為0.125g，采用擬靜力法進行分析，只考慮滑動方向的水平地震力的作用。地震荷載取值:

式中，U 為地震載荷矩陣;A 為單元質量矩陣;E 為地震加速度，取0.125g。

(2)考慮到研究區(qū)域常年被積雪覆蓋，且有凍土存在，邊坡表面均布積雪荷載(垂直于坡面)，可根據相關規(guī)范和式(2)進行計算:

式中，X 為積雪荷載，Pa;h 為積雪深度，m，取0.6 m;g為自由落體加速度，m/s2;ρ 為積雪的平均密度，kg/m3，取120 kg/m3;λ 為附加凍土系數，取1.1;α 為坡面角，取45°。經過計算，X =6 789 Pa。

3 開挖后邊坡計算結果分析

3.1 位移分析

邊坡開挖后，僅考慮原巖應力的情況下，得到了圖3 所示位移云圖(碎石土臺階坡面角為45°)，在4 729 ～4 759 m 高度范圍內出現(xiàn)了最大位移，約為84.6 mm。從上面的整體分析來看在碎石土區(qū)域臺階邊坡穩(wěn)定，未出現(xiàn)大的變形位移，且整體邊坡穩(wěn)定。

圖3 開采過程中位移云圖Fig.3 Displacement diagram in the mining process

3.2 剪應變增量分析

從圖4 可以看出，碎石土臺階最大剪應變增量小于0.009 36，巖石臺階面的剪應變增量不足0.001，整體邊坡以及局部碎石土區(qū)域臺階邊坡穩(wěn)定。

圖4 開采過程中剪應變增量云圖Fig.4 The shear strain increment cloud in the mining process

3.3 塑性區(qū)分析

從圖5 的特征來看，巖質臺階面部位以及碎石土邊坡區(qū)域無塑性區(qū)，說明在修改碎石土臺階坡面角情況下，整體邊坡以及局部碎石土區(qū)域臺階邊坡穩(wěn)定。

圖5 開采過程中塑性區(qū)分析圖Fig.5 The plastic zone in the mining process

4 積雪和地震載荷下計算結果分析

在初始地應力場條件下，按照一期設計開挖順序，進行逐步精細開挖計算，通常進行10 多次開挖計算，并把地震和積雪荷載和劣化因素加入到計算條件中，主要用于極端校核。

4.1 邊坡開挖后積雪和地震載荷下位移分析

從圖6 中可知(模型中碎石土臺階坡面角為45°):在4 729 ～4 759 m 高度范圍內出現(xiàn)了最大位移，約為366.8 mm，方向為指向坡體外面。從上面A－A 剖面邊坡開挖后積雪和地震載荷下整體位移分析來看在碎石土區(qū)域臺階邊坡欠穩(wěn)定，局部出現(xiàn)大的變形位移，巖質臺階面部位移量很小;A－A 剖面邊坡開挖后積雪和地震載荷下，碎石土臺階坡面局部有大的變形，但整體邊坡穩(wěn)定。

4.2 邊坡開挖后積雪和地震載荷下剪應變增量分析

從圖7 可以看出，巖質臺階面部位剪應變增量很小，為0.01 以下，碎石土臺階邊坡局部區(qū)域最大剪應變增量小于0.082，故碎石土區(qū)域出現(xiàn)很小的剪切應變，A－A 剖面邊坡開挖后積雪和地震載荷下局部碎石土區(qū)域臺階邊坡欠穩(wěn)定，但整體邊坡穩(wěn)定。

圖6 邊坡開挖后積雪和地震載荷下位移云圖Fig.6 Displacement diagram in the mining process under the heavy snow and the seismic load

圖7 邊坡開挖后積雪和地震載荷下剪應變增量云圖Fig.7 The shear strain increment cloud in the mining process under the heavy snow and the seismic load

4.3 邊坡開挖后積雪和地震載荷下塑性區(qū)分析

從圖8 的特征來看，巖質臺階面部位有少量塑性區(qū)，碎石土邊坡區(qū)域出現(xiàn)了較大的塑性區(qū)，說明在修改碎石土臺階坡面角情況下，A－A 剖面邊坡開挖后，在積雪和地震載荷下局部碎石土區(qū)域臺階邊坡欠穩(wěn)定，但整體邊坡穩(wěn)定。

圖8 邊坡開挖后積雪和地震載荷下塑性區(qū)分析圖Fig.8 The plastic zone in the mining process under the heavy snow and the seismic load

5 結 論

(1)西藏某露天礦邊坡可能同時面臨多種極端惡劣開采條件，采用ANSYS 軟件建立礦區(qū)邊坡A －A剖面的計算模型，同時加載積雪和地震荷載，運用FLAC3D軟件分析邊坡及臺階面的穩(wěn)定性。

(2)邊坡穩(wěn)定性計算參數，采用巖土力學室內試驗和經驗類比法確定，類比系數為0.82。根據國家抗震規(guī)范和地質詳勘報告獲得該礦區(qū)地震設防烈度按Ⅶ度考慮，采用擬靜力法計算的地震加速度按0.125g 計，考慮到研究區(qū)域常年被積雪覆蓋，且有凍土存在，邊坡表面均布積雪荷載X=6 789 Pa，采用安全平臺6 m，清掃平臺12 m，接渣平臺30 m;巖質臺階坡面角70°，碎石土覆蓋層坡面角45°，一期最終邊坡角38.5°的參數設計。

(3)由A－A 剖面在開采過程中各種工況下的計算結果可知，碎石土區(qū)域臺階邊坡及整體邊坡保持穩(wěn)定。在積雪和地震校核工況作用下，其穩(wěn)定性有所降低，但整體邊坡仍能保持穩(wěn)定。這說明最終臺階坡面角38.5°的設計能保證邊坡穩(wěn)定，且能減少剝離量，在保證安全的情況下能實現(xiàn)效益的最大化。

［1］ 胡 斌，王 偉，張 騰，等. 露天礦最終邊坡角優(yōu)化設計模擬分析［J］.金屬礦山，2013(9):14-18.

Hu Bin，Wang Wei，Zhang Teng，et al.Simulation analysis of the final sloping angle design and optimization on the open-pit mine［J］.Metal Mine，2013(9):14-18.

［2］ 劉宇飛.基于MATLAB 的土質邊坡穩(wěn)定可靠度計算及應用研究［D］.長沙:中南大學，2008.

Liu Yufei. Reliability Calculation and Application Research of Soil Slope Based on MATLAB［D］.Changsha:Central South University，2008.

［3］ 許玉娟.巖石凍融損傷特性及寒區(qū)巖質邊坡穩(wěn)定性研究［D］.長沙:中南大學，2012.

Xu Yujuan.Study on Degradation Characteristics of Rock in Condition of Freeze-thaw Cycle and Rock Slope Stability Analysis in Cold Region［D］.Changsha:Central South University，2012.

［4］ Zhang Zhaocai，Peng Chao，Yan Zhenxiong. Effects of rain infiltration on fissured slope stability［J］.Electronic Journal of Geotechnical Engineering，2014(19):2159-2164.

［5］ Zhang Yingbin，Chen Guangqi，Wu Jian，et al.Numerical simulation of seismic slope stability analysis based on tension-shear failure mechanism［J］.Geotechnical Engineering，2012，43(2)18-28.

［6］ Tan Bo，Zhang Rui，Lai Yanting.Analysis on stability of the expensive soil cutting slope based on ANSYS［J］.Applied Mechanics and Materials，2013，405/406/407/408:299-303.

［7］ Guo Fang，Gu Wei，Tang Jiejun，et al.Research on deformation stability of soft rock slope under excavation based on FLAC3D［J］.Applied Mechanics and Materials，2013，275/276/277:290-294.

［8］ 劉 春.邊坡工程中巖石力學參數隨機模糊選取研究［J］.巖土力學，2004，28(8):1327-1329.

Liu Chun.Research on random-fuzzy method for selecting rock parameters in slope engineering［J］.Rock and Soil Mechanics，2004，28(8):1327-1329.