溶质运移理论(一)水动力弥散的基本概念与弥散方程精选课件.ppt

1、地下水流溶质运移理论及模型地下水流溶质运移理论及模型第一章 水动力弥散的基本概念第二章 水动力弥散方程2一、流体类型一、流体类型3一、流体类型一、流体类型可混溶流体可混溶流体石油污染物在水体或含水层中的运移不可混溶流体不可混溶流体不同性质溶体之间无明显的突变界不同性质溶体之间无明显的突变界不同性质溶体之间有明显的突变界不同性质溶体之间有明显的突变界油油水污污染染物物水4二、示踪剂二、示踪剂5二、水动力弥散二、水动力弥散6水动力弥散分子扩散机械弥散由浓度高的方向向浓度底 的方向运动，趋于均一由于微观多孔介质中流速分布的不均一而引起的示踪剂（水质点）浓度在地下水含水层中不均匀分布的现象。三、三、水

2、动力弥散现象水动力弥散现象两部分两部分7四、机械弥散四、机械弥散（1）（2）（3）地下水质点运动速度的差异是产生水动力弥散的根本原因地下水质点运动速度的差异是产生水动力弥散的根本原因8四、机械弥散四、机械弥散平行于平行于平均流平均流速方向速方向上的弥上的弥散散垂直于垂直于平均流平均流速方向速方向上的弥上的弥散散纵向弥散纵向弥散横向弥散横向弥散Fick定律定律o 流体通量流体通量J J9五、流体参数五、流体参数dVdVmdVmNNNdmdd111Jo 多相分流体速度多相分流体速度 组分流速组分流速NNN111J10控制方程：控制方程：质量守恒定理质量守恒定理 在多组分组成的流体体系中任取一点P(

3、x,y,z),以P为中心取一微小的质量平衡体，其侧面分别平行与3个坐标面，边长分别为x、y、z。质量守恒原理:在时间t内，组分在这个单元体中的净流出（或流出）量（暂不考虑起内部有质量产生和消失），应等于这个单元中组分的质量变化11zyxtyxuutzxuutzyuuzzyxzzzyxzzyyxyzyyxyzyxxxzyxxx2,2,2,2,2,2zyxuuu,设设 分别表示分别表示组分密度、组分密度、x,y,z方向的速度方向的速度分量。分量。控制方程：控制方程：质量守恒定理（续）质量守恒定理（续）12其中：其中：经过经过t时间后，质量均衡体中时间后，质量均衡体中 的变化量。的变化量。将上式左右

4、两端同除以将上式左右两端同除以 得：得：tzyxtzuuyuuxuuzzyxzzzyxzzyyxyzyyxyzyxxxzyxxx/2,2,2,2,2,2控制方程：控制方程：质量守恒定理（续）质量守恒定理（续）13再对方程两端取极限，即令再对方程两端取极限，即令0,0,0,0tzyxtzuyuxuzyx即有即有：即即0divut若微小的质量均衡体内存在着若微小的质量均衡体内存在着组分的源汇项，则上式可组分的源汇项，则上式可改写为：改写为：Iutdiv多组分流体体系中多组分流体体系中组分的质量守恒方程组分的质量守恒方程控制方程：控制方程：质量守恒定理（续）质量守恒定理（续）14 I 是多组分组成的

5、流体中，是多组分组成的流体中，单位体积流体在单位体积流体在单位时间内，由于化单位时间内，由于化学反应或其它原因所产生（或消失）的学反应或其它原因所产生（或消失）的组分的质量。组分的质量。IutdivJutuutuututIdivdivdivdivdivdivuuu控制方程：控制方程：质量守恒定理（续）质量守恒定理（续）15ICDuCtCmgraddivdiv稀释的二元体系中稀释的二元体系中组分的对流组分的对流扩散方程扩散方程应用条件：1、二元体系二元体系；2、等温条件等温条件；3、低浓度低浓度；IutdivCDJmgrad控制方程：控制方程：质量守恒定理质量守恒定理+Fick定律定律16对于一

