煤层顶底板破坏深计算

1、煤层底板破坏深度计算目前，国内外对底板破坏深度的研究已经有许多种方法，本次研究主要是运用弹塑性力学方法结合莫尔库仑（Mohr-Coulomb）强度理论，依现场观测数据为依据，辅助进行计算机数值模拟，综合计算显德汪矿9#煤层底板岩体受采动影响的最大破坏深度，并提出该矿区9#煤层底板破坏深度的经验公式，为企业的安全开采提供科学依据。6.1底板岩体破坏带空间分布形态许多学者对煤层底板采动影响规律进行了研究，提出了煤层底板岩体采动带的空间分布形态。6.1.1近水平煤层煤层回采后，其顶板以冒落角向上冒落，最终形成顶板中部冒落的比较充分，采空区在中部充填较密实，而在采空区两侧顶板冒落得最不充分，充填不实（

2、图6-1）。煤层底板在采空区两侧有较大的自由空间，在地应力作用下，底板岩体能够充分膨胀，产生较多的采动裂隙，近水平煤层在采动边缘下方附近岩体的破坏深度最大。图6-1 煤层顶板岩体冒落示意图煤层底板中破坏带的形态也可用计算的方法得出。考虑到底板岩体的受力状态，以图6-2中的X1XI剖面作为计算模型，作用在弹性表面某一局部面积上的力系，被作用在同一局部面积上的另一静力等效力系所代替，则载荷的这种重新分布，只在离载荷作用很近的地方才使应力的分布发生显著变化，在离载荷较远处影响极小。图6-2中的X1XI剖面的应力分布图形可采用等效模型（图6-3）代替。图6-2中等效应力q(n+1)P0/2，作用宽度为

3、工作面端部至应力峰值距离（xa）的 图6-2 长壁工作面支承压力分布图图6-3 底板上应力简化示意图（P0原始应力）2倍，即L2 xa。煤层底板内岩体自重产生的应力为z，在平面应变状态中，底板岩体任意点M的主应力为： （6-1） （6-2） （6-3）在多向应力作用下，岩体发生破坏时服从Mohr-Coulomb破坏准则，即 1-K3Rc，将（6-16-3）式代入上式后，得： （6-4） 式中：Rc岩体的单轴抗压强度；岩体的容重。 （6-5）式中：底板岩体内摩擦角。 （6-6） 根据公式（6-4）计算得出的煤层底板岩体中破坏带的形态见图6-4所示，破坏带的范围按图中1 2 3 4的顺序发展。 图

4、6-4 岩体破坏带的发展过程（h1底板岩体最大破坏深度）6.1.2倾斜及急倾斜煤层（1971）采用石膏-硅藻土模拟了层状块状岩体在不同条件下受均布力作用的状况，得出了岩层中应力的传播规律（图6-5）。经分析知：在相同的外载作用下，当0°45°时，岩层在上山之间的应力传播深度大于下山方向；当45°90°时，则相反；当0°或90°时，岩层中的应力传播深度在各个方向相同。对于在煤层回采工作面，则煤壁正下方的底板岩体受集中应力的作用，其应力传播与图6-5有类似的规律。当煤层倾角0°45°时，煤层底板在下山方向的应力传播深度

5、大于上山方向6，从而导致煤层底板破坏带的深度在下山方向较大（图6-6a）。当煤层倾角45°90°时，煤层底板在上山方向的应力传播深度大于下山方向，在上山方向煤层底板破坏带的深度较大（图6-6b）。 图6-5 层状块状岩体底板内部应力等值线（为岩层倾角） a b 图6-6 煤层底板破坏形态 （a为倾斜煤层；b为急倾斜煤层）从受力角度分析，在缓倾斜及倾斜煤层的下出口附近，煤层底板承受的集中应力大于上出口附近的集中应力，煤层底板承受的水压力也是下出口附近的较大，煤层底板在下出口附近裂隙较发育。由于煤层的倾角不大，采空区冒落岩石的滑移起不了主要作用，然而急倾斜煤层不同，在自重力作用

6、下，冒落的采空区的岩石将向采空区下部滑移，充填坚实的冒落岩石限制了下出口附近煤层底板的膨胀，阻碍了裂隙的形成，所以，急倾斜煤层的底板破坏深度在上出口附近较大。显德汪矿9#煤层属于缓倾斜煤层，所以在下出口附近底板岩体的破坏深度值大于上出口。6.2底板岩体破坏深度理论计算6.2.1采场边缘破坏深度计算根据张金才等人的研究10，采场边缘的应力场为： （6-7） （6-8） （6-9） 根据弹性理论知，求解主应力的公式为： （6-10）把采场边缘的应力场计算公式（6-76-9）代入上式，可求得在平面状态下（即3=0）的采场边缘的主应力：（6-11） （6-12） （6-13）将（6-116-13）式代

