陳道志,楊素改,楊 松
(1.滕州郭莊礦業有限責任公司 錦丘煤礦, 山東 滕州市 277500;2.山東省天安礦業集團有限公司, 山東 曲阜市 273155)
煤層的開采必然會引起采場周圍巖層的變形、破斷和移動等現象,如果存在含水層,并形成導水通道,將會引發水災害。因此,在進行水體下壓煤開采時,就必須對煤層采動形成的裂隙是否相互連通以及相互連通的裂隙是否波及到水體等關鍵問題進行深入探討,為煤礦的安全生產提供保障[1-4]。
導水裂隙帶高度受到開采方法、覆巖巖性、時間等諸多因素的綜合影響[5],其值的確定便是一個復雜的問題。本文通過數值模擬計算得出導水裂隙帶高度,并通過現場實測進行對比分析,為煤礦安全生產提供一定的保障。
在實驗室內,將研究區域范圍內的煤系巖層根據相似原理用相似材料制成模型,或利用數值分析軟件建立數值模型。在滿足相似的初始條件、邊界條件等前提下,“開采”模型中的“煤層”,在相應的時間里模擬出相似的礦山壓力和礦山壓力顯現現象,推斷實際巖層發生的相應情況。
由于物理模擬需要大量相似的物質材料及繁瑣的工序,所以數值模擬顯得更為簡捷實用。通常使用的有離散元法、有限元法、三維快速拉格朗日法等。其中,三維快速拉格朗日法即FLAC3D法是一種基于三維顯示有限差分法的數值分析方法,它可以模擬巖土的三維力學特性,分析漸進破壞和失穩等,可以用來模擬采場大變形[6-7],為受采動影響圍巖變形破壞分析提供了便利。
理論計算法主要是解析法和數值法,兩者都以固體力學為基礎;經驗公式法是在分析確定了礦區覆巖結構類型的基礎上,根據《建筑物、水體、鐵路及主要井巷煤柱留設與壓煤開采規程》中給出的經驗公式來計算導水裂隙帶發育高度。
1.3.1 鉆孔分段注(放)水法
在采煤工作面周圍選擇合適的觀測場所,向采空區上方打仰斜鉆孔。鉆孔避開冒落帶、斜穿裂縫帶,達到預計的裂縫帶頂界后向上延伸一定高度。使用“鉆孔雙端封堵測漏裝置”系統,沿鉆孔分段封堵注水,由注水控制臺控制水壓、風壓,讀取注水流量。逐段觀測,直到測出整個鉆孔內各段的漏失量。根據漏失量變化情況確定裂縫帶的上界高度。
該方法具有精度高、工程量少、觀測效果明顯等優點,且觀測設備結構簡單、操作方便,易于現場應用。
1.3.2 鉆孔超聲成像探測法
不同的巖層具有不同的波阻抗,對聲波的反射能力也不同。通過向鉆孔孔壁發射超聲波脈沖并接收反射波,測量回波幅度和傳播時間,并由計算機處理成聲波幅度、時差圖像。聲波幅度圖像反映孔壁巖性硬度、巖石的結構,聲波時差圖像反映孔徑的大小。圖像顏色由深到淺分為16等級色,代表不同的硬度和孔徑。綜合分析圖像,可以看出裂縫的產狀及被充填的情況,即裂縫的開度或充填物性質[8]。
這種方法所獲得的資料比較直觀,同鉆孔沖洗液法相配合,可以很好地揭示巖層破壞特征,判定導水裂縫帶高度。
此外,還有鉆孔沖洗液耗量法[9]、三維地震技術[10]等實測方法。
某礦2702工作面位于-415 m水平,工作面按傾斜長壁布置,該面構造形態整體為一單斜構造,工作面煤層走向60°~90°,傾向北西,真傾角3°~10°。煤層平均煤厚1.25 m,厚度穩定,煤層可采指數為1.0,變異系數為10%。直接頂為泥巖,厚度5~7 m,硬度系數f=3;直接底為粉砂巖,厚度為1.1~1.2 m,硬度系數f=3。采用傾斜長壁后退式采煤法,全部垮落法管理頂板。采用綜合機械化采煤,工作面平均長度為160 m,傾斜長度為788 m,平均采高為1.2 m,循環進尺為0.6 m。
通過對上述地質條件進行分析,并考慮到現實情況,決定采用FLAC3D數值軟件對該面導水裂隙帶發育高度進行數值模擬,并利用鉆孔分段注水法進行實測。
2.1.1 數值模型的建立
利用FLAC3D軟件對自開切眼至充分采動全過程覆巖隨工作面推進時的變形、垮落情況進行模擬。模型尺寸為走向方向250 m,鉛直方向100 m,傾向200 m。模型中共有248750個單元,265251個節點,采空區部分用空單元模擬。在回采過程中,工作面每推進10 m進行一次運算,共開挖160 m。計算模型見圖1。

