开采损害学课程讲义2Word文档格式.doc

1、属于理论模型方法是建立在力学模型上的，以及建立在弹性或塑性理论基础上的计算方法。在这方面主要有以A.Salstowicz（1958）等为代表的固体力学理论；J.Litwiniszy（1963）等为代表的随机介质理论。建立在弹性或塑性理论基础上的计算方法如：有限单元法（FEM）；边界元法（BEM）；离散元（DEM）等方法；非线性力学（Nonlinear）等方法。目前应用情况简介2.2 概率积分法（重点）目前已成为我国乃至世界范围较为成熟、应用最广泛的预计方法之一。2.2.1 水平成层介质中的单元盆地开采沉陷的随机性随机介质理论为基础非连续介质单元模型，单元相互分离并发生相对运动。如图2.1在三维

2、问题中，地下（x0, y0, z0）处开采使地表点A（x, y, z）附近某一小块面积ds发生下沉这一事件的概率为：（2-10）d（x, y, z）为密度函数。在x-z剖面的z水平上x处的一段岩石条dx内有下沉发生，同时在y-z剖面的同一高度上y处的一段岩石条dy内有下沉发生。因此，ds 发生下沉这一事件的概率为发生上述两事件的概率之积，即：（2-11）f为密度函数，x2、y2对称性。过原点引另一组正交水平轴x、y, 使A点在这一系统中的坐标为（x, y, z ）。在新坐标系中，A点附近的微面ds发生下沉的概率为：（2-12）岩石水平成层f的形态皆一致。微面面积不变ds=ds ；开采点与被考虑

3、的微面相对位置不变P（ds）=P（ds从本质上讲，某一既定的微面在同一开采影响下的下沉概率与坐标轴方向的选择无关。（2-13）从而可得：（2-14）当采用下列方式选择坐标轴，使ox 经过A点，且：（2-15）代入（2-14）得：（2-16）将（2-16）两端对x2、y2取偏微分可得：（2-17）（2-18）由此可得：（2-19）移项得到：（2-20）上式左边为x的函数，右边为y的函数。此方程成立的条件是左右两端都不依赖于自变量x、y，故可令式（2-20）为常量k，从而有：（2-21）将x2看作自变量，解此常微分方程得到：（2-22）式中， P为积分常数。显然，远离中心两端的岩石下沉的概率小。因

4、此，从物理意义说k为负值，令k=-h2代入上式得：（2-23）同理（2-24）将上述结果代入式（2-13）有：（2-25）P（ds） z0煤层的开采以随机的方式传至z水平上的随机分布。由于，由此可得到顶板下沉盆地中的分块下沉体在z水平上造成的微小下沉盆地的表达式：（2-26）下沉盆地体积=（下沉概率采高开采面积） 单元开采分块下沉体，单元盆微下沉盆地。令开采面积Ds=1个微小单位，采高w（x0, y0, z0）=1个微小单位的开采为单元开采。单元开采在上覆岩层中造成的下沉盆地为单元下沉盆地we,于是得单元盆地下沉：（2-27）当然，孤立的下沉盆地并不存在，它们存在于统计意义之中。采出体积足够大

5、地表出现下沉；大面积开采=单元开采；原始的位置弹性变形逐渐垮落传递到地表（“三带”形成过程）下沉盆地。单元盆地是时间函数we=（x, y, z, t），对应的地表下沉盆地体积为： （2-28）Ve也是时间的函数。根据A.Salstowicz假设，下沉盆地体积的增长率与采空区域未压密的体积成正比，即：（2-29）式中 c比例系数，当t=0时，Ve=0；当t时，Ve=1 。解上述方程式可得：（2-30）盆地体积可由上述求得：（2-31）（换元ahx=l，dx=（1/h）dl；hy=l，dy=（1/h）dl；a）则：（与时间有关的概率）（2-32）从而：（2-33）t,（2-34）当宽与高均为一个微

6、小单位，长度为无穷大时的开采称为二维单元开采。当开采长轴方向平行于y轴，此时形成槽形盆地，槽长平行于y轴，而在平行于x轴的任何剖面上，此盆地的形状均相同，称为二维单元盆地，表示为：从而可得水平成层介质中二维单元下沉盆地的最终表达式为：y=0；（-d2）=d2；1/hdl=dd （2-35）2.2.2 单元盆地水平移动假定在单元开采影响下岩石发生的变形很小，从弯曲带直到地表的整个岩层中移动变形是连续的，岩石发生形变但体积不变，即：（2-36）又 （we轴与z轴的指向相反）从而有：式中k为积分常数，取决于边界条件。由上式将we先对z微分再对x积分即可求得ue。式中 ，再对x积分得： 由于对称的原因

