董良志 等

摘 要：放頂煤開采是一個破壞原巖平衡狀態(tài)的行為，這一行為將引起周圍巖體的移動，煤層開采過程中的覆巖移動規(guī)律是煤礦開采地壓控制的基礎(chǔ)，準(zhǔn)確判定煤層開采后上覆巖層的變形移動規(guī)律對于頂煤松動預(yù)裂技術(shù)優(yōu)化和礦壓防治工作意義重大，本文將巖土工程界運(yùn)用廣泛的數(shù)值模擬軟件FLAC3D用于模擬煤層開采引起上覆巖體的移動，通過計(jì)算機(jī)模擬得到了在采場中隨工作面推進(jìn)的上覆巖體移動規(guī)律。

關(guān)鍵詞：FLAC3D；覆巖；規(guī)律

1 概述

放頂煤開采是一個破壞原巖平衡狀態(tài)的行為，這一行為將引起周圍巖體的移動，放頂煤的開采將導(dǎo)致圍巖移動變形、巖層卸壓使得上部巖層地應(yīng)力降低、裂隙發(fā)育，一方面使煤層透氣性大幅度提高，為被保護(hù)層的卸壓瓦斯抽采提供了可能，另一方面削弱了煤與瓦斯突出的潛能，在對煤層開采工程以及巖石工程等開挖工程問題進(jìn)行研究時，人們普遍采用現(xiàn)場試驗(yàn)、數(shù)值模擬以及相似材料模型相結(jié)合的研究方法，其中采用計(jì)算機(jī)數(shù)值模擬研究方法能較全面的分析煤層開挖過程中上覆巖層的移動變形、應(yīng)力隨煤層開采的時空變化規(guī)律。放頂煤開采過程中會導(dǎo)致煤巖層頂?shù)装搴蜕细矌r層的移動變形，在頂?shù)装逡欢ǚ秶鷥?nèi)造成地應(yīng)力降低，研究煤層開采上覆巖層的移動變形規(guī)律有助于進(jìn)一步研究煤層壓力、滲透性變化及頂煤松動預(yù)裂技術(shù)優(yōu)化提供基礎(chǔ)。

2 FLAC3D簡介

FLAC起源于流體動力學(xué)，是快速拉格朗日差分分析的簡稱，該分析方法最早由用于固體力學(xué)領(lǐng)域，F(xiàn)LAC3D程序是快速拉格朗日差分分析方法的三維應(yīng)用，是FLAC模擬軟件向三維空間的擴(kuò)展。目前已經(jīng)成為采礦工程等巖石力學(xué)計(jì)算中的最重要數(shù)值方法之一，該程序軟件在邊坡穩(wěn)定性評價、動力學(xué)、地下洞室施工設(shè)計(jì)、地質(zhì)構(gòu)造及塊體運(yùn)動學(xué)、巷道支護(hù)設(shè)計(jì)及評價、隧道工程、河谷演化進(jìn)程再現(xiàn)、礦山工程等多個領(lǐng)域研究中得到廣泛應(yīng)用，特別適合求解巖土力學(xué)工程中非線性的大變形問題。

3 煤層賦存情況

孔莊煤礦井田內(nèi)主要含煤地層有太原組、山西組及下石盒子組，含煤地層總厚度約487.37m，含煤20多層，煤層平均總厚度約15.79m，含煤系為3.24%，其中可采或局部可采煤層5層（7、7下、8、17、21）總厚度約9.79m，主要煤層含煤系數(shù)為2.01%。井田內(nèi)的主要可采煤層為7煤層和8煤層，本礦井尚未開采17、21兩層煤。

7號煤層位于山西組中下部，上距下石盒子組底部分界砂巖約72m，下距8煤4.17～40.18m之間，西翼巖漿巖侵入嚴(yán)重，煤層兩極厚度為2.00～6.52m，一般厚度為4.54m，7煤結(jié)構(gòu)較簡單，只有少數(shù)工作面煤層中有夾矸，夾矸厚0.40～1.40m，另外在7251工作面出現(xiàn)了面積為60000m2含矸區(qū)域，夾矸最大厚度1.20m，夾矸將煤層分為上3.70m，下為0.80m，夾矸巖性多為砂質(zhì)泥巖、細(xì)砂巖。如圖1所示：

