概率积分法在煤矿采空区地表变形评价中应用.docVIP

概率积分法在煤矿采空区地表变形评价中应用 　　【摘 要】煤矿开采沉陷是由于煤矿地下煤层的开采导致的一种主要的地质灾害类型，对土地的破坏非常严重。本文通过对宁东某煤矿地质采矿条件的研究，利用概率积分法计算评价开采引起的地表变形及稳定性，为矿山企业安全生产及地质灾害防治提供科学参考依据。 【关键词】开采沉陷；地质灾害；概率积分法；地表变形；MathCAD软件 0 引言 煤炭从地下采出后，其上方覆盖岩层失去支撑，岩体内部的原有应力平衡状态受到破坏，引起岩体内应力的重新分布，使采空区周围的岩体产生位移、变形直至破坏。当采空区面积达到一定的范围后，引起采空区上方地表产生移动与变形，从而引发一种主要地质灾害类型，对矿区地质环境的破坏非常严重[1-2]，其主要表现为：地面塌陷及地裂缝使得地表耕地破坏；干旱山区地表水系受到破坏；平原地区大片地表移动形成盆地面积水；受采动影响的山地产生滑坡，危及工业及民用建筑物和生活生产设施的安全；道路塌陷、房屋变形破坏等[3-4]。 开采沉陷的预计方法主要有概率积分法、负指数函数法、典型曲线法和数值计算法（有限单元法、边界单元法和离散单元法等）。概率积分法源于上世纪50年代波兰学者李特威尼申引入岩层移动领域的随机介质理论，后由刘宝琛、廖国华等学者引入国内，经过几十年的研究，目前已成为我国较成熟的、应用最为广泛的预计方法之一[5]。 本文针对宁东某煤矿开发利用方案，通过研究煤矿采空塌陷与煤层厚度、倾角、采深、采厚、上覆岩层性质及开采方式等，建立计算模型并采用概率积分法计算走向及倾斜主断面上地表移动与变形，为该矿的安全生产和地质灾害防治提供可靠的科学依据。 1 地表移动变形预计[5-8] 1.1 半无限开时主断面上地表移动与变形的预计 地表移动观测的实践表明：当开采宽度达到相当大以后，再继续增加开采宽度，对地表下沉的最大值和停采线上方的地表变形几乎没有影响。将这种开采情况数学抽象，即认为停采线一侧煤层未被采动，而采区另一侧的煤层已全部采空，这就是所谓的半无限开采。 1.1.1 走向主断面 W（x）=■[erf（■）+1] T（x）=T■e■ K（x）=2？仔■（-■）e■ U（x）=bW■e■ ε（x）=2？仔b■（-■）e■ （1） 式中：W（x）、T（x）、K（x）、U（x）、ε（x）分别为主断面上下沉、倾斜、曲率、水平移动和水平变形值；W■、T■、K■、U■、ε■分别为最大下沉、最大倾斜、最大曲率、最大水平移动和最大水平变形，可按如下公式计算： W■=qmcosα、T■=W■/r、K■=1.52W■/r2、U■=bW■、ε■=1.52bW■/r；x――走向主断面上地表点至下沉曲线拐点的距离；e――自然对数的底；erf（x）――误差积分，erf（x）=■■e■dt；r-――主要影响半径；q――下沉系数；m――采厚；α――煤层倾角；b――水平移动系数，下同。 1.1.2 倾斜主断面 倾斜主断面上地表移动与变形的预计公式与式（1）相同，但需以y代替各式中的x。当预计下边界半个盆地移动与变形时，以下山方向的主要影响半径r1代替各式中的r；预计上边界半个盆地移动与变形时，用上山方向的主要影响半径r2代替各式中的r。其中，r1、r2按下式计算： r1=H1/tanβ r2=H2/tanβ 式中，H1――工作面下边界的开采深度；H2――工作面上边界的开采深度。 1.2 有限开采时主断面上地表移动与变形的预计 有限开采是指回采面积未能使地表出现充分移动的开采。 1.2.1 走向主断面 W（x）=Cy[W3（x）-W4（x-l）] T（x）=Cy[T3（x）-T4（x-l）] K（x）=Cy[K3（x）-K4（x-l）] U（x）=Cy[U3（x）-U4（x-l）] ε（x）=Cy[ε3（x）-ε4（x-l）] （2） 式中：W3（x）、T3（x）、K3（x）、U3（x）、ε3（x）分别为工作面从左边界向右侧半无限开采的下沉、倾斜、曲率、水平移动和水平变形；W4（x-l）、T4（x-l）、K4（x-l）、U4（x-l）、ε4（x-l）分别为工作面从右边界向左侧半无限开采的下沉、倾斜、曲率、水平移动和水平变形；Cy――倾斜方向采动程度系数，l――工作面走向方向的计算长度，L――工作面倾斜方向的计算长度，分别按下式计算： Cy =■erf■+1-■erf■+1 l=D■-2S■ L=（D■-S■-S■）■ 其中，D■、D■――工作面倾斜和走向的长度；S■、S■、S2――走向、下山和上山方向的拐点偏移距；θ――开采影响传播角，下同。 1.2.2 倾斜主断