束佳明 賈江鋒 張 臣 馬 寧
(山東科技大學礦業與安全工程學院,山東省青島市,266590)
隨著煤炭資源的不斷開發,礦井開采深度不斷增大,地質條件日趨復雜,礦井沖擊地壓等地質災害日益嚴重。根據現場實際生產情況及實測數據,在斷層、褶曲構造帶或煤層傾角變化帶等區域,易誘發沖擊地壓等地質災害。斷層作為礦井最為常見的地質構造,切割巖層的連續性,礦壓分布規律受斷層影響十分明顯。
現階段,針對斷層附近礦山壓力顯現特征及沖擊地壓的誘發,我國科學工作者已進行了大量的研究。代進等采用基本頂斷裂巖塊鉸接平衡理論,通過研究不同開采順序下,工作面采空區基本頂荷載在斷層兩側的傳遞規律和差異性,并結合實際現場,得出了斷層上盤先行開采時,斷層的載荷傳遞系數較大,上盤斷層煤柱壓力小,下盤煤體壓力較大;李志華等采用了斷層滑移失穩相似模擬試驗及數值模擬,研究了斷層附近工作面上下盤開采順序,通過對比分析得出工作面位于斷層下盤時沖擊礦壓危險性遠高于工作面位于上盤時的危險性。根據地質力學觀點,逆斷層構造附近更易發生沖擊地壓等動力災害,這是由于在水平方向的構造應力強烈擠壓作用下,斷層上下盤圍巖更容易發生彎曲,易積聚大量的彈性能。一旦斷層上下兩盤發生相對滑動將直接導致斷層帶圍巖積聚對的彈性能瞬間釋放,造成礦山傷亡事故。
本文以鮑店煤礦73上09工作面地質條件為基礎,通過對比不同尺寸的斷層煤柱,研究逆斷層下盤工作面開采時頂板垂直應力分布、煤壁下沉量及斷層的相對滑移量,進而分析誘導沖擊地壓的發生的規律。
UDEC是處理不連續介質的二維離散元程序,適用于模擬含節理或結構面不連續塊體集合問題,能夠滿足研究需要。以鮑店煤礦73上09工作面地質條件為基礎,針對逆斷層下盤工作面開采,建立UDEC數值計算模型。
考慮邊界效應,數值計算模型長550 m,高155 m,兩側各留設50 m邊界煤柱,采空區長度為120 m。數值計算模型如圖1所示。頂底板巖層物理力學參數見表1。斷層巖石物理力學參數如下:剪切模量1.52 MPa,體積模量2.30 MPa,容重2700 kg/m3,粘聚力0.05 MPa,內摩擦角16°,抗拉強度0.04 MPa。為減少模型運算時間,模型邊界附近網格劃分相對稀疏。

圖1 數值計算模型

巖性厚度/m剪切模量/GPa體積模量/GPa容重/kg·m-3粘聚力/MPa內摩擦角/(°)粉砂巖244.848.8227503.3030砂質泥巖101.002.1723601.3038粉砂巖124.848.8227503.3030細砂巖103.387.8727003.2628泥巖61.002.1723601.3038煤層82.173.6823602.1025粗砂巖83.36.8726903.1628砂質泥巖151.002.1723601.3038細砂巖123.387.8727003.2628粉砂巖304.848.8227503.3030
計算模型底部邊界采用全約束邊界條件,左右兩側邊界采用水平方向約束、豎直方向自由的邊界條件,頂端采用自由邊界條件。模擬工作面埋深400 m,模型頂板巖層高90 m,上部未模擬巖層自重以均布載荷垂直施加于模型頂端,因此補償載荷為7.75 MPa。
考慮到模型邊界的合理性,實際觀測表明,采空區的長度和寬度均達到和超過1.2~1.4H0(H0為平均開采深度)時,地表可達到充分采動。因此對于模型采空區長度為120 m時,地表達到充分采動。
在模型下盤煤層中分別設置應力監測線,監測工作面距斷層45 m、35 m、25 m、15 m和5 m時的應力狀態,得到工作面超前支承應力演化規律,煤柱垂直應力分布如圖2所示。
由圖2可以看出,當工作面距離斷層超過45 m時,工作面超前支承應力分布形式基本不變,表明斷層對工作面的支承應力影響不大。但由于工作面直接頂為硬厚粗砂巖,易形成大面積懸頂,工作面應力集中程度高。

