地下水向完整井的非稳定运动.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 ≤0.1时（即），相对误差不超过5%。一般生产上允许相对误差在2%左右。因此，当u≤0.01或u ≤ 0.05时，井函数可用级数的前两项代替，即： 于是，Theis公式可以近似地表示为下列形式： (4-13)式称为Jacob公式(1946)。 流量阶梯状变化时，当ui≤0.01时，即 (4-12)式可近似地表示为： 4 .对Theis公式和与之有关的几个问题的讨论1) Theis公式反映的降深变化规律将（4-11）式改写成无量纲降深形式，即，并给出曲线 [图4-3(a)]。曲线表明，同一时刻随径向距离r增大，降深s变小，当r→∞时，s→0，这一点符合假设条件。同一断面(即r固定)，s随t的增大而增大，当t=0时，s=0，符合实际情况。当t→∞时，实际上s不能趋向无穷大。因此，降落漏斗随时间的延长，逐渐向远处扩展。这种永不稳定的规律是符和实际的，恰好反映了抽水时在没有外界补给而完全消耗贮存量时的典型动态。图4-3反映了上述结论。 2) Theis公式反映的水头下降速度的变化规律将(4-9)式对t求导数，得：式(4-16)表明，抽水初期随着r的增大，值减小。因此，近处水头下降速度大，远处下降速度小。当r一定时，(4-16)式又表明，不同时刻的水头下降速度，由于和两个因素起着增、减两个方向相反的作用， 所以不是t 的单调函数； s-t曲线(图4-3b)不能沿着同一斜率变化，存在着拐点。可以利用，找出拐点的位置。为此有：所以 拐点出现的时间(此时u=1)为：图4-3的曲线也反映了上述结论，即每个断面