呂維赟,路 凱,任卓鑫,連 強
(1.潞安化工集團,山西 長治 046000;2.山西潞安環(huán)保能源開發(fā)股份有限公司 常村煤礦,山西 長治 046000;3.中國礦業(yè)大學 礦業(yè)工程學院,江蘇 徐州 221116;4.潞安化工集團 安全管理部,山西 長治 046000)
切頂卸壓圍巖控制是改善圍巖應力、煤礦實現(xiàn)可持續(xù)發(fā)展的重要技術之一。切頂卸壓圍巖控制技術被廣泛應用于各個礦井中,該技術改善了煤礦的圍巖應力環(huán)境,效果良好。近年來,許多專家學者針對該技術進行了大量研究工作。何滿潮院士及其團隊[1-3]提出以切頂短臂梁理論為基礎的切頂卸壓沿空留巷技術,在國內深部沿空留巷切頂卸壓圍巖控制中應用的效果較好,受到普遍認可;王炯等[4]通過相似材料模擬試驗,對比分析了切頂與非切頂情況下巷道圍巖變形規(guī)律; 鄭立軍[5]等人在古漢山礦1604綜放工作面沿空留巷中采用定向預裂爆破技術后,爆破效果較好,圍巖控制效果頗佳。王莊煤礦91 采區(qū)具有強度較低的泥巖直接頂和厚度較大的底板,給巷道的掘進和圍巖控制帶來了巨大的困難;此外,當開采的深度越來越大,其中的瓦斯含量也越來越多,增大了煤與瓦斯突出的風險。針對91-101 工作面的實際開采情況,現(xiàn)有的頂板處理技術還有許多不足:成本高、施工難度大、切頂裝備復雜,難以適應礦井復雜多變的工作環(huán)境,且切頂質量難以達到留巷要求。為此,以王莊煤礦91-101 工作面風巷為工程背景,通過力學分析、數(shù)值模擬,研究新型切頂卸壓沿空留巷圍巖豎向位移變化規(guī)律、應力及塑性區(qū)的分布特征,揭示切頂高度對圍巖變形影響規(guī)律[6-20]。
91-101 工作面位于91 采區(qū),東接540 運輸大巷,南為91-201 掘進工作面,北西為未采區(qū)域。地面標高為930~932 m,工作面標高為339~408 m,可采長度2 465 m。工作面回采可能會引起地面塌陷,對北康線與北栗至古城礦支線、693 縣道造成不同程度的破壞。工作面布置圖如圖1。
圖1 工作面布置圖Fig.1 Working face layout
91-101 工作面運巷長2 686 m,風巷長2 680 m,91-101 工作面風運巷可采長度均為2 465 m;工作面切眼長260 m,煤層總厚2.12 m,密度為1.42 t/m3,一次采全高,在91-101 工作面回風巷進行沿空留巷。
預裂爆破斷頂沿空留巷即采用定向預裂爆破技術切斷巷道頂板懸梁結構,減小懸臂長度,達到卸壓和留巷的目的。根據(jù)斷頂沿空留巷上覆巖體的結構特點,可以將模型進行簡化,簡化后的力學模型如圖2。l1段是預留巷道實體煤側,取l1長度為應力峰值區(qū)寬度的2 倍;b 段為預留巷道段,巷旁支撐載荷等效為集中力F1,左邊界固定約束等效為彎矩M1和力F2作用。
圖2 斷頂沿空留巷頂板懸梁力學計算模型Fig.2 Mechanical calculation model of roof suspension beam in goaf retaining with broken roof
由計算模型進行力學分析并計算,結合邊界和連續(xù)條件可以求得斷頂后預留巷道采空側巷旁位移量ω(x)為:
式中:E 為巖石彈性模型,GPa;I 為截面慣性矩,m4;q1為上覆巖體載荷;q2為實體煤支撐載荷邊值;q3為實體煤支撐荷載峰值。
比較斷頂和未斷頂?shù)捻敯鍛伊毫W模型,可得出未斷頂與斷頂?shù)南锏啦煽諅任灰屏恐睢鳓?為:
式中:l2為巷道寬度,m。
代入現(xiàn)場具體數(shù)值,可得出切頂后巷道采空側巷旁位移量比未斷頂時減小了46.7%。
工作面后方上覆巖層達到強度極限,在受到壓力后,巖層會垮落并在采空區(qū)形成堆積。針對91-101 工作面,上部巖層垮落是不規(guī)則的,經重新壓實后松散系數(shù)較小,針對巖石殘余破碎,殘余破碎膨脹系數(shù)Kp取1.1。
針對直接頂,取巖層的垮落厚度為Σh,則KpΣh便是堆積高度?!鳛槠渑c基本頂之間空隙,其值為:
式中:M 為回采高度,取1.6;Kp為巖石的殘余碎膨脹系數(shù),取1.1。
當M=Σh(Kp-1)時,△=0,采空區(qū)內填充滿垮落的直接頂巖層,支撐覆巖載荷。因此,若需要承載覆巖荷載,垮落的直接頂高度為:
將參數(shù)代入式(3)可得為19.2 m,即91-101 工作面風巷沿空留巷理論上需要斷頂垮落高度。
合理的切頂爆破應在考慮經濟適用性的情況下減小留巷圍巖應力水平,從而大幅減小巷旁支護阻力,減小頂板回轉變形,因此用數(shù)值模擬對切頂爆破過程中的切頂高度、切縫寬度、切頂角度等參數(shù)進行研究,揭示切頂沿空留巷圍巖變形規(guī)律。
為研究91-101 綜放工作面沿空留巷超前預裂爆破技術及參數(shù)優(yōu)化,利用FLAC3D模擬軟件建立數(shù)值計算模型。模型計算遵循莫爾-庫倫屈服準則,變形模式為大變形,上覆施加載荷為10 MPa,數(shù)值模型如圖3。此模型采用如下的邊界條件:周圍是滑動的邊界條件,模型的下面是固定的邊界條件,模型的上部施加相似模擬的巖層質量。工作面回采過后,在巷內支護和鋼管混凝土作用下共同維護巷道穩(wěn)定。
圖3 切頂爆破數(shù)值計算模型Fig.3 Numerical calculation model of roof cutting blasting
將切頂高度作為研究因素分析,研究切頂留巷后巷道的位移、應力以及塑性區(qū)的變化規(guī)律,針對91-101 工作面采區(qū),得出切頂高度的最優(yōu)解。
針對不同情況討論切頂爆破高度對留巷后巷道的影響,對這幾種情況下的模型進行圍巖應力、位移以及塑性區(qū)的分析,研究預裂切頂高度對圍巖變形的影響,并確定最佳預裂切頂高度。建立的幾種模型唯一變量為切頂爆破高度,分別取為17、18、19、20 m,其他因素不變,暫將切頂爆破角度定為0°,切頂爆破寬度為1 m。工作面處切頂后豎直位移圖如圖4,頂?shù)装逦灰屏咳鐖D5。
圖4 工作面處切頂后豎直位移圖Fig.4 Vertical displacement diagrams of working face after roof cutting
圖5 頂?shù)装逦灰屏縁ig.5 Displacement of roof and floor
結合4 個工況的位移云圖可知,切頂爆破高度分別為17、18、19、20 m 時,工作面巷道頂?shù)装逡平糠謩e為0.052、0.048、0.040、0.038 m。結合圖4、圖5 分析發(fā)現(xiàn),切頂爆破高度從17 m 增加至19 m 的時候,頂?shù)装逡平恐饾u減??;在切頂爆破高度由19 m 增加至20 m 的時候,巷道頂?shù)装逡平肯嗖畈淮螅f明在19 m 處,切頂爆破卸壓效果最佳,再增加切頂爆破高度不再增加卸壓效果。綜合分析對比4 種情況的位移云圖并結合經濟適用性,發(fā)現(xiàn)切頂爆破高度為19 m 時,卸壓效果最佳。
工作面處切頂后應力圖如圖6。
圖6 工作面處切頂后應力圖Fig.6 Stress diagrams after roof cutting at working face
結合4 種情況的應力云圖可知,經過切頂爆破卸載壓力后,回采工作面處應力重新分布,相較于周圍應力較小。切頂爆破高度為17、18、19 m 時,巷道圍巖應力以及墩柱的垂直應力隨著切頂高度的增加而減小,說明從17 m 到18 m 再到19 m,切頂爆破卸壓效果越來越好;切頂爆破高度為19 m 和20 m時,巷道圍巖應力以及墩柱的垂直應力相近,說明切頂爆破高度從19 m 往后卸壓效果不再大幅增加。綜合分析對比4 個情況的應力云圖并結合經濟適用性,發(fā)現(xiàn)切頂爆破高度為19 m 時,卸壓效果最佳。
工作面處切頂后塑性區(qū)圖如圖7。
圖7 工作面處切頂后塑性區(qū)圖Fig.7 Plastic zone diagrams after roof cutting at working face
結合4 個情況的塑性區(qū)云圖可知,切頂爆破高度分別為17、18、19、20 m 時,巷道周圍頂板及兩幫均受到大小不同的剪切和拉伸破壞。當切頂爆破高度從17 m 增加至19 m,巷道周圍的破壞程度逐漸減??;切頂爆破高度為19 m 和20 m 的時候,巷道周圍的破壞程度相近,說明切頂爆破高度為19 m時卸壓效果趨于最佳值,繼續(xù)增加切頂爆破高度,卸壓效果已不會大幅增加。綜合分析對比4 個工況的塑性區(qū)云圖并結合經濟適用性,發(fā)現(xiàn)切頂爆破高度為19 m 時,卸壓效果最佳。
針對工作面頂板下沉量大、巷道斷面收斂大等問題,需進行預裂切頂卸壓,考慮技術和經濟因素,對各項參數(shù)進行確定,切頂參數(shù)見表1。
表1 切頂參數(shù)表Table 1 Parameters of top cutting
現(xiàn)場設置3 個測量站點,每個站點用深基點位移計對圍巖的頂?shù)装逦灰埔约皟蓭臀灰七M行觀測測量。留巷期間巷道圍巖表面位移變化圖如圖8。
圖8 留巷期間巷道圍巖表面位移變化圖Fig.8 Surface displacement variation of roadway surrounding rock during roadway retention
由圖8 可知,在3#站點測量的頂?shù)装逦灰屏恳约皟蓭臀灰屏烤畲蟆?0 d 之后移近量基本不再變化。整體來看,3 個測站巷道頂?shù)装逦灰啤蓭臀灰谱兓口厔莼鞠嗤?,變形量相差很小,圍巖控制效果好。
1)通過建立力學模型,代入現(xiàn)場數(shù)據(jù)計算,對比未斷頂和斷頂時頂板的懸梁結構,得到沿空留巷采空側位移量表達式,切頂后巷道采空側位移比未斷頂時減小了46.7%。
2)采用數(shù)值模擬的方法計算得出,針對91-101工作面開采情況切頂高度為19 m 時效果最佳。
3)確定切頂?shù)母黜梾?shù),根據(jù)各項參數(shù)進行現(xiàn)場試驗,從分析出的數(shù)據(jù)可知圍巖控制效果好、可行。