2024年4月28日发(作者：宾得k52s)

第２０卷第１期　

安徽建筑工业学院学报（自然科学版）　

Ｖｏ１．２０　Ｎｏ．１　

２０１２年２月　

Ｊｏｕｒｎａｌ　ｏｆ　Ａｎｈｕｉ　Ｉｎｓｔｉｔｕｔｅ　ｏｆ　Ａｒｃｈｉｔｅｃｔｕｒｅ＆Ｉｎｄｕｓｔｒｙ　

Ｆｅｂ．２０１２　

近距离高瓦斯煤层群首采工作面顶板　

裂隙发育规律研究　

张德万，　肖峻峰，　卢　平　

（安徽建筑工业学院土木工程学院，合肥２３０６０１）　

摘要：确定顶板瓦斯运移规律及其富集区是综合治理高瓦斯工作面卸压瓦斯的关键。本文对我国西部地　

区典型的近距离高瓦斯煤层群４．２ｍ大采高工作面回采过程中顶板裂隙演化特征进行ＦＩ　ＡＣ如模拟，研究得　

出在工作面顶板覆岩ｌＯ倍采高处，纵向裂隙和横向裂隙发育最剧烈且相互导通，形成高瓦斯富集区。此对大　

采高工作面卸压瓦斯的安全高效治理具有理论和实际指导意义。　

关键词：ＦＬＡＣ３。；顶板裂隙；瓦斯运移通道；瓦斯富集区　

中图分类号：ＴＤ７１２．６　文献标识码：Ａ　文章编号：１００６—４５４０（２０１２）０１—０２５—０４　

Ｓｔｕｄｙ　ｏｎ　ｒｏｏｆ　ｃｒａｃｋｓ　ｄｅｖｅｌｏｐｍｅｎｔ　ｌａｗ　ｏｆ　ｔｈｅ　ｆｉｒｓｔ—ｍｉｎｅｄ　ｆａｃｅ　

ｄｕｒｉｎｇ　ｃｌｏｓｅ－ｄｉｓｔａｎｃｅ　ｈｉｇｈ　ｇａｓ　ｃｏａｌ　ｓｅａｍｓ　ｇｒｏｕｐ　

ＺＨＡＮＧ　Ｄｅ－ｗａｎ，ＸＩＡＯ　Ｊｕｎ—ｆｅｎｇ，ＬＵ　Ｐｉｎｇ　

（Ｓｃｈｏｏｌ　ｏＩ　ｃｉｖｉｌ　ｅｎｇｉｎｅｅｒｉｎｇ，Ａｎｈｕｉ　Ｕｎｉｖｅｒｓｉｔｙ　ｏｆ　Ａｒｃｈｉｔｅｃｔｕｒｅ，Ｈｅｆｅｉ，２３０６０１，Ａｎｈｕｉ，ＰＲ　Ｃｈｉｎａ）　

Ａｂｓｔｒａｃｔ：Ｔｈｅ　ｄｅｆｉｎｉｔｉｏｎ　ｏｆ　ｒｏｏｆ　ｃｒａｃｋｓ　ｄｅｖｅｌｏｐｍｅｎｔ　ｌａｗ　ａｎｄ　ｇａｓ　ｃｏｎｃｅｎｔｒａｔｉｏｎ　ｚｏｎｅ　ｉＳ　ｃｒｉｔｉｃａＩ　ｆｏｒ　ｃｏｒｎ—　

ｐｒｅｈｅｎｓｉｖｅ　ｃｏｎｔｒｏｌ　ｏｆ　ｇａｓ　ｗｉｔｈ　ｐｒｅｓｓｕｒｅ　ｒｅｌｉｅｆ　ｉｎ　ｈｉｇｈ　ｇａｓ　ｆａｃｅ．Ｉｎ　ｔｈｉｓ　ｐａｐｅｒ，ｒｏｏｆ　ｃｒａｃｋｓ　ｅｖｏｌｕｔｉｏｎ　

