稳定流抽水试验求取渗透系数K和影响半径R，基本理论依据就是裘布衣公式，下面具体介绍一下在单孔抽水试验中，（无观测孔）只能根据水井的出水量、水位降深等资料，则应消除抽水井附近所产生的三维流、紊流等影响，特别是在抽水时，水位降深很大的时候，应采用消除渗透阻力的方法。

当Q-S关系图，坐标横轴为小时出水量，纵轴为水位降深值，相对应点的连线呈直线时，地下水运动为平面流，承压水完整井采用下列公式：

• K=0.366Q(lgR-lgr)/Ms

式中，Q，涌水量，r井半径，M，承压含水层厚度，s，水位降深值，R影响半径(吉哈尔特提出):

• 完整计算代码：

import sympy
from decimal import *

while True:
    K = sympy.Symbol("K")
    R = sympy.Symbol("R")


    print("分别依次输入 流量或涌水量 Q0[L/s]、含水层厚度 M[m]、降深 S[m]、孔径或过滤管半径 r[m] 等数据！！！")
    Q_0 = input("Q0 = ")
    M = input("M = ")
    S = input("S = ")
    r = input("r = ")
    if Q_0 and M and S and r:
        try:
            Q_0 = float(Q_0) * 86.4
            M = float(M)
            S = float(S)
            r = float(r)
            if r == 0.1:
                r = 0.09999999
            results = sympy.solve([0.366 * Q_0 * sympy.log(R/r,10) / (M*S) - K,10*S*K**0.5-R], [K,R])
            print("K={0},R={1}".format(round(results[1][0],2),round(results[1][1],2)))
        except Exception as e:
            print("数据输入有误！请重新输入。\n")
    else:
        print("数据输入有误！请重新输入。\n")