韓東赟， 譚榮建

(昆明理工大學國土資源工程學院，云南昆明 650093)

?

地下采煤地表沉陷規(guī)律數(shù)值模擬研究

韓東赟， 譚榮建

(昆明理工大學國土資源工程學院，云南昆明 650093)

摘要為了研究煤層開采對地表的影響及變形規(guī)律，以柳樹溝煤礦一采區(qū)的113001工作面為研究背景，利用FLAC3D數(shù)值模擬軟件建立模型，選取合適的參數(shù)，對該工作面的開采沉陷進行模擬，生成地表下沉曲線、水平移動曲線，通過對下沉曲線和水平移動曲線的對比分析，得出柳樹溝煤礦一采區(qū)113001工作面開采地表移動變形規(guī)律。總體來說，由于煤層埋藏較深，開采厚度不大，其直接頂板巖層冒落巖塊體積膨脹，所以不會造成該工作面采煤完成后出現(xiàn)大規(guī)模的地表沉降。

關鍵詞開采沉陷；FLAC3D；地表變形；數(shù)值模擬

Numerical Simulation Study on the Laws of Surface Subsidence in Underground Mining

HAN Dong-yun, TAN Rong-jian (Faculty of Land Resource Engineering, Kunming University of Science and Technology, Kunming, Yunnan 650093)

Abstract To study the effects and deformation regularity of coal seam mining on surface, taking 113001 working face of the first district in Liushugou mine as research background, using numerical simulation software FLAC3Dto establish model, selecting appropriate parameters, the mining subsidence of the working surface was simulated, surface sinking curve and horizontal displacement curve were generated. Through the contrast analysis between the curves, basic laws of the surface deformation caused by coal mining were obtained. In general, due to deep buried coal sean and small mining thickness, the volume of roof caving rock mass increased large scale surface subsidence caused by coal mining could not occur.

Key wordsMining subsidence; FLAC3D; Ground deformation; Numerical simulation

地下礦層開采后，采空區(qū)周圍巖層的應力平衡狀態(tài)遭到破壞，應力重新分布，達到新的平衡，在此過程中，使巖層和地表產生連續(xù)的移動、變形和非連續(xù)的破壞(開裂、冒落等)，這種現(xiàn)象稱為“開采沉陷”。若地下開采的范圍較小、開采的礦物的埋藏深度較大，則開采沉陷波及的范圍往往只局限于開采區(qū)域周圍的巖體；若開采范圍較大、開采礦物的埋藏深度較小，則開采沉陷波及的范圍就會從巖體發(fā)展到地表，引起地表移動[1-2]。煤礦的大規(guī)模開采，既給人類帶來巨大的經濟效益和社會效益，也對人類生存環(huán)境帶來了一系列消極的影響，開采沉陷會導致地表耕地受到損毀；平原地區(qū)大片地表移動盆地積水；地表水系遭到破壞，影響居民生活用水和農田灌溉用水；使地表建筑物、鐵路、公路等生產生活設施遭到破壞，影響人們的生命財產安全[1-2]。因此，研究開采沉陷具有非常重要的意義。

對于開采沉陷的理論預測，國內外專家提出了多種方案，我國經過40多年的研究和實踐，開采沉陷研究逐漸成熟。目前，比較常用的地表移動變形計算方法有概率積分法、負指數(shù)函數(shù)法、威布爾函數(shù)法和典型曲線法等，其中概率積分法具有參數(shù)容易確定、實用性強等優(yōu)點，在各礦區(qū)使用比較廣泛，是最常用的方法，但概率積分法不能有效表現(xiàn)地層內部巖體的移動變形情況[3-4]。筆者采用FLAC3D數(shù)值模軟件對宣威市文興鄉(xiāng)柳樹溝煤礦一采區(qū)113001采煤工作面進行開采沉陷模擬，為治理該地區(qū)地表沉陷提供理論依據(jù)。

1研究區(qū)概況

