基于Midas GTS的地下開采對管狀皮帶影響的有限元分析

秦韶芬，胡海峰，張雪芹(太原理工大學礦業工程學院，山西 太原 030024)

?

基于Midas GTS的地下開采對管狀皮帶影響的有限元分析

秦韶芬，胡海峰，張雪芹
(太原理工大學礦業工程學院，山西 太原 030024)

隨著煤炭資源的大規模開采，越來越多的建(構)筑物不可避免的受到沉陷影響，這對建(構)筑物的安全構成了極大威脅。為保證鎮城底礦22620工作面上方興能電廠管狀皮帶的正常運轉，以鎮城底礦地質采礦條件為基礎，通過Midas GTS巖土結構有限元分析軟件，對管狀皮帶在地下開采過程的移動和變形進行了動態預計，并與現場實測數據作了對比。研究結果表明，所建立數值模型滿足工程精度需要，工作面推進至310m時，管狀皮帶的移動變形值大于規程規定的臨界變形值，皮帶將受到采動破壞。根據預計結果提出了皮帶走廊的修補保護措施，以保證電廠正常運營。

地下開采；數值模擬；移動變形；Midas GTS；管狀皮帶

地下煤炭資源采出后破壞了巖體內部原有的力學平衡狀態，使巖層發生移動、變形，當開采規模超過一定范圍，破壞和移動將會擴展到地表形成地表沉陷，對地表的建(構)筑物產生重大影響[1]。為了達到開采煤炭資源和保護地表建(構)筑物的平衡，有必要對地下開采引起的地表移動變形進行預測。從目前國內外的研究來看，預測方法主要有唯象法、力學法、數值模擬和相似模擬4類[2-6]。

近些年來隨著計算機和軟件的不斷發展，數值分析方法的不斷成熟，數值模擬逐漸成為預測開采沉陷問題的重要手段[7]，目前廣泛使用的數值方法主要有：有限差分法、離散單元法、彈塑性有限單元法等[2]，各類數值分析軟件也運用到開采沉陷中。Midas GTS是一種較新的大型巖土結構有限元分析軟件，具有直觀建模、自動劃分有限元網絡等特點，能提供強大的二維和三維模擬功能，為解決地表沉陷問題提供了有利工具[8]。

本文通過采用Midas GTS有限元分析系統，結合鎮城底22620工作面地質條件，對工作面上方興能電廠排灰管道在采動過程中的穩定性進行了分析，并據此提出了皮帶保護措施。

1 工程概況

1.1 工程背景

古交興能電廠位于鎮城底礦南六采區，電廠皮帶排灰管道從采區上方通過，如果在其范圍內全部留設保護煤柱，勢必造成煤炭資源的嚴重浪費。為了盡可能開采煤炭資源并保證興能電廠排灰通道的正常運行[9]，必須分析地下開采對排灰管道穩定性的動態影響，以便及時采取保護措施，保證皮帶走廊的安全。

1.2 工作面概況

22620工作面位于元家山村以東、王家坡村以南、雞兒溝西北部；井下位置屬于南六采區。工作面走向長1207m，傾向長170m，平均開采深度430m。所采煤層為2#、3#煤，煤層厚度3.5m，整體呈一向斜構造，煤層平均傾角為6°。工作面于2011年8月開始回采，開采為走向長臂式開采，頂板管理為全部垮落法。22620工作面上方建有地表巖移觀測站。觀測站共布置傾斜觀測線2條，走向觀測線1條，管狀皮帶觀測線1條，其中皮帶觀測線位置見圖1。

圖1 皮帶觀測線布置圖

2 模型建立

2.1 Midas GTS簡介

Midas GTS(Geotechnical and Tunnel analysis System)是將有限元分析與巖土結構的專業性有機結合的有限元分析軟件，具有基于CAD的三維幾何建模功能、自動劃分網格、映射網格等特點，其強大的前、后處理能力及分析能力為解決開采沉陷問題提供了優秀的平臺[10]。本次模擬結合開采實際對開采過程進行模擬，從而預測地面及管狀皮帶的動態移動與變形。

2.2 力學模型的選擇及巖層力學參數的選取