圖1 計算模型
2.1.2 數值模擬結果分析
(1) 應力場分析。工作面推進160 m時采場圍巖應力云圖見圖2。由圖2可知,受采動影響,重新分布的應力場呈拋物線形狀分布,且開切眼和煤壁前方出現明顯的應力集中現象。
(2) 位移分析。由圖3所示的工作面推進160 m時上覆巖層垂直位移云圖可見,采空區覆巖下沉特點表現為:下位巖層比上位巖層下沉值大,越往上,下沉值越小,下位巖層表現為垮落特性,向上依次表現為裂縫帶、彎沉帶;采空區中部比采空區邊緣處下沉值大,且下沉量隨與工作面煤柱的距離的減小而減小。

圖2 工作面推進160 m時垂直應力云圖

圖3 工作面推進160 m時垂直位移云圖
用數值模擬來判斷導水裂隙帶高度的判據有塑性區和應力兩種。應力判據是通過數值計算,得出各點的應力,根據該點的應力狀態來判斷其是否破壞[5]。據此,結合以上模擬結果分析,根據實驗室得出的該礦巖層力學參數,以各點應力狀態為判據,得出該礦覆巖導水裂隙帶高度為18.84 m。
本次現場實測選用鉆孔雙回路注(放)水法來實現?,F場共布置3個孔,其中采前孔一個,傾角70°,采后孔兩個,傾角分別為60°和70°。
2.2.1 鉆孔參數的確定
根據2702工作面的圍巖條件及巖石強度分級,確定選用中硬條件下裂隙帶高度的計算公式:
(1)
式中,M為采高,此處為1.2 m,代入(1)式計算得:H=16.1~27.3 m。因此,在考慮富余量的前提下,確定采前孔斜長為37 m,采后孔斜長分別為34 m和37 m。
2.2.2 實測結果分析
根據觀測數據,得出如圖4、圖5所示的采前孔、采后孔注水漏失量對比結果。

圖4 采后1#孔與采前孔漏失量對比

圖5 采后2#孔與采前孔漏失量對比
從圖4可以看出:在18.65 m位置處,采后1#孔漏失量曲線與采前孔漏失量曲線相交,當垂直深度超過18.65 m時,其漏失量與采前孔的基本一致,說明了采后1#孔的垂直深度18.65 m位置為該處導水裂隙帶發育高度。從圖5可以看出:采后2#孔與采前孔漏失量曲線在18.84 m位置處相交,在垂深大于18.84 m時,兩鉆孔漏失量基本吻合,說明了采后2#鉆孔的垂直深度18.84 m位置為該處導水裂隙帶發育高度。
通過兩個采后孔注水漏失量與采前孔注水漏失量的對比分析,可以確定該礦2702工作面覆巖導水裂隙帶發育高度為18.84 m。
通過數值模擬,得出某礦2702工作面的導水裂隙帶高度為18.84 m;通過現場實測,確定該面導水裂隙帶高度為18.84 m,兩者結果一致,表明數值模擬結果可靠,為覆巖導水裂隙帶發育高度的確定提供了簡單可行的方法。
參考文獻:
[1]黃炳香,劉長友,許家林.采動覆巖破斷裂隙的貫通度研究[J].中國礦業大學學報,2010,39(1):45-49.
[2]欒元重,李靜濤,班訓海,等.近距煤層開采覆巖導水裂隙帶高度觀測研究[J].采礦與安全工程學報,2010,27(1):139-142
[3]劉天泉.礦山巖體采動影響與控制工程學及其應用[J].煤炭學報,1995(2):1-5.
[4]張平松,胡雄武,劉盛東.采煤面覆巖破壞動態測試模擬研究[J].巖石力學與工程學報,2011,30(1):78-83.
[5]尹增德.采動覆巖破壞特征及其應用研究[D].青島:山東科技大學,2007.
[6]彭文斌.FLAC3D實用教程[M].北京:機械工業出版社,2009.
[7]陳育民,徐鼎平.FLAC/FLAC3D基礎與工程實例[M].北京:中國水利水電出版社,2008.
[8]韓緒山,張景考.聲波掃描成像在災害防治及環境保護中的應用[J].中國煤田地質,2003,15(2):59-62.
[9]張玉軍,李鳳明.高強度綜放開采采動覆巖破壞高度及裂隙發育演化監測分析[J].巖石力學與工程學報,2011,30(S1):2994-3001.
[10]潘冬明,程久龍,李德春,等.利用三維地震技術探測覆巖變形破壞研究[J].采礦與安全工程學報,2010,27(4):590-594.