急傾斜煤層深部開采對淺部斷層活化影響的數值模擬研究

摘要：淮南礦區某礦采區屬倒轉地層，灰巖含水層位于煤層頂板，斷層的活化極易導通上部含水層，引發突水事故。面對急傾斜煤層開采深度不斷增加的現狀，采用FLAC3D數值模擬的方法，對急傾斜煤層深部開采對頂板淺部斷層活化影響的規律進行了研究，總結了急傾斜煤層開采深度與斷層活化之間的時空關系，驗證了通過控制采高和隔離煤柱大小，可以有效降低急傾斜煤層開采對斷層的擾動作用，對急傾斜煤層后期深部開采斷層突水預測與防治具有重要的理論指導意義和工程應用價值。

關鍵詞：急傾斜煤層；深部開采；數值模擬；斷層活化

中圖分類號：TD163 文獻標識碼：A 文章編號：1674-098X（2015）04（a）-0000-00

Abstract： A certain ore mining area of Huainan mining area belongs to inversion formation， limestone aquifer located seam roof， the activation of fault easily conduct the upper aquifer， leading to water inrush disaster. Facing the continuously increased mining depth status of steep seam， using FLAC3D numerical simulation method， In steep seam mining deep effect on the fault activation of roof is studied， Summarizes the temporal relationship between the depth of steep seam and fault activation， verified by controlling the mining height and isolated coal pillar size， can effectively reduce the disturbance of steep seam faults， It has important theoretical guidance significance and engineering application value to the later stage of steep seam deep mining fault water inrush prediction and prevention.

Keywords： steep seam； deep mining； numerical； fault activation

中圖分類號：TD745

我國煤田水文地質條件復雜，在煤礦生產過程中，礦井突水事故頻發。據資料統計，其中80%的煤礦突水事故是由斷層導水引發的[1]。斷層活化是指斷層在原始不導水的狀態下，受地下工程的擾動作用，導水系數增大而成為導水通道。斷裂構造作為一種重要的突水通道，研究分析斷層在采動影響下的活化反應，對煤礦的安全生產具有重要意義。對于斷層活化的研究，學者們普遍認為，當斷層面的切應力τ大于斷層面的抗剪強度τn時，即τ≥τn時，斷層就會沿斷層面發生破裂，上下盤相對移動形成導水通道[2]。姜耀東等[3]把斷層接觸面以庫倫剪切模型代替，研究了采動過程中斷層接觸面上應力時空演化規律和斷層兩盤的運動規律；卜萬奎等[4]從斷層傾角變化的角度研究了斷層活化的規律。

雖然目前的研究已經對斷層活化的機理和規律有了一定的認識[5]，但是對于如何預測和監測斷層活化仍有很大的困難。急傾斜煤層的賦存特征和開采方式決定了其在開采過程中對斷層活化的影響與水平煤層具有較大差別[6]。急傾斜煤層進行淺部開采時，擾動作用強度和影響范圍有限，不具備誘發附近斷層活化的能力；但是隨著開采深度的增加，不平衡力積累到一定程度，就會導致淺部斷層活化，導通含水層，引發斷層突水事故。本文采用數值模擬軟件FLAC3D，以急傾斜煤層為研究對象，對煤層深部開采誘發淺部斷層活化進行了模擬和預測，并提出了預防和控制斷層活化的后期開采方案。

1FLAC3D計算原理簡介

FLAC3D采用了拉格朗日連續介質法（Fast Lagrangian Analysis of Continua），同時還采用了混合離散法和動態松弛法，可以模擬和記錄地質體在達到屈服極限時發生的破壞和塑性流動情況，分析漸進式破壞和模塊失穩，并且適合于模擬大變形；其根據已知應變增量，可以方便地求得應力增量、不平衡力并跟蹤系統的演化過程；由于計算過程中不會形成剛度矩陣，所以運算所需內存很小。除此之外，FLAC3D還具有內嵌的程序語言FISH，使用戶可以自行開發系統，增強了程序的靈活性，充分擴展了軟件的自身功能。[7]

