煤炭岩溶陷落柱突水预测三维数值模拟

《煤炭岩溶陷落柱突水预测三维数值模拟》由会员分享，可在线阅读，更多相关《煤炭岩溶陷落柱突水预测三维数值模拟（5页珍藏版）》请在金锄头文库上搜索。

1、 -1- 煤炭岩溶陷落柱突水预测三维数值模拟煤炭岩溶陷落柱突水预测三维数值模拟 王月庆，郑飞，周勇 河海大学地质及岩土工程系，南京（210098） E-mail：wangyueqing123 摘摘 要：要：岩溶陷落柱是中国北方型石炭二迭系煤田的一种特殊塌陷，广泛分布于 20 个煤田 45 个煤矿区。由于岩溶陷落柱属于隐伏垂向构造，其导致的突水具有隐蔽性、突发性，且 与岩溶水有天然联系等特点， 对煤矿安全及当地人民生活危害极大。 为研究陷落柱突水机理 及特征， 本文采用三维非稳定流数值模拟计算的理论， 以平朔矿区安家岭井工矿一号井田为 例，模拟预测了陷落柱突水时涌水量，并对其结果进行了分析总结。

2、 关键词：关键词：陷落柱；底板突水；三维数值模拟 1. 引言引言 煤矿突水被认为是安全生产的重大灾害之一， 也是世界产煤国家面对的一大安全开采难 题。我国煤田地质条件十分复杂，受水威胁的煤炭储量占探明储量的 27%，采矿中频繁发 生的突水事故严重威胁着煤矿的安全生产。中奥陶统灰岩层是北方煤田主要的充水含水层， 具有很高的承压水头，岩溶陷落柱的基底一般又均发育其中，若柱体充填物的压密、胶结程 度较差，在采动等外部因素的影响下，陷落柱很可能成为奥灰含水岩层的导水通道，危及矿 井安全1。 我国从“六五”计划将煤矿防治水列入国家重点攻关项目，到“七五”、“八五”煤矿防治水 一期和二期工业性试验的完成，

3、许多专家学者在煤矿涌(突)水量预测研究方面进行了不懈的 探索，取得了丰硕成果，为防治水工程的设计提供了科学依据。早期矿井突水量预测以解析 和统计预测方法为主，此后数值模型、地下水系统及管理在矿井涌水量预测中被广泛应用， 近年来， 广义双重或三重介质渗流模型、 广义双重介质渗流场与应力场耦合模型已经被应用 于矿井涌水量预测中。尽管理论方法及数学模型取得了较大进展，但矿井涌(突)水量预测的 准确性仍然不尽人意。 其中原因不是方法和模型不足以满足预测的需要， 而是基础地质背景 信息不能满足模型的要求， 地质条件的复杂性以及由于认识不清而导致人为模拟简化的不合 理性，是其主要原因1。因此，对地质条件的

4、深刻认识以及含水介质特征的掌握，是矿井涌 (突)水量准确预测的关键。本文仅以平朔矿区安家岭井工矿一号井田为例，系统的研究分析 了研究区水文地质条件并将其概化， 以矿井充水特征及规律的研究为目标， 采用三维非稳定 流模拟计算的理论建立了三维数值模型，预测了陷落柱内涌水量，并对其结果进行分析。 2. 工程概况工程概况 平朔矿区安家岭井工矿一号井位于山西省朔州市平鲁区, 北距平鲁 11km, 南距朔州 15km, 地表大部分被新生界覆盖,属典型的黄土丘陵地貌。地势上北、西、南较高，由透水 的石灰岩组成，东南低，总体为一向东南开口的 C 字型地形。地表水流向为北西南东 方向，地下水流向基本与地表相同。

