高偉程
(中國石油天然氣管道工程有限公司 沈陽分公司,遼寧 沈陽 110031)
隨著我國露天煤礦開采深度不斷增加,高陡邊坡日趨常態化,高陡邊坡問題日漸突出,高陡邊坡穩定性成為制約安全與經濟的關鍵因素之一[1]。而軟弱夾層是地層巖體中的非連續性面,其物理力學性質較差,往往是造成邊坡滑坡主要因素之一[2-3]。因此,探討含軟弱夾層高陡邊坡變形破壞機理,有利于高陡邊坡穩定性研究及危險邊坡治理。目前,國內外專家學者對露天邊坡已經進行了一系列研究。例如呂金星等[4]以勝利露天煤礦邊坡為例,基于數值模擬分析邊坡漸進性破壞,研究表明巖石的流變特性引起邊坡漸進破壞,最終導致滑坡;陶志剛等[5]基于高陡邊坡物理模擬試驗,確定邊坡滑坡時內部應力破壞規律及NPR 錨索力學特性,并證實NPR 錨索具有超前預警功能;何怡等[6]以四川某含弱層高邊坡為研究對象,基于應變軟化M-C 模型分析礦山邊坡漸進破壞過程,明確弱層含水量對邊坡破壞的影響;范昊天等[7]以蘇州某礦區滑坡區為研究對象,基于顆粒流離散元法建立PFC 邊坡模型,分析含軟弱夾層巖質邊坡破壞機理及滑動時應力狀態,為邊坡治理提供理論基礎。
基于平朔東露天礦以往勘探資料,得出露天礦南幫主要巖性有砂巖、泥巖、泥砂巖互層、粉質黏土、煤層及弱層。其中弱層由富含高嶺土等黏土礦物組成,在地下水的長期浸潤下,力學強度極低,嚴重影響邊坡穩定。同時根據現場情況,在底部1245 臺階上有明顯錯動跡象,在上部1 345 m 水平有裂縫顯現,存在邊坡發生大規模滑坡的危險。
根據礦區工程地質條件、以往地質鉆孔資料、滑坡反分析及節理化強度分析,以初步設計為依據,構建典型工程地質模型,其地層依次為粉質黏土、泥巖、泥砂巖互層、砂巖、弱層、泥巖、煤等,如圖1,相關巖土力學參數推薦值見表1。

圖1 典型工程地質剖面

表1 巖土力學參數推薦值
露天邊坡工程是一門復雜的涉及多學科的邊緣學科,其集采礦、安全、地質、力學等為一體[8]。露天邊坡失穩通常受多因素影響,其穩定性研究極為復雜而且十分困難,不僅寄托于理論學科研究,還應結合實踐的積累。通常采用物理模擬(相似實驗)、邊坡監測(GNSS、雷達監測等)和數值模擬(有限元、有限差分等)3 種研究方法,結合理論分析及現場實際情況得出邊坡變形破壞機制和規律[9]。
利用有限元軟件,構建典型剖面數值模擬模型,相關物理力學參數見表1。采用位移邊界條件進行模型約束,固定x=0 和x=474 2 個面x 方向位移,固定y=0 和y=100 2 個面y 方向位移,固定底面z=0豎直方向位移,頂面為自由邊界條件,介質的彈塑性狀態采用理想的摩爾-庫倫彈塑性模型描述,用以保持受力體系的平衡。
對典型剖面進行數值模擬計算,高陡軟弱夾層邊坡水平方向位移云圖、總位移云圖、剪應力云圖及剪切應變增量云圖如圖2~圖5。

圖2 水平方向位移云圖

圖3 總位移云圖

圖4 剪應力云圖

