條帶開采地表運移模擬與實際觀測對比分析

胡文鵬

(山西汾西礦業集團中興煤業有限責任公司，山西 交城 030500)

0 引言

條帶開采作為采出“三下”壓煤的重要途徑之一，由于其成本低廉，技術簡單，并能夠有效地控制地表沉陷和移動，保護地面建筑物，在我國礦區得到廣泛應用。條帶開采煤柱的穩定性及覆巖運動關系到地表建筑物等的安全，因此確定合理的采留比對于控制地表運移，保障地表設施、田地、河流、交通具有重要的意義[1-2]。

由于礦山地質的復雜性和巖體的特殊性，條帶開采引起的覆巖及地表移動和破壞機理比較復雜，上覆巖層的破壞和移動是自下而上逐漸發展的過程，并且受多方面因素的影響[3]。上覆巖層的破壞過程不但是一個連續移動變形的過程，更是一個巖體不斷破碎、斷裂發展的過程。可見，要揭示覆巖運動及地表移動規律，必須深入研究條帶開采過程中覆巖破壞及移動與地表移動的特征，而數值模擬方法可快捷方便地對條帶開采的各影響因素對覆巖破壞移動規律的作用關系進行研究，模擬結果直觀性強，在礦山開采沉陷工程覆巖破壞及地表移動規律研究中具有較好的適用性和廣泛性，但在數值模擬過程往往不能夠完全真實地反映煤層覆巖狀態，難免出現模擬結果與實際有誤差，本文圍繞某礦條帶開采條件(煤柱尺寸為80 m，采寬100 m)對地表的影響進行模擬與地表測量數據對比，為最終確定合理采留比提供參考[4-5]。

1 工程概況

工作面揭露主要煤層為二迭系山西組的3上煤層。根據工作面順槽實際揭露資料和鉆孔柱狀資料，二迭系山西組的3上煤層為黑色，以鏡煤、亮煤為主，平均煤厚5 m，煤層傾角平均5°，可采系數95%，條帶結構，層狀構造，屬于簡單結構。工作面內煤層賦存穩定，接近切眼段煤層傾角較大，本工作面與相鄰工作面較近，頂板裂隙比較發育。

煤層頂底板巖性，見表1。

表1 煤層頂底板情況表

2 條帶開采地表運移規律的數值模擬

2.1 模擬方案

為了研究不同采出率對條帶開采覆巖運動及地表移動規律的影響，為使各個方案的模擬結果有可比性，在進行某個影響因素的分析時，保持其他各影響因素相同。根據模擬的目的制定如下模擬方案。

設煤柱的寬度為80 m(礦方實際煤柱寬度)不變，改變采出條帶的寬度，采出條帶的寬度分別為60 m、80 m、100 m、120 m，共建立4個模型，這4個模型的采深均為500 m(煤層實際埋深)，采厚為5 m，煤層按水平煤層計算，每個方案均開采3個條帶，方案編號為：方案一、方案二、方案三和方案四。

各個方案的邊界條件，巖石力學參數等其他條件保持相同。

2.2 模擬結果分析

模擬結果分析見表2。

根據模擬結果，給出方案一～方案四地表下沉曲線圖及水平移動曲線圖，分別如圖1、圖2所示。

由地表下沉曲線圖1、水平移動值圖2可知，留設煤柱的寬度為80 m，在采寬分別為60 m和80 m時，采出率比較小，地表下沉曲線比較平緩，隨著采寬的增大，地表下沉曲線變的相對較陡。同樣，在采出率較小時，地表水平移動值較小，隨著采出率的增加，水平移動值明顯增大，并且隨著采出率的增加，地表下沉值及地表水平移動值的增幅更大。

圖2 采出率不同時地表水平移動曲線圖

采區寬度的不同代表著采動程度的不同，采區寬度較小時，往往表現為非充分采動，此時地表下沉值未達到最大值；隨著采區寬度的增加，采動程度逐漸接近充分采動，有時會達到超充分采動，這時的地表下沉值達到最大值，不再隨采區寬度的增加而增加。

