基于集合卡尔曼滤波的灌域尺度地下水多参数联合反演.pdfVIP

2018 年5 月 灌溉排水学报 第37 卷 第5 期 May 2018 Journal of Irrigation and Drainage No.5 Vol.37 文章编号：1672-3317（2018 ）05 - 0066 - 09 基于集合卡尔曼滤波的灌域尺度地下水多参数联合反演 * 张玉雪，朱 焱 ，杨金忠 （武汉大学水资源与水电工程科学国家重点实验室，武汉 430072 ） 摘 要：【目的】进行灌域尺度的地下水多参数联合反演。【方法】利用集合卡尔曼滤波（EnKF ）方法，基于二维潜水 运动模型SWF2D ，构建了地下水多参数联合反演模型SWF2D_DA 。以内蒙古河套灌区永济灌域为研究对象对地 下水参数进行了反演分析，该地区地下水运动共涉及8 个参数，分别为给水度Sy ，降雨入渗补给系数a ，两阶段的灌 P 溉入渗补给系数aI 1 、aI2 ，两阶段的潜水蒸发参数aA 1 、aA2 和两阶段的潜水蒸发极限埋深aS 1 、aS2 ，采用SWF2D_DA 模型 分别反演了2 待定参数情况和8待定参数情况。【结果】模型对2 参数反演平均相对误差MRE 为0.124% ，均方根误差 RMSE 为0.002 663 ，8 参数反演MRE 为0.376% ，RMSE 为0.003 283 ，反演结果均满足精度要求，同时反演8 参数会增 2 加得到理想参数所需的同化步。同时，还设置5 种观测误差方差，分别为0.01 、0.000 1、0.002 5、0.1、0.64 m ，讨论了 2 不同观测误差方差对反演结果的影响。当观测误差方差大于0.1 m 时，随着观测误差增大，MRE 、RMSE 、Spread 增 2 长较快，即观测误差过大会影响数据同化的反演精度；当观测误差方差在0.000 1～0.01 m 之间时，随着观测误差减 小，MRE 、RMSE 并无显著变化，即观测误差达到一定精度时，即使观测误差减小，模型也不能得到精度更高的解，反 而增加了观测成本。【结论】模型SWF2D_DA 可以实现大尺度复杂地区的多参数联合反演，且待求参数越少，得到 2 可靠的结果所需同化步越少；增加参数个数，需增加同化步。在目前的观测条件下，在观测误差方差为0.01 m 、相 对于地下水位最大变幅的相对观测误差在4.4% 以内的情况下，反演得到的地下水参数可满足精度要求，同时观测 成本较低。 关 键 词：多参数联合反演；数据同化；观测误差；集合卡尔曼滤波；地下水；模型 中图分类号：S274.1 文献标志码：A