基于RFPA2D下覆岩采动裂隙发育特征分析和应用.docVIP

基于RFPA2D下覆岩采动裂隙发育特征分析和应用 　　摘 要：为了更好地了解顶板岩体的工程地质特征，对研究区8煤层顶板岩石沉积环境进行了综合分析，提出研究区8煤层为一套高建设性浅水三角洲沉积体系及沉积相序，并利用全过程破裂软件RFPA2D对8煤层顶板开采下覆岩采动裂隙发育特征进行了模拟，这不仅为采区的初步设计提供了参考，也为后期巷道的合理设计支护提供了依据。 　　关键词：浅水三角洲；数值模拟；采动裂隙；顶板稳定性 　　1 序言 　　地下煤岩体被采动后，岩体周围的平衡应力状态被破坏，致使其发生变形甚至破坏。这些变化直接或间接的导致了顶板垮塌、煤与瓦斯突出、冲击地压、矿井透水、突水等问题的发生[1]。其中顶板事故的发生频率较高，约占事故总起数的54.42%[2]。这些问题对煤矿安全生产造成了很大的障碍，给我国人民的生产生活造成了很大损失。 　　芦岭煤矿位于安徽省宿州市东南20余公里处，随着经济的不断发展和浅层煤炭资源的逐渐枯竭，其开采深度也由浅水平向深水平延伸过度，煤层顶板问题的日趋严重直接困扰着该矿的正常生产。该矿8煤层位于下石盒子组下部，上距7煤层约20m，下距9煤层约3.5m，该段岩石主要为砂岩、粉砂岩、泥岩和煤。通过分析钻孔资料与生产资料知：该矿8煤层的直接顶板以泥岩为主，局部含细砂岩、粉砂岩、中砂岩及炭质泥岩，伪顶以砂质泥岩或泥岩为主，呈薄层状，局部含少量细砂岩，植物根茎化石碎片发育。 　　为了更加清楚的了解芦岭煤矿8煤层顶板工程地质特征，确保采取合理的顶板管理方法和科学的支护选型，在对沉积环境充分分析的基础上，对在采动效应下工作面顶板冒裂进行了数值模拟。这对于揭示顶板稳定性和相关工程力学问题的内在机制具有理论和实际应用意义，也为煤炭开采中顶板岩层控制提供了理论和实验依据[3]。 　　2 沉积环境分析 　　晚石炭世，华北含煤沉积盆地南带一度频繁地发生海进海退作用，沉积了一套含大量海相动物化石的浅海碳酸盐、潮坪及煤层互层的海陆交互相沉积。这很符合浅水三角洲形成的基本条件即地形十分平缓，基底整体一致下降[4-5]。盆地北部构造板块在早二叠世晚期快速上抬，致使南北高差明显加剧，由于重力的作用该浅海碳酸盐沉积盆地沉积了大量的陆源碎屑，形成了一套河道高建设性三角洲沉积建造。 　　在平面上砂体的主要分布形态为带状、分叉状分布，在剖面上的分布形态为透镜状、条带状分布（见图1-2），且顶板岩石中含少量植物化石，潮汐层理大量发育，结合测井曲线（见图3）分析知早二叠世晚期，8煤形成之后海水入侵使本区变为浅水海湾环境，在其上沉积了高建设性浅水三角洲沉积体系及其沉积相序，即8煤层的沉积环境。 　　在此基础上，总结了8煤层顶板五类不同的煤岩相组合类型，为后期顶板分类提供了依据，即：（1）煤-老顶组合类型；（2）煤-页岩-老顶组合类型；（3）煤-细砂岩-老顶组合类型；（4）煤-泥岩-老顶组合类型；（5）煤-粉砂岩、泥岩-老顶组合类型。 　　3 数值模型 　　大量生产数据表明，煤体经采空以后，煤层顶板附近的围岩会发生岩层弯曲、下沉、甚至坍落等现象。在连续介质、损伤介质力学原理的基础之上，开发了RFPA2D软件，它可以直观、有效地模拟岩石破裂全过程，并且具有应力分析和破坏分析两方面的功能。本次选用了第一类煤岩相组合类型中标志性较强、岩层分布明显、岩性描述准确、构造简单的2000-1孔，借助RFPA2D数值计算软件来模拟分析采动过程中围岩的变化规律，以及岩体变形、破坏特征。 　　3.1 顶底板岩石特征 　　模拟实例位于煤层顶板主要为伪顶发育区，其伪顶为0.1m～0.3m，老顶主要有泥岩及煤层组成，岩层呈层状分布，且容易发生离层现象，第一关键层厚度较薄。考虑到实际软件模拟需要，将岩层厚度统一四舍五入为整数，对于极薄的伪顶和其它岩层在模型中仅以相应特征的薄分层代替（见表1）。 　　3.2 计算模型与边界条件 　　根据8煤层分布特征模型尺寸为水平280m×垂 直130m，整个模型共划分为 36400 （280*130）个细观单元模型竖向范围。对应现场为：顶部距8煤层顶板116m，底部距8煤层底板6m，由于省略了现场顶部几百米厚覆岩产生的自重应力，因此在模型上方施加了一等效荷载。模型每一步的开挖长度为5m，高为2m。考虑到生产要求，需在模型左右两端分别预留一定的的煤柱。并使模型两端产生一定程度的水平约束，且底部固定。图形的灰度值代表岩层的弹性模量，岩层弹性模量越大，亮度越高。计算模型如图1所示。 　　4 结果分析 　　当工作面沿切眼开挖10m时，由于上覆岩层失去了下部岩层的支撑，使直接顶及上覆岩层之间皆出现了不同程度的裂隙，部分顶板发生了冒落。 　　当工作面推进35m时，第一关键层发生弯曲下沉，形成了垂直或斜交于岩层层面的