2. 渗透系数概念与现场试验物理基础

渗透系数是水文地质试验与计算中的核心参数,它把“水在含水层中流动的难易程度”转化为一个可以测量、可以代入公式的量。理解渗透系数,必须回到达西定律:v = K·i,其中 v 为比流量(单位时间通过单位过水断面的流量,量纲为 L/T),i 为水力坡度(无量纲),K 即渗透系数,量纲为 L/T,常用 m/d 或 cm/s。这个看似简单的线性关系,是现场抽水试验、注水试验、压水试验计算渗透系数的共同物理基础。

需要特别强调,达西定律中的 v 是比流量而非孔隙中的真实流速。含水层中实际过水断面只是岩土孔隙部分,真实平均流速 v' = v/n(n 为有效孔隙度),因此 v' 大于 v。现场试验测得的是宏观平均渗透性能,计算中使用的断面面积应为整个过水断面,而不是孔隙面积,这一点在整理试验数据时极易混淆。

达西定律的适用条件

达西定律并非在任何条件下都成立,现场试验方案设计与成果解释必须核查以下条件:

  • 层流条件:水流雷诺数较低,惯性力可忽略。粗颗粒地层(如卵砾石)中抽水降深过大时,井壁附近流速高,可能出现非线性流,此时按达西公式计算会低估渗透系数。
  • 均质各向同性或可等效处理:公式推导通常假定含水层均质、等厚、各向同性;实际层状或各向异性含水层需按等效或分层处理。
  • 稳定渗流场:抽水试验中要求流量与水位降深达到稳定,即 Q 与 s 不随时间变化,才能用稳定流公式(如裘布依公式)计算。
  • 水头损失仅为黏滞阻力:忽略紊流、惯性及化学堵塞等附加阻力。
  • 含水层饱和、水不可压缩:非饱和带或承压水头剧烈变化时需另行处理。

工程上常用一个经验判据辅助判断:当井中水位降深与含水层厚度之比过大、或抽水流量使井壁附近流速超过临界值时,应怀疑达西定律适用性,可通过多降深试验绘制 Q–s 曲线检验其线性。若 Q–s 呈明显上凸或下凹,说明存在非线性流或井损,需修正后再计算。

渗透系数的物理意义与量纲

渗透系数 K 综合反映了岩土孔隙结构(孔隙大小、连通性、颗粒级配)与流体性质(密度、黏度)对渗流的影响。从物理本质看,它可分解为渗透率 k(仅与介质结构有关,量纲 L²)与流体性质项:K = k·ρg/μ。因此同一地层,抽水(水)与抽气、抽油得到的“渗透系数”数值不同,现场试验报告必须注明流体种类与温度。

量纲上,K 为长度/时间,与流速同量纲,这提示我们:渗透系数在数值上等于水力坡度为 1 时的比流量。常见量纲换算关系为 1 m/d ≈ 1.157×10⁻³ cm/s,1 cm/s ≈ 864 m/d。实际资料整理中,黏土类弱透水层 K 常在 10⁻⁴~10⁻² m/d,粉细砂约 10⁻²~1 m/d,中粗砂及砾石可达 1~10² m/d,量级差异极大,因此试验方法与计算公式必须与地层渗透性匹配。

现场记录中应同时记录水温。水的运动黏度随温度变化,20℃ 与 10℃ 时同一天然介质的渗透系数可相差约 20%~30%,成果对比时宜统一换算到标准温度(如 10℃ 或 20℃)。

稳定流与非稳定流的基本假设

现场试验按水位随时间的变化分为稳定流与非稳定流两类,二者的假设和适用公式不同,下图概括了从达西定律到两类试验计算路径的物理基础关系。

达西定律 v = K·i 层流 · 饱和 · 均质等效 · 黏滞阻力主导 K = k·ρg/μ(量纲 L/T) 稳定流试验 Q、s 不随时间变化 假设:补给充分、降落漏斗定型 代表公式:裘布依 Q=2πKHM/(lnR/r) 非稳定流试验 Q 恒定、s 随时间下降 假设:弹性释水、无越流补给 代表方法:泰斯 s–lgt 直线段 现场判稳依据 连续 2~4 h 内 s 波动 ≤ 1~2 cm Q 波动 ≤ 5%,方可取稳定值 现场判据 s–lgt 曲线出现直线段 后期偏离提示边界补给或越流 注:两类试验均以达西层流为前提,非线性流或井损显著时须先修正

稳定流假设的核心是“补给与排泄达到动态平衡”,降落漏斗不再扩展,抽水量全部来自补给。裘布依公式 Q = 2πKHM/(ln(R/r)) 中,H 为含水层厚度、M 为水位降深、R 为影响半径、r 为井半径,计算时需先由观测孔或经验公式确定 R。现场判稳通常要求连续 2~4 小时水位波动不超过 1~2 cm、流量波动不超过 5%,取平均值代入。

非稳定流假设则承认含水层释放水量使水位持续下降,抽水量来自含水层的弹性释水(承压)或重力释水(潜水)。泰斯公式以 s–lgt 曲线的直线段斜率求导水系数 T,再由 T = K·M 反算 K。其关键假设包括:含水层无限延伸、等厚均质、无越流补给、井径无限小、释水瞬时完成。实际试验中若曲线后期上翘或下弯,往往说明存在边界补给、越流或井损,需要调整模型而非硬套公式。

现场操作中的关键控制点

以一个承压含水层单孔稳定流抽水试验为例(示例参数为教学假设):含水层厚度 M = 20 m,抽水流量 Q = 30 m³/h,稳定降深 s = 5 m,影响半径按经验取 R = 150 m,井半径 r = 0.1 m。代入裘布依公式:

K = Q·ln(R/r) / (2π·M·s)
  = (30/24)·ln(1500) / (2π×20×5)
  ≈ 1.25×7.313 / 628.3
  ≈ 0.0145 m/d

该结果属弱透水—中等透水区间,与粉细砂量级相符,可作为示例校核。实际操作中需注意:流量计与水位计应同步读数;动水位应在抽水稳定后每 30 min 记录一次;恢复水位阶段也应记录,用于检验参数一致性。若多降深试验的 Q–s 关系明显偏离直线,应先判断是井损还是非线性流,再决定是否修正公式。

不要把稳定流与非稳定流公式混用:稳定流公式要求 s 已稳定,非稳定流公式要求有随时间变化的水位序列。用未稳定的数据套裘布依公式,会系统性高估或低估 K,这是现场成果审查中最常见的错误之一。

综上,渗透系数的计算不是套公式,而是先核查达西定律适用条件、明确 K 的物理含义与量纲、再根据试验是否达到稳定选择对应假设与公式。只有物理基础清楚,后续的现场布置、数据整理与成果合理性判断才有可靠依据。