周永利，周 偉，陸 翔，才慶祥

(1.神華準能集團有限責任公司 科學技術研究院，內蒙古 準格爾 010300；2.中國礦業(yè)大學 礦業(yè)工程學院 煤炭資源與安全開采國家重點實驗室，江蘇 徐州 221116)

露天煤礦有限邊坡長度一般多在200m以下，對于完整巖體來說具有明顯的三維效應，三維邊坡結構對于邊坡穩(wěn)定性有著較大的影響。在巖體性質較弱或松散土石混合體邊坡一般采用易滑區(qū)煤炭回收技術回收煤炭，此時就要控制邊坡長度，通過回填和支撐形成擠壓壓力拱，防止邊坡滑坡和坍塌。邊坡太長就不易形成擠壓壓力拱而造成坍塌，局部坍塌會造成擠壓壓力拱失效，造成滑坡[1-4]。

我國蒙東地區(qū)10余個露天煤礦都存在明顯易滑區(qū)。在露天礦工作幫存在明顯弱面以及出現(xiàn)局部滑坡跡象時，根據(jù)地應力組合梁卸載模型，利用單斗挖掘機—卡車工藝工作線短的特點，采用“短工作線、高強度推進、快速回填”的開采技術，即采用短工作線來實現(xiàn)高強度推進，通過剝離物快速回填來及時加固邊坡，實現(xiàn)了露天采礦工程與邊坡治理一體化，安全地回采工作幫不穩(wěn)定區(qū)域的煤炭資源。

大型露天煤礦壓幫內排下端幫邊坡易滑區(qū)煤炭回收技術形成的基本特征是：①受采動影響明顯；②暴露面積變?。虎鄞嬖跁r間縮短。露天煤礦組合拱墻的特點：①由于露天煤礦采場工作幫、內排土場物料物理力學性質、采掘和排土工藝的差異，露天煤礦端幫組合拱墻具有典型非對稱性；②由于露天煤礦是沉積巖，各層拱墻之間具有粘結層，各層變形具有明顯協(xié)同性；③各層變形、破壞和失穩(wěn)的力學基礎是重力通過泊松效應轉化成的水平應力，因此，應力具有明顯的層次性。

我國大型近水平賦存露天煤礦在開采過程中，需要按照開采計劃進行剝離和采掘作業(yè)。原始地層的地應力經(jīng)過長期的卸載已經(jīng)處于平衡狀態(tài)，在露天開采過程中，由于開挖作業(yè)消除了原始地應力的有效支撐體，導致原巖應力二次卸載和重新分布，不同走向長度的邊坡，邊坡中的應力狀態(tài)不同，當邊坡走向長度較小時，應力會向兩側巖體中轉移，形成應力卸載拱，拱內的巖體處于受壓狀態(tài)，有較大的強度，邊坡的穩(wěn)定性較高[5]；如果邊坡走向長度較大時，應力重新分布后無法形成卸載拱，邊坡中部的巖體內會出現(xiàn)拉應力，由于巖體的抗拉強度很小，此時的邊坡穩(wěn)定性會降低。露天礦端幫邊坡采動過程中不平衡應力卸載對于邊坡整體穩(wěn)定十分不利，因此，研究采動過程中的卸荷規(guī)律，對于內排壓幫及端幫穩(wěn)定具有重要意義。

1 采動過程卸荷特點及應力模型

露天開采作業(yè)造成了原始地層的改變，破壞了地應力平衡，造成地應力的二次卸載。地應力卸荷會造成邊坡朝臨空面變形，破壞邊坡內部結構的完整性，但是，卸荷持續(xù)一段時間之后，巖質邊坡會形成自穩(wěn)結構，這種自穩(wěn)結構類似于自然界及工程中的多數(shù)穩(wěn)定結構，為類似“壓力拱”結構。相關研究表明，由于結構應力和自重應力的共同作用，在采掘區(qū)兩側形成壓力拱，且壓力拱的厚度和深度會隨著開挖深度的增大不斷增大。

壓力拱[6-8]的判斷方法主要有：①直接根據(jù)圍巖應力升高區(qū)來判斷；②基于壓力拱成拱系數(shù)來確定。這兩種判斷方法均是基于圍巖切向應力升高的特點而提出的。其中，基于成拱系數(shù)的判斷依據(jù)為：

式中，k為成拱系數(shù)；σθ、σθ0分別為開挖后、開挖前的切向應力，MPa。

