第二章 采动地表移动变形预计

重点：①预计理论体系概况；

②概率积分法。基本含义、基本概念、应用条件、应用方法、分布规律、特征值的确定方法，极值公式及计算、按特征值绘制移动变形分布图。

③半无限开采及半无限叠加方法； ④地表任一点移动变形预计方法；

⑤动态移动变形与静态方法的区别及其评价方法。 2.1 地表移动和变形预计理论方法概述

开采沉陷损害预计理论，可以概括为影响函数方法，理论模型方法，经验方法三大类型。 2.1.1 影响函数方法

①国内外学者及理论应用情况；

②假定开采单元矿层dv,其水平投影面积为dp,单元矿层开采引起地表点A的下沉表达式为：dwa?m?f(s)dp

（2-1）

③影响函数的可叠加性；

根据影响函数的叠加原理，对于开采范围为P的矿层开采引起地表点A的下沉量的通式表示为： 2.1.2 经验方法

①前苏联应用的负指数函数方法；②英国煤田方法（NCB.1975）；③波兰学者Z.Kowalczyk (1972)积分网格法；④中国学者何国清提出的威布尔分布法；⑤各矿区通过观测曲线拟合得出的适用本矿区的典型剖面曲线法等。 2.1.3 理论模型方法

属于理论模型方法是建立在力学模型上的，以及建立在弹性或塑性理论基础上的计算方法。在这方面主要有以A.Salstowicz（1958）等为代表的固体力学理论；J.Litwiniszy（1963）等为代表的随机介质理论。建立在弹性或塑性理论基础上的计算方法如：有限单元法（FEM）；边界元法（BEM）；离散元(DEM)等方法；非线性力学（Nonlinear）等方法。

目前应用情况简介 2.2 概率积分法（重点）

目前已成为我国乃至世界范围较为成熟、应用最广泛的预计方法之一。 2.2.1 水平成层介质中的单元盆地

开采沉陷的随机性随机介质理论为基础 ①非连续介质单元模型，②单元相互分离并发生相对运动。

如图2.1在三维问题中，地下（x0, y0, z0）处开采使地表点A(x, y, z)附近某一小块面积ds发生下沉这一事件的概率为：

wa?m???f(s)dp

P （2-9）

P(ds)??(x,y,z)ds

（2-10）

（x, y, z）为密度函数。

在x-z剖面的z水平上x处的一段岩石条dx内有下沉发生，同时在y-z剖面的同一高度上y处的一段岩石条dy内有下沉发生。因此，ds 发生下沉这一事件的概率为发生上述两事件的概率之积，即：

P(ds)?f(x2)dxf(y2)dy?f(x2)f(y2)ds

（2-11）

f为密度函数，x2、y2?对称性。

过原点引另一组正交水平轴x、y?, 使A点在这一系统中的坐标为(x坐标系中，A点附近的微面ds发生下沉的概率为：

, y?, z? )。在新

P(ds?)?f(x?2)f(y?2)ds?

（2-12）

岩石水平成层f的形态皆一致。微面面积不变ds=ds? ；开采点与被考虑的微面相对

位置不变P(ds)=P(ds? )。从本质上讲，某一既定的微面在同一开采影响下的下沉概率与坐标轴方向的选择无关。

P(ds)?P(ds?)?f(x2)f(y2)ds?f(x?2)f(y?2)ds?

（2-13）

从而可得：

f(x2)f(y2)?f(x?2)f(y?2)

（2-14）

当采用下列方式选择坐标轴，使ox? 经过A点，且：

x?2?x2?y2??? 2?y??0? （2-15）

代入（2-14）得：

f(x2)f(y2)?f(x2?y2)f(0)?Cf(x2?y2)

（2-16）

将（2-16）两端对x2、y2取偏微分可得：

df(x2)?f(x2?y2)?(x2?y2)?f(x2?y2)f(y)?C??C 2222d(x)?(x?y)?(x)?(x2?y2)2 （2-17）

