第6章地下水的非稳定渗流运动

第6章地下水的非稳定渗流运动随着工农业生产的不断发展，以及人口数量的不断增加，工业、农业及生活用水的需求量的不断增大，地下水作为重要的供水水源其开采量及开采规模迅速扩大，大多数地区普遍出现区域地下水的持续下降，而作为地下水运动要素均不随时间发生变化的稳定流理论及其裘布依（Dupuit）水量计算公式，无法解决和预测这一现象，以及未来地下水动态的变化趋势。本章主要讨论由于抽水而产生的非稳定渗流。

非稳定渗流理论所解决的主要问题1.评价地下水的开采量2.预报地下水位下降值3.确定含水层的水文地质参数泰斯以达西定律为基础，利用热传导理论提出了地下水非稳定井流的计算公式，称为泰斯公式。

泰斯非稳定流理论认为在抽水过程中地下水的运动状态是随时间而变化的，即动水位不断下降，降落漏斗不断扩大，直至含水层的边缘或补给水体，而且距抽水井越远，漏斗的曲率越小，扩展速度越来越缓慢。6.1非稳定渗流基本概念及其基本微分方程侧向边界离井很远，可不考虑其影响时，按无越流补给时处理，此时越流补给强度e=0。6.1.1轴对称二维不稳定潜水井流基本微分方程本节所要研究的问题是：在均质、各向同性、隔水底板水平的无限含水层中，单个完整井进行抽水的情况（考虑为二维流动）。渗流遵守达西线性定律，渗入强度为e。取一以井轴为中心的单元环柱体作为均衡地段，以dt为均衡时段。设断面r的流量为Q，断面r+dr的流量为Q+dQ，则均衡方程为：根据达西定律V=kJ可得上式的负号，是表示Q与h/r的方向相反，有：简化为

将上式代入式（6.1）:（6.2）

（6.1）

使式（6.2）线性化的方法，常用的有下列两种。第一种线性化的方法，是将式（6.2）左端部分中作为乘数的h用平均值hm代替，并视为常量，则式（6.2）可改写为当无渗入时（e=0），方程可写为对于水平二维无压流动，令则式（6.3）和（6.4）可写成为（6.3）

（6.4）

式中a为水位传导系数，m2/d；（6.5）

第二种线性化的方法，是在式（6.2）的两端均乘以h，并令势函数得：

再以平均值hm代替h，，并将式（6.5）代入上式，得当e=0时：令：T=kh——导水系数，表示含水层的导水性能；将T、a代入上式则得潜水完整井非稳定流的微分方程：或6.1.2不稳定承压井流基本概念及其基本微分方程1.承压含水层的弹性水量首先分析：在承压含水层中抽水（假定含水层的顶底板是不透水的，而且抽水时保持承压状态），抽出的水是哪里来的？从潜水含水层中抽水，它导致含水层的疏干，表现为地下水位——自由液面——的下降，抽出的水量正是含水层被疏干部分的水量（当e=0时）。但是，从承压含水层中抽水，周围形成的降落漏斗并不是对含水层的疏干，而只是构成水头（压力）的降低。

