蘆嶺煤礦8煤層覆巖移動規(guī)律數(shù)值模擬

鄒立強，陳仕剛，丁宏偉

（山東省煤田地質局第二勘探隊，山東濟寧272000）

蘆嶺煤礦8煤層覆巖移動規(guī)律數(shù)值模擬

鄒立強*，陳仕剛，丁宏偉

（山東省煤田地質局第二勘探隊，山東濟寧272000）

煤層開采將引起原有應力平衡的破壞，導致覆巖的變形與破斷，不利于礦井的安全生產(chǎn)?；趲r石損傷力學的基本原理，分析了覆巖變形失穩(wěn)的條件。結合礦區(qū)鉆孔99-4揭露的開采地質條件和巖石力學性能等參數(shù)，運用巖石破裂過程分析系統(tǒng)（RFPA2D）建立了力學模型。結果表明：覆巖中的離層面和縱向破斷面的張開和閉合是一個動態(tài)演化過程；細中砂組合覆巖采動中直接頂基本是在工作面推進40m時發(fā)生離層破壞，而基本頂則是在工作面推進約65m時整體垮落，之后周期垮落步距10～15m。

煤層開采；覆巖失穩(wěn)；移動規(guī)律；動態(tài)演化；數(shù)值模擬

煤層開采將引起上覆巖層應力重新分布，使之產(chǎn)生變形、破斷和滲透變化等一系列復雜的工程地質作用過程。采用數(shù)值模擬的研究方法對上覆巖層的裂隙發(fā)育情況與其破壞規(guī)律進行的分析，目前國內外的研究角度較為全面，研究成果也較為豐碩［1-5］。其中對上覆巖層的結構研究中，各種梁理論的建立在時間上較早［6-8］。陳忠輝、馬占國等人借助損傷力學的理論建立一個簡單的力學模型或從時變邊界的力學分析總結覆巖的移動破壞機理［9-10］?；诟矌r采動失穩(wěn)條件的分析，并且利用鉆孔對巖芯進行巖石物理、力學性質參數(shù)測定，在分析蘆嶺煤礦8煤層頂部巖石結構及力學性質特征研究的基礎上，建立了以此為依據(jù)的力學模型，預期演繹出8煤層在真實開采的條件下上覆巖層移動與破斷規(guī)律，依此利用既得規(guī)律特征有效規(guī)避采場頂板事故，為采場進行科學設計提供依據(jù)。

蘆嶺煤礦位于宿縣礦區(qū)東部，宿東向斜西南翼的東南段。8煤層屬于下石盒子組沉積，與上部煤層底板間距平均20m，下部煤層頂板平均3.5m，沉積巖則以砂巖、粉砂巖和泥巖為主。本次模型數(shù)據(jù)源為三水平一采區(qū)上的99-4鉆孔，主要模擬位于井田中部F7～11-12線之間的區(qū)域，該區(qū)發(fā)育主要有斜切斷層F40、FD17，根據(jù)斷層分形理論計算［11］，其值為0.79248，相關系數(shù)為0.9776，因此相對構造比較簡單。研究區(qū)頂板粒度自下而上變粗，基本頂細中砂為主，基本構成則是長石以及石英砂巖，其中大型楔形交錯層理時有出現(xiàn)。

1　采動后覆巖失穩(wěn)垮落的條件分析

采場覆巖破壞有很多類型，例如離層、冒落、層間錯動，剪切破壞和塑性變形，塊體移動、沿軟弱結構面的滑動，巖爆、煤爆和片幫等，可以分為4種破壞機制：張破壞、剪破壞、結構體的滑落和轉動［12-13］。

張破壞是煤巖張應變達到或者超過巖體的損傷閾值造成的，其理論公式為：

式中：εi（i=1，2或3）——巖體在重分布應力作用下i個方向上的張應變；

εi，0——巖體的極限張應變。

剪破壞的判據(jù)可用Mohr-Coulomb強度理論：

式中：τ——抗剪強度；

σ——法向應力；

C、φ——巖體的粘聚力和內摩擦角。

巖體應變的同時巖體內力會重新分布，然而巖體結構并不會隨即變化，而是產(chǎn)生“梁”這種結構模型，其穩(wěn)定狀態(tài)的持續(xù)關鍵可能決定于巖層移動中形成離層區(qū)附近的幾個巖塊，在文獻［13］中的“S-R”穩(wěn)定條件分析即是對砌體梁的滑落與轉動變形失穩(wěn)條件的詳細論述。

覆巖失穩(wěn)垮落條件的建立與選用，比較系統(tǒng)的反映了覆巖在采動中的應力場的變化，為覆巖的失穩(wěn)變形類型提供了依據(jù)，有利于對覆巖失穩(wěn)破壞機理的分析與研究。

2　數(shù)值模擬模型

