穩(wěn)定流抽水試驗(yàn)求取滲透系數(shù)K和影響半徑R呀邢，基本理論依據(jù)就是裘布衣公式，下面具體介紹一下在單孔抽水試驗(yàn)中，（無(wú)觀測(cè)孔）只能根據(jù)水井的出水量诵叁、水位降深等資料，則應(yīng)消除抽水井附近所產(chǎn)生的三維流钦椭、紊流等影響拧额，特別是在抽水時(shí)碑诉，水位降深很大的時(shí)候，應(yīng)采用消除滲透阻力的方法侥锦。
當(dāng)Q-S關(guān)系圖联贩，坐標(biāo)橫軸為小時(shí)出水量，縱軸為水位降深值捎拯，相對(duì)應(yīng)點(diǎn)的連線呈直線時(shí)泪幌，地下水運(yùn)動(dòng)為平面流，承壓水完整井采用下列公式：
- K=0.366Q(lgR-lgr)/Ms
式中署照，Q祸泪，涌水量，r井半徑建芙，M没隘，承壓含水層厚度，s禁荸，水位降深值右蒲，R影響半徑(吉哈爾特提出):
- 完整計(jì)算代碼：
import sympy
from decimal import *
while True:
K = sympy.Symbol("K")
R = sympy.Symbol("R")
print("分別依次輸入 流量或涌水量 Q0[L/s]、含水層厚度 M[m]赶熟、降深 S[m]瑰妄、孔徑或過(guò)濾管半徑 r[m] 等數(shù)據(jù)！Ｓ匙间坐！")
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("數(shù)據(jù)輸入有誤！請(qǐng)重新輸入邑退。\n")
else:
print("數(shù)據(jù)輸入有誤竹宋！請(qǐng)重新輸入。\n")
