Data assimilation methods for estimating a heterogeneous conductivity field by assimilating transient solute transport data via ensemble Kalman filter

被引:12
作者
Tong, Juxiu [1 ,2 ,3 ]
Hu, Bill X. [1 ,2 ,3 ]
Yang, Jinzhong [4 ]
机构
[1] China Univ Geosci, Key Lab Groundwater Cycle & Environm Evolut, Minist Educ, Beijing 100083, Peoples R China
[2] China Univ Geosci, Collage Water Resources & Environm Sci, Beijing 100083, Peoples R China
[3] Florida State Univ, Dept Earth Ocean & Atmospher Sci Geol Sci, Tallahassee, FL 32306 USA
[4] Wuhan Univ, State Key Lab Water Resources & Hydropower Engn S, Wuhan 430072, Peoples R China
关键词
data assimilation; ensemble Kalman filter; hydraulic conductivity; solute transport; mixed Neumann; Dirichlet boundary; observation error; inflation factor; model error; MODEL; FLOW; CALIBRATION;
D O I
10.1002/hyp.9523
中图分类号
TV21 [水资源调查与水利规划];
学科分类号
081501 ;
摘要
An ensemble Kalman filter (EnKF) is developed to identify a hydraulic conductivity distribution in a heterogeneous medium by assimilating solute concentration measurements of solute transport in the field with a steady-state flow. A synthetic case with the mixed Neumann/Dirichlet boundary conditions is designed to investigate the capacity of the data assimilation methods to identify a conductivity distribution. The developed method is demonstrated in 2-D transient solute transport with two different initial instant solute injection areas. The influences of the observation error and model error on the updated results are considered in this study. The study results indicate that the EnKF method will significantly improve the estimation of the hydraulic conductivity field by assimilating solute concentration measurements. The larger area of the initial distribution and the more observed data obtained, the better the calculation results. When the standard deviation of the observation error varies from 1% to 30% of the solute concentration measurements, the simulated results by the data assimilation method do not change much, which indicates that assimilation results are not very sensitive to the standard deviation of the observation error in this study. When the inflation factor is more than 1.0 to enlarge the model error by increasing the forecast error covariance matrix, the updated results of the hydraulic conductivity by the data assimilation method are not good at all. Copyright (c) 2012 John Wiley & Sons, Ltd.
引用
收藏
页码:3873 / 3884
页数:12
相关论文
共 45 条
[1]  
[Anonymous], 1989, FLOW TRANSPORT POROU, DOI DOI 10.1007/978-3-642-75015-1
[2]  
Bear J., 1972, Dynamics of Fluids in Porous Media
[3]   Data assimilation for transient flow in geologic formations via ensemble Kalman filter [J].
Chen, Yan ;
Zhang, Dongxiao .
ADVANCES IN WATER RESOURCES, 2006, 29 (08) :1107-1122
[4]  
Constantinescu EM, 2006, Q J ROYAL METEOROLOG, DOI 10.1256/qj.yy.n
[5]  
Dai ZX, 2010, HYDROGEOL J, V18, P607, DOI 10.1007/s10040-009-0543-y
[6]   Inverse problem of multicomponent reactive chemical transport in porous media: Formulation and applications [J].
Dai, ZX ;
Samper, J .
WATER RESOURCES RESEARCH, 2004, 40 (07) :W074071-W0740718
[7]   Calibration framework for a Kalman filter applied to a groundwater model [J].
Drécourt, JP ;
Madsen, H ;
Rosbjerg, D .
ADVANCES IN WATER RESOURCES, 2006, 29 (05) :719-734
[8]   An Ensemble Kalman filter with a 1-D marine ecosystem model [J].
Eknes, M ;
Evensen, G .
JOURNAL OF MARINE SYSTEMS, 2002, 36 (1-2) :75-100
[9]   Sampling strategies and square root analysis schemes for the EnKF [J].
Evensen, G .
OCEAN DYNAMICS, 2004, 54 (06) :539-560
[10]  
Evensen G., 2006, Data Assimilation: The Ensemble Kalman Filter