非饱和土壤水的集合卡尔曼滤波Ⅱ不一致性问题解决方法.pdf

水 利 学 报 2015年12月 SHUILI XUEBAO 第46卷 第12期 文章编号：0559-9350（2015）12-1470-09 非饱和土壤水的集合卡尔曼滤波Ⅱ： 不一致性问题解决方法比较 张秋汝，史良胜，林 琳，杨楚慧 （武汉大学 水资源与水电工程科学国家重点实验室，湖北 武汉 430072） 摘要：集合卡尔曼滤波（EnKF）能够便捷地根据观测信息改进土壤水的预测精度，但在非线性问题中存在不一致 性问题：当同时更新变量（水头或含水率）和非饱和水力参数时，EnKF会导致参数和变量之间不再服从Richards 方程关系。本文以Borden试验场的饱和渗透系数数据为基础，构造土柱试验，研究了不一致性对土壤水数据同 化带来的破坏以及相应的迭代型解决方案。研究结果表明：在非均匀土壤中，EnKF可能会引发强烈的不一致 性，对含水率预测和参数估计的精度造成破坏；不一致性的峰值位置与水流锋面保持一致，且受边界条件影响； 对于同类型的观测与待求参数，不一致性一般是随着观测点数量与待求参数数量的比例的增大而减小；当观测与 水头或含水率之间具有强烈的非线性时（如强烈的干湿交替），推荐采用能在保证计算效率的前提下有效降低不一 致性影响的CEnKF或MREnKF方法；当观测信息相对充足，不一致性可忽略时，EnKF方法更具优势。 关键词：非饱和流；集合卡尔曼滤波；不一致性；迭代型方法 中图分类号：S152.7 文献标识码：A doi ：10.13243/ki.slxb 1 研究背景 ［1］ 针对非饱和土壤水数据同化问题，在“非饱和土壤水的集合卡尔曼滤波 Ⅰ” 中，我们已对状态向 量和非饱和流算法的选择进行了研究，提出联合更新参数和变量（水头）的必要性并推荐采用非迭代 型非饱和流算法—Ross算法。但由于集合卡尔曼滤波（EnKF）采用了高斯分布的假设，在线性问题中 ［2］ 该方法是完全严格的，而在强烈非线性问题中存在理论和方法上的缺陷 。随着数据同化技术在石油 工程、气象学和水文等领域的持续应用，不少研究者注意到非线性可能带来的破坏。例如，Gu等［3］ 发现EnKF可能更新出没有物理意义的数值，出现饱和度超过1或者小于0的情况。 非饱和土壤水运动受非线性Richards方程的控制，同时土壤水分运动往往受到快速变化的气象条 ［4］ 件的影响，因此非饱和水分运动通常呈现出非常强烈的非线性特征 。非线性给土壤水数据同化方法 带来的主要问题是不一致性：当应用增广状态向量时，在每一个同化步，EnKF方法需要同时对变量 （水头或含水率）和非饱和水力参数进行更新，更新时依赖的是由统计样本得到的观测与变量、参数 的相关关系，相当于在观测时间点对变量和参数同时做一次线性校正；由于这个校正完全由统计样 本得到，并没有考虑到参数和变量之间的非线性关系，更新后的参数和变量之间不再满足Richards方 ［5］ ［6］ 程 。在其它领域中，不一致性问题也得以关注。Wen等 在石油库藏模拟中发现EnKF会导致静态 ［7］ 参数与动态变量之间的不一致性；Hendricks等 比较了EnKF和迭代型EnKF在饱和流中的计算性 能，发现迭代型方法略优于原始EnKF，但对于非饱和流问题，迭代型方法的潜在效益仍有待深入研 究。为缓解或消除不一致性带来的不利影响，多种迭代型EnKF得以提出，一般可划分为如下两类： 收稿日期：2014-10-13 基金项目：国家自然科学基金项目51279141） 作者简介：张秋汝（1990-），女，湖北随