饱和-非饱和土壤中吸附性溶质水动力弥散实验研究(全文).docVIP

饱和-非饱和土壤中吸附性溶质水动力弥散实验研究 张富仓，康绍忠，潘英华(西北农林科技大学 教育部旱区农业水土工程重点实验室) 摘 要：本文通过室内土柱对轻壤土和重壤土饱和、非饱和水动力弥散实验，得到了吸附性溶质水动力弥散参数的计算方法。根据一维饱和弥散实验结果，由相应的解析解，用简化计算方法求得饱和水动力弥散系数Dsh和阻滞因子R；对于一维非饱和弥散实验，在分析溶质的水动力弥散理论的基础上，推导出求解非饱和土壤水动力弥散系数Dsh(θ)的计算公式。 关键词：吸附性溶质；水动力弥散系数；阻滞因子 基金项目：国家杰出青年基金研究项目的部分内容作者简介：张富仓(1962-)，男，陕西人，副研究员，主要研究方向：土壤水动力学和溶质迁移。 溶质在多孔介质中运动是一个复杂的物理化学过程，自70-80年代以来，国外许多学者对此进行了较多的研究[1～3]。近年来，国内也开展了这方面的工作，取得了一些进展。然而大多数的研究都是对非吸附性溶质如Cl-、Br-等离子进行的[4～5]，而对吸附性溶质的水动力弥散试验研究，大多数是在饱和条件下进行的，其溶质的水动力弥散系数的确定方法依据实测浓度的最小二乘法、最大拟然法或者用近似的解析公式求得[6～8]。对吸附性溶质在饱和和非饱和土壤水动力弥散系数及其相关参数的测定和计算国内至今没有统一方法，国外基本形成了以Laryea(1982)为代表的测定和计算方法。本文在国内外学者研究的基础上，以陕西省2种典型土壤为例，对吸附性溶质在饱和与非饱和土壤的水动力弥散系数及其有关参数的室内测定和计算方法进行了试验研究。 1 基本原理 1.1 饱和水动力弥散系数的测定原理 一维饱和-非饱和土柱中，水流为稳态、溶质为吸附性的CDE，即： (1) 式中：C为土柱溶质浓度(g/L)；S为单位质量土壤中吸附溶质的量；γd为土壤容重；ν为孔隙平均水流速(cm/min)；x为距离变量(cm)；t为时间变量(min)。 对于大多数吸附性离子如K+和重金属离子，溶质的吸附符合线性吸附，即S=KdC，则上式可改写为： (2) 令 R=(1+γd/θKd) (3) (4) 式中：R为阻滞因子(量纲为1)；式(4)表达了在稳态水流条件下反应性溶质的运移方程，溶质在土壤中吸附、交换阻滞、散布等都可用同一参数阻滞因子R综合表示。Dsh为饱和水动力弥散系数或称离弥散系数。 对于一维饱和流土壤水动力弥散系数Dsh，可用垂直土柱法测定，对于式(4)，其定解条件可以描述为 C(x,0)=0,0≤x＜∞ (5a) C(0,t)=C0,t＞0 (5b) C(∞,t)=0,x→∞，t＞0 (5c) C/C0=1/2{erfc[x-(νt)/R/2+exp(xv/Dsh)·erfc[x+(νt)/R/2]]} (6) 当D/(νx)0.075时，上式可近似表示为： C/C0=1/2erfc[Rx-νt/2] (7) 式中：erfc为余误差函数，Abramowitz和Stegun(1970)定义为 (8) 根据实测资料，利用式(7)可求得Dsh和R.因一维饱和流入渗条件下，土壤孔隙平均水流速ν与时间t关系表达式为： t=TL/ν (9) 式中：T为孔隙体积(V/V0)；V为出流体积；V0为土壤空隙体积；L为土柱长度。 将式(9)代入式(7)可得： (RLν/4DshT)0.5-(LνT/4DshR)0.5=erfc-1(2C/C0) (10) 式中：erfc-1(2C/C0)为2C/C0的补余误差函数的逆。其中C/C0和T由实验测定可得，设n为实验测定点数，用最小二乘法可得到R和Dsh关系为： R=a/b (11) Dsh=Lν/(4ab) (12) 其中： a=n∑[Ti0.5erfc-1(2Ci/C0)]－∑T\-i∑[erfc-1(2Ci/C0)/T i0.5]/n2－∑Ti∑(１/Ti) (13) b=(a/n)∑(1/Ti)-(1/n)∑[erfc-1(2Ci/C0)/Ti0.5] (14) 　　图1表示了相对浓度Ci/C0从0.05～0.95变化时erfc-1(2Ci/C0)的关系曲线。根据实验得到的Ci/C0、孔隙体积Ti以及实验点数n，由式(13、14)可求得a和b值，进而得到R和Dsh的值。 1.2 非饱和水动力弥散系数的测定原理 若溶质是吸附性溶质，则单位土壤体积某种溶质的数量(μg/cm3)用下式表示， 图1 erfc-1(2C/C0)与C/C0关系曲线 M=θC+γdS (15) 　　则非饱和土壤溶质运移方程为： (16) (17) 以KNO3溶液为例，在水平非饱和入渗条件下，式(17)的初始和边界条件为： 初始：θ=θi,C-=C-i，C+=C