7、，在开采中心线上的点不发生水平移动，即：， 故 由此得：（2-37）2.2.3 半无限开采时地表移动盆地走向主断面的移动变形预计所谓半无限开采，如图2.2所示，是指沿工作面推进方向在x =+ ,0区间已被开采，而沿垂直工作面推进方向的开采尺寸足够大，使之达到充分采动。1. 移动变形计算表达式（1）下沉预计表达式 在图2.2中，采深为H，开采厚度为 m，取坐标原点通过开采边界。由于垮落的顶板岩石碎胀的原因，顶板最大下沉量一般小于开采厚度，为：（2-38）式中 m采高，mm； h下沉系数，取决于顶板管理方法。当开采ds一段矿层时，地表产生的下沉盆地为：（2-39）显然，当矿层自s=0 开采至 s=

8、,地表稳定后的下沉盆地表达式可写为： （3-40）令h（x-s）=l, ds=-dl/h，相应的积分限变为hx及-，代入上式得：（2-41）当 x=-时，地表点应有最大值wmax,故（2-42） 故上式 （2-43）将上式代入式（2-41）得：（2-44）在上式中，当wmax 给定时，取不同的h（见式2-48）值，作出w（x）对应曲线见图1.23表明，对于不同的h值所得到的理论下沉盆地具有不同的横向“发育”，下沉曲线也具有不同的最大斜率，其可通过对式（2-44）的微分得出：（2-45）当x=0，盆地具有最大的斜率：（2-46）又，从几何关系知（2-47）式中 r 主要影响半径， 。其中：H采深

9、，tgb主要影响角正切（见图2.62.7）。将式（2-47）代入式（2-44）得：（2-49）在这里引入误差积分函数：则上式可写成如下形式：（2-50）（2）倾斜、曲率表达式 在式（2-35）对x取一阶、二阶导数，便可得到下沉盆地倾斜分布及曲率分布表达式：（2-51）（2-52）（3）水平移动、水平变形表达式 由式（2-47）对z求导，得h， 并将其代入式（2-37）中，得单元水平移动表达式：（2-53）将式（2-47）代入式（2-35）中并对 x求导得单元盆地的倾斜表达式：（2-54）比较式（2-53）与（2-54）得：式中，是不依赖于x而依赖于z的参量，令其等于B（z），则有：（2-55）

10、对于地表来说，z等于开采深度H，B（z）为常数，并可令它等于B，则根据上式可得半无限开采条件下的水平分布表达式：式中 。由于当 x=0 时，上式具有最大值 令 b被称之为水平移动系数，则：（2-56）在上式中对 x微分得水平变形表达式：（2-57）2. 极值及其位置在式（2-51）、（2-56）不难看出，当x=0时，地表倾斜和水平变形达到最大值：（2-58）（2-59）在式（2-52）、（2-57）中对x微分，并令其等于零，同时当 时得：（2-60）（2-61）根据上述分析，将概率积分法的半无限开采下沉盆地主断面上任意点的移动变形预计公式列于表2-1，其移动变形分布形态如图2.3所示。为了便于

11、计算，在表2-2中以为引数给出了无因次量、的数值。用移动变形的最大值，乘以表2-2给出的移动变形的相对值（分布系数），就得到下沉盆地主断面内各点的移动和变形值。表2-1 概率积分法的半无限开采下沉盆地主断面内任意点移动变形预计公式指 标主断面内任意点计算公式极值计算公式极值出现的位置下 沉倾 斜曲 率水平变形水平移动2.2.4有限开采时地表移动盆地主断面的移动变形预计1. 走向主断面的移动变形预计设矿层沿倾斜方向已达到充分采动，沿走向方向没有达到充分采动，在移动盆地走向主断面（图2.4）的计算宽度l。根据叠加原理可求得此条件下的地表下沉分布预计式为：（2-62）令： 得：，在上式中作变元，换元有： （2-63）由此可见，有限开采实际上是由两个半无限开采叠加的结果，即w（x）