吸附性土壤溶质运移参数识别的粒子群-差分算法.pdfVIP

第50卷第3期 土 壤 学 报 V01．50．No．3 2013年5月 ACTAPEDOLOGICASINICA May，2013 吸附性土壤溶质运移参数识别的粒子群．差分算法半 任长江1 白 丹” 周蓓蓓1 梁 伟2 陈 (1西安理工大学水利水电学院，西安710048) (2陕西师范大学旅游与环境学院，西安710062) 摘 要 非饱和土壤水分和溶质运移参数(导水率、扩散率、水动力弥散系数)的识别是进行数值模拟 的关键。基于垂直一维非饱和吸附性土壤水分-溶质运移方程，以溴化钾为人渗溶液，分别以土壤含水率和溶 质溶度的实测值与计算值标准差最小为优化目标，建立水分一溶质运移参数识别的多目标优化模型。应用权 重系数法，将这一多目标优化问题转化为单目标优化问题，采用动态权重、异步时变学习因子的粒子群算法 对模型求解。通过土柱实验，测定了220、380和780rain三组时刻的土壤含水率以及钾离子浓度的空间分布 值，以前两组数据识别参数，第三组数据进行验证。结果表明人渗780min权重系数为0．5时含水率以及钾 离子浓度的计算值和实测值相关系数分别为0．977和0．952，标准差分别为0．0074、2．369，实测值和计算值 吻合较好，这表明粒子群识别水分和溶质运移参数是可行的。数值模拟表明，权重越大含水率实测和计算值 的相关系数越大而浓度实测值和计算值的相关系数越小。 关键词 水分运移；溶质运移；粒子群算法；有限差分法；参数识别；多目标函数 中图分类号$278 文献标识码 A 土壤水分以及溶质运移方程参数的确定是土 Li_”1、崔凯等。1“、Li等’1引利用同伦正则化法和参数 壤水分和溶质运移数值模拟的基础。目前对土壤 HY等¨4。采用遗传 迭代优化算法，杨坤等…。，Sun 水分和溶质参数的求解主要有：借助现代测量设 算法对土壤水分参数进行了数值反演，数值反演在 备、模型模拟软件、数值计算、数值反演等方法。 解决此类问题中应用较广，但也存在算法不收敛或 Kandasamy等…利用高光谱影像技术通过对土壤光者陷入局部最优解等问题。因此寻求设备简单、方 谱的分析，Mohamed和Marwan口1利用TDR(Time法简便、计算精确、算法收敛的计算水分和溶质运 Domain Reflectometry)对土壤介电常数的测定，Can—移参数方法是本文要研究的内容。 demir等旧3通过对土壤物理特性的测定等方法间接 地计算出了土壤水分运移相关参数。现代测量设 1 非饱和水分一溶质运移方程及其 备能够快速测量水分以及溶质运移参数，但高昂的 参数 设备费用使其应用受到限制；刘建军等H1利用 Hydrus一1D软件，Zhang口1采用Mike2一D软件包，张定水头定浓度条件下一维垂直非饱和吸附性 嘉和工明玉¨1利用MT3DMS软件分别对一维、二水分．溶质运移方程…1： 维、三维等水分、溶质运移参数进行了求解。模拟 软件能够快速的对研究问题给出模拟结果，但模拟 面00=小㈩外丁OK(O) 的精度易受软件运行参数设置的影响。Duan和 Cui1基于小波函数用Wavelet—Galerkin法，杨文涛 (·+等]掣=小寸掣㈩ 等’81利用Gm—Cxtfit方法，王康等p1采用粒子平均移 0=0。，c=C。 (t=0，髫) 动和分形理论等数值方法计算土壤溶质运移参数。 (t，z=0) 数值法求解此类问题具有一定的精度，但也