基于三维离散元法的煤层综放开采动态应力分析.docxVIP

基于三维离散元法的煤层综放开采动态应力分析 1 工作面开采中的覆岩运动及其分布规律 近年来，随着我国一些煤矿逐渐向大型煤矿、大采高采等高强度开采方式的转型，以及采空区的快速形成，导致了大面积采空区的快速形成。此外，采空区顶板的开发设计和开采顺序不足。尤其是具有深部影响倾向的坚硬顶板运动及其动力负荷是充分考虑屋顶的主要动力因素[1.3]。例如，2007年兖州南屯煤矿的“3.13”屋顶损坏，当采空区的高度接近工作面长度时，由于采空区高大岩石的突然崩溃，导致工作面中道路的荷载破坏。为此,国内外学者[4~11]对于煤层开采顶板岩层的运动以及应力分布规律开展了相关研究工作,并且取得了一定的成果,但是由于受研究手段的限制,对于采空区覆岩运动形态以及采空区岩体中应力变化规律研究得不够深入。理论研究和生产实践都证明,冲击地压、矿震以及顶板事故几乎都与采掘工作面推进过程中覆岩运动破坏和采动围岩应力场应力大小及分布密切联系在一起。然而,煤层的开采,尤其是工作面不断推进过程中,采空区顶板岩体在形态和结构上均呈现出强烈的不连续性,其运动演化过程是几何或材料非线性的结果,离散元法是20世纪70年代初兴起的一种数值计算方法,特别适用于非连续介质、节理岩体的应力分析,在采矿工程、隧道工程、边坡工程以及放矿力学等方面都有一定的应用[13~17]。本文选取国内典型的冲击地压矿井华丰煤矿四层煤的煤岩地质赋存条件,采用美国Itasca公司研制的大型三维离散元软件3DEC建立煤层综放开采模型,实现综放开采采场、采空区顶板岩层在工作面推进过程中的运动规律以及采动围岩(包括不同层位)中采岩层块体应力动态演化特征的研究,其研究结果形成规律性认识以指导生产。 2 数值模拟 2.1 岩块的离散元分析 离散元法引入我国相对较晚,在岩土工程界的有些学术论文中将离散元法(discrete element method)、个别元素法(distinct element code)、不连续变形分析法(discontinuous deformation analysis)等混为一谈。1992年Cundall和Hart认为只有满足下面2个条件的计算程序才可称之为离散元法: (1)允许块体产生有限位移和旋转,并允许块体间完全分离。 (2)在运算过程中,必须能自动判别各块体间的接触点。 同年,Cundall和Hart指出符合离散元定义的程序主要有: (1)个别元素法程序(distinct element programs)。其代表程序有TRUBAL,UDEC,3DEC,DIBS,3DSHEAR和PFC。这些程序使用显式的,与时俱进的方式直接解答运动等式。块体可以是刚性的,也可以是变形的(子划分成单元),接触也可以变形。 (2)模态法。其代表程序为CICE。这种方法与刚性体的离散元很相似,但是对于变形块体使用了模态叠加。 (3)不连续变形分析法。其对应的计算程序便为DDA。这种方法,接触是刚性的,块体可以是刚性的也可以是变形的,通过迭代格式来取得非相互贯通的条件,块体变形来自与应变模式的叠加。 (4)惯量–交换法。这种方法接触和块体都是刚性的,在瞬间的碰撞过程中2个接触的块体交换惯量,可以表示摩擦滑动。 3DEC将岩体视为由许多完整岩块所组成,各完整岩块间由岩体中的不连续面分隔,而各完整岩块与岩块间的接触面视为岩块的边界。完整岩块可被模拟成刚体(rigid block)或可变形体(deformable block),3DEC在模拟可变形岩块时,将岩块自动分割成许多次级块体(sub-block),每个次级块体可配合所选用的材料组成律及外力作用情况,计算岩块的受力及应力分布情况。在节理的模拟方面,主要根据位移–作用力法则,计算岩块在节理面上的剪应力及正应力,作为岩块的应力边界条件,因此可模拟岩块大位移与转动的情况。 2.2 煤岩体力学参数的选取 数值模拟计算所用煤岩物理力学参数主要来自于实验室试验。同时,由于自然界岩体富含节理和弱面,其力学参数往往比实验室所用岩块的力学参数小很多。因此,煤岩体的力学参数取试验所得数据的1/2~1/20。数值模拟计算所用参数如表1,2所示。 2.3 煤岩块的运动 3DEC真正体现了自然界煤岩体处于真三轴的应力状态。模型具有6个自由面,水平方向和底部边界采用取薄层的方法进行速度为0固定。模型为单一煤层开采,煤层埋深为670 m,模型顶部之上顶板按等效载荷施加到模型顶部边界上,等效载荷按下式计算: 式中:q为等效载荷(N/m2);h为未模拟煤岩层厚度(m);ρ为未模拟煤岩层密度(kg/m3);g为重力加速度,一般取9.8 m/s2。 各单元在运动时一般不发生弹跳,这是由于运动时其动能转化为热能而耗散掉的缘故,因此煤岩块的运动是不可逆的过程。为避免岩块在平衡位置振动,应加入