基于Geo-Studio的某邊坡穩定性分析

周忠志　陳換換

摘 要：本文以準東露天煤礦首采區非工作幫軟弱煤層底板邊坡作為研究對象，通過分析該礦區工程地質條件和邊坡特征，得出堆土高度為0 m、20 m、40 m、60 m、80 m、100 m和120 m七種工況的等效附加荷載值。采用二維有限元數值模擬軟件Geo-Studio對不同工況進行模擬分析，結果表明，在坡頂排土場堆積土體高度不斷增加的情況下，在非工作幫470 m平臺以上的坡體剪應力集中嚴重，邊坡失穩的可能性增大，得出坡頂堆載對影響軟弱煤層底板邊坡失穩的規律，對準東煤礦非工作幫邊坡做出科學評價。

關鍵詞：應力應變分析;邊坡穩定性;穩定系數;Geo-Studio

中圖分類號：TU457 ? ? 文獻標志碼：A ? ? 文章編號：1003-5168（2022）3-0063-04

DOI：10.19968/j.cnki.hnkj.1003-5168.2022.03.015

Analysison Stability of Slope Based on Geo-Studio Software

ZHOU Zhongzhi? ? CHEN Huanhuan

（1.Qingdao Housing Security Center，Qingdao 266000，China;2.Qingdao Real Estate Registry，Qingdao 266000，China）

Abstract：In this paper， the floor slope of non working side weak coal seam in the first mining area of Zhundong open pit coal mine is taken as the research object. By analyzing the engineering geological conditions and slope characteristics of the mining area， the equivalent additional load values of seven working conditions with soil stacking heights of 0 m， 20 m， 40 m， 60 m， 80 m， 100 m and 120 m are obtained. The two-dimensional finite element numerical simulation software Geo-Studio is used to simulate and analyze different working conditions. The results show that with the increasing height of the accumulated soil in the slope top dump， the shear stress concentration of the slope above the 470 m platform of the non working slope is serious， and the possibility of slope instability is increased. The law of slope top surcharge affecting the slope instability of the floor of weak coal seam is obtained， The non working slope of Zhundong coal mine is scientifically evaluated.

Keywords：stressstrainanalysis;slopestability;safetyfactor;Geo-Studio

0 引言

20世紀，許多露天礦煤礦主急于取得經濟效益，并且相關規劃設計人員不顧及煤礦長遠發展，因此，部分露天礦存在設計上的不合理。其中，以排土場緊鄰露天礦邊坡問題最為突出。隨著露天礦開采量的日益增長，排土場上堆積的排棄物逐漸增多，緊鄰排土場的坡體下滑力增大，邊坡失穩的可能性上升，從而增大滑坡災害的可能性。因此，研究坡頂堆載對邊坡穩定性的影響具有重大意義[1-4]。

本文以準東露天煤礦首采區非工作幫邊坡為例，在分析邊坡巖土體力學特性的基礎上，采用Geo-Studio軟件中的模擬分析，得到坡體在不同高度的坡頂堆載作用下剪應力云圖，從而得出了邊坡失穩的內因[5-7]。再通過Geo-Studio軟件對非工作幫邊坡進行穩定性分析，計算得到坡體在不同高度的坡頂堆載作用下的穩定系數，得出引起邊坡失穩時的極限堆積土體高度[8-10]。為準東露天煤礦的邊坡治理提供理論依據，并提出具體的預防和治理措施。

1 工程地質條件

準東露天煤礦位于新疆維吾爾自治區吉木薩爾縣，距吉木薩爾縣的距離為130 km，距離烏魯木齊市的距離為220 km[11]。礦區處于準噶爾盆地腹地內，地貌類型主要以沖積湖平原、剝蝕殘余丘陵、風積沙漠和剝蝕波狀平原為主。礦區年平均降水量為106 mm，年蒸發量為1 202～2 382 mm，5—8月偶有雷陣雨，冬季積雪稀少。區內常年多風，風力一般為4～5級，經常有7～8級大風，最大可達10級，多以西北風為主。降雨量稀少，屬于大陸性沙漠氣候[12-13]。

