不同数学模型模拟非均质土壤溶质大尺度运移对比分析.pdfVIP

水 利 学 报 2010年2月 SHUILI XUEBAO 第41卷 第2期 文章编号：0559-9350（2009）02-0164-09 不同数学模型模拟非均质土壤溶质大尺度运移的对比分析 1 1 1，2 3 高光耀 ，冯绍元 ，黄冠华 ，詹红兵 （1. 中国农业大学 中国农业水问题研究中心，北京 100083； 2. 中国农业大学 中国—以色列国际农业研究培训中心，北京 100083； 3. 美国得克萨斯农工大学 地质与地球物理系, 美国 德州 77843） 摘要：利用长度为1250cm的非均质土柱中NaCl运移的试验资料，分别应用对流－弥散方程（CDE）、两区模型 （TRM）、对流对数正态传递函数模型（CLT）、分数微分对流 弥散方程（FADE）和连续时间随机游走理论（CTRW） - 等数学模型对长土柱中溶质的运移过程进行模拟，对比各个模型的模拟结果，并从模型机理上分析模拟结果差异 的原因，探讨各种模型描述非均质土壤中溶质大尺度运移的适用性。结果表明：CDE、CLT和FADE难以描述非 均质土柱中的不规则穿透曲线，但FADE对穿透曲线尾部的模拟精度比CDE和CLT有所提高；与CDE、CLT和 FADE相比，两者的模拟结果与实测值吻合更好，在穿透曲线尾部CTRW的模拟值略大于TRM的模拟结果；两者 不仅较好地描述了穿透曲线的提前穿透和拖尾等不规则特征，而且捕捉到了非均质土柱中溶质不规则运移程度随 着尺度增加而减弱的变化趋势。尽管TRM和CTRW均能很好地模拟长土柱试验所研究的饱和非均质土壤中保守 性溶质的一维运移问题，但对于复杂的溶质运移问题的模拟，TRM和CTRW各有其优缺点。 关键词：溶质运移；模型比较；非均质土壤；大尺度；数值模拟 中图分类号：S152 文献标识码：A 1 研究背景 溶质在土壤中的运移过程非常复杂，溶质的分布不仅与水流运动和土壤结构有关，而且与溶质 ［1］ 的性质以及溶质随水流运动过程中所发生的一系列物理、化学和生物化学反应过程有密切关系 。数 学模型能够定量地描述溶质运移过程的本质规律，而且具有预测功能，是研究土壤溶质运移问题的 有效手段。长期以来，许多学者通过对土壤结构和溶质运移特征进行不同程度的概化，提出各种理 论方法和数学模型对土壤溶质运移过程进行模拟研究。 - 基于水动力弥散理论的对流 弥散方程（Convection-DispersionEquation，CDE）是常用的研究溶质 ［2］ 运移问题的数学模型 。利用CDE来模拟均质土壤中溶质运移过程往往能够得到较好的结果，但 CDE对于非均质土壤中溶质运移经常出现的提前穿透和拖尾等不规则现象却难以较好地描述,影响了 ［3-5］ 溶质运移预测的可靠性 。在CDE基础上建立起来的土壤溶质运移两区模型（Two-RegionModel， TRM）考虑非均质土壤孔隙中存在的不动水体在溶质运移过程中所起的作用，更真实地反映了土壤孔 隙流速分布特征，对上述的两种不规则现象能够进行比较精确的描述，被广泛用来研究室内的饱和 ［6-9］ 与非饱和土壤中溶质运移规律 。相比于CDE和TRM等以溶质运移机理和过程为构造基础的确定性 数学模型，许多学者认为土壤中溶质的运移现象是土壤内在的不稳定过程，适合用描述随机过程的 方法来研究其数量特征［10-11］。Jury提出的传递函数模型（TransferFunctionModel，TFM）就是一种简单 ［11］ 且常用的随机模型 。TFM采用溶质运移时间的概率密度函数表示溶质在土壤中