第七章 地下水运动中的若干专门问题 10页

第七章地下水运动中的若干专门问题§1非饱和带的地下水运动一、关于非饱和带水分的基本知识1.含水率，饱和度和田间持水量包气带中的空隙，一部分被水充填，另一部分被空气充填。含水率(θ)：表示单位积中水所占的体积，式中：（Vw）0——典型单元体中水的体积；V0——典型单元体的体积饱和度：岩石的空隙空间中被水占据部分所占的比例。式中：（V0）0——典型单元体中的空隙体积含水率与饱和度的关系：θ=nSw式中：n——孔隙度。田间持水量：在长时间重力排水后仍然保留在土中的水量。2.毛管压力毛管压强：在多孔介质的孔隙中，液体和气体接触是，二者存在压力差，这个压力差称毛管压强。用pc表示pc=pa-pw式中：pa——空气的压强；pw——水的压强毛管压强取决于界面的曲率，曲率愈大（液面愈弯曲，毛管压强愈大。以上毛管压强是以绝对压强为基准，如果以相对压强为基准，这时有：pc=pa-pw–pa∴pc=-pw毛管压强相对大气压强为负值。即，非饱和带孔隙中的水处于小于大气压强的情况下。非饱和带水流中任何点的水头式中：z——位置水头；hc=pc/r——毛管压力水头；∴H=z-hc3土壤水分特征的曲线\n水分特征曲线：反映毛管压力水头（或毛管压强）和土壤含水率或饱和度关系的曲线。如图：随着含水率的减少，毛管压力增加，当含水率减小到某一值时，压强继续增大时，含水率不在减小。相应的饱和度为：影响特征曲线的因素：（1）不同质地的土壤，其水分特征曲线不同。一般说，土壤的粘粒含量愈高。同一负压条件下土壤的含水率愈大，或者同一含水率下其负压愈高。这是因为，粘粒含量增多。使土壤中细小孔隙发育的缘故。（2）土壤结构。如图，为一砂壤土不同干容重的水分特征曲线，在同一负压下，土壤愈密实，（大），相应的含水率一般也大。原因，土壤愈密实，大孔隙数量减少，中孔隙增多。（3）温度的影响。温度升高，水的粘滞性下降，所以表面张力降低，在同样的负压下，含水率要低一些。（4）土壤水分变化过程的影响。对于同一土壤，土壤脱湿（由湿变干）过程测得的水分特征曲线不同，如图，在相同的负压下，排水（脱湿）时的含水率要大于吸湿时的含水率。这种现象称为滞后现象。（5）容水度：毛管压力水头变化一个单位时，从单位体积土中释放出的水体积。数值上等于，水分特征曲线的斜率的负倒数。它是含水率和毛管压强的函数，可用或表示。二、非饱和带水运动的基本方程非饱和带地下水的运动，也可以用达西定律描述，这时的渗透系数是变化的，与含水率和毛细压力水头有关，是和的函数，其关系如图，随着含水率的增大，渗透系数增大，随毛管压力水头的减小，渗透系数增大。在非饱和带中，定律的表达式为：在三个坐标轴的分量为：2基本微分方程第一章推的渗流连续性方程，如下：在饱水带中，全部孔隙被水充满，等式右端用孔隙度，在非饱和带中，部分孔隙被水充满，所以用含水率取代，并两边除（近似为常数），得：\n将vx、vy、vz代入上式，得：二式为非饱和流的基本方程3基本方程的几种形式（1）以含水率为因变量的表达式将H换成θ，将H=z-hc代入上（1）式，得：上式进一步变换定义K（θ）/C（θ）=D（θ）为扩散系数，得：对于垂向一维流动，去掉前两项，得：轴向上取正值，轴向下取负值。（2）以毛管压力水头为因变量得表达式：将H=z-hc代入（2）式，得：代入上式：对于垂向一维流动，去掉前两项，得：轴向上取正值，轴向下取负值。§2水动力弥散理论用来模拟地下水中污染物和化学成分得运移过程，预测地下水污染得发展趋势。一水动力弥散现象及其机理\n例1.在一口井中注入一种示踪剂，示踪剂在随地下水向前流动得过程中，向外围扩散，形成一个以中心点浓度最大，向四周浓度逐渐减小的过渡带，并且随示踪剂迁移的距离增大，过渡带也越来越宽。如图（书中）例2.在均匀流的砂柱中，用含有示踪剂浓度为的水去替代，在砂柱另一端测量示踪剂浓度，得曲线如图（书）。水在流动过程中并非一个突变界面，而是一个过液带。这种现象称水动力弥散。水动力弥散是机械弥散和分子扩散所引起的。1机械弥散液体在多孔介质中运动的三种情况：（1）由于液体粘性的作用和结合水的摩擦阻力，使得靠近孔隙壁的水流速度趋于零。孔隙中心部位流速最大。孔隙大小不一，造成不同孔隙之间沿轴部的最大流速有差异；（1）由于空隙的弯弯曲曲，水流方向也随之不断地改变。由于上述三种情况，造成了地下水质点运动速度，在大小和方向上的不均一，造成了示踪剂有的运动快，有的运动慢，从而形成了上述过渡带。这种由于速度不均一所造成的这种物质运移现象称为机械弥散。2分子扩散一般溶质都有由浓度高向浓度低的地方运移的性质，以求浓度趋于均一。这种由于液体中所含溶质的浓度不均一而引起的物质运移现象叫分子扩散。分子扩散服从定律式中：Is——单位时间内通过单位面积的溶质的质量；dc/ds——溶质在溶液中的浓度c沿s方向变化的浓度梯度；Dd—扩散系数。机械弥散和分子扩散是同时出现的，当流速较大时，机械弥散是主要的；当流速甚小时，分子扩散的作用就变得明显。水动力弥散还分为沿水流方向和垂直与水流方向的弥散，沿水流方向的弥散称纵向弥散，垂直水流方向的弥散称横向弥散。二水动力弥散系数分子扩散服从定律：式中：D″——为分子扩散系数；I″——为由于分子扩散在单位时间内通过单位面积的溶质质量。机械弥散也服从定律：式中：D′——为机械扩散系数；I′——为由于机械扩散在单位时间内通过单位面积的溶质质量。由于水动力弥散是分子弥散和机械扩散共同作用的，定义水动力弥散系数：水动力弥散定律如下：式中：I——单位时间内通过面积的溶质质量；D——水动力弥散系数；Dc/ds——浓度梯度。\n如果我们取方向与流速方向一致，轴和轴与流速方向垂直，上式可用下式表示：三对流弥散方程及其定解条件如图，以渗流区内任一点为中心，取一无限小的六面体单元，各边长为Δx、Δy、Δz，选择x轴与P点处的平均流速方向一致。（即纵向弥散方向为轴方向）在对流弥散问题中，包括两个子问题：其一，溶质随地下水的流动或流出单元体；其二，溶质通过自身的弥散流入或流出单元体。水动力弥散引起的物质运移：设，沿轴方向溶质的质量变化率为dIx/dx，如果假设Ix为在abcd面上，单位时间内通过单位面积溶质的质量，那么，时间内通过面流入单元体的溶质质量为：IxnΔyΔzΔt因为沿轴方向溶质的质量变化率为dIx/dx，所经距离Δx后，变化了（dIx/dx）Δx，所以，Δt时间内，通过a′b′c′d′面流出单元体的溶质质量为：所以，沿x轴方向流入与流出单元体的溶质质量差为：同理，沿y轴方向和z轴方向流入与流出单元体的溶质质量差为：所以，通过弥散单元体内溶质质量的变化为：随地下水流的物质运移：设沿x轴方向，在abcd面地下水的流速为vx，则单位时间通过abcd面单位面积流入单元体的水量为vx×1×1，流入的溶质质量为vx×c：（为溶质的浓度），那么，在Δt时间内流入abcd面溶质量为：vxcΔyΔzΔt设沿方向，通过单位面积溶质质量的变化率为：经Δx距离的变化量为：在a′b′c′d′面，单位时间单位面积流出单元体的溶质质量为：在Δt时间内流出面的溶质质量为：\n所以，沿轴方向流入与流出单元体的溶质的质量差为：同理，沿轴和轴方向流入与流出单元体的溶质质量差为：所以，随地下水流流入与流出单元体的溶质质量差为：Δt时间内，流入和流出单元体总的溶质质量差为：另外，设单元体内溶质浓度随时间的变化率为dc/dt，那么，Δt时间内单元体内体积溶质浓度变化量为：所以，Δt时间内单元体内溶质质量变化量为：上述二量应相等，并消去ΔxΔyΔzΔt，得：代入上式，并两边同除以，则得：式中，u——实际流速。上式为对流弥散方程。如果，有其它源汇项时，并设单位时间单位体积含水层内由源汇引起的溶质质量的变化量为f。在上式的左边加一项f即可。关于溶质运移的数学模型除微分方程外，还应有定解条件：初始条件：初始时刻的浓度分布。表达式如下：C(x，y，z，0)=C0(x，y，z)边界条件：有两类。一类是已知浓度的边界条件，表示如下：式中：Γ1——表示一类边界。另一类是通量边界，即单位时间内通过单位边界面积的溶质质量已知。如：隔水边界：\n补给边界：所以，要确定一个水动力弥散问题的解，即求得浓度的分布，要给出下列信息：（1）微分方程；（2）研究空间区域和时间区域；（3）研究区域水头场的分布；（4）有关参数，如弥散度和等；（4）定解条件。四一维弥散问题解设投放示踪剂前，含含水层中示踪剂的浓度为0，然后在河渠中连续注入浓度为C0的示踪剂。在均匀流情况下，ux=u为常数（1）数学模型：定解问题的解：其中当时，上解可近似为：（2）利用实验资料求纵向弥散系数有一个观测孔时，观测孔距河渠距高已知，在观孔中可测得不同时刻的，从而可求得。求参步骤：①据实验资料作Ci/C0—t关系曲线，如图②在图上找出Ci/C0为0.84和0.16二点，并读出其横坐标t0.84和t0.16。③代入下式求纵向弥散系数式中：u——地下水实际流速§3海岸带含水层中的咸淡水界面天然条件下，在海岸带含水层中的地下水一般是流向海的，由于海水比淡水的比重大，海水体将位于淡水体的下方，是楔型，如图。处于平衡状态。当在海岸边抽取淡水时，这时淡水的水位下降，打破了原来的平衡，引起海水向内陆的入侵，以达到新的平衡，这时界面向陆地推进。该现象为海水入侵。海水与淡水是可以溶混的，由于水动力弥散，在海水与淡水之间形成了一个过渡带，在过渡带中地下水的矿化度由小变大，直到海水矿化度。\n过渡带的宽度在不同的地区，其宽度不同，当其宽度较小，且与含水层的厚度相比较小时，可以以为海水与淡水之间是一突变界面；否则宽度较宽时，则作为水动力弥散问题加以研究。一作突变界面处理静止界面的近似解当淡水和海水处于一种平衡状态时，界面是静止的，如图，假设淡水的容重为rf，海水的容重为rs。在界面上的点，受海水的压力为：γshs在界面上的点，受到淡水的压力为：γf（hs+hf）此二压力应相等：γf（hs+hf）=γshs解得：令则：hs=δhf一般海水密度为1.025g/cm3，容重γs=10045N/m3，淡水密度为1.000g/cm3，容重γs=9800N/m3。代入可求得：δ=40∴hs=40hf说明：在离海岸任一距离上，稳定界面在海面以下得深度为该处淡水高出海面得倍。以上仅是一种近似解法。二确定界面的形状及海水入侵的范围（一）厚度固定的水平承压含水层中的界面问题水流是稳定流，如图，设，原点位于坡脚（点），x轴的正向指向海。含水层厚度为M，承压水头为H，假设地下水为水平流，忽略垂向分速度。由Darcy定律，有：K=Kf（Kf含水层中淡水的渗透系数）H=Hf（Hf含水层中淡水的水头）由前面知，hs=d+h(x)∵hs=δH∴d+h(x)=δH两边对求导，得：（与无关）代入上式得：对上式整理得：两边积分得：\n当x=0时，h(x)=M，代入上式得：代入上式得：上式表明界面得形状是一条抛物线。利用此式可确定X处的h（x）。另外，由d+h(x)=δH得：h(x)=δH-d代入（1）式得：两边积分，得：当x=0时，H=H0，代入上式，得：代入，得：此式表明承压水含水层淡水的水头面的形状也是一条抛物线。当x=L时，h=0，这时d+h=d=δH∴H=d/δ代入上式，得：此式表示出海水入侵深度与流向海的淡水流量q0和界面坡脚以上测压水头H0之间的关系。又∵在x=0处，有d+M=Δh0，代入上式，得：此式中，均为常数，所以等式的右端为常数，这样可得出：与是反比。当增大时，减小。（一）潜水含水层中的界面问题含水层上部有均匀入渗补给，水流是稳定流，且假设水流为水平流动。如图。处的单宽流量为：前面得出了h=δhf代入上式，得：\n分离变量，得：两边积分得：当x=0时，hf=H0，代入上式得：则：此式表明：界面形状为一条椭圆线。当x=L时，hf=0，有：此式表示与的关系。即越大，越小，所以，可以通过控制来控制海水入侵的范围。返回上页