非饱和土中溶质迁移参数反演的HISR方法  

HISR Method for Parameter Inversion of Solute Transport through Unsaturated Soils

作　　者：崔凯[1] 李兴斯[1] 李宝元[2] 

机构地区：[1]大连理工大学工业结构装备分析国家重点实验室,辽宁大连116023 [2]大连理工大学工程力学系,辽宁大连116023

出　　处：《地球科学（中国地质大学学报）》2004年第2期247-252,共6页Earth Science-Journal of China University of Geosciences

基　　金：国家重点基础性研究发展规划项目 (G19990 3 2 80 5 )

摘　　要：以非饱和土中溶质迁移参数反演问题为背景 ,依据正则化方法的思路 ,以Itakura Saito距离作为同伦函数中的平凡问题 ,将同伦方法引入非线性参数反演问题的求解 ,进而提出一种求解非线性参数反演问题的大范围收敛 (HomotoyItakura SaitoRegularization ,HISR)方法 .为保证迭代稳定性 ,并同时削弱观测噪声的影响 ,同伦参数的修正采用了连续化修正方法 .本文将HISR方法应用于求解带有平衡及非平衡吸附效应的一维非饱和土中溶质迁移参数反演问题 ,计算结果表明HISR方法具有大范围收敛性及计算稳健性 ,同时有较强的抵抗观测噪声的能力 .Taken the parameter inversion of solute transport through unsaturated soils as the application background and homotopy method imported for solving the nonlinear parameter inversion problem, this paper presents a global convergence method, HISR (homotopy Itakura-Saito regularization) method. The Itakura-Saito distance between the iterative parameter value and the estimated parameter value was adopted as the trivial problem in the homotopy function based on the idea of regularization. Moreover, a continuous method is employed to modify the homotopy parameter during iteration in order to ensure the ability of the iteration and the capability of the noise resistance. Problems of parameter inversion of solute transport coupled with equilibrium and non-equilibrium effects through one-dimensional unsaturated soils are carried out as numerical examples, and the computational results clearly demonstrated the features of global convergence and stability of the HISR method. Besides, a favorable solution is obtained even though experimental quantities were contaminated heavily by noise by HISR method.

关 键 词：非饱和土 溶质迁移 非线性参数反演 同伦方法 Itakura-Saito距离 正则化 

分 类 号：P64[天文地球—地质矿产勘探]