专门水文地质学 (71).doc

吉林大学精品课专门水文地质学教材专门水文地质学 10.4矿坑涌水量预测 一、矿坑涌水量预测的内容、方法、步骤与特点 （一）矿井涌水量预测的内容及要求 矿坑涌水量预测是一项重要而复杂的工作，是矿床水文地质勘探的重要组成部分。 矿坑涌水量是指矿山开拓与开采过程中，单位时间内涌入矿坑包括井、巷和开采系统的水量。通常以m3/h表示。它是确定矿床水文地质条件复杂程度的重要指标之一，关系到矿山的生产条件与成本，对矿床的经济技术评价有很大的影响。并且也是设计与开采部门选择开采方案、开采方法，制定防治水疏干措施，设计水仓、排水系统与设备的主要依据。因此，在矿床水文地质调查中，要求正确评价未来矿山开发各个阶段的涌水量。其内容与要求包括可概括为以下四个方面 （1）矿坑正常涌水量指开采系统达到某一标高水平或中段时，正常状态下保持相对稳定的总涌水量，通常是指平水年的涌水量。 （2）矿坑最大涌水量是指正常状态下开采系统在丰水年雨季时的最大涌水量。对某些受暴雨强度直接控制的裸露型、暗河型岩溶充水矿床来说，常常还应依据矿山的服务年限与当地气象变化周期，按当地气象站所记录的最大暴雨强度，预测数十年一遇特大暴雨强度产生时，可能出现暂短的特大矿坑涌水量，作为制订各种应变措施的依据。 （3）开拓井巷涌水量指包括井筒立井、斜井和巷道（平、平巷、斜巷、石门在开拓过程中的涌水量。 （4）疏干工程的排水量是指在规定的疏于时间内，将一定范围内的水位降到某一规定标高时，所需的疏干排水强度。 对于地质勘探阶段来说，主要是进行评价性的计算，以预测正常状态下矿坑涌水量及最大涌水量为主。至于开拓井巷的涌水量预测和专门性疏干工程的排水量的计算，由于与矿山的生产条件密切相关，一般均由矿山基建部门或生产部门承担。 （二）矿坑涌水量预测的方法 根据当前矿床水文地质计算中常用的各种数学模型的地质背景特征极其对水文地质模型概化的要求，可作如下类型的划分 （三）矿坑涌水量预测的步骤 矿坑涌水量预测是在查明矿床的充水因素及水文地质条件的基础上进行的。它是一项贯穿矿区水文地质勘探全过程的工作，一个正确预测方案的建立，是随着对水文地质条件认识的不断深化，不断修正、完善而逐渐形成的，一般应遵循如下三个基本步骤 1. 选择计算方法与相应的数学模型 详勘阶段均要求选择2个或2个以上的计算方法，以相互检验，映证。选择时必须考虑三个基本要素 （1）矿床的充水因素及水文地质条件复杂程度。如位于当地侵蚀基准面之上，以降水入渗补给的矿床，应采用水均衡法；水文地质条件简单或中等的矿床，可采用解析法或比拟法；水文地质条件复杂的大水矿床，要求采用数值方法； （2）勘探阶段对矿坑涌水量预测的精度要求； （3）勘探方法、勘探工程的控制程度与信息量如水均衡法，要求不少于一个水文年的完整均衡域的补给与排泄项的动态资料；Qs曲线方程外推法，要求具抽水试验的水位降达到预测标高水柱高度的1/21/3；解析法，要求勘探工程全面控制含水层的非均质各向异性、非等厚的结构特征及其边界条件与补给、径流与排泄，并提供数值模型的建立、识别、预测所需的完整信息数据，这些数据的获取，只有采用大型抽、放水试验对渗透场进行整体控制与揭露才可能做到。 因此，计算方法与相应数学模型类型的选择，与矿床的充水因素及水文地质条件复杂程度、勘探方法勘探工程的控制程度及信息量是相互关联的，统一在最佳技术经济条件这一原则下，所以数学模型类型选择是否合理，可以用以下标准衡量 一是对矿床水文地质条件的适应性指能否正确刻画水文地质条件的基本特征； 二是对勘探方法勘探工程控制程度的适应性指是否最充分的利用勘探工程提供的各种信息，即信息的利用率；同时，也可理解为所选数学模型要求的勘探信息是否有保证，即信息的保障率。 2. 构造水文地质模型 矿坑涌水量预测中数学模型的作用，是对水文地质条件进行量化，因此预测精度主要取决于对充水因素与水文地质条件判断的准确性，由于不同数学模型类型对水文地质条件的刻画形式与功能各异，因此必须按数学模型的特点构造水文地质模型，称水文地质条件概化。概化后的水文地质模型称水文地质概念模型，它在地质实体与数学模型之间起中介桥梁作用。下面以最基本的预测方法解析法与数值法为例作一讨论。 （1） 概化已知状态下的水文地质条件； （2） 给出未来开采状态下的内边界条件； （3） 预测未来开采状态下的外边界条件。 解析法将复杂的含水层结构与内外边界，以理想化模式构造理论公式，因此必须按解析解要求进行概化。如含水层均质等厚，内外边界几何形态规则，边界供水条件简单、确定。 数值法以近似分割原理对复杂的含水层结构、内外边界条件进行量化“逼真”，概化时要求以控制水文地质条件与内外边界的节点参数、水位与流量来构造水文地质概念模型。 随着数学模型研究的不断进展，现代水文地质计算对水文地质模型的要求越来越高。目前，对复杂的大水矿床来说，一个可靠的水文地质模型的建立，必须贯穿整个勘探过程，并大致经历三个阶段。即 第一阶段通过对以往资料的整理，提出水文地质模型的“雏型”，作为下一步勘探设计的依据。尤其对大型抽（放）水试验的设计具有特殊重要的意义； 第二阶段根据进一步勘探提供的各种信息数据。特别是大型抽放水资料，通过流场分析或数值模拟，完成对“雏型”模型的调整，建立水文地质模型的“校正型”； 第三阶段在“校正型”的基础上，按开采方案给出疏干工程的内边界条件，根据勘探资料预测不同疏干条件下的外边界条件，建立水文地质概化模型的“预测型”。 3. 计算数学模型，评价预测结果 应该指出，不能把数学模型的解 仅仅看作是一个单纯的数学计算，而应看作是对水文地质模型和数学模型进行全面验证识别过程，也是对矿区水文地质条件从定性到定量在回到定性的不断深化的认识过程。 （四）矿坑涌水量预测的特点 虽然矿坑涌水量预测的原理方法与供水水资源评价类同，但其预测条件、预测要求与思路各有不同。如 （1）供水水资源评价，以持续稳定开采确保枯水期安全开采量为目标；而矿坑涌水量预测则以疏干丰水期的最大涌水量为目标； （2）矿床大多分布于基岩山区。含水质的非均质性突出。参数代表性不易控制；边界条件复杂、非确定性因素多，常出现紊流、非连续流与管道流。定量化难度大。 （3）矿山井巷类型及其分布千变万化，开采方法、开采速度与规模等生产条件复杂且不稳定，与供水的取水建筑物简单、分布有序、生产稳定形成显明对比，给矿坑涌水量预测带来诸多不确定性因素。 （4）矿坑涌水量预测多为大降深下推。此时开采条件对水文地质条件的改变难以预料和量化，这与供水小降深开采有明显差异。 （5）矿床水文地质勘探从属于矿产地质勘探，与专门性的供水水文地质勘探对比，前者一般投入小、工程控制程度低，预测所需的信息量相对少而不完整。 以上特点，决定了矿坑涌水量预测中存在诸多产生误差的客观条件。因此属于评价性计算，为矿山设计及采前进一步专门性补充勘探提供依据。 二、QS曲线方程外推法 （一）原理与应用条件 指用稳定井流条件下抽水试验的Qf（s）方程，外推未来疏干水位降的涌水量。实质上也是一种相似条件下的比拟法。应用时的前提条件是一、抽水试验建立Qf（s） ，应符合稳定井流条件；二、抽水试验的各种条件应与预测对象的疏干条件接近。因此，必须重视试验的技术条件，包括 1.应将抽水试验孔布置在预测对象的分布地段，保证水文地质条件的一致性； 2.采用大口径（或孔组）试验，计算时为消除井径对涌水量的影响，需做井径换算； 3.抽水降深应大于疏干水位水柱高度的1/21/3，计算时的外推疏干降深不应超过1.75倍的抽水降深，主要考虑疏干状态下的补给条件； 4.用枯季抽水试验预测正常涌水量，根据雨季试验预测季节性最大涌水量； 5.要排除抽水过程中一切自然和人为随机影响因素的干扰。 Qs曲线方程法的优点是回避各种水文地质参数求参过程中的失真，计算简单易行。适用于建井初期的井筒涌水量预测。上水平疏干资料外推下水平的涌水量，以及矿床规模小、矿体分布集中、边界条件和含水结构复杂的涌水量预测。 （二）计算方法与步骤 1. 鉴别Qs曲线类型 （1）曲度法即用曲度n值进行鉴别，当n＝1时，为直线QqS；1＜n＜2时，为幂曲线Q；n＝2时，为抛物线SaQbQ2；n＞2时，为半对数曲线SablgS。如果n＜1时，表明抽水试验不正确。 2. 确定方程参数 （1）最小二乘法应根据Qf（s）类型选用最小二乘法，如常见的幂函数型 （2）图解法即利用直角坐标的图解，a为图解中纵坐标上所切的截距线段；b为直线对水平倾角的正切。 其它类型详见地下水动力学。 3. 井径换算 （1）根据稳定井流理论中，不同地下水运动（层流、紊流）状态下井径与涌水量关系进行换算，如层流 ；紊流。试验表明井径对涌水量的影响，一般比对数大，而小于平方根。 （2）用二次或二次以上不同孔径的抽水试验资料建立井径d的换算公式Qmdn，m、n值可用最小二乘法求得。 4. 水均衡论证 将预测结果进行疏干条件下的水均衡论证，分析所得的最大涌水量是否符合疏干状态下的边界补给条件。 （三）实例 1. 广东某金属矿，曾用Qs曲线法预测50m水平的涌水量为14450m3／d，与坑道放水试验的外推涌水量14000m3/d接近，而解析法预测的结果12608m3／d则偏小12％。 2. 湖南叶花香铜矿，1970年采用试生产坑道的排水量与水位降资料。分别用两种曲线类型对未来四个中段的涌水量进行外推预测，后经实际检验见表14-2。 表14-2 计算方法 不同降深的水量 52.45m 73.33m 99.47m 水量（m3/d） 误差 （） 水量（m3/d） 误差 （） 水量（m3/d） 误差 （） 计算 实际 计算 实际 计算 实际 SaQbQ2 15145 16937 -10.60 18223 23644 -22.90 21514 31155 -31.00 Q 14840 -12.40 18093 -23.50 21582 -30.70 实际检验表明预测结果与实际情况较接近，产生误差的主要原因是其间曾发生两次大突水淹井，大量泥沙溃入矿坑，地面出现岩溶坍陷和裂缝所致。 三、相关分析法 （一）原理与应用条件 相关分析是根据涌水量与主要影响因素之间相关关系的密切程度建立回归方程，利用抽水试验或开采初期的疏干资料，预测矿坑涌水量或外推开采后期下水平的涌水量。其原理已在供水中介绍。根据实际资料的统计，多元复相关预测远比单相关效果好，其回归方程表达的内容丰富，可反映除降深外的各种影响因素。它的应用条件与Qs关系方程类同，但对原始数据的采集有严格要求 1代表性（规范）要求不少于一个水文年（包括丰、平、枯季节）的动态观测数据，同时数据（择本）量不少于30个； 2一致性指应与预测对象上条件相一致； 3独立性与相关性即多自变量有独立的变化规律，相互间关系不大；而与涌水量之间均存在密切的相关关系，（规范）要求相关系数不低于0.7。 （二）实例与计算方法 1. 利用勘探阶段抽水试验资料预测矿坑涌水量 如广东沙洋矿通过在勘探阶段设计相距6m的两个抽水孔和十余个不同距离的观测孔组成的群孔抽水试验，取得了复相关计算所需的涌水量Q与井径r（是将距抽水孔不同距离观测孔的位置概化为疏干状态下的坑道系统不同面积的作用半径）、水位降S（即不同作用半径的水位降，以模拟疏干水位降）有关资料，（见表3）通过求参建立了复相关幂函数预测方程 其复相系数达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 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-4。 表14-4 相关系数r lgQ0 lgP0 lgH0 lgF0 lgT0 lgQ0 0.664 0.451 0.593 0.175 lgP0 0.664 0.340 0.680 0.323 lgH0 0.451 0.340 0.559 0.523 lgF0 0.593 0.680 0.559 0.778 lgT0 0.175 0.323 0.523 0.778 根据判别得知，生产时间T0对的影响不大（相关系数为0.175）。用多元复相关计算，求得四元复相关曲线回归方程 其复相关系数。在此基础上建立了比拟公式 预测结果与传统的单位涌水量法相比，使误差减少1.4倍（式中P、H、F为设计值）。 四、解析法 （一）解析法的应用条件 解析法是根据解析解的建模要求，通过对实际问题的合理概化，构造理想化模式的解析公式，用于矿坑涌水量预测。具有对井巷类型适应能力强、快速、简便、经济等优点，是最常用的基本方法。解析法预测矿坑涌水量时，以井流理论和用等效原则构造的“大井”为主，后者指将各种形态的井巷与坑道系统，以具有等效性的“大井”表示，称“大井”法。因此说矿坑涌水量计算的最大特点是“大井法”与等效原则的应用，而供水则以干扰井的计算为主。 稳定井流解析法应用于矿坑疏干流场处于相对稳定状态的流量预测。包括①在已知某开采水平最大水位降条件下的矿坑总涌水量；②在给定某开采水平疏干排水能力的前提下，计算地下水位降深（或压力疏降）值。 非稳定解析法用于矿床疏干过程中地下水位不断下降，疏干漏斗持续不断扩展，非稳定状态下的涌水量预测。包括①已知开采水平水位降（s）、疏干时间（t），求涌水量（Q）；②已知Q、s，求疏干某水平或漏斗扩展到某处的时间（t）；③已知Q、t，求s，以确定漏斗发展的速度和漏斗范围内各点水头函数隨时间的变化规律，用于规划各项开采措施。在勘探阶段，以选择疏干量和计算量最大涌水量为主。 （二）计算方法 如上所述，应用解析法预测矿坑涌水量时，关键问题是如何在查清水文地质条件的前提下，将复杂的实际问题概化。它可概括为如下三个重要方面分析疏干流场的水力特征，合理概化边界条件，正确确定各项参数。 1. 分析疏干流场的水力特征 矿区的疏干流场是在天然背景条件下，迭加开采因素演变而成。分析时，应以天然状态为基础，结合开采条件作出合理概化。 （1）区分稳定流与非稳定流 矿山基建阶段，疏干流场的内外边界有受开拓井巷的扩展所控制，以消耗含水层储量为主，属非稳定流；进入回采阶段后，井巷输廊大体已定，疏干流场主要受外边界的补给条件控制，当存在定水头（侧向或越流）补给条件时，矿坑水量被侧向补给量或越流量所平衡，流场特征除受气候的季节变化影响外，呈现对稳定状态。基本符合稳定的“建模”条件，或可以认为两者具等效性；反之，均属非稳定流范畴。 如河北开滦煤矿，其矿坑涌水量随坑道走向的延展而增加，但这种延展暂停时，涌水量立即出现相对的稳定。此时仅表现有受降水的季节变化在一定变幅范围内上下波动，并出现强出水点和边缘出水点成袭夺中间出点现象，而总涌水量不变。又如辽宁复州粘土矿，其涌水量随采深增加，但某一水平进入回采时，其涌水量就逐渐稳定，并保持到下一水平突水进止，在此阶段虽然也出现下水平突水点袭夺上水平突水点现象，但总涌水量却保持不变。由此可见，在某些矿区的疏干过程中，不仅存在疏干流场的相对稳定阶段，而且隨矿山工程的进展而不断相互转化。 但选用稳定流解析法时要慎重，必须进行均衡论证，判断疏干区是否真正存在定水头供水边界或定水头的越流系统。此外，对于地下水储存量较大的矿区，要单独计算疏干过程中不同阶段含水层储存量的 放强度，大量生产实际表明，矿坑最大涌水量，并非出现在疏干过程的稳定阶段。 （2）区分达西流与非达西流 在矿坑涌水量计算时，常遇到非达西流问题，它涉及解析法的应用条件，在宏观上可概括有两种情况 一是暗河管道岩溶充水矿床，地下水运动为压力管道流与明渠流；此外，分水岭地段的充水矿床，矿坑涌水量直接受垂向入渗降雨强度控制，与水位降深无关。两者均与解析法的“建模”条件相距甚大，矿坑涌水量预测应选择水均衡法或各种隨机统计方法。 二是，局部状态的非达西流，常发生在大降深疏干井巷附近与某些特殊构造部位，它只对参数计算与参数的代表性产生影响。在宏观上，它是一个流态概化问题，不存在解析法的应用条件问题。 （3）区分平面流与空间流 严格讲，在大降深疏干条件下，地下水运动的垂向速度分量不能忽略，均为三维空间流（包括非完整井巷的地下水运动）问题，其分布范围仅限于井巷附近，均为含水层厚度的1.54.75倍。因此，在矿坑涌水量预测中，大多将其纳入二维平面流范畴，在宏观上不影响预测精度。计算时应根据井巷类型作出不同的概化 如竖井的涌水量计算，可概化为平面径向流问题，以进流公式表达；计算水平巷道涌水时，以剖面平面流近似，采用单宽流量解析公式，但其两端上往往也产生辐射流（见图），需要考虑它的存在，并采用平面径向流公式补充计算巷道端部的进水口。 坑道系统则复杂得多，根据“大井法”原理，一般以近似的径向流概化，但当坑道系统近于带状的狭长条形时，也可概化为剖面流问题。 对于倾斜坑道，根据阿勃拉莫夫有关水电比拟法的研究，证明坑道的倾斜对涌水量影响不大，可根据坑道的倾斜度，分别按竖井或水平巷道进行近似。即若坑道倾斜度＞450时，视其与竖井近似，用井流公式计算；若坑道倾斜度＜450时，则视其与水平巷道近似，用单宽流量公式计算。 根据解析解的存在条件，一些简单的非完整井巷涌水量计算，可以运用三维空间问题予以解决。此时，可根据非完整井的特点，运用地下水动力学中映射法与分段法的原理来求解。通常用平面分段法解决完整竖井的涌水量计算，用剖面分段法解决非完整平巷的涌水量计算。 （4）区分潜水与承压水 与供水不同，在降压疏干时，往往出现承压水转化为潜水或承压－无压水。此外，在陡倾斜含水层分布的矿区，还可能出现坑道一侧保持原始承压水状态，而另一侧却由承压水转化无压水或承压－无压水的现象。概化时，需从宏观角度作等效的近似处理。 2. 边界条件的概化 图14-2 恩口矿区边界条件转化示意图 1－T1sP2d 下叠大冶组；2－P2l上叠龙潭组隔水层； 3－P1mP1q4 下二叠茅口组与栖霞组岩溶含水层； 4－P1q3 下二叠栖霞组李子塘段隔水层；5－Ⅰ，Ⅱ， Ⅲ，疏干水平。 图14-3 铜录山矿区边界概化图 1－Mb 大理岩； 2－γδ岩浆岩。 边界条件概化的失误是导致解析解失真的主要原因之一。由于理想化要求常与实际条件相差甚远，成为解析法应用中的难点，也是解析法预测矿坑涌水量的重要环节。 （1）侧向边界的概化 解析法要求将复杂的边界补给条件概化为隔水与供水两种进水类型；同时，将不规则的边界形态，简化为规则的。但实际问题中一般难以具有上述理想条件，其进水条件 常常既不完全隔水，又不具有无限补给能力，它的分布也极不规则。为此，必须通过合理的概化，缩小理论与实际的差距，满足近似的计算要求。其要点是 ①立足于整体概化效果； ②以均衡为基础，用好等效原则。等效原则是边界概化中的无奈之举，即通过对概念（如相对隔水边界、近似定水边界）寻找近似处理的途径；或根据等效原则将垂向越流补给和侧向补给共同构造定水头边界，将局部进水口概化为区域进水边界等。但这些等效原则的应用，必须建立在区域水均衡条件论证的基础上，并涉及参数的优化处理。 ③充分考虑开采因素。疏干流场始终处于补给量与疏干量不断变化的动平衡状态，随着开采条件的变化，边界的位置及其进水条件常发生转化，如湖南恩口煤矿的东部边界（见图2），在Ⅰ水平疏干时东部壶天河不起作用；开采延伸至Ⅱ水平时，因排水量增大漏斗扩展到壶天河，成为茅口灰岩的定水头供水边界；当疏干达到Ⅲ水平时，排水量随降深继续增加，当壶天河的补给能力无法与其平衡时，其定水头供水边界已不复存在，漏斗扩展至由隔水层构造的隔水边界，但壶天河仍以变水头集中补给形式平衡疏干漏斗的发展。概化时，应与西部边界的供水条件作统一的整体殾考虑，如仅就东部边界而言，可用等效原则按第一类越流边界处理，但须从均衡出发，确定一个相当于第一类越流作用的“引用越流系数”取代；此外，也可单独计算壶天河的渗漏量，作为矿坑涌水量的一部分。 ④边界几何形态的概化也需认真对待。如湖北铜汞山铜矿的露天矿涌水量预测矿坑充水来自围岩大理岩，与东西两侧岩浆岩隔水层呈似以300交角，向南敞开（见图3）所。60年代勘探时，概化为东侧直线隔水的环状供水边界，采用非完整井稳定井流域，预测矿坑涌水量为59587985m3/d，而实际涌水量仅3790 m3/d，误差57111。70年代回访调查验证计算时，采用300扇形补给边界的稳定流近似计算，得涌水量3685 m3/d，周期实际涌水量为3416 m3/d，误差仅7.8。证明边界形态概化的重要性。 ⑤边界概化应把重点放在主要供水边界上。孙纳正教授运用数值法，与解析法对理想化模型的对比验证计算表明简化供水边界的形状往往会带来较大的误差，但简化隔水边界的形状影响一般不大。 （2）各种类型侧向边界条件下的计算方法 ① 映射法。即根据地下水动力学中的映射迭加原理，获得矿坑涌水量预测的描述各种特定边界条件下的解析公式。可采用如下一般形式表示 稳定流 非稳定流 式中，与分别为稳定流与非稳定流的边界类型条件系数。详见地下水动力学及有关文献。 图14-4 某矿区辐射流计算图 ②分区法。也称卡明斯基辐射流法。它是从研究稳定状态下的流网入手，根据疏干流场的边界条件与含水层的非均质性特点，沿流面和等水压面将其分割为若干条件不同的扇形分流区见图14-4，每个扇形分流区内其地下水流都呈辐射流，其沿流面分割所得的各扇形区边界为阻水边界，而沿等水压面分割所得的扇形区边界为等水头边界。常用卡明斯基平面辐射流公式分别计算各扇形区的涌水量 潜水 承压水 式中，与为分流区辐射状水流上下游断面的宽度；与为与断面隔水底板上的水头高度；为与断面的间距。 然后，按下式求各分区流量的总和 每个扇形区内的下游断面，是以直接靠近井巷的硫干漏斗等水头线的一部分为准；而上游断面则以远离井巷的供水边界上等水头线面一部分为准。 （3）垂向越流补给边界类型的确定及其计算 当疏干含水层的顶底板为弱透水层时，其垂向相邻含水层就会通过弱透水层对疏干层 产生越流补给，出现所谓的越流补给边界。越流补给边界分定水头和变水头二类，解析法对 后者的研究还尚待解决。 产生定水头垂向越流补给的矿坑涌水量计算，可用增加越流参数项的形式来表示 稳定流 非稳定流 式中越流参数；垂向弱透水层渗透系数；垂向弱透水层厚度；零阶二类修正贝塞尔函数。 3. 参数确定 （1）渗透系数（K）值 渗透系数是解析公式中的主要参数。我国矿山大多为分布于基岩山区的裂隙、岩溶充水矿床，充水含水层的渗透性具明显不均匀性，根据解析计算要求，应作均值概化，同时这也是保证渗透系数具有代表的措施之一。矿坑涌水量预测中常用的方法有两种。 ① 加权平均值法又可分为厚度平均、面积平均、方向平均法等。如厚度平均，则公式为 式中－承压（潜水）含水层各垂向分段厚度； －相应分段的渗透系数。 ②流场分析法有等水位线图时，可采用闭合等值线法 或据流场特征，采用分区法 式中、－任意两条（上、下游）闭合等水位线的长度； －两条闭合等水位线的平均距离； －两条闭合等水位线的平均距离； －含水层的平均厚度； －涌水量； 、－辐射状水流上、下游断面上的宽度； 、－和断面隔水底板以上的水头高度； －和断面之间的距离。 （2）大井引用半径的确定 矿坑的形状极不规则，龙其是坑道井巷系统，分布范围大，形状千变万化，构成了复 杂的内边界。根据解析法计算模型的特点，要求将它理想化。经观测，坑道系统排水时，其周边逐渐形成了一个统一的降落漏斗。因此，在理论上可将形状复杂的坑道系统看成是一个理想“大井”在工作，此时整个坑道面积，看成是相当于该“大井”的面积。整个坑道系统的涌水量，就相当于“大井”的涌水量，这样就使一般的井流公式能适应于坑道系统的涌水量计算。这种方法，在矿坑涌水量预测中称为“大井法”。“大井”的引用半径，在一般情况下用下式表示 式中，为坑道系统分布范围所圈定的面积。确切地说，近似等于为保证并田设计生产率所必需的坑道所圈定面积的大小，或者以降落漏斗距坑道最近处的封闭等水位线所围起 来的面积。如果开采面积近于圆形、方形时，采用上式较准确，对于形状特别的面积，可采用其它专门公式计算。 （3）引用半径的确定 对于在实际问题中是否存在裘布依关于园形定水头的假定条件，以及齐姆模型中影响半径概念是否合理，学术界存在众多争议，特别是非稳定流理论与越流理论的产生，人们认识到某种稳定状态的出现可能是越流作用所造成，这样影响半径的概念离开了原有的含义，于是引出了所谓“引用影响半径”。从稳定井流理论的实际应用出发，上述概念是可以被接受的，即根据等效原则，将疏干量与补给量相平衡时出现的稳定流场，其边界用一个引用的园形等效外边界进行慨化，其与“大井”中心的水平距离称为引用影响半径，也称为补给半径。即。同理，在用割面平面流解析公式计算使表水中坑道涌水量，也就有了引用影响带宽度（）的概念，即疏干坑道中心与外边界之间的距离。 在稳定流条件下，引用影响半径为一个常量，也称补给半径；在非稳定流条件下则是一个不断变化着的变量，这样在理论上解决了稳定井流理论及其引用影响半径计算公式的 实用问题。 矿山疏干实际表明，矿坑排水的影响范围，总是随时间的延长、排水量的增加以及坑道 的推进而不断扩大，直到天然边界为止，它不可能被限制在一个不是边界的理想“半径”之 内。此外，对比计算表明，若确定影响半径的误差为23倍，则矿坑涌水量的计算误差可达 30％60％；此外，若取偏低值其误差远比取偏高值要大。因此，矿坑涌水量预测时，能否 用解析公式及常见的经验公式来近似的确定影响半径值得探讨；对开拓井巷的捅水量预测， 最好采用抽水试验外推法，即根据多落程的抽水试验，确定降深与影响半径或流量与影响半 径的线性关系，外推某琉干水位或某疏于量的相应疏干半径值。如 ；或 对坑道系统的涌水量预测，应根据疏干中心天然水文地质边界线之间距离的加权平均值计算，即塞罗瓦特科公式 式中“大井”的引用半径； 井巷轮廓线与各不同类型水文地质边界间的平均距离； －各类型水文地质边界线的宽度。 （4）最大疏干水位降深的确定 在理论上，目前解析解还无法处理承压区与无压区同时并存与大降深的潜水问题，对于 实际问题，则是矿床硫干时最大可能水位降是多少，如何近似确定最大硫于水位降深值。 爱尔别尔格尔在实验中取得的潜水最大水位降深等于潜水含水层一半的结论，即扩大应用到承压含水层时，，一直是水文地质计算中所遵循的概念。近年来，我国通过渗流槽及野外抽水试验．证明这一结论是保守的。可以超过，在矿坑涌水量计算中，通常不考虑这一概念。据观测，在长期硫于条件下的大截面井巷系统外缘，动水位一般不超过12m，它所引起的涌水量计算偏大值一般为0.5％1％。因此，矿坑涌水量预测时，最大疏干水位降一般取。 另一个理论问题，即最大水位降时的最大涌水量计算问题。众所周知，当时，裘布依公式在理论上就会“失真”，这正是稳定井流理论的最大缺陷之一；而泰斯公式则是从承压水含水层建立起来的，扩大到无压含水层使用时，作最大降深疏干时，承压含水层均转化为无压水层，常把隨时间变化的含水层厚度作线性处理，即取不变的平均值，这种线性化处理必然带来误差，据研究当降深超过含水层厚度30时，非稳定井流公式要偏实际情况，出现明显误差，更不用说是作最大水位降的计算了。如上所述，不难看出矿坑涌水量预测时，作最大水位降的最大疏干量计算，对解析法来说不是很合科适宜的。 （三）实例 1最佳疏干量（）的确定 某铁矿地处灰岩区，裂隙岩溶发育较均匀，地下水运动符合达西定律，矿区内有部分 地下水动态长期观测资料，其它地质条件略。 （1）要求 ①当疏干水平或中段的水位降深确定后，则疏干量是时间的函数。这样，疏干量就是与疏干时间有关的一组数据。某水平的正常疏干量，应是该水平预测的矿坑涌水量值。设计部门要在一组具不同疏干强度及与其相应的时间的对比中，选出最佳疏干方案，即选择排水能力要求不过大，而疏干时间又不长的方案。 ②疏干时间通常要求控制在两个雨季之间，否则的计算则无意义。 （2）任务 给定的条件是①疏干中段水位降（）确定为零米标高；②疏干时间要求在两个雨季间完成。 （3）最佳疏干量的计算与分析 第一步初选疏干时间段。根据第二项任务，在现有地下水动态曲线（图）上初选三个时间段，即270 d、210 d、150d，供计算分析。 第二步确定相应的值。根据给定的零米标高，从动态曲线图上确定出各霎时间段相对应的值，即 d 270 210 150 m 100 90 80 图5 某矿区地下水水位动态曲线图 第三步求相应的值，利用公式（符号为常用地下水动力学符号） 在已知、、、、、的条件下，求得相应的、、，作为第四步分析的初值。 第四步绘制不同疏干强度条件下的曲线。在初值、、的范围内，通过内插给出一组供进一步分析的疏干量数据。其公式为 图6 不同疏干量条件下Sft曲线 分析不同疏干量时的随的变化规律，并绘制不同疏干量条件下的曲线。（图） 第五步绘制不同定降深条件下的曲线。根据图作出不同降深条件下的疏干量与时间的关系曲线 图，进行不同条件下，疏干量与疏干时间的对比分析。 第六步绘制降深与最佳疏干量的关系曲线。根据图中各曲线的拐点，求出不同降深条件厂的最佳疏干强度，即拟稳定疏干量与降深的关系曲线图148。 第七步确定最佳硫干量，并检验其可行性。根据图148取得的不同降深的最佳 图8 降深与疏干量关系曲线 图7 不同降深条件下Qft曲线 疏干量检验它们达到时所需的时间，是否满足任务要求，即是否能在两个雨季之间完成疏干任务。如符合需要，预测就算完成；不符合，则还要重复进行，直至所选取的 最佳疏干量满足任务要求的与时为止。 从图148取＝90m，则为20000m3／d从图中求得200d；可行性检验200＜210d，故符合技术要求。 继之，求雨季最大疏干量。 雨季地下水位上升，如以表示雨季的时段长，以表示水位上升幅度，为保证开采水平中段的正常生产，必须将雨季特别是丰水年雨季抬高的水头降下去。因此，雨季的最大疏干量应为开采水平正常疏干量即正常涌水量 ，亦即在前面所确定的最佳疏干量，再加雨季时段拾高所增加的疏干量，称疏干增量。则 图 9 广东曲塘矿区水文地质示意图 1-相对隔水边界；2-河流；3-断层；4-上升泉；5-强烈岩溶化灰岩；6-砂页岩、泥灰岩 上述计算，关键是雨季及其时段内地下水位上升幅度的确定。一般按动态观测资料给出抬高的平均值，较为可靠。将所得、代入前面所列公式，则可计算出雨季增加的疏干量。 2. 稳定涌水量的计算 广东曲塘多金属矿，位于一构造盆地边缘，地势平缓，雨量充沛，地表水系发育。矿体位于当地侵蚀基准面以下，赋存于含水性差的晚泥盆系天子岭组泥灰岩中。由于地层缺失，在某些地段使强烈岩溶化的中晚石炭系壶天群灰岩直接复盖其上，构成矿区主要充水层。此外，地表分布有弱含水的第四系冲积粘土夹砾石。矿区北、西有隔水层，东、南开阔，有地下水与地表水联系密切。根据边界的概化（图）选择直交隔水边界的稳定流“大井”公式，计算各开采中段稳定涌水量 将－40m中段的矿坑涌水量计算的参数和结果列入表 表 曲塘矿－40米中段矿坑涌水量计算表 计算参数 预计矿坑涌水量计算表（m3/d） 实际排水量 m3/d 误差 m m m m/d m m m m 壶天群 天子岭组 第四系 总计 149.77 151.59 47.80 2.374 4345 455 391.50 55078.50 28120 4620 910 33650 39250 16.6 表中的，其中为比例系数，系根据多落程水位降低的抽水试验求得。 五、数值法 （一）应用条件 数值计算用近似分割原理摆脱解析法处理实际问题时的严格理想化要求，使其更接近实际，因此主要用于水文地质条件复杂的大水矿床，并依据大流量抽放水试验为其对水文地质条件整体暴露，并提供建模、模型识别、大降深疏干预测的各种信息资料。 矿坑涌水量数值计算其原理方法虽与供水水资源评价完全一致，但由于矿床所处自然环境复杂，开采条件变化大，不确定因素多，又要求作大降深下推预测。因此，矿坑涌水量数值计算的最大特点是模型识别的条件差、任务重、难度大。不仅要为原始状态下水文地质模型的各项未知条件与不确定因素，通过定量化过程得到识别与校正；同时，还要为大降深数值预测建立内边界的互动机别，即隨内边界（面积、降深）变化外边界的下推规律及其水均衡条件。 （二）计算方法与步骤 1. 数学模型的选择 数学模型的建立，既要考虑需要，又要分析其可能与效果。即实际问题的复杂程度是否 具有所选模型相应的资料，一般来说，，平面二维数学模型已能满足解决实际问题的基本要求。但对于由弱透水层连接的多层层状含水层层组结构，可以从实际出发选择准三维模型；即用平面二维问题刻划含水层的基本特性一，以垂向一维流描述含水层之间的作用；对于在垂向上具明显非均 质特征的巨厚含水层，在较大降深的开采量和水位预报时，为避免失真最好采用三维流数学模型。 下面以非均质二维非稳定流地下水模型为例进行讨论，它由泛定方程和定解条件（初始条件和边界条件）组成 泛定方程 初始条件 边界条件 式中导水系数。潜水，；承压水， 水头； 承压含水层厚度； 潜水为给水度，承压水为储水系数； 单位时间、单位面积上的垂向转化量； 水位边界； 流量边界； 计算域； 平面坐标； 时间。 泛定方程是一个描述地下水渗流场收、支平衡的水均衡方程，其对水量转化规律的刻画是以达西定律为依据的，它由两部分组成一是均衡