2.1數(shù)值模擬軟件的選取

采場上覆巖層的運動是一個非常復雜的動態(tài)發(fā)展過程。由于覆巖材料存在著非勻質性，破壞的非連續(xù)性，力學行為的非線性等特點［14-16］，從而以往使用的模擬系統(tǒng)都存在著或多或少的局限性。東北大學巖石破裂與失穩(wěn)研究中心開發(fā)的巖石破裂過程分析系統(tǒng)RFPA2D是利用有限元的理論基礎，充分考慮了巖體的各向力學性能的不一致性，以各微小基元的參數(shù)通過宏觀方式表現(xiàn)出來，對巖體的本構關系模擬較為具體。RFPA2D分析系統(tǒng)假定細觀基元的力學性質服從某種分布，以物理統(tǒng)計方法描述其離散性，建立了細觀與宏觀介質力學性能的聯(lián)系。因此采用巖石破裂過程分析軟件RFPA2D系統(tǒng)會更真實的模擬覆巖采動的動態(tài)發(fā)展過程，更好的分析其變形破斷規(guī)律、機理。

2.2覆巖巖石力學參數(shù)的測定

本次模型采用鉆孔99-4的巖芯材料作為基本數(shù)據(jù)來源。表1是嚴格按照巖石力學試驗規(guī)程，對8煤層上覆關鍵巖層進行力學性能進行測試的結果。根據(jù)測試結果泥巖抗壓強度在7.6～32.5MPa之間，平均在25MPa，彈性模量在10.2～25.1GPa之間，平均在17.65GPa。粉砂巖抗壓強度在51.9～107.1MPa之間，平均在79.5MPa，彈性模量在33.1～50.6GPa之間，平均在41.85GPa。粗砂巖抗壓強度為131.4MPa，彈性模量為27.633GPa。

表1　數(shù)值模擬模型中各巖層的主要力學參數(shù)

2.3模擬模型基本數(shù)據(jù)

模型建立以無限接近相應區(qū)域的原始尺寸為原則，在水平方向上選用了280m，而在上下方向的選用上，采取了采動過程中短期可能夠影響到的130m厚度，這樣按1m×1m的基元劃分將模型歸為36400個基元的系統(tǒng)。另外，短期影響不到的上覆巖層根據(jù)地應力化為與之相應的載荷。為防止邊界效應的影響，在模型的兩端分別預置保護煤柱50m。近似模擬采煤過程，以2m采高、5m為步距分步分階段向前推進。在邊界的兩端采用水平約束，而在邊界的上部巖層可以存在水平滑動。巖層組合模型見如圖1所示。

3　模型模擬煤層開采過程中上覆巖層的移動特點

在切眼開拓的瞬間，伴隨著應力的重新分布（應力大小在模型中由亮度表示，亮度越大表示應力越大），當開挖近10m時，在采空區(qū)兩端應力出現(xiàn)明顯集中，而在采空區(qū)上部有限范圍內亮度變暗，說明地應力有所減弱，底板影響不明顯。當采空區(qū)跨度達30m時，直接頂內出現(xiàn)黑點，說明直接頂移動破壞開始。這是由于層間滑動引起直接頂?shù)膹埿云茐?，發(fā)生了直接頂?shù)膶娱g錯動—離層—冒落（圖2a-b）。

采空區(qū)跨度為40m時，直接頂陸續(xù)產(chǎn)生多處黑點，隨著工作面的推進黑點進一步擴大，而且與上位巖層之間產(chǎn)生線狀黑色區(qū)域。說明隨著直接頂內部張力的不斷發(fā)展促成了離層裂隙的出現(xiàn)，同時基本頂內部也出現(xiàn)了應力集中，逐漸發(fā)育有豎直方向的裂隙與離層裂隙，裂隙帶繼而出現(xiàn)。采空區(qū)跨度不斷擴大，直接頂出現(xiàn)大面積的暴露，煤壁位置剪切應力持續(xù)加大，其上部出現(xiàn)明顯裂隙。開挖至55m時，直接頂大面積的垮落，顯然此時基本頂以“砌體梁”形式支撐著上覆巖層。另外，基本頂也逐漸移動下沉，裂隙進一步向上發(fā)展（圖2c-d）。

圖1　巖層組合數(shù)值模擬模型

圖2　數(shù)值模擬結果

