徐旭東，南瑩浩，邵小朋

(華北科技學(xué)院安全工程學(xué)院，北京東燕郊 065201)

傾斜煤層淺部開采裂隙帶發(fā)育高度數(shù)值模擬研究

徐旭東，南瑩浩，邵小朋

(華北科技學(xué)院安全工程學(xué)院，北京東燕郊 065201)

為保證傾斜煤層安全開采的需求，以許廠煤礦330采區(qū)3336工作面為研究地質(zhì)單元，在以往的工作面實際情況和研究成果的基礎(chǔ)上，以傾斜煤層圍巖力學(xué)為研究角度，利用FLAC3D建立三維開采模型，通過數(shù)值模擬研究了開挖距離、采高對導(dǎo)水裂隙帶發(fā)育高度的影響，結(jié)果表明：隨開挖距離和采高的增加，覆巖導(dǎo)水裂隙帶高度也增加，但開挖距離對覆巖導(dǎo)水裂隙帶高度起主導(dǎo)作用。當(dāng)開挖距離為200m時，頂板破壞深度達到最大值43m。

傾斜煤層;許廠煤礦;數(shù)值模擬;裂隙帶;破壞深度

Abstract：In order to guarantee the safe exploitation of inclined coal seam，taking the 3336 working face of No.330 mining area of Xuchang Coal Mine as research geological unit，on the basis of the actual situation and research results in the past，taking the surrounding rock mechanics of inclined coal seam as the research angle，using FLAC3D to build 3D mining model，the influence of excavation distance and mining height is studied on the development of fractured zone through numerical simulation.The results show that with the increase of excavation distance and mining height，the height of water induced fracture zone increases，but the excavation distance plays a dominant role in the height of fracture zone of overlying rock.When the excavation distance is 200m，the roof failure depth reaches the maximum of 43m.

Key words：inclined seam;Xuchang Coal Mine;numerical simulation;fracture zone;failure depth

0 引言

對于我國水下采煤的專題研究，主要通過經(jīng)驗公式、相似模擬試驗結(jié)合現(xiàn)場測試以及數(shù)值模擬等方法確定導(dǎo)水裂隙帶高度［1－2］。通過大量現(xiàn)場試驗，得出了適用于分層綜采及普采的覆巖“兩帶”高度計算的經(jīng)驗公式。但在淺部傾斜煤層對“兩帶”高度發(fā)育的影響方面，目前研究的還較少。“兩帶”經(jīng)驗公式已經(jīng)不能滿足許廠煤礦330采區(qū)傾斜煤層安全開采需求。在其他學(xué)者的研究成果上，利用FLAC3D數(shù)值模擬軟件，綜合運用包括水文地質(zhì)學(xué)、巖石力學(xué)、地質(zhì)學(xué)、采礦學(xué)等方法和理論，對3煤層頂板導(dǎo)水裂隙帶進行研究分析及數(shù)值模擬，對于煤礦安全生產(chǎn)具有不可小覷的指導(dǎo)意義［3－5］。

1 井田概況

許廠煤礦地處濟寧城區(qū)東北部，礦井平面面積約56.587km2。濟寧煤田許廠煤礦地層區(qū)屬華北地層區(qū)——魯西地層分區(qū)，區(qū)域上主要發(fā)育早古生代、中生代、新生代地層，含煤巖系為古生代石炭——二疊紀(jì)月門溝群山西組、太原組。井田內(nèi)地層自下而上依次發(fā)育奧陶紀(jì)馬家溝群，石炭——二疊紀(jì)月門溝群本溪組、太原組和山西組，二疊紀(jì)石盒子群，侏羅紀(jì)淄博群三臺組及第四系。研究區(qū)3336工作面位于330東翼采區(qū)，所采煤層為3下煤，影響掘進的含水層主要為3下煤層頂板砂巖含水層、第四系孔隙含水層，粉細砂巖和粉砂巖為煤層頂?shù)装宓闹饕獛r性。

2 模型的建立及邊界條件

2.1 巖石物理力學(xué)性質(zhì)

