土壤水盐运移的解析解及其水动力弥散系数D的确定.pdfVIP

农业高效用水理论与技术 ， 土壤水盐运移的解析解及 水动力弥散系数D的确定 肖娟1 雷廷武2江培福3 (1．太原理工大学水利科学与工程学院，太原030024；2．中国农业大学水利与土木工 程学院，北京 l00083；3．水利部综合事业局，北京100053) 摘要：通过室内小土柱试验，我们得到饱和土壤孔隙体积数(时间)和对应的流出液浓 度的数据，即可绘出一条穿透曲线。采用一维对流一扩散方程，在一定的边界条件下求出了 该方程的解析解。．把解析解与实验所得的穿透曲线加以比较，用黄金分割法搜索，最后取得 最佳拟合时的水动力弥散系数D。结果表明，经典的一维对流一扩散方程能较成功地预报本 实验结果。但当水力传导度很小，盐分不能完全穿透土柱时，拟合效果较差，另外，对相同的 钠吸附比SAR，溶液浓度越高，则拟合效果越好；对相同的盐分浓度，SAR越小，拟合效果越 好。水动力弥散系数随着孔隙流速的增加而增加，对于粉质壤土，指数函数能最好地拟合计 算得到的数据，而对于粉质粘壤土，直线关系能最好地拟合数据。 关键词：解析解；水动力弥散系数；穿透曲线；孔隙平均流速 引 言 土壤水力传导度是随着土壤可交换钠的百分比(ESP)的增加和土壤溶液电导率(Ec)的 在盐碱土壤上用咸水灌溉时，为保证土壤可持续运行及考虑地下水不受污染，盐分在土壤中 的运移的研究受到广泛的重视。而穿透曲线或者说混合置换实验中的流出液浓度随时间变 化的测定则是研究溶质在土壤中运移的常用方法。 土壤中溶质的运动是十分复杂的。溶质运移主要是通过对流和水动力弥散两种作用实 现的(雷志栋，1988)，两者之间相互联系。有许多数学模型可以描述土壤溶质的运移(E． 溶质运移模型分析水盐运动规律时，水动力弥散系数D是在数值计算中需要确定的一个重 and Bennett，1995)。 要参数(zheng 国内外测定土壤中水动力弥散系数的计算方法有很多。用穿透曲线法求解水动力弥散 ·691· 现代节水高效农业与生态灌区建设 uwe 作都只局限于流速恒定的情况(Bresler，1982)。Kai 据同时估算孔隙平均流速伽和水动力弥散系数风的斜率法及等斜率法，对金属在土壤运移 进行估计，结果表明，该方法具有计算快速而简单，对饱和及非饱和流均适用的优点。任理等 (1999)在上述基础上通过BTC的实验数据，采用“新陈代谢”思想在抛物插值条件下，应用斜 率法求得水动力弥散系数，在实际应用中取得良好的结果。朱红霞等(1998)采用一维对流一 扩散方程研究了非反应性化学物质在土壤中的迁移行为，并确定土壤不同团聚体粒径和不 同孔隙流速时的水动力弥散系数，参数的确定是采用方程的解析解，通过最小二乘法，最佳 拟和穿透曲线而得到。其试验值和模拟值有很好的相关性。在众多的研究中，无论是同时确 定№和风，还是在已知蜘时估算风，都是采用的孔隙平均流速，如做出盐分通过所研究 土体不同时间的概率，其中以50％的概率所对应的时间作为盐分平均迁移时间(又称中值 时间)，再计算迁移的平均速度(任理，1999b)。并且不同水质时的水动力弥散系数的研究也 较少。 本研究主要以两种土壤为研究对象。 研究饱和土壤不同水质条件下的水盐运 移情况和水动力弥散系数的确定。参数D 叫¨ 可通过用对流一弥散方程最佳拟合穿透 曲线而得到，在拟合过程中不采用平均流 速，而是以每一实测时间的实际流速为基 础，消除了以孔隙平均流速代替实际流速 带来的误差。 1．材料和方法 1．1实验设备 实验用土柱高为15cm，直径为5cm， 规格相同，都为有机玻璃制成。供水设备 图1试验装置示意图 为马氏瓶，规格为6升。为保持水位恒定， 采用马氏瓶供水。出口用软管和一量杯连接。实验装置如图1。 表1土壤的物理化学性质 EC CEC