2. 承压含水层抽水物理过程解析
要正确使用Theis公式，先要弄清抽水时承压含水层里到底发生了什么。承压含水层与潜水含水层的根本区别在于：含水层顶板被一层不透水的隔水层（如黏土、页岩）覆盖，含水层内部的水承受着上覆地层与静水压力，处于承压状态。井未抽水时，测压水位高于含水层顶板；抽水后，井中水位下降，但含水层本身并未被疏干，只是水头（测压水位）降低。这一区别决定了承压含水层的释水机制与潜水含水层完全不同。
一、承压含水层概念模型
把实际地层抽象为概念模型，是建立Theis方程的前提。一个典型的承压含水层概念模型包含以下要素：
- 含水层：厚度为 M（m）的均质、各向同性水平层，渗透系数为 K（m/d），导水系数 T = K·M（m²/d）。
- 隔水顶板与底板：上下均为不透水边界，抽水过程中无越流补给，含水层与外界无垂向水量交换。
- 完整抽水井：井贯穿整个含水层厚度，以定流量 Q（m³/d）抽水，井径相对于影响范围可忽略（视为线汇）。
- 初始水头：抽水前水头面水平，初始水头 H₀ 处处相等，降深 s = H₀ − H(r,t) 初始为零。
- 侧向边界：含水层在水平方向无限延伸，抽水的影响尚未到达边界。
这一模型把复杂的层状地层简化为“等厚、均质、无限、无越流”的理想承压系统，后续所有解析解都建立在此假设之上。实际工作中，若含水层厚度变化大、存在明显越流或边界较近，就需要修正模型或改用其他方法。
二、抽水引起的降深时空演化
抽水开始后，井周围水头下降，形成以井为中心的降落漏斗。降深 s 随径向距离 r 和时间 t 变化，其演化具有明显的阶段性：
- 早期：抽水初期，水量主要来自井附近含水层的弹性释放，降深快速增大，影响范围很小，漏斗陡而窄。
- 中期：随着抽水持续，降落漏斗不断向外扩展，降深增速逐渐变缓，漏斗形态趋于平缓。
- 后期：在无限含水层中，降深随时间近似按对数规律缓慢增大，漏斗持续但缓慢地扩展；若存在定水头边界或越流补给，则降深最终趋于稳定。
从空间上看，同一时刻降深随 r 增大而减小，井壁处降深最大，远处趋近于零。Theis解正是用井函数 W(u) 描述这种“时间—距离”双重衰减关系，其中 u = r²S/(4Tt)，S 为储水系数。理解 u 的物理含义很关键：u 越小，表示该点受抽水影响越显著；u 越大，影响越微弱。
三、弹性释水机制
承压含水层抽水时，含水层并未被疏干，水量从何而来？答案是弹性释水。当水头下降时，含水层内孔隙水压力减小，有效应力增大，导致两方面的压缩与膨胀：
- 水的膨胀：孔隙水压力降低，水体体积略微膨胀，释放出少量水。
- 骨架压缩：有效应力增大使含水层骨架被压缩，孔隙体积减小，把孔隙中的水“挤”出来。
两者共同作用的结果，用储水系数 S（无量纲）表示：单位面积、单位水头降深下，从整个含水层厚度中释放的水量。承压含水层的 S 通常很小，约为 10⁻⁵～10⁻³ 量级，远小于潜水含水层的给水度（0.1～0.3）。这正是承压含水层抽水时降深发展快、影响范围大、水位恢复慢的物理根源。工程上常用储水率 Sₛ = S/M（1/m）描述单位厚度含水层的释水能力。
四、Theis模型的假设条件
Theis公式是承压含水层非稳定流的标准解析解，但其成立依赖一组严格的假设。实际计算前必须逐条核对，否则结果可能严重偏离：
| 假设条件 | 物理含义 | 不满足时的后果 |
|---|---|---|
| 含水层均质、各向同性、等厚 | T、S 为常数 | 降深曲线畸变，参数反演失真 |
| 含水层水平、无限延伸 | 无侧向边界影响 | 边界附近降深偏大或偏小 |
| 上下隔水层不透水 | 无越流补给 | 后期降深被高估 |
| 完整井、定流量抽水 | 井贯穿含水层，Q 恒定 | 井附近流态三维化，解不适用 |
| 水头下降即瞬时释水 | 忽略释水滞后 | 早期段拟合偏差 |
| 井径无限小 | 视为线汇 | 井壁附近误差增大 |
下面用一张图把“概念模型—抽水响应—弹性释水—假设条件”四者的关系串起来，帮助建立整体物理图像：
五、从物理过程到计算准备
把上述物理图像落到实际计算上，需要先完成三件事。第一，确定含水层参数：通过抽水试验或经验值获得 T 与 S，例如某承压含水层 K = 8 m/d、M = 25 m，则 T = 200 m²/d；S 取 5×10⁻⁴。第二，明确抽水条件：定流量 Q = 1000 m³/d，观测井距抽水井 r = 50 m。第三，核对假设：若含水层在 500 m 外有河流定水头边界，则Theis解在后期会高估降深，需改用叠加法或镜像法处理。
以示例参数估算：当 t = 1 d 时，u = r²S/(4Tt) = 50²×5×10⁻⁴/(4×200×1) ≈ 0.00156，此时 u 很小，W(u) ≈ −0.5772 − ln u ≈ 5.88，降深 s = Q·W(u)/(4πT) = 1000×5.88/(4π×200) ≈ 2.34 m。当 t 增至 10 d，u 减小到约 1.56×10⁻⁴，W(u) ≈ 8.18，s ≈ 3.25 m。可见降深随时间增大但增速递减，这正是承压含水层弹性释水 S 极小、漏斗缓慢扩展的定量体现。上述数值为教学示例假设，实际工程应以现场试验参数为准。
理解了概念模型、降深演化、弹性释水与假设条件这四层内容，就为下一步套用Theis公式、查算井函数 W(u) 和反求参数打下了物理基础。计算时每一个符号背后都对应着具体的含水层行为，而不是纯粹的数学代入。