根據(jù)許廠礦現(xiàn)有的鉆孔資料和巖土測試數(shù)據(jù)，如表1所示，本次模擬的輸入?yún)?shù)，需要選取包括體積模量、切變模量、抗拉強度、抗剪強度、內(nèi)摩擦角、內(nèi)聚力、密度等各煤巖層的巖石物理力學(xué)性質(zhì)，進行網(wǎng)格劃分。

表1 巖石物理力學(xué)性質(zhì)

2.2 邊界條件

根據(jù)3336工作面的實際情況，確定該模型的邊界約束條件：

(1)上部邊界條件：通常情況下，數(shù)值計算時，將煤層上覆巖層的重量設(shè)定為上部邊界條件

式中：γ—上覆巖層的平均體積力，kN/m3;

H—煤層的埋深，m。

單位體積上覆巖層重力取27 kN/m3，3#煤層垂直深度平均值約為200m，代入式得：

(2)下部邊界條件：本模型的底板作為下部邊界條件，簡化成為位移邊界條件，z方向簡化為固定鉸支座，在x、y方向視為可以運動，即v=0。

(3)兩側(cè)邊界條件：本模型中的實體煤巖體作為兩側(cè)邊界條件，簡化成為位移邊界條件，其他方向簡化為固定鉸支座，在z方向視為可以運動，即u=w=0。

2.3 模型建立

FLAC3D模擬計算采用摩爾—庫倫本構(gòu)模型，運用動態(tài)運動方程進行求解，保證其模擬物理上的運動過程沒有數(shù)值障礙［6－7］。以 3336工作面為原型，如圖1，長(x方向)寬(y方向)高(z方向)分別為300 m、300 m、180 m，施加上覆巖層自重應(yīng)力在模型頂部，3煤層采高為3m。

與巖石的抗壓強度相比，其抗拉與抗剪強度要低很多，因此巖石內(nèi)部剪張力作用是巖石破壞的主要原因，“塑性破壞”可被認(rèn)為是巖體破壞的主要形式。研究中，把彈性變形和塑性變形的分界點作為巖石的屈服極限，巖石在達到屈服極限之前，巖體變形為彈性變形，屈服極限之后，巖體變形為塑性變形。通常采用屈服準(zhǔn)則判斷某點的應(yīng)力是否進入塑性狀態(tài)，本次模型采用摩爾—庫侖屈服準(zhǔn)則進行數(shù)值模擬。建立計算模型如圖1所示：

圖1 三維計算模型

2.4 模擬計算過程

本次模擬的計算過程如下：

(1)對3#煤層進行開挖時，按照每20m為開挖的步距，沿煤層的走向進行一次性開采;

(2)為保護煤柱，兩側(cè)應(yīng)留設(shè)50 m的邊界，以消除邊界的影響;

(3)承受的水平壓力為垂向應(yīng)力乘以側(cè)壓系數(shù)0.75。

2.5 數(shù)值模擬結(jié)果與分析

在數(shù)值模擬的計算中，對3煤層進行開挖時，按照每20m為開挖的步距，模擬出分步開挖后煤 巖層塑性破壞區(qū)分布圖，如圖2所示。

圖2 工作面推進不同距離時塑性區(qū)變化圖

由于煤層的開采導(dǎo)致采空區(qū)頂板受到自重及上覆巖層重力的作用，巖層向下位移和彎曲，巖層受拉應(yīng)力作用，當(dāng)拉應(yīng)力大于其極限強度時，直接頂首先發(fā)生斷裂，然后破碎，最終發(fā)生冒落的過程。

由于煤層的開采導(dǎo)致采空區(qū)周圍巖石發(fā)生位移和變形，該過程較復(fù)雜，對于煤層上覆巖層，從煤層的直接頂巖開始，其移動過程由上而下依次發(fā)生冒落、斷裂、離層、裂隙等。

根據(jù)數(shù)值模擬結(jié)果，可得煤層在開挖的各個階段的裂隙帶高度，如表2所示

表2 煤層在開挖各個階段的頂板破壞深度