3 地表變形觀測與對比分析

為了解決數值模擬所采用地質資料的礦井正確的選擇條帶的采留、寬問題以及觀測開采過程中及開采結束后地表移動變形情況，合理優化開采方案，檢驗地表移動數值模擬值的正確性，根據該礦的開采實際及要求，在該采區上方布置地表巖移觀測站，沿條帶采區正上方沿走向布置一條觀測線，分析總結該礦開采引起的地表變形規律，以指導礦井安全生產，依據相關規程規定進行了比較細致的觀測；通過對測區地表觀測數據進行分析，指導礦井生產設計，得出煤礦開采引起的地表變形對地表建(構)筑物的損壞程度等級，以便對井田范圍及周邊的地面設施進行合理保護或搬遷避讓。

3.1 地表移動觀測站布設及觀測方法

本次地表移動觀測線設在條帶開采采區工作面的近正上方，為走向觀測線，觀測線長度為990 m，全線共布設34個測點，點間距為30 m，觀測點布置圖如圖3所示。觀測點編號從東自西起連續編號，此觀測線基本為一條直線，觀測線橫跨3條工作面，根據已有的資料表明相鄰工作面的開采對該觀測站產生的影響很小，因此觀測站位置較為理想。

圖3 觀測點布置圖

3.2 觀測結果分析

觀測數據：每次觀測結束后，檢查觀測數據，使數據符合《煤礦測量規程》的有關規定，然后進行改正計算和平差計算，使觀測結果比較準確。在該采區，所采用的開采設計是留設煤柱80 m，采出條帶寬度為100 m，即模擬中方案三所采用的留寬和采寬，共采出3個條帶，觀測站對各觀測點進行了10次觀測。此采區工作面于2015年12月已經回采結束，到2017年1月18日，地表下沉進入穩定階段，由此可得到最終的地表移動規律。根據觀測數據計算出地表下沉值與水平移動值，其中第10次測量得到的地表下沉值與水平移動值，見表3。

m次觀測時n點的下沉：

Wn=Hn0-Hnm

(1)

式中：Wn—n點的下沉值，mm；Hn0、Hnm—分別為首次和m次觀測時點的高程，mm。

n號點的水平移動：

Un=[(Xnm-Xn1)×cosα+(Ynm-Yn1)×sinα]×1 000

(2)

式中：Un—n號點的水平移動，mm；(Xnm，Ynm)—m次觀測時n點的坐標；(Xn1，Yn1)—第一次全面觀測時n點的坐標；α—水平移動方向的方位角。

表3 觀測點觀測數據

曲線對比：將第10次所測數據導入EXCEL，得到兩者的地表下沉曲線和水平移動曲線，與模擬所得數據作對比，如圖4，圖5所示。

圖4 地表下沉曲線

圖5 地表水平移動曲線

3.3 模擬結果與觀測數據對比

在采寬為100 m，留寬為80 m時，數值模擬地表最大下沉值為323 mm，最大水平移動值為89 mm，而實際觀測最大下沉值為386 mm，最大水平移動值120 mm，數值模擬結果比實際數據要小，誤差主要由建模過程中對巖層參數簡化造成的。數值模擬在對地質條件作了部分簡化，去除了部分較薄的軟弱巖層，使得上覆巖層中的堅硬巖層所占比例增加，有利于控制地表的下沉與移動，形成數值模擬結果比實際數據要小的這種情況，但差距不大。

條帶開采工作面觀測站走向線下沉水平移動和水平變形曲線的發育規律符合開采沉陷的一般規律。根據地表移動規律分析，截止工作面結束工作面開采未達到充分采動。

4 結論

(1)觀測數據和模擬數據相對比，采用FLAC3D預計地表下沉和水平移動值與現場監測數據相差不大，能夠較準確地預計開采引起的地表下沉與變形，證明數值模擬是一種可靠的研究手段。

(2)通過對觀測數據的分析，得出地表下沉和移動的一般規律，地表下沉曲線以最大下沉點對稱，地表變形曲線各點的移動指向最大下沉點。