首采區現狀：露天坑南北向長約為2.2 km，東西向寬約為1.4 km。準東露天煤礦首采區的西側為工作幫，在工作幫上設計臺階。其中，542 m水平以上，臺階高度為10 m;542 m水平以下，臺階高度為16 m。非工作幫不設計臺階，邊坡角為8～13°，首采區最大采深320 m。本文的研究對象為準東露天煤礦非工作幫邊坡，經過不斷露天開采之后，首采區的Bm煤層已經得到充足開采，非工作幫邊坡即經已開挖完煤層之后得到的Bm煤底板，主要是以泥巖為主，厚度僅為0.5～5 m。

2 應力應變數值模擬分析

為了更好展現坡體在不同高度的堆積土體作用下應力應變變化，采用二維有限元模擬軟件Geo-Studio中的SIGMA/W模塊進行模擬，得到坡體剪應力場的變化特征。

2.1 計算參數的選取

計算參數主要來源于神華新疆公司提供的《新疆準東煤田-吉木薩爾縣五彩灣礦的勘探報告》，見表1。

2.2 計算工況的確定

通過實地踏勘，排土場堆積物以碎石、開采的煤層塊和煤砂泥巖為主，計算得出土體容重為25 kN/m，在進行模擬分析時，視為堆積物近似均勻地分布于排土場上。依據準東煤田五彩灣礦區三號露天煤礦一期工程設計，排土場設計最大堆土高度為120 m，分析時取堆土高度為0 m、20 m、40 m、60 m、80 m、100 m和120 m，得到了不同高度下坡頂堆載的等效附加荷載值，見表2。

2.3 建立模型

因為研究對象為非工作幫邊坡即，經已開挖完煤層之后得到的Bm煤底板，主要以泥巖為主，在建模型時僅考慮了泥巖層。依據研究對象底部X、Y方向無位移，左側和右側X方向無位移，Y方向可自由移動，定義研究對象邊界條件，建立模型。

2.4 數值模擬結果及分析

為了研究坡體在不同高度的堆積土體作用下的變化特征和穩定性發展趨勢，同時考慮滑坡發生的潛在模式，模擬堆土高度為0 m、60 m、120 m時三種工況下的剪應力云圖，如圖1至圖3所示。

在坡頂排土場堆積土體高度不斷增加的過程中，坡頂、便道493 m臺階和首采區坑底位置出現了明顯的應力集中現象。隨著每增加20 m高堆積土體，坡頂處的剪應力值由50 kPa增加到400 kPa;便道493 m臺階位置剪應力值由100 kPa增加到150 kPa;首采區坑底位置處的剪應力值由200 kPa增加到250 kPa。應力集中現象在坡頂位置處最為明顯。因此，隨著排土場堆積土體高度的增加，在470 m平臺以上的邊坡極易發生失穩。

通過以上SIGMA/W模塊對坡體進行應力應變的模擬分析得到：坡頂排土場堆積土體荷載主要引起470 m臺階以上坡體應力集中嚴重和應變明顯增大，造成上部土體的滑動變形，從而擠壓下部坡體產生變形而失穩，從滑坡受力類型分析，屬于后推式滑坡。

3 穩定性計算及結果分析

按照《巖土工程勘察規范》取Fs=1.20作為邊坡安全系數的參考值。為了更加明確邊坡在不同高度堆積土體下穩定狀況，再采用SLOPE/W模塊進行穩定性計算，分別得到7種不同堆積高度下的安全系數，為1.339、1.314、1.275、1.238、1.110、1.014、0.971，整理得到穩定系數曲線圖，如圖4所示。

從圖4可以得出，坡體在堆積土體為0 m時的穩定系數時1.339，與選定的邊坡安全系數的參考值Fs=1.20比較，坡體處于穩定狀態。隨著坡頂堆積土體高度由20 m增加到120 m后，坡體的穩定系數由1.314下降到0.971。通過分析可以得到，隨著坡頂排土場堆積土體高度的不斷增加，坡體的穩定性系數與土體堆積高度呈負相關，即堆積土體高度增加時坡體的穩定安全系數隨之減小。