ｃｈａｒａｃｔｅｒ　ｉｎ　４．２ｍ　ｔｈｉｃｋ　ｆｉｒｓｔ—ｍｉｎｅｄ　ｆａｃｅ　ｄｕｒｉｎｇ　ｃｌｏｓｅ－ｄｉｓｔａｎｃｅ　ｈｉｇｈ　ｇａｓ　ｃｏａｌ　ｓｅａｍｓ　ｇｒｏｕｐ　ｉｎ　ｗｅｓｔｅｒｎ　ｒｅ—　

ｇｉｏｎ　ｏｆ　ｏｕｒ　ｃｏｕｎｔｒｙ　ｗａｓ　ｓｔｕｄｉｅｄ　ｂｙ　ＦＬＡＣａ。ｓｏｆｔｗａｒｅ．Ｒｅｓｕｌｔｓ　ｓｈｏｗ　ｔｈａｔ　ｉｔ　ｉｓ　ｉｎ　１　０　ｔｉｍｅｓ　ｍｉｎｉｎｇ　ｈｅｉｇｈｔ　

ｏｆ　ｒｏｏｆ　ｓｔｒａｔａ　ｔｈａｔ　ｔｈｅ　ｌｏｎｇｉｔｕｄｉｎａｌ　ａｎｄ　ｌａｔｅｒａｌ　ｃｒａｃｋｓ　ｄｅｖｅｌｏｐｅｄ　ｖｉｏｌｅｎｔｌｙ。ｃｏｍｍｕｎｉｃａｔｅｄ　ｍｕｔｕａｌｌｙ　ａｎｄ　

ｆｏｒｍｅｄ　ａｂｕｎｄａｎｔ　ｇａｓ　ｃｏｎｃｅｎｔｒａｔｉｏｎ　ｚｏｎｅ．Ｉｔ　ｃａｎ　ｈｅｌｐ　ｔｏ　ｓａｆｅｌｙ　ａｎｄ　ｅｆｆｉｃｉｅｎｔｌｙ　ｃｏｎｔｒｏ１　ｇａｓ　ｗｉｔｈ　ｐｒｅｓｓｕｒｅ　

ｒｅｌｉｅｆ　ｆｒｏｍ　ｌａｒｇｅ　ｍｉｎｉｎｇ　ｈｅｉｇｈｔ　ｆａｃｅ　ｉｎ　ｔｈｅｏｒｅｔｉｃａｌ　ａｎｄ　ｐｒａｃｔｉｃａｌ　ａｓｐｅｃｔ．　

Ｋｅｙ　ｗｏｒｄｓ：ＦＬＡＣⅢ；ｒｏｏｆ　ｃｒａｃｋｓ；ｇａｓ　ｆｌｏｗ　ｐａｓｓａｇｅ；ｒｉｃｈ　ｇａｓ　ｃｏｎｃｅｎｔｒａｔｉｏｎ　ｚｏｎｅ　

我国煤矿大多为高瓦斯煤层群赋存条件，特　围岩应力场重新分布，上覆岩层划分为“三　

征为“三高一低”，８０　９／６以上重大恶性事故与瓦斯　

带”［３　］，根据“砌体梁”理论＿５　］，采空区顶板岩　

有关，目前瓦斯治理仍是世界性难题¨１　］。对于　

层卸压，产生裂隙Ｌ７　］，规则垮落带和断裂带中顶　

４ｍ以上大采高工作面而言，上、下向扰动卸压范　板岩层产生卸压膨胀，存在竖向裂隙发育区，该区　

围更大，邻近层瓦斯大量涌人工作面，瓦斯控制难　域离层裂隙和竖向破断裂隙发育，横向和竖向裂　

度大，因此在厚煤层大采高条件下探寻瓦斯富集　

隙贯通，并和不规则垮落带相连通，为围岩卸压瓦　

区域及其运移规律尤为重要。煤层回采后工作面　

斯和本煤层工作面采空区积聚的瓦斯提供良好的　