5、地貌形态为低山丘陵，井田位于神头泉域水文地质单元 的补给迳流区，七里河河床为井田最低侵蚀基准面。研究区范围北起白家辛窑村, 南到 太西村,东起东易村、上窑子旧村,西至西部党家山边煤层露头线, 井田面积 20.21km2,石炭系 上统太原组为本区主要含煤地层, 其中 4#、9#煤层为主要可采煤层2。 井田发育陷落柱 8 个，其中 X1、X3、X4、X5、X7 控制程度较高，被认定为可靠陷落 柱。其特征如表 1 所示2： -2- 表 1 陷落柱特征一鉴表 长轴 短轴 编号 方向 长度 (m) 方向 长度 (m) 面积 (m2) 发育 煤层 X1 NE 50 SW 40 2090 4、9、11 X3

6、 NE 180 NW 60 103904、9、11 X4 SN 135 EW 75 8742 4、9、11 X5 EW 80 SN 50 3841 4、9、11 X7 SW 55 EN 80 229914、9、11 为研究上述陷落柱对煤炭井工矿开采所造成影响的范围及大小， 确保煤田井工矿安全生 产，特展开本次模拟研究。 3. 研究区水文地质概念模型研究区水文地质概念模型 本次模拟计算将平朔矿区安家岭井工矿一号井作为一个统一的水文地质系统， 研究区以 煤田开采范围为界，井田面积约 20.21km2。全区在平面上剖分成 100100 个的矩形网格单 元，垂向上从上往下剖分成八层，山西组为一层；太原

7、组煤系地层剖分成四层，其中 4 煤、 9 煤、各自单独为一层，4煤与 9煤之间砂岩为一层，9煤以下砂岩与 11煤合并成一层； 本溪组为一层；奥陶纪灰岩剖分成二层。每层有效单无格 3176 个，无效单元格 6284 个。整 个模型共计 25408 个有效单元格。 计算目的层为二叠纪山西组砂岩含水层、 石炭纪太原组煤 系含水层、本溪组弱含水层和奥陶纪灰岩含水层，各层均概化为非均质各向异性。煤系含水 层系统和奥灰含水层系统初始流场由现场实测数据给出， 本溪组弱含水层初始流场由煤系含 水层与奥灰含水层初始流场数据插值所得。 整个地下水流态概化为三维非稳定流。 四周均按 通用水头边界处理，系统的顶、底部

8、均为隔水边界。 4. 研究区地下水系统数学模型研究区地下水系统数学模型 根据区内水文地质概念模型，建立下列与之相适应的数学模型34： ()()() 2 0 0 2 ( ,) ( , )( ,) ( ,) ( , ) ( ,) (2) (1) (3) xxyyzz t xxyyzz xyz hhhEh kkkW xxyyzzTt x y z h x y z thx y z x y z hhh kkkq x y z t nnn x y z = += = += 式中： y S E S = 承压含水层 潜水含水层 ； M T B = 承压含水层 潜水含水层 ； M S Ss= ； S 为储水系数；Sy

9、 为给水度；M 为承压含水层单元体厚度（m） ；B 为潜水含水层单元体地 下水饱和厚度（m） ；Kxx，Kyy，Kzz 分别为各向异性主方向渗透系数（m/d） ；h 为点（x， -3- y，z）在 t 时刻的水头值（m） ；W 为源汇项（1/d） ；h0为计算域初始水头值（m） ；h1为第 一类边界的水头值（m） ；Ss为储水率（1/m） ；t 为时间（d） ；为计算域； 1 为第一类边 界； 2 为第二类边界；q（x，y，z，t）为第二类边界上单位面积的侧向补给量（m/d） ；nx 为 x 轴上的单位法向量；ny为 y 轴上的单位法向量；nz为 z 轴上的单位法向量。 5. 模型及计算结果模

10、型及计算结果 有效单元格 无效单元格 图 1 研究区平面剖分图 （从上到下依次为 XL7、XL5、XL4、XL3、XL1） 图 2 研究区纵向剖分图（第 37 列，陷落柱 XL3） 图 3 研究区横向剖分图（第 95 行，陷落柱 XL2） 平朔矿区安家岭井工矿一号井数值模拟如图 1 所示，其长（x 方向）宽（y 方向）高（z 方向）分别为 22370m、17750m、925m。其中，井田内部为有限单元格，外部为无效单元 格， 模型边界概化为通用水头边界。 陷落柱内部水文地质参数由专门针对陷落柱内部进行的 实测抽水试验数据计算而来，如表 2 所示。井田内水文地质参数采用几次抽水试验均值， Kxx