圖5 剪切應變增量云圖
由圖2 和圖3 可知,邊坡呈現上部第四系沉降,下部沿弱層剪出,總位移最大約0.256 m,發生在第四系黏土層區域,水平最大位移量約0.218 m,發生在弱層上方+1 355~1 380 m 水平。由圖4 可知,弱層上部1 355~1 380 m 水平出現明顯應力集中,表現為拉應力,總體呈應力剪出狀。由圖5 可知,邊坡剪應變增量較大的位置主要集中在+1 355 m 水平弱層位置,尚未見貫通。綜上,南幫邊坡潛在滑移面是由切層巖體和軟弱夾層面組成的復合滑面,即南端幫邊坡的破壞模式主要為沿弱層滑動的“坐落-滑移式”。
選取極限平衡理論體系分析邊坡穩定性[10]。依據“露天煤礦工程設計規范”相關規定,礦區非工作幫臨時性邊坡穩定系數Fs取值范圍為1.00~1.20,根據礦區地質資料、邊坡邊坡現狀及南幫邊坡對露天礦的作用,本次選取安全儲備系數為1.10。根據所建立的外排土場典型工程地質剖面圖,巖土體物理力學指標推薦值見表1,利用Morgenstern-Price 法分析外排土場邊坡穩定性,穩定分析結果如圖6。

圖6 穩定分析結果
根據典型剖面邊坡穩定計算結果可知,邊坡穩定系數為1.01,邊坡處于極限平衡態,正在蠕動變化,易發生滑坡事故。
鑒于南幫上部1 345 m 水平有裂縫顯現的現狀,為防止出現坡頂裂隙與弱層貫通,造成邊坡滑動。考慮露天礦地質條件、礦區現狀、治理的迫切性、施工費用經濟性及施工安全性,選削坡減載治理為主,GNSS 監測為輔的治理設計方案。
治理工程涉及+1 320~+1 380 m 水平范圍,削坡減載治理方案設計主要原則為:在開采境界內,保障上部運輸道路完整,保障上部供電系統、其它設備位置不變,保障+1 350 m 水平運輸道路寬度不小于40 m,同時確保上部松散第四系黃土穩定。
以現狀邊坡特征為基礎,參考礦區初步設計,同時結合邊坡現場實際條件,選取單臺階角度為55°,設計邊坡整體南移10、15、20、25、30 m,治理設計曲線如圖7。

圖7 治理設計曲線
根據圖7 計算結果及南幫安全儲備系數要求,確定南幫最優削幫方案為南擴25 m,單坡面角55°,穩定系數1.110,滿足安全儲備系數要求。同時可知,該設計方案下南幫單剖面開挖量約455 m2,其中土方量273 m2,巖石182 m2。
為確保南幫邊坡的長治久安,保證礦區正常生產作息安全,需建立健全“礦山邊坡監測控制系統”,保證邊坡實時在線監測與控制,避免“為了監測而監測”或“只監不控”局面的繼續,故在邊坡變形區建立在線邊坡監控系統(GNSS 邊坡監控系統),對監測點進行數據采集分析,繪制速度監測成果,治理前位移曲線圖如圖8。

圖8 治理前位移曲線圖
由圖8 可知,變形區范圍內DLT-01~DLT-04監測點日均位移處于20.36~29.43 mm/d。未出現明顯加速階段,邊坡處于蠕動變形階段與現場實際情況相符。GNSS 監測系統可進行有效解算、處理與分析,能夠形成相應的位移數據曲線,為邊坡“臨滑預警”提供基礎數據支撐。
1)基于數值模擬分析可知,平朔東露天礦邊坡變形破壞機理為沿軟弱夾層的“坐落-滑移”式破壞模式。水平最大位移量約0.218 m,發生在弱層上方+1 355~1 380 m 水平。
2)基于現場勘查及極限平衡法可知,南幫現狀邊坡處于極限平衡狀態,并正在蠕動變形,亟需開展邊坡治理工程。
3)綜合考慮各種治理措施,針對南幫現狀提出:以削幫減載為主,GNSS 監測為輔。治理設計邊坡穩定系數為1.10,滿足安全儲備系數要求;GNSS 監測系統,監測現狀邊坡變形速率約為0.025 m/d,需要盡快推進治理方案。