第2章 水动力弥散方程.

1、 第二章 水动力弥散方程用来描述地下水系统当中溶质运移规律的数学方程（微分方程）。本章主要内容有：2-1. 水动力弥散方程的有关参数1、流体的密度、浓度； 2、多组分流体的流速； 3、流体的通量。2-2. 溶液中组分的质量守恒方程2-3. 组分的的对流扩散（Fick方程）2-4. 多孔介质中水动力弥散方程2-5. 源汇项2-6. 初始条件与边界条件2-1 水动力弥散方程的有关参数2-1-1 流体的密度（）所谓的流体密度指的是单位流体体积的质量，常用 表示，量纲ML-3。多组分流体的密度实际上对于非均质的多组分流体而言，其密度是随着组成它的各种组分的浓度不同而变化的。假设某多组份流体共有N种组分

2、其某一组分称为 ，取该液体中一体积为dv的微元，其质量为dm，该液体中在dv微元中组分的质量为dm则 组分的质量密度：dvdm若将所有N种组分的质量密度进行求和： dvdmdvdmdvdmNNN1121就等于该溶体的体系密度。某一组分的质量的密度：实际上就是水化学中学过的某一组分的浓度。浓度定义为单位体积流体某种溶质的质量。2-1-2 多组分流体的流速 u组分的质点流速u是指在dv内组分的各个分子的统计平均速度，也就是各个分子的速度之和除以分子的个数。对每种多组分流体来看溶液中各种组分的速度是不相等的。流体体系的质点流速:流体体系中各组分的质量平均速度 u一般情况下， 组分的质点流速 与流体体

3、系的质量平均流速 是不相等的，两者存在一个偏差：uuuuuuuu 或u称为组分质点相对于质量平均速度 的扩散速度。u2-1-3 流体的通量流体的质量通量流体的质量通量 ：流体在单位时间内通过单位面积的流体的质量：流体在单位时间内通过单位面积的流体的质量uJ组分的质量通量组分的质量通量 单位时间内通过与流体方向垂直的单位单位时间内通过与流体方向垂直的单位面积上的面积上的组分的质量。组分的质量。uJ组分相对与溶体质量平均流速组分相对与溶体质量平均流速 的的质量扩散通量质量扩散通量 ：对流体体系来说，显然有：对流体体系来说，显然有：JJJu()JuuuJ10NJ111110NNNNNJuuuuuu这

4、是因为这是因为2-2 溶液中组分的质量守恒方程（连续介质） 在多组分组成的流体体系中任取一点P(x,y,z),以P为中心取一微小的质量平衡体（如图2-1），其侧面分别平行与3个坐标面，边长分别为 x、y、 z，质量守恒原理:在时间t内，组分在这个单元体中的净流出（或流出）量（暂不考虑起内部有质量产生和消失），应等于这个单元中组分的质量变化用方程的形式可表示为：zyxtyxuutzxuutzyuuzzyxzzzyxzzyyxyzyyxyzyxxxzyxxx2,2,2,2,2,2质量守恒方程（连续介质）zyxuuu,设设 分别表示分别表示组分密度、组分密度、x,y,z方向的速度方向的速度分量。分量

5、。其中：其中： 经过经过t时间后，质量均衡体中时间后，质量均衡体中 的变化量。的变化量。将上式左右两端同除以将上式左右两端同除以 得：得：tzyxtzuuyuuxuuzzyxzzzyxzzyyxyzyyxyzyxxxzyxxx2,2,2,2,2,2再对方程两端取极限，即令再对方程两端取极限，即令0, 0, 0, 0tzyxtzuyuxuzyx即有即有：即即0divut若微小的质量均衡体内存在着若微小的质量均衡体内存在着组分的源汇项，则上式可改写为：组分的源汇项，则上式可改写为：Iutdiv多组分流体体系中多组分流体体系中组分组分的质量守恒方程的质量守恒方程 多组分组成的流体中，多组分组成的流体

6、中，单位体积流体在单位时间内，由于化单位时间内，由于化学反应或其它原因所学反应或其它原因所产生产生（或（或消失消失）的）的组分的组分的质量质量。I 上述质量守恒方程中，至少包括上述质量守恒方程中，至少包括 4个未知变个未知变量量 有时可以独立给出（如抽、注、示踪剂的速率），但有时也有时可以独立给出（如抽、注、示踪剂的速率），但有时也与与 有关，如吸附作用、溶解作用，不能简单的有关，如吸附作用、溶解作用，不能简单的 给定，因此上述给定，因此上述方程不能单独求解，还必须引入方程不能单独求解，还必须引入通量通量与与驱动力驱动力之间的关系式，即之间的关系式，即质质量通量量通量与与组分密度组分密度间的关