覆巖破碎的抗剪強度一般遵循Mohr-Coulomb準則，所以本次模擬采取摩爾-庫倫(Mohr-Coulomb)模型。建模時按力學性質特點對不同地層進行一些合并，模型所包含的巖層從上到下依次為：黃土、砂質泥巖、砂巖、砂質泥巖、灰巖、泥巖、煤層、細砂巖。所采用的巖層力學參數見表1。

2.3 幾何模型建立

本次模型以工作面走向為X方向，以傾向方向為Y方向，豎直方向為Z方向。根據彈塑性力學理論，開挖后應力變化的影響范為所開挖范圍的3～5倍。此次模擬的開挖尺寸為采長480m，采寬170m。為了滿足計算需要并保證計算精度，取模型走向長度為1500m，傾向長度為550m。模型盡量符合地表實際，最高點473m，最低點436m。詳細的三維模型及其網格劃分見圖2。

表1 巖石物理力學參數

圖2 三維計算模型及網格劃分圖

2.4 邊界條件設定

本模型采用的邊界約束條件：①模型的上邊界地面為自由邊界；②模型的前后、左右邊界均施加水平法向約束；③模型底邊界施加固定約束。

3 數值模擬結果及分析

3.1 數值模擬結果

數值模型工作面采長480m，本次模擬結合觀測情況，對開采過程中五個階段的地表移動變形進行分析。數值模擬首先根據地層物理力學參數、模型邊界條件等計算初始應力場，然后在初始應力場的基礎上按順序模擬開挖，每開挖一次為一個施工階段。進行施工階段分析后，可輸出任意斷面的水平和豎向位移云圖[10]。圖3、圖4為工作面推進230m時，走向主斷面上的水平和豎向位移圖。

觀測點各階段下沉曲線如圖5所示。由圖5可知，隨著開采的進行，位于皮帶走廊的各觀測點產生不同程度的下沉，下沉量自皮帶走廊觀測線中心向兩側逐漸減小。工作面推進230m時，位于回采工作面以外的G35、G38、G61、G63下沉量均在10mm以下，基本不受開采影響。隨著工作面的推進，各觀測點下沉量不斷增加，不均勻沉降問題也愈加明顯。工作面推進至480m時，預計最大沉值位于G50，達到1280mm；最小為63mm。為保證皮帶支架結構的安全，建議解開支架固定螺栓，卸載地表沉陷引發的結構應力。

地表點的傾斜、曲率和水平變形值由節點的下沉差值和水平移動差值計算而得。

圖3 走向主斷面豎向位移云圖

圖4 走向主斷面水平位移云圖

圖5 皮帶觀測點下沉曲線

傾斜曲線見圖6。由圖6可知，隨著工作面的推進，皮帶走廊傾斜變形顯著增加。最大值由0.95mm/m增加至11.7mm/m，大于規程規定的臨界值3mm/m；軌道巷一側最大傾斜變形預計出現在G44～G46之間，皮帶巷一側出現在G54～G55之間。傾斜變形會造成皮帶支架歪斜，進而破壞皮帶穩定性。為減緩傾斜影響，建議在皮帶的支架側面加斜撐。

模擬結果顯示水平變形曲線與曲率曲線形狀相似。曲率曲線見圖7，位于工作面邊界外側的支架將產生拉伸變形和正曲率變形，而工作面中心部位支架受到壓縮變形和負曲率變形的影響。工作面推進到380m時，最大正負曲率變形分別達0.25mm/m2和-0.35mm/m2，超出規程臨界值0.2mm/m2。拉伸變形易引起皮帶支架產生豎向裂縫，而壓縮變形易使皮帶支架產生剪切和水平擠壓裂縫。建議在壓縮變形嚴重部位挖變形縫以減緩壓縮變形對皮帶的影響。

圖6 皮帶觀測點傾斜曲線

圖7 皮帶觀測點曲率曲線

3.2 模擬值與觀測值對比分析

本文采用Midas GTS軟件進行采空區施工過程的動態數值模擬，得出五種工況下管狀皮帶觀測點的移動變形值。表2給出了模擬結果與觀測值的對比。

表2 模擬與實測結果對比

