Application of the bias aware Ensemble Kalman Filter with Confirming Option (Bias-CEnKF) in groundwater flow data assimilation

YANG Yun; WU Jichun; LUO Qiankun; QIAN Jiazhong

doi:10.16030/j.cnki.issn.1000-3665.202112017

2022 Vol. 49, No. 6

Article Contents

Article navigation > Hydrogeology & Engineering Geology > 2022 > 49(6): 13-23

Next Article Previous Article

YANG Yun, WU Jichun, LUO Qiankun, QIAN Jiazhong. Application of the bias aware Ensemble Kalman Filter with Confirming Option (Bias-CEnKF) in groundwater flow data assimilation[J]. Hydrogeology & Engineering Geology, 2022, 49(6): 13-23. doi: 10.16030/j.cnki.issn.1000-3665.202112017

Citation:

YANG Yun, WU Jichun, LUO Qiankun, QIAN Jiazhong. Application of the bias aware Ensemble Kalman Filter with Confirming Option (Bias-CEnKF) in groundwater flow data assimilation[J]. Hydrogeology & Engineering Geology, 2022, 49(6): 13-23. doi: 10.16030/j.cnki.issn.1000-3665.202112017

Application of the bias aware Ensemble Kalman Filter with Confirming Option (Bias-CEnKF) in groundwater flow data assimilation

1.
Department of Hydrosciences, School of Earth Sciences and Engineering, Nanjing University, Nanjing, Jiangsu　210023, China
2.
Huai River Water Resources Commission, Bengbu, Anhui　233001, China
3.
School of Resources and Environmental Engineering, Hefei University of Technology, Hefei, Anhui　230009, China

More Information

Corresponding author: WU Jichun, jcwu@nju.edu.cn

Received Date: 09 December 2021

Revised Date: 23 January 2022

Publish Date: 15 November 2022

Abstract

Abstract

The Ensemble Kalman Filter (EnKF) has been widely applied for real-time simulation of groundwater flow and solute transport. The majority of previous studies tend to assume no bias in forecast models, therefore ignoring the model uncertainties. This assumption, however, may be invalid when a conceptual model is not accurately generalized. As a result, forecast bias will lead to incorrect estimation of the system parameters or states. In this work, a bias aware Ensemble Kalman Filter with Confirming Option (Bias-CEnKF) is proposed to take into account the forecast bias by the model during the filtering process. The proposed method is tested in a real-time groundwater simulation considering model uncertainties, by setting inaccurate boundary conditions, initial conditions and recharge items. The results show that the standard EnKF may lead to filter divergence and assimilation failure, if the forecast model is not accurately generalized. Instead, Bias-CEnKF not only achieves better performances, but also reduces the inconsistency caused by the nonlinear relationship among the parameters, variables and bias corrections. Four scenarios are investigated, with the results showing the aquifer hydraulic conductivities and heads obtained by Bias-CEnKF are close to the real values, and the prediction results are also more reliable. This study further improves the applicability of the EnKF under the uncertain condition of model generalization, and provides a reliable method for groundwater data assimilation under complex field conditions.
- data assimilation /
- Ensemble Kalman Filter /
- conceptual model uncertainty /
- forecast bias /
- inconsistency

FullText(HTML)

References

