巖溶煤礦部分開采合理采高的研究*

陽志成，朱川曲

(湖南科技大學 資源環境與安全工程學院， 湖南 湘潭市 411201)

巖溶煤礦部分開采合理采高的研究*

陽志成，朱川曲

(湖南科技大學 資源環境與安全工程學院， 湖南 湘潭市 411201)

為了探究巖溶煤礦部分開采方法的合理采高，防止巖溶塌陷，以豐城礦務局某煤礦的地質條件和生產條件為背景，采用FLAC3D數值模擬軟件，模擬出了不同采高時上覆溶洞周圍塑性區的發育情況，確定了礦井的最合理采高以及巖溶塌陷易發性的評價指標。研究結果表明，地下開采會使巖溶地層和表層產生拉應力，形成拉伸破壞區，對于淺埋溶洞，其周圍拉伸破壞區易發展至表層破壞區而導致巖溶塌陷。由數值模擬結果可知，充分開采后上覆溶洞周圍塑性區演化至非常接近貫通但又不貫通至表層塑性區時的采高為1.4 m，因此，該礦最合理采高為1.4 m。通過對數值模擬監測數據的計算分析及地表移動的計算分析可知，地表水平拉伸變形量大于1.89 mm/m時，地表可能發生巖溶塌陷。研究結論可為該巖溶煤礦安全、經濟、綠色開采及相關研究提供依據。

煤礦；合理采高；水平拉伸量；FLAC3D；巖溶塌陷

0 引 言

我國的巖溶面積大約有346萬km2，占全國總面積的近1/3，是世界上巖溶最發育的國家之一[1]。巖溶地區進行煤炭開采易造成巖溶塌陷，據研究，巖溶塌陷主要包括[2￣6]：重力塌陷、潛蝕塌陷、真空吸蝕塌陷、沖爆塌陷、振動塌陷、荷載塌陷、溶蝕塌陷、滲壓致塌陷等。長期以來，在研究采煤誘發巖溶地表塌陷機理時，人們普遍認為，在巖溶地層受到采動影響時，產生的裂隙與導水裂隙帶聯通，是采煤引起煤礦巖溶塌陷的重要原因，但湯伏全、胡炳南[7￣8]等在相關研究中發現，一些巖溶煤礦垮落帶和導水裂縫帶最大的高度相對于上覆巖層的厚度來講是很小的，采空區不可能直接溝通至含水地層引起地表塌陷，而采動垂直離層和水平拉伸裂縫可能是巖溶塌陷的基礎原因；王明立[9]在研究中發現地表拉伸變形量的數值達到1.8 mm/m時可能發生巖溶塌陷。采用部分開采、留煤柱開采、條帶開采及充填開采等開采方法都能有效減少和控制這兩種巖溶塌陷，基于部分開采相對其他幾種采煤方法具有經濟性和簡單有效性等優點，研究部分開采合理的采高，對巖溶煤礦的安全開采、防止巖溶塌陷及水資源保護有著重要意義。

1 數值模擬

1.1 礦井概況

該礦位于江西省境內，年生產能力約55萬t，地面平均標高為+27 m。平均煤層厚度為2.5 m。平均開采深度520 m。某工作面開采后，地表發生巖溶塌陷，塌陷坑的直徑達20 m，塌陷深度達5.0 m。根據地質報告，發生塌陷處的采區地層從上往下依次為松散層、長興巖層、粒砂巖層、細砂巖層、粉砂巖層、泥巖層、煤層及細砂巖層，共8層。表層松散層厚20 m，長興巖層厚140 m，該巖層中巖溶分布較廣，溶洞數目眾多、大小不一且深淺不同，造成該處塌陷的溶洞為此巖層中的一個淺埋溶洞。長興巖層有4組巖層共計厚度為450 m，前3層自上而下分別是粒砂巖、細砂巖、粉砂巖，而煤層處于長興巖層下的第四層巖石中，厚度為2.5 m，巖土力學參數值見表1。

1.2 模型建立