收稿日期：２０１１—１２—２８　

基金项目：安徽建筑工业学院硕博科研启动项目（Ｎｏ．２００６１２０１２２）、国家自然科学基金（５０９７４００２、４１１７２１４７）、安徽省教育厅自然科　

学重点项目（ＫＪ２０ｌ０ＡＯ５８０）。　

作者简介：张德万（１９８４一），男，在读硕士研究生，主要研究方向为矿山安全。　

２６　安徽建筑工业学院学报（自然科学版）　第２Ｏ卷　

储集场所ｌ＿９］。　

针对我国西部地区典型的近距离高瓦斯煤层　

群４．２ｍ大采高工作面，本文利用ＦＬＡ　。软件进　

行模拟，通过分析煤层开采后顶板覆岩纵向位移　

场、纵向应力场和横向应力场的变化，得出煤层顶　

板覆岩移动规律和裂隙发育规律。　

１　２４２０７工作面概况　

２４２０７工作面是我国西部地区典型的高瓦　

斯、大采高、松软顶板工作面，隶属山西省华晋焦　

煤有限责任公司沙曲矿。工作面底板标高在＋　

３６Ｏｍ一４－４５０ｍ之间，地面标高为＋８６２ｍ一＋　

１００７ｍ。工作面回采３　＋４　合并煤层，厚度在　

３．８５ｍ一４．３６ｍ之间，平均４．２０ｍ，瓦斯含量　

１１．７２　／ｔ。２　煤层上距３＃＋４　煤层平均２０ｍ，　

平均厚度为１．０４ｍ，瓦斯含量１０．６５ｍ。／ｔ；５　煤层　

下距３　＋４　煤层平均５．６ｍ，平均厚度为３．３ｍ，　

瓦斯含量１２．０６ｍ。／ｔ。工作面伪顶不发育，局部　

有０．２ｍ的泥岩；直接顶为灰黑色中细砂质泥岩，　

厚度５．５ｍ；老顶为灰白色中砂岩，厚度５．５９ｍ；直　

接底为１．１ｍ的灰色中砂岩；老底为２．５ｍ的黑色　

粉砂岩。工作面设计走向长度１６６３ｍ，倾向长度　

２２０ｍ，工作面整体呈单斜构造，煤层走向３３０。，倾　

向ＳＷ，倾角４。一７。，平均倾角５。。　

２　ＦＬＡＣ。Ｄ基本原理　

ＦＬＡＣ。。（３一Ｄｉｍｅｎｓｉｏｎａｌ　Ｆａｓｔ　Ｌａｇｒａｎｇｉａｎ　

Ａｎａｌｙｓｉｓ　ｏｆ　Ｃｏｎｔｉｎｕａ）＿１　ｕ］是用于工程力学计算　

的显式有限差分程序。该程序可模拟土、岩石等　

材料的力学行为。快速拉格朗日分析采用混合离　

散方法，将区域离散为常应变六面体单元的集合　

体，又将每个六面体看作以六面体角点的常应变　

四面体的集合体，应力、应变、节点不平衡力等变　

量均在四面体上进行计算。快速拉格朗日分析用　

速率来求某一时步的单元应变增量，即：　

△％一丢（　，＋　（１）　

有了应变增量，即可由本构方程求出应力增　

量，各时步的应力增量叠加即可得到总应力。作　

用在可变形固体上的方程组为本构关系，或者称　

为应力一应变准则。首先由速度梯度得出应变速　

率，公式如下：　

一　

［警＋警］　（２）　

式中，　为应变率分量；　为速度分量。　

本构关系的形式如下：　

：一Ｍ（ａ　，　，　）　（３）　

式中，Ｍ（）为本构关系的函数形式，ｋ为可能　

出现也可能不出现的一个历史参数，依赖于特定　

的定律；“：一”为“由…代换”。　

一

般说来，非线性的本构定律以增量的形式　

表示，因为应力和应变之间的对应关系并非唯一　

的，上式给出了以前的应力张量和应变（应变增　

