水动力弥散方程

第三章 求解水动力弥散方程的有限差分法在符合定解条件的情况下，水动力弥散方程的解析解，可求得计算区域任一点的溶质浓度及其随时间的变化，可分析多种因素对溶质运移过程的影响，利用实测资料可反求溶质运移有关参数，还可用以验证数值方法的可靠程度等。

但对于较复杂的初边值条件或非稳定水流运动中的溶质运移问题，一般很难求得解析解。

即使通过一定简化求得解析解，也由于公式过于繁琐以致难以实际应用。

实践中，对这类问题，常用数值方法求解。

求解水动力弥散方程的数值方法，主要有有限差分法和有限单元法。

考虑到前面有关章节对这两种数值方法已做过详细介绍，三、四两章将主要讨论几种对求解水动力弥散方程较为有效和实用的方法第一节 有限差分计算中出现的跳动和数值弥散根据前面章节所述差分法的一般概念，可以直接写出水动力弥散方程的差分方程。

一般情况下，一维弥散问题隐式格式可用追赶法求解，二维弥散问题常用交替方向隐式差分格式（ADI 法）求解。

具体解法与一般地下水运动差分方程的解法类同，此处不再赘述。

但值得注意的问题是，由于对流项的存在和差分过程中截断误差的影响，使水动力弥散问题的一般差分解法常遇到两个数值困难，即跳动（过量）和数值弥散。

所谓跳动是指在浓度峰面处，用数值法求得的相对浓度大于或小于l ；而数值弥散是指在浓度峰面处，数值法求得浓度的曲线形状比实际界面平缓。