由表2可知，數值模擬預計的移動變形與實際觀測數據存在差距。產生這種結果的原因主要有：①所建模型與真實地表有一定差異；②上覆巖層復雜多變，而所建模型經過一系列簡化，不能完整表達研究區域的地層構造；③地表移動變形模擬的準確性很大程度取決于巖土力學參數的選取。模擬所采用巖層物理力學參數與真實值之間存在差異。

經計算，此次模擬下沉預測誤差在2.4%～7.6%之間，水平移動在1.5%～7.1%之間，傾斜誤差在2.7%～10.4%之間，水平變形誤差在3.3%～10.8%之間，曲率誤差在0%～20%之間。除曲率外，其它移動變形值預測誤差均在11%以內，與實測數據誤差較小。因此，本文認為采用Midas GTS軟件能較好地進行地表移動變形預測。

4 結 論

1)本文采用Midas GTS軟件進行了開采過程的數值模擬，經與皮帶實測移動變形資料對比，除曲率外的移動變形預測誤差均在11%以內，表明Midas GTS作為一種較新的有限元分析軟件，用來預測地表開采沉陷規律是可行的。

2)由于上覆巖層及地形的復雜性，所建三維模型很難完全表達實際情況，為了用數值模擬方法更準確地分析開采沉陷問題，需進一步深入研究。

3)分析結果表明，隨著工作面的推進電廠皮帶將受到采動破壞。為保證電廠的正常運行，建議及時采取保護措施。

[1] 何國清，楊倫．礦山開采沉陷學[M]．徐州：中國礦業大學出版社，1994：35．

[2] 梅松華，盛謙，李文秀.地表及巖體移動研究進展[J].巖石力學與工程學報，2004，23(S1)：4535-4539.

[3] 何萬龍，康建榮．山區地表移動與變形規律的研究[J].煤炭學報，1992，17(4)：1-15.

[4] 秦世界，張和生，李國棟.基于FLAC3D的煤礦開采沉陷預計及與概率積分法的對比分析[J].煤炭工程，2014，46(6)：96-98.

[5] 胡海峰，康建榮.基于有限元法的采動坡體穩定性計算[J].太原理工大學學報，2000，31(4) ：343-345，353.

[6] 郭春穎，李云龍，劉軍柱.UDEC在急傾斜特厚煤層開采沉陷數值模擬中的應用[J].中國礦業，2010，19(4)：71-74.

[7] 劉恩彥，張欽禮，吳鵬，等.后觀音山鐵礦采場結構參數優化研究[J].金屬礦山，2011(9):33-35.

[8] 王新民，柯愈賢，胡威，等.露天轉地下開采地表沉陷預計及安全性分析[J].科技導報，2012，30(25)：27-31.

[9] 白利明.井工開采對興能電廠地表管狀皮帶構筑物的影響研究[J].山西煤炭，2013，33(6)：32-34.

[10] MIDAS Information Technology Co.，Ltd..MIDAS User Guide[R].北京:北京邁達斯技術有限公司，2006.

[11] 金愛兵，明世祥，孫金海.采空區地表沉陷及移動規律的數值模擬[J].采礦技術，2008，8(4)：70-72.

Finite element analysis of effect of underground mining on tubular belt based on Midas GTS

QIN Shao-fen，HU Hai-feng，ZHANG Xue-qin

(School of Mines，Taiyuan University of Technology，Taiyuan 030024，China)

With the mining of coal resources，more and more building are influenced by the subsidence.For the sake of the normal operation of ash discharge pipe of Gujiao xingneng plant，which pass over the 22620 face of Zhenchengdi colliery，the movement and deformation are predicted based on the finite element analysis software Midas GTS.The forecasts and measured data are compared.The results show that，the numerical model established meet with engineering need.When the working face was advanced to 310 meters，the movement and deformation of belt is greater than the critical deformation value of the regulations.The belt will suffer from mining damage.Protection measures for belt are put forward to ensure the normal operation.

underground mining；numerical simulation；movement and deformation；Midas GTS；tubular belt

2015-01-11

山西自然科學基金項目“采煤沉陷區地表破壞及對生態環境的影響研究”資助(編號：2014011001-3)

秦韶芬(1990-)，女，山西晉城人，太原理工大學碩士研究生，主要從事開采沉陷與變形監測方面的研究。E-mail：qsftyut@163.com。

TD5

A

1004-4051(2015)09-0141-04