6、维流动二维水动力弥散：对于一维流动二维水动力弥散：IuCxyCDxCDyyCDxCDxtCxyyyxxyxx对于一维流动一维水动力弥散：对于一维流动一维水动力弥散：IuCxxCDxtCxxx对流对流扩散方程扩散方程简化成简化成17六、水动力弥散方程六、水动力弥散方程多孔介质中溶质的分子扩散通量多孔介质中溶质的分子扩散通量多孔介质分子扩散系数，数值上小于多孔介质分子扩散系数，数值上小于溶质的对流量溶质的对流量机械弥散通量机械弥散通量联立上述两式，得联立上述两式，得（1）（2）将所有平均号将所有平均号“-”略去略去18六、水动力弥散方程六、水动力弥散方程19六、水动力弥散方程六、水动力弥散方程20

7、七、源汇项七、源汇项1）放射性密度与化学、放射性密度与化学、生物生物化学反应化学反应2）吸附与解吸吸附与解吸3）抽水与注水抽水与注水21七、源汇项：放射性生物化学反应七、源汇项：放射性生物化学反应teCC0ceCtCt0cuCxxCDdivtCijij 若研究对象是地下水中某种放射性物质作为示踪剂，若研究对象是地下水中某种放射性物质作为示踪剂，则它的浓度分布受对流和弥散的影响外，还将受到其自身则它的浓度分布受对流和弥散的影响外，还将受到其自身的放射性衰变的影响。的放射性衰变的影响。若由于化学反应或生物化学反应而使示踪剂在单位体积溶液中的消耗速率若由于化学反应或生物化学反应而使示踪剂在单位体积溶

8、液中的消耗速率或产生速率与其浓度成正比，也可以用上述式子表示。或产生速率与其浓度成正比，也可以用上述式子表示。22七、源汇项：吸附与解吸七、源汇项：吸附与解吸(汇)水吸附解吸离子交换固体(源)(源、汇)在一定条件下，溶液中某些溶质在多孔介质的固相表在一定条件下，溶液中某些溶质在多孔介质的固相表面产生吸附、解吸或者离子交换等物理化学作用。如果这面产生吸附、解吸或者离子交换等物理化学作用。如果这些溶质属于我们的研究对象，则这些作用的结果应该综合些溶质属于我们的研究对象，则这些作用的结果应该综合到源汇项中，如果固相表面吸附示踪剂，视为汇，否则，到源汇项中，如果固相表面吸附示踪剂，视为汇，否则，称为解

9、吸，视为源，而离子交换即可视为汇也可视为源。称为解吸，视为源，而离子交换即可视为汇也可视为源。23七、源汇项：吸附与解吸七、源汇项：吸附与解吸溶解相与吸附之间的吸附解析作用溶解相与吸附之间的吸附解析作用,往往是一个可往往是一个可逆的过程逆的过程:AAA吸附解析tCtCIb汇多孔介质含水率多孔介质含水率固体骨架密度固体骨架密度溶解项溶解项A的密度的密度吸附相的密度吸附相的密度24七、源汇项：吸附与解吸七、源汇项：吸附与解吸 （i）对于非均衡吸附作用：CKCKIrf吸附作用常数AA解析作用常数AA非均衡吸附作用非均衡吸附作用均衡吸附作用均衡吸附作用吸附作用吸附作用25对于均衡吸附作用：对于均衡吸附

10、作用：CKCKKCddf平衡常数对于饱水多孔介质对于饱水多孔介质=n=C（孔隙率）（孔隙率）tCKtCtCdbbtCKntCdb七、源汇项：吸附与解吸七、源汇项：吸附与解吸26iiiijidbCuxxCDxtCKn 1dbdKnR1diijdijiRuCxxCRDxtC只是用只是用 去除以水动力弥散系数去除以水动力弥散系数 和流速和流速u，由于，由于 ，因，因此吸附作用产生的后果，相对于此吸附作用产生的后果，相对于 和和 均减小均减小 ，起到减缓，起到减缓弥散的作用。所以把弥散的作用。所以把 称为：减缓因子。称为：减缓因子。DudR1dRDdR1dR七、源汇项：吸附与解吸七、源汇项：吸附与解吸

