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 的上下限区间,形成涌水量的区间预测而非单点值。

解析法参数获取:现场试验 → 参数 → 解析计算 现场试验 / 室内试验 单孔稳定流抽水 2~3 个降深,稳定 8 h Q—s 曲线检验井损 非稳定流抽水 + 观测孔 水位—时间全过程记录 含停抽后水位恢复段 注水 / 压水试验 弱富水层、基岩裂隙层 室内给水度试验 原状砂样重力疏干 连通试验 / 流速流向 识别优势径流通道 水文地质参数 渗透系数 K (m/d) 砂砾 10~100|粉砂 0.1~1 导水系数 T = K·M 解析公式的主控参数 (m²/d) 储水系数 S 承压:10⁻⁵~10⁻³ 给水度 μ 无压:砂砾 0.15~0.30 影响半径 R (m) 随降深增大而增大 解析法涌水量预测 大井法 / 裘布依公式 Q = f(K, T, S/μ, R, s) 矿坑尺度降深外推 非稳定流解析解 泰斯 / 博尔顿 / 纽曼 降深随时间发展过程 敏感性分析 S 变幅 3~5 倍 R 取上下限区间 输出:涌水量区间 而非单点值 Q_min ~ Q_max 供疏干设计取值 ×M 经验值仅用于初勘或资料匮乏阶段;实测参数优先,且必须做敏感性分析,避免单点取值误导疏干设计

四、参数取值的工程校核与常见误区

把上述参数代入解析公式之前,建议做三项校核。第一,量纲与厚度一致性: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 给出上下限,不取单点
常见误区:①直接用单井小降深试验的 R 预测矿坑大降深涌水量,系统性偏小;②把承压含水层的 S 误用于无压含水层,导致早期降深预测严重失真;③用经验表 K 值代替实测值却不做敏感性分析;④忽略井损,用抽水井水位直接代入裘布依公式,使 K 偏小、涌水量预测偏保守甚至失真。
实操建议:每个参数至少保留一条独立校核路径——K 用抽水试验与注水/压水试验互校,S 与 μ 用现场非稳定流与室内试验互校,R 用实测外推与经验公式互校。三条路径结果落在同一数量级内,参数方可进入解析计算;否则应补充试验或扩大取值区间。

综上,解析法参数获取的关键在于“试验定值、经验定界、敏感性定区间”。现场试验给出参数的主体量级,经验公式与方法用于外推和校核,最终以区间形式输入解析模型,才能让矿坑涌水量预测既有物理依据,又能反映勘察阶段固有的不确定性。