D0I:10.13374/j.issn1001-053x.2004.06.028 第26卷第6期 北京科技大学学报 Vol.26 No.6 2004年12月 Journal of University of Science and Technology Beijing Dec.2004 灵新煤矿西天河下安全开采技术综合分析 蔡美峰》任奋华”来兴平12) 1)北京科技大学士木与环境工程学院,北京1000832)西安科技大学能源学院,西安710054 摘要为了研究灵新煤矿河下开采引起的覆岩破坏规律和地表移动规律,采用工程地质 调查、彩色钻孔孔电视、钻孔冲洗液漏失量法、物理与数值模拟、地应力测量及地表监测等多 种研究手段,获得了矿区地应力分布规律、工作面开采的垮落带和导水裂缝带的高度、岩层 移动规律及矿压显现规律,实现了L3414综采工作面安全开采和环境保护的双重目标 关键词地应力测量:覆岩;破坏高度;开采方式;控水开采 分类号TD325 西天河是宁夏回族自治区灵武矿区内惟一 1地应力测量 的常年地表径流,它阻隔毛乌素沙漠的扩展,对 区域生态环境有重要的调节作用.它流经灵武矿 采用实现完全温度补偿并考虑岩体非线性 区8387m,给矿井生产带来巨大的安全隐患,造 的地应力解除测量技术,首次进行了灵新矿区地 成大量煤炭储量的滞压.本文采用了现场调查、 应力现场实测,获得了矿区地应力分布规律,建 室内实验、理论分析、物理与数值模拟、现场试验 立了矿区地应力场模型.为了测定整个矿区区域 等多种研究手段,对开采引起的覆岩破坏规律、 内地应力的数值、方向和空间分布规律,在不同 地表移动规律进行了多方面的综合性研究:在维 位置、不同深度布置了3个测点,进行地应力测 持西天河现状和灵新煤矿现有的采煤方法和开 量的应力解除实验.测点的具体位置如下:1号孔 拓方式的前提下,提出并优化了L3414综采工作 位于一采区1160m水平的轨道石门,方位196°: 面西天河下控水采煤方案:通过L3414综采工作 2号孔位于三采区1186m水平的轨道石门,方位 面的开采试验,进行现场矿压、岩移和水文监测, 194°:3号孔位于三采区1186m水平的风巷绕道, 取得了较全面、系统的技术成果和开采经验. 钻孔方位14°.其测量结果见表1. 表】地应力测量结果 Table 1 Results of in-situ stress measurement 最大主应力 中间主应力 最小主应力 测点 数值MPa方向/()倾角)数值MPa方向)倾角/()数值MPa方向/)倾角/() 4.05 287.5 -2.9 1.97 183.8 -77.9 1.34 198.1 11.7 2 3.69 300.1 -11.2 1.91 272.5 77.6 1.60 29.0 5.7 3 3.36 94.1 6.2 2.18 352.8 11.8 1.82 201.2 76.6 2覆岩破坏高度预测 量法的观测结果,首次准确地确定了灵新煤矿十 四煤的垮落带和导水裂缝带的高度.十四煤综采 在L3314综采工作面采空区上方地表布置了 的垮落带高度为15.3120.65m(对应采厚为2.7 两个覆岩破坏观测钻孔,采用钻孔冲洗液法、彩 2.2m),冒采比为5.79.4:导水裂缝带的高度为 色钻孔电视系统等进行了全面观测",通过全面 41.42-56.03m(对应采厚为2.2-2.7m),裂高采厚 对比分析彩色钻孔电视系统与钻孔冲洗液漏失 比为18.9-20.8(见表2). 收稿日期20040508蔡美蜂男,60岁,教授,博士 *教育部博士点基金资助项目No.20020008021)
第 2 6 卷 第 6 期 2 0 0 4 年 1 2 月 北 京 科 技 大 学 学 报 J o u r n a l o f U n vi e招 iyt o f S e ci n e e a n d Tc c h n o ot yg B e ij i n g V心 1 . 2 6 刊 0 . 6 D ec 。 2 0 4 灵新煤矿西天河下安全开采技术综合分析 蔡 美峰 ” 任奋 华 ` , 来兴 平 ` ,2) l) 北京科 技大 学土 木与环 境工 程学 院 , 北京 10 0 0 8 3 2) 西 安科技 大学 能源 学 院 , 西 安 71 0 0 54 摘 要 为 了研 究灵 新煤 矿河 下 开采 引起 的覆岩 破 坏规 律和 地表 移动 规律 , 采用 工程 地质 调查 、 彩色钻 孔孔 电视 、 钻 孔冲洗 液漏 失 量法 、 物理 与数值 模 拟 、 地 应力 测量及 地表 监 测等 多 种研 究手 段 , 获得 了矿 区 地 应 力分布 规律 、 工作面 开采 的垮 落 带和 导水 裂缝 带 的高度 、 岩层 移动规 律 及矿 压显 现规 律 , 实现 了 L 3 4 14 综 采工作 面 安全 开采 和环 境保 护 的双重 目标 . 关键 词 地应 力测 量 ; 覆 岩 ; 破坏 高度 ; 开采 方式 ; 控 水开采 分类 号 T D 3 2 5 西天 河 是 宁夏 回族 自治 区灵 武 矿 区 内惟 一 的常年 地表 径 流 , 它 阻隔 毛 乌素 沙 漠 的扩 展 , 对 区域 生态 环境 有 重要 的调 节 作用 . 它流经 灵 武矿 区 8 3 8 7 m , 给 矿 井 生产 带 来 巨 大 的安全 隐患 , 造 成 大量 煤 炭储 量 的滞 压 . 本 文采 用 了现 场 调 查 、 室 内实验 、 理论 分析 、 物 理 与数值 模拟 、 现场 试验 等 多种 研 究手 段 , 对 开采 引起 的覆 岩破 坏 规 律 、 地表 移动 规 律进 行 了多 方面 的综 合性 研 究 ; 在维 持 西 天河 现 状 和 灵 新煤 矿 现 有 的采 煤 方法 和 开 拓 方式 的 前提 下 , 提 出 并优 化 了 L 3 41 4 综 采 工作 面西 天 河 下控 水 采煤方 案 ; 通 过 3L 4 14 综 采 工作 面 的开采 试验 , 进 行现 场 矿压 、 岩 移和 水 文监 测 , 取 得 了较 全面 、 系 统 的技 术成 果 和 开采 经验 . 1 地 应 力测 量 采 用 实现 完 全 温 度 补 偿 并 考 虑 岩 体 非 线 性 的地 应 力解 除测 量技 术 , 首 次进 行 了灵新 矿 区地 应 力现 场 实测 , 获 得 了矿 区地 应 力 分布 规律 , 建 立 了矿 区地 应力 场模 型 . 为 了测 定整 个矿 区 区域 内地 应 力 的数 值 、 方 向和 空 间分 布 规律 , 在不 同 位 置 、 不 同深 度 布置 了 3 个测 点 , 进 行 地 应 力测 量 的应 力解 除实验 . 测点 的具 体位 置 如下 : 1 号孔 位 于一 采 区 1 16 0 m 水平 的轨道 石 门 , 方 位 19 60 ; 2 号 孔位 于 三采 区 1 186 m 水 平 的轨 道石 门 , 方 位 194 “ ; 3 号 孔位 于 三采 区 1 186 m 水平 的风巷 绕道 , 钻 孔方 位 14 “ . 其 测量 结 果 见表 1 . 表 1 地应 力测 量 结果 1恤b l e 1 R e s u lst o f i n 一 s i t u S t er s s m e a s u 代 m e n t 测 点 最 大主应 力 al 中 间主应 力氏 最小 主应 力仍 数值迎讨P a 方 向o(/ ) 倾角 o(/ ) 数 值/州田a 方 向o(/ ) 倾 角 o(/ ) 数值 / M P a 方 向o(/ ) 倾 角(/ o) 4 0 5 3 . 6 9 3 . 3 6 2 8 7 . 5 3 0 0 . 1 9 4 . 1 1 . 9 7 1 . 9 1 2 . 1 8 1 8 3 . 8 2 72 . 5 3 5 2 . 8 1 . 3 4 1 6 0 1 . 8 2 1 1 7 5 . 7 7 6 . 6 2 覆 岩 破坏 高度预 测 在 L 3 14 综采 工作 面 采 空区上 方 地表 布置 了 两个 覆岩 破 坏 观测 钻 孔 , 采 用 钻 孔冲 洗 液法 、 彩 色 钻 孔 电视 系统 等进 行 了全 面 观测` , 周 . 通 过全 面 对 比 分 析 彩 色钻 孔 电视 系 统 与 钻 孔冲 洗 液 漏 失 收稿 日 期 2 0 4 -() 5刁8 蔡美峰 男 , 60 岁 , 教 授 , 博士 * 教 育部 博士 点基金 资助 项 目 (N 。 2 0 0 2 0 0 0 8 02 1) 量 法 的观 测 结果 , 首 次准 确地 确 定 了灵新 煤矿 十 四煤 的垮 落 带和 导水 裂缝 带 的高度 . 十 四煤综 采 的垮 落带 高度 为 15 . 3 1一 20 . 65 m (对 应 采 厚 为 2 . 7 - .2 2 m ) , 冒采 比 为 5 . 7一 .9 4 ; 导 水裂 缝 带 的 高度 为 4 1 . 4 2一 5 6 . 0 3 m (对 应采 厚 为 2 2一 2 . 7 m ) , 裂高 采 厚 比 为 18 . 9一 2 0 . 8 ( 见表 2 ) . DOI: 10. 13374 /j . issn1001 -053x. 2004. 06. 028
Vol.26 No.6 蔡美峰等:灵新煤矿西天河下安全开采技术综合分析 573 表2覆岩破坏高度实测结果 Table 2 Measuring results of the damage height of overburden rock 冒落带 导水裂缝带 钻孔编号 钻孔位置 采厚/m 基岩柱垂尚/m 高度/m 冒采比 高度m裂采比 D 上边界 2.7 15.31 5.7 56.03 20.8 107.56 D, 上边界 2.2 20.65 9.4 41.42 18.8 105.97 3物理模拟与数值模拟分析 厅米力向 3.1相似材料模拟 200m 325m 141m 采用相似材料模拟方法研究L3414综采工作 面开采过程中岩层移动规律,,即在模型的开采 煤层上覆岩层及地表布设不同的观测点,用于观 700m 测研究岩层及地表的移动变形规律,图1为地表 图2数值模拟计算平面模型 下沉规律曲线 Fig.2 Plane model for numerical modelling 现象:在工作面推进过程中,工作面的两端始终 会发生应力集中现象,由于工作面前端受到扰动 的影响较大,因此工作面前端的应力值普遍大于 1作面后端的应力值. 综合分析相似材料模拟和数值计算结果,工 作面开采后岩层下沉呈非连续性,彼此间存在着 9 13 17 31 测点 缝隙,在靠近断裂边界的岩层中存在着永久性的 图1地表下沉曲线 离层缝隙:在初次采动条件下岩层断裂边界以外 Fig.1 Regular curve of surface subsidence 没有出现离层. 相似材料模拟实验开采的结果为:在采厚为 4开采现场观测 2.38m的条件下,工作面上覆岩层的垮落带高度 为10-14m:导水裂缝带高度为32-46m:周期来 4.1矿压观测 压步距为7~14m:裂隙多以斜向甚至纵向为主, 矿压观测的目的是认识和总结L3414工作面 纵横连接,交错贯通. 开采中矿山压力显现的基本规律,包括来压显现 3.2数值模拟分析 特点和周期来压步距和来压强度,为水体下采煤 运用美国明尼苏达大学和美国Itasca Consult- 的参数分析提供技术数据和参考依据.沿工作面 ing Group Inc..开发的二维有限差分计算程序 倾向设了5条观测线,分别距运输顺槽7.5,45,87, FLAC进行计算,该程序主要适用模拟计算地质 135和174m.图3为平行切眼段矿压显现规律图, 材料和岩土工程的力学行为.根据观测结果, 它反映了支架阻力与工作面推进度的相互关系. L3414综采工作面横向变形很小,因此本次计算 矿压观测结果表明,工作面支架阻力、初撑力和 采用平面应变模型假设,即垂直于计算剖面方向 周期来压步距在开采过程中变化不大,周期来压 的变形为零,充分模拟研究了工作面开采过程中 步距为7~14m,其值与相似材料模拟实验所得结 岩层的力学状态、破坏场的分布以及位移场的变 果一致. 化.图2为数值模拟计算的平面模型. 4.2地表移动规律观测研究 由数值模拟分析结果可知:L3414综采工作 在L3414综采工作面地表沿其走向和倾向布 面开采中最大主应力场的分布服从由上至下逐 置了两条观测线,全程观测L34141作面开采过 渐增加的趋势,由于受到工作面开采扰动的影 程主断面上地表移动与变形情况(见图4). 响,预留煤柱和采空区交接区域会发生应力集中 地表移动观测结果表明:在整个开采过程
从 〕1 . 2 6 N o . 6 蔡美 峰等 : 灵新煤矿西 天 河 下 安全开 采技 术综合 分析 . 5 7 3 - 表 2 覆 岩破 坏 高度实 测结 果 Ta b l e 2 M e a s u r i n g re s u l t s o f t h e d a m a g e h e ig h t o f o v e r b u r d e n r o c k 钻孔 编 号 钻 孔位 置 采 厚/ m 冒落 带 导水 裂缝 带 高度/m 冒采 比 高度加 裂采 比 基 岩柱 垂高/ m 卜边 界 卜边 界 巧 . 3 1 之0 石 5 5 . 7 9 . 4 5 6 . 0 3 4 1 . 4 2 2 0 . 18 . 1 0 7 . 5 6 10 5 . 9 7 D 3 物理 模 拟 与数 值模 拟 分 析 3 . 1 相似材 料 模 拟 采 用相 似材料 模拟 方法 研 究L 3 4 14 综采 工 作 面开 采 过程 中岩 层移 动 规律 ` 34J , 即 在模 型 的开 采 煤层 _ [ 覆岩层 及地 表布 设不 同 的观测 点 , 用于 观 测 研究 岩层 及地 表 的移 动变 形规 律 , 图 1 为地 表 下 沉规 律 曲线 . 开 采方 向 L〕 41 4 综采 l 作面 图 2 数 值模 拟计 算平 面模型 P l a n e m o d e l fo r n u m e r i e a l m o d e l li n g 篡巡二晦 4 崔耳ì /日日 9 13 1 7 2 1 测 点 图 l 地表 下 沉 曲线 F i g . l R e g u l a r e u vr e o f s u r fa e e s u b s id e n e e 现 象 ; 在 工 作面 推进 过 程 中 , _ _ L 作 面 的两端 始 终 会 发 生 应 力集 中现象 , 由于 工 作 面前端 受到 扰动 的影 响较大 , 因 此 工 作面 前 端的应 力值 普遍 大于 1 _ 作 面后 端 的应 力值 . 综合 分 析相似 材 料模 拟和 数 值计 算 结果 , 1 _ 作 面开 采 后岩 层 下 沉呈 非连 续性 , 彼此 间存在 着 缝 隙 , 在靠近 断 裂边 界 的岩层 中存 在着 永久 性 的 离层 缝 隙 ; 在 初 次采动 条件 下 岩层 断裂 边界 以 外 没有 出 现离 层 . 相似 材料 模拟 实验 开 采 的结果 为 : 在 采厚 为 .2 38 m 的条件 下 , 工 作面 上 覆岩 层 的垮 落带 高度 为 10一 14 m ; 导水 裂缝 带 高度 为 3 2一 46 m ; 周 期 来 压 步距 为 7一 14 m ; 裂 隙 多 以斜 向甚至 纵 向为主 , 纵 横连 接 , 交 错 贯通 . 3 .2 数 值模 拟 分 析 运用 美 国明尼 苏达 大 学和 美国 It as ca C o ns ul t - in g G or uP nI o . 开 发 的 二 维 有 限 差 分 计 算 程 序 F L A C DZ 进行 计算 , 该 程序 主要适 用模拟 计 算地质 材 料 和 岩 土工 程 的力 学行 为 15 一 8] . 根据 观 测 结果 , L 34 14 综 采 工 作 面横 向变 形 很小 , 因 此本 次计算 采 用平面 应变 模 型假 设 , 即垂直 于计算 剖 面方 向 的变形 为零 , 充分 模拟 研究 了工 作面 开采 过程 中 岩层 的力学状 态 、 破坏 场 的分布 以及位 移 场 的变 化 . 图 2 为 数值 模拟 计 一 算 的平 面模 型 . 由数 值模 拟 分析 结果 可 知 : 3L 4 14 综 采 工 作 面 开 采 中最 大 主 应 力 场 的分 布 服 从 由上 至 下逐 渐 增 加 的趋 势 , 由 于 受 到 工 作 面 开 采 扰 动 的 影 响 , 预 留煤 柱和采 空 区 交接 区 域 会发 生应 力集 中 4 开 采 现 场 观 测 .4 1 矿 压观 测 矿 压观 测 的 目的是认 识和 总结 L 341 4 工 作 面 开采 中矿 山压 力显现 的 基本 规律 , 包 括来 压显 现 特 点和 周 期来压 步距 和 来 压 强度 , 为水体 下采煤 的参 数分析 提 供技 术数据 和 参考 依据 . 沿 工 作 面 倾 向设 了 5 条观 测线 , 分别 距运 输顺 槽 7 . 5 , 45 , 87, 13 5和 1 74 m . 图 3 为平行 切眼 段矿 压显 现规律 图 , 它反 映 了支 架 阻力 与工 作 面推 进度 的相互 关系 . 矿 压 观测 结 果表 明 , 工 作 面支 架 阻 力 、 初撑 力和 周期 来压步 距在 开采 过程 中变 化 不大 , 周 期来压 步距 为 7一 14 m , 其 值 与相似 材 料模 拟实 验所 得结 果 一 致 . .4 2 地 表移 动 规律 观测 研 究 在 L 3 4j 4 综 采工 作 面地表 沿 其走 向和 倾 向布 置 了 两条 观测 线 , 全 程 观测 L 34 14 止 _ 作 而开采 过 程 主 断面 上地 表 移动 与变 形 情况 ( 见 图 4) . 地 表 移 动 观 测 结 果表 明 : 在 整 个 开 采过 程
·574· 北京科技大学学报 2004年第6期 0 有效填堵 30 5 结语 20 通过灵新煤矿西天河下L3414综采工作面试 移梁前初准力 采,安全采出水下压煤55万t,直接经济效益 玉一移架后术阻力 10 1574.1万元,实现了河下压煤的安全开采及环境 116 139162185208231 保护的双重目标,L3414综采工作面在西天河下 工作面推进度/m 图3L3414综采工作面平行切眼段矿压显现规律 最小防水煤岩柱垂高小于79m条件下的成功开 Fig.3 Regularities of ground pressure in L3414 mechaniz- 采,为指导西天河下其它类似条件煤层的安全合 ed mining face parallel to the opening 理开采提供了可靠依据,具有广泛的应用前景. 参考文献 1任奋华,蔡美峰,来兴平,等.采空区覆岩破坏高度监 测分析.北京科技大学学报,2004,26(2):116 走向观测 开切 L41)综采T作面 2康永华,王济忠,孔凡铭,等.覆岩破坏的钻孔观测 方法).煤炭科学技术,2002,30(1226 3北京科技大学,宁夏煤业集团有限责任公司.宁夏 灵新煤矿西天河下安全开采技术研究[R].北京:北 京科技大学,2004 图4地表移动观测布置 4李云鹏,王芝银.回采过程的动态模拟及其应用门. 煤炭学报,1997,22(2):138 Fig.4 Scheme of surface subsidence monitoring 5谢兴华,高延法,基于数值计算的覆岩移动离层量 中,工作面上方及河床均出现不同程度的下沉, 计算方法[J.中国地质灾害与防治学报,2002,13(3): 沿走向观测线上的最大下沉值为1081mm,此处 48 在西天河北岸50m,通过井下涌水量及河流上下 6谢文兵,邓建明,孙厚涛,等.数值模拟在工程实践 中的应用分析[J).矿山压力与顶板管理,2001(1:83 游流量观测,涌水量和河水流量均没有出现异 7刘红元,刘建新,唐春安,采动影响下覆岩垮落过程 常,由此可证明地表水没有遭受破坏,分析其原 的数值模拟).岩土工程学报,2001,23(2):201 因有两个:一是泥岩和粉砂岩类隔水岩层在全基 8朱建明,徐秉业,朱峰,等.FLAC有限差分程序及其 岩柱中所占比例超过60%,其整体隔水性能较 在矿山工程中的应用J.中国矿业,2000,9(4):78 好;二是对地表裂缝及河床下沉区域及时采取了 Comprehensive Analysis of Safe Mining Underneath the Xitian River at Lingxin Coal Mine of China CAI Meifeng,REN Fenhua",LAI Xingping 1)Civil and Environmental Engineering School,University of Science and Technology Beijing.Beijing 100083,China 2)Xi'an University of Science and Technology,Xi'an 710054,China ABSTRACT In order to study the damage of overburden rock and the surface subsidence caused by water-control mining at Lingxin Coal Mine,the in-situ stress distribution law,the damage height and crack height,the overburden rock displacement law and the ground pressure law were abtained by means ofengineering geological investigation, color bore-bole TV inspecting,leakage monitoring in drilling fluid,simulation experiment,numerical computing, in-situ stress measurement and subsidence measurement.The result shows that safe mining at L3414 working face underneath the Xitian River and environmental protection can be implemented ultimately. KEY WORDS ground stress measurement;overburden rock;damage height;mining method;water-control min- ing
5 7 4 . 北 京 科 技 大 学 学 报 2 0 4 年 第 6期 二 、 立 \ 厂~ 移 入 . 八 冬年蕊 一 之 苏声 ~ 移架前 初撑力 ~ 移架后 初撑 力 一 移架 前末阻力 一盆一 移架后 末阻力 盯 二 1 1 1 有 效填 堵 . . `叼贬只国默侧 1 16 13 9 1 6 2 1 8 5 2 0 8 2 3 1 工 作 面推进 度 m/ 图 3 L3 4 14 综 采工 作面 平 行切 眼段矿 压 显现 规律 F gi · 3 R e g u l a r 血i s o f g or u n d P re s s u代 ni L 3 4 1 4 m e c h a n七 - de m ni ni g af c e P a ar l l e l t o tb e 0 P e n in g 一一 暇 塑 走 向观 . 了八 二。 , 砚 鬓 L 34 13 “ 黔 作 , 嘟 硫评丽 ~ 、 、 、 、 谊 擎 尸一一 一一 图 4 地 表 移动 观测 布 置 F ig . 4 S e h e皿 e o f s u r fa ce s u b s id e n ce m o u iot r i n g 中 , 工 作面 上 方 及河 床 均 出现 不 同程 度 的 下沉 , 沿 走 向观 测 线上 的最 大下 沉值 为 1 0 l8 m m , 此 处 在 西天 河 北岸 50 m . 通 过井 下 涌水 量及 河 流上 下 游 流 量 观 测 , 涌水 量 和 河 水 流 量 均 没有 出现 异 常 , 由此 可 证 明地 表 水 没有 遭受破 坏 . 分 析其 原 因 有两 个 : 一 是泥 岩和 粉 砂岩 类 隔水 岩层 在全 基 岩 柱 中所 占 比 例 超 过 60 % , 其 整 体 隔 水 性 能较 好 : 二是 对地 表裂 缝及 河 床下 沉 区 域 及 时采取 了 5 结 语 通 过灵 新煤 矿 西 天河 下 L 3 4 14 综 采 工作面 试 采 , 安 全 采 出水 下 压 煤 5 万 t , 直 接 经 济效 益 15 74 . 1 万 元 , 实 现 了河下 压煤 的安 全 开采 及环 境 保 护 的双 重 目标 . L 3 41 4 综采 工 作 面 在西 天河 下 最 小 防 水煤 岩 柱 垂 高小 于 79 m 条 件 下 的成 功开 采 , 为指 导西 天 河下 其它 类 似条 件煤 层 的安 全合 理 开 采提 供 了可 靠依 据 , 具有 广 泛 的应 用 前 景 . 参 考 文 献 1 任奋 华 , 蔡美 峰 , 来 兴 平 , 等 . 采 空区 覆 岩 破坏 高度 监 测 分析 [J ] . 北 京科 技 大学 学报 , 2 0 0 4 , 26 (2 ) : 1 16 2 康 永华 , 王 济 忠 , 孔 凡铭 , 等 . 覆 岩破 坏 的钻 孔观 测 方法 [J] . 煤 炭 科学 技术 , 2 0 0 2 , 3 0 ( 12 ) : 2 6 3 北京 科技 大学 , 宁夏 煤业 集 团有 限责 任 公司 . 宁夏 灵新 煤矿 西天 河 下安 全 开采技 术 研究 [R] . 北京 : 北 京 科技 大学 , 2 0 04 4 李 云 鹏 , 王芝 银 . 回 采过 程 的动态 模 拟及 其应 用 1[J 煤 炭学 报 , 1 9 9 7 , 2 2 ( 2 ) : 1 3 5 5 谢兴 华 , 高延 法 . 基 于数值 计算 的 覆岩 移动 离层 量 计 算方法 1J[ 、 中国地 质灾 害与 防治 学报 , 2 0 02 , 1 3 (:3) 4 8 6 谢 文兵 , 邓 建 明 , 孙 厚涛 , 等 . 数 值模 拟 在工 程 实践 中的应 用分 析 J[] . 矿 山 压 力与 顶板 管理 , 2 0 0 1 (l ) : 83 7 刘 红元 , 刘 建 新 , 唐 春安 . 采动影 响 下覆岩垮 落过 程 的 数值 模拟 [J ] . 岩 土 工程 学报 , 20 0 1 , 2 3 (2 ) : 2 0 - 8 朱建 明 , 徐 秉业 , 朱 峰 , 等 . F L A C 有 限差 分程序 及 其 在矿 山 工程 中 的应用 IJ] . 中 国矿业 , 20 0 , 9( 4) : 78 C o m P r e h e n s i v e A n a ly s i s o f S a fe M i n i n g U n d e m e a ht ht e X it ian R i v e r a t L i n g x i n C o a l M i n e o f C h i n a CA I 人企沙 ” 岁 , , 尺五厅 矛七月 h u a ,) LA I iX n gP i心 ,2) l ) C i v il an d nE v l r o n ln e n t a l E gn in e e r l n g S c h o o l , U in v e rs ity o f s e i e n e e an d eT e hn o l o g y B e ij l l l g , B e ij ign l 0 0 0 8 3 , C h in a 2 ) X i , an U n i v esr iyt o f S e ien e e an d eT c hn o 1 o g y l 兀 , an 7 1 00 5 4 , C h in a A B S T R A C T nI o r d e r t o s tu dy ht e d am ag e o f o v e r b ur d e n r o e k an d ht e s u ir 益e e s ub s i d e n e e e au s e d b y w a t -er e o n tr o l m i n in g at L i n g x i n C o a l M ine , ht e i n 一 s iut s tr e s s di s itr b u t i o n l a v 叽 het d aJ . ag e h e l g h t an d e r a c k h e i gh t , t h e vo e br u r d e n r o e k d i s P l a e e m ent l aw an d t h e gr o u n d P r e s s uer l aw w e r e ab t alne d b y m e ans o f e n g l n e e r l n g ge o l o g i e al vln e s it g iat on , e o l o r b o er 一 b o l e T V i n s P e e it n g , l e a k a g e m o n iot inr g i n 面l li n g fl u i d , s im u l at i o n e XP e ir m e in , n 切 nr e ir e a l e o m P ut i n g , in 一 s iut s etr s s m e as uemr e in an d s ub s i d e n e e m e as u丁e 幻。 e in . hT e er s u lt s h o w s ht at s a fe m iin n g at L 3 4 14 w o ikr n g acf e un d e m e aht th e X it i an 死v e r an d e vn ior nm e n ta l P fot e e t i o n e an b e im P l e m e in e d u lt im at e l.y K E Y WO R D S gr o un d s tr e s s m e a s uer m e lt : o v e br u r d en r o e k : d am ag e h e i g h t ; m i n i n g m e t h o d ; w at er 一 e o n otr l m i n - in g