11、27七、源汇项：抽水与注水七、源汇项：抽水与注水 如果含水层当中有抽水或注水井，含水层中示踪剂如果含水层当中有抽水或注水井，含水层中示踪剂的质量就会发生变化，而且抽水与注水导致源汇项的变的质量就会发生变化，而且抽水与注水导致源汇项的变化不同。当抽水时：化不同。当抽水时：假设单位时间内从单位体积含水层中的抽水量为假设单位时间内从单位体积含水层中的抽水量为W。CnWI孔隙率孔隙率为抽水点处的溶质浓度为抽水点处的溶质浓度表示失去的溶质表示失去的溶质28七、源汇项：抽水与注水七、源汇项：抽水与注水当注水时当注水时若向含水层中注入含有示踪剂的水（示踪剂浓度若向含水层中注入含有示踪剂的水（示踪剂浓度C0）

12、0CnWI0CnWCuxxCDxtCiijiji29八、初始条件与边界条件八、初始条件与边界条件 水动力弥散方程揭示了溶质在地下水中运移的一般水动力弥散方程揭示了溶质在地下水中运移的一般规律，对于一个具体问题，我们必须知道其初始的状态，规律，对于一个具体问题，我们必须知道其初始的状态，以及边界条件，才能达到地下水中溶质的空间分布规律以及边界条件，才能达到地下水中溶质的空间分布规律及其随时间的变化。及其随时间的变化。初始条件初始条件边界条件边界条件30八、初始条件与边界条件八、初始条件与边界条件 描述综合初始时刻，研究区描述综合初始时刻，研究区D内各点（内各点（x,y,z)处的浓度处的浓度分布状

13、态的条件（数学表达式）分布状态的条件（数学表达式）zyxzyxCC,00,初始条件确指原始状态；初始时刻可以任意选定，只要已知那一时初始条件确指原始状态；初始时刻可以任意选定，只要已知那一时刻研究区各点的浓度即可。初始条件的如何选取，应该根据研究问刻研究区各点的浓度即可。初始条件的如何选取，应该根据研究问题的需要、资料状况及计算与模拟方法等因素确定。例如：题的需要、资料状况及计算与模拟方法等因素确定。例如：t=0t=0时向时向某区域注入含示踪剂的水，若在此之前研究区某区域注入含示踪剂的水，若在此之前研究区D D不含该示踪剂，则不含该示踪剂，则C C（x x，y y，z z）=0=0。如：在弥散

14、试验时，可将示踪剂注入前浓度分布视。如：在弥散试验时，可将示踪剂注入前浓度分布视为初始状态。又如设计治理地下污水方案时，可将现状污染物分布为初始状态。又如设计治理地下污水方案时，可将现状污染物分布视为初始条件。视为初始条件。31 边界条件指的是研究区边界上的溶质浓度分布和变化边界条件指的是研究区边界上的溶质浓度分布和变化情况或边界上流入（或流出）研究区的浓度分布和变化情情况或边界上流入（或流出）研究区的浓度分布和变化情况。况。主要有以下三类情况：第一类、第二类与第三类。主要有以下三类情况：第一类、第二类与第三类。第一类边界条件第一类边界条件：也叫给定浓度边界，即已知边界：也叫给定浓度边界，即已

15、知边界上浓分布的边界。上浓分布的边界。tzyxfCB,110,1tBzyx为B1上的已知函数为研究区D的第一类边界八、边界条件八、边界条件32八、边界条件八、边界条件第二类边界条件：第二类边界条件：给定弥散通量边界给定弥散通量边界,指已知边界弥散通指已知边界弥散通量随时间变化规律的边界条件，或者称之为量随时间变化规律的边界条件，或者称之为Neumann边界条件，边界条件，tzyxfngradCDB,220,2tBtzyx为边界为边界B2上某点（上某点（x，y，z）处的外法线方）处的外法线方向上的单位向量向上的单位向量水动力弥散系数水动力弥散系数研究区上的第二类边界研究区上的第二类边界已知函数，已知函数，定义在定义在B2上上33八、边界条件八、边界条件第三类边界条件第三类边界条件:指已知边界上溶质通量随时间变化指已知边界上溶质通量随时间变化规律的边界条件，或称之为规律的边界条件，或称之为Cauchy边界条件。边界条件。zyxfnDgradCuCB,33已知函数已知函数孔隙平均流速孔隙平均流速