非饱和土中溶质地表径流迁移模型及解析模拟 Model for chemicals transfer to surface runoff in unsaturated soil and its analytical modeling.pdfVIP

第20卷第1期 水科学进展 、瑚．20．No．1 2009年1月 ADVANCESINWATERScIENCE Jan．。2009 非饱和土中溶质地表径流迁移模型及解析模拟 童菊秀，杨金忠，暴入超 (武汉大学水资源与水电工程科学国家重点实验室，湖北武汉430072) 摘要：针对非饱和农田在降雨引起肥料溶质地表径流迁移流失的实际问题，建立了二层结构溶质地表径流迁移模 型，根据溶质质量守恒定律和水量平衡原理，求得模型的解析解，并采用室内非饱和土壤实验验证解析模型的合理 性。通过敏感性分析．探讨了模型参数和土壤物理参数对溶质地表径流迁移流失过程的影响。研究结果表明：非饱 和土的吸附性增大将增强土壤溶质迁移流失到地表径流水的能力；土壤干容重增大将减少土壤溶质的地表径流迁 移；土壤的初始体积含水率越高，土壤溶质的地表径流运移作用越强。 关键词：非饱和土壤；溶质迁移；地表径流；解析模型 中图分类号：P641．2 文献标识码：A 文章编号：1001-6791(2009)0l瑚l咖8 在农田灌溉或降雨过程中，土壤肥料地表径流迁移流失的研究多数假定初始土壤处于饱和状态，在灌溉或 降雨刚开始就产生径流【卜4|，且土壤侵蚀比较严重。然而在一般情况，刚开始灌溉或降雨时，地下水位往往在 地表以下，此时地表以下一定深度内的土层处于非饱和状态。农田地面比较平坦，土壤侵蚀不甚严重，且农田 中有田埂阻挡，田面形成一定的积水才产生径流。研究表明，在降雨过程中，土壤在雨滴打击及径流冲刷作用 下，形成一定厚度的扰动层，称之为“混合层”【卜6I。混合层内溶质参与径流迁移，而此层以下的溶质不参与 大，而且随着深度增加，进入径流几率呈指数递减。因此他们提出了有效混合深度(后也称为土壤混合层深度) J发现，完全、不完全混合理论模型不能很好 的概念，即在这个深度内溶质均匀分布且完全混合。王全九等【7,8 反映非饱和土壤溶质径流迁移的过程，进一步发展了非完全混合理论，将径流与入渗作用分别考虑，提出等效 径流迁移深度的概念，并应用到非饱和黄土高原土壤中。 本文运用非完全混合理论，对非饱和土壤进行模拟降雨，考虑到积水层的存在，并将土壤混合层和积水． 径流混合层看作一个整体，整个抽象为混合区．混合区以下的二层结构体系。将混合区作为主要研究对象，以 溶质质量守恒定律和水量平衡方程为基础，建立了溶质地表径流迁移控制方程，考虑到非饱和土壤中入渗率的 变化特征，得到问题的解析解。通过对已有实验数据的对比分析，证明了所提出模型的正确性和有效性。最后 通过分析模型和实验中的参数，探讨了非饱和土壤的物理特性对溶质地表径流迁移流失的影响。 1二层结构数学模型． 1．1溶质质量守恒方程 溶质以吸附性和溶解性存在于土壤中，在地表附近的土壤混合层内的溶解性溶质浓度为Q(图1)，采用下 列的平衡线性吸附方程考虑土壤混合层中土颗粒对溶质的吸附特性[6．9]： cl=I|}dG (1) 收稿日期：2008-01．16 100205； 2006AA啪207) 作者简介：童菊秀(1981一)，女，湖北武汉人，博士研究生，主要从事地下水、土壤水运动理论及水资源研究。 E-mail：tongjxlll@163．cofll 万方数据 第1期 童菊秀，等：非饱和土中溶质地表径流迁移模型及解析模拟 式中 kd为溶质在固．液相间的吸附分布系数；G为土壤混合层中吸附性溶质浓度；G为土壤混合层中溶解 性溶质浓度。 土壤混合层内的溶质是入渗水、径流水溶质的来源，假定无侧向壤 中流，仅考虑土壤中垂向的溶质运移[6]。忽略径流水中微量土壤颗粒， 即不考虑径流水中吸附性溶质。为了简化，在整个混合区及其以下土壤 的交界处，人渗和扩散运动可视为以垂直向下