认证主体：于**（实名认证）
IP属地：天津
下载本文档
1、文档供参考，可复制、编制，期待您的好评与关注!吉林大学精品课 >> 专门水文地质学 >> 教材 >> 专门水文地质学§ 10.4 坑涌水量预测一、矿坑涌水量预测的内容、方法、步骤与特点（一）矿井涌水量预测的内容及要求矿坑涌水量预测是一项重要而复杂的工作，是矿床水文地质勘探的重要组成部分。矿坑涌水量是指矿山开拓与开采过程中，单位时间内涌入矿坑（包括井、巷和开采系统）的水量。通常以 m3/h表示。它是确定矿床水文地质条件复杂程度的重要指标之一，关系到 矿山的生产条件与成本，对矿床的经济技术评价有很大的影响。并且也是设计与开采部门 选择开采方案、开采方法
2、，制定防治水疏干措施，设计水仓、排水系统与设备的主要依据。因此，在矿床水文地质调查中，要求正确评价未来矿山开发各个阶段的涌水量。其内容与要求包括可概括为以下四个方面：（1）矿坑正常涌水量：指开采系统达到某一标高（水平或中段）时，正常状态下保持相对稳定的总涌水量，通常是指平水年的涌水量。（2）矿坑最大涌水量：是指正常状态下开采系统在丰水年雨季时的最大涌水量。对某 些受暴雨强度直接控制的裸露型、暗河型岩溶充水矿床来说，常常还应依据矿山的服务年 限与当地气象变化周期，按当地气象站所记录的最大暴雨强度，预测数十年一遇特大暴雨 强度产生时，可能出现暂短的特大矿坑涌水量，作为制订各种应变措施的依据。（3）
3、开拓井巷涌水量：指包括井筒 （立井、斜井）和巷道（平、平巷、斜巷、石门 ）在开 拓过程中的涌水量。（4）疏干工程的排水量：是指在规定的疏于时间内，将一定范围内的水位降到某一规 定标高时，所需的疏干排水强度。对于地质勘探阶段来说，主要是进行评价性的计算，以预测正常状态下矿坑涌水量及 最大涌水量为主。至于开拓井巷的涌水量预测和专门性疏干工程的排水量的计算，由于与 矿山的生产条件密切相关，一般均由矿山基建部门或生产部门承担。（二）矿坑涌水量预测的方法根据当前矿床水文地质计算中常用的各种数学模型的地质背景特征极其对水文地质模型概化的要求，可作如下类型的划分：1 / 25文档供参考，可复制、编制，期待您
4、的好评与关注!非确定性统计模型Q s曲线方程 回归方程解析解 井流方程稳定井流公式非稳定井流公式数学模型分类确定模型数值解有限元法有限差法水均衡法混合型模型(三)矿坑涌水量预测的步骤矿坑涌水量预测是在查明矿床的充水因素及水文地质条件的基础上进行的。它是一项 贯穿矿区水文地质勘探全过程的工作，一个正确预测方案的建立，是随着对水文地质条件 认识的不断深化，不断修正、完善而逐渐形成的，一般应遵循如下三个基本步骤：1 .选择计算方法与相应的数学模型详勘阶段均要求选择 2个或2个以上的计算方法，以相互检验，映证。选择时必须考 虑三个基本要素：(1)矿床的充水因素及水文地质条件复杂程度。如：位于当地侵蚀基
5、准面之上，以降 水入渗补给的矿床，应采用水均衡法；水文地质条件简单或中等的矿床，可采用解析法或 比拟法；水文地质条件复杂的大水矿床，要求采用数值方法；(2)勘探阶段对矿坑涌水量预测的精度要求；(3)勘探方法、勘探工程的控制程度与信息量：如：水均衡法，要求不少于一个水文年的完整均衡域的补给与排泄项的动态资料；Qs曲线方程外推法，要求具抽水试验的水位降达到预测标高水柱高度的1/21/3;解析法，要求勘探工程全面控制含水层的非均质各向异性、非等厚的结构特征及其边界条件与补给、径流与排泄，并提供数值模型的建立、识 别、预测所需的完整信息数据，这些数据的获取，只有采用大型抽、放水试验对渗透场进 行整体控
6、制与揭露才可能做到。因此，计算方法与相应数学模型类型的选择，与矿床的充水因素及水文地质条件复杂 程度、勘探方法勘探工程的控制程度及信息量是相互关联的，统一在最佳技术经济条件这 一原则下，所以数学模型类型选择是否合理，可以用以下标准衡量：一是对矿床水文地质条件的适应性：指能否正确刻画水文地质条件的基本特征；二是对勘探方法勘探工程控制程度的适应性：指是否最充分的利用勘探工程提供的各 种信息，即信息的利用率；同时，也可理解为所选数学模型要求的勘探信息是否有保证， 即信息的保障率。2 .构造水文地质模型矿坑涌水量预测中数学模型的作用，是对水文地质条件进行量化，因此预测精度主要 取决于对充水因素与水文地
7、质条件判断的准确性，由于不同数学模型类型对水文地质条件 的刻画形式与功能各异，因此必须按数学模型的特点构造水文地质模型，称水文地质条件 概化。概化后的水文地质模型称水文地质概念模型，它在地质实体与数学模型之间起中介 桥梁作用。下面以最基本的预测方法一一解析法与数值法为例作一讨论。2 / 25文档供参考，可复制、编制，期待您的好评与关注!(1)概化已知状态下的水文地质条件；(2)给出未来开采状态下的内边界条件；(3)预测未来开采状态下的外边界条件。解析法将复杂的含水层结构与内外边界，以理想化模式构造理论公式，因此必须按解 析解要求进行概化。如含水层均质等厚，内外边界几何形态规则，边界供水条件简单
8、、确 定。数值法以近似分割原理对复杂的含水层结构、内外边界条件进行量化“逼真”，概化时要求以控制水文地质条件与内外边界的节点参数、水位与流量来构造水文地质概念模型。随着数学模型研究的不断进展，现代水文地质计算对水文地质模型的要求越来越高。目前，对复杂的大水矿床来说，一个可靠的水文地质模型的建立，必须贯穿整个勘探过程， 并大致经历三个阶段。即：第一阶段：通过对以往资料的整理，提出水文地质模型的“雏型”，作为下一步勘探设计的依据。尤其对大型抽(放)水试验的设计具有特殊重要的意义；第二阶段：根据进一步勘探提供的各种信息数据。特别是大型抽(放)水资料，通过流场分析或数值模拟，完成对“雏型”模型的调整，
9、建立水文地质模型的“校正型”；第三阶段：在“校正型”的基础上，按开采方案给出疏干工程的内边界条件，根据勘探资料预测不同疏干条件下的外边界条件，建立水文地质概化模型的“预测型”。3.计算数学模型，评价预测结果应该指出，不能把数学模型的解仅仅看作是一个单纯的数学计算，而应看作是对水文地质模型和数学模型进行全面验证识别过程，也是对矿区水文地质条件从定性到定量在回 到定性的不断深化的认识过程。(四)矿坑涌水量预测的特点虽然矿坑涌水量预测的原理方法与供水水资源评价类同，但其预测条件、预测要求与思路各有不同。如：(1)供水水资源评价，以持续稳定开采确保枯水期安全开采量为目标；而矿坑涌水量 预测则以疏干丰水
10、期的最大涌水量为目标；(2)矿床大多分布于基岩山区。含水质的非均质性突出。参数代表性不易控制；边界 条件复杂、非确定性因素多，常出现紊流、非连续流与管道流。定量化难度大。(3)矿山井巷类型及其分布千变万化，开采方法、开采速度与规模等生产条件复杂且 不稳定，与供水的取水建筑物简单、分布有序、生产稳定形成显明对比，给矿坑涌水量预 测带来诸多不确定性因素。(4)矿坑涌水量预测多为大降深下推。此时开采条件对水文地质条件的改变难以预料 和量化，这与供水小降深开采有明显差异。(5)矿床水文地质勘探从属于矿产地质勘探，与专门性的供水水文地质勘探对比，前 者一般投入小、工程控制程度低，预测所需的信息量相对少而
11、不完整。以上特点，决定了矿坑涌水量预测中存在诸多产生误差的客观条件。因此属于评价性 计算，为矿山设计及采前进一步专门性补充勘探提供依据。二、QS曲线方程外推法3 / 25文档供参考，可复制、编制，期待您的好评与关注!(一)原理与应用条件指用稳定井流条件下抽水试验的Q=f (s)方程，外推未来疏干水位降的涌水量。实质上也是一种相似条件下的比拟法。应用时的前提条件是：一、抽水试验建立Q=f (s),应符合稳定井流条件；二、抽水试验的各种条件应与预测对象的疏干条件接近。因此，必 须重视试验的技术条件，包括：1 .应将抽水试验孔布置在预测对象的分布地段，保证水文地质条件的一致性；2 .采用大口径(或孔
12、组)试验，计算时为消除井径对涌水量的影响，需做井径换算；3 .抽水降深应大于疏干水位水柱高度的1/21/3,计算时的外推疏干降深不应超过1.75倍的抽水降深，主要考虑疏干状态下的补给条件；4 .用枯季抽水试验预测正常涌水量，根据雨季试验预测季节性最大涌水量；5 .要排除抽水过程中一切自然和人为随机影响因素的干扰。Qs曲线方程法的优点是：回避各种水文地质参数求参过程中的失真，计算简单易行。 适用于建井初期的井筒涌水量预测。上水平疏干资料外推下水平的涌水量，以及矿床规模 小、矿体分布集中、边界条件和含水结构复杂的涌水量预测。(二)计算方法与步骤1 .鉴别Qs曲线类型(1)曲度法：即用曲度 n值进行
13、鉴别：，当n=1时，为直线Q=qS; 1vnv2时，为哥 曲线Q=; n=2时，为抛物线 S=aQ+bQ2; n>2时，为半对数曲线 S=a+blgS。如果nv 1时， 表明抽水试验不正确。2 .确定方程参数(1)最小二乘法：应根据 Q=f (s)类型选用最小二乘法，如常见的募函数型：(2)图解法：即利用直角坐标的图解，a为图解中纵坐标上所切的截距线段；b为直线对水平倾角的正切。其它类型详见地下水动力学。3 .井径换算(1)根据稳定井流理论中，不同地下水运动(层流、紊流)状态下井径与涌水量关系 进行换算，如层流：；紊流：。试验表明：井径对涌水量的影响，一般比对数大，而小于平方根。(2)用
14、二次或二次以上不同孔径的抽水试验资料建立井径d的换算公式：Q=mdn, m、n值可用最小二乘法求得。4 .水均衡论证将预测结果进行疏干条件下的水均衡论证，分析所得的最大涌水量是否符合疏干状态 下的边界补给条件。(三)实例1 .广东某金属矿，曾用Qs曲线法预测+50m水平的涌水量为14450m3/d,与坑道放 水试验的外推涌水量(14000m3/d)接近，而解析法预测的结果 (12608m3/d)则偏小12%。4 / 25文档供参考，可复制、编制，期待您的好评与关注!2 .湖南叶花香铜矿，1970年采用试生产坑道的排水量与水位降资料。分别用两种曲线 类型对未来四个中段的涌水量进行外推预测，后经实
15、际检验见表 14-2。5 / 25文档供参考，可复制、编制，期待您的好评与关注!表 14-2计算方法不同降深的水量52.45 m73.33 m99.47 m水量（m/d）误差（为水量（m-）误差（%水量（nVd ）误差（%计算实际计算实际计算实际S=aQ+bQ1514516937-10.601822323644-22.902151431155-31.00Q=14840-12.4018093-23.5021582-30.70实际检验表明预测结果与实际情况较接近，产生误差的主要原因是其间曾发生两次大 突水淹井，大量泥沙溃入矿坑，地面出现岩溶坍陷和裂缝所致。三、相关分析法（一）原理与应用条件相关分析
16、是根据涌水量与主要影响因素之间相关关系的密切程度建立回归方程，利用 抽水试验或开采初期的疏干资料，预测矿坑涌水量或外推开采后期下水平的涌水量。其原 理已在供水中介绍。根据实际资料的统计，多元复相关预测远比单相关效果好，其回归方 程表达的内容丰富，可反映除降深外的各种影响因素。它的应用条件与 Qs关系方程类同,但对原始数据的采集有严格要求：1代表性：（规范）要求不少于一个水文年（包括丰、平、枯季节）的动态观测数据， 同时数据（择本）量不少于 30个；2一致性：指应与预测对象上条件相一致；3独立性与相关性：即多自变量有独立的变化规律，相互间关系不大；而与涌水量之间均存在密切的相关关系，（规范）要求
17、相关系数不低于0.7。（二）实例与计算方法1.利用勘探阶段抽水试验资料预测矿坑涌水量如广东沙洋矿通过在勘探阶段设计相距6m的两个抽水孔和十余个不同距离的观测孔组成的群孔抽水试验，取得了复相关计算所需的涌水量Q与井径r （是将距抽水孔不同距离观测孔的位置概化为疏干状态下的坑道系统不同面积的作用半径）、水位降S （即不同作用半径的水位降，以模拟疏干水位降）有关资料，（见表3）通过求参建立了复相关募函数预测方程：其复相系数达 0.9468,复相关机误仅 0.0721,完全可用于未来矿山各设计水平与面积的矿坑涌水量预测。经实际排水资料检验，预测误差偏小3856%,主要与开采导函大量地面岩溶坍陷有关。6
18、 / 25文档供参考，可复制、编制，期待您的好评与关注!34.49155.1681.1472.0531.7052.0331.5671.7841.4271.39369.1452.9842.9022.1162.4742.0051.8611.3082.利用矿山观测资料外推预测可充分考虑矿坑涌水量的增长和各项生产因素间的关系，并根据它们之间的密切程度来建设涌水量方程。在原苏联顿巴斯煤矿的某些涌水量预测中，首先，在30个矿井中建立了 320个观测点，获得了涌水量（Q2）与各生产因素（包括矿产量 Po、开采深度Ho、开采 面积F。、生产时间To等）之间的相关关系，以及其密切程度，见表 14-4。表 14-
19、4相关系数r、XlgQ0lgP0lgH0lgF0lgT0lgQ0一0.6640.4510.5930.175lgP00.664个0.3400.6800.323lgH00.4510.340一0.5590.523lgF00.5930.6800.559一0.778lgT00.1750.3230.5230.778一根据判别得知，生产时间To对的影响不大（相关系数为 0.175）。用多元复相关计算,求得四元复相关曲线回归方程：其复相关系数：。在此基础上建立了比拟公式:预测结果与传统的单位涌水量法相比，使误差减少1.4倍（式中P、H、F为设计值）。四、解析法（一）解析法的应用条件解析法是根据解析解的建模要求
20、，通过对实际问题的合理概化，构造理想化模式的解 析公式，用于矿坑涌水量预测。具有对井巷类型适应能力强、快速、简便、经济等优点， 是最常用的基本方法。解析法预测矿坑涌水量时，以井流理论和用等效原则构造的“大井” 为主，后者指将各种形态的井巷与坑道系统，以具有等效性的“大井”表示，称“大井” 法。因此说：矿坑涌水量计算的最大特点是“大井法”与等效原则的应用，而供水则以干 扰井的计算为主。稳定井流解析法：应用于矿坑疏干流场处于相对稳定状态的流量预测。包括在已知 某开采水平最大水位降条件下的矿坑总涌水量；在给定某开采水平疏干排水能力的前提 下，计算地下水位降深（或压力疏降）值。非稳定解析法：用于矿床疏
21、干过程中地下水位不断下降，疏干漏斗持续不断扩展，非 稳定状态下的涌水量预测。 包括：已知开采水平水位降 （s）、疏干时间（t）,求涌水量（Q）;7 / 25文档供参考，可复制、编制，期待您的好评与关注!已知Q、s,求疏干某水平或漏斗扩展到某处的时间( t);已知Q、t,求s,以确定漏 斗发展的速度和漏斗范围内各点水头函数随时间的变化规律，用于规划各项开采措施。在 勘探阶段，以选择疏干量和计算量最大涌水量为主。(二)计算方法如上所述，应用解析法预测矿坑涌水量时，关键问题是如何在查清水文地质条件的前 提下，将复杂的实际问题概化。它可概括为如下三个重要方面：分析疏干流场的水力特征， 合理概化边界条件
22、，正确确定各项参数。1.分析疏干流场的水力特征矿区的疏干流场是在天然背景条件下，迭加开采因素演变而成。分析时，应以天然状 态为基础，结合开采条件作出合理概化。(1)区分稳定流与非稳定流矿山基建阶段，疏干流场的内外边界有受开拓井巷的扩展所控制，以消耗含水层储量 为主，属非稳定流；进入回采阶段后，井巷输廊大体已定，疏干流场主要受外边界的补给 条件控制，当存在定水头(侧向或越流)补给条件时，矿坑水量被侧向补给量或越流量所 平衡，流场特征除受气候的季节变化影响外，呈现对稳定状态。基本符合稳定的“建模” 条件，或可以认为两者具等效性；反之，均属非稳定流范畴。如河北开滦煤矿，其矿坑涌水量随坑道走向的延展而
23、增加，但这种延展暂停时，涌水 量立即出现相对的稳定。此时仅表现有受降水的季节变化在一定变幅范围内上下波动，并 出现强出水点和边缘出水点成袭夺中间出点现象，而总涌水量不变。又如辽宁复州粘土矿， 其涌水量随采深增加，但某一水平进入回采时，其涌水量就逐渐稳定，并保持到下一水平 突水进止，在此阶段虽然也出现下水平突水点袭夺上水平突水点现象，但总涌水量却保持 不变。由此可见，在某些矿区的疏干过程中，不仅存在疏干流场的相对稳定阶段，而且随 矿山工程的进展而不断相互转化。但选用稳定流解析法时要慎重，必须进行均衡论证，判断疏干区是否真正存在定水头 供水边界或定水头的越流系统。此外，对于地下水储存量较大的矿区，
24、要单独计算疏干过 程中不同阶段含水层储存量的放强度，大量生产实际表明，矿坑最大涌水量，并非出现在疏干过程的稳定阶段。(2)区分达西流与非达西流在矿坑涌水量计算时，常遇到非达西流问题，它涉及解析法的应用条件，在宏观上可 概括有两种情况：一是暗河管道岩溶充水矿床，地下水运动为压力管道流与明渠流；此外，分水岭地段 的充水矿床，矿坑涌水量直接受垂向入渗降雨强度控制，与水位降深无关。两者均与解析 法的“建模”条件相距甚大，矿坑涌水量预测应选择水均衡法或各种随机统计方法。二是，局部状态的非达西流，常发生在大降深疏干井巷附近与某些特殊构造部位，它 只对参数计算与参数的代表性产生影响。在宏观上，它是一个流态概
25、化问题，不存在解析 法的应用条件问题。(3)区分平面流与空间流严格讲，在大降深疏干条件下，地下水运动的垂向速度分量不能忽略，均为三维空间8 / 25文档供参考，可复制、编制，期待您的好评与关注!流(包括非完整井巷的地下水运动)问题，其分布范围仅限于井巷附近，均为含水层厚度的1.54.75倍。因此，在矿坑涌水量预测中，大多将其纳入二维平面流范畴，在宏观上不影响预测精度。计算时应根据井巷类型作出不同的概化：如：竖井的涌水量计算，可概化为平面径向流问题，以进流公式表达；计算水平巷道 涌水时，以剖面平面流近似，采用单宽流量解析公式，但其两端上往往也产生辐射流(见 图)，需要考虑它的存在，并采用平面径向
26、流公式补充计算巷道端部的进水口。坑道系统则复杂得多，根据“大井法”原理，一般以近似的径向流概化，但当坑道系 统近于带状的狭长条形时，也可概化为剖面流问题。对于倾斜坑道，根据阿勃拉莫夫有关水电比拟法的研究，证明坑道的倾斜对涌水量影响不大，可根据坑道的倾斜度，分别按竖井或水平巷道进行近似。即：若坑道倾斜度450时，视其与竖井近似，用井流公式计算；若坑道倾斜度V450时，则视其与水平巷道近似，用单宽流量公式计算。根据解析解的存在条件，一些简单的非完整井巷涌水量计算，可以运用三维空间问题 予以解决。此时，可根据非完整井的特点，运用地下水动力学中映射法与分段法的原理来 求解。通常用平面分段法解决完整竖井
27、的涌水量计算，用剖面分段法解决非完整平巷的涌 水量计算。(4)区分潜水与承压水与供水不同，在降压疏干时，往往出现承压水转化为潜水或承压-无压水。此外，在 陡倾斜含水层分布的矿区，还可能出现坑道一侧保持原始承压水状态，而另一侧却由承压 水转化无压水或承压-无压水的现象。概化时，需从宏观角度作等效的近似处理。2 .边界条件的概化边界条件概化的失误是导致解析解失真的主要原因之一。由于理想化要求常与实际条 件相差甚远，成为解析法应用中的难点，也是解析法预测矿坑涌水量的重要环节。(1)侧向边界的概化图14-3 铜录山矿区边界概化图1-Mb 大理岩；2丫6岩浆岩。图14-2 恩口矿区边界条件转化示意图1T
28、is+P2d 下叠大冶组；2P21上叠龙潭组隔水层； 3 Pim+Piq4下二叠茅口组与栖霞组岩溶含水层； 4-P1q3下二叠栖霞组李子塘段隔水层；5- I , 口，hi ,疏干水平。解析法要求将复杂的边界补给条件概化为隔水与供水两种进水类型；同时，将不规则的边界形态，简化为规则的。但实际问题中一般难以具有上述理想条件，其进水条件常9 / 25文档供参考，可复制、编制，期待您的好评与关注!常既不完全隔水，又不具有无限补给能力，它的分布也极不规则。为此，必须通过合理的概化，缩小理论与实际的差距，满足近似的计算要求。其要点是：立足于整体概化效果；以均衡为基础，用好等效原则。等效原则是边界概化中的无
29、奈之举，即：通过对概 念（如相对隔水边界、近似定水边界）寻找近似处理的途径；或根据等效原则将垂向越流 补给和侧向补给共同构造定水头边界，将局部进水口概化为区域进水边界等。但这些等效 原则的应用，必须建立在区域水均衡条件论证的基础上，并涉及参数的优化处理。充分考虑开采因素。疏干流场始终处于补给量与疏干量不断变化的动平衡状态，随着开采条件的变化，边界的位置及其进水条件常发生转化，如湖南恩口煤矿的东部边界（见图2）,在I水平疏干时东部壶天河不起作用；开采延伸至n水平时，因排水量增大漏斗扩 展到壶天河，成为茅口灰岩的定水头供水边界；当疏干达到出水平时，排水量随降深继续 增加，当壶天河的补给能力无法与其
30、平衡时，其定水头供水边界已不复存在，漏斗扩展至 由隔水层构造的隔水边界，但壶天河仍以变水头集中补给形式平衡疏干漏斗的发展。概化 时，应与西部边界的供水条件作统一的整体骰考虑，如仅就东部边界而言，可用等效原则 按第一类越流边界处理，但须从均衡出发，确定一个相当于第一类越流作用的“引用越流 系数”取代；此外，也可单独计算壶天河的渗漏量，作为矿坑涌水量的一部分。边界几何形态的概化也需认真对待。如湖北铜汞山铜矿的露天矿涌水量预测：矿坑充水来自围岩大理岩，与东西两侧岩浆岩隔水层呈似以300交角，向南敞开（见图 3）所。60年代勘探时，概化为东侧直线隔水的环状供水边界，采用非完整井稳定井流域，预测矿坑涌水
31、量为59587985m3/d,而实际涌水量仅 3790 m3/d,误差57111%。70年代回访调查 验证计算时，采用 300扇形补给边界的稳定流近似计算，得涌水量 3685 m3/d,周期实际涌 水量为3416 m3/d,误差仅7.8%。证明边界形态概化的重要性。边界概化应把重点放在主要供水边界上。孙纳正教授运用数值法，与解析法对理想 化模型的对比验证计算表明：简化供水边界的形状往往会带来较大的误差，但简化隔水边 界的形状影响一般不大。（2）各种类型侧向边界条件下的计算方法映射法。即根据地下水动力学中的映射迭加原理，获得矿坑涌水量预测的描述各种特定边界条件下的解析公式。可采用如下一般形式表示
32、：稳定流非稳定流式中，与分别为稳定流与非稳定流的边界类型条件系数。详见地下水动力学及有关文献。分区法。也称卡明斯基辐射流法。它是从研究稳定状态下的流网入手，根据疏干流场的边界条件与含水层的非均质性特点，沿流面和等水压面将其分割为若干条件不同的扇形分流区（见图14-4）,每 图14-4某矿区辐射流计算图个扇形分流区内其地下水流都呈辐射流，其沿流面分割所得的各扇形区边界为阻水边界，而沿等水压面分割所得的扇形区边界为等水头边界。常用卡明斯基平面辐射流公式10 / 25文档供参考，可复制、编制，期待您的好评与关注!分别计算各扇形区的涌水量Qi潜水承压水式中，与为分流区辐射状水流上下游断面的宽度；与为与
33、断面隔水底板上的水头高度; 为与断面的间距。然后，按下式求各分区流量的总和：每个扇形区内的下游断面，是以直接靠近井巷的硫干漏斗等水头线的一部分为准；而 上游断面则以远离井巷的供水边界上等水头线面一部分为准。(3)垂向越流补给边界类型的确定及其计算当疏干含水层的顶底板为弱透水层时，其垂向相邻含水层就会通过弱透水层对疏干层 产生越流补给，出现所谓的越流补给边界。越流补给边界分定水头和变水头二类，解析法 对后者的研究还尚待解决。产生定水头垂向越流补给的矿坑涌水量计算，可用增加越流参数项的形式来表示： 稳定流非稳定流式中：一一越流参数TM垂向弱透水层渗透系数；一一垂向弱透水层厚度;Kg 零阶二类修正贝
34、塞尔函数。3 .参数确定(1)渗透系数(K)值加权平均值法: 式为：K cp渗透系数是解析公式中的主要参数。我国矿山大多为分布于基岩山区的裂隙、岩溶充 水矿床，充水含水层的渗透性具明显不均匀性，根据解析计算要求，应作均值概化，同时 这也是保证渗透系数具有代表的措施之一。矿坑涌水量预测中常用的方法有两种。又可分为厚度平均、面积平均、方向平均法等。如厚度平均，则公nMi(Hi)Ki i 1 nMi(Hi) i 1式中：-承压(潜水)含水层各垂向分段厚度;Ki-相应分段的渗透系数。流场分析法：有等水位线图时，可采用闭合等值线法:11 / 25文档供参考，可复制、编制，期待您的好评与关注!或据流场特征
35、，采用分区法：KCPQ n 3 b2）h； * i i Inbllnb2 2L式中：、一任意两条（上、下游）闭合等水位线的长度；，'两条闭合等水位线的平均距离；两条闭合等水位线的平均距离；-含水层的平均厚度；涌水量；、-辐射状水流上、下游断面上的宽度；、一和断面隔水底板以上的水头高度；和断面之间的距离。（2）大井引用半径的确定矿坑的形状极不规则，龙其是坑道（井巷）系统，分布范围大，形状千变万化，构成了复杂的内边界。根据解析法计算模型的特点，要求将它理想化。经观测，坑道系统排水时， 其周边逐渐形成了一个统一的降落漏斗。因此，在理论上可将形状复杂的坑道系统看成是 一个理想“大井”在工作，此
36、时整个坑道面积，看成是相当于该“大井”的面积。整个坑 道系统的涌水量，就相当于“大井”的涌水量，这样就使一般的井流公式能适应于坑道系 统的涌水量计算。这种方法，在矿坑涌水量预测中称为“大井法”。“大井”的引用半径,在一般情况下用下式表示：式中，为坑道系统分布范围所圈定的面积。确切地说，近似等于为保证并田设计生产 率所必需的坑道所圈定面积的大小，或者以降落漏斗距坑道最近处的封闭等水位线所围起 来的面积。如果开采面积近于圆形、方形时，采用上式较准确，对于形状特别的面积，可 采用其它专门公式计算。（3）引用半径的确定对于在实际问题中是否存在裘布依关于园形定水头的假定条件，以及齐姆模型中影响 半径概念
37、是否合理，学术界存在众多争议，特别是非稳定流理论与越流理论的产生，人们 认识到某种稳定状态的出现可能是越流作用所造成，这样影响半径的概念离开了原有的含 义，于是引出了所谓“引用影响半径” 。从稳定井流理论的实际应用出发，上述概念是可以 被接受的，即根据等效原则，将疏干量与补给量相平衡时出现的稳定流场，其边界用一个 引用的园形等效外边界进行慨化，其与“大井”中心的水平距离称为引用影响半径，也称 为补给半径。即。同理，在用割面平面流解析公式计算使表水中坑道涌水量，也就有了引 用影响带宽度（L0）的概念，即疏干坑道中心与外边界之间的距离。在稳定流条件下，引用影响半径为一个常量，也称补给半径；在非稳定
38、流条件下则是 一个不断变化着的变量，这样在理论上解决了稳定井流理论及其引用影响半径计算公式的 实用问题。12 / 25文档供参考，可复制、编制，期待您的好评与关注!矿山疏干实际表明，矿坑排水的影响范围，总是随时间的延长、排水量的增加以及坑 道的推进而不断扩大，直到天然边界为止，它不可能被限制在一个不是边界的理想“半径” 之内。此外，对比计算表明，若确定影响半径的误差为23倍，则矿坑涌水量的计算误差可达30%60%;此外，若取偏低值其误差远比取偏高值要大。因此，矿坑涌水量预测时，能否 用解析公式及常见的经验公式来近似的确定影响半径值得探讨；对开拓井巷的捅水量预测， 最好采用抽水试验外推法，即根据
39、多落程的抽水试验，确定降深与影响半径或流量与影响 半径的线性关系，外推某琉干水位或某疏于量的相应疏干半径值。如：;或对坑道系统的涌水量预测，应根据疏干中心天然水文地质边界线之间距离的加权平均 值计算，即塞罗瓦特科公式：式中：一一“大井”的引用半径；井巷轮廓线与各不同类型水文地质边界间的平均距离；-各类型水文地质边界线的宽度。（4）最大疏干水位降深的确定在理论上，目前解析解还无法处理承压区与无压区同时并存与大降深的潜水问题，对 于实际问题，则是矿床硫干时最大可能水位降是多少，如何近似确定最大硫于水位降深值。爱尔别尔格尔在实验中取得的潜水最大水位降深等于潜水含水层一半的结论，即（扩大应用到承压含水
40、层时，），一直是水文地质计算中所遵循的概念。近年来，我国通过渗流槽 及野外抽水试验.证明这一结论是保守的。可以超过，在矿坑涌水量计算中，通常不考虑这一概念。据观测，在长期硫于条件下的大截面井巷系统外缘，动水位（）一般不超过12m,它所引起的涌水量计算偏大值一般为0.5%1 %。因此，矿坑涌水量预测时，最大疏干水位降一般取。另一个理论问题，即最大水位降时的最大涌水量计算问题。众所周知，当时，裘布依 公式在理论上就会“失真”，这正是稳定井流理论的最大缺陷之一；而泰斯公式则是从承压 水含水层建立起来的，扩大到无压含水层使用时，（作最大降深疏干时，承压含水层均转化为无压水层），常把随时间变化的含水层厚
41、度作线性处理，即取不变的平均值，这种线性化 处理必然带来误差，据研究当降深超过含水层厚度30%时，非稳定井流公式要偏实际情况，出现明显误差，更不用说是作最大水位降的计算了。如上所述，不难看出矿坑涌水量预测 时，作最大水位降的最大疏干量计算，对解析法来说不是很合科适宜的。（三）实例1最佳疏干量（）的确定某铁矿地处灰岩区，裂隙岩溶发育较均匀，地下水运动符合达西定律，矿区内有部分13 / 25文档供参考，可复制、编制，期待您的好评与关注!地下水动态长期观测资料，其它地质条件略。（1）要求当疏干水平（或中段）的水位降深（S）确定后，则疏干量（）是时间（）的函数。这样，疏干 量就是与疏干时间有关的一组数
42、据。某水平的正常疏干量，应是该水平预测的矿坑涌水量值。设计部门要在一组具不同疏干强度及与其相应的时间的对比中，选出最佳疏干方 案，即选择排水能力要求不过大，而疏干时间又不长的方案。疏干时间通常要求控制在两个雨季之间，否则的计算则无意义。（2）任务给定的条件是：疏干中段水位降（ S）确定为零米标高；疏干时间要求在两个雨 季间完成。（3）最佳疏干量的计算与分析第一步：初选疏干时间段。根据第二项任务，在现有地下水动态曲线（图）上初选三 个时间段，即270 d、210 d、150d,供计算分析。第二步：确定相应的值。根据给定的零米标高，从动态曲线图上确定出各霎时间段相 对应的S值，即：(d)t3 27
43、0210150S(m)S31009080第三步：求相应的值，利用公式（符号为常用地下水动力学符号） 在已知、t3、S3的条件下，求得相应的、，作为第四步分析的初值。第四步：绘制不同疏干强度条件下的曲线。在初值、的范围内，通过内插给出一组供进一步分析的疏干量数据。其公式为：图6不同疏干量条件下 S=f（t）曲线分析不同疏干量时的 S随的变化规律，并绘 制不同疏干量条件下的曲线。（图）第五步：绘制不同定降深 S条件下的曲线。根据图作出不同降深 S条件下的疏干量14 / 25文档供参考，可复制、编制，期待您的好评与关注!与时间的关系曲线 （图），进行不同S条件下，疏干量与疏干时间的对比分析。第六步：
44、绘制降深 S与最佳疏干量的关系曲线。根据图中各曲线的拐点，求出不同降深条件厂的最佳疏干强度，即拟稳定疏干量与降深的关系曲线（图148）。第七步：确定最佳硫干量，并检验其可行性。根据图148取得的不同降深的最佳图8降深与疏干量关系曲线图7 不同降深条件下 Q=f（t）曲线疏干量检验它们达到 S时所需的时间，是否满足任务要求，即是否能在两个雨季之间 完成疏干任务。如符合需要，预测就算完成；不符合，则还要重复进行，直至所选取的 最佳疏干量满足任务要求的S与时为止。从图148取S = 90m,则为20000m3/d从图中求得=200d ;可行性检验：200v210d, 故符合技术要求。继之，求雨季最大
45、疏干量。雨季地下水位上升，如以表示雨季的时段长，以 S表示水位上升幅度，为保证开采水 平（中段）的正常生产，必须将雨季（特别是丰水年雨季）抬高的水头S降下去。因此，雨季的 最大疏干量应为开采水平正常疏干量（即正常涌水量），亦即在前面所确定的最佳疏干量，再加雨季时段拾高 S所增加的疏干量，称疏干增量。则：区1 区2国3区4日 5 | P2C1+D 6响9设爰H H遮Utt”1-相对隔水边界；壑-潮f3缄层；4-上升泉；5- 强烈岩溶化灰岩；6-砂页岩、泥灰岩15 / 25上述计算，关键是雨季及其时段内地 下水位上升幅度的确定。一般按动态观测 资料给出抬高S的平均值，较为可靠。将 所得、S代入前面
46、所列公式，则可计算出 雨季增加的疏干量。2,稳定涌水量的计算广东曲塘多金属矿，位于一构造盆地 边缘，地势平缓，雨量充沛，地表水系发 育。矿体位于当地侵蚀基准面以下，赋存 于含水性差的晚泥盆系天子岭组泥灰岩 中。由于地层缺失，在某些地段使强烈岩 溶化的中晚石炭系壶天群灰岩直接复盖文档供参考，可复制、编制，期待您的好评与关注!其上，构成矿区主要充水层。此外，地表分布有弱含水的第四系冲积粘土夹砾石。矿区北、 西有隔水层，东、南开阔，有地下水与地表水了解密切。根据边界的概化（图）选择直交 隔水边界的稳定流“大井”公式，计算各开采中段稳定涌水量：R4Ln2 , 28r b1b2 b1 b2将-40m中段
47、的矿坑涌水量计算的参数和结果列入表表 曲塘矿40米中段矿坑涌水量计算表计算参数预计矿坑涌水量计算表（m3/d）实际排水量（m3/d）误差%(m)(m)(m)(m/d)(m)(m)(m)(m)壶天群天子岭组第四系总计149.77151.5947.802.3744345455391.5055078.50281204620910336503925016.6表中的，其中为比例系数，系根据多落程水位降低的抽水试验求得。五、数值法（一）应用条件数值计算用近似分割原理摆脱解析法处理实际问题时的严格理想化要求，使其更接近 实际，因此主要用于水文地质条件复杂的大水矿床，并依据大流量抽放水试验为其对水文 地质条件
48、整体暴露，并提供建模、模型识别、大降深疏干预测的各种信息资料。矿坑涌水量数值计算其原理方法虽与供水水资源评价完全一致，但由于矿床所处自然 环境复杂，开采条件变化大，不确定因素多，又要求作大降深下推预测。因此，矿坑涌水 量数值计算的最大特点是：模型识别的条件差、任务重、难度大。不仅要为原始状态下水 文地质模型的各项未知条件与不确定因素，通过定量化过程得到识别与校正；同时，还要 为大降深数值预测建立内边界的互动机别，即随内边界（面积、降深）变化外边界的下推 规律及其水均衡条件。（二）计算方法与步骤1.数学模型的选择数学模型的建立，既要考虑需要，又要分析其可能与效果。即实际问题的复杂程度是 否具有所
49、选模型相应的资料，一般来说，平面二维数学模型已能满足解决实际问题的基本要 求。但对于由弱透水层连接的多层层状含水层层组结构，可以从实际出发选择准三维模型； 即用平面二维问题刻划含水层的基本特性一，以垂向一维流描述含水层之间的作用；对于 在垂向上具明显非均质特征的巨厚含水层，在较大降深的开采量和水位预报时，为避免失真最好采用三维流数 学模型。下面以非均质二维非稳定流地下水模型为例进行讨论，它由泛定方程和定解条件（初 始条件和边界条件）组成：泛定方程：16 / 25文档供参考，可复制、编制，期待您的好评与关注!-Tx-Hx X初始条件:边界条件:_ HTx cos(n,x) x_ HTy cos(
50、n, y) | 2 yq(x, y,t)式中Tx 导水系数。潜水，；承压水，Tx k m水头；承压含水层厚度；S潜水为给水度，承压水为储水系数；单位时间、单位面积上的垂向转化量； 水位边界； 流量边界；G 计算域；平面坐标；时间。泛定方程是一个描述地下水渗流场收、支平衡的水均衡方程，其对水量转化规律的刻画是以达西定律为依据的，它由两部分组成：一是均衡基本项（、 S项），指方程带有水头 函数的偏导项，表征渗流场各均衡单元内部及相互间的水量交换。其中含的水量渗透基本 项，指渗流场水量的侧向交换条件，反映了含水层介质的渗透性、非均质性、含水层的几 何形态、渗流运动状态；面含S的水量储存与释放基本项，
51、指渗流水量的储存与消耗。二是垂向交换项（项），包括源、汇项（即计算域内各井的抽水或注水强度），垂向入渗补给和消耗以及越流项。在模型中应是一个给定的已知函数，但在实际中某些垂向交换量常常 是未知的，因此它也可引入参数（如降水入渗垂向越流系数等）在模型中参与求参。初始条件：是指开采初始条件的地下水水头，为已知条件。边界条件：在二维模型中仅指侧向边界条件。当已知边界水头变化规律时，可按已知 水位边界表示（），称一类边值问题；当已知边界的流量变化规律时，可用已知流量边界表 示（），称二类边值问题，其强度以单宽流量表示。由一类边界和二类边界共同组成的混合 边界，称混合边值问题。矿坑涌水量数值计算中，大多
52、采用混合边值，即以勘探工程控制矿区主要水量交换边 界，用表现资料给定一类边值，以免模型识别失真，解决数学模型求解的唯一性问题；同 时，在水均衡研究基础上，以流量边界模拟次要边界，参于调参与识别，解决工程量不是 的困难。2.水文地质条件的概化17 / 25文档供参考，可复制、编制，期待您的好评与关注!水文地质条件概化是数值计算中的一个重要环节。要求根据勘探资料按数值方法对实 际问题的特点进行概化。它反映了勘探信息的利用率和保证率，以及对水文地质条件的研 究程度，直接关系计算精度。（1）含水层结构的概化：包括含水层的空间形态与结构参数分区的概化。含水层的空 间形态，是利用含水层顶、底板标高等值线图
53、，给出每一剖分节点（离散点）坐标（）上 的含水层顶、底板标高，由模型自动识别含水层的厚度，完成几何形态的概化。含水层的 非均质结构参数分区，是在水文地质分区的基础上（即依据、S的分布特点，结合岩性和松散沉积物的成因类型、基岩的构造条件、岩溶地区的水动力条件，进行水文地质分区） 按水文地质条件的宏观规律和渗流运动的特点，在空间上渐变地进行参数分区及参数分级， 给出各分区参数的平均值及其上、下限，作为模型调试的依据。对取水层与相邻含水层相 互作用概化，一般要求地质模型给出与相邻含水层的连接位置与坐标，其连接方式可以是 断层，“天窗”或通过弱透水层的越流补给。（2）地下水流态的概化：当水位降较大时，
54、在开采井附近常出现复杂的非达西流与三 维流，此外某些局部的构造部位或岩溶发育地段，甚至出现非渗流或非连续流状态。但这 些复杂水流状态的分布范围一般不大，因此在宏观上仍可考虑用二维达西流进行概化。（3）边界条件的概化：数值法能较真实地模拟边界复杂的边界条件，它与数理统汁模 型相结合，可以处理解析法无能为力的各种非确定边界问题。概化时，要求根据边界分布 的空间形态，给出边界的坐标，确定边界作用的性质，有无水量交换及其交换方式，并根 据动态观测或抽水试验资料，用数理统计方法概化水位或流量的变化规律，并按不同时段 给出边界节点的水位或单宽流量。计算边界的选择与确定对数值计算的精度及其工程量的投资关系极
55、大。操作时应遵循 两个基本原则：一是在经济上要求以最小的工程控制边界条件；二是在技术上要求所确定 的主要边界，具有一定的工程控制，能为模型的识别、校正和预测提供可靠的计算数据。 具体表现在，首先，尽可能的取自然边界和确定性边界，以节约勘探工程和提高模型的可 靠性；其次，应避免置计算边界于源、汇项附近，并远离供水中心，以缩小边界条件概化 不当对计算结果的不良影响；此外，模型识别与预测的边界必须一致，否则模型识别的成 果将失去意义。在二维地下水模型中，垂向水量交换是作为水量附加项（项）列入方程中的，因此在 概化时应特别慎重。同时要求给出含水层中的人工抽（注）水井的坐标、类型及其抽（注） 水强度。（
56、4）初始条件的概化：按初始时刻各控制节点实测水位资料绘制的等水位线图，给出 各节点的水位作为初始条件。由于控制节点的数量有限，等水位线图的整理难免在一定的 随意性，在含水层结构或边界条件较复杂的情况下，最好利用模型的小步长运行进行校正。3 .计算区域的离散数值法根据分割近似原理，将一个反映实际渗流场的光滑连续的水头曲面，用一个由 若干彼此衔接无缝不重叠的三角形（有限元法）或方形、矩形（有限差分法）拼凑成的连 续但不光滑的水头折面代替，将非线问题简化为线性问题求解。按离散化要求部分时，首 先要选好控制性节点，它是具有完整水位资料的观测孔。由于观测孔的数量有限，要有许18 / 25文档供参考，可复制、编制，期待您的好评与关注!多插值点平补充，完
0/150
联系客服
本站为文档C2C交易模式，即用户上传的文档直接被用户下载，本站只是中间服务平台，本站所有文档下载所得的收益归上传人(含作者)所有。人人文库仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私，请立即通知人人文库网，我们立即给予删除！