7、系。间的关系。zyxuuu,I 在多组分组成的溶体体系中，一种组分的运移受两个因素的驱动在多组分组成的溶体体系中，一种组分的运移受两个因素的驱动：2-3 组分的对流扩散方程（连续介质）一一是受流体的流动的控制是受流体的流动的控制，即该组分按平均流速随这个流体体系，即该组分按平均流速随这个流体体系的运移，的运移， 即即对流对流；二二是该组分的自身分子扩散是该组分的自身分子扩散，即由，即由浓度梯度浓度梯度引起的相对于平均流速引起的相对于平均流速运移的运移的分子扩散分子扩散。uutuututIdivdivdivdiv下面在下面在组分质量守恒方程基础上建立组分质量守恒方程基础上建立组分的组分的对流对流

8、扩散方程扩散方程：引入引入组分的质量扩散通量组分的质量扩散通量 则上式可写成：则上式可写成：JIJutdivdiv 是是组分的质量通量组分的质量通量 的的对流对流分量。分量。 是是组分的质量通量组分的质量通量 的的扩散扩散分量。分量。uuJu 对于上有溶质、溶剂两种组分构成的二元体系，对于上有溶质、溶剂两种组分构成的二元体系，组分在等温组分在等温条件（忽略热扩散）相对于质量平均速度条件（忽略热扩散）相对于质量平均速度 的扩散通量的扩散通量 可依可依Fick定律得出：定律得出：uJgradmDJ 表示溶质的分子扩散系数：表示溶质的分子扩散系数：mDICDuCtCmgraddivdiv二元体系中二

9、元体系中组分的组分的对流对流扩散方程扩散方程 对于低浓度溶液，浓度对于低浓度溶液，浓度C的改变并不明显地影响的改变并不明显地影响 ，于是，于是 可视为常量，可视为常量， 也可视为常数。则：也可视为常数。则：mDCDJmgradICDuCtCmgraddivdiv稀释的二元体系中稀释的二元体系中组分的对流组分的对流扩散方程扩散方程 将上述对流扩散方程加上适当的边界条件和初始条件。即可用来解决流动的地表水中组分的分布及变化规律（例如地表水体中污染物质的迁移）。应用条件：1、二元体系二元体系；2、等温条件等温条件；3、低浓度低浓度；因此，必须对上述方程的各变量在典型单元体上取平均值，也就是从微观水平

10、上的研究过渡到比较粗的宏观水平上来研究多孔介质中所发生的现象。方程中微观变量C、 ，都是相对于流体的质点而言的，而实际工作中都是取它们在典型单元体上的平均值。u2-4 多孔介质中水动力弥散方程 上述对流扩散方程是对流体连续介质建立的，若从这种微观水平上来研究多孔介质中的溶质输运，则需把多孔介质的骨架作为问题的边界。 将速度将速度 和浓度和浓度 在典型单元体的空隙体积在典型单元体的空隙体积 上取平上取平均均 值值 和和 ：CuVV。uCvoVVVodVuVu1voVVVoCdVVC1速度速度 和浓度和浓度 可分别用平均值可分别用平均值 、 和偏差和偏差 、 之和之和来表示：来表示：uuCCuC显