為簡化計算，對分析系統作如下假定：分析中不考慮巖溶堆積物的影響，按空洞進行計算；溶洞為漫長時期形成的，在煤炭開挖前，整個模型處于自然平衡狀態；開采期間，不考慮溶洞的尺寸等的物理力學性質的變化；模型為連續和均勻、各向同性的彈塑性體；不考慮水的作用、地震荷載、采動震動等對溶洞穩定性的影響。

表1 巖土力學參數取值

劉之葵、劉鐵雄等對巖溶空間形態的統計表明，在溶洞尺度較大(溶洞洞徑大于15 m)時，一般溶洞空間形態近似大廳狀，在溶洞洞徑不大時，溶洞斷面一般為圓形或似圓形[10￣11]。因此，將造成上述塌陷的溶洞的形狀設置成圓形。再根據最大影響角正切為1.89，設計模型寬度為1100 m。因此，建立模型長、寬、高分別為1100，1，542.5 m(見圖1)。

1.3 模擬方案

(1) 模擬過程。建立模型，本構關系和材料特性(本次模擬為摩爾-庫侖材料)；用不同采高開挖煤層，得出充分采動后溶洞周圍塑性區剛好聯通地表時的采高；提取采高最合理時的監測點數據，并對模擬數據進行分析。

圖1 數值模擬模型

(2) 采高的確定。模擬以充分采動后溶洞周圍塑性區演化至非常接近貫又不貫通至表層塑性區時的采高為最合理的采高。利用二分法來找出合理的采高，即從中點出發確定含合理采高的上下限，不斷縮小含合理采高的采高范圍直至模擬出最合理采高。

(3) 監測點與監測內容的確定。在模型頂部X軸方向的中間線上，沿X軸方向每隔110 m取一個監測點，兩端加中間共計11個監測點，監測監測點X方向上的位移(見圖2)。

圖2 監測點布置

1.4 模擬過程及結果分析

1.4.1 模擬過程

(1) 采高取煤層厚度的一半即1.25 m時，發現充分采動后塑性區未演化至表層塑性區，模擬結果見圖3。

(2) 采高為1.25 m時充分采動后塑性區未演化至表層塑性區，因此采高取1.25 m與2.5 m的中間值即1.875 m，繼續模擬時發現未達到充分采動時溶洞周圍塑性區已演化至表層塑性區，模擬結果見圖4。

(3) 采高為1.875 m時溶洞周圍塑性區已演化至表層塑性區，因此采高取1.25 m至1.875 m的中間值即1.5625 m，受模擬精度的影響，此時采高取1.55 m，繼續模擬時發現未達到采動時塑性區已演化至表層塑性區，模擬結果見圖5。

(4) 采高為1.55 m時溶洞周圍塑性區已演化至表層塑性區，因此采高取1.25 m至1.55 m的中間值即1.4 m，繼續模擬時發現充分采動后塑性區演化至非常接近貫通而又未貫通至表層塑性區，模擬結果見圖6。

(5) 確定最合理采高在1.4 m至1.55 m之間，取1.45 m時模擬發現未充分采動時塑性區演化已貫通至表層塑性區(見圖7)，因此，確定采高1.4 m為最合理采高。

圖3 采高為1.25 m時的塑性區演化結果

圖4 采高為1.875 m時的塑性區演化結果

圖5 采高為1.55 m時的塑性區演化結果

圖6 采高為1.4 m時的塑性區演化結果

圖7 采高為1.45 m時的塑性區演化結果

1.4.2 模擬結果分析

合理采高即采高為1.4 m時，充分采動后監測點的水平移動數據見表2。

采用圖2所示的坐標系，擬合出地表水平位移公式如下：

u(x)=6.394×10-9x3-1.047×10-5x2+0.003824x+0.003876

擬合度0.9853，符合精度要求。再由上式得地表水平變形公式為：

ε(x)=u'(x)=1.9182×10-8x2-2.094×10-5x+0.003824

計算得采高為1.4 m時溶洞處地表水平拉伸量為1.89 mm/m。

模擬結果表明，地下開采會使巖溶地層和表層產生拉應力，受拉應力的作用，溶洞周圍產生了拉伸破壞區，表層也產生了拉伸破壞區，二者一旦貫通，可能造成巖溶塌陷。部分開采能有效控制地表水平變形量，當采高控制在1.4 m以內時，溶洞周圍破壞區與表層破壞區不會貫通，采礦相對安全。