當k<0，表明開挖引起切向應力下降；當k>0，表明開挖造成了切向應力升高，當k=0，則表明切向應力沒有變化。很多研究證明，壓力拱邊界的準確確定直接影響工程質量和作業(yè)環(huán)境的安全性。

根據(jù)彈塑性理論，求得極坐標中的應力分量關系式為：

式中，σx為x方向的正應力，kN；σy為y方向的正應力，kN；τxy為xy方向的剪應力，kN。

當θ=0°時，σθ=σy；當θ=90°時，σθ=σx，根據(jù)成拱系數(shù)的零值點，便可以圈定壓力拱范圍。

2 支撐拱模型

滑坡“支撐拱”的計算模型如圖1所示。在圖1(a)中，弧為滑坡主斷壁，弧CYD為滑坡破壞面趾，折線AM′MC，BN′ND分別為滑坡的兩翼。根據(jù)“支撐拱”模型實驗分析，當滑坡兩翼之間的土體達到塑性極限平衡狀態(tài)，且土體沿兩翼邊界面及滑面發(fā)生向下滑移時，M，N兩點之間土體的最大主應力軌跡將形成一條連續(xù)的拱線——最大主應力拱，亦即平臥支撐拱。拱圈厚度為b，拱頂至主斷壁距離為t，拱的跨度(亦即M，N兩點之間的距離)為l，拱的矢高為f。M，N兩點處的H，V分別為拱座水平推力和豎向反力。在拱座處，最大主應力σ1，作用平面與滑坡側翼面夾角為θ。在平臥支撐拱上方一定高度的M′，N′處，由于土體產(chǎn)生了向坡體下方的連續(xù)滑動，使得滑坡兩翼對土體的摩阻力得到充分的發(fā)揮，從而形成了拱頂凸向下方的最小主應力拱，在該拱座處，最大主應力σ1作用平面與滑坡側翼面(破裂面)夾角為θ′。

圖1 滑坡平臥“支撐拱”模型

材質均勻的滑動土體中，在平臥拱上方取距滑坡主斷壁y處的土條記為i，土條高為h，寬為dy，長為L，Wi是其自重；q與q+dq分別為平行于滑面方向作用于土條兩側單位長度單位高度的均布荷載；E及E+ΔE分別為作用于土條兩側的條間力合力，方向與土條底部平行；土條底部滑面傾角為α；Ni及Ti分別為作用于土條底部的總法向反力和切向阻力；Fl為土條在滑坡兩側的阻力。

若Pi為土條由于本身自重而產(chǎn)生的下滑力，則有：

Pi=Wisinα=rhLsinαdy

(5)

式中，γ為土體容重，N/m3。

根據(jù)莫爾-庫侖準則，土條在滑面上的平均切向阻力為：

Ti=Wicosαtanφs+csLdy=γhcos2αtanφsdy+csLdy

(6)

式中，cs，φs分別為滑面抗剪強度指標。

而土條在兩側的阻力為：

式中，cf為滑坡側翼面抗剪粘聚力，MPa；φf為滑坡側翼面摩擦角，(°)；K0，Kw為不同情況下土的側壓力系數(shù)；γhK0/2代表高度為h的土體作用在滑坡側翼面上的靜止土壓力，其中K0=l-sinφ；qKw代表平行于滑面方向的土壓力q對側翼面形成的水平壓力；Kw為當滑坡土體中最小主應力拱形成以后，由于應力傳遞而增大了滑坡側翼面上的正應力時，該處土的側壓力系數(shù)。同時，根據(jù)定義有：

E=qhL

(8)

假定滑坡土體為剛塑性體，則當土條i處于極限平衡狀態(tài)時：

E+Pi-(E+ΔE)-Ti-2Fi=0

(9)

求得實際作用在拱上方的均布載荷為：

由于拱結構和荷載的對稱性，圖1(a)表示的拱結構可簡化為圖2(a)所示的半拱結構。建立新坐標系xo′y，將坐標原點設在o′點。o′點只作用有水平軸向力T，均布荷載par與τar分別為拱由于自重而產(chǎn)生的平行于滑面方向單位面積的下滑力以及拱在滑面上單位面積的抗滑力。

par=γhcosαsinα

(11)

設平臥支撐拱拱軸為合理拱軸線，對于拱軸上—點m(x，y)建立力矩平衡方程，即：

