Mine Engineering 矿山工程, 2020, 8(3), 361-367
Published Online July 2020 in Hans. http://www.hanspub.org/journal/me https://doi.org/10.12677/me.2020.83046
Study on Deformation and Failure Characteristics of Roof and Floor Mining Based on FLAC3D Model Ubiquitous
Yuan Zhang, Haifeng Lu, Duoxi Yao
School of Earth and Environment, Anhui University of Science and Technology, Huainan Anhui
Received: Jun. 28th, 2020; accepted: Jul. 13th, 2020; published: Jul. 20th, 2020
Abstract
In the deep coal seam mining, the mining face of roof and floor strata sedimentary structure affect the strength of the roof and floor, roof and floor when bedding plane development, are more likely to happen damage deformation, causing working face roof and floor of floor heave of caving accidents, serious threat to the mine staff safety of life and destruction of mining equipment of working face. In this paper, Model ubiquitous and Model mohr of FLAC3D numer-ical simulation are used to compare and analyze the failure and deformation of the working face floor during the bedding development, and the following results are obtained: The struc-tural strength of the roof and floor of the coal mining face is lower than that of the area where the bedding structure is developed. The top and bottom plates are more prone to failure and deformation. When exploiting coal resources in the area with bedding structure development, attention should be paid to strengthening the prevention of water damage to roof and floor and the support of goaf.
Keywords
Coal Seam Floor, FLAC3D, Model Ubiquitous and Model Mohr, Failure Deformation Characteristics
基于FLAC3D多节理模型的顶底板采动变形破坏特征研究
张 元，鲁海峰，姚多喜
安徽理工大学地球与环境学院，安徽 淮南
文章引用: 张元, 鲁海峰, 姚多喜. 基于FLAC3D多节理模型的顶底板采动变形破坏特征研究[J]. 矿山工程, 2020, 8(3): 361-367. DOI: 10.12677/me.2020.83046
张元 等
摘 要
收稿日期：2020年6月28日；录用日期：2020年7月13日；发布日期：2020年7月20日
在对深部煤层进行开采时，采煤工作面顶底板的地层沉积结构会对顶底板的强度产生影响，当顶底板层理面发育时，更易发生破坏变形，进而引发工作面顶底板的垮落底鼓事故，严重威胁矿井工作人员的生命安全以及破坏工作面开采设备。本文通过FLAC3D数值模拟的多节理和库仑摩尔本构模型对比分析工作面底板在层理发育时的破坏变形情况，得出：层理构造发育区域，采煤工作面顶底板的结构强度较层理不发育地区要低。导致顶底板更容易发生破坏变形。在层理构造发育地区开采煤资源时，应注意加强对顶底板的水害防治与采空区的支护。
关键词
煤层底板，FLAC3D，多节理模型，破坏变形特征
Copyright ? 2020 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/
1. 引言
Open Access 深部煤资源开采时，由于开采扰动会对采煤工作面顶底板产生扰动破坏，由于顶底板破坏变形，可能会引起工作面的垮塌事故，严重威胁矿井工作人员的生命安全以及破坏工作面开采设备。目前大量学者针对采煤工作面顶底板破坏变形特征进行了深入研究[1]-[13]，多数学者利用数值模拟方法对采煤工作面顶底板破坏变形特征从多角度进行研究[1]-[10]，其中部分学者基于多节理本构模型模拟分析了工作面顶底板层状地层的破坏变形特征[1] [2] [3]。部分学者从实际开采工作面实例出发，结合数值模拟分析了工作面布设规格，顶底板岩性对顶底板破坏变形特征的影响[4]-[10]。另有部分学者基于相似材料模拟对采煤工作面顶底板支撑应力分布特征以及不同构造环境中顶底板破坏变形特征进行分析[11] [12] [13]。然而对于层理构造对采煤工作面顶底板的破坏变形研究比较少，对于水平层理构造发育地区，多数情况下人们往往会忽视其对工作面顶底板的破坏变形的影响。为探究地层沉积层理构造对采煤工作面顶底板破坏变形的影响，本文基于FLAC3D数值模拟方法，采用多节理模型模拟地层层理构造，分析采煤工作面顶底板的破坏变形特征，并与摩尔库仑模型进行对比分析，探究层状地层对采煤工作面顶底板的破坏变形的影响规律。
2. 数值模型
2.1. 模型参数
为模拟层理发育时的工作面顶底板围岩的层状构造，本文采用多节理本构模型来模拟沉积层理，工作面顶板为细砂岩，底板为砂质泥岩。模型参数如表1所示。
DOI: 10.12677/me.2020.83046
362
矿山工程
张元 等
Table 1. Multi-joint model parameters 表1. 多节理模型参数
密度 地层 抗拉强度MPa 摩擦角(?) 内聚力MPa 剪切模量GPa 体积模量GPa 节理抗拉强度MPa 节理摩擦角(?) 节理内聚力MPa
(kg/m3) 顶板 2487 煤层 1380 底板 2483
2.1 1.2 2.2
38 32 32
5.8 3.5 4.8
6.01 2.01 7.35
5.97 4.91 9.97
1.4 0.6 1.3
18 12 15
2.2 1.5 2.2
同时为对比分析两种本构模型下的工作面底板破坏变形情况，库仑摩尔模型的模型除节理参数，其余参数均与多节理模型一致。
2.2. 数值模型及工作面布置
本次数值工作面采动数值模型的尺寸为300 * 300 * 200 m (图1)，由100 m厚顶板，97 m厚底板及3 m厚的煤层组成(图1中黑色部分)，模型边界为位移与应力组合的混合边界，底部为全位移约束，四周为侧向位移约束。工作面处于?400 m水平，等效在模型顶部施加10 MPa的垂向荷载，四周边界施加2.5 MPa的水平应力。
Figure 1. Numerical model and meshing
图1. 数值模型及网格划分
工作面开采方向平行煤层，沿x轴方向开采，工作面跨度150 m，模拟时取开采步长20 m，分5步开采，总推进长度为100 m。
3. 模拟结果分析
3.1. 破坏特征分析
为探究工作面底板破坏特征，对比分析在两种本构模型下的每一步开采对工作面底板破坏的影响，如图2所示。
由图2可知，在多节理模型中，破坏类型主要为张拉破坏和剪切破坏，张拉破坏区域主要分布在顶底板中部区域，剪切破坏主要分布在工作面顶底板的两侧。并且随着开采进度顶底板塑性破坏区分布范围逐渐扩大，工作面顶底板破坏深度逐渐向顶底板围岩延伸，最后一步开挖达到最大破坏深度30 m，破坏区垂向分布范围40~60 m。开挖过程中，先采区域顶底板破坏垂向分布范围较后采区域要广。由此可知在层理发育的地层环境中，开采工作面顶底板容易发生破坏，且先采区域破坏程度要大于后采区域。为对比分析，导出库仑摩尔模型逐步开挖底板的塑性区分布。
DOI: 10.12677/me.2020.83046
363
矿山工程
张元 等
Figure 2. Distribution of plastic zone in XZ section of multi-joint model 图2. 多节理模型XZ切面塑性区分布
Figure 3. Distribution of plastic zone in the XZ section of Coulomb's model 图3. 库仑摩尔模型XZ切面塑性区分布
由图3分析可知，在库仑摩尔模型中，随着工作面开采进尺，顶底板的塑性区分布纵向延伸范围较多节理模型的破坏范围要小，垂向破坏分布范围在20 m左右，并且破坏类型以顶底板的张拉破坏为主，仅在第三步开挖时在底板右侧出现剪切破坏。并且先采区与后采区的纵深破坏范围分布受开挖进尺影响较小。
对比分析两种模型下工作面底板的破坏特征分析可知，多节理模型中工作面顶底板更容易发生破坏，随着开挖底板破坏深度逐渐增大，且先采区域的底板破坏程度较后采区域破坏大。同时，库仑摩尔模型中底板破坏程度较多节理模型要小，且先采区与后采区的破坏程度差别不大。综上，可知在沉积层理构
DOI: 10.12677/me.2020.83046
364
矿山工程
张元 等
造发育区域，开挖扰动对采煤工作面顶底板破坏区域的分布范围影响较大，为预防开采次生裂隙导通含水层，应注意加强顶底板突水水害的防治。
3.2. 位移特征分析
为分析工作面顶底板的位移情况，导出采区中部XZ横截面的位移云图。对比分析在两种本构模型下的每一步开采对工作面顶底板变形的影响，如图4所示。
Figure 4. Vertical displacement of XZ section of multi-joint model 图4. 多节理模型XZ切面垂直位移
由图4可知，工作面顶底板最大位移量分布于开挖总步长中部区域，且随着开采逐步进行顶底板最大垂直位移量逐步增大，为进一步分析每步开挖顶底板的位移特征，导出每步开挖后顶底板的最大垂向位移，如表2所示。
Table 2. Multi-joint model gradually excavates the maximum vertical displacement of roof and floor 表2. 多节理模型逐步开挖顶底板最大垂直位移
开挖步 顶板位移 底板位移 相对位移变量
20 m ?0.4347 0.2965 0.7312
40 m ?3.957 2.8254 6.7824
60 m ?5.6557 3.6914 9.3471
80 m ?7.4118 4.3487 11.7605
100 m ?9.0443 4.8527 13.897
同理导出库仑摩尔模型逐步开挖垂直位移云图和最大位移量，如图5和表3。
Table 3. The maximum vertical displacement of the roof and floor of the Coulomb Moore Model 表3. 库仑摩尔模型逐步开挖顶底板最大垂直位移
开挖步距 顶板位移/cm 底板位移/cm 相对位移变量/cm
20 m ?0.4347 0.2965 0.7312
40 m ?3.6561 2.5692 6.2253
60 m ?5.3426 3.4508 8.7934
80 m ?7.2147 4.2655 11.4802
100 m ?8.6664 4.7566 13.423
DOI: 10.12677/me.2020.83046
365
矿山工程
基于FLAC3D多节理模型的顶底板采动变形破坏特征研究 - 图文
