3. 解析法计算所需水文地质参数获取
解析法(如大井法、集水廊道法、裘布依公式等)预测矿坑涌水量的可靠性,首先取决于参数取值是否贴近实际含水层。解析公式中真正起作用的往往只有几个参数:渗透系数 K、导水系数 T、储水系数 S(或弹性释水系数 μe)、给水度 μ 与影响半径 R。它们之间并非彼此独立:T=K·M(M 为含水层厚度),S 与 μ 分别控制承压与无压含水层的非稳定流过程,R 则把抽水试验的局部信息外推到矿坑尺度。因此本节的核心不是“查表填数”,而是弄清每个参数的物理含义、可用的现场试验手段、以及经验取值的适用边界。
一、渗透系数 K 与导水系数 T:现场抽水试验是主渠道
渗透系数 K 表征含水介质的透水能力(m/d),导水系数 T=K·M 表征整个含水层厚度的导水能力(m²/d)。对矿坑涌水量预测而言,T 往往比 K 更直接,因为解析公式中的降深—流量关系多由 T 控制。
现场获取的首选方法是单孔稳定流抽水试验,按裘布依公式反算:
承压完整井:K = 0.366·Q·(lg R − lg r) / (M·s)
无压完整井:K = 0.732·Q·(lg R − lg r) / [(2H − s)·s]
其中 Q—稳定流量(m³/d),s—降深(m),r—井半径(m),
M—承压含水层厚度(m),H—无压含水层初始水位至底板厚度(m)
实操要点:抽水前需洗井至水清砂净,正式抽水前做一次试验抽水确定合理降深;稳定标准一般取连续 8 h 内水位波动不超过平均降深的 1%、流量波动不超过 3%。至少做 2~3 个降深(如 1/3、2/3、最大降深),用 Q—s 曲线检验是否符合线性或幂函数关系,若出现明显弯曲,说明井损或三维流影响显著,应改用带观测孔的资料计算。有观测孔时优先用观测孔水位,可避开井损:K=0.366Q(lgr2−lgr1)/[M(s1−s2)]。
当含水层富水性弱、单井出水量小、难以形成稳定降深时,可改用注水试验(常水头或变水头)或压水试验(适用于基岩裂隙含水层),按相应规范公式换算 K。对于矿坑疏干涉及的强富水岩溶裂隙带,还可结合连通试验与钻孔流速流向测定判断优势径流通道,避免把局部高渗带误当作全区平均 K。
经验取值只在初勘阶段或资料匮乏时使用:砂砾石 K 约 10~100 m/d,中细砂约 1~10 m/d,粉砂约 0.1~1 m/d,黏性土常小于 0.01 m/d;基岩则强烈受裂隙发育程度控制,同一岩性可相差 2~3 个数量级。经验值必须标注来源与适用条件,并在解析计算中做敏感性分析,而不是当作确定值直接代入。
二、储水系数 S 与给水度 μ:区分承压与无压
承压含水层的储水系数 S 由水的压缩与含水层骨架压缩共同贡献,量级通常为 10⁻⁵~10⁻³,意味着承压含水层“弹性释放”的水量很小,但水头传导快。无压含水层的给水度 μ 是重力疏干孔隙率,砂砾层约 0.15~0.30,中细砂约 0.10~0.20,粉砂约 0.05~0.10,量级比 S 大 2~4 个数量级。这一差别直接决定非稳定流解析解(如泰斯公式、博尔顿模型)中降深随时间发展的形态:承压含水层早期降深增长快、后期平缓;无压含水层存在明显的“延迟给水”段。
现场获取 S 与 μ 的主要手段是非稳定流抽水试验:布置 1~2 个观测孔,记录抽水全过程(尤其停抽后的水位恢复)的水位—时间序列,用配线法或直线图解法(s—lgt、s—lg(t/r²))求参。承压含水层用泰斯公式配线求 T 与 S;无压含水层用博尔顿模型或纽曼模型,可同时识别早期弹性释水与后期重力给水两个阶段。若仅有单孔资料,可用恢复试验(剩余降深法)估算 T,S 的精度会明显下降,此时宜结合室内试验或经验值校核。
室内试验可测给水度:取原状砂样做饱和后重力排水试验,测定单位体积疏干水量;但原状样扰动、夹层缺失会使室内 μ 偏大或偏小,需与现场非稳定流结果对照。经验上,矿坑疏干预测中若无实测 S,承压段可按 10⁻⁴~10⁻³ 量级试算,无压段按岩性取 μ,并对 S 做 3~5 倍变幅的敏感性检验。
三、影响半径 R:最易被误用、也最需要论证的参数
影响半径 R 是解析公式把“无限含水层”截断为有限圆域的等效边界,并非真实存在的水位零降深点。现场确定方法有三类:一是抽水试验实测,用观测孔资料按 s—lgr 直线外推至 s=0 求 R,这是最可靠的做法;二是经验公式,如库萨金公式 R=2s√(K·H)(无压)或吉哈尔特公式,适用于初步估算;三是影响半径经验表,按岩性给出粗砂数百米、中细砂百余米、粉砂数十米的量级。
需要特别提醒:R 随降深增大而增大,且与 K、含水层厚度正相关,把一个小降深试验得到的 R 直接用于矿坑大降深(数十米至上百米)预测,会系统性低估涌水量。工程上更稳妥的做法是:用实测 R 标定模型后,按矿坑疏干的实际降深用经验公式外推,并在解析计算中给出 R 的上下限区间,形成涌水量的区间预测而非单点值。
四、参数取值的工程校核与常见误区
把上述参数代入解析公式之前,建议做三项校核。第一,量纲与厚度一致性:T=K·M 中的 M 必须与解析公式假定的含水层厚度一致,承压段用承压厚度,无压段用饱和厚度,矿坑疏干后水位降至含水层内部时,应按剩余饱和厚度分段计算。第二,参数来源标注:每个参数应注明是实测、经验公式还是经验表,实测值注明试验日期、降深与稳定标准,便于复核。第三,敏感性排序:一般 T 与 R 对涌水量影响最大,S 主要影响时间过程而非最终稳定流量,μ 在无压疏干预测中影响显著。可按下表组织参数取值与依据:
| 参数 | 首选获取方法 | 示例取值(教学假设) | 取值依据与备注 |
|---|---|---|---|
| 渗透系数 K | 单孔稳定流抽水,2~3 个降深 | 砂砾石层 25 m/d | Q—s 曲线线性良好,无明显井损 |
| 导水系数 T | T=K·M,M 取 18 m | 450 m²/d | 与抽水试验反算 T 对照偏差<15% |
| 储水系数 S | 非稳定流抽水配线法 | 2×10⁻⁴ | 承压段;按 3~5 倍变幅做敏感性 |
| 给水度 μ | 室内原状样 + 博尔顿模型反演 | 0.18 | 无压段;与室内试验值对照 |
| 影响半径 R | 观测孔 s—lg r 外推,经验公式校核 | 试验 180 m,外推至疏干降深 350 m | 给出上下限,不取单点 |
综上,解析法参数获取的关键在于“试验定值、经验定界、敏感性定区间”。现场试验给出参数的主体量级,经验公式与方法用于外推和校核,最终以区间形式输入解析模型,才能让矿坑涌水量预测既有物理依据,又能反映勘察阶段固有的不确定性。