表2 采高為1.4 m充分采動后監測點的水平位移量

2 地表水平變形預測

采用概率積分法[13￣14]預測巖溶塌陷坑位置的水平移動變形，公式如下：

式中，Wcm=mqcosα；r為等價計算工作面的主要影響半徑；m是煤層開采的法向高度1.4m；q為下沉系數0.6；α為開采煤層的傾角10°；L為等價計算工作面各邊界的直線段；Hd為等價開采影響深度520m；β為主要影響角，其正切取值1.89；η、ζ為積分變量；Ucm= bWcm，b為平移動系數0.3；θ0為開采影響傳播角85°；ε、i分別為水平變形、傾斜變形。計算出部分開采采高為1.4m時的溶洞處地表水平拉伸變形量為1.06mm/m。

可見，預測公式計算比數值模擬計算的水平拉伸變形量要小，如果以數值模擬的結果來作參考，則預測公式中采高為1.4m時溶洞周圍的塑性區不會演化至表層破壞區，采高取1.4m也是安全的。由于現場實際的地層結構復雜，采用數值模擬和地表移動變形的預測理論來計算地表變形情況可能都會有所偏差。因此，綜合二者的結果，確定該礦該地區采高為1.4m，水平拉伸變形量大于1.89mm/m時巖溶塌陷發生的可能性較大。

3 結 論

(1) 部分開采采煤法能有效控制巖溶煤礦巖溶層和表層的拉伸破壞，設計合理的采高能防止巖溶塌陷，從而實現安全、綠色開采。實驗礦井合理的采高為1.4m。

(2) 地表拉伸變形量大于1.89mm/m時，溶洞周圍塑性區演化至表層破壞區而導致地表塌陷的可能性比較大，為防止采動引起的巖溶塌陷，采取合理的措施是必要的。

[1]袁道先.新形勢下我國巖溶研究面臨的機遇和挑戰[J].中國巖溶,2009,28(4):330￣331.

[2]康彥仁.論巖溶塌陷形成的致塌模式[J].水文地質工程地質,1992,19(4):32￣34.

[3]王 濱,李治廣,董 昕，等.巖溶塌陷的致塌力學模型研究[J].自然災害學報,2011,20(4):119￣125．

[4]張麗芬,曾夏生,姚運生，等.我國巖溶塌陷研究綜述[J].中國地質災害與防治學報,2007,18(3):126￣130．

[5]王衛軍,趙延林,李青鋒，等.礦井巖溶突水災變機理[J].煤炭學報,2010,35(3):443￣448．

[6]程 星,彭世壽.巖溶區地下水位下降致塌的數值模擬研究[J].地球與環境,2005,33(S1):119 ￣123．

[7]湯伏全,田志洲.采動影響下松藻礦區近地表溶洞的變形破壞規律[J].西安礦業學院,1993(3):249￣250.

[8]胡炳南,張慎勇,陳佩佩.江西坪湖礦控制巖溶塌陷的壓煤開采優化設計研究與實踐[C]//2010全國"三下"采煤與土地復墾學術會議,2010:11￣12.

[9]王明立.煤礦開采誘發巖溶地表塌陷的機理分析[J].煤炭科學技術,2014,42(8):5￣6.

[10]劉之葵,梁金城,朱壽增，等.巖溶區含溶洞巖石地基穩定性分析[J].巖土工程學報，2003,25(5):629￣633.

[11]劉鐵雄.巖溶頂板與樁基作用機理分析與模擬試驗研究[D].長沙:湖南大學,2004.

[12]高延法,張順峰，等.覆巖導水裂縫與巖層拉伸變形量的關系研究[J].采礦與安全工程學報,2012,29(3):303￣304.

[13]武 雄,高明顯,黃萬鵬，等.棗莊巖溶塌陷形成條件及主要影響因素[J].地學前緣,2007,14(6):227￣233.

[14]張華興.井筒煤柱開采理論與技術研究[D].北京:煤炭科學研究總院,1991.

國家自然科學基金項目(51474104).

2016￣11￣07)

陽志成(1987-)，男，湖南攸縣人，在讀碩士，主要從事保水開采方向的研究，Email：286105181@qq.com。