基于UDEC的水庫下壓煤開采安全性分析

張　進

(國投晉城能源投資有限公司)

?

基于UDEC的水庫下壓煤開采安全性分析

張進

(國投晉城能源投資有限公司)

摘要為分析裴溝煤礦魔洞王水庫下壓煤開采的安全性，利用UDEC數值模擬軟件建立了走向回采工作面的開采模型，模擬計算得出不同開采厚度條件下上覆巖層破壞發展過程及對導水裂隙帶發育高度的影響，并與經驗公式計算的結果對比分析，驗證數值模擬的正確性。根據導水裂隙帶發育高度和地表裂隙發育情況，分析得出魔洞王水庫下開采安全防水巖柱尺寸，結果表明，魔洞王水庫下回采是安全可行的。首采的31071工作面已經安全通過水庫區，從而驗證了分析正確性。

關鍵詞UDEC覆巖破壞導水裂隙帶高度

影響水體下壓煤開采安全性的因素主要是導水裂隙帶發育高度，當保安煤柱高度不夠或留設不合理時，有可能發生水體或者泥沙導通井下或工作面涌水量超限等問題[1]，嚴重影響礦井的安全生產工作。目前國內外對裂隙帶發育高度的研究方法包括理論分析、現場鉆(物)探實測、相似模擬分析和經驗公式計算等[2]。隨著科學技術的發展，數值模擬分析以其計算快捷方便，能夠直觀顯示覆巖破壞的發展變化過程而應用越來越廣泛。利用UDEC建立數值模型，分析魔洞王水庫下采煤覆巖破壞特征和導水裂隙帶發育狀況，對該礦區進行水庫下采煤安全性分析和合理選擇采高具有重要的參考意義。

1工程概況

裴溝煤礦始建于1960年，是鄭煤集團的主力生產礦井。目前礦井開采的31采區位于井田東部，主采的二1煤層平均厚約7.8 m。31采區地表的魔洞王水庫始建于上世紀50年代，為小(二)型，水庫平均深約3.13 m，水域面積為600.6萬m2，覆蓋31采區中北部大片區域，壓煤量預計達500余萬t。31采區共布置7個工作面，大多數工作面要穿過水庫。工作面與水庫相對位置關系見圖1。

2數值模擬

UDEC是針對處理不連續介質的二維離散元程序，該程序基于離散單元法DEM開發，模型內劃分的塊體單元及其節理面可以被壓縮、擠壓而產生滑動、變形、分離等破壞，能夠形象地顯示節理巖體的非線性大變形特征。

圖1　31采區工作面布置與魔洞王水庫相對位置關系

2.1模型建立

計算模型根據裴溝煤礦的實際采礦地質條件，并考慮計算模型的邊界效應對于模擬開采區域的影響，確定模型尺寸為1 000 m×300 m(長×高)，見圖2。判斷巖體模型破壞采用mohr-coulomb屈服準則(式1)和Hoek-Brown破壞準則(式2)。

圖2　開采模型

(1)

(2)

式中，σ1，σ3分別為最大、最小主應力，MPa；c為黏聚力，MPa；φ為內摩擦角，(°)；σci為單軸抗壓強度，MPa；m，s分別為反映巖石的軟硬程度和破碎程度的經驗參數。

2.2模型力學參數

根據裴溝礦31采區地質鉆孔柱狀圖資料確定模型各上覆巖層及其層位關系。參考該礦實測巖石力學參數，依據Hoek-Brown破壞準則確定計算模型巖體及節理的物理力學參數[3]，見表1、表2。

表1　巖體力學參數

表2　節理力學參數

2.3模擬結果分析

計算模型均以20m為計算步距，分別模擬采高為3,5和7.5 m時上覆巖層的破壞狀況。用豎向位移云圖表示覆巖垮落和裂隙發育，根據劃分的離散元塊體垮落、離層情況，判斷覆巖發展破壞和裂隙發育狀況[4],見圖3。

圖3　不同采高時上覆巖層堅向位移變化云圖　　由圖3可知：

(1)當采高為3 m時，工作面從80 m推進到100 m，覆巖破壞高度基本保持在42 m；當采高為5 m 時，工作面從100 m推進到140 m，覆巖破壞高度基本保持在55 m；而當采高為7.5 m時，工作面從140 m推進到180 m，覆巖破壞高度基本保持在74 m。這表明采高分別為3，5，7.5 m時，導水裂隙帶發育高度分別為42，55，74 m。這也說明覆巖破壞有時間和空間的發展過程，在覆巖充分破壞以前，裂隙發育隨著工作面的推進不斷向上發展，而當覆巖充分破壞以后，裂隙發育高度基本保持不變。

(2)不同采高條件下覆巖破壞形態均由上凸形逐漸向馬鞍形過渡，且在覆巖破壞高度達到充分后，馬鞍形的破壞形態更加明顯。這是由于隨著工作面的推進，覆巖先后發生垮落破壞，當回采過的空間達到一定程度后，上部巖層不再破斷而產生彎曲下沉，并作用在已垮落的巖體上，下部已垮落巖體受上部巖層的壓實作用，覆巖破壞過程趨于穩定，部分裂隙逐漸閉合，裂隙高度略有下降。而在工作面端部，受超前支撐壓力作用，覆巖破壞劇烈，裂隙發育高度較大，從而在形態上表現為裂隙發育中間低兩端高的馬鞍形[5]。

(3)根據裂隙高度與采高的關系，采高為3 m時，裂采比(裂隙帶發育高度與采高之比)為14；采高為5 m時，裂采比為11；采高為7.5 m時，裂采比為9.87，即隨著采高的增加，裂隙帶發育高度與采高之比逐漸減小。

3水庫下壓煤回采安全性分析

