地下水动力学-第五讲课件.pptx

地下水向完整井的非稳定运动（续第四章）三、考虑弱透水层弹性释水补给和越流补给的完整井流前述在研究越流补给时忽略了弱透水含水层的弹性释水补给，当弱透水层较厚时，这种补给量是可观的，不能忽略不计。（一）模型1、物理模型、基本假设与数学模型（1）物理模型1960年，M. S. Hantush提出的三层结构模型，根据弱透水层弹性释水与相邻含水层关系分三种情况进行了研究，如地下水动力学_第四章4-16.JPG图4-16所示。1）与两弱透水相邻的越流含水层为定水头含水层；2）与两弱透水层相邻为隔水层；3）上弱透水含水层与定水头含水层相邻，下弱透水含水层与隔水层相邻。地下水动力学讲稿_第五讲地下水向完整井的非稳定运动（续第四章）（2）基本假设1）抽水含水层与弱透水含水层均质，各向同性且等厚（K1、K、K2，μ*1、μ*、μ*2，m1，M，m2）。2）含水层产状水平，无限分布；天然水力坡度为零。3）单井定流量抽水。4）抽水含水层抽水时能得到弱透水含水层弹性释水的补给；在弱透水层中水流垂直运动，抽水含水层水流水平运动。5）与弱透水层相邻条件： a) 定水头含水层； b) 隔水含水层； c) 定水头含水层与各水含水层。地下水动力学讲稿_第五讲初始条件：a）定水头边界条件与越流含水层相邻顶板与主含水层相邻底板b）隔水边界条件与越流含水层相邻顶板与主含水层相邻底板地下水向完整井的非稳定运动（续第四章）（3）数学模型1）上弱透水含水层（仅有垂向运动）地下水动力学讲稿_第五讲a）定水头边界条件与主含水层相邻顶板与越流含水层相邻底板b）隔水边界条件与主含水层相邻顶板与越流含水层相邻底板地下水向完整井的非稳定运动（续第四章）2）下弱透水含水层（仅有垂向运动）初始条件：地下水动力学讲稿_第五讲初始条件：边界条件地下水向完整井的非稳定运动（续第四章）3）抽水（主）含水层（仅有水平运动） 式中“正负”号表示上下弱透水层向主抽水层补给弹性释水和越流补给。地下水动力学讲稿_第五讲地下水向完整井的非稳定运动（续第四章）（二）数学模型的求解1、求解思路上述数学模型为偏微分方程组，求解思路：第一步进行Laplace变换，化为常微分方程组；第二步进行Hankel变换，求解Bessel方程；第三步进行Hankel逆变换；第四步进行Laplace逆变换。可求出s(r,t), s1(r,z,t), s2(r,z,t)。详细求解过程可参考：张蔚榛：地下水非稳定流计算和地下水资源评价，pp. 81-109；Hantush, M. S. , Modification of the theory of leaky aquifiers; J. Geophy Research, Vol. 65, No. 11地下水动力学讲稿_第五讲地下水向完整井的非稳定运动（续第四章）2、数学模型的解（1）三层结构，相邻含水层为定水头情况的解考虑弱含水层弹性释水时的简化近似解1）抽水初期 当 和三种情况有相同的近似解，主抽水层的降深公式为：式中考虑弱透水层弹性释水的越流井函数 可查表4-8，p. 119，曲线地下水动力学_第四章4-17.JPG见图-17。地下水动力学讲稿_第五讲地下水向完整井的非稳定运动（续第四章）余误差函数误差函数（概率积分） 计算系数 越流因素：地下水动力学讲稿_第五讲地下水向完整井的非稳定运动（续第四章）2）抽水时间较久时的近似解a）相邻含水层为定水头情况下 降深公式可近似表示为 式中 W(u1,α)为不计弱透水层弹性释水的越流系统井函数。当 和地下水动力学讲稿_第五讲地下水向完整井的非稳定运动（续第四章）b）两相邻层为隔水层 降深公式可近似表示为 式中：W(u2)为无越流含水层的井函数；当 和地下水动力学讲稿_第五讲地下水向完整井的非稳定运动（续第四章）c）上相邻含水层为定水头含水层，下相邻层为隔水层 住含水层的降深近似公式 式中： W(u3,r/B1)为不计弱含水层弹性释水越流系统井函数；当 和地下水动力学讲稿_第五讲地下水向完整井的非稳定运动（续第四章）3、公式讨论（1）当B1=B2→∞，或μ*1=μ*2=0时，有β=0，即： erfc（0）=1； 主含水层降深公式简化为Theis公式（2）对第一、三情况（定水头）下，当K2=0， μ*1=μ*2=0，此时即为不考虑弱透水层弹性释水的情况。（3）由地下水动力学_第四章4-17.JPG图4-17，H(u,β)~1/u曲线，反映了S与t（1/u）和β的关系，总的来看，s随β的增大而减小，当β=0时，曲线即为Theis曲线；可以得出： 1）当μ*1、μ*2增大（