式中，Sx為o′與m(x，y)兩點之間拱圈(圖2(b)中陰影部分)的面積。

計算求解后得有：

可見拱軸線為二次拋物線。

圖2 半拱結構分析

露天礦采場邊坡大多數(shù)都是有限邊坡，近水平露天煤礦實現(xiàn)內排，端幫邊坡就是采場一側的支撐和內排土場一側的回填動態(tài)變化的過程，此時端幫邊坡如果受到破壞，其經(jīng)歷兩個過程：完整巖體形成的支撐拱提供抗滑力，拱中的應力超過巖體的強度時，巖體裂隙開始擴張，開始擴容，此時擴容作用又形成相互擠壓的抗滑拱，兩個拱形成共同作用的組合拱；排土場邊坡是松散土石混合體組成的有限邊坡，松散土石混合體相互擠壓形成抗滑拱，同時，巖體的內摩擦力也作為抗滑力，在長時間堆載作用下松散土石混合體重塑，物理力學性質有所強化，提高了邊坡的抗滑力。

3 支撐-回填有限邊坡數(shù)值模擬

為了研究露天開采時的壓力拱效應，借助數(shù)值模擬工具，對易滑區(qū)掏槽開采進行研究，確定其自穩(wěn)結構，揭示成拱規(guī)律，為易滑區(qū)的陡幫開采和資源回收，提供技術保障。端幫邊坡開挖以后受到地應力的卸載作用，邊坡不斷變形，假設在露天礦邊坡開采的范圍以內邊坡的地應力是恒定的，且方向一致。露天煤礦邊坡大多是層狀沉積巖，各層之間交結，在滑動體內部也形成了層狀梁，端幫邊坡地應力卸載可以視為非對稱組合拱墻的卸載，采場和內排土場作為兩個固支點。

FLAC3D數(shù)值模擬軟件是通過有限差分法來計算模型的應力狀態(tài)，該軟件在巖土工程問題中是明顯優(yōu)于極限平衡法的，它能夠從三維角度對邊坡進行分析，能夠更真實的反映邊坡的穩(wěn)定性狀態(tài)。小龍?zhí)兜V務局布沼壩露天煤礦西端幫屬于典型的易滑邊坡，按照原始設計方案執(zhí)行開采作業(yè)，形成的端幫邊坡穩(wěn)定性僅能達到1.1左右。該邊坡若受到地下水或震動的影響，極易形成大規(guī)模的滑坡。為了回收端幫煤炭資源，采用條帶式掏槽靠幫開采，按照采掘槽底部寬度D為50m、100m、150m的順序進行數(shù)值模擬，得到采掘槽兩側的應力分布表現(xiàn)為明顯的“壓力拱”結構，具體如圖3所示。研究分析的西幫邊坡高度為180m，垂直邊界范圍是邊坡高度的4～5倍，設置單元數(shù)72695個，網(wǎng)格節(jié)點數(shù)89560個。

圖3 掏槽開采時壓力拱分布

由圖3可知，在端幫易滑區(qū)執(zhí)行條帶式掏槽靠幫開采時，邊坡地應力朝采掘槽中心卸載，并表現(xiàn)為標準的“壓力拱”形狀，其壓力拱具有以下基本特征：①隨著采掘槽寬度的不斷增大，壓力拱的半徑也不斷增大，當采掘槽寬度增大到一定程度之后，壓力拱效用僅作用在采掘兩側的邊坡上，在采掘槽底部的應力效應變??；②壓力拱具有明顯的時間效應，隨著掏槽寬度及推進位置的發(fā)展順序，壓力拱集中于長期存在的邊坡，且存在時間越長，壓力拱對其作用越明顯，因此，縮短邊坡暴露時間能降低地應力卸載形成的壓力拱在邊坡上的滑移驅動效應；③壓力拱具有非對稱性，對于開采形成的對稱采掘槽，壓力拱效用基本對稱，對著開采作業(yè)沿端幫推進，打破了原有的平衡，工作幫一側的壓力拱效用較小，非工作幫一側的壓力拱效用較強；④當采掘槽延伸到最終位置時，剝離物沿非工作幫一側開始內排，此時非工作幫轉變?yōu)閮扰磐翀龉ぷ鲙?，內排作業(yè)消弱了壓力拱效用，并促使整個壓力拱結構沿推進方向向前移動，始終包圍在采掘槽兩側；⑤當采掘槽的深度不同時，壓力拱的形狀也有差異；采掘槽越深，壓力拱的拱頂曲率越大，反之越小。