压力降低为什么能释放出水来？物体均具有可压缩性，只是程度不同而已。当作用在物体上的压力增大时，物体的体积缩小，密度增大；反之，当压力减小时，其体积增大，密度减小。对于承压含水层（取一处于平衡状态的地层柱体来研究，见图6.2），含水层上覆岩体外部荷载的重量和大气压力由两部分力与其平衡，一是含水层多孔介质对它的反力ps，另一是承压水作用在隔水顶板上的浮托力p（p=hpg，其中hp是承压含水层顶面的测压高度；g是水的重率）。这是抽水前的平衡状态。如果发生水头降低。也即含水层中每点地下水的压力p减小，它将引起下列作用：（1）由于水压的降低，地下水的体积发生膨胀，从而释放出部分地下水；（2）水的压力p的降低，即地下水对上覆岩体的浮托力降低，为了维持平衡，这部分力将转嫁到含水层多孔介质上，从而压缩含水层，其结果使含水层的空隙率n变小和含水层厚度变薄，这两个因素均使得从含水层中释放出部分地下水；（3）由于压力的降低，组成含水层骨架的固体部分将会膨胀，而这又引起含水层厚度和空隙率的变化，其关系比较复杂。考虑到含水层固体部分的压缩性一般比水和含水层要小得多，因此，建立微分方程时可以忽略固体部分的压缩性，将它视为刚体。如果承压含水层测压水头上升，则发生相反的过程。上述分析说明：假如水头降低，承压含水层会释放出部分地下水；如果水头升高，承压含水层也会储存部分地下水，这就是通常所说的“弹性储量”。弹性储量提供承压抽水井水量的概念，是与齐姆稳定井流的设想不相同的。后者假设，承压抽水井全靠“水平补给”。可以想像，如果没有弹性储量，依据水流连续性原理，则在抽水开始的一刹那，各断面(包括r→∞)的流量均等于抽水井的流量Q。或者为了把矛盾暴露得更突出些，考虑承压含水层中沟流的情况，则在刚抽水的一瞬间，各断面(包括r→∞)的流速均相等。这显然不符合实际情况。因此，弹性储量必须加以考虑。承压含水层由于水的来源是含水层的弹性压缩与水的弹性膨胀，因此其基本微分方程的建立除根据水均衡原理和渗流基本定律外，还应与水及含水层的状态方程（体积与压力间的关系）有关。2.水的状态方程假定水近似地符合弹性变形，依虎克定律，有因为V随p的增大而减小，即dV/dp<0，而bw规定为正值，所以上两式右侧有一负号。bw的物理意义是：当压力改变一个单位时，单位体积水的增量。bw的单位：压力的单位通常采用kg/cm2（大气压），故bw的单位采用cm2/kg。对方程（6.12）进行积分依照马克劳林级数:在压力变化不大时上式可近似取头两项,

（6.13）

式（6.13）可写成（6.14）

（6.12）

压力p的变化引起水的体积V的变化，但是水的质量m和重量G是不变的。由Vr=m和Vg=G的关系可知：若体积V增大，则密度r和重度g相应减小，有∴由式（6.12）得3.岩层（多孔介质）的状态方程4.轴对称二维不稳定承压井流基本微分方程这里讨论均质、各向同性、等厚的承压含水层中完整井的抽水情况。考虑含水层底板（或顶板）为弱透水层，抽水时通过它有越层渗透，其越流强度为e。设断面r的重量流量为gQr，断面r+dr的重量流量为gQr+d（gQr），单元环柱体中水的重量为G，则其均衡方程为按达西定律则

均衡段内地下水的重量为式中n——空隙率，

可推得：此式与潜水井微分方程式（6.3）对比，从形式上看，承压含水层中的Mg

(nbw+bs)起着潜水含水层的给水度m的作用。自己看P84~P86，可推出：

（6.32）当无越流时e

=0时

（6.33）上两式就是轴对称二维非稳定承压井流的基本微分方程。对比式（6.33）和（6.10），如果令j=M(H－h)=Ms（6.34）则式（6.33）可写成（6.35）

式（6.35）和（6.10）的形式完全相同，只是其中的势函数j不同。对非完整井，可推导出均质各向异性介质中地下水三维流动的微分方程为：或

6.2无越流含水层中的单个定流量完整井流因为无越流故无垂向渗流所以：e

=0。6.2.1无限承压含水层中单个定流量完整井流含水层均是有限的。如果含水层是如此之大，以致边界对于含水层研究区段的水头分布没有明显的影响，则可称它为无限含水层。对于压力传导系数小的含水层进行短时间抽水的情况，可视为无限的含水层。1.承压含水层定流量抽水时的Theis（泰斯）公式。承压含水层中单个井定流量抽水的数学模型是在下列假设条件下建立的：（1）含水层为均质各向同性、等厚、侧向无限延伸，产状水平、导水系数T为常数；（2）抽水前地下水的水头面水平；（3）定流量抽水井，井径无限小；（4）含水层中水流服从达西定律；（5）抽水后，水头下降引起地下水从贮存量中的释放是瞬时完成的，储水系数m*为常数储水系数是常数，故承压含水层的压力传导系数a=T/m*是常数，其中m*是释水系数。所以：导水系数T=kM，式中M为含水层的厚度。（6）承压井流要保持承压状态，则水位降深s不得大于（H-M）。（7）承压井按泰斯理论，抽出的水来自含水层储存量的弹性释放，并且是瞬时完成，所以T=kM为常数。抽水井抽水时，抽水量完全是来自含水层中的储存量，随着抽水时间的延续，以井轴为对称的下降漏斗不断的扩展，水流始终是非稳定渗流状态。在上述假设条件下，将坐标原点放在含水层底板抽水井的井轴处，井轴为z轴，如右图所示。根据上述假设条件(1)和(5)可以应用不稳定承压井流基本微分方程（6.6）,(2)是初始条件,(3)和含水层侧向无限延伸是内外边界条件。令j=kMH+C=M1(H－h)=M1s因此，该定解问题可写成为：（e=0时）