在坡頂排土場堆積土體的高度低于60 m時，坡頂排土場堆積土體對邊坡穩定性影響較小。如在堆積土體高度為20 m時，經與0 m時的邊坡進行分析對比研究，發現此種工況下邊坡的穩定系數與0 m下基本相同。在堆積土體高度由20 m到60 m逐漸增大的過程中，邊坡穩定系數雖然在逐漸降低，但都是高于選定的邊坡安全系數的參考值Fs=1.20，表明坡體在穩定狀態。分析得出原因，在堆積土體未達到一定高度時，邊坡的失穩與坡頂排土場堆積土體相關性較小，其他影響因素在對邊坡的穩定性起主要作用，如開采爆破震動和軟弱夾層。在坡頂排土場堆積土體高度由60 m增大80 m時，邊坡穩定系數由1.238快速下降到1.110，低于選定的邊坡安全系數參考值Fs=1.20，坡體處于不穩定狀態。當土體堆積高度增大到80 m時，整個邊坡的滑面上移，造成滑體由坡體的中部滑出。分析得出原因，滑面的上移造成了抗滑力降低，滑體由坡體中部剪出，此時坡頂排土場堆積土體對邊坡失穩起到主要作用。

隨著坡頂排土場土體堆積高度的增大，邊坡的穩定系數不斷降低，兩者是負相關。并且可以看出，在堆積土體的高度為0～60 m時，邊坡穩定系數下降幅度較小;當堆積土體高度由60 m增加到80 m時，邊坡穩定系數下降幅度較大，并且已經低于選定的邊坡安全系數參考值Fs=1.20。分析其原因，隨著坡頂排土場堆積高度的增加，滑面明顯上移，坡頂荷載直接將滑體由邊坡中部壓出，此時，滑體上下滑力增大，抗滑力減小。

4 結語

采用SIGMA/W和SLOPE/W模塊相結合的數值模擬手段，從剪應力場和穩定系數對坡體穩定性進行研究，結果表明。

①坡頂排土場堆積土體荷載主要引起470 m臺階以上坡體應力集中嚴重和應變明顯增大，造成上部土體的滑動變形擠壓下部坡體產生變形而失穩，屬于后推式滑坡。

②隨著坡頂排土場堆積土體高度的增加，造成邊坡穩定系數從1.339減小到0.971，穩定系數逐漸降低，邊坡失穩的趨勢可能性增大。

③當坡頂排土場堆積土體高度達到60 m時，穩定系數已接近極限平衡值，當達到80 m以上時，穩定系數已經低于選定的安全系數參考值Fs=1.20，因此，當堆積高度達到60 m時，要預防邊坡失穩的發生。

④因堆土荷載的存在，需要對準東露天煤礦非工作幫邊坡不穩定的區域進行加固治理，通過布置抗滑樁進行加固處理，并引入監測預警預報系統，對邊坡實時監控。

參考文獻：

[1]丁參軍，張林洪，于國榮，等.邊坡穩定性分析方法研究現狀與趨勢[J].水電能源科學，2011（8）：112-114.

[2]黃昌乾，丁恩保.邊坡工程常用穩定性分析方法[J].水電站設計，1999（1）：53-58.

[3]李銀海.邊坡工程常用穩定性分析方法[J].中國水運，2005（5）：94-95.

[4]李曉敏，宋靜.邊坡工程常用穩定性分析方法綜述[J].巖土工程與地下工程，2011（3）：88-92.

[5]姜清輝，王笑海.三維邊坡穩定性極限平衡分析系統軟件SLOPE3D的設計及應用[J].巖石力學與工程學報，2003（7）：1121-1125.

[6]韋樓，蔣廷球.高邊坡穩定性分析對比研究：簡化BISHOP法和簡化JANBU法[J].中國西部科技，2011（12）：8-10.

[7]曾繁文.基于GEO-SLOPE軟件對排土場邊坡穩定性的分析[J].2013（2）：182-183.

[8] BREBBIA C A ，TELLES J C F ，WROBEL L C，Boundaryelementtechnique[J].BerlinHeidelbergNewYorkTokyo，1984.

[9]涂杰楠，李劍.Geo-Studio在黃土邊坡中的應用[J].山西建筑，2007（35）：359-360.

[10]田洪寧，牛彥博，羅剛.滑坡穩定性分析中Geo-studio軟件的應用[J].2012（16）：32-33.

[11]陸衛國.準東露天礦采場軟巖邊坡穩定性研究[J].露天采礦技術，2010（5）：53-54.

[12]韓勇.準東露天礦排土場邊坡穩定性分析研究[J].露天采礦技術，2010（5）：63-64.