量）下对应战略的新估计值。最简单的本构定律　

为各向同性的弹性本构关系：　

：一　４－｛岛（Ｋ一专Ｇ）垂ＫＫ＋２Ｇｅｏ｝Ａｔ　

（４）　

式中，岛为Ｋｒｏｎｅｃｋｅｒ符号；Ａｔ为时间步；Ｇ、Ｋ为　

剪切模量和体积模量。　

３模型的建立及计算　

本模型选择以沙曲矿２４２０７工作面为原型，　

建立受覆盖岩层重力作用即在重力场作用下的三　

维空间力学模型。为便于建模和剖分，同时充分　

体现各岩层组合特征，将研究区内力学性质相近　

的岩层归并为一组。　

坐标系按如下规定：ｚ轴为煤层走向方向，Ｙ轴　

为煤层倾向方向，　轴为竖直方向即重力方向。沿ｚ　

轴方向长６００ｍ，Ｙ轴方向长４００ｍ，　轴方向高１５０ｍ，　

模型底边即ｚ＝０处到地面设定高度为６００ｍ，煤层　

倾角５。，模型自下到上分为１６组岩层，共划分３００００　

个单元体，３３２０１个节点，３　＋４　煤层为模拟的开采　

煤层，且划分为模型的第７层（ｓｑ７）。据上述参数及　

条件，模型网格划分图如图１所示。　

图１模型网格图　

工作面为仰采面，模型从ｚ一５０ｍ处开切眼，　

采高为４．２ｍ，采面倾向长度为２２０ｍ，模拟开采分　

第１期　张德万，等：近距离高瓦斯煤层群首采工作面顸板裂隙发育规律研究　２７　

为６个阶段，开采长度分别为ｌＯｍ，２０ｍ，２０ｍ，　

５０ｍ，５０ｍ，５０ｍ，即从ｘ一５０处开采到ｘ一２５０处，　

累计开采２００ｍ。模型计算采用莫尔一库仑（Ｍｏ－　

ｈｒ－－Ｃｏｕｌｏｍｂ）屈服准则。　

４模拟结果分析　

（１）横向离层裂隙发育区分析　

当工作面向前推进时，采空区上方上覆岩层　

将产生离层裂隙，在ＦＬＡＣ。。数值计算分析过程　

中，纵向位移可真实地反映岩层沉降而产生离层　

裂隙发育状况。图２不同累计进尺下　一１００面　

上的竖向位移云图。　

由图２可看出，在距离切眼５０ｍ（ｚ—ｉ００）的切　

面上，纵向位移场均呈不均匀沉降状态分布；每个阶　

段的开采过程中，该面上的竖向位移最大值分别是　

３．５２ｍｍ、ＩＩ．９８ｍｍ、２７．４７ｍｍ、７７．８６ｍｍ、１０７．５５ｍｍ　

和１４７．４８ｍｍ。随着工作面的推进和采动影响，在竖　

直方向上出现不均匀下沉，沉降量不同的岩层将伴　

随着横向裂隙发育。位移变化趋势与回采推进距离　

直接相关，而位移变化的大小与距开采煤层顶板的　

高度有关。工作面顶板裂隙的发育可以通过工作面　

顶板以上不同倍采高的纵向位移随工作面推进的变　

化进行考察。图３为工作面顶板以上不同倍采高竖　

向位移随工作面推进的变化曲线。　

震ｌ圈｛　

Ｊ　

（ｅ）累计进）￣１５０ｍ，ｘ＝ｌＯ０切面竖向位移云图（Ｄ累计进尺２ｏｏｍ　＝ｌ００切面烃向位移云图　

圈２不同累计进尺下ｚ—１００面上的竖向位移云图　

开采到距离切眼的位置　

冀一　

ｍ　５ｏｍ】００ｍ　１５ｏｍ　２ｏ０ｍ　

－●・５倍聚商　

ｔ７　

ＩＯ倍采商　

ｌ３ｆ　聚高　

