swd5土壤溶质运移

土壤水动力学(SoilWaterDynamics)

第5章土壤溶质运移

(SoilSoluteTransport)尚松浩清华大学水利水电工程系水文水资源研究所Email:shangsh@第5章土壤溶质运移土壤溶质运移及其基本方程土壤溶质运移的其它分析方法土壤溶质运移的数值模拟参考文献：李韵珠，李保国.土壤溶质运移.北京：科学出版社,19985.1土壤溶质运移及其基本方程土壤溶质土壤溶质运移的对流与弥散溶质运移的对流－弥散方程土壤中溶质运移与水分运动的关系5.1.1土壤溶质土壤液相是包含各种溶质的水溶液有机溶质：可溶性有机物：氨基酸、腐殖酸、糖类、有机－金属离子络合物有机污染物：酚类、苯胺类、硝基苯类、有机氯农药、有机磷农药、阿特拉津、丙烯腈和丙烯醛…无机溶质：八大离子：Ca2+、Mg2+、Na+、K+、Cl-、SO42-、CO32-、HCO3-其它微量离子：Fe3+、Cu2+…无机污染物：铅、镐、汞等重金属元素；氟、砷和硒等非金属元素；无机氮溶质形态：离子态、水合态、络合态土壤溶液中溶质的来源：自然来源：降水带入的可溶盐：滨海地区，火山喷发，工业废气土壤基质中盐分溶解土壤生物活动：游离的氨基酸、糖类等可溶性有机物地下水中的盐分随潜水蒸发的迁移人工来源：灌溉、排污、农药及化肥的施用…土壤溶液的性质：浓度与活度：浓度:质量浓度(kg/m3)、摩尔浓度(mol/m3)…离子活度：实际溶液中离子的有效浓度或热力学浓度离子强度：导电性酸碱性氧化还原性土壤溶质的物理、化学、生物变化过程：溶解与沉降吸附与交换水解与络合氧化还原生物化学过程：有机质分解，氮素转化…土壤溶质的运移对流——随土壤水分运动的运移扩散——土壤溶液浓度梯度所引起的溶质运移弥散——土壤孔隙中水的微观流速的不均匀性所产生的运移土壤溶质运移研究的意义土壤盐碱化发生规律与防治土壤污染——粮食、生态地下水污染5.1.2土壤溶质运移机理

——对流、分子扩散与机械弥散对流(Convection)：溶质随土壤水分的运动而运移溶质对流通量：Jc=qc=vθcq：土壤水通量c：土壤水溶液的质量浓度v=q/θ：土壤水溶液的平均孔隙流速对流作用下溶质的穿透时间：tL=L/v=Lθ/q分子扩散(Moleculardiffusion)：不同浓度的溶质间在浓度梯度的作用下产生的质量交换原因：分子不规则热运动（布朗运动）趋势：高浓度→低浓度溶质的分子扩散符合Fick第一定律：自由水溶液：土壤水溶液：Jd——分子扩散通量；Ds——扩散系数溶质扩散系数Ds

：随土壤含水率的减小Ds(θ)=D0θαDs(θ)=D0aexp(bθ)b=10,a=0.001~0.005机械弥散(mechanicaldispersion)：又称对流弥散(convectiondispersion)，由于土壤水流和土壤孔隙结构相互作用使溶质在土壤孔隙中运移和分散的过程产生原因：土壤孔隙不均匀，流速方向和大小不同，使得溶质在水流运动中分散溶质的机械弥散通量：机械弥散系数Dh(v)：Dh(v)=λ|v|经验常数λ与土壤质地、结构有关机械弥散系数还与空间尺度有关水动力弥散(hydrodynamicdispersion)：综合考虑分子扩散与机械弥散溶质的水动力弥散通量：水动力弥散系数(综合扩散－弥散系数)：Dsh(v,θ)当土壤水流速度较大时，机械弥散作用超过扩散作用；土壤溶液静止时，只有分子扩散作用溶质的混合置换：置换流体(Displacingfluid)被置换流体(Displacedfluid)穿透曲线(BTC:Breakthroughcurve)5.1.3溶质运移的对流－弥散方程土壤溶质通量：溶质运移连续方程(一维，不考虑源汇项)一维溶质运移方程：

根据水流连续方程：对流－弥散方程：CDE稳定流情况下：土壤含水率θ

、水流通量q、平均孔隙流速v均为常数，Dsh(v,θ)为常数，CDE简化为线性PDE，在一些简单情况下可以求出解析解：考虑源汇项：Sc——单位时间、单位体积土壤中由于化学、生物作用所生成（消失减少）的溶质质量考虑溶质的动态储存(沉淀与溶解)：流动区、非流动区双区模型：θm、θim—土壤水流动区和非流动区的含水率cm、cim—土壤水流动区和非流动区的溶质浓度vm—土壤水流动区的平均孔隙流速流动区和非流动区之间溶质迁移：β——质量迁移系数5.1.4土壤溶质运移与水分运动的关系土壤水分运动对溶质运移的影响水流通量q影响溶质对流土壤含水率θ

、水流通量q（土壤平均孔隙流速v=q/θ）影响溶质的水动力弥散土壤溶质的源汇项、动态储存一般与θ有关溶质对土壤水分运动的影响溶质势：ψs=-cRT/(ρwgM)

水流通量：

f0——选择系数：＝1：完全选择，溶质不能通过(半透膜)；＝0：溶质能通过，溶质势对水流无影响粘土层存在一定的半透膜作用，非粘土中一般可以忽略溶质势的影响溶质的存在导致土壤导水率的降低土壤溶质运移模型的其它问题：二维、三维：多组分溶质反应运移参数：尺度效应纵向弥散系数、横向弥散系数土壤水盐耦合运移：水、溶质方程联合求解忽略溶质对水分运动的影响：先求解水流方程，然后求解溶质运移方程5.2土壤溶质运移的其它分析