目前7#煤層采用走向長壁后退式綜采放頂煤采煤方法，其中機(jī)采2.5m，放頂煤層厚度2.3m，由于頂板處煤層層理、裂隙較發(fā)育，易破碎垮落，局部揭露煤層底板泥巖，層理、節(jié)理較發(fā)育，泥質(zhì)膠結(jié)，穩(wěn)定性差，為防止冒頂事故發(fā)生，需要對放頂過程中的上覆巖層移動變形規(guī)律進(jìn)行研究，以便為開采設(shè)計(jì)提供科學(xué)依據(jù)?，F(xiàn)采用數(shù)值模擬研究的方法研究走向長壁后退式綜采放頂煤采煤時對上覆巖體的影響。

4 模擬模型的建立

4.1模型模擬參數(shù)選取 本數(shù)值模擬試驗(yàn)建立工作面沿走向模型，通過數(shù)值模擬計(jì)算，從理論上研究放頂煤層開采對上覆煤巖體移動變形影響作用。

FLAC3D數(shù)值模擬軟件中包含了11種材料本構(gòu)模型，其中包括七種塑性模型和三種彈性材料模型以及一種空單元模型，巖石工程界常用的摩爾-庫倫模型就屬于七種塑性模型中的一種，針對開采保護(hù)層的物理模型，我們采用摩爾-庫倫模型。

考慮到使用摩爾-庫倫模型時剪切模量、剪脹角、密度、粘聚力、體積模量、摩擦角等

是摩爾-庫倫模型所需輸入的參數(shù)，該模型的破壞包絡(luò)線對應(yīng)于摩爾-庫倫判據(jù)加上拉伸分離點(diǎn)。由于本次模擬研究研究具有一般性，目的在于得到各個因素對應(yīng)力變化的影響情況，所以本次數(shù)值模擬將所需的參數(shù)取實(shí)測的平均值，各個參數(shù)取值如表1：

4.2 模擬區(qū)域范圍確定 一般來說煤層開采后頂板覆巖移動遠(yuǎn)大于底板巖層移動，本文主要研究煤層放頂開采后上覆煤巖移動規(guī)律?？紤]到煤層、頂板及底板三者之間的變形協(xié)調(diào)規(guī)律以及煤層采出后在包括頂?shù)装逶趦?nèi)的圍巖原始應(yīng)力要發(fā)生變化、產(chǎn)生變形、移動乃至破壞，其波及范圍往往很大，但是按照圣維南原理，其影響范圍仍有一定的限度。為最大限度模擬本區(qū)域煤層開采后的上覆巖層移動變形規(guī)律，模擬范圍的確定是一個難點(diǎn)，模擬區(qū)域過小往往不能準(zhǔn)確的再現(xiàn)開挖工程中周圍巖體的移動、變形、破壞規(guī)律，模擬區(qū)域過大則計(jì)算時間長、工作量大。本次數(shù)值模擬模型的范圍：煤層頂板以上3.9m砂質(zhì)泥巖外加30m的間接頂板，煤層底板以下2.5砂質(zhì)泥巖外加30m的間接底板以及4.65m厚的煤層，考察停采線和開切眼各向外擴(kuò)展50m作為模型的左、右邊界，傾向模型的上、下、左、右邊界取法同走向模型。傾向模型和走向模型都采用類似的約束邊界條件，煤層邊界以上的巖層作為外載荷施加在模型的上邊界上，其他方向采用應(yīng)力邊界條件。

4.3 模擬網(wǎng)格的劃分 煤層在放頂開采后所導(dǎo)致的上覆煤巖移動、破壞是一個漸進(jìn)的演變過程，隨著采空區(qū)范圍的不斷擴(kuò)大，應(yīng)力的不斷調(diào)整而發(fā)展變化的，必須通過對開采過程進(jìn)行數(shù)值模擬才能夠使之轉(zhuǎn)化成實(shí)際的動態(tài)問題來掌握其發(fā)生、發(fā)展演化史。據(jù)此結(jié)合工作面的推進(jìn)進(jìn)度，將走向模型分解成一連串的4個模型分別模擬開采10m，30m，50m，100m進(jìn)行分步開采。