柳樹溝煤礦位于云南省宣威市35°方位，平距47 km，地處宣威市文興鄉(xiāng)境內。礦井主要可采煤層為C30煤層，煤層位于宣威組第一段第一亞段(P2x1-1)中下部，下距峨眉山組玄武巖約8 m。煤層傾向北西，傾角8°～9°，屬緩傾斜煤層。煤層厚1.75～2.38 m，平均厚2.06 m，結構單一，一般不含夾矸。煤層頂板巖性主要為細砂巖，局部為粉砂巖、泥質粉砂巖，底板巖性主要為粉砂巖，局部為泥質粉砂巖、粉砂質泥巖。C30煤層在礦區(qū)層位及厚度均穩(wěn)定，屬穩(wěn)定的中厚煤層。礦山開采方式為地下開采，采用斜井開拓；采煤方法采用長壁式采煤法，后退式回采，一次采全厚；采用全部垮落法管理頂板，礦井投產時一采區(qū)的113001工作面長度為80 m，回采工作面推進495 m。

2模型建立及數(shù)值計算

2.1FLAC3D簡介FLAC3D是美國ITASCA咨詢集團公司開發(fā)的三維快速快速拉格朗日分析程序，該程序能較好地模擬地質材料在達到強度極限或屈服極限時發(fā)生的破壞或塑性流動的力學特性，特別適用于分析漸進破壞失穩(wěn)以及模擬大變形。目前，F(xiàn)LAC3D已廣泛應用于巖土工程、采礦工程、隧道工程、道路與鐵路工程等領域的科學研究中[5-7]。

2.2主要地層礦區(qū)出露的地層從老至新分為二疊系上統(tǒng)峨眉山玄武巖組(P2β)、二疊系上統(tǒng)宣威組(P2x)、第四系(Q)。峨嵋山玄武巖組(P2β)巖性主要為灰、深灰色、灰綠色、黑綠色致密塊狀玄武巖，厚度大于50 m。宣威組第一段第一亞段(P2x1-1)底部為淺灰色、灰色薄層狀泥質粉砂巖、含凝灰質泥巖；中下部為C30煤層，煤層厚1.75～2.38 m，平均厚2.06 m；中上部為淺灰色中厚至厚層狀凝灰?guī)r屑，玄武巖屑砂礫巖、砂巖夾泥質粉砂巖、泥巖和煤線。宣威組第一段第二亞段(P2x1-2)底部為灰綠色中厚層狀細粒玄武巖屑砂巖；中上部為灰色薄層狀菱鐵質粉砂巖、細粒玄武巖屑砂巖與泥質粉砂巖、粉砂質泥巖、泥巖、薄煤呈不等厚互層。第四系(Q)巖性為殘坡積物、沖、洪積物及耕植土，結構松散。

2.3模型建立采用采礦工程運用FLAC3D解決問題時常用的單元類型矩形建模。為了消除開挖時模型邊界效應的影響，模型建立盡量大些。模型選取沿煤層走向為Y軸方向，長度取800 m；煤層傾向為X軸方向，長度取300 m；鉛垂方向為Z軸方向，高度為187 m。模型共計33 600個網格和37 023個節(jié)點，模型的煤層厚2.06 m，工作面沿煤層傾向布置，即在模型中沿X軸方向布置[8-11]。建立的初始模型見圖1。

圖1　 FLAC3D建立的模型Fig. 1　Model constructed with FLAC3D software

2.4本構模型及主要力學參數(shù)本構模型是對巖土材料力學性質的經驗性描述，表達的是外載條件下巖體、土體的應力-應變關系，因此，本構模型的選擇是數(shù)值模擬的一個關鍵性步驟。摩爾-庫倫模型代表的材料類型為松散或膠結的粒狀材料，如土體、巖石、混凝土，為巖土力學通用模型，在邊坡穩(wěn)定和地下開挖中應用廣泛。故此次數(shù)值計算的本構模型采用摩爾-庫倫準則。選擇計算模型本構關系后，需確定模擬計算的各巖層的主要力學參數(shù)，包括體積模量、剪切模量和內聚力、內摩擦角、抗拉強度。此次模擬計算所用的基本力學參數(shù)根據(jù)礦井地質勘探資料提供的地質柱狀圖、生產地質報告巖石力學化驗成果及相關文獻給出的力學參數(shù)綜合確定，具體見表1。

2.5確定邊界條件及初始應力狀態(tài)根據(jù)現(xiàn)場模型的分析和研究，研究范圍內的煤巖體系處于一個無窮大地殼平面內，在其邊界處由于鄰近地層的作用，其變形可認為是0。因此對模型邊界變形的設定情況為X軸正負兩個方向上約束為0，Y軸正負兩個方向上約束為0，Z軸負方向上約束為0，Z在正方向上可自由變形。由于實際工程模型處于重力場中，所以在模型進行運算時，施加重力場作用，重力加速度為9.8 m/s2。