为了解释数值跳动和数值弥散产生的原因，现以均质各向同性含水层一维弥散问题为例加以说明]7[。

此时水动力弥散方程为x cv xc D t c ∂∂-∂∂=∂∂22 （3-3-1） 根据差分法的基本概念，可直接写出其相应隐式差分方程为xC C vx C C CD tC C k i k i k i k i k i k i k i∆--∆+-=∆-+-++++++-+2)2(11112111111（3-3-2）式中，1+k iki C C 、分别为时段始、末的结点i 的溶液浓度。

令2x tr ∆∆=，上式整理后为 ki k i i k ii k i i C rC D C B C A 111111=++++++- （3-3-3） 式中，。

水动力弥散方程水体中的物质运移和扩散往往会受到水流的影响，因此涉及到水动力弥散方程。

水动力弥散方程是描述物质在水动力作用下在水体中弥散和扩散的方程。

在环境保护、污染防治、水资源利用和水力工程等领域中，水动力弥散方程非常重要。

弥散的基本概念在介绍水动力弥散方程之前，需要先了解一些基本概念。

弥散弥散是指物质在水中因为分子热运动而发生的无规则传递过程。

在水中，物质均呈现出弥散的现象，即物质会沿着水流的方向不断扩散。

扩散扩散是指物质在稳定均匀的介质中自发地运动，使得物质的浓度分布趋于均匀的传递过程。

对流对流是指流体中由于温度等差的非均匀性而引起的流动。

水动力域中，对流一般指水流的流动。

分子扩散分子扩散是指物质在介质中因分子热运动而发生的扩散过程。

水动力弥散方程的构建在水动力弥散方程中，要考虑物质的对流和扩散。

如果仅考虑扩散，则十分简单，其方程为：$$\\frac{\\partial c}{\\partial t}=D\ abla^2c$$其中，c表示物质的浓度，t表示时间，D为扩散系数。

但实际上，流体内部还会存在对流影响，所以在含有对流的情况下，水动力弥散方程为：$$\\frac{\\partial c}{\\partial t}+v\ abla c=D\ abla^2c$$其中，v表示水流的速度。

这个方程告诉我们随着时间的推移，浓度c会发生变化。

变化是由扩散和对流两种机制引起的，从而影响水体中物质的分布情况。

水动力弥散方程的本质意义是用数学语言描述了物质在水动力作用下如何弥散和扩散。

物理解释物理上，扩散作用是由分子的玻尔兹曼方程描述的，而对流作用是由沃滕变换描述的。

弥散过程是扩散和对流两种作用的综合体现。

在弥散过程中，对流所起的作用是将物质从一处地方迅速“输送”到其他地方，从而影响弥散的速率。

对流作用越强，同样的物质浓度分布会更快地发生变化；反之，扩散作用相对于对流影响变弱，则物质的浓度分布变化更缓慢。

地下水动力学概念总结---- King Of Black Spider 说明：带下划线的是重点，重点116个，次重点22个，共138个。

第0章地下水动力学：Groundwater dynamics研究地下水在孔隙岩石、裂隙岩石和岩溶（喀斯特）岩石中运动规律的科学，它是模拟地下水流基本状态和地下水中溶质运移过程，对地下水从数量上和质量进行定量评价和合理开发利用，以及兴利除害的理论基础。

主要研究重力水的运动规律。

第1章渗流：Seepage flow是一种代替真实地下水流的、充满整个岩石截面的假想水流，其性质（密度、粘滞性等）与真实地下水相同，充满整个含水层空间（包括空隙空间和岩石颗粒所占据的空间），流动时所受的阻力等于真实地下水流所受的阻力，通过任一断面及任一点的压力或水头均与实际水流相同。

越流：Leakage 当承压含水层与相邻含水层存在水头差时，地下水便会从水头高的含水层流向水头低的含水层的现象。

对于指定含水层来说，水流可能流入也可能流出该含水层。

贮水系数：storativity又称释水系数或储水系数，指面积为一个单位、厚度为含水层全厚度M的含水层柱体中，当水头改变一个单位时弹性释放或贮存的水量，无量纲。

μ* = μs M。

既适用于承压含水层，也适用于潜水含水层。

导水系数：Transmisivity 是描述含水层出水能力的参数；水力坡度等于1时，通过整个含水层厚度上的单宽流量；亦即含水层的渗透系数与含水层厚度之积，T=KM。

它是定义在一维或二维流中的水文地质参数。

单位：m2/d。

非均质介质：如果在渗流场中，所有点不都具有相同的渗透系数，则称该岩层是非均质的。

各向异性介质：渗流场中某一点的渗透系数取决于方向，渗透系数随渗流方向不同而不同。

达西定律：Darcy’s Law 是描述以粘滞力为主、雷诺数Re< 1~10的层流状态下的地下水渗流基本定律，指出渗流速度V与水力梯度J成线性关系，V=KJ，或Q=KAJ，为水力梯度等于1时的渗流速度。

地下水水质的数学模擬(三)——水动力弥散方程的解析解法及其应
用
地下水水质的数学模拟是地下水地下水水质保护的重要方法之一。

在地下水水质模拟中，水动力弥散方程是一个重要的方程，可以用来模拟地下水的水流和污染物的扩散。

下面是水动力弥散方程的解析解法及其应用:
一、水动力弥散方程的解析解法
1. 欧拉法
欧拉法是一种经典的求解水动力弥散方程的方法。

该方法的主要思想是将水动力弥散方程转化为一个积分方程，然后通过欧拉方法来求解积分方程。

欧拉法的基本思路是将时间域问题转化为频域问题，并使用频率分析方法来求解。

2. 拉格朗日法
拉格朗日法是一种基于拉格朗日平动理论的解析方法。

该方法的主要思想是将水动力弥散方程转化为一个拉格朗日方程，然后通过拉格朗日方程来求解水动力弥散方程。

拉格朗日法适用于求解非线性水动力弥散方程。

二、水动力弥散方程的应用领域
1. 地下水污染控制
水动力弥散方程可以用来模拟地下水的水流和污染物的扩散，从而帮助人们掌握地下水的污染状况，并为地下水污染控制提供科学的决策支持。

2. 水文地质勘探
水动力弥散方程也可以用来求解水文地质勘探中的勘探参数，从而帮助人们掌握地下水的分布情况，为水文地质勘探提供科学的决策支持。

水动力弥散系数一、基本概念在研究地下水溶质运移问题中，水动力弥散系数是一个很重要的参数。

水动力弥散系数是表征在一定流速下，多孔介质对某种污染物质弥散能力的参数，它在宏观上反映了多孔介质中地下水流动过程和空隙结构特征对溶质运移过程的影响。

水动力弥散系数是一个与流速及多孔介质有关的张量，即使几何上均质，且有均匀的水力传导系数的多孔介质，就弥散而论，仍然是有方向性的，即使在各向同性介质中，沿水流方向的纵向弥散和与水流方向垂直的横向弥散不同。

一般地说，水动力弥散系数包括机械弥散系数与分子扩散系数。

当地下水流速较大以致于可以忽略分子扩散系数，同时假设弥散系数与孔隙平均流速呈线性关系，这样可先求出弥散系数再除以孔隙平均流速便可获取弥散度。

分子扩散系数D '与介质的性质有关。

经验证明：T D D d ⋅=' (6-25)式中 d D ——溶质在静水中的分子扩散系数，它主要取决于溶质分子的特性和温度；T ——多孔介质的弯曲度。

机械弥散系数D ''是一个与地下水流速有关的量。

在各向同性介质中，经试验证明为：UU U U D ji T L ij T ij )(ααδα-+='' (6-26)式中 ij δ——Kronecker 记号，当j i =时，1=δ，当j i ≠时，0=δ；L α——纵向弥散度； T α——横向弥散度；U ——地下水实际速度，i U 、j U 为实际速度的分量；二、水动力弥散系数确定的试验方法水动力弥散系数可通过室内或现场弥散试验确定。

弥散系数的计算方法一般分两类：一是利用解析公式直接或间接求解；二是采用标准曲线对比法。

1.一维室内弥散试验测定水动力弥散系数 （1）试验原理以人工配制的均质各向异性岩样，进行示踪剂注入实验。

具体假设及要求如下： ①. 试验流场为均质不可压缩的稳定的一维流场，渗流为定水头补给的一维弥散； ②. 多孔介质是均质的，渗透系数，孔隙度和弥散系数都是常数； ③. 流体是不可压缩的均质液体，密度、粘滞度为常数，温度不变； ④. 试验土柱（或砂柱）及其中之流体，示踪剂的初始浓度为一定值。