Comparison of Two Ensemble Kalman-Based Methods for Estimating Aquifer Parameters from Virtual 2-D Hydraulic and Tracer Tomographic Tests

Abstract

We compare two ensemble Kalman-based methods to estimate the hydraulic conductivity field of an aquifer from data of hydraulic and tracer tomographic experiments: (i) the Ensemble Kalman Filter (EnKF) and (ii) the Kalman Ensemble Generator (KEG). We generated synthetic drawdown and tracer data by simulating two pumping tests, each followed by a tracer test. Parameter updating with the EnKF is performed using the full transient signal. For hydraulic data, we use the standard update scheme of the EnKF with damping, whereas for concentration data, we apply a restart scheme, in which solute transport is resimulated from time zero to the next measurement time after each parameter update. In the KEG, we iteratively assimilate all observations simultaneously, here inverting steady-state heads and mean tracer arrival times. The inversion with the dampened EnKF worked well for the transient pumping-tests, but less for the tracer tests. The KEG produced similar estimates of hydraulic conductivity but at significantly lower costs. We conclude that parameter estimation in well-defined hydraulic tests can be done very efficiently by iterative ensemble Kalman methods, and ambiguity between state and parameter updates can be completely avoided by assimilating temporal moments of concentration data rather than the time series themselves.