df(y2)?f(x2?y2)?(x2?y2)?f(x2?y2)f(x)?C??C

d(y2)?(x2?y2)?(y2)?(x2?y2)2 （2-18）

由此可得：

2df(x2)2df(y) f(y)?f(x)22d(x)d(y)2 （2-19）

移项得到：

1df(x2)1df(y2) ???2222f(x)d(x)f(y)d(y) （2-20）

上式左边为x的函数，右边为y的函数。此方程成立的条件是左右两端都不依赖于自变

量x、y，故可令式（2-20）为常量k，从而有：

df(x2)2?kf(x) 2)d(x （2-21）

将x2看作自变量，解此常微分方程得到：

f(x2)?Pekx

2 （2-22）

式中， P为积分常数。显然，远离中心两端的岩石下沉的概率小。因此，从物理意义说k

为负值，令k=-h2代入上式得：

f(x2)?P?exp(?h2x2)

（2-23）

同理

f(y2)?P?exp(?h2y2)

（2-24）

将上述结果代入式（2-13）有：

P(ds)??(x,y,z)?P2?exp?h2(x2?y2ds

?? （2-25）

P(ds) z0煤层的开采以随机的方式传至z水平上的随机分布。由于

w(x,y,z)?P(ds)w(x0,y0,z0)，由此可得到顶板下沉盆地中的分块下沉体在z水平上造成

的微小下沉盆地的表达式：

w(x,y,z)?P2?exp?h2x2?y2w?x0,y0,z0)?ds0

???? （2-26）

下沉盆地体积=（下沉概率×采高×开采面积） 单元开采分块下沉体，单元盆微下沉盆地。 令开采面积s=1个微小单位，采高w(x0, y0, z0)=1个微小单位的开采为单元开采。单

元开采在上覆岩层中造成的下沉盆地为单元下沉盆地we,于是得单元盆地下沉：

we(x,y,z)?P2exp?h2x2?y2

???? （2-27）

当然，孤立的下沉盆地并不存在，它们存在于统计意义之中。

①采出体积足够大地表出现下沉； ②大面积开采=∑单元开采；

③原始的位置→弹性变形→逐渐垮落→传递到地表（“三带”形成过程）→下沉盆地。 单元盆地是时间函数we=(x, y, z, t)，对应的地表下沉盆地体积为：

Ve???????we(x,y,z,t)dxdy

???? （2-28）

Ve也是时间的函数。根据A.Salstowicz假设，下沉盆地体积的增长率与采空区域未压密

的体积成正比，即：

dVe?c??1?Ve? dt （2-29）

式中 c—比例系数，当t=0时，Ve=0；当t??时，Ve=1 。 解上述方程式可得：

Ve?1?e?ct

（2-30）

盆地体积Ve可由上述求得：

????2222Ve???dxdy?P2????Pexp?hx?y?????h2 （2-31）

????22Ve???Pexp?hxdxp?h2y2dy?P????Pex?????????h?P?h ?2（①换元hx=，dx=（1/h）d；hy=，dy=（1/h）d；②?e??d???） ??则： 从而：

P?2h2??1?e?

?ct （与时间有关的概率） （2-32）

we(x,y,z,t)?h2??1?e?exp??h?x?ct22?y2

?? （2-33）

t??,

we(x,y,z)?h2?exp??h2?x2?y2??

（2-34）

当宽与高均为一个微小单位，长度为无穷大时的开采称为二维单元开采。当开采长轴方向平行于y轴，此时形成槽形盆地，槽长平行于y轴，而在平行于x轴的任何剖面上，此盆地的形状均相同，称为二维单元盆地，表示为：

we(x,z,t)??1?e?ct????h2???exp?h2x2??y???d?

2????从而可得水平成层介质中二维单元下沉盆地的最终表达式为：y=0；（- 2）=2；1/hd=d