Python 稳定流抽水试验承压水完整井计算渗透系数K和影响半径R

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

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

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

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

• 完整计算代码：

# ************************************************************************
# _*_coding:utf-8_*_
# Author: Pandas120 （微信）
# Desription: 稳定流抽水试验承压水完整井求取渗透系数K和影响半径R。
# ************************************************************************
import sympy
from decimal import *

while True:
    K = sympy.Symbol