從表2可以看出，3煤頂板的破壞深度隨著開挖距離的不斷增加而增加，但當(dāng)開挖距離達到150m后，頂板破壞深度基本不變，變化不大，當(dāng)開挖距離為200 m時，頂板破壞深度到達最大值43m。因此可推測頂板最大破壞深度在43m。

3 結(jié)論

通過對許廠煤礦3336工作面采用FLAC3D數(shù)值模擬軟件模擬3煤層開采時導(dǎo)水裂隙帶發(fā)育情況，并根據(jù)模擬結(jié)果，得出如下結(jié)論：

(1)工作面開挖距離150m是導(dǎo)水裂隙帶發(fā)育的一個轉(zhuǎn)折點，是一個從發(fā)育充分、快速向發(fā)育緩慢、循序漸進的分界點;當(dāng)工作面開挖距離為200m時得到最大的頂板破壞深度43 m，由此可推斷出該深度為最大破壞深度。

(2)裂隙帶呈現(xiàn)“馬鞍型”，表現(xiàn)為中間較低、兩邊較高的形態(tài)，導(dǎo)水裂隙帶起始邊界位于煤壁附近;從裂隙帶的分布區(qū)域來看，采場周邊為裂隙帶高度發(fā)育的較大區(qū)域，且該區(qū)域受到拉伸破壞的影響導(dǎo)致更易產(chǎn)生導(dǎo)水裂隙，因此采場邊界的突水風(fēng)險性應(yīng)值得重視。

(3)僅從工作面的長度以及采高兩個角度研究對首采面導(dǎo)水裂隙帶發(fā)育高度的影響程度，其他可能的影響因素，還有待進一步研究，如區(qū)段煤柱寬度、大面積開采、開采時間、采空面積、埋深等。

(4)采用數(shù)值模擬不僅可很大程度上減少實際測量所需的人力物力財力，且能給予煤礦企業(yè)實際開采以指導(dǎo)，對煤礦安全生產(chǎn)有著不可小覷的意義。

［1］ 王錦山.孤島工作面上覆巖層變形破壞規(guī)律相似試驗研究［J］.山東科技大學(xué)學(xué)報(自然科學(xué)版)，2011，30(5)：6－11.

［2］ 王慶照，蔣升，司馬俊杰.厚煤層重復(fù)采動覆巖破裂發(fā)育規(guī)律研究［J］.山東科技大學(xué)學(xué)報(自然科學(xué)版)，2010，29(4)：67－72.

［3］ 杜元洲，杜何辛.覆巖破壞機理和導(dǎo)水裂隙帶發(fā)育高度研究廠［J］.煤礦安全，2011，31(04)：70－72.

［4］ 李學(xué)良.基于FLAC3D的采動區(qū)覆巖破壞高度數(shù)值模擬研究［J］.煤炭技術(shù)，2012，31(10)：83－85.

［5］ Pariseau W G.App lication of finite element analysis to miming engineering［M］.Hudson J A.Comp rehensive Rock Engineering.Oxford：Pergamon Press，1993：491－522.

［6］ 馬亞杰，武強，章之燕，等.煤層開采頂板導(dǎo)水裂隙帶高度預(yù)測研究［J］. 煤炭科學(xué)技術(shù)，2008，36(5)：59－62.

［7］ 劉增輝，楊本水.利用數(shù)值模擬方法確定導(dǎo)水裂隙帶發(fā)育高度［J］.礦業(yè)安全與環(huán)保，2006，33(5)：16－19.

Numerical simulation study on development height of fracture zone in shallow coal seam

XU Xu－dong，NAN Ying－h(huán)ao，SHAO Xiao－peng
(College of Safety Engineering，North China Institute of Science and Technology，Yanjiao，065201，China)

TD741

A

1672－7169(2017)03－0009－05

2017－04－17

中央高?；究蒲袠I(yè)務(wù)費資助(3142017036)

徐旭東(1993－)，男，江蘇南京人，華北科技學(xué)院在讀碩士研究生，研究方向：礦井水害防治。E－mail：1256057641@qq.com