[1]	EVENSEN G. Sequential data assimilation with a nonlinear quasi-geostrophic model using Monte Carlo methods to forecast error statistics[J]. Journal of Geophysical Research:Oceans,1994,99(5):10143 − 10162. doi: 10.1029/94JC00572 CrossRef Google Scholar
[2]	EVENSEN G. The Ensemble Kalman Filter:Theoretical formulation and practical implementation[J]. Ocean Dynamics,2003,53(4):343 − 367. doi: 10.1007/s10236-003-0036-9 CrossRef Google Scholar
[3]	沈晔,李海涛,黎涛,等. 地下水位预测: 集合卡尔曼滤波(EnKF)应用概述[J]. 水文地质工程地质,2014,41(1):21 − 24. [SHEN Ye, LI Haitao, LI Tao, et al. Groundwater level forecast: Overview of application of the Ensemble Kalman filter (EnKF)[J]. Hydrogeology & Engineering Geology,2014,41(1):21 − 24. (in Chinese with English abstract) Google Scholar
[4]	兰天, 康学远, 施小清, 等. 基于EnKF综合水头和浓度观测数据推估地下水流模型参数[J]. 水文地质工程地质,2017,44(5):6 − 13. [LAN Tian, KANG Xueyuan, SHI Xiaoqing, et al. Joint assimilation of heads and concentrations for estimating parameters of groundwater flow models using the Ensemble Kalman Filter[J]. Hydrogeology & Engineering Geology,2017,44(5):6 − 13. (in Chinese with English abstract) Google Scholar
[5]	刘琨, 黄冠华. 基于集合卡尔曼滤波法的二维土壤水流状态变量和参数联合估计[J]. 水利学报,2019,50(3):399 − 408. [LIU Kun, HUANG Guanhua. Joint state and parameter estimation of two-dimensional soil water flow model based on Ensemble Kalman Filter method[J]. Journal of Hydraulic Engineering,2019,50(3):399 − 408. (in Chinese with English abstract) Google Scholar
[6]	XU T, GÓMEZ-HERNÁNDEZ J J, CHEN Z, et al. A comparison between ES-MDA and restart EnKF for the purpose of the simultaneous identification of a contaminant source and hydraulic conductivity[J]. Journal of Hydrology,2021,595:125681. doi: 10.1016/j.jhydrol.2020.125681 CrossRef Google Scholar
[7]	HOUTEKAMER P L, ZHANG F Q. Review of the ensemble Kalman filter for atmospheric data assimilation[J]. Monthly Weather Review,2016,144(12):4489 − 4532. doi: 10.1175/MWR-D-15-0440.1 CrossRef Google Scholar
[8]	JORDI A, WANG D P. Estimation of transport at open boundaries with an ensemble Kalman filter in a coastal ocean model[J]. Ocean Modelling,2013,64:56 − 66. doi: 10.1016/j.ocemod.2013.01.002 CrossRef Google Scholar
[9]	EMERICK A, REYNOLDS A. Combining sensitivities and prior information for covariance localization in the ensemble Kalman filter for petroleum reservoir applications[J]. Computational Geosciences,2011,15(2):251 − 269. doi: 10.1007/s10596-010-9198-y CrossRef Google Scholar
[10]	崔凯鹏, 吴吉春. 观测数据时空密度对集合卡尔曼滤波计算精度的影响[J]. 水利学报,2013,44(8):915 − 923. [CUI Kaipeng, WU Jichun. Effect of observation data time/spatial density on Ensemble Kalman Filter[J]. Journal of Hydraulic Engineering,2013,44(8):915 − 923. (in Chinese with English abstract) Google Scholar
[11]	曹少华, 曾献奎, 蒋建国, 等. 基于局域化集合卡尔曼滤波的含水层弥散度场识别研究[J]. 南京大学学报(自然科学版),2016,52(3):429 − 437. [CAO Shaohua, ZENG Xiankui, JIANG Jianguo, et al. An approach based on localizaed ensemble Kalman filter to identify groundwater dispersivity field[J]. Journal of Nanjing University (Natural Sciences),2016,52(3):429 − 437. (in Chinese with English abstract) Google Scholar
[12]	康学远, 施小清, 邓亚平, 等. 基于EnKF融合地球物理数据刻画含水层非均质性[J]. 水科学进展,2018,29(1):40 − 49. [KANG Xueyuan, SHI Xiaoqing, DENG Yaping, et al. Assimilation of hydrogeophysical data for the characterization of subsurface heterogeneity using Ensemble Kalman Filter (EnKF)[J]. Advances in Water Science,2018,29(1):40 − 49. (in Chinese with English abstract) Google Scholar
[13]	南统超, 吴吉春. 集合卡尔曼滤波估计水文地质参数的局域化修正[J]. 水科学进展,2010,21(5):613 − 621. [NAN Tongchao, WU Jichun. Localization corrections for the estimation of hydrogeological parameters using ensemble Kalman filter[J]. Advances in Water Science,2010,21(5):613 − 621. (in Chinese with English abstract) Google Scholar
[14]	张秋汝, 史良胜, 林琳, 等. 非饱和土壤水的集合卡尔曼滤波Ⅱ:不一致性问题解决方法比较[J]. 水利学报,2015,46(12):1470 − 1478. [ZHANG Qiuru, SHI Liangsheng, LIN Lin, et al. Ensemble Kalman filter for unsaturated soil water flowⅡ:Comparison of methods to deal with inconsistency[J]. Journal of Hydraulic Engineering,2015,46(12):1470 − 1478. (in Chinese with English abstract) Google Scholar
[15]	周海燕. 基于集合卡尔曼滤波法的非高斯含水层参数识别[D]. 北京: 中国地质大学（北京）, 2012 Google Scholar ZHOU Haiyan. Characterizing non-Gaussian aquifer model parameters based on the ensemble Kalman filter[D]. Beijing: China University of Geosciences, 2012. （in Chinese with English abstract） Google Scholar
[16]	XU T, GÓMEZ-HERNÁNDEZ J J. Simultaneous identification of a contaminant source and hydraulic conductivity via the restart normal-score ensemble Kalman filter[J]. Advances in Water Resources,2018,112:106 − 123. doi: 10.1016/j.advwatres.2017.12.011 CrossRef Google Scholar
[17]	钟乐乐, 曾献奎, 吴吉春. 基于高斯过程回归的地下水模型结构不确定性分析与控制[J]. 水文地质工程地质,2019,46(1):1 − 10. [ZHONG Lele, ZENG Xiankui, WU Jichun. Quantification and reduction of groundwater model structural uncertainty based on Gaussian process regression[J]. Hydrogeology & Engineering Geology,2019,46(1):1 − 10. (in Chinese with English abstract) Google Scholar
[18]	ROJASR, FEYENL, DASSARGUES A. Conceptual model uncertainty in groundwater modeling: combining generalized likelihood uncertainty estimation and Bayesian model averaging[J]. Water Resources Research,2008,44(12):W12418. Google Scholar
[19]	ELSHALLAS, YE M, FINKEL M. Evaluating two multi-model simulation-optimization approaches for managing groundwater contaminant plumes[J]. Journal of Hydrology,2020,590:125427. doi: 10.1016/j.jhydrol.2020.125427 CrossRef Google Scholar
[20]	孙晓卓, 曾献奎, 吴吉春, 等. 一种改进的地下水模型结构不确定性分析方法[J]. 水文地质工程地质,2021,48(6):24 − 33. [SUN Xiaozhuo, ZENG Xiankui, WU Jichun, et al. An improved method of groundwater model structural uncertainty analysis[J]. Hydrogeology & Engineering Geology,2021,48(6):24 − 33. (in Chinese with English abstract) Google Scholar
[21]	ERDAL D, NEUWEILER I, WOLLSCHLÄGER U. Using a bias aware EnKF to account for unresolved structure in an unsaturated zone model[J]. Water Resources Research,2014,50(1):132 − 147. doi: 10.1002/2012WR013443 CrossRef Google Scholar
[22]	XUE L, ZHANG D X. A multimodeldata assimilation framework via the ensemble Kalman filter[J]. Water Resources Research,2014,50(5):4197 − 4219. doi: 10.1002/2013WR014525 CrossRef Google Scholar
[23]	REFSGAARD J C, VAN DER SLUIJS J P, BROWN J, et al. A framework for dealing with uncertainty due to model structure error[J]. Advances in Water Resources,2006,29(11):1586 − 1597. doi: 10.1016/j.advwatres.2005.11.013 CrossRef Google Scholar
[24]	KOLLAT J B, REED P M, RIZZO D M. Addressing model bias and uncertainty in three dimensional ground water transport forecasts for a physical aquifer experiment[J]. Geophysical Research Letters,2008,35(17):L17402. doi: 10.1029/2008GL035021 CrossRef Google Scholar
[25]	WEN X H, CHEN W H. Real-time reservoir model updating using ensemble Kalman filter with confirming option[J]. SPE Journal,2006,11(4):431 − 442. doi: 10.2118/92991-PA CrossRef Google Scholar
[26]	崔凯鹏. 集合卡尔曼滤波在地下水流及溶质运移数据同化中的应用探讨[D]. 南京: 南京大学, 2013 Google Scholar CUI Kaipeng. The application and discussion of ensemble Kalman filter method in groundwater flow and solute transport data assimilation[D]. Nanjing: Nanjing University, 2013. （in Chinese with English abstract） Google Scholar
[27]	康学远, 施小清, 史良胜, 等. 基于集合卡尔曼滤波的多相流模型参数估计—以室内二维砂箱中重质非水相污染物入渗为例[J]. 吉林大学学报(地球科学版),2017,47(3):848 − 859. [KANG Xueyuan, SHI Xiaoqing, SHI Liangsheng, et al. Inverse multiphase flow simulation using ensemble Kalman filter:Application to a 2D sandbox experiment of DNAPL migration[J]. Journal of Jilin University (Earth Science Edition),2017,47(3):848 − 859. (in Chinese with English abstract) Google Scholar
[28]	ZHANG D X, LU Z M. An efficient, high-order perturbation approach for flow in random porous media via Karhunen-Loève and polynomial expansions[J]. Journal of Computational Physics,2004,194(2):773 − 794. doi: 10.1016/j.jcp.2003.09.015 CrossRef Google Scholar
[29]	HARBAUGH A W, BANTA E R, HILL M C, et al. MODFLOW-2000, the U. S. Geological survey modular ground-water model-user guide to modularization concepts and the ground-water flow process[R]. Reston: USGS, 2000. Google Scholar