11、然显然0u0CICCDuuCCtCCmgraddivdivuuuCCC于是连续性水动力方程可以写成：于是连续性水动力方程可以写成：注意到注意到：在典型单元体上的液相体积中取平均值，得在典型单元体上的液相体积中取平均值，得0uCuC0uCuC并且：并且：梯度的平均梯度的平均等于等于平均的梯度；散度的平均平均的梯度；散度的平均等于等于平均的散度；平均的散度；对时间导数的平均值对时间导数的平均值等于等于平均对时间求导平均对时间求导。可得：。可得：ICDuCuCtCmgraddivdivdivuCuCCuC )( uCuC ICCDuuCCtCCmgraddivdiv)(展开，得展开，得ICDCDuC

12、uCuCuCtCtCmmgraddivgraddivdivdivdivdivCDCDCDmmmgradgradgrad 在溶液连续体中的分子扩散系数在溶液连续体中的分子扩散系数 在多孔介质典型单元在多孔介质典型单元体的空隙体积上取平均之后变为体的空隙体积上取平均之后变为mD mD多孔介质的分子扩散系数多孔介质的分子扩散系数它是一个张量，一般讲，他在数值上要它是一个张量，一般讲，他在数值上要小于小于 。mDuC在平均过程中而引入的宏观水平上的附加变量为在平均过程中而引入的宏观水平上的附加变量为 机械弥散变量机械弥散变量，它是由于速度偏差，它是由于速度偏差 而产生的弥散而产生的弥散通量。如果不存在

13、速度偏差通量。如果不存在速度偏差 ，并忽略分子扩散，并忽略分子扩散，则溶质呈现则溶质呈现“活塞式推进活塞式推进”或迁移。或迁移。uuCDuCgrad 机械弥散系数（张量）机械弥散系数（张量） 实验表明实验表明：机械弥散通量类似介质中的分子扩散通量，也服从：机械弥散通量类似介质中的分子扩散通量，也服从于类似于类似FickFick扩散定律的形式：扩散定律的形式：DIuCCDtCdivgraddivDDD 水动力弥散系数水动力弥散系数上式称为水动力弥散方程。上式称为水动力弥散方程。CDJDgrad水动力弥散通量多孔介质的分子扩散系数机械弥散系数水动力弥散系数水动力弥散系数ICDDuCtC grad)

