《2022年2022年矿坑涌水量的常用预测方法汇总 .pdf》由会员分享，可在线阅读，更多相关《2022年2022年矿坑涌水量的常用预测方法汇总 .pdf（27页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、吉林大学精品课 专门水文地质学 教材 专门水文地质学10.4矿坑涌水量预测一、矿坑涌水量预测的内容、方法、步骤与特点（一）矿井涌水量预测的内容及要求矿坑涌水量预测是一项重要而复杂的工作，是矿床水文地质勘探的重要组成部分。矿坑涌水量是指矿山开拓与开采过程中，单位时间内涌入矿坑(包括井、巷和开采系统)的水量。通常以m3/h 表示。它是确定矿床水文地质条件复杂程度的重要指标之一，关系到矿山的生产条件与成本，对矿床的经济技术评价有很大的影响。并且也是设计与开采部门选择开采方案、开采方法，制定防治水疏干措施，设计水仓、排水系统与设备的主要依据。因此，在矿床水文地质调查中，要求正确评价未来矿山开发各个阶段
2、的涌水量。其内容与要求包括可概括为以下四个方面：（1）矿坑正常涌水量：指开采系统达到某一标高(水平或中段 )时，正常状态下保持相对稳定的总涌水量，通常是指平水年的涌水量。（2）矿坑最大涌水量：是指正常状态下开采系统在丰水年雨季时的最大涌水量。对某些受暴雨强度直接控制的裸露型、暗河型岩溶充水矿床来说，常常还应依据矿山的服务年限与当地气象变化周期，按当地气象站所记录的最大暴雨强度，预测数十年一遇特大暴雨强度产生时，可能出现暂短的特大矿坑涌水量，作为制订各种应变措施的依据。（3）开拓井巷涌水量：指包括井筒(立井、斜井 )和巷道（平、平巷、斜巷、石门)在开拓过程中的涌水量。（4）疏干工程的排水量：是指
3、在规定的疏于时间内，将一定范围内的水位降到某一规定标高时，所需的疏干排水强度。对于地质勘探阶段来说，主要是进行评价性的计算，以预测正常状态下矿坑涌水量及最大涌水量为主。 至于开拓井巷的涌水量预测和专门性疏干工程的排水量的计算，由于与矿山的生产条件密切相关，一般均由矿山基建部门或生产部门承担。（二）矿坑涌水量预测的方法根据当前矿床水文地质计算中常用的各种数学模型的地质背景特征极其对水文地质模型概化的要求，可作如下类型的划分：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 2
4、7 页 - - - - - - - - - 混合型模型水均衡法有限差法有限元法数值解非稳定井流公式稳定井流公式井流方程解析解确定模型回归方程曲线方程非确定性统计模型数学模型分类sQ（三）矿坑涌水量预测的步骤矿坑涌水量预测是在查明矿床的充水因素及水文地质条件的基础上进行的。它是一项贯穿矿区水文地质勘探全过程的工作，一个正确预测方案的建立，是随着对水文地质条件认识的不断深化，不断修正、完善而逐渐形成的，一般应遵循如下三个基本步骤：1. 选择计算方法与相应的数学模型详勘阶段均要求选择2 个或 2 个以上的计算方法，以相互检验， 映证。 选择时必须考虑三个基本要素：（1）矿床的充水因素及水文地质条件复
5、杂程度。如：位于当地侵蚀基准面之上，以降水入渗补给的矿床，应采用水均衡法；水文地质条件简单或中等的矿床，可采用解析法或比拟法；水文地质条件复杂的大水矿床，要求采用数值方法；（2）勘探阶段对矿坑涌水量预测的精度要求；（3）勘探方法、勘探工程的控制程度与信息量：如：水均衡法，要求不少于一个水文年的完整均衡域的补给与排泄项的动态资料；Qs 曲线方程外推法， 要求具抽水试验的水位降达到预测标高水柱高度的1/21/3；解析法，要求勘探工程全面控制含水层的非均质各向异性、非等厚的结构特征及其边界条件与补给、径流与排泄， 并提供数值模型的建立、识别、预测所需的完整信息数据，这些数据的获取，只有采用大型抽、放
6、水试验对渗透场进行整体控制与揭露才可能做到。因此，计算方法与相应数学模型类型的选择，与矿床的充水因素及水文地质条件复杂程度、勘探方法勘探工程的控制程度及信息量是相互关联的，统一在最佳技术经济条件这一原则下，所以数学模型类型选择是否合理，可以用以下标准衡量：一是对矿床水文地质条件的适应性：指能否正确刻画水文地质条件的基本特征；二是对勘探方法勘探工程控制程度的适应性：指是否最充分的利用勘探工程提供的各种信息， 即信息的利用率；同时， 也可理解为所选数学模型要求的勘探信息是否有保证，即信息的保障率。2. 构造水文地质模型矿坑涌水量预测中数学模型的作用，是对水文地质条件进行量化，因此预测精度主要取决于
7、对充水因素与水文地质条件判断的准确性，由于不同数学模型类型对水文地质条件的刻画形式与功能各异，因此必须按数学模型的特点构造水文地质模型，称水文地质条件概化。概化后的水文地质模型称水文地质概念模型，它在地质实体与数学模型之间起中介桥梁作用。下面以最基本的预测方法解析法与数值法为例作一讨论。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 27 页 - - - - - - - - - （1） 概化已知状态下的水文地质条件；（2） 给出未来开采状态下的内边界条件；（3） 预测未来
8、开采状态下的外边界条件。解析法将复杂的含水层结构与内外边界，以理想化模式构造理论公式，因此必须按解析解要求进行概化。如含水层均质等厚，内外边界几何形态规则，边界供水条件简单、确定。数值法以近似分割原理对复杂的含水层结构、内外边界条件进行量化“逼真”，概化时要求以控制水文地质条件与内外边界的节点参数、水位与流量来构造水文地质概念模型。随着数学模型研究的不断进展，现代水文地质计算对水文地质模型的要求越来越高。目前，对复杂的大水矿床来说，一个可靠的水文地质模型的建立，必须贯穿整个勘探过程，并大致经历三个阶段。即：第一阶段：通过对以往资料的整理，提出水文地质模型的“雏型”，作为下一步勘探设计的依据。尤
9、其对大型抽（放）水试验的设计具有特殊重要的意义；第二阶段：根据进一步勘探提供的各种信息数据。特别是大型抽(放)水资料，通过流场分析或数值模拟，完成对“雏型”模型的调整，建立水文地质模型的“校正型”；第三阶段：在“校正型”的基础上，按开采方案给出疏干工程的内边界条件，根据勘探资料预测不同疏干条件下的外边界条件，建立水文地质概化模型的“预测型”。3. 计算数学模型，评价预测结果应该指出，不能把数学模型的解仅仅看作是一个单纯的数学计算，而应看作是对水文地质模型和数学模型进行全面验证识别过程，也是对矿区水文地质条件从定性到定量在回到定性的不断深化的认识过程。（四）矿坑涌水量预测的特点虽然矿坑涌水量预测
10、的原理方法与供水水资源评价类同，但其预测条件、 预测要求与思路各有不同。如：（1）供水水资源评价，以持续稳定开采确保枯水期安全开采量为目标；而矿坑涌水量预测则以疏干丰水期的最大涌水量为目标；（2）矿床大多分布于基岩山区。含水质的非均质性突出。参数代表性不易控制；边界条件复杂、非确定性因素多，常出现紊流、非连续流与管道流。定量化难度大。（3）矿山井巷类型及其分布千变万化，开采方法、开采速度与规模等生产条件复杂且不稳定，与供水的取水建筑物简单、分布有序、生产稳定形成显明对比，给矿坑涌水量预测带来诸多不确定性因素。（4）矿坑涌水量预测多为大降深下推。此时开采条件对水文地质条件的改变难以预料和量化，这
11、与供水小降深开采有明显差异。（5）矿床水文地质勘探从属于矿产地质勘探，与专门性的供水水文地质勘探对比，前者一般投入小、工程控制程度低，预测所需的信息量相对少而不完整。以上特点， 决定了矿坑涌水量预测中存在诸多产生误差的客观条件。因此属于评价性计算，为矿山设计及采前进一步专门性补充勘探提供依据。二、 QS 曲线方程外推法（一）原理与应用条件名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 27 页 - - - - - - - - - 指用稳定井流条件下抽水试验的Q=f （s）
12、方程，外推未来疏干水位降的涌水量。实质上也是一种相似条件下的比拟法。应用时的前提条件是：一、抽水试验建立Q=f （s），应符合稳定井流条件；二、抽水试验的各种条件应与预测对象的疏干条件接近。因此，必须重视试验的技术条件，包括：1.应将抽水试验孔布置在预测对象的分布地段，保证水文地质条件的一致性；2.采用大口径（或孔组）试验，计算时为消除井径对涌水量的影响，需做井径换算；3.抽水降深应大于疏干水位水柱高度的1/21/3，计算时的外推疏干降深不应超过1.75倍的抽水降深，主要考虑疏干状态下的补给条件；4.用枯季抽水试验预测正常涌水量，根据雨季试验预测季节性最大涌水量；5.要排除抽水过程中一切自然和
13、人为随机影响因素的干扰。Qs 曲线方程法的优点是：回避各种水文地质参数求参过程中的失真，计算简单易行。适用于建井初期的井筒涌水量预测。上水平疏干资料外推下水平的涌水量，以及矿床规模小、矿体分布集中、边界条件和含水结构复杂的涌水量预测。（二）计算方法与步骤1. 鉴别 Qs 曲线类型（1）曲度法：即用曲度n 值进行鉴别：1212lglglglgQQSSn，当 n1 时，为直线 Q=qS；1n2 时，为幂曲线Q=bSa；n2 时，为抛物线S=aQ+bQ2；n2 时，为半对数曲线S=a+blgS 。如果 n1 时，表明抽水试验不正确。2. 确定方程参数（1）最小二乘法：应根据Q=f （s）类型选用最小
14、二乘法，如常见的幂函数型：QSQSNSSblglg)lg(lg)lg()(lg22NSbQalg1lglg（2）图解法：即利用直角坐标的图解，a 为图解中纵坐标上所切的截距线段；b 为直线对水平倾角的正切。其它类型详见地下水动力学。3. 井径换算（1）根据稳定井流理论中，不同地下水运动（层流、紊流）状态下井径与涌水量关系进行换算，如层流：）（）（井井孔孔孔井rRrRQQlglglglg； 紊流：孔井孔井rrQQ。试验表明：井径对涌水量的影响，一般比对数大，而小于平方根。（2）用二次或二次以上不同孔径的抽水试验资料建立井径d 的换算公式： Q=mdn，m、n 值可用最小二乘法求得。4. 水均衡论
15、证将预测结果进行疏干条件下的水均衡论证，分析所得的最大涌水量是否符合疏干状态下的边界补给条件。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 27 页 - - - - - - - - - （三）实例1. 广东某金属矿， 曾用 Qs 曲线法预测 +50m 水平的涌水量为14450m3d，与坑道放水试验的外推涌水量(14000m3/d)接近，而解析法预测的结果(12608m3d)则偏小 12。2. 湖南叶花香铜矿，1970 年采用试生产坑道的排水量与水位降资料。分别用两种曲线
16、类型对未来四个中段的涌水量进行外推预测，后经实际检验见表14-2。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 27 页 - - - - - - - - - 计算方法不同降深的水量52.45m73.33m99.47m水量（m3/d）误差（% ）水量（m3/d）误差（% ）水量（m3/d）误差（% ）计算实际计算实际计算实际S=aQ+bQ215145 16937 -10.60 18223 23644 -22.90 21514 31155 -31.00 Q=bSa14840
17、 -12.40 18093 -23.50 21582 -30.70 实际检验表明预测结果与实际情况较接近，产生误差的主要原因是其间曾发生两次大突水淹井，大量泥沙溃入矿坑，地面出现岩溶坍陷和裂缝所致。三、相关分析法（一）原理与应用条件相关分析是根据涌水量与主要影响因素之间相关关系的密切程度建立回归方程，利用抽水试验或开采初期的疏干资料，预测矿坑涌水量或外推开采后期下水平的涌水量。其原理已在供水中介绍。根据实际资料的统计，多元复相关预测远比单相关效果好，其回归方程表达的内容丰富，可反映除降深外的各种影响因素。它的应用条件与Qs 关系方程类同，但对原始数据的采集有严格要求：1 代表性：（规范）要求不
18、少于一个水文年（包括丰、平、枯季节）的动态观测数据，同时数据（择本）量不少于30 个；2 一致性：指应与预测对象上条件相一致；3 独立性与相关性：即多自变量有独立的变化规律，相互间关系不大；而与涌水量之间均存在密切的相关关系，（规范）要求相关系数不低于0.7。（二）实例与计算方法1. 利用勘探阶段抽水试验资料预测矿坑涌水量如广东沙洋矿通过在勘探阶段设计相距6m 的两个抽水孔和十余个不同距离的观测孔组成的群孔抽水试验，取得了复相关计算所需的涌水量Q 与井径 r（是将距抽水孔不同距离观测孔的位置概化为疏干状态下的坑道系统不同面积的作用半径）、水位降 S（即不同作用半径的水位降，以模拟疏干水位降）有
19、关资料，（见表 3）通过求参建立了复相关幂函数预测方程：536.11843.3189.11SrQ其复相系数达0.9468，复相关机误仅0.0721，完全可用于未来矿山各设计水平与面积的矿坑涌水量预测。经实际排水资料检验，预测误差偏小3856%，主要与开采导函大量地面岩溶坍陷有关。表 3 S （m）r （m）Q（L/s）50.53 54.76 125.50 127.90 150.90 202.20 216.50 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 27 页 -
20、- - - - - - - - 34.491 1.147 1.705 55.168 2.053 2.033 1.567 1.784 1.427 1.393 69.145 2.984 2.902 2.116 2.474 2.005 1.861 1.308 2. 利用矿山观测资料外推预测可充分考虑矿坑涌水量的增长和各项生产因素间的关系，并根据它们之间的密切程度来建设涌水量方程。在原苏联顿巴斯煤矿的某些涌水量预测中，首先，在30 个矿井中建立了 320 个观测点，获得了涌水量（Q2）与各生产因素（包括矿产量P0、开采深度H0、开采面积 F0、生产时间T0等）之间的相关关系，以及其密切程度，见表14-
21、4。表 14-4 相关系数 r lgQ0lgP0lgH0lgF0lgT0lgQ00.664 0.451 0.593 0.175 lgP00.664 0.340 0.680 0.323 lgH00.451 0.340 0.559 0.523 lgF00.593 0.680 0.559 0.778 lgT00.175 0.323 0.523 0.778 根据判别得知，生产时间T0对的影响不大（相关系数为0.175） 。用多元复相关计算，求得四元复相关曲线回归方程：11. 0024. 0051. 00072.0FHPQ其复相关系数： 。在此基础上建立了比拟公式：11.0024. 0051.000FF
22、HHPPQQ预测结果与传统的单位涌水量法相比，使误差减少1.4 倍（式中P、H、F 为设计值）。四、解析法（一）解析法的应用条件解析法是根据解析解的建模要求，通过对实际问题的合理概化，构造理想化模式的解析公式，用于矿坑涌水量预测。具有对井巷类型适应能力强、快速、简便、经济等优点，是最常用的基本方法。 解析法预测矿坑涌水量时，以井流理论和用等效原则构造的“大井” 为主，后者指将各种形态的井巷与坑道系统，以具有等效性的“大井”表示，称“大井”法。因此说：矿坑涌水量计算的最大特点是“大井法”与等效原则的应用，而供水则以干扰井的计算为主。稳定井流解析法：应用于矿坑疏干流场处于相对稳定状态的流量预测。包
23、括在已知某开采水平最大水位降条件下的矿坑总涌水量；在给定某开采水平疏干排水能力的前提下，计算地下水位降深（或压力疏降）值。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 27 页 - - - - - - - - - 非稳定解析法： 用于矿床疏干过程中地下水位不断下降，疏干漏斗持续不断扩展，非稳定状态下的涌水量预测。包括：已知开采水平水位降（s） 、疏干时间 （t） ，求涌水量 （Q） ；已知 Q、s，求疏干某水平或漏斗扩展到某处的时间（t） ；已知 Q、t，求 s，以确定
24、漏斗发展的速度和漏斗范围内各点水头函数隨时间的变化规律，用于规划各项开采措施。在勘探阶段，以选择疏干量和计算量最大涌水量为主。（二）计算方法如上所述， 应用解析法预测矿坑涌水量时，关键问题是如何在查清水文地质条件的前提下，将复杂的实际问题概化。它可概括为如下三个重要方面：分析疏干流场的水力特征，合理概化边界条件，正确确定各项参数。1. 分析疏干流场的水力特征矿区的疏干流场是在天然背景条件下，迭加开采因素演变而成。分析时， 应以天然状态为基础，结合开采条件作出合理概化。（1）区分稳定流与非稳定流矿山基建阶段， 疏干流场的内外边界有受开拓井巷的扩展所控制，以消耗含水层储量为主，属非稳定流；进入回采
25、阶段后，井巷输廊大体已定，疏干流场主要受外边界的补给条件控制，当存在定水头（侧向或越流）补给条件时，矿坑水量被侧向补给量或越流量所平衡，流场特征除受气候的季节变化影响外，呈现对稳定状态。基本符合稳定的“建模”条件，或可以认为两者具等效性；反之，均属非稳定流范畴。如河北开滦煤矿，其矿坑涌水量随坑道走向的延展而增加，但这种延展暂停时，涌水量立即出现相对的稳定。此时仅表现有受降水的季节变化在一定变幅范围内上下波动，并出现强出水点和边缘出水点成袭夺中间出点现象，而总涌水量不变。又如辽宁复州粘土矿，其涌水量随采深增加，但某一水平进入回采时，其涌水量就逐渐稳定，并保持到下一水平突水进止，在此阶段虽然也出现
26、下水平突水点袭夺上水平突水点现象，但总涌水量却保持不变。由此可见， 在某些矿区的疏干过程中，不仅存在疏干流场的相对稳定阶段，而且隨矿山工程的进展而不断相互转化。但选用稳定流解析法时要慎重，必须进行均衡论证，判断疏干区是否真正存在定水头供水边界或定水头的越流系统。此外，对于地下水储存量较大的矿区，要单独计算疏干过程中不同阶段含水层储存量的放强度， 大量生产实际表明，矿坑最大涌水量， 并非出现在疏干过程的稳定阶段。（2）区分达西流与非达西流在矿坑涌水量计算时，常遇到非达西流问题，它涉及解析法的应用条件，在宏观上可概括有两种情况：一是暗河管道岩溶充水矿床，地下水运动为压力管道流与明渠流；此外， 分水
27、岭地段的充水矿床， 矿坑涌水量直接受垂向入渗降雨强度控制，与水位降深无关。两者均与解析法的“建模”条件相距甚大，矿坑涌水量预测应选择水均衡法或各种隨机统计方法。二是，局部状态的非达西流，常发生在大降深疏干井巷附近与某些特殊构造部位，它只对参数计算与参数的代表性产生影响。在宏观上， 它是一个流态概化问题，不存在解析法的应用条件问题。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 27 页 - - - - - - - - - （3）区分平面流与空间流严格讲， 在大降深疏干条件
28、下，地下水运动的垂向速度分量不能忽略，均为三维空间流（包括非完整井巷的地下水运动）问题，其分布范围仅限于井巷附近，均为含水层厚度的1.54.75 倍。因此，在矿坑涌水量预测中，大多将其纳入二维平面流范畴，在宏观上不影响预测精度。计算时应根据井巷类型作出不同的概化：如：竖井的涌水量计算，可概化为平面径向流问题，以进流公式表达；计算水平巷道涌水时， 以剖面平面流近似，采用单宽流量解析公式，但其两端上往往也产生辐射流（见图），需要考虑它的存在，并采用平面径向流公式补充计算巷道端部的进水口。坑道系统则复杂得多，根据“大井法”原理，一般以近似的径向流概化，但当坑道系统近于带状的狭长条形时，也可概化为剖面
29、流问题。对于倾斜坑道， 根据阿勃拉莫夫有关水电比拟法的研究，证明坑道的倾斜对涌水量影响不大， 可根据坑道的倾斜度，分别按竖井或水平巷道进行近似。即：若坑道倾斜度450时，视其与竖井近似，用井流公式计算；若坑道倾斜度450时，则视其与水平巷道近似，用单宽流量公式计算。根据解析解的存在条件，一些简单的非完整井巷涌水量计算，可以运用三维空间问题予以解决。此时，可根据非完整井的特点，运用地下水动力学中映射法与分段法的原理来求解。通常用平面分段法解决完整竖井的涌水量计算，用剖面分段法解决非完整平巷的涌水量计算。（4）区分潜水与承压水与供水不同，在降压疏干时，往往出现承压水转化为潜水或承压无压水。此外，在
30、陡倾斜含水层分布的矿区，还可能出现坑道一侧保持原始承压水状态，而另一侧却由承压水转化无压水或承压无压水的现象。概化时，需从宏观角度作等效的近似处理。2. 边界条件的概化边界条件概化的失误是导致解析解失真的主要原因之一。由于理想化要求常与实际条件相差甚远，成为解析法应用中的难点，也是解析法预测矿坑涌水量的重要环节。（1）侧向边界的概化lPq34qPl+ml Pl2Pd2P+S1T壶天河图 14-2恩口矿区边界条件转化示意图1T1s+P2d 下叠大冶组； 2P2l 上叠龙潭组隔水层；3P1m+P1q4下二叠茅口组与栖霞组岩溶含水层；b号矿体矿坑名师资料总结 - - -精品资料欢迎下载 - - -
31、- - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 27 页 - - - - - - - - - 解析法要求将复杂的边界补给条件概化为隔水与供水两种进水类型；同时，将不规则的边界形态， 简化为规则的。 但实际问题中一般难以具有上述理想条件，其进水条件常常既不完全隔水，又不具有无限补给能力，它的分布也极不规则。为此，必须通过合理的概化，缩小理论与实际的差距，满足近似的计算要求。其要点是：立足于整体概化效果；以均衡为基础，用好等效原则。等效原则是边界概化中的无奈之举，即：通过对概念（如相对隔水边界、近似定水边界） 寻找近似处理的途径
32、；或根据等效原则将垂向越流补给和侧向补给共同构造定水头边界，将局部进水口概化为区域进水边界等。但这些等效原则的应用，必须建立在区域水均衡条件论证的基础上，并涉及参数的优化处理。充分考虑开采因素。疏干流场始终处于补给量与疏干量不断变化的动平衡状态，随着开采条件的变化，边界的位置及其进水条件常发生转化，如湖南恩口煤矿的东部边界（见图14-2） ，在水平疏干时东部壶天河不起作用；开采延伸至水平时，因排水量增大漏斗扩展到壶天河， 成为茅口灰岩的定水头供水边界；当疏干达到水平时，排水量随降深继续增加，当壶天河的补给能力无法与其平衡时，其定水头供水边界已不复存在，漏斗扩展至由隔水层构造的隔水边界，但壶天河
33、仍以变水头集中补给形式平衡疏干漏斗的发展。概化时， 应与西部边界的供水条件作统一的整体殾考虑，如仅就东部边界而言，可用等效原则按第一类越流边界处理， 但须从均衡出发，确定一个相当于第一类越流作用的“引用越流系数” 取代；此外，也可单独计算壶天河的渗漏量，作为矿坑涌水量的一部分。边界几何形态的概化也需认真对待。如湖北铜汞山铜矿的露天矿涌水量预测：矿坑充水来自围岩大理岩，与东西两侧岩浆岩隔水层呈似以300交角，向南敞开（见图14-3）所。60 年代勘探时，概化为东侧直线隔水的环状供水边界，采用非完整井稳定井流域，预测矿坑涌水量为59587985m3/d，而实际涌水量仅3790 m3/d，误差 57
34、111%。70 年代回访调查验证计算时，采用300扇形补给边界的稳定流近似计算，得涌水量3685 m3/d，周期实际涌水量为 3416 m3/d，误差仅 7.8%。证明边界形态概化的重要性。边界概化应把重点放在主要供水边界上。孙纳正教授运用数值法，与解析法对理想化模型的对比验证计算表明：简化供水边界的形状往往会带来较大的误差，但简化隔水边界的形状影响一般不大。（ 2）各种类型侧向边界条件下的计算方法 映射法。即根据地下水动力学中的映射迭加原理，获得矿坑涌水量预测的描述各种特定边界条件下的解析公式。可采用如下一般形式表示：稳定流RQcrR/)(2非稳定流rRKUQ/4式中，R与rR分别为稳定流与
35、非稳定流的边界类型条件系数。详见地下水动力学及有关文献。分区法。也称卡明斯基辐射流法。它是从研究稳定状态下的流网入手，根据疏干流场的边界条件与含水层的AAllh2h1AA5575图 14-4某矿区辐射流计算图名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 27 页 - - - - - - - - - 非均质性特点，沿流面和等水压面将其分割为若干条件不同的扇形分流区(见图 14-4)，每个扇形分流区内其地下水流都呈辐射流，其沿流面分割所得的各扇形区边界为阻水边界，而沿等水
36、压面分割所得的扇形区边界为等水头边界。常用卡明斯基平面辐射流公式分别计算各扇形区的涌水量iQ潜水LhhbbbbQi2lnln)(22212121?承压水LhhMbbbbQi)(lnln)(212121?式中，1b与2b为分流区辐射状水流上下游断面的宽度；1h与2h为1b与2b断面隔水底板上的水头高度；L为1b与2b断面的间距。然后，按下式求各分区流量的总和：nniiQQQQQ211每个扇形区内的下游断面，是以直接靠近井巷的硫干漏斗等水头线的一部分为准；而上游断面则以远离井巷的供水边界上等水头线面一部分为准。（3）垂向越流补给边界类型的确定及其计算当疏干含水层的顶底板为弱透水层时，其垂向相邻含水
37、层就会通过弱透水层对疏干层产生越流补给， 出现所谓的越流补给边界。越流补给边界分定水头和变水头二类，解析法对后者的研究还尚待解决。产生定水头垂向越流补给的矿坑涌水量计算，可用增加越流参数项B的形式来表示：稳定流)(/200BrKkMSQ非稳定流)(/40BruWTSQ式中：B越流参数KTM；K垂向弱透水层渗透系数；M垂向弱透水层厚度；0K零阶二类修正贝塞尔函数。3. 参数确定（1）渗透系数（ K）值渗透系数是解析公式中的主要参数。我国矿山大多为分布于基岩山区的裂隙、岩溶充水矿床， 充水含水层的渗透性具明显不均匀性，根据解析计算要求，应作均值概化， 同时这也是保证渗透系数具有代表的措施之一。矿坑
38、涌水量预测中常用的方法有两种。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 27 页 - - - - - - - - - 加权平均值法：又可分为厚度平均、面积平均、方向平均法等。如厚度平均，则公式为：niiiniiiiCPHMKHMK11)()(式中：)(iiHM承压（潜水）含水层各垂向分段厚度；iK相应分段的渗透系数。流场分析法：有等水位线图时，可采用闭合等值线法：hLLMrQKcpCP)(221或据流场特征，采用分区法：niCPLhhbbbbQK122212121
39、2lnln)(式中：1L、2L任意两条（上、下游）闭合等水位线的长度；r两条闭合等水位线的平均距离；h两条闭合等水位线的平均距离；CPM含水层的平均厚度；Q涌水量；1b、2b辐射状水流上、下游断面上的宽度；1h、2h1b和2b断面隔水底板以上的水头高度；L1b和2b断面之间的距离。（2）大井引用半径0r的确定矿坑的形状极不规则，龙其是坑道(井巷 )系统，分布范围大，形状千变万化，构成了复杂的内边界。根据解析法计算模型的特点，要求将它理想化。经观测，坑道系统排水时，其周边逐渐形成了一个统一的降落漏斗。因此，在理论上可将形状复杂的坑道系统看成是一个理想“大井”在工作，此时整个坑道面积，看成是相当于
40、该“大井”的面积。整个坑道系统的涌水量， 就相当于 “大井”的涌水量，这样就使一般的井流公式能适应于坑道系统的涌水量计算。这种方法，在矿坑涌水量预测中称为“大井法”。 “大井”的引用半径0r，在一般情名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 12 页，共 27 页 - - - - - - - - - 况下用下式表示：FFr565.00式中，F为坑道系统分布范围所圈定的面积。确切地说，近似等于为保证并田设计生产率所必需的坑道所圈定面积的大小，或者以降落漏斗距坑道最近处的封闭等水位
41、线所围起来的面积。如果开采面积近于圆形、方形时，采用上式较准确，对于形状特别的面积，可采用其它专门公式计算。（3）引用半径0R的确定对于在实际问题中是否存在裘布依关于园形定水头的假定条件，以及齐姆模型中影响半径概念是否合理，学术界存在众多争议，特别是非稳定流理论与越流理论的产生，人们认识到某种稳定状态的出现可能是越流作用所造成，这样影响半径的概念离开了原有的含义，于是引出了所谓“引用影响半径”。从稳定井流理论的实际应用出发，上述概念是可以被接受的，即根据等效原则，将疏干量与补给量相平衡时出现的稳定流场，其边界用一个引用的园形等效外边界进行慨化，其与“大井”中心的水平距离称为引用影响半径，也称为
42、补给半径。即00rRR。同理，在用割面平面流解析公式计算使表水中坑道涌水量，也就有了引用影响带宽度（0L）的概念，即疏干坑道中心与外边界之间的距离。在稳定流条件下，引用影响半径0R为一个常量，也称补给半径；在非稳定流条件下则是一个不断变化着的变量，这样在理论上解决了稳定井流理论及其引用影响半径计算公式的实用问题。矿山疏干实际表明，矿坑排水的影响范围，总是随时间的延长、排水量的增加以及坑道的推进而不断扩大，直到天然边界为止，它不可能被限制在一个不是边界的理想“半径”之内。此外， 对比计算表明， 若确定影响半径的误差为23 倍，则矿坑涌水量的计算误差可达30 60；此外，若取偏低值其误差远比取偏高
43、值要大。因此，矿坑涌水量预测时，能否用解析公式及常见的经验公式来近似的确定影响半径值得探讨；对开拓井巷的捅水量预测，最好采用抽水试验外推法，即根据多落程的抽水试验，确定降深与影响半径或流量与影响半径的线性关系，外推某琉干水位或某疏于量的相应疏干半径值。如：maSR1；或maQR1对坑道系统的涌水量预测，应根据疏干中心天然水文地质边界线之间距离的加权平均值计算，即塞罗瓦特科公式：LLbrRCPCP0式中：0r“大井”的引用半径；cpb井巷轮廓线与各不同类型水文地质边界间的平均距离；名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名
44、师精心整理 - - - - - - - 第 13 页，共 27 页 - - - - - - - - - L各类型水文地质边界线的宽度。（4）最大疏干水位降深maxS的确定在理论上， 目前解析解还无法处理承压区与无压区同时并存与大降深的潜水问题，对于实际问题，则是矿床硫干时最大可能水位降是多少，如何近似确定最大硫于水位降深m axS值。爱尔别尔格尔在实验中取得的潜水最大水位降深等于潜水含水层一半的结论，即HS2/1max(扩大应用到承压含水层时，HHS2/1max)，一直是水文地质计算中所遵循的概念。近年来，我国通过渗流槽及野外抽水试验证明这一结论是保守的。maxS可以超过H8.0，在矿坑涌水量
45、计算中，通常不考虑这一概念。据观测，在长期硫于条件下的大截面井巷系统外缘，动水位(h)一般不超过12m，它所引起的涌水量计算偏大值一般为0.51。因此，矿坑涌水量预测时，最大疏干水位降一般取HSmax。另一个理论问题，即最大水位降HSmax时的最大涌水量计算问题。众所周知，当HSmax时，裘布依公式在理论上就会“失真”，这正是稳定井流理论的最大缺陷之一；而泰斯公式则是从承压水含水层建立起来的，扩大到无压含水层使用时，(作最大降深疏干时，承压含水层均转化为无压水层)，常把隨时间变化的含水层厚度作线性处理，即取不变的平均值，这种线性化处理必然带来误差，据研究当降深超过含水层厚度30%时，非稳定井流
46、公式要偏实际情况，出现明显误差，更不用说是作最大水位降的计算了。如上所述，不难看出矿坑涌水量预测时，作最大水位降的最大疏干量计算，对解析法来说不是很合科适宜的。（三）实例1 最佳疏干量（最Q）的确定某铁矿地处灰岩区，裂隙岩溶发育较均匀，地下水运动符合达西定律，矿区内有部分地下水动态长期观测资料，其它地质条件略。（1）要求当疏干水平(或中段 )的水位降深 (S)确定后，则疏干量(Q)是时间 (t)的函数。这样，疏干量Q就是与疏干时间t有关的一组数据。 某水平的正常疏干量，应是该水平预测的矿坑涌水量值。设计部门要在一组具不同疏干强度Q及与其相应的时间t的对比中，选出最佳疏干方案，即选择排水能力要求
47、不过大，而疏干时间又不长的方案。疏干时间通常要求控制在两个雨季之间，否则Q的计算则无意义。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 14 页，共 27 页 - - - - - - - - - （2）任务给定的条件是： 疏干中段水位降（S）确定为零米标高；疏干时间要求在两个雨季间完成。（3）最佳疏干量的计算与分析第一步：初选疏干时间段t。根据第二项任务，在现有地下水动态曲线（图）上初选三个时间段，即270 d、210 d、150d，供计算分析。第二步： 确定相应的Q值。根据给定的
48、零米标高，从动态曲线图上确定出各霎时间段相对应的S值，即：t(d) 3t270 2t210 1t150 S(m) 3S100 2S90 1S80 第三步：求相应的Q值，利用公式（符号为常用地下水动力学符号）：)4(/42TtrWTSQ在已知1t、1S、2t、2S、3t、3S的条件下，求得相应的1Q、2Q、3Q，作为第四步分析的初值。第四步： 绘制不同疏干强度Q条件下的)(tfS曲线。在初值1Q、2Q、3Q的范围内，通过内插给出一组供进一步分析的疏干量数据。其公式为：)4(42TtrWTQS分析不同疏干量时的S随t的变化规律，并绘制不同疏干量条件下的)(tfS曲线。（图 6）第五步：绘制不同定降
49、深S条件下的图 5 某矿区地下水水位动态曲线图721210t (d)150210270(m)h100806040200图 5 某矿区地下水水位动态曲线图图 6 不同疏干量条件下曲线S=f （ t）(m3/d)Q5=40000Q4=36000Q3=30000Q2=23000Q1=15000t (d)EDCBAdcba（m）S1009080名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 15 页，共 27 页 - - - - - - - - - )(tfQ曲线。根据图作出不同降深S条件下
50、的疏干量Q与时间t的关系曲线)(tfQ(图 7)，进行不同S条件下，疏干量Q与疏干时间的对比分析。第六步：绘制降深S与最佳疏干量Q的关系曲线。根据图中各)(tfS曲线的拐点，求出不同降深Q条件厂的最佳疏干强度Q，即拟稳定疏干量与降深的关系曲线(图 148)。第七步：确定最佳硫干量，并检验其可行性。根据图148 取得的不同降深的最佳疏干量最Q检验它们达到S时所需的时间t，是否满足任务要求，即是否能在两个雨季之间完成疏干任务。如符合需要，预测就算完成；不符合，则还要重复进行，直至所选取的最佳疏干量满足任务要求的S与t时为止。从图 148 取S90m，则佳Q为 20000m3d 从图中求得t=200