1、1地下水运动中的专门问题6.l 非饱和带的地下水运动在地下水面以上的非饱和带（即包气带）也有水的运动。在许多情况下，研究非饱和带的地下水运动具有很大的意义。例如，在地下水资源评价中，必须研究“三水”(即大气水、地表水相地下水)的相互转化，而非饱和带的地下水运动是其转化的重要环节。入渗的水必须经过非饱和带才能到达潜水面，故研究水在非饱和带的运动，对于入渗的计算很重要。其次，各种施加在地表的污染物将随入渗的水一起运动，经过非饱和带进入地下水中。因此研究地下水污染时，也必须研究非饱和带中水的运动。6.1.1 非饱和带水分的基本知识1. 含水率、饱和度和田间持水量在非饱和带中，空隙空间的一部分充填了水
2、，其余部分充填了空气。水分和空气的相对份量是变化的。可以用二个变量来表示水分含量的多少。为含水率 ，表示单位体积中所占的体积：(6-1)式中， 为含水率，无量纲；(V w)0 为典型单元体中水的体积； V0 为典型单元体的体积；另个为饱和度 Sw，表示岩石的空隙空间中水所占据部分所占的比例：(6-2)式中，S w 为饱和度，无量纲；(V 0) 0 为典型单元体中的空隙体积。显然，含水率 不能大于空隙度 n。而饱和度 Sw 不能大于 1。两者之间有下列关系：nS w (6-3)因为利用了典型单元体的概念上述定义对于任一点都是适用的。在长时间重力排水后仍然保留在土中的水量称为田间持水量。此时，水以
3、簿膜水的形式和在颗粒接触点附近以孤立的悬挂环形式存在。从图 6-1 可以看出，空隙度减去田间持水量，相当于排水空隙度，即排水时的有效空隙度。2毛管压力当多孔介质空孔隙中有两种不相混溶的流体(如水和空气) 接触时，这两种液体之间的压力存在着不连续性。此压力差的大小取决于该点界面的曲率(它又取决于饱和度) ，这个压力差 pc 称为毛管压强： wa(6-4)式中， 空气的压强， p水的乐压强。如假设孔隙中的空气是在 101325Pa(一个大气压)下、并取大气压强作为测量流休压强的基准，则 ap0，于是：wcp(6-5)故非和带孔隙中的水处于小于大气压强的情况下。正如在毛细管现象中见到的一样，在周围水
4、面以上的毛管内的压强是负的。2和饱和带的情况一样，可以定义非饱和带水流中任何点的水头(毛管水头) ：（6-6）式中，为水的容重； wcph(6-7)称为毛管压力水头。某些作者用符号 表示压力水头的负值，即：0(6-8)zH(6-9)对于饱和非饱和流动，可以写出统一的水头表达式； p(6-10)式中，压强少可正可负。在饱和带中，p 为水的压强，取正值；在非饱和带，p 为毛管压强的负数，取负值。其余符号同前。图 6-1 非饱和带的含水量曲线图 6-2 土壤水分特征曲线（据 Richards 和 Weaner）3土壤水分特征曲线反映毛管压强 pc或毛管压力水头 hc和土壤含水率 或饱和度 Sw 关系
5、的曲线，称为水分特征曲线(图 6-2)。它表示非饱和带中水分的能量和数量之间的关系，反映了包气带中水的基本特征。从曲线上还可以看出，即使在相当高的压强下，土样中仍保持一定的水，含水率不再进一步减小。这个含水率记作 0，相应的饱和度为：nSw0不同土的水分特征曲线是不同的。在同样条件下，粘性土要比砂保持更多的水分，具有更高的含水率。土的颗粒级配，对持征曲线的形状也有影响，如图 6-2 的曲线 I 和 II。温度的变化对它也有影响。温度升高时，表面张力降低，在同样吸力下含水率要低一些。水分持征曲线斜率的负倒数称为容水度，记作 C：3cdhHC （6-11）容水度不是常数，它随含水率或毛管压强而变化
6、，记作 C()或 C(hc)。它表示毛管压力水头变化一个单位时从单位体积土中释放出的水体积，是计算非饱和带水运动的重要参数。图 6-3 吸湿和排水情况下的水分特征曲线 (据 J. Bear)实验表明：同一土样在同样的温度下，排水过程和吸湿过程的水分特征曲线是不同的(图 6-3)。在同一 pc 或 hc 下，排水时的含水率要大于吸湿时的含水率。这种现象称为滞后现象。土壤从饱和到干燥或从干燥到饱和的水分持征曲线称为主线。土样从部分湿润到开始排水或从半干燥状态重新润湿时，水分特征曲线是顺着一些中间曲线由一条主线移至另一条主线，这些中间曲线称为扫描曲线。因此，水分持征曲线随土壤的干、湿历史的不同而变化
7、。故容水度 C()不是含水率 的单值函数。4非饱和流动中的给水度概念已经介绍过给水度的概念。给水度是单位体积含水层中所排出的重力水的体积。但实际上，当潜水面下降时，其间的水并未全部排出，只是由饱和带的水变成非饱和带的水，水分分布曲线发生相应的改变。实际排出的水体积只相当于排水前后两条水分分布曲线间的那一部分面积。为此，需要这样来定义给水度：从地表一直延伸到含水层底板的一个单位水平面积垂直土柱，当潜水面降低一个单位时，由重力所排出的水的体积。由于重力排水的迟后，给水度 也是时间 t 的函数。只有当长时间排水后才趋近于某一常数值。6.1.2 非饱和带水运动的基本方程1. 运动方程1931 年，Ri
8、chards 提出，Darcy 定律可引伸应用于非饱和带水的运动。但此时的渗透率 k 和渗透系数 K 不再是常数，而与土壤的含水率有关。当含水率( 或饱和度)减小时，一部分空隙为空气充填，因而过水断面减小，渗流途径的弯曲程度增加，导致渗透率或渗透4系数减小。因此，该情况下 k 和 K 可记作含水率 或饱和度 Sw 的函数 k()，K( )或 k(Sw)，K(Sw)。这样，非饱和带中的 Darcy 定律表达式为：v=K()J (6-12)如用渗透率来表达时，则有： )()zpSkvw)()zpSkwr(6-13)式中：k 一饱和土的渗透率；k(Sw)非饱和土的渗透率，为饱和度 Sw 的函数；kr
9、(Sw)一 相对渗透率， kkwr)()(；水的动力粘滞系数。相对渗透率为非饱和土的渗透率和同一种土饱和时的渗透率的比值，为含水率 或饱和度 Sw的函数。非饱和砂的相对渗透率 和饱和度 Sw的关系表示在图 6-4 中。当饱和度 (含水率)减少时，大孔隙首先开始排水，渗透在较小的孔隙中进行，过水断面减小，渗流途径的弯曲度增加，相对渗透率急剧减小。到达 A 点，孔隙中的水变得不连续了，相对渗透率等于零。此时的饱和度为身 Sw0，相应的含水率为 0n Sw0。图 6-4 非饱和砂的相对渗透率与饱和度的关系(据 Wyckoff 和 Botset, 1936)2基本微分方程在第一章中，我们已经得到了渗流
10、的连续性方程(1-65)式。对于非饱和流动，把等式右端的空隙度 n 换成含水率 ，方程仍然是适用的。在非饱和带中，一般不考虑介质的变形，即单元体体积xy z 不随时间而变化。于是可以约去等式两端的 xyzt；同时，在非饱和流动中，水的密度变化很小，可当作常数。于是，相应的连续性方程为：(6-14)将运动方程(6-12)式代入(6-14)式中，得：(6-15)式(6-15) 即为非饱和流的基本微分方程，称为 Richards 方程。上述方程中，既含有含水率 ，又含有水头 H，为解决问题方便起见，可以把基本微分方程化成以下几种表达形式。(1) 以含水率 为因变量的表达式：前已述及，非饱和带的水头
11、Hz- hc，又由水分特征曲线、毛管压力水头 hc 和含水率 之间存在着函数关系，因此(6-15)式可改写为：5（7-16）由于水分持征曲线各处的斜率不同，C 不是常数而是随含水率 而变的变数，即CC( )。令：参数 D()是渗透系数和容水度的比值，称为扩散系数，量纲为L 2T-1。它是一个重要的参数。引入 D()以后，(6-16)式变为：（6-18）这是二阶的非线性偏微分方程。对于一维的垂直流动，可简化为：（6-19）z 轴向上取正值，z 轴向下取负值。(2) 以毛管压力水头为因变量的表达式：从水分特征曲线可知，毛管压力水率之间存在着函数关系。因此，非饱和土的渗透系数同样是毛管压力水头的函数
12、，即 K=K(hc)或K()，C C(h c)或 C()。于是，（6-20）（6-21）考虑到（6-22）(3) 饱和非饱和流的表达式：在饱和非饱和流动中，常以压强 p 或水头 H 为因变量，有 pzH。如果不忽略密度的变化，连续性方程(6-14)可写为：（6-23）再将 v 用运动方程(6-13)代入，容重 g，含水率 nS w，则得：(6-24)该方程中的某些参数的取值范围如下：6以上考虑的模型都是单相流模型，只研究水的运动，即凡是水流到的地方，空气自然被排走。实际上，岩石空隙是既存在空气也存在水的二相系统，也必然是更复杂的模型，这里就不介绍了，请读者参考有关的专著。思考题：1为什么对于饱
13、和流动，不透水边界的边界条件为0nH，而对于垂直入渗的非饱和流动，不透水边界的边界条件为 1zhc？2图 6-5 为一个垂直入渗情况下的饱和非饱和流动模型。地面入渗率为R(t)，潜水面的埋深 s(t)随时间而变化，初始埋深为 s(0) L，潜水含水层底板隔水。试写出该情况下的饱和非饱和流动的数学模型。图 6-5 垂直入渗情况下的饱和-非饱和流动模型他在建立微分方程时应用了下列假定：(1) 具有原生孔隙的岩层中广泛发育有随机分布的裂隙，二者都充满着整个研究区，形成两个重叠的连续系统。也就是说，孔隙和裂隙的分布彼此都是连续的。即所谓的“二重性”假定。双重介质的名称即由此而来。根据这一假定，在渗流区
14、中的每个点上都有两个水头，一个是孔隙水头 H，另一个是裂隙水头 Hf；(2) 孔隙以贮水为主，裂隙以导水为主，水自孔隙经裂隙流向别处，其总的渗透性决定于裂隙的渗透性；水在裂隙中的流动服从 Darcy 定律；(3) 孔隙和裂隙的初始水头相等，它们之间交换的水量与其水头差成正比，即：Qpf = C(H-Hf) (6-25)式中，Q pf 为单位体积的含水层在单位时间内从孔隙流入裂隙的水量，C 是比例常数；(4) 含水层骨架可以压缩，但其固体颗粒的压缩性忽略不计，看成是刚性的。6.2.2 微分方程的建立根据上述假定和水流连续性原理，可以建立裂隙承压含水层的微分方程。仍然考虑图1-29 中那个单元体的
15、水均衡，所取坐标轴和主渗透方向一致。根据第二个假定，这个单元7体中的水通过裂隙流入、流出，并服从 Darcy 定律。因此，在t 时段内沿 x 轴方向流入这个单元体与流出这个单元体的水量之差(净流入这单元体的水量) 为：式中， fxK为裂隙的主渗透系数(与 x 轴平行) 。同理可得， t 时段内沿 y 轴方向和 z 轴方向净流入这个单元体的水量分别为：和式中，fyK、 fz为裂隙的主渗透系数(分别与 y 轴和 z 轴平行 )。根据假设，孔隙中释放出的水要进入裂隙，其量为：因此，在t 时间内单元体中总的水量变化(净流入量) 为：（6-25）这个时间内，单元体内由于贮存的变化所引起的水量变化为：式中
16、 fs为裂隙贮水率。单元体内水量变化必然引起贮存的变化，两者应相等。于是得（6-26）同理，对孔隙有（6-27）式中， xK、 y、 z为沿 x，y，z 轴方向孔隙的主渗透系数， s为岩块的贮水率。由于孔隙中的水力坡降很小，可以认为：故有（6-28 ）式中 sC，称为承压水迁移系数。为了建立关干 Hf 的方程，把 (6-28)式改写为这是关于（HH f）的一阶线性微分方程。按已知公式，可求得其道解为：8（6-29）再利用初始条件：带入上式，可得 C=0。有此求得其解为（6-30）把它代入(6-28)式可得：（6-31）把它代入(6-26)式，便可得只包括裂隙水头 Hf 的方程：（6-32）上式
17、是描述承压双重介质裂隙水流的基本微分方程。在二维情况下可简化为：（6-33）式中 fxT、fy为裂隙的主导水系数 (分别与 x, y 轴平行)。 *,f为裂隙和岩块的贮水系数。把上述方程和描述多孔介质中渗流的基本微分方程比较，其不同之处只是多了一项：（ 二维情况下为 ）从前面的讨论中不难发现，它表示：因此，它的物理意义为单位时间内单位体积含水层(二维情况下为单位面积的柱体) 中从孔隙流入裂隙的水量。它是一个和时间有关的量，抽水早期，即 t 值很小时，它很小，抽出的水主要来自裂隙内水的释放，从而造成裂隙水头的迅速下降。随着时间的增长，这一水量相应地增大，裂隙水头的下降速度也随之减缓。可见，多孔岩
18、块中的水是逐渐释放出来的。这是由该数学表达式的性质决定的，从而造成孔隙水头的下降落后于裂隙水头的下降，在时间上存在着迟后。因此，这部分水量也可以称为延迟弹性释水量。随着时间的增长，迟后效应逐渐变小，孔隙中释放的水量逐渐跟得上裂隙中水位的下降，最后迟后效应小到可以忽略不计。和潜水 Boulton 方程比较，不难发现，两者在形式上是相似的，都包含有延迟效应项，即延迟弹性释水量项和潜水迟后重力排水项。但前音有明确的物理含意。延迟弹性释水项中包含一个新的参数承压水迁移系数。根据定义，其中比例常数C 是反映孔隙和裂隙之间水量交换特征的参数，与孔隙的渗透系数 K 及多孔岩块的几何特征有关。T DStrel
19、tsova(1976)认为，CK/L ，L 为岩块的特征长度，用岩块的平均大小9或岩块中心到它表面的平均距离来表示。所以值取决于 K，L， s等，它是反映孔隙、裂隙发育情况及其连通程度的特征量。越小，从孔隙向裂隙运移的水量越少，延迟时间越长；反之，越大，从孔隙向裂隙运穆的水量越多，延迟的时间越短。对无压含水层，相应近似的裂隙水流基本方程，在二维情况下为：（6-34）式中， *Ch，称为水迁移系数，h 为含水层厚度。方程右端第二项的物理含意为延迟弹性释水量与延迟重力排水量之和。对于完整抽水井，如果含水层是均质、等厚的承压含水层，抽水前，所有裂隙和孔隙中的静水压力相等，以定流量抽水，则有：（6-3
20、5 ）式中，s f 为裂隙中的降深，此时有定解条件：s (r，0)=0s (，t)=0 t 0ffrTQ2lim0t 06.3 水动力弥散理论随着近年来地下水遭到不同程度的污染，地下水中溶质运移理论可以用来模拟地下水中污染物的运移过程，预测地下水污染的发展趋势，控制地下水污染等方面。6.3.1 水动力弥散现象及其机理先考察一个实例。通过它们可以大致了解水动力弥散现象是怎么回事。例：若在一口井中瞬时注入某种浓度的一种示踪剂，则在附近观测孔中可以观察到示踪剂不仅随地下水流一起位移，而且逐渐扩散开来，超出了仅按平均实际流速所预期到达的范围，并有垂直于水流方向的横向扩散，不存在突变的界面。上述事实说明
21、，存在一种特殊的现象。因为如果不存在这种现象，示踪剂应按水流的平均流速移动；含示踪剂和不含示踪刑的水的接触界面应该是突变的；示踪剂也不应公横向扩展开来，即有一个以实际平均流速移动的直立锋面。以上事实说明，在两种成分不同的可以混溶的液体之间存在着一个不断加宽的过渡带。这种现象称为水动力弥散。因此，所谓水动力弥散就是多孔介质中所观察到的两种成分不同的可混溶液体之间过渡带的形成和演化过程。这是一个不稳定的不可逆转的过程。水动力弥散是由溶质在多孔介质中的机械弥散和分子扩散所引起的。兹分述如下。（1）机械弥散在多孔介质中，无论液体运动速度的大小还是方向，都是很不均一的。这主要和下列情况有关：由于液体有粘
22、滞性以及结合水对重力水的摩擦阻力，使得最靠近隙壁部分的(重力)水流速度趋近于零，向轴部流速逐渐增大，至轴部最大，孔隙的大小不一，造成不同孔10隙间轴部最大流速有差异，孔隙本身弯弯曲曲，水流方向也随之不断改变，因此对水流平均方向而言，具体流线的位置在空间是摆动的。这几种现象是同时发生的，由此造成开始时彼此靠近的示踪剂质点群在流动过程中不是一律按平均流速运动，而是不断向周围扩展，超出按平均流速所预期的扩展范围。沿平均速度方向和垂直它的方向上，都可以看到这种扩展现象。液体通过多孔介质流动时，由于速度不均一所造成的这种物质运移现象称为机械弥散。（2）分子扩散分子扩散是由于液体中所含溶质的浓度不均一而引
23、起的一种物质运移现象。浓度梯度使得物质从浓度高的地方向浓度低的地方运移，以求浓度趋向均一。因此，即使在静止液体中也会发生分子扩散，使示踪剂扩散到越来越大的范围。分子扩散使同一流束内的浓度趋于均一，而且相邻流束间在浓度梯度的作用下也有物质交换，导致横向浓度差的减小。物理学的知识告诉我们，分子扩散服从 Fick 定律。该定律揭示了溶液中溶质的扩散，在单位时间内通过单位面积的溶质质量入与该溶质的浓度梯度成正比，即：式中： 为该溶质在溶液中的浓度 c 沿方向 s 变化的浓度梯度，比例系数 Dd 称为扩散系数，量纲为L 2T-1。不同溶质的扩散系数各不相同，同一物质在不同温度下的扩散系数也不同。在浓度低
24、的情况下，可以认为它是一个与浓度无关的常数。由于扩散是沿着浓度减小的方向进行的，而扩散系数总是正的，所以式中要加一负号。液体在多孔介质中流动时，机械弥散和分子扩散是同时出现的，事实上也不可分。这种划分带有某种人为的性质。事实上， “纯”机械弥散不可能存在。因为当示踪剂质点沿着微小的流管运移时，分子扩散不仅使流管中的浓度趋于拉平，而且还使示踪剂质点从一条流管移向相邻的另一条流管，导致横向浓度差的减小。但分子扩散，即使在没有水流运动的情况下也能单独存在。当流速较大时，机械弥散是主要的；当流速甚小时，分子扩散的作用就变得很明显。显然，机械弥散和分子扩散都会使溶质既沿平均流动方向扩展又沿垂直于它的方向
25、扩展。前者称为纵向弥散，后者称为横向弥散。除了机械弥散和分子扩散外，某些其它现象也会影响多孔介质中溶质的浓度分布，如多孔介质中固体颗粒表面对溶质的吸附、沉淀，水对固体骨架的溶解及离子交换等。此外，液体内部的化学反应也可导致溶质浓度的变化。6.3.2 水动力弥散系数由于多孔介质几何结构的复杂性，从微观水平上研究一个点的运动规律实际上足不可能的；同样，从微观水平宋研究弥散也是困难的。因此，和定义渗流速度一样，也从宏观上来描述弥散现象，亦即将其定义在典型单元体(REV)上的平均值。分子扩散服从 Fick 定律，通过实验和理想模型的研究，证实机械弥散也能用这个定律来描述。根据 Fick 定律，多孔介质
26、中的分子扩散可用下式描述：I Dgradc 式中，D”为多孔介质中的分子扩散系数，量纲为L 2T-1 ，是二秩张量；c 为该溶质在溶液中的浓度；I为由于分于扩散在单位时间内通过单位面积的溶质质量，对于机械弥散有：11I一 Dgradc 式中，D为机械弥散系数，量纲为L 2T-1 ，也是二秩张量； I为由于机械弥散造成的个单位时间内通过单位面积的溶质质量，c 的含义同前。D 和 D的量纲相同，由此定义水动力弥散系数 D：D 也是二秩张量。由于水动力弥散在单位时间内通过单位面积的溶质的质量则为 II 十I一 Dgradc，如果我们选择 x 轴与该点处的平均流速方向一致，y 轴和 z 轴则与平均流速
27、方向垂直，则上式也可以写成下列更容易被我们理解的形式：或此时水动力弥散系数张量：坐标轴方向称为弥散主轴。D xx 称为纵向弥散系数，D yy， Dzz 称为横向弥散系数。由于弥散主铀的方向依赖于流速方向，即使在均质各向同性介质中，各点弥散主轴的方向也会随着水流方向的改变而各不相同。水动力弥散系数在研究地下水物质运移问题中是一个很重要的参数。通过大量在末固结的多孔介质中的实验，得到了如图 7-10 所示的曲线。图中纵坐标是从实验室得到的纵向弥散系数 DL 与溶质在所研究的液相中的分子扩散系数 Dd 的比值，横坐标是一个无量纲的量：12图 6-10 分于扩散和水动力弥散间的关系 (据 J. Bea
28、r)称为 Peclet 数。其中，u 为实际平均流速，d 为多孔介质的某种特征长度，如多孔介质的平均粒径等。该无量纲数表示实际流速和分子扩散系数相比的相对大小，Pe 数愈大，表示流速相对愈大。根据这条曲线的变化情况，大致上可以分五个区。第 I 区：实际流速很小，以分子扩散为主，相当于曲线上寻接近于常数的一段。第 II 区：对应的 Peclet 数 Pe 约在 0.4 到 5 之间，曲线开始向上弯曲，机械弥散已达到和分子扩散相同的数量级。因此，应当研究两者的和，而不应忽略其中的任何一个。第 III 区：物质运移主要由机械弥散和横向分子扩散相结合而产生。横向分子扩散往往会削弱纵向的物质运移，实验结
29、果得出 DL/Dd (Pe)m， =0.5，1M1.2。第 IV 区：以机械弥能为主，分子扩散的作用已经可以忽略不计，但流速尚未达到偏离 Darcy 定律的程度。本区相当于图中的直线部分。实验给出于 DL/DdPe，=1.8。第 V 区；仍属于机械弥散为主的区域，与第 IV 区的区别在于水流速度已达到越出Darcy 定律适用的范围。惯性力和紊流的影响造成纵向物质运移的减少，曲线斜率减缓。横向弥散试验得到了和纵向弥散相类似的结果。上述曲线说明，弥散系数和水流速度、分子扩散有关。它们间的关系如下式所示：式中：ijD机械弥散系数，为一个二秩对称张量，这是它的一个分量；kmij,多孔介质的弥散度，为一
30、四秩张量；在饱和流动中它反映多孔介质固体骨架 的几何性质，量纲为L；13u实际平均流速，u k，u m 分别为它在坐标轴 xk、x m 上的分量；表示水流通道形状持征的系数，无量纲；在微观水平上考虑相邻流线之间内分子扩散所引起的对物质运移影响的因数，这个影响和机械弥散是不可分的。Pe 较大时，由 f(Pe, )的表达式可以看出，f(Pe, )1。也就是说，分子扩散对机械弥散系数的影响就变得微不足道了。这时机械弥散系数和实际平均流速之间呈线件关系。对于大多数实际问题来说，都属于这种情形，总是假定，f(Pe, )1。如果在某一点上选择坐标轴，使得其中一个坐标油(如 f 轴 )祁该点处的平均流速方向
31、一致(即弥散主轴)，并忽略分子扩散， ，f(Pe, )1，则：式中， L， T 分别称为纵向弥散度和横向弥散度。纵向机械弥散系数xD和横向机械弥散系数yD，及 z称为弥散系数的主值。由于弥散主轴依赖于水流方向，所以除了均匀流（u x常数，u y =uz =0）以外，一般说来即使在各向同性介质中各点的弥散系数也各不相同，随空间位置而变化。6.3.3 对流弥散方程及其定解条件考虑由某种溶质和溶剂组成的二元体系。以充满液体的渗流区内任点 p 为中心，取一无限小的六面体单元，各边长为x、 y 和z，选择 x 轴与 p 点处的平均流速方向一致，来研究该单元中溶质的质量守恒。先研究由水动力弥散所引起的物质
32、运移。在t 时间内，由于弥散和水流运动所引起的单元体内总的溶质质量变化为在t 时间内，单元体内溶质的浓度发生了 tc的变化，单元体内的液体体积为 nxy z，则由它所引起的该单元体中溶质质量的变化为： tzyxtcn根据质量守恒定律，上述两者应该相等，当坐标轴与水流平均流速方向一致时，经过整理、简化后得到：上式称为对流一弥散方程(水动力弥散方程 )。它右端后三项表示水流运动 (习惯地把它喻为对流)所造成的溶质运移，前三项表示水动力弥散所造成的溶质运移。如果还有化学反应或其它原因所引起的溶质质量变化，且单位时间单位体积含水层内由此而引起的溶质质量的变化为 f，则应把它加到方程式的右端，有：要确定
33、一个水动力弥散问题的解，即求得浓度的分布，还要给出下列信息：研究空间 和时间区间 0，T；研究区域水头场的分布；有关参数，如弥散度 L 和 T 等；14定解条件。初始条件给出初始时刻(t0)区域 上的浓度分布，即：c(x,y,z,0)= c0(x,y,z)c0 是已知函数。边界条件通常有二种类型。一种是已知浓度的边界条件，即： 11 z)y,(x|t,( |t)z,y(x式中， 1为研究区的边界， 是已知函数。另一种是已知单位时间内通过边界单位面积的溶质质量的边界条件。在三维条件下，形式复杂，不易理解。兹以一维问题的几种常见例子具体说明如下。(1) 多孔介质 a 的边界外为另一多孔介质 b，根
34、据单位时间通过边界的溶质的质量要保持连续的原则，当渗透速度为 v 时有：(2) 如边界为隔水边界，则通过边界的流量和溶质的量均为零，由上式0xcDvL及 v0 得边界 2上有边界条件：|16.3.4 一维弥散问题的解考虑流速方向与 x 轴方向一致的半无限一维均匀流的情况，示踪剂连续注入，纵向弥散系数 DxxD L 在均匀流情况下不随坐标 x 而变化，u x u 为常数，一维情况下(749)式化为：同时有定解条件：当麦足够大时，该定解问题的解为利用(757) 式可以求得任意时刻，任意距离 f 处的相对浓度 cn。因为示踪剂浓度c。是已知的，即可求得该处的浓度 f(2， )。反之，也可利用实验室或
35、野外的一维弥散的实际观测资料，求出纵向弥散系数几，因为流速“已知，也可以算出纵向弥散度吨。根据对流弥散方程，在适当的初始条件、边界条件下求得的解，可以用来预报地下水中污染物的时、空分布。其结果和实验室的实验结果，一般也拟合得很好。但应用于野外试验时，却发现利用对流弥散方程反求得的弥散度值要比实验室实验所得的值大几个数量级，而且弥散度值看来和污染物分布的范围有关，随着它的增大而增大(称为尺度效应) 。15对产生这种观象的较普遍的看法是，是受岩层非均质性影响的结果。非均质性引起复杂的速度分布。这种分布导致在一定程度上与均质介质中因粒间速度变化引起的机械弥散现象相类似的污染物分布。虽然非均质性是造成上述现象的主要原因已没有多少疑问，但这种影响的方式还没有了解得很清楚。总之，虽然过去 20 多年中对地下水中溶质运移的研究已经取得了很大成绩，但这个领域的研究还只是处于它的“幼年阶段” 。还需要大量的室内和野外试验，以便积累在野外尺度上溶质运移过程的知识，为发展其数学模型提供基础。没有这些，就难以充分地评价人类活动对地下水环境质量的影响。
