1. 课程导入与实战目标

承压含水层抽水试验的降深预测，是水文地质计算中最经典也最容易出错的一类问题。说它经典，是因为Theis在1935年提出的非稳定流公式至今仍是承压含水层参数反演与降深预报的理论基石；说它容易出错，是因为实际计算中涉及单位换算、井函数查算、参数反演、叠加原理等多个环节，任何一处疏忽都会让结果偏离数倍。本课程不打算停留在公式推导层面，而是围绕"给定一套抽水试验数据，如何算出降深、如何反求参数、如何判断结果是否可信"这条主线，把Theis公式真正用起来。

为什么Theis公式是承压含水层计算的核心

承压含水层与潜水含水层的本质区别在于：承压含水层的释水来自含水层骨架压缩与水的弹性膨胀，释水系数S通常在10⁻⁵～10⁻³量级，远小于潜水的给水度；同时承压含水层顶板隔水，抽水时不存在重力疏干过程。Theis公式正是建立在这一弹性释水机理之上，其降深表达式为：

s(r,t) = Q / (4πT) · W(u)
u = r²S / (4Tt)
W(u) = ∫ᵤ^∞ (e^(-x)/x) dx   （井函数）

式中：s为降深（m）；Q为抽水流量（m³/d）；T为导水系数（m²/d）；S为储水系数（无量纲）；r为观测孔到抽水井的距离（m）；t为抽水时间（d）；W(u)为井函数。这个公式把"抽多少水、抽多久、在多远处、降多深"四件事用一个井函数串了起来，是后续所有配线法、直线法、叠加计算的原点。

本课程的实战任务总览

整门课程围绕一个贯穿始终的实战任务展开：某承压含水层场地开展单井定流量抽水试验，主井以Q=1200 m³/d抽水，距主井r=50 m处布设观测孔，测得不同时刻的降深序列；要求利用Theis理论完成参数反演，并预测抽水30天时不同距离处的降深。这是一个典型的"正演—反演—再正演"闭环任务，示例参数均为教学假设，不代表任何实际场地实测值。

围绕这个任务，课程将依次解决以下问题：

• 正演计算：已知T、S，如何用井函数算出任意时刻、任意距离的降深；井函数如何查算或近似。
• 参数反演：已知抽水试验的s–t数据，如何用配线法、Jacob直线法反求T和S，两种方法各自适用条件是什么。
• 边界与叠加：遇到隔水边界、补给边界或多井干扰时，如何用镜像法与叠加原理扩展Theis公式。
• 结果校核：如何判断反演参数是否合理，降深预测的误差来源在哪里，如何做敏感性分析。

学习路径与成果要求

课程的学习路径遵循"先建立物理图像，再掌握计算工具，最后完成综合任务"的顺序。前段重点是把Theis公式的每个符号与现场对象对应起来——抽水井、观测孔、隔水顶底板、弹性释水；中段进入计算实操，包括井函数表的读取、W(u)–u曲线的使用、Jacob法的半对数直线拟合；后段则把单井计算扩展到边界条件和多井系统，并完成贯穿任务的完整计算报告。

成果要求方面，学完本课程你应当能够独立完成三件事：其一，给定Q、T、S、r、t，手算或借助表格算出降深，并说明井函数近似式的误差范围；其二，给定一组抽水试验降深数据，用至少两种方法反求T、S并比较结果；其三，针对含边界或多井的场地，建立叠加计算方案并给出预测降深。这三项能力对应课程的三类核心练习，也是最终综合任务的评分依据。

需要特别强调的是，Theis公式的适用前提是：含水层均质、等厚、无限延伸，抽水流量恒定，井径无限小（完整井），释水瞬时完成。实际场地很难完全满足这些条件，因此课程在讲解计算的同时，会反复提醒你关注适用性判断——比如抽水初期井储效应显著时数据不可用，抽水时间足够长（u<0.01）时才能用Jacob近似，隔水边界附近降深会偏大等。这些判断不是附加知识，而是决定计算结果能不能用的关键。