3.1導水裂隙帶發育高度計算

分析31采區內鉆孔柱狀圖可知，該采區上覆巖層組成主要為中粒砂巖、細粒砂巖、粉砂巖、砂質泥巖和泥巖等。綜合分析各巖層組分，砂巖、粉砂巖、泥巖組成比例約為0.34∶0.08∶0.58，覆巖類型為軟弱偏中硬。

利用經驗公式(表3)計算采高分別為3，5和7.5 m時導水裂隙帶發育高度，結果見表4。

表3　緩傾斜厚煤層分層開采導水裂隙帶計算公式

注：ΣM為累計采厚，ΣM<15 m；M為單層采厚，M<3 m；±為中誤差。

表4　導水裂隙帶高度計算結果

3.2模擬結果與經驗公式計算結果對比

對比分析不同采高時數值模擬與經驗公式計算得出的導水裂隙帶發育高度，可以看出，在采高為3和5 m時，數值模擬與經驗公式計算結果非常接近，而采高為7.5 m時，數值模擬計算值高于經驗公式計算值。由于經驗公式是前人基于以往分層開采和普采的統計資料分析確定，對于綜采工作面有一定的局限性。實測和理論分析都表明，隨著采高的增大，綜采工作面覆巖破壞強度和裂隙發育高度與分層開采和普采相比明顯增大。這說明利用UDEC數值模擬研究覆巖破壞特征比利用經驗公式具有更廣泛的適用性，分析計算的結果也更真實可靠。

3.3防水巖柱高度分析

防水煤巖柱指開采煤層頂板至地面水體底界之間的巖層整體。其防隔水性能與該段巖柱的含水性、隔水性及抵抗采動裂隙產生的能力有關，防水煤巖柱的安全高度是判斷水體下采煤是否安全的重要指標。防水煤巖柱高度應大于導水裂隙帶與保護層厚度之和，計算公式為

(3)

式中,Hsh為防水巖柱高度,m；Hli為導水裂隙帶高度,m；Hb為保護層厚度,m；h為地表裂縫深度,m。

利用建立完備的開采沉陷預計模型[6]，推導得出回采后表土層裂縫臨界水平變形值和裂縫發育極限深度的估算公式：

(4)

(5)

式中,εJ為計算裂縫臨界水平變形值,mm;h為計算裂縫極限深度,m;c為黏聚力,MPa；φ為內摩擦角,(°);μ為泊松比;E為彈性模量,GPa；γ為巖土層容重,kN/m3。

根據31采區表土層的物理力學性質[7]，E=13.1 MPa,c=0.029 MPa,μ=0.3,φ=18°,γ=15.8 kN/m3,經計算，回采后地表裂縫臨界水平變形值可達5.5 mm，裂縫發育深度約3.7 m，考慮安全因素，裂縫發育深度取6 m。

依據裴溝礦軟弱偏中硬的巖層性質，保護層厚度Hb取一次開采厚度的6倍[8]，分析31采區鉆孔柱狀圖可知，煤層埋深在水庫上部約204 m，中部約252 m，下部約394 m。不同采高時防水巖柱高度與上覆巖層厚度的對比關系見表5。

由表5可知，開采煤層上覆巖層的最小厚度均遠大于需要留設的防水巖柱的高度，且水庫底有約18 m厚的第四系黏土作為隔水層，表明魔洞王水庫下采煤是安全可行的。

表5　防水巖柱與煤層上覆巖層厚度對比關系

431071工作面試采

目前31071試采工作面已經回采完畢，采高為5 m，通過分析井下和地表的觀測資料，在回采期間工作面沒有出現涌水量超限現象，地表水庫水體也沒有受到回采的影響，這表明分析結果是準確的、可靠的。

5結論

根據裴溝煤礦的覆巖賦存特點，利用UDEC數值計算軟件模擬分析得出不同回采厚度對上覆巖層破壞和導水裂隙帶發育高度的影響，并實現了覆巖發展破壞過程的可視化；對比分析數值模擬結果和經驗公式計算結果，表明UDEC數值模擬具有更廣泛的適用性和準確度；分析得出魔洞王水庫下采煤防水煤巖柱尺寸，結合31采區的巖層賦存特征，表明魔洞王水庫下開采是安全可行的；首采的31071工作面已安全通過水庫庫區，從而驗證了分析是正確的、可靠的。

參考文獻

[1]康永華，靳仁昌.水體下放頂煤開采研究現狀及其發展趨勢[J].煤礦開采,2003,8(1):15-18.

[2]盛佳，李向東.基于Hoek-Brown強度準則的巖體力學參數確定方法[J].采礦技術,2009,9(2):12-14.

[3]曾先貴，李文平，李洪亮等.綜放開采近斷層導水裂隙帶發育規律研究[J].采礦與安全工程學報,2006,23(3):306-310.

[4]張玉軍，李鳳明.高強度綜放開采采動覆巖破壞高度及列席發育演化監測分析[J].巖石力學與工程學報,2011,30(S1):2994-3001.

[5]國家煤炭工業局.建筑物、水體、鐵路及主要井巷煤柱留設與壓煤開采規程[M].北京：煤炭工業出版社,2000.

[6]吳侃,黃珍珍,王欣.礦山開采沉陷的完備預計模型[J].煤礦開采,2008,11(4):4-6.

[7]王金啟.裴溝煤礦水體下放頂煤開采安全性分析[J].中州煤炭,2008(1):18-19.

[8]陳俊杰，郭文兵，鄒友峰.大型水體下頂水安全開采的可行性研究[J].中國安全科學學報，2011,21(2):57-62

(收稿日期2015-11-23)

張進(1987—)，男，助理工程師，碩士，048026 山西省晉城市鳳臺東街1188號。