西部干旱矿区采煤引起潜水损失量的定量评价方法.pdf

第 42 卷第 1 期煤炭学报Vol42No1 2017 年1 月 JOUNAL OF CHINA COAL SOCIETYJan2017 赵春虎， 虎维岳， 靳德武西部干旱矿区采煤引起潜水损失量的定量评价方法 J 煤炭学报， 2017， 42 1 169174doi 1013225/j cnkijccs20165001 Zhao Chunhu， Hu Weiyue， Jin Dewu of quantitative uation on amount of groundwater loss from unconfined aquifer caused by mining disturbance in the arid area of Western China J Journal of China Coal Society， 2017， 42 1 169174doi 1013225/jcnkijccs 20165001 西部干旱矿区采煤引起潜水损失量的定量评价方法 赵春虎 1， 2， 虎维岳1， 2， 靳德武1， 2 1. 中国煤炭科工集团西安研究院有限公司， 陕西 西安710054; 2. 陕西省煤矿水害防治技术重点实验室， 陕西 西安710177 摘要 针对西部干旱矿区煤炭资源开发对地下水扰动的相关研究定量化程度不高问题， 通过地下 水动力学分析方法， 得出采掘扰动形成的覆岩垮落带、 裂隙带、 弯曲带及地面沉陷区的空间尺度是 控制地下水损失的主要因素， 其扰动面积越大、 冒裂带高度越高、 弯曲带保护层渗透能力越强则地 下水渗漏量越大; 提出了覆岩导水裂缝边界化， 弯曲带渗透能力分区， 地面沉陷区重新剖分的采掘 扰动影响下地下水系统数值评价模型的构建方法; 并通过构建案例矿井的评价模型， 计算得出神东 矿区补连塔煤矿由于煤矿开采引起松散层潜水损失量为 1. 90104m3/d。 关键词 地下水; 松散含水层; 定量评价; 数值模拟; 导水裂缝带 中图分类号 TD823文献标志码 A 文章编号 02539993 2017 01016906 收稿日期 20160711修回日期 20161109 责任编辑 许书阁 基金项目 国家自然科学基金资助项目 41402220 ; 国家 “十三五” 重点研发计划资助项目 2016YFC0501102 ; 中煤科工西安研究院自筹基 金资助项目 2014MS009 作者简介 赵春虎 1980 ， 男， 陕西扶风人， 副研究员， 博士。Email zhaochunhu cctegxian. com of quantitative uation on amount of groundwater loss from uncon- fined aquifer caused by mining disturbance in the arid area of Western China ZHAO Chun- hu1， 2， HU Wei- yue1， 2， JIN De- wu1， 2 1Xi an esearch Institute of China Coal Technology Engineering Group Corp， Xi an710054， China; 2Shaanxi Key Laboratory of Prevention and Con- trol Technology for Coal Mine Water Hazard， Xi an710077， China Abstract The quantitative degree of coal resources exploitation on water environment disturbance is insufficient in the coal mining area of western ChinaIn this paper， it concludes that the spatial scale is the main factor to control the loss of groundwater by the of groundwater dynamics analysis， which includes the caving zone， the fracture zone， the bending zone and the ground subsidence areaThe of numerical uation model of groundwater system is pro- posed， which is based on the boundary of the hydraulic conductivity of overlying rock， the partition of bending zone， and the re- subdivision of surface subsidence areaAccording to the quantitative uation of groundwater dynamic sys- tem disturbance in the Bulianta coal mine， the simulation results show that the loss of groundwater resources is about 1 90104m3/d due to the mining disturbance Key words groundwater; loose aquifer; quantitative uation; numerical simulation; fractured zone 地下水是西部干旱地区水资源主要组成部分， 具 有极其重要的资源功能和生态功能。而随着煤炭资 源开发西进战略的实施， 我国煤炭基地分布与水资源 丰富程度已呈逆向分布。尤其西部干旱矿区在规模 化、 现代化的煤炭资源开发过程中不可避免地造成了 对地下水影响和破坏， 加剧了地区资源性缺水与生态 中国煤炭期刊网 w w w . c h in a c a j . n et 煤炭学报 2017 年第 42 卷 环境破坏， 是西部煤炭工业可持续发展面临的主要问 题之一。由于能源基础结构的不同， 国外对地下水水 动力扰动问题的研究极少， 主要以矿区水污染风险评 价、 关闭矿井修复、 复垦以及相关法规制定等为主要 研究方向 14 。国内众多学者采用不同技术方法针 对西部地区煤炭开采造成的地下水位下降、 地下水水 量损失等水动力响应问题进行了研究56 ， 但是尚未 形成公认的评判方法， 直到 2011 年 2 月， 环保部第一 次正式颁布了环境影响评价技术导则地下水环 境 ， 对于指导我国开展地下水资源与环境保护工作 具有重要的现实意义。 目前水文地质界常用大井法、 降落漏斗法 78 ， 以及煤炭生产管理部门的吨煤排水量法等均以宏观 分析为主， 评价的定量程度不足。近年来数值法在水 资源评价、 预测和管理中得到迅速的发展， 使地下水 资源研究从传统的研究方法转到模型研究， 大大提高 了地下水资源评价定量化程度。如刘怀忠 9 、 张 伟 10 、 李莹11 、 白乐12 、 吴喜军13 等基于数值模拟 手段将矿井涌水量或围岩渗透性乘以一定的比例系 数 如 1. 2， 1. 5 倍等 ， 通过变化该类比例系数来预 测地下水动态变化， 达到定量评价采煤对地下水影响 程度的目的， 其中涌水量系数和围岩渗透系数的给定 具有较强的主观性， 忽视了采煤导致地下水响应机制 问题。本文将采矿覆岩破坏与地下水系统模拟研究 结合， 通过分析采煤对地下水扰动的控制因素， 探索 数值化的处理方法， 建立定量化程度高的地下水数值 评价模型， 为我国西部干旱矿区煤水资源协调开发 提供科学依据。 1地下水环境扰动机理 1. 1采动覆岩的破坏规律 采煤活动导致围岩体含、 导水能力的变异是地下 水水动力变化的根本原因。目前， 采矿与地下水研究 的相关学者、 技术人员对采动覆岩破坏有较为一致的 认识， 以经典的采矿 “覆岩分带” 14 理论为研究基础， 认为煤层开采后采动覆岩由下至上依次分为垮落带、 裂隙带和弯曲带， 同时在地表形成地面沉陷区。因 而， 典型采动覆岩扰动 图 1 可以总结为“三带一 区” 的变形损伤规律。 1. 2对地下水影响机制分析 陕北与神东煤炭基地位于我国西北干旱地区， 松 散层地下潜水是区域惟一具有供水意义和重要生态 价值的水资源15 。而采矿扰动使赋存地下水的岩层 发生“三带一区” 损伤变形， 由此产生的松散层地下 水向采掘空间的渗漏流失。因此， 须以采矿覆岩破坏 图 1采动覆岩影响示意 Fig. 1Influence of mining overburden rock 为研究基础， 分析采矿引起松散含水层地下水的渗漏 损失问题。 煤层开采导致的地下水损失 渗漏 量主要是由 松散层地下水从垂向上向采掘空间的渗漏量 Q1 、 冒 裂带周边揭露含水层侧向排泄量 Q2和地面沉陷积水 引起无效蒸发量 Q3三部分组成。 1 地下水渗漏量 Q1 深埋煤层开采冒裂带一般未发育至松散含水层 底部， 当余留一定的基岩保护层厚度时， 为了便于分 析扰动因素对地下水渗漏的影响， 认为采掘扰动形成 的导水裂缝为一倒梯形， 如图 2 所示， 定义覆岩体顶 部截面积为 A m2 ; 冒裂带 导水裂缝带 发育高度 为 L1 m ; 松散孔隙含水层与煤层之间基岩厚度 为 M m ; 其渗透系数为 K m/d ; 松散含水层与保 护层界面地下水水压力为 P MPa ; u 为保护层中地 下水的实际流速; r 为地下水容重 N/m3 。 图 2地下水环境扰动示意 Fig. 2Disturbance of groundwater system 定义煤层顶板为 “0” 基准面， 根据经典地下水动 力学分析方法， 得导水裂缝带顶部 L1高度处 的水 头值为 H1 L 1 P/r u2/2g 由于井下排水， 冒裂带内直接与大气联通， 孔隙 压力为大气压， 即 P0， 则 H1 L 1 u 2 /2g 保护层顶部 L2高度处 的水头值为 071 中国煤炭期刊网 w w w . c h in a c a j . n et 第 1 期赵春虎等 西部干旱矿区采煤引起潜水损失量的定量评价方法 H2 L 2 P/r u2/2g 则松散含水层向采掘空间渗漏量 Q1为 Q1 KA L2 P/r L1 M L1 1 由于 L2 L 1ML1 ， 则 Q1 KA 1 P M L1 r 2 根据式 2 可以看出， 渗漏量 Q1与采动导水裂缝 带高度 L1、 弯曲带渗透能力 K、 导水裂缝带揭露范围 A 以及含水层水压 P 呈正相关关系， 与余留的保护层 厚度 ML1 呈负相关。即扰动揭露的面积越大、 冒 裂带高度越高、 弯曲带保护层渗透能力越强则地下水 渗漏量越大。 2 侧向排泄量 Q2 浅埋煤层开采冒裂带一般直接发育至松散含水 层内或地表， 即保护层厚度为 0， 采掘空间顶部含水 层地下水直接疏干。地下水资源流失以侧向排泄为 主， 根据基于稳定流分析的“大井法” ， 含水层地下水 向采掘空间的侧向排量泄 Q2为 Q2 1. 366K 2H S S lg 0 lg r0 3 式中， K 为渗透系数， m/d; H 为水头高度， m; S 为由 于矿井排水而产生的水位降深值， m; 0为引用影响 半径， m; r0为假想 “大井” 的半径 即为导水裂缝揭露 区域的引用半径 ， m。 假设采动裂缝揭露的松散含水层均质无限分布， 其天然水位近似水平， 因此引用影响半径 0可采用 下式计算 0 r 0， 10槡 S K， r0F/ 槡 ， 其中， F 为 导水裂缝影响的含水层面积， m2。从式 3 可以得 出， 导水裂缝影响的含水层面积 F 越大， 地下水侧向 排泄量 Q2越大。 3 地下水无效蒸发量 Q3 深埋煤层开采一般在地面形成沉陷区， 当松散层 地下水埋深较浅时可形成地面积水， 地下水由潜水蒸 发转化成蒸发能力极强的水面蒸发， 则产生了地下水 无效蒸发量 Q3 ， 定义松散层潜水的蒸发强度为 E1， 水面蒸发为 E2， 假设地面已形成地面沉陷积水， 积水 面积为 A2 m2 ， 则 Q3 E2 E 1 A2 4 根据式 4 可以看出， 无效蒸发量 Q3与地区水面 蒸发能力 E2和积水面积 A2正相关。从采掘扰动角度 分析， 积水面积与采掘煤层具体地质条件、 含水层水 文地质条件相关， 即采煤厚度越大、 面积越大、 水位埋 深越小则更易形成地面积水。 综上以上 3 种形式松散层地下水损失情况， 导水 裂缝带影响的高度、 平面范围、 弯曲带保护层厚度及 渗透能力、 地面沉陷范围等是控制地下水资源漏失的 主要因素， 即采掘扰动形成的“三带一区” 的空间尺 度是地下水扰动程度的控制因素。 2采煤对地下水扰动的评价模型构建 由前面分析可知， 采掘扰动形成的导水裂缝带、 弯曲带、 地面沉陷区的空间尺度是地下水损失的控制 要素， 下面就“三带一区” 在地下水系统数值模型数 值化处理技术展开研究。 1 导水裂缝带数值化处理技术 从地下水系统角度分析， 当导水裂缝发育到某一 高度 层位 ， 在冒裂带发育范围内的含水层会被直 接沟通， 采空区顶部含水层地下水体沿垂向采动裂缝 迅速涌入采空区， 在无垂向无越流或大气降水直接补 给时顶部含水层被迅速疏干。 如图 3 所示， 在导水裂缝带外围边界处， 形成了 切割含水层的地下水排泄“带” ， 地下水沿排泄“带” 涌入采空区形成了较为稳定的矿井正常涌水， 造成地 下水水位在排泄“带” 处下降至该含水层底板， 其水 位高度等于该处的含水层底板标高。 图 3地下水与导水裂缝带关系示意 Fig. 3elationship between groundwater and water flowing fractured zone 在导水裂缝带上方外围边界处， 当导水裂缝未揭 露上覆的含水层， 上部的含水层在垂向上可能形成稳 定的越流补给， 由于导水裂缝直接与大气联通， 其孔 隙水压力为大气压力， 即水头高度 Hc近似等于该处 的位置水头标高 冒裂带高程 Hm 为 Hc H m 从水文地质条件概化的角度， 导水裂缝切割含水 层的接触带构成了地下水系统的“内边界” ， 可以概 化为地下水运动的一类水头边界条件， 且水头边界水 头值即为该处含水层底板标高。以此为约束条件， 在 地下水系统数值模型中将导水裂缝切割含水层的接 触带数值化处理成地下水系统的一类内部边界， 通过 识别含水层底板标高， 来定义水头边界的水位值， 以 171 中国煤炭期刊网 w w w . c h in a c a j . n et 煤炭学报 2017 年第 42 卷 计算模拟含水层地下水的变化特征。 2 “弯曲带” 覆岩渗透能力分析与数值化处理 技术 浅埋煤层开采导水裂缝一般直接发育至地表， 覆 岩为垮落带与裂隙带“两带型” 采动裂缝， 深埋煤层 开采易形成典型的“三带一区型” 结构， 弯曲带虽未 形成具有较强导水能力的扰动裂缝， 地应力状态变化 使覆岩发生不同程度与形式变形。如图 4 所示， 位于 采掘空间上方弯曲带的压缩区其岩体渗透能力相对 减小， 两侧拉伸区岩体渗透能力相对减小。根据西部 矿区相关研究成果， 文献 16在神东矿区利用地球 物理手段探测出采前采后弯曲带覆岩含水性整体较 为稳定， 文献 1718 利用原位钻孔压水试验得出位 于采动拉伸区的离石组黄土渗透能力较采前增大了 12 个数量级， 而同位于压缩区与拉伸区的保德组红 土， 其采前采后渗透能力相当的结论。可见弯曲带渗 透能力变异对地下水漏失影响小。 图 4弯曲带渗透能力分区 Fig. 4Curve of the capacity of the bending zone 在地下水系统建模时， 由于采掘扰动后弯曲带覆 岩仍为连续介质， 可以根据渗透能力在采掘前后的变 化趋势， 进行渗透能力参数分区， 即划定渗透能力变 化趋势一致的区段， 进行参数分区， 分区内的渗透能 力参数进行重新赋值， 达到采掘扰动下渗透能力变化 的数值化处理目的。 3 地面沉陷区数值化处理技术 地面沉陷降低了潜水水位埋深， 易形成地面积 水， 从而加大了地下水资源的无效蒸发， 因而在建模 过程中应予以考虑。如图 5 所示， 根据地面沉陷预测 或实测结果， 通过改变地面高程值， 对模型的上边界 进行重新剖分。 需要注意的是在重新改变了上边界条件后， 须对 模型进行试运行， 分析沉陷区最低点的地面标高 Hmin 与潜水位标高 Hw的关系， 当 Hmin H w， 地下水出露于 地表形成积水， 此时需对地面进行蒸发能力的重新分 区， 积水区范围按照水面蒸发处理， 其他区域仍按照 潜水蒸发处理; 当 HminHw， 地下水未出露， 全区仍为 图 5地面沉陷数值化处理示意 Fig. 5Numerical treatment of ground subsidence 潜水蒸发， 但是由于埋深的降低， 潜水蒸发量可能增 大。 综上， 以地下水系统数值模拟为基础， 通过将采 煤扰动形成的覆岩导水裂缝边界化、 弯曲带渗透能力 分区， 地面沉陷区重新剖分的数值化处理手段是地下 水系统数值评价模型构建的基本思路。 3评价案例 3. 1案例矿井概况 补连塔井田位于西部神东矿区中部， 主采煤层埋 深小于 150 m， 地表松散层厚度约为 30 m， 采用综合 机械化采煤工艺。如图 1 所示， 井田内松散层类孔隙 水主要接受大气降水的入渗补给， 在沟谷地下水出露 于地表， 形成常年或季节性的地表水体 补连沟、 活 鸡兔等 ， 乌兰木伦河是区内所有地表水的排泄区。 如图 6 所示， 模型在垂向上按照岩性变化剖分为 4 层， 研究区水文地质单元较为完整， 分别以常年流 水的乌兰木伦河、 呼和乌苏沟、 活鸡兔沟为模型的一 类水头边界。 图 6研究区 补连塔井田 范围示意 Fig. 6Sketch map of study area 3. 2采掘扰动影响下地下水系统模拟 据矿井实测， 补连塔井田采动覆岩为典型的“两 271 中国煤炭期刊网 w w w . c h in a c a j . n et 第 1 期赵春虎等 西部干旱矿区采煤引起潜水损失量的定量评价方法 带” 破坏形式， 不存在弯曲带渗透能力变异与地面沉 陷积水问题， 因此在评价模型中定义导水裂缝接触带 为模型的“内边界” ， 概化为含水层地下水运动的一 类水头边界条件， 松散含水层水头边界水头值为该处 松散含水层底板标高， 降雨补给以研究区平水年数据 为准。 1 水位 由图 7 a 分析可知 蓝色箭头为地下水流向， 红 色箭头为地下水转地表水流向 ， 煤层未开采前， 表 现出大气降水和地下水补给地表水的三水转化特征， 乌兰木伦河为井田的地下水排泄点。局部由于地形 控制， 以地表分水岭为界限地下水分别向两侧的补连 沟、 呼和乌素沟和活鸡兔沟等地表水体排泄。 图 7地下水流场变化特征 Fig. 7Groundwater level contour map under mining disturbance 补连塔井田规模化开采开始于 2000 年左右， 如 图 7 b 所示， 2012 年井田基本以 1 2 煤回采完毕的 三、 四盘区为中心， 形成较为明显的降落漏斗， 水位下 降了 20 m 以上， 采空区顶部含水层直接被疏干 黄 色区域 。 2 地下水均衡影响 表 1 中， 未采掘前大气降水是地下水补给的主要 来源， 侧向补给次之， 河流仅在局部地段补给地下水 存在补给。地下水向河流等地表水体的排泄是主要 排泄形式， 蒸发排泄次之， 由于模拟区水文地质单元 较为完整， 地下水向区外侧向排泄量为 0。 表 1天然条件下地下水均衡 Table 1Groundwater balance under natural conditions 源汇项数值占比/合计 降雨补给量/ 104m3d 1 6. 7171. 69 河流补给量/ 104m3d 1 0. 707. 469. 36 侧向补给量/ 104m3d 1 1. 9520. 85 蒸发排泄量/ 104m3d 1 3. 0732. 82 9. 36 向河流排泄量/ 104m3d 1 6. 2967. 18 如图 8 所示， 随着井田采掘范围逐年增加， 受采 掘扰动范围扩大， 散含水层向采掘区段排泄的地下水 水量逐年增大， 以 2012 年为例， 地下水向采掘区段的 侧向排泄量 即地下水的漏失量 为1. 9104m3/d， 较 之未采掘前， 其排泄量主要是袭夺了河流排泄减量和 潜水蒸发的减量。 图 8采掘影响下水均衡变化 Fig. 8Water balance under mining disturbance 3 对比分析 2012 年补连塔井田四盘区的实测矿井正常排水 量约为 0. 852104m3/d 不含采空积水量 ， 由于矿井 排水量不包括补给至采空区积水和井下复用等水量， 其结果相对偏小。根据井田水文地质特征， 通过大井 法计算地下水向采掘区的侧向排泄量， 其中松散含水 层 K 取 2. 385 m/d， 水头高度、 水位降深均取含水层 厚度值 30 m， F 取截止 2012 年已采掘的三四盘区面 积约为 2. 49107m2， 计算可得松散含水层向采掘区 的地下水排泄量为 4. 4104m3/d， 由于大井法为四面 汇水条件， 无限含水层， 规则影响半径导致其计算结 果明显偏大， 因此数值评价模型计算结果 1. 9 104m3/d 能正确、 定量地反映出采掘扰动引起地下 水的损失量。 4结语 1 “三带一区” 是地下水赋存介质扰动的主要 表现形式 由采动 “覆岩分带” 理论可知， 采掘扰动形 371 中国煤炭期刊网 w w w . c h in a c a j . n et 煤炭学报 2017 年第 42 卷 成了覆岩的垮落带、 裂隙带及弯曲带。 2 “三带一区” 空间尺度是控制地下水损失的 主要因素 通过地下水动力学分析可知， 扰动影响的 面积越大、 冒裂带高度越高、 弯曲带保护层得到渗透 能力越强则地下水渗漏量越大。 3 提出了 “三带一区” 模型化处理技术 在地下 水系统数值模型中将导水裂缝切割含水层的接触带 数值化处理成地下水系统的一类内部边界; 根据弯曲 带渗透能力在采掘前后的变化趋势， 进行渗透能力参 数分区; 根据地面沉陷区预测或实测结果， 通过改变 地面高程值， 对模型的上边界进行重新剖分。 4 通过构建神东矿区补连塔井田地下水扰动 的计算机评价模型， 模拟得出由于采掘扰动使井田地 下水水位下降 20 m 以上 至 2012 年 ， 地下水资源 损失量可达 1. 90104m3/d。 参考文献 eferences 1Chugh Y P， Behum P TCoal waste management practices in the USA An overview J International Journal of Coal Science Tech- nology， 2014， 1 2 163176 2Skousen J， Zipper C E Post- mining policies and practices in the Eastern USA coal region J International Journal of Coal Science Technology， 2014， 1 2 135151 3胡振琪， 龙精华， 王新静论煤矿区生态环境的自修复、 自然修复 和人工修复 J 煤炭学报， 2014， 39 8 17511757 Hu Zhenqi， Long Jinghua， Wang XinjingSelf- healing， natural resto- ration andartificialrestorationofecologicalenvironmentfor coal miningJ Journal of China Coal Society， 2014， 39 8 1751 1757 4Manca P P， Desogus P， Orr GThe reuse of abandoned Acquare- si mine voids for storage of the Masua flotation tailingsJInterna- tional Journal of Coal Science Technology， 2014， 1 2 213220 5范立民， 马雄德， 蒋辉， 等西部生态脆弱矿区矿井突水溃沙危险 性分区J煤炭学报， 2016， 41 3 531536 Fan Limin， Ma Xiongde， Jiang Hui， et alisk uation on water and sand inrush in ecologically fragile coal mine JJournal of Chi- na Coal Society， 2016， 41 3 531536 6王力， 卫三平， 王全九榆神府煤田开采对地下水和植被的影响 J煤炭学报， 2008， 33 12 14081414 Wang Li， Wei Sanping， Wang Quanjiu Effect of coal exploitation on groundwater and vegetation in the Yushenfu Coal MineJJour- nal of China Coal Society， 2008， 33 12 14081414 7潘天杭用降落漏斗体积法求潜水含水层中抽水时之影响半径 J水文地质工程地质， 1957 9 4142 Pan TianhangThe influence radius of the water pumping in the aqui- fer by the of the volume of the descent funnel J Hydrogeology and Engineering Geology， 1957 9 4142 8张子文稳定降落漏斗中心最大下降及其体积计算方法探讨 J工程勘察， 1989 3 3436 Zhang Ziwen Discussion on the maximum descent and its vol- ume calculation of the stable descend funnel center in JGeotechnical Investigation Surveying， 1989 3 3436 9刘怀忠煤矿开采对矿区地下水系统扰动的定量评价研究D 徐州 中国矿业大学， 2011 Liu Huaizhong Quantitative uation of groundwater system dis- turbance in mining areaD Xuzhou China University of Mining and Technology， 2011 10李杨浅埋煤层开采覆岩移动规律及对地下水影响研究D北 京 中国矿业大学 北京 ， 2012 Li YangStudy on the movement law of overlying strata in shal- low coal seam mining and its influence on groundwaterDBei- jing China University of Mining and Technology Beijing ， 2012 11李莹陕北煤炭分布区地下水资源与煤炭开采引起的水文生态 效应D西安 长安大学， 2008 Li YingThe hydrological and ecological effects caused by groundw- ater resources and coal mining in the coal distribution area of Northern Shaanxi ProvinceDXi an Chang an University， 2008 12白乐， 李怀恩， 何宏谋， 等煤矿开采区地表水地下水耦合模拟 J煤炭学报， 2015， 40 4 931937 Bai Le， Li Huaien， He Hongmou， et alIntegrated simulation of sur- face water and groundwater in a high intensive coal mining area JJournal of China Coal Society， 2015， 40 4 931937 13吴喜军煤炭开采地区河道径流变化与生态基流研究以陕北 窟野河流域为例D西安 西安理工大学， 2013 Wu XijunStudy on river runoff change and ecological base flow in coal mining areaDXi an Xi an University of Technology， 2013 14钱鸣高， 石平五， 许家林矿山压力与岩层控制M徐州 中国 矿业大学出版社， 2010 98101 Qian Minggao， Shi Pingwu， Xu Jialin Mining pressure and stra- ta controlMXuzhou China University of Mining and Technology Press， 2010 98101 15王双明， 黄庆享， 范立民， 等生态脆弱矿区含 隔 水层特征及 保水开采分区研究J煤炭学报， 2010， 35 1 714 Wang Shuangming， Huang Qingxiang， Fan Limin， et al Study on overburden aquclude and water protection mining regionazation in the ecological fragile mining areaJJournal of China Coal Soci- ety， 2010， 35 1 714 16Du W， Peng S， Zhu G， et alTime- lapse geophysical technology- basedstudyonoverburdenstratachangesinducedby modern coal miningJInternational Journal of Coal Science Technology， 2014， 1 2 184191 17隋旺华开采沉陷土体变形工程地质研究M徐州 中国矿业 大学出版社， 1999 Sui Wanghua The engineering geology of the soil deation of mining subsidenceMXuzhou China University of Mining and Technology Press， 1999 18李涛陕北煤炭大规模开采含隔水层结构变异及水资源动态研 究D徐州 中国矿业大学， 2012 Li TaoStudy on the structure variation and water resources of the large- scale mining of coal in the coal of Northern ShaanxiD Xuzhou China University of Mining and Technology， 2012 471 中国煤炭期刊网 w w w . c h in a c a j . n et