（常数）

该定解问题可用积分变换法，分离变量法或博尔兹门变换法求解，它的解是其中令∵

j

=M(H－h)=Ms式（6.45）、式（6.46）和式（6.47）称为泰斯公式。W（u）称为承压水定流量的井函数。右端为指数积分，将其展成幂级数并逐项积分。W（u）可表示为下列无穷级数的形式因为0.577216=ln1.78107，当抽水时间t较长、抽水井u≤0.01或观测井u≤0.05时，式（6.48）第二项以后（不包括第二项）可忽略不计，则式（6.45）、式（6.46）和式（6.47）可简化为见P89表6.1、P91[例6.1]。2.对泰斯公式的简要分析见P91,属于了解的范畴。6.2.2无限潜水含水层中单个定流量完整井流在潜水含水层中抽水时，潜水面是随时间不断变化的上界面。1.地下水向潜水井的非稳定渗流的主要表现以及与承压井的非稳定渗流的比较：（1）潜水井的导水系数T=kh是随时间和距离而变化的，而承压井的T=kM为常数；（2）潜水井降深较大时，垂直分速度不可忽略，在井附近为三维流。水平含水层中的承压井可作为二维流处理；（3）潜水井中抽出的水量主要来自含水层的重力疏干。考虑疏干的滞后性时，虽然潜水面下降了，而潜水面以上新形成的饱和带中仍有水继续向下排水，补给潜水层。这时，潜水层的给水度是变化的，故潜水含水层的水位传导系数a=T/m是变化的。给水度m是随着抽水时间的增加而逐渐趋向于稳定的最大值，一般提供的给水度值，便是这个最终值。导水系数，式中：。从实测的水头降深s和t时间的关系曲线分析，潜水含水层存在滞后疏干的现象，s－t曲线反映出抽水过程三阶段：抽水初期s－t曲线与承压井的泰斯曲线一致。这时主要为弹性释放，可以储水系数表示。潜水位下降了，重力疏干因滞后反应所引起的作用还很小。所以，含水层的作用相当于贮水系数小的承压含水层。这阶段的时间很短，也许只有几分钟。抽水中期s－t曲线的变化很象有越流补给时半承压含水层的情况，明显偏离泰斯曲线，曲线斜率变小，甚至出现短时间的假稳定。此时，重力疏干的作用逐渐明显，贮存的重力水逐步释放出来，起到连续再补给的作用。弹性释放的作用依然存在，但所占比例已逐渐减弱。此时给水度的数值逐级增大而趋向潜水含水层的最大值。这个阶段，由于潜水层性质的不同，可能是几分钟，也可能是几天。抽水后期s－t曲线又与泰斯曲线相一致，水头下降速度增大，降落漏斗扩展，这时主要为重力疏干作用，由于抽水时间的增长，重力排水已跟得上水位的下降，滞后作用可忽略不计，此时的给水度达到最大值。（4）当降深不大时，含水层为无限含水层时潜水完整井单井抽水非稳定流运算模型参照承压水完整井的方式进行一系列代换导出时，将式（6.5）代入式（6.43），则得计算潜水井流的基本方程：上三式也称为泰斯公式。

所以泰斯公式有6个式子，即：公式（6.45）、式（6.46）、式（6.47）、式（6.58）、式（6.59）和式（5.60）。式中W（u）的求解与承压水完整井非稳定流时相同。当抽水时间t较长、抽水井u≤0.01或观测井u≤0.05时，式（6.48）第二项（不包括第二项）以后可忽略不计，则式（6.58）、式（6.59）和式（6.60）可简化为2.博尔顿法数学模型及其解*井函数的数值解在双对数坐标纸上绘得曲线簇,曲线，适用于抽水初期；

它包括两组曲线，

左边为A组

右边为B组曲线，适用于抽水后期。因式（6.63）的条件是h

→∞，两组曲线的中间部分为一水平线，可用式（6.65）表示。当h