掊粟高　

（ａ）距离切ＵＩ￣１０ｍ（ｘ＝６０的切面）　

开采到距离切眼的位置　

ｍ　５ｒ）ｍ　１ｏｎｍ１５０ｍ　２０Ｉ　Ｉｍ　

、　

ｆ：　薷丽　

●７　采　

１０倍聚‘　

一Ｉ３＊采，　

・＿Ｉ５＊　０　

Ｃｏ）距离切Ｉ］￣５０ｍ（ｘ＝１Ｏ０的切面）　

图３工作面顶板以上不同倍采高竖向位移随　

＂１－作面推进的变化曲线　

由图３（ａ）可以看出，在距离切眼ｌＯｍ切面（ｚ　

一６０）上，当工作面累计推进ｌＯｍ时，５—１５倍采　

高处沉降量比较均匀，几乎重合为一条曲线，也相　

对较小，沉降不明显，离层裂隙发育也不明显；当　

工作面累计推进３０ｍ时，不同采高处开始呈不均　

匀沉降，随着工作面的推进，１３、１５倍采高沉降曲　

线趋于重合，位于所有曲线上面，而５、７倍采高的　

下沉曲线趋于重合且，位于所有曲线的下方，１０　

倍采高处纵向位移曲线位于所有曲线中间。图３　

（ｂ）可以看出类似于图３（ａ）的变化趋势。当工作　

面累计进尺到切眼５０ｍ切面处及之前，不同倍才　

高处纵向位移曲线趋于重合；采过切面后，曲线分　

离，不同倍采高下沉不均匀性加剧，５、７倍接近重　

合，１３、１５倍接近重合，１０倍处于两条重合曲线　

中间。　

通过对比图３（ａ）、（ｂ）曲线的变化，可以判断　

在工作面顶板以上１０倍采高处覆岩离层最明显，　

横向裂隙发育最剧烈；当工作面采过切面后，不同　

倍采高处下沉的不均匀性才开始加剧，是因为采　

空区顶板覆岩受采动影响后应力场重新分布而逐　

渐被压实。图４是煤层累计开采２００ｍ时在ｘ一　

６０倾向切面纵向应力云图和ｘ＝＝＝６０倾向切面上　

２８　安徽建筑工业学院学报（自然科学版）　第２０卷　

不同倍采高的竖向应力（ｓｚｚ）的变化情况。由图４　

可看出，随着工作面的逐步推进，采空区顶板覆岩　

应力先释放再回复的变化规律，表明采空区顶板　

・■－ｘ＝６０　

覆岩随开采逐步被压实。　

（ａ）ｘ＝６０倾向切面纵向应力云图　

０．００Ｅ＋００　

开采到距离切眼的位置　

ｍ　３０ｍ　５０ｍ　ｌ００ｍ　ｌ５０ｍ　２０　

．

２．００Ｅ＋０７　

．．－ｌ倍采高　

邑．４００Ｅ＋０７　

＿．卜３倍采高　

６０ｏＥ＋０７　

・　

．．Ｉ５倍采高　

．

ｊ　

画｛．ＯＯＥ＋０７　

／　＝＝Ｉ｝％　一　

・●－７倍采高　

厂　■　—　

・■－ｌ０倍采高　

．

１，００Ｅ＋０７　

￣ｌｍＤｌ３倍采高　

１　２０Ｅ＋０７　

－—－－ｌ５倍采高　

．

１．４０Ｅ＋Ｏ７　

（ｂ）ｘ＝６０ｔ￣向切两纵向应力变化曲线　

图４　ｘ＝６０倾向切面上随工作面推进不同倍采高竖向应力变化　

（２）纵向裂隙发育区分析　

煤层开采后，卸压区内岩层以拉应力为主，当　

其超过岩体的极限抗拉强度时，便出现纵向裂隙，　

因此在顶板覆岩周期断裂时，煤壁前方顶板岩层　

内产生开口向上的纵向裂隙，在煤壁后方顶板岩　