表1　各巖層力學參數(shù)

2.6求解及開挖模擬該案例通過監(jiān)控點(150，400，187)在Z軸方向的位移迭代變化和模型中最大不平衡力來判斷模型是否達到平衡狀態(tài)。達到平衡狀態(tài)后進行開挖，工作面長80 m (即X軸方向上110～190 m之間)，沿煤層走向開挖，從模型Y軸方向150 m處開始，到模型Y方向645 m結束，每步開挖100 m，一共開挖5次。

3模擬結果分析

通過三維有限差分數(shù)值模擬軟件FLAC3D，計算位移場分布特征，獲得其位移等值線云圖，開挖后的位移變形見圖2～4。煤層采出后，底板在水平方向受壓，在垂直方向上減壓，造成底板向采空區(qū)方向隆起；其直接頂板巖層由于失去支撐而破碎、冒落、堆積在采空區(qū)，由于冒落巖塊體積膨脹，冒落帶達到一定程度后自行終止，經過一定的沉降變形，最終達到穩(wěn)定狀態(tài)。隨著回采工作面的不斷推進，地表下沉盆地在地表的影響范圍不斷擴大，下沉盆地呈碗狀。

圖2　回采工作面推進495 m沉降云圖Fig.2　Advancing 495 m settlement cloud image of working face

圖3　回采工作面推進495 m沿X軸方向位移云圖Fig.3　Advancing 495m along X axis direction displacement cloud of mining face

圖4　回采工作面推進495 m沿Y軸方向位移云圖Fig.4　Advancing 495m along Y axis direction displacement cloud of mining face

為了探究煤層開采時地表移動變形規(guī)律，模型計算結束后采用一種FLAC3D常用的后處理方法，利用“print”輸出命令配合log文件記錄的方式對計算結果進行進一步處理，提取某些關鍵節(jié)點的變形值，運用Excel繪制沉降曲線和水平移動曲線。在未設置歷史監(jiān)測變量的情況下，可通過手動查找的方式找出關鍵節(jié)點對應的ID。該研究在模型地表中心，沿傾向和走向布置兩條觀測線，傾向觀測線每隔15 m布置1個測點，共計21個；走向觀測線每隔20 m布置1個測點，共計41個[12-13]。

3.1傾向觀測線上地表移動傾向觀測線上的地表移動情況見圖5和圖6。由圖5并結合圖2可知，在傾向觀測線上，地表下沉量最大值出現(xiàn)在觀測線的中心處，遠離觀測線中心，地表下沉量逐漸減小。當回采工作面推進100 m時，煤層開采對地表基本沒有影響，回采工作面推進到200 m時，地表最大下沉量增加到5 mm，隨著回采工作面的推進，地表下沉量顯著增加，下沉速度逐漸加快。煤層開采完后，地表最大下沉量達24.8 mm。

由圖6及圖3可知，地表水平移動量分為兩個區(qū)段，正值表示測點向礦層的上山方向移動，負值表示測點向礦層的下山方向移動，這是由于以下山方向的觀測線控制點作為計算的起始點這一規(guī)定引起的。水平移動發(fā)生在盆地的邊緣區(qū)，指向盆地的中心，由于傾角的影響，下山方向一側的水平移動量與上山方向一側的水平移動量不對稱，下山方向一側的水平移動量大于上山方向一側的水平移動量。在正水平移動區(qū)，在煤柱邊緣處水平變形量達到最大值，向兩側逐漸減小，向右側在地表最大下沉處水平移動量為0；在負水平移動區(qū)，煤柱邊緣處水平移動量達最大值，兩側逐漸減小，向左側在地表最大下沉處水平移動量為0。隨著工作面的推進，地表水平移動量變化的特點保持不變，只是最大水平移動量在增加。

圖5　傾向觀測線上地表下沉曲線Fig.5　Surface subsidence curve of inclination observation line

圖6　傾向觀測線上地表水平移動曲線Fig.6　Surface horizontal movement curve of inclination observation line