7、入Mohr-Coulomb破坏准则方程（1-K3=Rc）后，可得到平面应力的采场边缘破坏区的边界方程： （6-14）根据（6-14）方程，可绘制出采场边缘由于应力集中而形成的破坏区形态（图6-7）。根据该图形可知，垂直于开采层的岩体破坏深度h为： （6-15） 即： （6-16） 图6-7 采场边缘岩体破坏形态为求解煤层底板岩体最大破坏深度，令dh/d0，得： （6-17） 求解以上三次方程，得到有效解约为74.84°，也就是当为74.84°时，采场边缘底板岩体的破坏深度为最大值（hm）： （6-18）岩体最大破坏深度距工作面端部的距离Lm为： （6-19）根据以上分析可知

8、，煤层底板岩体最大破坏深度与工作面倾斜长度成线性关系，与岩体原始应力的平方成正比，与岩体抗压强度成反比。6.2.2塑性理论计算岩体最大破坏深度如果煤层底板主要由软岩构成，煤体边缘一定范围内的底板岩体，当作用在其上的支承压力达到或超过其临界值时，岩体中将产生塑性变形，形成塑性区；当支承压力达到导致部分岩体完全破坏的最大载荷时，支承压力作用区域周围的岩体塑性区将连成一片，致使采空区内底板隆起，已发生塑性变形的岩体向采空区内移动，并形成一个连续的滑移面。6.2.2.1煤层屈服区长度的计算煤层屈服区长度Xa，可以通过现场实际测量获得，也可以通过计算得到。（1） 利用煤层内聚力（Cm）的计算公式： （6

9、-20）式中：煤层内摩擦角；Cm煤层内聚力；m煤层采高。显德汪矿1192工作面9#煤层参数为：煤层的内摩擦角45°；煤层的内聚力Cm2.8MPa；煤层的采高m 3.50m；平均埋深H=385m；26 KN/m3；n为最大应力集中系数，利用经验公式n1+0.23Lx0.47，得n3.266。 （6-21）屈服区长度xa： （6-22） （2） 利用国外学者得出的煤层屈服长度的经验公式： （6-23）（3）A.H.Wilson提出的煤层屈服区长度计算公式： （6-24）显德汪矿1192工作面煤层参数为：m3.50m，K15.828，代入下式 （6-25）计算得F6.708，代入公式（62

10、5），得 xa6.00m研究过程中，由于实验室9#煤层内聚力数值无直接测试数据，计算所采用数据为同一地区周围矿井的相似材料参数数据，所得到的屈服区长度结果误差较大。经验公式的计算，考虑的影响因素较少，计算结果存在误差相对也较大，经过对三种方案的对比，并考虑实际情况，选取第三种计算方案作为煤层屈服长度，即xa6.00m。6.2.2.2底板最大破坏深度的塑性解魏西克（A.S.Vesic）通过压膜试验及现场实际经验，提出了岩土产生塑性滑移时的极限承载力的计算公式。结合煤层支撑压力作用的特点，张金才对极限承载力计算公式进行了修改与补充，即底板岩体的承载极限： （6-26）式中：xa煤体屈服区的长度；C

11、底板岩体的内聚力。煤层底板岩体的滑移线（即塑性区的边界）由三个区域组成：主动极限区aab、过渡区abc及被动极限区acd（图6-8），主动区和被动区的滑移线各由两条直线组成；过渡区的滑移线一组由对数螺线组成，另一组为自a为起点的放射线。其对数双螺线（图6-9）方程为： 图6-8 煤层底板中塑性破坏区 图6-9 对数螺线示意图 （为主动区，过渡区，被动区）煤矿生产过程中煤层底板岩体发生底鼓的现象可以用塑性区的形成和发展过程加以解释。煤层开采后，在采空区周围的底板岩体上产生支承压力，当支承压力作用区域的岩体（图6-8中区）所承受的应力超过其极限强度时，岩体将会产生塑性变形，由于这部分岩体在垂直方向

12、上受到压缩，在水平方向上必然会膨胀，膨胀的岩体挤压过渡区（图6-8中区）的岩体，并且将应力传递到这一区域。过渡区的岩体受到挤压后将继续挤压被动区的岩体（图6-8中区）。由于过渡区和被动区有临空面（采空区），在主动区传递来的力的作用下，过渡区和被动区的岩体将向采空区内膨胀（即底板的压延作用）。图6-10 煤层底板最大破坏深度计算简图根据张金才对魏西克（A.S.Vesic）公式进行修正的岩土层极限承载力的综合计算公式，可得到底板岩体的极限载荷，从而得出极限支承压力条件下破坏区的最大深度和长度计算公式。煤层底板岩体最大破坏深度h1确定：根据底板最大破坏深度计算简图，在aba中，在aef中，而 （6-