层内产生开口向下的纵向裂隙＿１　。图５为顶板　

以上不同倍采高处水平应力（ｔＯｔ向）随工作面推进　

的变化曲线。可以看出，１３和１５倍采高处，采前　

和采后水平应力（Ｙ方向）始终保持受压状态，不　

会出现纵向裂隙；７倍和１Ｏ倍采高区域，随着工　

作面的推进，水平应力值（倾向）均由负变为正，说　

明应力从原始的受压变状态成采动后的受拉状　

态，且最大为３．８０×１Ｏ　ｐａ，模拟分为１６层，顶板　

７—１０倍各层抗拉强度极限最大为１．８×１Ｏ　ｐａ，　

固会产生纵向裂隙，考虑重力场重新分配，由图４　

可知７采高陪以下将随着工作面的推进而逐渐被　

压实，判定顶板以上７陪到１０倍采高出处于纵　

向裂隙最发育区。　

．．●－　＝８Ｏ　

—●　１００　

・　；１５Ｏ　

－　＝２０Ｏ　

（ａ）顶板以上７倍采高　

开采到距离切眼的位置　

２　Ｏ０Ｅ＋Ｏ６　—————————————————一　

・●・，ｃ＝６０　

．．Ｏａ－＝８０　

。●＿ｒ＝１００　

—●■＝１５０　

￣ｔ－＝２００　

（ｂ）顶板￣ｌ－１０倍采高　

开采到距离切眼的位置　

・・ｘ＝６０　

●ａ＝８０　

‘．－ｘ＝ｌ００　

一　１５０　

・一ｘ＝２００　

（ｃ）顶板￣Ｊ：１３倍采高　

开采到距离切眼的位置　

０　００Ｅ＋００　…～一　一　—　——　——　一…一一——～・　——　

．

５．００Ｅ＋０５１￣ｋｎ　３０ｍ　５０ｍ　１００ｍ　１５０ｍ　２００ｍ　

￣ｘ＝６０　

，＂Ｉｔ－ｘ＝８０　

－●・ｘ＝ｌ０ｏ　

．

■ｒ　

－ｕ－ｘ＝ｌ５０　

＂＇ｘ＝２００　

一～

二二　

（ｄ）项板以上１５倍采商　

图５顶板以上不同倍采高处水平应力（倾向）随＂１－作面　

推进的变化曲线　

５结　论　

按照沙曲矿２４２０７工作面的实际地质条件，　

基于ＦＬＣ３。软件的基本原理，建立了相似模型，进　

行了采动上覆岩层移动规律及裂隙发育的三维快　

速拉格朗日分析，得到如下主要结论：　

（１）竖向位移场均呈现不均匀沉降状态分布，　

不同倍采高处的位移采前下沉基本一致，采后下　

沉不均匀性加剧，并且随着工作面的逐步推进，采　

空区逐步被压实。　（下转第５５页）　

第１期　何俊，等：基坑降水及综合措施在安徽颍上县耿楼枢纽工程中的应用　５５　

专人负责，每日定期观测，做好记录。发现特殊情　

况及时报告。　

降水施工的效果是分不开的。不同的工程有不同　

的地质条件和工程工况，只要我们针对性的，通过　

合理的计算和施工，采用不同的降水方案，就能取　

得很好的效果。　

参考文献　

１龚晓南．深基坑工程设计施工手册Ｉ－Ｍ］．北京：中国建　

测压管管底不高于１３ｍ高程，测压管上部周　

围用粘土封堵，防止漏气。待施工至地面时，拆除　

测压管，孔眼用粘土填实。　

进行降水时在周围设置一定数量的沉降观　

测点。　

４结束语　

本工程施工期为１６个月，经历１个洪水期，　

发生１０年一遇标准以上洪水的可能性是存在的，　

为此必须考虑发生１０年一遇以上洪水的特殊情　

况。虽然节制闸位于旱地施工，但节制闸紧邻沙　

