条带开采对煤层稳定性影响的数值模拟研究.docxVIP

条带开采对煤层稳定性影响的数值模拟研究 0 地表沉陷控制策略 煤层地下开采破坏了岩体现有的应力平衡，并将岩屑表面周围的岩石移动并变形。随着开采的进行,采空区围岩某些部位达到弹塑性极限应变状态并且进入流动状态,此时可称之为破坏或失稳,导致井巷工程及地表建筑物的破坏。目前控制岩层及地表沉陷的方法主要有充填开采、部分开采(条带或房柱式开采)、覆岩离层注浆等。尽管条带开采的资源回收率相对较低,但由于它能有效地控制上覆岩层和地表沉陷,不需要增加生产成本,生产管理简单,已在我国“三下”采煤中被广泛应用,并取得了显著效果。在保证煤层开采整体稳定的前提下尽可能提高采出率是当前迫切需要解决的实际问题。研究煤层开采中不同开采方案对地表沉陷、采空区围岩的影响,并用安全系数进行稳定性的评价具有重要意义。 1 基于强度法确定系数的原理和实现 1.1 抗剪储备验算 自1975年,Zienkiewicz等首次提出了抗剪强度折减系数概念以来,由此确定的强度储备安全系数得到了岩土界的广泛认同,并在边坡稳定性分析中得到普遍采用。强度折减的原理是在岩土弹塑性数值计算中折减岩土体的抗剪切强度参数,使岩土体达到极限破坏状态。抗剪强度折减系数定义如下:在外荷载保持不变的情况下,岩土材料所能发挥的最大抗剪强度与外荷载在岩土体内所产生的实际剪应力之比,其公式如下: 式中:τf为岩土材料的最大抗剪强度;τ为实际剪应力。 外荷载所产生的实际剪应力应与抵御外荷载所发挥的最大抗剪强度相等,即按照实际强度计算所确定的,实际中得以发挥的抗剪强度相等。当假定岩土材料各部分抗剪强度的发挥程度相同时,这种抗剪储备安全系数通常就称之为评价工程的整体稳定的安全系数Fs。因此,该方法所确定的安全系数实际上反映的是岩体结构的强度裕度。 在煤矿开采过程中,煤层开挖卸荷的影响使得围岩反复进行应力调整,形成一种在采空区周围岩体径向应力相对较小而切向应力相对较大的二次应力状态,在采空区围岩内部形成一个与采空区空间结构密切相关的应力场。煤层开挖实践表明,采空区围岩失稳最直观的体现是采空区围岩出现大变形并发展成岩体垮落,其力学机制是不利应力状态使得岩体发生剪切屈服并出现塑性流动。故式(1)给出的强度储备安全系数同样可适用于描述煤层开采的整体稳定性。 1.2 采空区至地表岩体塑性区域连通的数值问题 理论上讲,边坡或隧道以及矿山开采等地下工程,当采用强度折减法时,围岩达到弹塑性极限应变状态并且进入无限流动状态时方可称之为破坏或失稳。目前,关于强度折减法中破坏的判别标准大致分为以下2类:(1)采用数值计算不收敛作为破坏或失稳的标志;(2)采用广义剪应变、广义塑性应变或者等效塑性应变区域贯通全断面作为破坏或失稳的标志。 在边坡工程中经常采用数值计算不收敛作为破坏或失稳的标志,这在数值实现上非常方便,在均质边坡的整体稳定计算中也被证明是符合实际情况的。但是,通过笔者的实际计算,在煤层地下开采中,当计算结果显示采空区之间及采空区至地表岩体塑性区域连通,并且地表最大沉降和采空区顶板最大位移突然增大时,数值计算收敛性仍然很好。因而,在煤层地下开采的模拟中,采用数值计算不收敛作为破坏或失稳的标志是不合理的。 另外,通过研究发现,塑性区域贯通仅是破坏的必要条件而非充分条件。在强度折减法计算过程中,塑性区的发展是一个连续的过程。开始时,塑性区贯通全断面,塑性应变值较小。随着强度参数的继续折减,塑性应变值继续增大,当塑性应变值和位移突然增大时,表示此时岩土体沿破坏面发生无限流动,即围岩破坏状态。因此,笔者定义煤层条带开采整体失稳的临界状态判断标准为采空区之间及采空区至地表岩体塑性区域连通,并以地表最大沉降和采空区顶板最大位移突然增大作为临界状态判断标准的验证条件。 值得注意的是,地下采煤工程和边坡工程具有不同的工程特点:首先,在多数煤层地下开采过程中,岩土体虽然会出现贯通的塑性区域,但沿破坏面发生的流动不会是无限的。其次,煤层地下开采中岩层多为层状分布,各岩层的物理、力学性质往往差别很大,无法简化为均质岩体。因为地下采煤工程和边坡工程不同的工程特点,使强度折减法在地下采煤工程中的应用表现出自身的特点。 1.3 mohs-coolomb准则 FLAC3D采用拉格朗日算法,考虑变形对结点坐标的影响,适合于建立大变形非线性模型,而且它以运动方程为基本差分方程。这样就能够模拟材料和结构逐渐稳定,进而发生塑性破坏,最后又逐渐稳定的动态过程,能够反映采矿工程中岩体受采动影响发生变形与破坏的过程。 屈服准则描述了不同应力状态下材料某点进入塑性状态,并使塑性变形继续发展所必须满足的条件。Mohr-Coulomb准则作为一种描述摩擦类岩土体的剪切强度准则,是目前岩土力学研究中应用最为广泛的屈服准则,其表达式为 式中:I1,J2分别为应力张量第1不变量