隨著工作面的向前推進，基本頂內部裂隙逐漸貫通。采空區(qū)跨度達到65m時，基本頂向下移動速度加快，承載能力加速下降并最終破斷，形成懸臂梁結構。此時采空區(qū)后方煤壁上部剪應力集中明顯，隨后整體垮落。采空區(qū)上部巖層協(xié)同垮落，裂隙急速向上擴展至75m左右（圖2e）?；卷敵醮慰迓浜?，工作面繼續(xù)推進時上部巖層移動速度和規(guī)模則相對穩(wěn)定，然而工作面前方時變煤壁上方應力卻逐漸變大，向前推進至75m時，基本頂上方又一次出現(xiàn)破斷，但離層裂隙向上發(fā)育極為有限。之后基本頂變化又處于短暫穩(wěn)定狀態(tài)，并逐漸形成了穩(wěn)定的冒落帶（圖2f）?；卷斂迓浜?，冒落帶由于壓實作用在采動過程中形成的裂隙又逐漸閉合。

4　結論

（1）煤巖開采過程中，采空區(qū)上方應力的不斷變化造成局部應力集中，使上覆巖層發(fā)生移動和破壞。直接頂是由于張力作用產(chǎn)生豎向裂隙與離層，進而冒落?；卷攧t是由于豎向裂隙的不斷發(fā)展并逐漸貫通而形成的“砌體梁”結構，達到承載力上限引起破壞；在采動過程中，巖體的破壞結構面隨著垮落穩(wěn)定又逐漸閉合，因此離層面和縱向破斷面的張開和閉合是一個動態(tài)演化過程。

（2）模擬結果顯示了細中砂組合覆巖采動中存在的特點：8煤層模擬區(qū)域內直接頂每步開挖都會發(fā)生應力擾動，在40m步距時離層裂隙極為發(fā)育，達55m步距時冒落變得劇烈?；卷?shù)某醮慰迓洳骄酁?5m，周期垮落步距為10～15m。離層裂隙發(fā)育高度約75m左右，形成穩(wěn)定的冒落帶。

［1］Dziurzynski W，Krach A.Mathematicalmodel of Methane Emis?sion Caused by a Collapse of Rockmass Crump［J］.Archives of Mining Sciences，2001，46（4）：433-449.

［2］石必明，劉澤功.保護層開采上覆煤層變形特性數(shù)值模擬［J］.煤炭學報，2008，33（1）：17-22.

［3］王戈，鞠巍，石興龍，等.工作面上覆巖層活動規(guī)律的數(shù)值模擬分析［J］.煤礦安全，2012，43（4）：146-149.

［4］乙體成，鄭樹權，王社新，等.東保衛(wèi)煤礦大傾角煤層頂板活動規(guī)律與綜采試驗研究［J］.礦業(yè)安全與環(huán)保，2013（6）：50-53.

［5］陳愛和.綜采工作面覆巖破斷規(guī)律試驗及高位鉆孔參數(shù)優(yōu)化［J］.礦業(yè)安全與環(huán)保，2014（5）：47-50.

［6］錢鳴高，劉聽成.礦山壓力及控制［M］.北京：煤炭工業(yè)出版社，1991：79-107.

［7］劉俊.基于板模型的頂板力學分析及應用［D］.徐州：中國礦業(yè)大學力學與建筑工程學院，2008.

［8］浦海，黃耀華，陳榮華.采場頂板X-O型斷裂形態(tài)力學分析［J］.中國礦業(yè)大學學報，2011，40（6）：835-840.

［9］陳忠輝，謝和平，林忠明.綜放開采頂煤冒放性的損傷力學分析［J］.巖石力學與工程學報，2002，21（8）：1136-1140.

［10］馬占國，黃偉，郭廣禮，等.覆巖失穩(wěn)破壞的時變邊界邊界力學分析［J］.遼寧工程技術大學學報，2006，25（4）：515-517.

［11］汪宏志，孫林，胥翔，等.基于分形理論的蘆嶺煤礦8煤層頂板斷裂構造復雜程度研究［J］.中國煤炭地質，2012，24（6）：22-25.

［12］隋旺華.開采覆巖破壞工程地質預測的理論與實踐［J］.工程地質學報，1994，2（2）：29-37.

［13］錢鳴高.采場圍巖控制理論與實踐［J］.礦山壓力與頂板管理，1999，3（4）：12-15.

［14］宋世杰，趙曉光，王雙明，等.覆巖內砂巖層數(shù)對開采沉陷的影響分析與數(shù)值模擬［J］.礦業(yè)安全與環(huán)保，2014（5）：9-12.

［15］TANG C A，KAISER P K.Numerical Simulation of Cumulative Damage and Seismic Energy Release During Brittle Rock Failure［J］.Fundamental J.Int.Rock And Min.Sci.：Part I，1998 （2）：113-121.

［16］黃明利，唐春安.非均勻因素對I型裂紋擴展、相互作用影響的數(shù)值分析［J］.巖石力學與工程學報，2002，21（8）：1111-1114.

TD325

A

1004-5716（2016）04-0126-04

2015-03-30

鄒立強（1986-），男（漢族），山東濟南人，助理工程師，現(xiàn)從事地質工程和礦井地質方向研究工作。