颍河大堤，仍有必要做好安全度汛工作。　

通过理论计算和结合现场实际地质条件和工　

程工况，颍上耿楼枢纽的深井降水施工取得了良　

好的效果，确保了基坑开挖的安全和开挖后的干　

施工。应该说颍上耿楼枢纽的顺利完工，跟深井　

（上接第２８页）　

（２）工作面顶板以上５倍和７倍采高沉降均　

匀，１３倍和１５倍趋于一致，顶板以上１Ｏ倍采高　

处离层明显，横向裂隙最发育。　

（３）工作面顶板以上７—１０倍采高处以拉应　

力为主，在顶板覆岩周期断裂时顶板岩层内产生　

纵向裂隙。　

（４）工作面顶板以上１０倍采高处，纵横向裂　

隙发育剧烈且相互导通，成为良好的瓦斯运移通　

道，属于高瓦斯富集区域。　

参考文献　

１袁亮．高瓦斯矿区复杂地质条件安全高效开采关键　

技术［Ｊ］．煤炭学报，２００６，３１（２）：１７４—１７８．　

２卢平，袁亮，程桦等．低透气性煤层群高瓦斯　

采煤工作面强化抽采卸压瓦斯机理及试验［Ｊ］．煤炭　

学报，２０１０，３５（４）：５８Ｏ一５８５．　

３黄志安，童海方，张英华等．采空区上覆岩层“三带”　

的界定准测和仿真确定［Ｊ］．北京科技大学学报，　

２００６，２８（７）：６０９－６１２．　

４钱鸣高，缪协兴，许家林等．岩层控制的关键层理论　

筑工业出版社，１９９８．　

２刘　勇．浅谈轻型井点降水［Ｊ］．当代建设，２００３　

（５）：４３．　

３　张清国．建筑工程机械［Ｍ］．重庆：重庆大学出版　

社，２００４．　

４　中华人民共和国水利部．水闸施工规范（ＳＬ２７—９１）　

Ｅｓ３．　

５周志芳，郭耿新，汪北华，等．深基坑降水的双层结构　

模型及有限元计算［Ｊ］．河海大学学报（自然科学版），　

２００４（２）：１７９—１８３．　

６郭建堂．水利水电工程施工手册Ｉ－Ｍ］．北京：中国电力　

出版社，２００４．　

Ｉ－Ｍ］．徐州：中国矿业大学出版社，２００３：１７—２８．　

５钱鸣高，刘听成．矿山压力及其控制ＥＭ－］．北京：煤炭　

工业出版社，１９９１：１５０—１５５．　

６钱鸣高，缪协兴，何富连．采场“砌体梁”结构的关键　

块分析ＥＪ］．煤炭学报，１９９４，１９（６）：５５７—５６２．　

７刘泽功，袁亮，戴广龙等．采场覆岩裂隙特征研究及　

在瓦斯抽放中应用［Ｊ］．安徽理工大学学报，２００４，２４　

（４）：１０—１５．　

８刘伟韬，武强，李献忠等．覆岩裂隙带发育高度的　

实测与数值仿真方法研究［Ｊ］．煤炭工程，２００５，　

（１１）：５５～５８．　

９刘泽功，袁亮，李平等．开采煤层顶板环形裂隙　

圈内走向长钻孔法抽采瓦斯研究ＥＪ］．中国工程科学，　

２００４，６（５）：３２—３８．　

１Ｏ陈育民，徐鼎平．ＦＬＡＣ／ＦＬＡＣ３Ｄ基础与工程实例　

［Ｍ］．北京：中国水利水电出版社，２００９：７４—１１１．　

１ｌ李围．隧道及地下工程ＦＬＡＣ解析方法Ｉ－Ｍ］．北　

京：中国水利水电出版社，２００９：９９—２５９．　

１２卢平，刘泽功，廖光煊等．高瓦斯综采面顶板覆岩　

卸压抽放瓦斯试验研究ＥＪ］．力学与实践，２００３，２５　

（４）：５３—５６．