11、，Kyy 为 0.024m/d，Kzz 为 0.0024m/d；储水率为 1.510-8m-1。 -4- 表 2 陷落柱水文地质参数一览表 主轴方向渗透系数（m/d） 储水率 （m-1） 陷落柱 kxx kyy kzz ss XL1 0.8 0.8 0.8 1.010-7 XL3 2.4 2.4 2.4 1.010-7 XL4 0.5 0.5 0.5 1.010-7 XL5 0.8 0.8 0.8 1.010-7 XL7 0.8 0.8 0.8 1.010-7 利用上述模型，采用布置抽水井方法来模拟煤炭井工矿巷道排水，分两个方案分别。将地 下水水位分别疏降到 4#煤底和 9#煤底进行预测。运行

12、模型，得到了各个方案情况下，每个 陷落柱从某一底板向上涌水量数值。如表 3、表 4 所示。 表 3 将地下水疏降到 4#煤底板时各陷落柱底板突水数值一览表(m3/d) 4#煤 砂岩 9#煤 11#煤 本溪组 XL1 5.4624 45.466 53.813 66.402 750.55 XL3 2.3025 12.96 19.872 28.268 302.37 XL4 - 7.1314 7.2657 9.0981 254.19 XL5 - 3.9844 7.3296 9.6452 149.37 XL7 - 3.0775 3.5402 4.1252 49.748 表 4 将地下水疏降到 9#煤底板

13、时各陷落柱底板突水数值一览表(m3/d) 9#煤 11#煤 本溪组 XL1 30.942 72.206 443.55 XL3 15.652 34.756 211.47 XL4 - 6.9853 170.12 XL5 - 10.146 100.93 XL7 - 3.4148 60.742 6. 结论及分析结论及分析 （1）三维数值模拟计算是将整个平朔矿区一号井工矿作为一个统一的水文地质系统，地下 水流态概化为三维非稳定流，四周均按通用水头边界处理，系统的顶、底部均为隔水边界， 模拟合理，计算结果切合实际。说明将三维数值模型模拟法能取得较好效果。 （2）陷落柱的充水受地质构造运动、地下水径流条件、

14、柱体内物质组成、压实和胶结情况 以及承受水压大小等多种条件与因素的控制和影响， 而各种因素又彼此促进和相互制约， 只 有处在现代岩溶水强径流带和集中排泄带并隐伏埋藏在地下水水位以下者， 才能构成突水的 潜在威胁。 由模型预测结果看， 平朔矿区安家岭井工矿一号井中陷落柱对煤炭开采所造成的 水灾害问题不严重。 （3）由于煤层开采过导致应力集中，往往容易产生卸压带。周边围岩又因卸荷、松动等原 因，产生新的裂隙，破坏原有的平衡，加大了陷落柱引起突水的可能性。若岩溶陷落柱一旦 突水，水量大且迅猛，预测及防治难度极大。所以对陷落柱突水问题仍需引起足够重视。 -5- 参考文献参考文献 1尹尚先,武强等北方岩

15、溶陷落柱的充水特征及水文地质模型J 岩石力学与工程学报, 2005, 24 (1) :77-82 2卢世威,许霞平朔矿区安家岭井工矿一号井充水因素分析J 中国煤田地质, 2007,19 (3) :37-69 3 Hsieh P A.Documentation of a com-puter program to simulate horizontal-flow barrers using the U.S.Geologicalsurvey modular three-dimensional fin-ite-difference ground-water flow mod-elR. Washingt

16、on: U.S.Geological Su-rvey Open-File Report,1993:92-477. 4 McDonald MG.Amethod of converting no-flow cells to variable-head cellsfor the U.S.Geological Survey modularfinite-differecne groundwater flowmo-del R. Washington: U.S.Geological Survey Open-File Report,1992:91-536. Three-dimensional numerical simulation of Water inrush forecast in coal