3. Theis公式推导与数学结构

Theis公式是承压含水层非稳定流抽水试验分析的基石。它把“抽水过程中水位随时间持续下降”这一现象，转化为一个可以用数学精确描述的解析解。理解这个公式，关键不在于背下表达式，而在于弄清它是从什么物理假设出发、经过怎样的数学变换得到、每个符号在真实含水层中对应什么，以及它在什么条件下才成立。

一、从质量守恒到井函数：推导的物理起点

Theis解建立在三个核心假设之上：含水层为均质、各向同性、等厚的承压含水层；水流服从达西定律；抽水前水位面水平，抽水井为定流量完整井。在这些前提下，取以井轴为中心的极坐标，对含水层微元体建立水量平衡：单位时间内流入微元的水量减去流出的水量，等于含水层弹性释放的水量。由此得到轴对称承压含水层的非稳定流控制方程：

∂²s/∂r² + (1/r)·∂s/∂r = (S/T)·∂s/∂t

其中 s 为降深，r 为距井轴距离，t 为抽水时间，T 为导水系数，S 为储水系数。初始条件为 s(r,0)=0，无穷远处降深为零，井壁处满足定流量边界。Theis通过引入相似变量 u = r²S/(4Tt)，把偏微分方程化为常微分方程，最终得到降深表达式：

s(r,t) = Q/(4πT) · W(u)
u = r²S/(4Tt)
W(u) = ∫ᵤ^∞ (e⁻ˣ/x) dx

这里的 W(u) 就是井函数，它是一个指数积分函数。物理上，W(u) 表示在给定无量纲时间 u 下，含水层已经“贡献”出的无量纲降深总量。u 越小，说明抽水时间越长或距离越近，W(u) 越大，降深越大；u 越大，W(u) 越小，降深趋于零。

二、井函数的级数展开与实用计算

井函数无法用初等函数闭式表达，但可以展开为收敛级数。对积分式作分部积分并逐项展开，得到：

W(u) = -0.5772 - ln u + u - u²/(2·2!) + u³/(3·3!) - u⁴/(4·4!) + …

其中 -0.5772 为欧拉常数。这个级数在 u 较小时收敛很快，通常 u < 0.01 时只取前两项即可获得很高精度，即 W(u) ≈ -0.5772 - ln u。这正是Jacob近似法的数学来源。当 u 较大时，级数收敛变慢，工程上多直接查井函数表或采用数值积分。实际计算中，若 u 在 0.01～0.1 之间，建议保留前三到四项，避免近似误差影响参数反演。

三、降深公式中各参数的物理意义

公式 s = Q/(4πT)·W(u) 中，Q 为抽水井定流量，单位与时间一致；T 为导水系数，反映含水层传递水量的能力，T 越大，同样流量下降深越小、影响范围扩展越快；S 为储水系数，反映含水层释放弹性储水的能力，S 越小，水位下降越迅速、影响传播越远；r 为观测点到井轴的距离；t 为自抽水开始计的时间。需要特别强调符号约定：本文及图示中降深 s 定义为抽水前水位减去抽水后水位，即水位下降为正，这与常规水头定义方向相反，全文保持一致。

四、适用条件与局限性

Theis公式的成立依赖一组严格的理想化条件，实际应用中必须逐条核对，否则参数反演结果会系统性偏离。主要条件包括：含水层为承压、均质、各向同性、等厚且无限延伸；抽水井为完整井且定流量；水流为达西流；储水系数在抽水过程中保持不变；忽略井筒储存效应和井壁水跃。当抽水初期井筒储存占主导时，实测降深曲线会偏离Theis解，此时应剔除早期数据或采用考虑井储的模型。

局限性方面，最突出的是对边界条件的敏感性。若抽水影响范围内存在定水头边界（如河流），降深会提前趋于稳定，Theis解会高估降深；若存在隔水边界，降深会持续增大，Theis解会低估降深。此外，当含水层存在越流补给时，长时间抽水后降深不再按Theis规律增长，需改用Hantush-Jacob等越流模型。对于非完整井、非均质含水层或变流量抽水，Theis解只能作为近似参考。

综上，Theis公式的价值不仅在于给出一个降深数值，更在于它建立了一套“流量—时间—距离—含水层参数”之间的定量关系框架。掌握井函数的级数展开，可以在没有查表条件时快速估算；理解各参数的物理意义，才能在反演结果异常时判断问题出在数据还是模型假设；而清楚适用条件与局限性，则是避免把理想解误用于真实复杂含水层的关键。