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)直接决定叠加结果的可靠性。
五、回顾中容易出错的几个细节
- 单位统一: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。框架一旦选定,所有单井项必须同框架计算。
- 降深参照一致:所有降深均以同一初始水位为基准,正号表示水位下降。若某井为回灌井,其"降深"为负值,叠加时按代数和处理。