4. 单井抽水试验降深计算回顾

井群干扰抽水试验的降深叠加计算,本质上是在单井降深计算的基础上做线性叠加。因此在进入井群问题之前,必须先把单井条件下"给定流量求降深"和"给定降深求流量"这两类计算彻底理清。本节回顾裘布依公式与泰斯公式的适用条件、计算路径和结果判读方式,为后续叠加计算奠定基础。

一、稳定流:裘布依公式的计算逻辑

裘布依公式描述的是承压完整井在稳定流条件下的降深分布。其基本形式为:

s = Q / (2πKM) · ln(R / r)
Q = 2πKM · s / ln(R / r)

式中:s 为抽水井降深(m);Q 为抽水流量(m³/d);K 为含水层渗透系数(m/d);M 为承压含水层厚度(m);R 为影响半径(m);r 为计算点到井轴的距离(m)。对于潜水完整井,需用 (H² − h²) 替代 2Ms,即 Q = πK(H² − h²) / ln(R/r),其中 H 为初始水位至隔水底板的厚度,h 为抽水后井中水深。

裘布依公式的关键特征是:降深与流量成正比,与距离呈对数关系。这意味着距离井越远,降深衰减越慢——这正是井群干扰能够波及较大范围的根本原因。实际计算中,影响半径 R 通常由经验公式或抽水试验实测确定,例如承压井常用 R = 10s√K 的经验式(示例性经验关系,具体系数应结合场地资料选取)。

裘布依公式要求含水层等厚、各向同性、完整井、稳定流。若抽水尚未达到稳定,或含水层边界明显,直接套用会产生系统偏差,此时应改用泰斯公式或考虑边界条件的解析解。

二、非稳定流:泰斯公式的计算逻辑

泰斯公式描述承压含水层中定流量抽水时的非稳定降深过程:

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

式中:T 为导水系数(m²/d),S 为储水系数(无量纲),t 为自抽水开始起算的时间(d),W(u) 为井函数。当 u 较小时(一般 u < 0.01),可用雅可比近似 W(u) ≈ ln(2.25Tt / r²S),此时降深与时间的对数近似呈直线关系,这也是配线法和直线图解法求参的基础。

泰斯公式与裘布依公式的关系值得注意:当抽水时间足够长、影响范围趋于稳定时,泰斯解在形式上可退化为对数形式,与裘布依公式结构一致。因此在井群叠加计算中,若采用稳定流叠加,本质上是把每个单井的裘布依降深逐点相加;若采用非稳定流叠加,则把每个单井在对应时刻的泰斯降深相加。

三、单井降深计算的完整操作路径

以某承压完整井为例(以下参数为教学示例假设值):含水层厚度 M = 20 m,渗透系数 K = 5 m/d,抽水流量 Q = 500 m³/d,影响半径 R = 300 m。计算距井轴 r = 50 m 处的稳定降深:

T = K·M = 5 × 20 = 100 m²/d
s = 500 / (2π × 100) × ln(300 / 50)
  = 0.7958 × ln(6)
  = 0.7958 × 1.7918
  ≈ 1.43 m

若改用泰斯公式估算同一位置在 t = 2 d 时的降深,取 S = 0.001:

u = 50² × 0.001 / (4 × 100 × 2) = 2.5 / 800 ≈ 0.003125
u < 0.01,用雅可比近似:
s ≈ 500 / (4π × 100) × ln(2.25 × 100 × 2 / (50² × 0.001))
  ≈ 0.3979 × ln(18)
  ≈ 0.3979 × 2.8904
  ≈ 1.15 m

两个结果量级一致但数值不同,原因在于稳定流假设影响范围已达 R,而非稳定流在 2 d 时影响范围尚未扩展到 300 m。这一差异提醒我们:叠加计算前必须明确采用稳定流还是非稳定流框架,两者不能混用。

四、单井计算与井群叠加的衔接关系

单井计算得到的是"一个井在任意点产生的降深"这一基本单元。井群干扰抽水时,若各井同时以定流量抽水,且含水层满足线性叠加条件(承压含水层、降深远小于含水层厚度、无越流或边界影响可忽略),则任意点的总降深等于各单井在该点产生降深的代数和:

s_total(P) = Σ s_i(P)

其中每个 s_i(P) 都用本节回顾的单井公式计算,只是把距离 r 换成第 i 口井到计算点 P 的距离 r_i。因此单井计算的准确性、参数取值的一致性(同一 K、M、S、R)直接决定叠加结果的可靠性。

单井降深计算 → 井群叠加的衔接关系 单井降深计算(基本单元) 稳定流:裘布依公式 s = Q/(2πKM)·ln(R/r) 适用:完整井、等厚、稳定流 参数:Q、K、M、R、r 特征:s 与 Q 成正比    s 与 ln r 呈线性衰减 非稳定流:泰斯公式 s(r,t) = Q/(4πT)·W(u) u = r²S/(4Tt) 参数:Q、T、S、r、t u<0.01 时用雅可比近似 特征:s 随时间对数增长 输出:单井在计算点 P 的降深 s_i(P) (距离 r 换成第 i 井到 P 的距离 r_i) 逐点 求和 井群降深叠加 叠加前提 承压含水层、线性叠加成立 各井参数取值一致(K、M、S、R) 叠加公式 s_total(P) = Σ s_i(P) r_i = 第 i 井到 P 的距离 结果判读 总降深大于任一单井降深 井距越小,干扰叠加越显著 框架一致性原则 稳定流与泰斯解不可混用叠加 注:图中参数与公式为教学示例,实际取值应依据场地抽水试验资料确定

五、回顾中容易出错的几个细节

  • 单位统一:Q 用 m³/d、K 用 m/d、M 用 m 时,T 的单位为 m²/d,降深为 m。若 Q 用 L/s,需先乘 86.4 换算为 m³/d,否则降深会差近两个数量级。
  • r 的取值:计算抽水井自身降深时,r 应取井半径 r_w(并考虑井损);计算观测孔降深时取观测孔到井轴距离。叠加计算中每个 r_i 都必须按计算点重新量取,不能沿用单井试验中的固定距离。
  • 稳定流与非稳定流不可混用:若井群抽水尚未稳定,用裘布依公式叠加会低估早期降深;若已稳定,用泰斯公式叠加则需取足够长的 t。框架一旦选定,所有单井项必须同框架计算。
  • 降深参照一致:所有降深均以同一初始水位为基准,正号表示水位下降。若某井为回灌井,其"降深"为负值,叠加时按代数和处理。
叠加法成立的前提是含水层对降深的响应是线性的。当降深接近或超过承压含水层顶板、含水层转为无压流时,线性叠加不再严格成立,此时应改用无压流公式或数值方法,不能简单相加。