≥100时，曲线中间部分仍趋近于一水平线，如h<100时，曲线中间部分就不是水平线了，而是一条比抽水初期和后期的斜率小得多的曲线，可用两组曲线间它们的共同切线连接。曲线簇说明：抽水初期以弹性释放为主，水头降深与左边的泰斯曲线吻合；中期滞后重力排水的影响曲线簇偏离泰斯曲线，水头下降速度变小，并随不同的r/D，以不同方式向水平线趋近；后期滞后重力排水影响逐渐减弱，水头下降速度又由小变大。滞后重力排水影响基本结束时，曲线簇与右边的泰斯曲线趋近。；6.3有越流补给时承压含水层中的单个定流量完整井流如果抽水层的顶板或底板不是隔水层，而是弱透水层或弱含水层，那么当抽水层抽水时，水位下降，抽出的水除了抽水层本身的弹性释放之外，还得到弱透水层的弹性释放补给和相邻含水层通过弱透水层的补给。这种含水层系统称为越流系统。它包括抽水含水层和相邻含水层。见下图。越流系统可分为三种类型：（1）弱透水层的弹性储量很小，可忽略不计，而且在抽水含水层抽水期间，相邻含水层的水位不变。（2）弱透水层中释放出来的水量相当大，甚至是越流补给的主要来源，相邻含水层的水位保持不变；所以应考虑弱透水层的弹性储量。（3）补给层的水头随主含水层的抽水情况而变化，这种类型的计算很复杂，本书不讨论此类问题。本书只介绍第一类及第二类特定条件下的计算公式。越流系统中常用的三个系数（1）弱透水层越流系数b：其含义为，水头降低一个单位时，单位时间内，单位水平面积上，相邻含水层通过弱透水层补给抽水含水层的水量。（2）阻越流系数B阻越流系数B表示越流补给量大小，B愈大则渗透系数就愈小，故垂直补给量愈小。（3）越流补给强度e越流补给强度是单位时间通过单位水平面积补给抽水含水层的水量，（m/d）

（6.68）式中s——抽水含水层水头降深。一、第一类越流系统地下水流向承压完整井的非稳定运动1.汉土斯假设：（1）越流系统中每层都均质各向同性、产状水平、等厚、侧向无限延伸；（2）抽水含水层和相邻含水层初始水位水平且相等，抽水后，抽水层中水流为平面径向流；（3）抽水后相邻含水层越流补给抽水层，但其中水位保持不变；（4）抽水层中水流服从达西定律；（5）水和含水层均为弹性体，储水量的释放是瞬时完成的；（6）抽水井以定流量抽水为井径无限小的完整井；（7）弱透水层的弹性释放水量可忽略不计，通过其中的水流为垂向一维流。在上述水文地质条件下抽水时，井的抽水量Q由二部分组成：一是由于水位降低，抽水含水层本身的弹性释放量；二是在抽水含水层与相邻含水层之间的水位差作用下，来自相邻含水层的越流补给量。随着抽水时间的延续，水位差增大，降落漏斗扩展，使越流补给量在井的抽水量中所占比重也逐渐增大。当井的抽水量与越流补给量相等时，抽水含水层的降落漏斗达到稳定，抽水含水层不再释放储存量。2.汉土斯数学模型及其解在上述假设条件下，可应用承压二维渗流的微分方程，结合相应的初始条件和边界条件，构成一个理想的数学模型：s（r，0）=0=0这一数学模型有非稳定流和稳定流的解。（1）有越流补给时非稳定流的解将上式积分变换，得：

