为一阶第二类虚宗量Bessel函数.ppt

第三章 地下水向完整井的稳定运动 主要内容 第三章 地下水向完整井的稳定运动 §3.1 水井的分类及井流特征 §3.2 地下水向承压水井和潜水井的稳定流动 §3.3 越流含水层中地下水向承压井的稳定流动 §3.4 流量和水位降深关系的经验公式 §3.5 地下水向干扰井群的稳定运动 §3.6 均匀流中的井 §3.7 井损与有效井径的确定方法 §3.3 越流含水层中地下水向承压井的稳定流动 3.3.1 地下水运动的数学方程及其解 图3-9表示有越流补给时，在无限承压含水层中的一口完整井。 因从井中抽水，造成水头降低，和相邻含水层(图中为潜水含水层)之间产生水头差或将原有的水头差扩大，相邻含水层中的水通过弱透水层越流补给抽水含水层。当抽水延续一定时间后，进入抽水含水层降落漏斗范围内的越流量和抽水量平衡时，水流达到稳定状态。 此时假设: 发生越流的潜水含水层,有足够的补给量维持初始水位不变;弱透水层的弹性释放量很小,可以忽略不计,且流向井的水流基本上仍保持水平流动。 在此假设条件下,抽水含水层内的水头满足方程（1-84）。对于稳定流动,与该方程相应的以柱坐标表示的方程为： （3-29） 把水头改用降深表示,令H0- H = s,并代入上式,经变量代换后得: （3-30） 相应的边界条件为: s=0， 当r→∞时; 当r=rw时 方程(3-30)是零阶虚宗量Bessel方程。其通解为： （3-31） 式中，α，β为待定系数， 分别为零阶第一类和第二类虚宗量Bessel函数。 由边界条件知，当r→∞时， =0，而 ≠0，把它们代入（3-31）式可得a=0。因而有： 再考虑井壁边界条件： 得： 最后得： （3-32） 式中的 K1(rw/B) 为一阶第二类虚宗量Bessel函数。 在一般越流含水层中,越流因素B都有相当大的值,故实际上rw/B1。对于Bessel函数函数,当x1时，xK1(x)=1(如当x0.02时,误差小于1%)。因此，（3-32）式可简化为 （3-33） 式(3-33)式称为Hantush-Jacob 公式。虚宗量Bessel函数有表可查。 在抽水井附近，rw/B1。对于第二类虚宗量Bessel函数，当x1时，xK1(x) =ln(1.123/x)。故(3-33)式又可简化为: （3-34） 采用此式计算，误差不大。当r/B0.35时，误差小于5%；当r/B0.1时，误差小于1%。 越流量占抽水量的比例: 设Qr代表径向距离r处的侧向流入量，则有 井的流量为： 取二式比值，并注意当rw/B《1时， ，因而有： （3-35） 表明侧向流入量占抽水井流量的比例，与径向距离r与越流因素B的比值有关。 图3-10 Qr/Q与r/B关系曲线（据J. Bear） 图3-10为Qr/Q与r/B的关系曲线。由图可以看出,当r=4B时，Qr/Q=0.05，表示侧向流入量（来自该断面到无穷远处的越流量）只占抽水井流量的5%，而95%的抽水井流量是来自r4B地段的越流量。 3.3.2 含水层参数计算方法 可利用式（3-33）和式（3-34）,根据稳定流抽水试验资料求参数。此时,要求有距抽水井不同距离r的若干个观测孔。测得各观测孔的水位降深后,可用下述方法求导水系数T，越流因素B和