Operational real-time modeling with ensemble Kalman filter of variably saturated subsurface flow including stream-aquifer interaction and parameter updating

被引:71
作者
Franssen, H. J. Hendricks [1 ]
Kaiser, H. P. [2 ]
Kuhlmann, U. [3 ]
Bauser, G. [1 ]
Stauffer, F. [1 ]
Mueller, R. [2 ]
Kinzelbach, W. [1 ]
机构
[1] ETH, Inst Environm Engn, CH-8093 Zurich, Switzerland
[2] Water Supply Zurich, CH-8023 Zurich, Switzerland
[3] TK Consult, Zurich, Switzerland
关键词
MONTE-CARLO METHODS; DATA ASSIMILATION; STEADY-STATE; GROUNDWATER-FLOW; HYDRAULIC CONDUCTIVITY; TRANSIENT; TRANSPORT; FIELD; TRANSMISSIVITY; CALIBRATION;
D O I
10.1029/2010WR009480
中图分类号
X [环境科学、安全科学];
学科分类号
08 ; 0830 ;
摘要
Urban groundwater is frequently contaminated, and the exact location of the pollution spots is often unknown. Intelligent monitoring of the temporal variations in groundwater flow in such an area assists in selectively extracting groundwater of drinking water quality. Here an example from the city of Zurich (Switzerland) is shown. The monitoring strategy consists of using the ensemble Kalman filter (EnKF) for optimally combining online observations and online models for the real-time characterization of groundwater flow. We conducted numerical simulation experiments for the period January 2004 to December 2007 with a 3-D finite element model for variably saturated groundwater flow. It was found that the daily assimilation of piezometric head data with EnKF results in a better characterization of piezometric heads than does a model which is inversely calibrated with historical data but not updated in real time. The positive impact of model updating with observations can still be observed 10 days after the update. These simulations also suggest that parameters (hydraulic conductivity and leakage) are successfully updated: 1 and 10 day piezometric head predictions are better with than without updating of parameters. Additional experiments with a synthetic model for the same site, in which the only difference is that certain parameter values are selected as the unknown "true" conditions, show that EnKF also successfully updates unknown parameters. However, this is only the case if spatially distributed hydraulic conductivities and leakage coefficients are jointly updated and if a damping parameter is used. The mean absolute error of estimated log leakage coefficients decreased by up to 63%; for log hydraulic conductivity a decrease of up to 27% was observed. From January 2009 the method has been operational at the Water Works Zurich and showed a remarkable performance until present (October 2010).
引用
收藏
页数:20
相关论文
共 46 条
[1]   Pilot points method incorporating prior information for solving the groundwater flow inverse problem [J].
Alcolea, Andres ;
Carrera, Jesus ;
Medina, Agustin .
ADVANCES IN WATER RESOURCES, 2006, 29 (11) :1678-1689
[2]  
Allen R.G., 1998, FAO Irrigation and Drainage Paper 56
[3]   Real-Time Management of an Urban Groundwater Well Field Threatened by Pollution [J].
Bauser, Gero ;
Franssen, Harrie-Jan Hendricks ;
Kaiser, Hans-Peter ;
Kuhlmann, Ulrich ;
Stauffer, Fritz ;
Kinzelbach, Wolfgang .
ENVIRONMENTAL SCIENCE & TECHNOLOGY, 2010, 44 (17) :6802-6807
[4]  
Bear J., 1979, Hydraulics of Groundwater
[5]  
Burgers G, 1998, MON WEATHER REV, V126, P1719, DOI 10.1175/1520-0493(1998)126<1719:ASITEK>2.0.CO
[6]  
2
[7]   ESTIMATION OF AQUIFER PARAMETERS UNDER TRANSIENT AND STEADY-STATE CONDITIONS .1. MAXIMUM-LIKELIHOOD METHOD INCORPORATING PRIOR INFORMATION [J].
CARRERA, J ;
NEUMAN, SP .
WATER RESOURCES RESEARCH, 1986, 22 (02) :199-210
[8]   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
[9]   Assessment of model uncertainty for soil moisture through ensemble verification [J].
De Lannoy, Gabrielle J. M. ;
Houser, Paul R. ;
Pauwels, Valentijn R. N. ;
Verhoest, Niko E. C. .
JOURNAL OF GEOPHYSICAL RESEARCH-ATMOSPHERES, 2006, 111 (D10)
[10]  
*DELT H ING MBH, 2006, SPRING 3 2 SOFTW