解得：为不考虑弱透水层弹性释水时越流系统的井函数，其值见P101表6.2。（2）有越流补给时非稳定流的解的分析①有越流时的降深比无越流时小将有越流时的解式（6.71）和泰斯的解式（6.45）对比可看出，当u值相同时，因为恒为正值，故积分较为小。因此，该含水层的降深s比无越流的承压含水层的降深值要小。这是因为在越流时，水井中抽出的水，一部分来自越流补给，抽水含水层可以少释放一些弹性储量造成的。当k/=0时，弱透水层变为不透水层，越流因素B→∞，没有越流补给。式中的→0，则变成泰斯公式。时，就是泰斯曲线。②s—t曲线的形状按表6.2有越流井函数绘制—标准曲线，见下图。此曲线相当于s—t曲线。抽水初期，降深较小，越流尚未进入抽水层，井中抽水的水量几乎全部是消耗抽水层的贮存量，可看出这时的降深曲线与泰斯曲线一致。在理论上即越流量等于0，或B→∞，→W(u)，降深曲线与泰斯曲线一致。值时的曲线形状。在其他条件和越小，与泰斯曲线一致的过程越长，相邻含水层的分析，小，即B大，则透弱水层的渗透性小，厚度大，标准曲线中表示了不同r相同时,越流量进入抽水含水层的时间越迟。从所以越流发生得迟。抽水中期，降深曲线变缓，偏离泰斯曲线。这说明越流开始进入抽水层，这时抽水量由两部分组成：一是来自抽水层的弹性释放，二是来自相邻含水层的越流补给。其他条件和r相同时,大，越流补给量大，消耗抽水层的储存量小，所以抽水层的降深也小，偏离泰斯曲线早。抽水后期，降深曲线趋向水平直线，也就是说水头不再下降，抽水量等于越流补给量，流水由非稳定流变为稳定流。③水头下降速度当存在越流时，含水层中水头下降速度为

与泰斯公式的水头下降速度对比看出，恒为正值，故有越流时水头下降速度比无越流的承压含水层为小。t足够大时，井周围的降落漏斗等速下降。（3）有越流补给时的稳定流解当t足够大时，u便很小。实用上只要u≤导，井函数部分为：

，经数学推一般越流含水层的阻越流系数B都相当大，故《1，上式要求[注：见表5.5，大，则小]

。稳定时的降深即最大降深，可表示为：式中

为零阶第二类虚宗量塞尔函数，见P103表6.3。

能满足当≤0.05时，在抽水井附近，则按贝塞尔函数性质：则稳定流中求得的裘布依公式为

对比上两式得:R=1.12BR就是裘布依公式中的影响半径，在一定水文地质条件下B是常数，因而R也是定值。此时，R的意义是在抽水效果不变的条件下，把面积形式补给的越流量转化为具有圆柱形侧向补给时的圆周半径。见P104[例6.2]

第二类越流系统水流向承压完整井的非稳定运动的计算公式当弱透水层上下为两个相邻含水层，抽水时间又很短的第二类越流系统（弱透水层中释放出来的水量相当大，不能忽略时）可用下列公式计算：

erfc(x)——误差函数的补函数。

——井函数自变量；——第二类越流系统井函数；6.4集中开采区定流量完整干扰井群在集中开采区，往往在较小范围内布置数量很多的单井，形成相当密切的干扰井群。由于单井的数量很多，用前面讲的叠加方法进行计算是相当繁琐的。为了简化计算，可以把集中开采区的干扰井群近似概化为具有均匀开采（或补给）强度的开采地区，变成假想的连续型干扰井群。若井群的总开采量为Q，开采区面积为A，则开采强度e=Q/A。实际的干扰井群无论怎样密集总是离散型的，概化为连续型只是一种近似，把井群作为一个整体看待。在平面上按井的分布应用较多的连续型开采地区有圆形和矩型。在含水层稳定展布的平原区，各种井群往往不均匀地分布着，难以概化为单一的具有均匀开采强度的圆形或矩型开采地区。因此，可以按井孔的分布和流量的情况、概化成若干个不同开采强度、不同几何形状的开采地区的组合。还可加上若干个难以包括的单井,见上图。然后将若干个概化地区和单井叠加起来，就可求解。6.4.1圆形开采区右图表示一平面无限延伸的承压含水层，厚度为M，压力传导系数a=T/m*，圆形开采地区的半径为R，均匀开采强度为e，总开采流量为Q，开采后区域降落漏斗、概化如图中所示。水流为轴对称的，可列出以降深表示的数学模型：可推得开采中心区最大降深

当u0≤0.05，即抽水时间较长，时

时，

时，

时，

此式表示在距离比较远处，圆形开采区的作用与单井的作用相同，可以用单井公式计算。见P106[例6.3]。6.4.2矩形开采区图6.12表示一平面无限延伸的承压含水层，厚度为M，压力传导系数a=T/m*，矩形开采地区的边长为2L和2b，总开采流量为Q，开采后区域降落漏斗、概化如图中所示，均匀开采强度为e

=Q/4Lb，