14、(divdiv记则建立建立水动力弥散方程，水动力弥散方程，我们涉及了我们涉及了3 3个水平，即分子水平、微观水个水平，即分子水平、微观水平和宏观水平。地下水动力学中，一般仅涉及宏观水平。但对于弥平和宏观水平。地下水动力学中，一般仅涉及宏观水平。但对于弥散问题，必须涉及到此散问题，必须涉及到此3 3个水平。这是因为：不讲分子水平，就无个水平。这是因为：不讲分子水平，就无法理解分子扩散、不讲微观水平，就搞不清机械弥散；法理解分子扩散、不讲微观水平，就搞不清机械弥散；但为了解决问题，我们最终不得不上升到宏观水平上来。但为了解决问题，我们最终不得不上升到宏观水平上来。 习惯上去掉“”，水动力弥散方程即

15、表示为：IuCxxCDtCijijdivIuCCDtCdivgraddiv写成微分的形式：3 , 2 , 1,ji I多组分组成的流体中，多组分组成的流体中，单位体积流体在单位时间内，由于化单位时间内，由于化学反应或其它原因所学反应或其它原因所产生产生（或（或消失消失）的某组分的）的某组分的质量质量。展开：展开：IuCzuCyuCxzCDyCDxCDzzCDyCDxCDyzCDyCDxCDxtCzyxzzzyzxyzyyyxxzxyxx对于一维流动二维水动力弥散：IuCxyCDxCDyyCDxCDxtCxyyyxxyxx对于一维流动一维水动力弥散：IuCxxCDxtCxxx2-5 源汇项源汇项

16、源汇项 系指在单位时间内、单位液相体积中由于化学反应、生系指在单位时间内、单位液相体积中由于化学反应、生物化学作用或抽注水等产生减少物化学作用或抽注水等产生减少组分的质量。组分的质量。2-5-1 放射性密度与化学、生物化学反应放射性密度与化学、生物化学反应设其变化规律为：设其变化规律为：teCC0CeCtCt0即：衰变速率与当时浓度成正比CuCxxCDtCijijdiv3 , 2 , 1,ji 若由于化学反应或生物化学反应而使示踪剂在单位体积溶液中的消耗速率或产生速率与其浓度成正比，也可以用上述式子表示。2-5-2 吸附与解吸 吸附与解吸：在一定条件下，溶液中某些溶质在多孔介质的固相表面产生吸

17、附、解吸等物理化学作用。 这些作用的结果应该综合到源汇项中：如果固相表面吸附示踪剂，称为吸附，视为汇；否则，称为解吸，视为源。(汇 )水吸 附解 吸离 子 交 换固 体(源 )(源 、 汇 )溶解相溶解相与吸附之间的吸附解析作用与吸附之间的吸附解析作用往往是一个可逆的过程：AAA吸附解析0)1()(tCntnC多孔介质多孔介质孔隙度孔隙度固体骨固体骨架体积架体积溶解相溶解相A A的密度的密度吸附相吸附相的密度的密度对于饱水多孔介质对于饱水多孔介质 n= =常数常数（孔隙率），则（孔隙率），则tCntCn)1 (tCnnI)1 (常数CnnC)1 ()(则则在体积为在体积为1 1的含水层中，如果

18、体系中某种溶质的吸附的含水层中，如果体系中某种溶质的吸附- -解吸为平衡，解吸为平衡，则有则有 I 表示单位时间内在单位液相体积由于这些作用增加或减少的示踪剂的质量。 （i）对于非均衡吸附作用：）对于非均衡吸附作用：0CKCKrf 吸附作用常数吸附作用常数解析作用常数解析作用常数rK （ii）对于均衡吸附作用：）对于均衡吸附作用：CKCKKCdrf平衡常数平衡常数tCKnntCnnId)1 ()1 (tCKnnId10CKCKrffKiiiijidCuxxCDxtCKnn 113 , 2 , 1,ji令ddKnnR11则上式可以写为：则上式可以写为：diijdijiRuCxxCRDxtC3 ,

19、 2 , 1,ji该方程形式上也不再含有源汇项。只是用只是用 去除以水动力弥散系数去除以水动力弥散系数 和流速和流速 ，由于，由于 ，因，因此吸附作用产生的后果，相对于此吸附作用产生的后果，相对于 和和 均减小均减小 ，起到减缓，起到减缓弥散的作用。所以把弥散的作用。所以把 称为：减缓因子。称为：减缓因子。DudR1dRDudR1dR将其代入对流将其代入对流弥散方程中整理得到：弥散方程中整理得到：2-5-3 抽水与注水抽水与注水 如果有抽水或注水井，含水层中示踪剂的质量就会发生变化：如果有抽水或注水井，含水层中示踪剂的质量就会发生变化：（i）当抽水时：）当抽水时：若假设单位时间内从单位体积含水

20、层中的抽水量为若假设单位时间内从单位体积含水层中的抽水量为W。则。则CnWI孔隙率孔隙率为抽水点处的溶质浓度为抽水点处的溶质浓度表示失去的溶质表示失去的溶质则，水动力弥散方程可写为：则，水动力弥散方程可写为：CnWCuxxCDxtCiijiji3 , 2 , 1,ji（ii）单位时间内向单位体积含水层中注入含有示踪剂的水（示踪）单位时间内向单位体积含水层中注入含有示踪剂的水（示踪剂浓度剂浓度C0）)(0CCnWI水动力弥散方程为：水动力弥散方程为：)(0CCnWCuxxCDxtCiijiji3 , 2 , 1,ji0CnWI0CnWCuxxCDxtCiijiji或或 若含水量若含水量不等于常数

21、（饱气带）。不等于常数（饱气带）。（iii）含水层中注入浓度）含水层中注入浓度C0的放射性示踪剂，该示踪剂又与固体颗的放射性示踪剂，该示踪剂又与固体颗粒发生均衡吸附作用，此时的水动力弥散方程可表示为：粒发生均衡吸附作用，此时的水动力弥散方程可表示为：0CnWCRuCxxCRDxtCdiiidiji3,2, 1,ji则：则：0WCCRuCxxCRDxtCdiiidiji3,2,1,ji研究非饱和带中的溶质运移研究非饱和带中的溶质运移非饱和水动力弥散方程。2-6 初始条件与边界条件水动力弥散方程揭示了溶质在地下水中运移的一般规律，对于一个具体问题，我们必须知道其初始的状态，以及边界条件，才能达到地

22、下水中溶质的空间分布规律及其随时间的变化。2-6-1 初始条件初始条件 描述综合初始时刻，研究区D内各点（x,y,z)处的浓度分布状态的条件（数学表达式） 初始条件如),(|00zyxCCtDzyx,初始条件确指原始状态；初始时刻可以任意选定，只要已知那一时刻研究区各点的浓度即可。 例如：t=0时向某区域注入含示踪剂的水，若在此之前研究区D不含该示踪剂，则C(x,y,z)=0。如：在弥散试验时，可将示踪剂注入前浓度分布视为初始状态。设计地下污水治理方案时，可将现状污染物分布视为初始条件。2-6-2 边界条件 边界条件指的是研究区边界上的溶质浓度分布和变化情况或边界上流入（或流出）研究区的浓度分

23、布和变化情况。主要有以下三类情况 1、第一类边界条件：给定浓度边界，即已知边界上的浓度分布。 tzyxfCB,110,1tBzyx为B B1上的已知函数为研究区D的第一类边界 2、第二类边界条件 给定弥散通量边界 指已知边界弥散通量随时间变化规律的边界条件，或者称之为Neumann边界条件，tzyxfnCDB,grad220,2tBtzyx水动力弥散系数为研究区上的第二类边界为边界B2上某点（x，y，z）处的外法线方向上的单位向量已知函数，定义在B2上。 3、第三类边界条件 给定溶质通量边界 指已知边界上溶质通量随时间变化规律的边界条件，或称之为Cauchy边界条件。zyxfnCDuCB,gr

24、ad330,3tBtzyx孔隙平均流速已知函数为研究区上的第三类边界数学模型：描述实际问题的函数或数学方程。数学模型：描述实际问题的函数或数学方程。地下水溶质运移数学模型：微分方程地下水溶质运移数学模型：微分方程+初始条件初始条件+边界条件。边界条件。一般可以写成：一般可以写成：2-7 数学模型IuCxxCDxtCxxx),(|00zyxCCtDzyx,zyxfnCDuCB,grad330,3tBzyx0,tDzyx例例1 1： 无限长多孔介质柱体，初始示踪剂呈阶梯函数分布无限长多孔介质柱体，初始示踪剂呈阶梯函数分布 设有一无限长均质砂柱，以速度设有一无限长均质砂柱，以速度u做稳定流动，且初始

25、浓度呈阶做稳定流动，且初始浓度呈阶梯状分布（梯状分布（图图48）。）。解：该问题属于一维稳定流动一维水动力弥散问题，取坐标系如图解：该问题属于一维稳定流动一维水动力弥散问题，取坐标系如图4-8，则数学模型可写成，则数学模型可写成xCuxCDtCL2200)0 ,(10 xCxCxC0),(0tCtC0),(1tCtC例例2 2： 半无限长多孔介质柱体，一端为定浓度边界半无限长多孔介质柱体，一端为定浓度边界 设有一半无限长均质砂柱，一维稳定流动，孔隙平均流速为设有一半无限长均质砂柱，一维稳定流动，孔隙平均流速为u，其一端为定浓度边界，求其浓度分布。其一端为定浓度边界，求其浓度分布。 解：取坐标系

26、如图解：取坐标系如图410所示。该问题的数学模型可描述如下所示。该问题的数学模型可描述如下 xCuxCDtCL220)0 ,(xC0), 0(CtC0),( tC例3： 一维稳定流动二维水动力弥散问题 在均质各向同性、等厚的承压含水层中存在着一维稳定流在均质各向同性、等厚的承压含水层中存在着一维稳定流动，孔隙平均流速为动，孔隙平均流速为u，取，取x坐标轴平行地下水流向（图坐标轴平行地下水流向（图411）。）。这时，对于水动力弥散系数而言却是均质各向异性的，在这时，对于水动力弥散系数而言却是均质各向异性的，在x方向方向为纵向弥散系数为纵向弥散系数DL，而在，而在y方向为横向弥散系数方向为横向弥散

27、系数DT。 xCuyCDxCDtCxTL22220)0 ,(yxC0), 0(CtyC0),(tyC0),( txC 习题： 1.设计Fick试验仪，验证Fick定律，计算多孔介质分子扩散系数及水动力弥散系数。2.写出一维砂柱溶质运移数学模型 例一、如下图多孔介质A的边界外为另一种介质B则在A-B的边界上，溶质通量应该保持BLxALxxCDCunxCDCun21n为孔隙率为孔隙率例二、如右图多孔介质为无孔隙介质，则通过该边界的流量和溶质的质量都为0。000BxLxxCuxCDu0LD例三、多孔介质外围为Cr的河水。0BdrLxdrLxxCDuCxCDCunxCDuCxCDCun若忽略分子扩散：若忽略分子扩散：0BrLxuCxCDCu当是湖泊，水流静止：当是湖泊，水流静止：0BLxxCDCu例四、多求介质边界为空气，此时两侧的浓度相等：0BxC浓度不变，变化率为0。方向导数：方向导数：coscoscoszuyuxulu梯度：梯度：sMusMulu)()(lim定义：若在数量场中的一点定义：若在数量场中的一点M处，存在这样的矢量处，存在这样的矢量G，其方，其方向为函数向为函数u(M)在点在点M处变化率最大的方向，其模也正好是处变化率最大的方向，其模也正好是这个最大变化率的数位