4地下水向完整井的非稳定运动解决方案.ppt

在上述假设条件下，抽水后将形成以井轴为对称轴的下降漏斗，将坐标原点放在含水层底板抽水井的井轴处，井轴为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）式，有： 两边积分得： 令 ，则有： 故： 利用r=0时的F（r）值，由（4-6）可以确定C值： 但由(4-7)式，有： 把上式代入(4-5)式，有： 为计算方便，对（4-8）式进行变量代换，令： 同时更换积分上下限，当τ=0 时， 当τ=t时，y = 于是， 其中， 在地下水动力学中，采用井函数W(u)代替(4-9)式中的指数积分式： 则(4-9)式可改写成： 式中，s——抽水影响范围内，任一点任一时刻的水位降深； Q——抽水井的流量；T——导水系数； t——自抽水开始到计算时刻的时间； r——计算点到抽水井的距离；μ*——含水层的贮水系数。 (4-9)式为无补给的承压水完整井定流量非稳定流计算公式，也就是著名的Theis公式。 为了计算方便，通常将W(u)展开成级数形式： 并制成数值表（表4-1），只要求出u值，从表4-1中就可查出相应的W（u）值；反之亦然。 2. 流量变化时的计算公式 Theis公式是在假定流量固定不变的情况下导出的。这种情况通常只有在抽水试验时才能做到。实际上，很多生产井的流量是季节性变化的。如农用井在灌溉季节抽水量大，非灌溉季节抽水量小。工业用水也有类似情况，常随需水量而变化。在这种情况下，怎样应用Theis公式? 首先需要绘出生产井的Q=f（t）关系曲线，即流量过程线。 然后将流量过程线概化，用阶梯形折线代替原曲线，坐标选择如图4-2所示。概化原则是矩形面积等于曲线于横坐标所围成的面积。其中，每一个阶梯都可视为定流量，应用Theis公式。把各阶梯流量产生的降深，按叠加原理叠加起来，即得流量变化时水位降深的计算公式。 当0tt1时，水位降深为： t时刻经历若干个阶梯流量后所产生的总水位降深为： 式中，设t0=0，相应的Q0=0。 (4-12)式为流量变化时，经概化呈阶梯状变化后的计算公式。 3 . Theis公式的近似表达式 如前述，Theis公式中的井函数，可以展开成无穷级数形式，即： 前三项之后的级数是一个交错级数。根据交错级数的性质可知，这个级数之和不超过u。也就是说，当u很小，井函数W(u)用级数前两项(-0.577216-lnu)代替时，其舍掉部分不超过2u。因此， 当u ≤0.01(即 )，井函数用级数前两项代替时，其相对误差不超过0.25%； 当u≤0.05时(即 )，相对误差不超过2%； 当u