在开采中心x=0，y=0处降深最大，上式简化为：见P106[例6.4]。6.4.3潜水含水层在降深不大时可近似应用承压井的公式（6.89）见P107[例6.5]。6.5无越流含水层中水流向完整干扰井群的非稳定渗流6.5.1.1承压完整干扰井群根据叠加原理，若n个井同时抽水，t时刻在A点的水位降深为：当时，其近似公式为：6.5.1.2潜水完整干扰井群对潜水含水层在降深不大时可近似用下式：当时，其近似公式为：6.5.2边界附近地下水向井的非稳定渗流6.5.2.1直线补给边界承压含水层一侧为定水位的河流，另一侧无限延伸。有一抽水井在工作。设以河流边界为镜面，在另一侧对称位置上有一注水井，以Q注=Q抽工作。根据叠加原理可写出直线补给边界附近的单井非稳定抽水时，t时刻任意一点A的降深计算公式：当u1≤0.01，u2≤0.01时，可用下列近似公式：

上式的右端没有时间t的变量，说明sA不随时间而变化，和稳定流的方程式一样。所以有补给边界存在的条件下，抽水井工作一段时间后，便可达到稳定状态。见P109[例6.6]。6.5.2.2直线隔水边界如工作井是在承压含水层中的抽水井，虚构井可看成是抽水井，流量同真实的抽水井一样。根据叠加原理可写出直线隔水边界附近的单井非稳定抽水时，t时刻任意一点A的降深计算公式：当u1≤0.01，u2≤0.01时，可用下列近似公式：6.6地下水向非完整井的非稳定流渗流运动基本方程式对于承压水非完整井：对于潜水非完整井：6.7水文地质参数的确定利用稳定流抽水试验计算水文地质参数

一、单井稳定抽水试验计算渗透系数k

利用裘布依型稳定流公式进行渗透系数计算时，若没有观测孔而只能根据抽水井的出水量、水位下降等数据，则应消除抽水井附近产生的三维流、紊流的影响；特别是在抽水井水位下降值较大的情况下，最好采用下列消除渗透阻力的方法：首先根据单井内水位下降值s与相应的出水量Q绘制Q—s关系曲线，如右图所示，再按所得曲线类型选择适当的计算公式。1.当承压井Q—s(或潜水井Q—)呈直线关系时（见右图中的a型曲线），地下水运动为平面流，可直接采用第五章的公式计算渗透系数k。(注：潜水含水层在自然情况下：

)。见P112[例6.7]。2.当Q～s(或)关系呈曲线关系时（见图6.14中的b型曲线），抽水井壁及其附近含水层中，已产生三维紊流，不符合裘布依的基本假定条件，因此不能直接用稳定流公式进行计算。为了消除三维流、紊流的影响，在计算时应采用消除阻力法。方法如下：首先绘制或关系曲线，若根据三次水位下降的Q、s

P113[例6.8]值所做的承压水的或潜水的则可将直线在纵轴上的截距a值代入第五章的公式计算渗透系数，这种方法称为截距法。即用a代替承压井公式中的；用a代替潜水井公式中的，然后再求渗透系数k。

呈直线时，关系曲线式中

s0=H-h0

；h0——井中水深。二、带观测孔的单井稳定抽水试验计算渗透系数k为了避免抽水井附近的三维紊流影响，所选的最近观测孔距主井的距离一般为含水层厚度的一倍，而所选的最远观测孔距第一个观测孔的距离也不宜太远，以保证各观测孔内有一定的水位下降值，并使各观测孔的水位下降值在的直线段内。1.有一个观测孔

—lgr曲线承压水完整井潜水完整井当观测孔距抽水井较近时，易受三维紊流影响，如采用6.7.1.1中的方法，则渗透系数偏小；当远离抽水井时，如采用6.7.1.1中的方法，则渗透系数偏大。2.有两个观测孔时承压水完整井：潜水完整井：

承压水非完整井（过滤器紧接含水层顶板，L<0.3M，r1=0.3r2，观测井和抽水井深度相等）：潜水非完整井（抽水井过滤器淹没）：见P115[例6.9]

三、影响半径R的确定裘布依公式是在抽水井位于岛状含水层中心位置的条件下建立的，因此公式的特定含意为：在距抽水井R远处有一个实际的定水头补给；R表示实际的影响范围；R处的水位下降值为零；R值是一个不受抽水下降值s及抽水量Q影响的常数值。生产实践表明，抽水的影响范围是随抽水时间t的延长、流量Q的增加而扩大的，不可能把抽水的影响范围限定在一个“半径”内。同时，在天然条件下，降落漏斗多不对称，边界也不明显，单井抽水影响范围实际上不是一个圆。不能简单地用一个“半径”来确定。目前确定影响半径的方法主要有如下几种。1.无观测孔（1）不考虑地下水流向潜水完整井承压水完整井