煤層在采動過程中，破壞了周圍巖石的應力平衡，在重力和地應力的作用下，巖石向采空區移動，在不發生破壞的情況之下，相鄰巖石的位移變化相差不大，而斷層帶作為軟弱接觸面，當采動造成的不平衡力達到臨界值時，斷層就會在巖石移動過程中沿斷層面發生斷裂，使斷層上、下兩盤發生錯動和滑移，增大了斷層上、下兩盤的位移差，即斷層進入活化狀態。本文基于FLAC3D強大的數值模擬和監測記錄功能，監測在急傾斜煤層開采到不同深部時斷層上、下盤的相對位移量，研究預測在不同煤層開采強度下的斷層活化情況。

2計算模型建立

2.1 工程概況

淮南礦區某礦為二疊系急傾斜煤層，地面標高+20m。受推覆構造作用，煤田地層發生倒轉，總體呈傾伏倒轉褶曲形態，主要可采煤層地層傾角為50°～80°，煤層頂底板巖石主要為砂巖、泥巖和砂質泥巖，第四系厚度140m。倒轉翼灰巖含水層成為煤層頂板，可采煤層均位于倒轉褶曲帶下翼，灰巖不僅上覆于煤系地層上方，而且是礦井主要水害威脅含水層。該礦急傾斜煤層自上而下分水平進行開采，每一水平的段高為30米，水平之間留設5m的隔離煤柱，開采方法為偽傾斜柔性掩護支架采煤法，沿偽傾斜布置采煤工作面，工作面偽傾斜角度為30°。受區域推覆構造影響，研究區域小斷層較發育，其中有正斷層位于煤層頂板，切割倒轉翼太原組灰巖含水層。距離煤層約220m。煤層目前已開采至-480m水平，繼續進行深部開采可能會引發該斷層的活化，導通太原組灰巖含水層，引發突水事故。本文選擇了該煤層底板至太原組灰巖，埋深140m～740m的地層范圍，建立了包含斷層的三維地質模型。

2.2 模型建立及參數選擇

數值模擬過程中共建立了三組模型，通過模擬不同的工作面段高和隔離煤柱大小，研究煤層深部采動下對淺部斷層活化的影響。三組模型的工作面參數分別為：工作面段高30m，保留5m煤柱；工作面段高30m，保留6m煤柱；工作面段高25m，保留5m煤柱。三維模型根據研究區域范圍內的勘探線地質剖面圖建立。模型尺寸為x=600m，y=200m，z=600m，最大模擬深度740m，每組模型的煤層厚度3m，東西走向，煤層傾角70°。斷層落差10m，方位角205°，傾角70°，延伸到-130m～-480m。模型頂部壓力用等效載荷代替，頂部為第四系地層，取平均密度 ，平均厚度 ，則 （負號代表方向向下）。三維模型如圖1所示。

結合地層綜合柱狀圖，為了簡化網格剖分的工作量，同時充分體現各層之間的力學性質的差異，將巖性、力學性質和分布厚度相近或相同的巖層劃歸為一個層組，合并后巖層自上到下可劃為13個層組，根據該煤礦巖石力學試驗資料以及鄰近煤礦的資料，獲取巖層物理力學參數，如表1所示。假設巖石各向同性、均質且符合摩爾-庫倫（Mohr—Coulomb）彈塑性模型介質，模型底部限制垂直移動，前后和左右限制水平方向位移。模型共劃分了201560個單元，233751個節點。

數值模型選取沿煤層走向的中線y=100m處取xz平面垂直切片，作為數值模擬研究的分析斷面，為了監測整個斷層的活化情況，在斷層上、下盤以垂直深度間隔50m選取8對監測點（見圖2），通過記錄監測點的位移量，計算每對監測點的位移差，來研究斷層的活化規律。

3數值計算與分析

3.1 模擬過程

按照煤礦實際開采情況，模擬開采標高-200m至-700m的煤層，監測和記錄煤層往深部開采過程中每個監測點的水平位移量和垂直位移量，計算斷層上盤與下盤的位移差，繪制位移差曲線，以此來判斷斷層的活化情況。

為了方便后文敘述，圖1中標明了斷層帶附近的監測點在模型中的位置及標高。

3.2 斷層活化規律分析

1）工作面段高30m，保留5m煤柱