在掏槽開采作業(yè)過程中，應力集中區(qū)域的分布具有明顯特征，其集中位置與采掘槽尺寸，開挖先后順序密切相關，具體表現(xiàn)為以下幾個規(guī)律性特征：①尺寸對稱的采掘槽，即工作幫與非工作幫邊坡角度相同時，兩側應力集中效果相同，應力集中的程度與采掘槽邊坡角度有關，角度越大，應力越大，反之較小；也與邊坡高度成正比，高度越大，應力越大，反之越?。虎诜菍ΨQ采掘槽，最大應力集中在邊坡角度大的一側坡頂線位置，邊坡角度小的一側，應力集中效應較弱；③開采擾動對于應力集中位置同樣產(chǎn)生一定的影響，當采掘槽左右對稱時，工作幫一側受開采作業(yè)的擾動，成為整個邊坡的應力集中區(qū)域；④應力集中區(qū)域與邊坡存在(暴露)時間有關，存在時間越長，應力集中效應越明顯，存在時間越短，應力集中效應越小，主要是地應力卸載時間不同造成的，采用內排壓幫，壓覆暴露一段時間的端幫，可以有效的阻止端幫變形的繼續(xù)發(fā)展，并迫使應力集中區(qū)域迅速轉移到內排臺階表面。

從應力集中區(qū)域的變化規(guī)律可知：掏槽開采作業(yè)破壞了端幫的完整性，造成了不平衡地應力的再次卸載，地應力卸載時間的長短會直接影響應力在邊坡表面的分布情況及邊坡巖體的蠕變位移，卸載時間越長，應力和蠕變位移越大。采掘槽工作幫不斷推進對于邊坡造成的擾動使最大應力主要集中在采掘槽工作幫一側的邊坡坡頂線位置，最大應力會隨著采掘槽的推進不斷向前移動。在采掘槽非工作幫壓幫內排之后，內排臺階阻止巖石蠕變的繼續(xù)發(fā)展，并提供與地應力卸載方向相反的反作用力，此時非工作幫一側的應力集中區(qū)域便向前移動至內排臺階表面。由此可見，隨著采排作業(yè)的動態(tài)發(fā)展，應力、應變都在不斷的發(fā)生變化，并與采掘槽推進速度之間存在一定的匹配關系。

小龍?zhí)兜V物局布沼壩露天煤礦西端幫屬于多斷層、軟巖端幫，若受到地下水或爆破作業(yè)的影響，極易發(fā)生滑坡。該礦端幫邊坡只能實行條帶式掏槽靠幫開采，掏槽寬度D按50m、100m、150m的順序遞增，對不同掏槽寬度時的邊坡穩(wěn)定性進行模擬分析，結果如圖4—圖6所示。

圖4 布沼壩露天礦西端幫邊坡穩(wěn)定性分析結果

圖5 掏槽開采時邊坡穩(wěn)定性分析結果

圖6 掏槽開采50m壓幫時邊坡穩(wěn)定性分析結果

從圖5可以看出，布沼壩露天礦西端幫邊坡穩(wěn)定隨著掏槽寬度的增加，成比例遞減，二者的線性關系明顯，當采用內排壓幫時，邊坡穩(wěn)定系數(shù)出現(xiàn)大幅度提高，增加的幅度達到5%左右。

5 結 論

1)掏槽寬度與邊坡穩(wěn)定性呈現(xiàn)負相關性，采掘槽寬度越大，易滑區(qū)的暴露面積越大，不平衡地應力的有效卸載空間越大，巖體發(fā)生蠕變的時間越長，整體結構松弛度越大，穩(wěn)定狀態(tài)逐漸惡化；

2)壓幫內排及時跟進可以提高端幫邊坡的穩(wěn)定性，其主要機理是內排臺階提供了有效支撐力，阻止端幫邊坡的蠕變發(fā)展，限制了地應力的卸載空間，掩埋了變形較大的坡腳位置，并使應力集中區(qū)域轉移到內排臺階表面。

3)對于不同巖性，不同尺寸的邊坡，實行掏槽式靠幫開采時，都存在可執(zhí)行的最大掏槽，對于掏槽開采方案的可執(zhí)行性，要結合礦區(qū)既有采掘設備的具體參數(shù)來綜合論證。