（2）考虑地下水流向①承压含水层地下水上游方向地下水下游方向

地下水流向垂直方向

抽水时井的降落漏斗为椭圆形，三者之间的关系为R1=7.3R2=2.7R3②潜水含水层式中③经验公式承压含水层（集哈尔特公式）潜水含水层（库萨金公式）2.有一个观测孔（1）承压含水层：r0≤r1≤0.178R时

（2）潜水含水层

式中r1——抽水井至观测孔之间的距离（m）；

s1——观测孔内水位下降值（m）。其它符号同上。3.有两个或两个以上观测孔（1）承压含水层（2）潜水含水层

式中s1、s2——在s－lgr曲线上任意两点的水位下降值；、——在－lgr曲线上任意两点的纵坐标值；r1、r2——s()－lgr曲线上，纵坐标为s1、s2或、的两点至抽水井的距离。4.图解法确定影响半径（1）浸润曲线外推法在制图纸上，将主井抽水时，最大一次水位下降稳定后，测得的主井及各观测井的稳定水位绘在图上，以光滑的曲线按自然趋势连接，并外推使其与自然水位线相交于一点，该点至抽水井的距离，即为影响半径。如右图所示。（2）s—lgr直线交绘法此方法需观测孔不少于三个。在半对数坐标纸上，按抽水井至观测孔距离标出抽水井、观测孔位置，然后根据抽水同一时刻测得的抽水井和观测孔的稳定水位，点在图上。其观测孔各点连线，应为直线，并将各条直线外推，与lgr轴交于一点，此点的值即为影响半径。如右图所示。无越流含水层中利用非稳定流抽水试验计算水文地质参数在无越流含水层中进行非稳定流抽水试验，主要是为了确定含水层的导水系数T、储水系数u*或压力传导系数a。计算这些参数常用的方法有：试算法、配线法（也叫标准曲线对比法或量板法）、直线图解法、恢复水位法等等。一、试算法此法的优点为计算简单易懂，缺点为计算结果精度不高。1.当有一个抽水井和一个观测井时使抽水过程中流量Q保持不变，若取两个时段t1、t2及相应的观测井中的水位下降值s1、s2，则有：；

上式中只有压力传导系数a是未知数，但由于它居于井函数之中，一般情况下提不出来，不能直接求解。通常采用试算的方法，根据经验给定一个a值，代入式上看能否满足，如果不满足就再另选一个a值直到满足为止。当给定一个a值之后，则为已知数，即可在表6.1中查得W（u）值，如果之比值等于时，则为满足。上述试算的过程还可用作图的方法来代替，如下图所示。设以a作为横坐标，给以任意的a1、a2、a3……an值代入式（6.125）可得到相应的B1、B2、B3

……Bn值，把这些数值点在右图上，即可得曲线1，然后再根据实测的值，在曲线1上确定A点，A点在横坐标上的投影点就是所求的a值。a值确定以后，可按右式计算导水系数T：，以B作为纵坐标，如果没有观测井，亦可粗略地把抽水井本身当作观测井。2.当有两个观测井时可取任意同一时间t的两个观测井的s1、s2和r1、r2进行计算，这样得：

二、配线法这是通过实测（试验）曲线与理论曲线（标准曲线）对比确定含水层水文地质参数的方法。此种方法可分为降深—时间距离、降深—时间及降深—距离共三种配线方法。当只有一个观测井资料时应采用前两种配线法，若有两个以上观测井资料时可采用降深—距离配线法。1.配线法的优缺点(1)优点：因能充分利用抽水试验的全部资料，故能避免个别资料的偶然误差，所以此法的最大优点是计算精度高。(2)缺点：有以下两个缺点①抽水初期实测曲线常与理论曲线不符。其原因主要是推导公式时应用了一些假设，这些假设与实际不符，例如假设贮存量瞬时释放，而实际上总有一个过程。此外抽水初期涌水量不易稳定与理论要求不符合，故非稳定抽水试验时间不宜过短。②抽水后期曲线比较平缓时，同理论曲线不容易拟合准确，常因个人判断不同引起误差。故在确定抽水延续时间和观测精度时，应考虑所得资料能画出s—t或s—曲线的弯曲部分，便于拟合。如果后期实测数据偏离标准曲线，则可能是含水层外围边界的影响或含水层岩性发生了变化等。这就需要把试验数据和具体水文地质条件结合起来分析。2.降深~时间距离配线法（1）原理承压完整井的非稳定流公式可表示为对上两式变换形式并取对数得上两式右端的第二项在同一次抽水试验中都是常数，所以在双对数坐标系内，对于定流量抽水s—曲线和标准曲线