3.2走向觀測線上地表移動走向觀測線上的地表移動情況見圖7和圖8。由圖7可知，在走向觀測線上，地表最大下沉量出現(xiàn)在采空區(qū)上方中心處，地表下沉盆地的范圍和最大下沉量隨著回采工作面的推進而增大。當回采工作面推進100 m時，對地表的下沉影響較??；當回采工作面推進到495 m時，地表最大下沉量達24.8 mm，在最大下沉值處的兩側下沉量逐漸減小。隨著煤層的開挖，地表下沉速度逐漸增大，當回采工作面推進到一定距離時，下沉速度逐漸減小。

圖7　走向觀測線上地表下沉曲線Fig.7 Surface subsidence curve of the observation line

由圖8可知，走向觀測線上的地表水平移動量也分2個區(qū)段，在開切眼側的正水平移動區(qū)和在停采線側的負水平移動區(qū)。采空區(qū)上方的地表水平移動量零點處對應地表最大下沉值處，隨著回采工作面推進，地表水平移動0點向停采線方向移動?；夭晒ぷ髅嫱七M100 m時，對地表水平移動影響較小，隨著回采工作面的繼續(xù)推進，地表水平移動量逐漸變化，最大水平移動量增加。

圖8　走向觀測線上地表水平移動曲線Fig.8　Surface horizontal movement curve of the observation line

4結語

該研究采用FLAC3D數(shù)值模擬軟件模擬柳樹溝煤礦一采區(qū)113001采煤工作面開采對地表變形的影響，結果表明：隨著回采工作面的不斷推進，地表下沉盆地在地表的影響范圍不斷擴大，下沉盆地呈碗狀。當C30煤層回采完后，地表最大沉降值為24.8 mm，最大值發(fā)生在采空區(qū)上方地表略偏向傾向方向，但由于傾角比較小，這種趨勢不是特別明顯；從傾向方向地表水平位移看，下山方向位移最大值為-9.7 mm，上山方向位移最大值為7.8 mm，可以看出在傾向方向水平位移不對稱，下山方向大于上山方向。走向方向地表水平位移正最大值為8.4 mm，負移最大值為-8.3 mm，正負位移基本上呈對稱分布?？傮w來說，由于煤層埋藏較深，開采厚度不大，其直接頂板巖層冒落巖塊體積膨脹，所以柳樹溝煤礦一采區(qū)113001工作面采煤完成后不會造成大規(guī)模的地表沉降。

參考文獻

[1] 何國清，楊倫，凌庚娣，等．礦山開采沉陷學[M]．徐州：中國礦業(yè)大學出版社，1991．

[2] 侯長祥，馮濤，熊仁欽，等.礦床“三下一上”開采[M].北京：煤炭工業(yè)出版社，2001．

[3] 國土資源部土地整理中心． 土地復墾方案編制實務[M]．北京：中國大地出版社，2011．

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

[5] 彭文斌.FLAC3D實用教程[M].北京：機械工業(yè)出版社，2007.

[6] 陳育民，徐鼎平.FLAC/FLAC3D基礎與工程實例[M].北京：中國水利水電出版社，2013.

[7] 孫書偉，林杭，任連偉.FLAC3D在巖土工程中的應用[M].北京：中國水利水電出版社，20011.

[8] 武崇福，劉東彥，方志.FLAC3D在采空區(qū)穩(wěn)定性分析中的應用[J].河南理工大學學報，2007,26(02):136-140．

[9] 朱慶偉，蔣軍.煤層開采沉陷對橋梁建筑物的影響分析[J].測繪科學，20014,39(4):80-91．

[10] 朱強，高明中，孫超，等.急傾斜煤層開采地表移動規(guī)律數(shù)值模擬研究[J].安徽理工大學學報，2012,32(3):71-74．

[11] 盧志剛，劉興權，唐義宏.基于GIS和FLAC3D礦山地表沉陷可視化[J].現(xiàn)代礦業(yè)，2011(11):44-46．

[12] 侯利斌，翟英達，韓偉.李村煤礦開采沉陷規(guī)律數(shù)值模擬研究[J].煤礦安全，2013,44(4):58-61．

[13] 劉劍，郝兵元，黃輝．大采寬工作面開采沉陷規(guī)律的FLAC3D模擬研究[J].煤炭技術，2014,33(11)：132-134．

收稿日期2015-12-28

作者簡介韓東赟(1990- )，女，云南紅河人，碩士研究生，研究方向：土地整治、礦區(qū)土地復墾、地表沉陷。

中圖分類號S 289；TD 82

文獻標識碼A

文章編號0517-6611(2016)03-243-04