地下水向完整井的非稳定运动-吉林大学.PPT

第四章 地下水向完整井的非稳定运动 第四章 地下水向完整井的非稳定运动 §4-1 承压含水层中的完整井流 §4-2 有越流补给的完整井流 §4-3 有弱透水层弹性释水补给和越流补给的完整井流 §4-4 潜水完整井流 天地不可一日无和气， 人心不可一日无喜神。 著名科学家 §4-1承压含水层中的完整井流 当承压含水层侧向边界离井很远，边界对研究区的水头分布没有明显影响时，可以把它看作是无外界补给的无限含水层。 4.1.1 定流量抽水时的Theis公式 承压含水层中单井定流量抽水的数学模型是在下列假设条件下建立的： (1) 含水层均质各向同性，等厚，侧向无限延伸，产状水平； (2) 抽水前天然状态下水力坡度为零； (3) 完整井定流量抽水，井径无限小； (4) 含水层中水流服从Darcy定律； (5) 水头下降引起的地下水从贮存量中的释放是瞬时完成的。 The Theis solution assumes the following: The aquifer is confined and has an apparent infinite extent; The aquifer is homogeneous, isotropic, and of uniform thickness over the area influenced by pumping; The piezometric surface was horizontal prior to pumping; The well is fully penetrating and pumped at a constant rate; The well diameter is small, so well storage is negligible； Water removed from storage is discharged instantaneously with a decline in head。 Data requirements: Drawdown vs. time at an observation well Finite distance from the pumping well to observation well Pumping rate (constant) 在上述假设条件下，抽水后将形成以井轴为对称轴的下降漏斗，将坐标原点放在含水层底板抽水井的井轴处，井轴为Z轴，如图4-1所示。 单井定流量的承压完整井流，有如下的数学模型： 式中，s=H0-H。 下边研究如何求降深函数s (r, t)。为此，利用Hankel变换，将方程式(4-1)两端同乘以rJ0(βr)，并在(0,∞)内对r积分。 设导压系数 ，则有： 方程式右端 方程式左端，利用分部积分，同时注意到边界条件式 （4-3）与式(4-4)，有： 按Bessel函数的性质，有： 因此，有： 上述定解问题，经过Hankel变换，消去了变量r，转变为常微分方程的初值问题，即： 其解为： 再通过Hankel逆变换由 求s，即： 先计算方括号内的积分，为此设： 将(4-6)式对r求导数，有： 根据（4-6）式，有： 两边积分得： , 令C1=lnC ，则有： 故： 利用r=0时的F(r) 值，由（4-6）可以确定C值： 由(4-7)式，有： 把上式代入(4-5)式，有： 为计算方便，对（4-8）式进行变量代换，令： 同时更换积分上下限，当τ=0 时， ; 当τ=t时，y=∞ , 于是， 其中， (4-9)式为无补给的承压水完整井定流量非稳定流计算公 式，也就是著名的Theis公式。 在地下水动力学中，采用井函数W(u)代替(4-9)式中的指数 积分式： 则(4-9)式可改写成： 式中，s——抽水影响范围内，任一点任一时刻的水位降深(m)； Q——抽水井的流量(m3/d, m3/h)； T——导水系数(m2