接著要將網(wǎng)格離散，網(wǎng)格的離散原則上是劃分得越細(xì)越好，但是由于計(jì)算機(jī)容量和計(jì)算速度的限制，不可能將網(wǎng)格劃分得過細(xì)。本次模擬試驗(yàn)研究將走向、傾向模型均劃分成132800個單元，同時考慮到邊界效應(yīng)，網(wǎng)格劃分采用等間距與不等間距相結(jié)合。

5 上覆巖體位移變化規(guī)律模擬

通過對上述模型進(jìn)行煤層的開挖模型模擬，分別模擬了煤層開采的工作面推進(jìn)10m、30m、50m、100m時上覆巖體的豎向的移動變形規(guī)律，模擬結(jié)果如圖2～5所示。

6 模擬結(jié)果分析

圖6即為煤層工作面上覆巖體觀測線垂直移動量隨開采長度的變化規(guī)律圖，當(dāng)煤層開采的工作面推進(jìn)到30m時，采空區(qū)上部觀測線移動量左右對稱，最大向下移動量為11.8mm，工作面煤柱邊沿下部觀測線未出現(xiàn)下移，基本呈倒V字形。

隨著工作面的進(jìn)一步推進(jìn)，采空區(qū)越來越大，設(shè)置在上覆巖層中的觀測線移動量也相應(yīng)增加，移動量向工作面?zhèn)绕?，由左右對稱逐漸向非對稱轉(zhuǎn)移，而且隨著工作面的推進(jìn)，下移量由小變大，在煤柱邊沿下部也逐漸出現(xiàn)了向下移動變形的狀況，當(dāng)煤層工作面開采到100m時，觀測線最大下移量為89.0mm。

7 結(jié)論及下一步工作

通過數(shù)值模擬得到采場中隨工作面推進(jìn)的上覆巖體移動規(guī)律，工作面推進(jìn)后，采空區(qū)上部巖層隨之向下移動，且移動呈左右對稱，基本呈倒V字形，煤柱底板側(cè)幾乎不移動；隨著推進(jìn)距離加大，上部巖層移動量加大，且底板側(cè)移動量大于頂板側(cè)，移動中心整體向底板側(cè)偏移。本次模擬的數(shù)據(jù)量不夠大，下一步一方面要以10m為一個單元詳細(xì)對比分析上部巖層隨推進(jìn)距離影響變化量，計(jì)算得到頂板最活躍范圍；另一方面，將模擬推進(jìn)距離增大到500m，計(jì)算得到具體推進(jìn)多少米后覆巖不再移動，重新穩(wěn)定達(dá)到平衡，從而指導(dǎo)生產(chǎn)。

參考文獻(xiàn)：

[1]閆長斌，徐國元，李夕兵.爆破震動對采空區(qū)穩(wěn)定性影響的FLAC 3D分析[J].巖石力學(xué)與工程學(xué)報，2005，24（16）：2894-2899.

[2]彭文斌.FLAC 3D實(shí)用教程[M].北京，機(jī)械工業(yè)出版社，2008：110.

[3]何學(xué)秋.含瓦斯煤巖流變動力學(xué)[M].徐州：中國礦業(yè)大學(xué)出版社，1995.

[4]蔣承林.煤與瓦斯突出的球殼失穩(wěn)機(jī)理及防治技術(shù)[M].徐州：中國礦業(yè)大學(xué)出版社，1998.

[5]劉寶安.下保護(hù)層開采上覆煤巖變形與卸壓瓦斯抽采研究[D].合肥：安徽理工大學(xué)，2006.

作者簡介：董良志（1982-），江蘇沛縣人，本科學(xué)歷，2013年1月畢業(yè)于中國礦業(yè)大學(xué)采礦工程專業(yè)，助理工程師，歷任上海大屯能源股份有限公司孔莊煤礦技術(shù)員等職務(wù)，從事采掘技術(shù)管理工作。