标准曲线在形状上是相同的，只是纵横坐标平移了（2）具体步骤①从井函数表6.1中选5—6对以上的u~W(u)的对应值，在双对数坐标纸上绘制（理论曲线）；和记下对应的坐标值，代入承压水完整井的非稳定流公式即可确定有关参数。的距离。只要将两条曲线重合，任选一匹配点，

②在与理论曲线张双对数坐标纸上选择适当的坐标系，根据实测点据绘制s—曲线，使理论曲线与s—曲线重合理论曲线基本上在实测点据的点群中央]。

[注：③任选一匹配点（在曲线上或外均可），记录匹配点的对应坐标值W（u）、、s和T和储水系数。

代入泰斯公式分别计算为计算方便，尽量取数简易的点作为匹配点；如匹配点选在W（u）=1，上。

上。同一导水系数3.降深~时间配线法的具体步骤（此时r为定值）（1）由实测点据在双对数坐标纸上绘制s—t曲线；在与s—t曲线同一张双对数坐标纸上选择适当的坐标系绘理论曲线，使s—t曲线和

曲线基本重合[注：上在实测点据的点群中央]，曲线基本只是纵坐标和横坐标平移了、、的距离。

对承压完整井的非稳定流公式变换形式并取对数得：（2）任选一匹配点，由匹配点读出W（u）、、s和t；（3）计算出T和u*。见P122[例6.10]。4.降深—距离配线法的思路在实测值中任取一个t为定值，利用所有观测孔的降深值，在双对数坐标纸上将s—r实际资料曲线与标准曲线拟合，最后求出水文地质参数。直线图解法

本章只介绍利用一个观测井的资料来确定含水层参数的方法。将式（5.37）换算成常用对数

在某一抽水过程中，上式中T、a、r、Q均为常数，故s—lgt在半对数坐标中呈直线关系，所以上式可变换为如下形式：

故

将直线延长并与横轴s=0处交于t0点，1.计算原理:（1）将坐标（t0，0）代入（1）式则可得：；

；

则：故：

2.计算步骤（1）求s0、t0和直线的斜率m以lgt为横坐标、s为纵坐标，在观测井资料的点群中央作一条直线交纵坐标得s0，交横坐标得t0。并求出直线的斜率m，见右图。（2）用式（6.135）计算T，用式（6.136）计算a值。见P124[例6.11]

恢复水位法1.计算原理假如某井以定流量Q进行抽水，持续进行了tP时间之后停止抽水测定恢复水位，则时间tP之后的剩余水位下降值s，可以考虑为该井仍以流量Q继续抽水，并从停止抽水的时间起有一个流量Q的注水井开始工作，这样正负流量相抵消，得到了停止抽水的效果，如右图所示。根据势的叠加原理，停止抽水后的剩余水位下降值（6.45）得到:——恢复水位持续时间。tP——抽水持续时间；可按公式当时:

从上式可以看出：以实测值s为纵坐标，点群中央绘制一条s—为横坐标，可在直线；则该直线的斜率m为：只适用于抽水的非稳定流过程2.计算步骤见P126。[例6.12]越流系统中利用非稳定流抽水试验计算水文地质参数表示越流补给的承压含水层的水文地质参数有：导水系数T、储水系数用的是标准曲线对比法、拐点法等。下面重点介绍第一越流系统有关参数的几种求法。、越流系列数B等。确定参数的方法很多，最常一、标准（理论）曲线对比法1.计算原理标准曲线对比法是从式（6.71）出发，即：由上两式可以看出：与（或u）的关系和s与t（或r2）

同的，据此则可将井函数的关系是一致的，说明两者在双对数坐标系中的曲线形状是相及自变量

数表，并绘制成标准双对数曲线。

各值列成函u及值。

然后用实际观测的资料制成模数相同的双对数曲线，再用重叠法使两曲线达到最大的吻合，最后从相应的理论曲线上读出所需要的计算数值——

、2.计算步骤（1）在表6.3即