3. 降深叠加原理的数学基础

井群干扰抽水试验中，观测孔所测得的降深并不是某一口井单独作用的结果，而是多口井同时抽水时各自产生降深的叠加。要把这一直觉变成可计算、可校核的工具，就必须先弄清叠加原理的数学前提：它本质上是把描述地下水运动的偏微分方程当作线性方程来处理，从而允许把多个解直接相加。承压含水层与潜水含水层在这一点的适用性差别很大，这正是本节要解决的核心问题。

一、叠加原理成立的数学条件

承压含水层中，若假定含水层均质、等厚、各向同性，且水头变化引起的水量释放由弹性释水系数 μ_e 描述，则定流量抽水满足的方程为：

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

其中 s 为降深，T 为导水系数，S 为储水系数。该方程对 s 是线性的：若 s₁、s₂ 分别是两口井单独抽水时的解，则 s₁+s₂ 仍是同一方程的解。边界条件（定流量、定水头）与初始条件（初始降深为零）也都是线性的，因此叠加在数学上严格成立。这就是井群干扰计算可以“逐井求降深、再代数相加”的根本依据。

需要强调的是，线性成立依赖三个前提：含水层参数不随水头变化、抽水引起的渗流仍满足达西定律、各井的边界条件互不改变。一旦其中任一条件被破坏，叠加就只能作为近似。

二、承压含水层中的叠加表达

设井群有 n 口井，第 i 口井以流量 Q_i 抽水，观测点 P 到该井的距离为 r_i，抽水历时为 t。按泰斯公式，P 点由第 i 口井引起的降深为：

s_i = (Q_i / (4πT)) · W(u_i)，  u_i = r_i²·μe / (4T·t)

则 P 点的总降深为各井贡献之和：

s_P = Σ s_i = (1/(4πT)) · Σ [ Q_i · W(u_i) ]

若各井流量相同为 Q、且观测点距各井距离相近，可近似写成 s_P ≈ n·Q·W(u)/(4πT)，这正是“干扰使降深增大”的定量来源。实际计算时，通常先对每口井单独查 W(u) 表或计算井函数，再逐项相加，最后与观测孔实测降深对比，用差值反推参数或检验模型。

三、潜水含水层的适用性差异

潜水含水层不能直接套用上述线性叠加。原因在于潜水含水层的导水系数 T = K·h 随饱和厚度 h 变化，而 h 本身又随降深改变，方程中出现 h·∂h/∂r 这类非线性项。工程上常用的处理办法是降深代换法：令 H₀ 为初始水位，则把潜水问题转化为承压问题：

s' = (H₀² - h²) / (2H₀)

在降深相对于饱和厚度较小（一般 s/H₀ 小于 0.1～0.2）时，代换后的 s' 近似满足线性方程，此时可对 s' 逐井叠加，再反算实际水位 h。若降深较大，非线性显著，叠加结果会系统性偏小或偏大，必须改用非线性数值模型或分段线性化处理。因此，“潜水能否叠加”不是绝对的，而取决于降深与含水层厚度的比值。

四、物理意义：干扰不是“抵消”而是“累加”

叠加的物理含义是：每口井抽水都在含水层中形成一个独立的降落漏斗，多个漏斗在空间上相互嵌套、在数值上相加。因此井群中心的降深总是大于任一单井在该点的降深，井距越近、流量越大、抽水时间越长，叠加后的降深越大。这也解释了干扰抽水试验中“总降深大于各单井降深之和的观测值”往往意味着参数估计偏大或存在越流、补给等额外因素，需要用叠加残差去排查。

五、示例：三口井干扰的叠加计算

假设某承压含水层 T = 200 m²/d，μ_e = 2×10⁻⁴，三口井流量均为 500 m³/d，抽水 1 d，观测孔到三口井距离分别为 30 m、50 m、80 m。按 u_i = r_i²μ_e/(4Tt) 计算得 u₁≈3.4×10⁻⁵、u₂≈9.4×10⁻⁵、u₃≈2.4×10⁻⁴，对应 W(u) 分别约为 9.6、8.6、7.5（示例取值，实际应查表）。则各井降深分量约为 1.91 m、1.71 m、1.49 m，总降深约 5.11 m。这一结果说明：即使最远井的贡献也有近 1.5 m，井群干扰不可忽略。上述数值为教学示例假设，非实测结果。

掌握上述数学基础后，后续井群干扰资料整理时，就能明确每一步“为什么可以相加、在什么条件下相加、相加后如何还原”，而不是机械套用公式。