圖2為煤層采動過程中斷層兩盤觀測點的水平方向和垂直方向位移差曲線圖。由圖2可以看出，煤層進行淺部開采時，在水平方向和垂直方向上斷層上、下盤位移差很小且無明顯變化，基本未受到采動影響，當煤層開采至-570m標高時，不同位置的觀測點曲線開始在水平方向和垂直方向同時出現分化，其中-130m至-230m的觀測點仍保持較低的變化率，而深部-230m以下的觀測點水平方向位移差和垂直方向位移差出現快速增長，說明斷層底部位置最先發生活化反應，并隨著煤層采深的增加斷層兩盤距離不斷增大，煤層開采結束后，水平方向最大位移差為28.28mm，垂直方向最大位移差為47.10mm，根據斷層傾角計算，位移差法向分量為42.39mm，垂直法相分量為34.94mm，斷層破碎帶的拉破壞大于剪切破壞。

通過計算可以看出，急傾斜煤層深部開采對于斷層活化的影響程度顯著大于淺部煤層開采。計算結果表明，當煤層開采至-570m標高時，斷層帶同時出現剪切破壞和拉破壞，斷層發生活化，而且隨著采深的增加活化程度不斷加劇。在實際生產中，煤層預計深部開采將超過-700m標高，斷層作為安全隱患極易發生活化引發突水事故。通過降低煤層回采率，來減弱對斷層的擾動作用，需要增加隔離煤柱的高度和減小工作面的高度。因此在原模型的基礎上，降低煤層回采率，驗證上述兩種方法是否能起到控制斷層活化的作用。

2）工作面段高30m，保留6m煤柱

圖3為工作面高度保持30m不變，隔離煤柱高度增加到6m后斷層兩盤觀測點的水平方向和垂直方向位移差曲線圖。由圖3可以看出，當煤層開采至-600m標高時，斷層發生活化反應，煤層開采結束后，水平方向最大位移差為26.04mm，垂直方向最大位移差為29.35mm，斷層的活化反應速度在時空分布上均有所降低。

2）工作面段高25m，保留5m煤柱

圖4為隔離煤柱高度保持5m不變，工作面高度減小到25m后斷層兩盤觀測點的水平方向和垂直方向位移差曲線圖。由圖4可以看出，當煤層最早開采至-670m標高時，斷層發生活化反應，煤層開采結束后，水平方向最大位移差為7.88mm，垂直方向最大位移差為12.69mm，活化程度微弱，斷層的活化反應速度在時空分布上明顯降低。

因此，通過增加隔離煤柱的高度以及減小工作面的高度，可以減弱斷層的活化程度，延遲斷層的活化時間，對控制斷層的活化具有明顯效果。

4結論

1）在急傾斜煤層開采過程中，隨著開采深度的增加，雖然工作面距離淺部斷層越來越遠，但是由于持續的擾動作用和不平衡力的積累，當煤層開采至-570m標高時，斷層帶開始出現剪切破壞和拉破壞，斷層發生活化反應，并隨著開采深度的增加活化程度不斷加劇，斷層破碎帶上拉破壞大于剪切破壞。

2）通過增大隔離煤柱高度和減小工作面高度可以有效的減小急傾斜煤層深部開采對淺部斷層的擾動作用，因此煤礦可以綜合考慮兩種方法，在保證煤層回采率的基礎上，提前做好頂板探放水工作，確保煤礦的安全生產。

參考文獻

[1]宋振騏，郝建，湯建泉，石永奎. 斷層突水預測控制理論研究[J]. 煤炭學報，2013，09：1511-1515.

[2]許進鵬，張福成，桂輝，張同俊. 采動斷層活化導水特征分析與實驗研究[J]. 中國礦業大學學報，2012，03：415-419.

[3]姜耀東，王濤，趙毅鑫，王辰. 采動影響下斷層活化規律的數值模擬研究[J]. 中國礦業大學學報，2013，01：1-5.

[4]卜萬奎，茅獻彪. 斷層傾角對斷層活化及底板突水的影響研究[J]. 巖石力學與工程學報，2009，02：386-394.

[5]錢鳴高， 石平五. 礦山壓力與巖層控制[M]. 徐州：中國礦業大學出版社，2003：107-115.

[6]李棲鳳. 急傾斜煤層開采[M]. 北京：煤炭工業出版社，1984.

[7]魏學勇，武強，趙樹賢. FLAC3D在礦井防治水中的應用[J]. 煤炭學報，2004，06：704-707.