圖2 下盤工作不同寬度斷層煤柱應力分布
當工作面距離斷層25~45 m時,受斷層切割煤體影響,斷層與工作面之間形成大煤柱,煤巖體應力集中明顯,并隨著工作面向斷層推進而逐漸升高,頂板與煤體均處在高應力狀態下,煤柱最大應力達到26.7 MPa(距離斷層25 m);應力峰值區間縮小,應力峰值與煤壁間的距離縮短。斷層下盤工作面前方煤巖體積聚了大量的彈性能,斷層煤柱區域工作面巷道發生高應力沖擊危險性較高。
下盤工作面距離斷層5~15 m時,工作面進入小煤柱支承階段,工作面的開采對斷層的構造應力的影響也達到最大,超前支承應力影響區轉移至上盤,下盤應力全部降低,高應力區轉移至上盤。由于小煤柱發生較為劇烈的塑性破壞,承載能力顯著降低,小煤柱應力降低,應力峰值降低6.5~18.95 MPa,峰值位置距煤壁位置也降低至2~7 m。此階段,斷層下盤煤柱若突然發生脆性破壞,極易誘發動載型沖擊災害。
在煤壁上方5 m處布置監測點,采用History命令,記錄工作面每次開挖后煤壁上方頂板下沉量的歷史值,得出其運動曲線圖,曲線的斜率反映頂板的下沉速率。提取下盤工作面每次開挖計算8000步的數據(運行5000步左右時已基本達到平衡狀態),繪制出下盤工作面開挖頂板運動曲線,如圖3所示。
由圖3可以看出,下盤工作面距斷層25~55 m時,受斷層切割影響,大煤柱支承狀態下,煤巖體處于高應力狀態,但大煤柱的支承能力強,煤壁上方的最大位移僅由1.38 m增至1.68 m,位移及下沉速率均較少,易誘發高應力型沖擊災害。當工作面距斷層15 m時,進入小煤柱護巷階段,煤柱所承受的上覆載荷減少,且承載能力顯著降低,煤柱破壞嚴重,此時的頂板位移激增至2.02 m,下沉速率明顯增大,頂板運動釋放較多能量,易誘發煤巖體整體失穩型沖擊災害。當工作面距離斷層5 m時,煤柱進一步破壞,下沉量仍較大,易誘發頂板沉降型沖擊災害。

圖3 下盤開采頂板運動曲線
為研究斷層面的滑移量,在平行于斷層面的上盤和下盤布置等高的監測點,采用History命令,記錄工作面每次開挖后監測點下沉量的歷史值,取斷層上盤和下盤的位移的差值,可作為斷層面滑移量,如圖4所示。

圖4 下盤開采斷層滑移量變化曲線
由圖4可知,斷層面的相對滑移量隨著工作面距斷層的距離增大而減小,斷層與工作面之間的距離越近,斷層面相對滑移量越大。當下盤工作面距斷層25~55m時,斷層的滑移量從0.416 m增加至0.495 m,斷層滑移速率及滑移量較小。工作面距斷層15 m時,斷層的滑移量與滑移速率均有明顯的變化,滑移增加0.05 m,煤柱受上覆巖層垮落影響而發生破壞,釋放較多的能量,易發生煤巖體整體失穩型沖擊災害。當工作面距離斷層5 m時,煤柱進一步破壞,此時的斷層滑移量也達到最大值0.61 m,易誘發頂板沉降型沖擊災害。
當工作面與斷層之間留煤柱超過15 m時,煤巖體處于高應力狀態下,頂板下沉量及其下沉速率、斷層滑移量均較小,易誘發高應力型沖擊災害;當工作面與斷層之間留設小煤柱時,煤巖體發生破壞,下盤應力降低,高應力區全部轉移至上盤,而頂板下沉量及其下沉速率、斷層滑移量較大,易誘發頂板沉降型沖擊災害。
(1)隨著開挖不斷接近斷層,斷層帶附近應力不斷集中,進入小煤柱支承階段,下盤工作面應力降低,高應力區轉移至上盤;
(2)頂板位移及其沉降速率、斷層滑移量隨煤柱的增大而減小,表明工作面距離斷層近,越易發生沖擊地壓;
(3)逆斷層下盤開采時,大煤柱支承階段,煤巖體處于高應力狀態,主要以應力型沖擊災害為主,可采取煤層注水、大直徑鉆孔卸壓等措施進行處理;進入小煤柱階段,煤柱承載能力降低,易誘發頂板沉降型沖擊災害,此階段應時刻監測礦壓數據,對端頭處加強支護護。