13、27）因此： （6-28）由 ，可以求出破坏区的最大破坏深度h1。所以，得到 （6-29）将（6-29）式和r0代入（6-28）式，即可得到底板的最大破坏深度h1。 （6-30）式中底板岩体的内摩擦角（底板岩体以泥岩、细砂岩为主，综合考虑取39°），将以上参数代入6-30式，求得： （1）煤层底板岩体最大破坏深度距工作面端部的水平距离L1： （6-31） （2）采空区内底板破坏区沿水平方向的最大长度L2：（6-32） 上述计算结果，可绘制出按塑性理论求得的显德汪矿9#煤层底板破坏带的分布形态（图6-11）。 6.00 10.92 60.32 13.48 图6-11 9#煤底板破坏带形

14、态（理论计算结果图，图中单位m）6.3基于MATLAB的底板岩体破坏经验公式6.3.1 MATLAB简介在科学与工程计算中，为了把握某些规律，经常需要研究和探寻一些变量之间的关系。而变量之间的关系有时是确定的，有时又是非确定的。研究这些确定性与非确定性之间关系的一个重要工具就是回归分析方法。回归分析是数理统计学中重要的一部分，在科学实验、数学模型的建立、决策预测及自动控制中都有着广泛的应用。MATLAB作为一种强大的科学计算工具，已受到各专业人员的广泛重视。MATLAB最突出的特点就是简洁。它用更直观的，符合人们思维习惯的代码代替C和Fortran语言的冗长代码，使更多的工程技术人员摆脱了烦琐

15、的编程工作，能够将注意力集中在专业技术研究的核心问题上。 MATLAB（Matrix Laboratory）是美国Mathworks公司自80年代中期推出的数学软件。经过近20年的发展，它己经成为科学计算、视图交互系统、动态系统仿真等的基本工具。现在，MATLAB可运行在Windows系列，OS2，Unix，Linux等十几个平台上。MATLAB的优势主要体现在以下几个特点上：（1）语言简洁紧凑，使用方便灵活，库函数极其丰富，压缩了一切不必要的编程工作；（2）运算符丰富。MATLAB是由C语言编写的，它提供了和C语言几乎一样的运算符；（3）语法限制不严格，程序设计自由度大；（4）程序可移植性很

16、好，基本上不用做修改就可在各种型号计算机和操作系统上运行；（5）图形功能强大，数据的可视化非常简单；（6）功能强劲的工具箱是MATLAB的另一重大特色，常用的工具箱有：分析与综合工具箱（Analysis and Synthesis Toolbox），神经网络工具箱（Neural Network Toolbox），最优化工具箱（Optimization Toolbox）等；（7）原程序的开放性，除内部函数以外，所有的MATLAB的核心文件和工具箱文件都是可读可改的源文件，用户可通过对源文件的修改以及加入自己的文件构成新的工具箱。6.3.2 MATLAB的程序实现回归分析是数理统计中最常用的方法之

17、一，一般用最小二乘法确定回归方程中的系数，其矩阵计算过程颇为复杂。而用MATLAB实现则使问题大大简化，MATLAB中有四个函数可以用于回归分析和拟合；polyfit（x，y，n）、leastsq（/function/，x）、regress（y，x）和 pinv（A）*ypolyfit（x，y，n）只能用于多项式线性回归，leastsq（/function/，x）可用来做非线性回归，regress（y，x）可用于多元线性回归，pinv（A）*y可用于求解线性方程组。6.3.3 经验公式确定通过对1192工作面的现场观测分析、数值模拟、理论计算结果分析，获得了工作面回采过程中的最大破坏深度，结合

18、工作面突水实例及室内分析计算表明，显德汪矿深部水平底板破坏深度与工作面埋深、倾斜长度、采高之间有一定的相关关系，但将1192工作面各因素的相关数据带入显德汪矿浅部经验公式中，发现计算结果与实际误差很大，说明1291工作面研究中得出的经验公式已不适合深部水平开采底板破坏深度的预测，因此本次研究基于MATLAB数学软件，结合本矿及邻区典型工作面相关数据，利用回归分析方法重新得出埋深在650m之内的煤层底板破坏深度与上述诸因素的关系式（经验公式）： （6-31） 式中：h煤层底板破坏深度（m）；H煤层埋藏深度（m）；L工作面倾斜长度（m）；M工作面回采高度（m）。 6.3.4各因素对经验公式的作用将

19、有关数值代入经验公式，使有关参数从下限变到上限，可求出增量h值，以分析各参数的作用。根据显德汪矿的实际情况，各参数的上下限取值如下：埋 深：Hmin206m Hmax636m倾斜长度：Lmin90m Lmax135m 回采高度：Mmin3m Mmax5m则： 计算结果表明，工作面倾斜长度对h的影响最大，其次是煤层埋深。表6-2 实测底板破坏深度与经验公式计算结果对比表地 点对底板破坏深度实测结果计算破坏深度（m）误 差埋 深（m）斜 长（m）采 高（m）破 坏深 度（m）相对误差（）绝对误差（m）邯郸王凤矿1930面118802.5109.009.9-0.99峰峰三矿3707面1301351.51213.815.31.84邢台矿7607窄面320605.49.78.89.3-0.91邯郸王凤矿1830面123704