Dynamic inversion for hydrological process monitoring with electrical resistance tomography under model uncertainties

Abstract

[1] We propose an approach for imaging the dynamics of complex hydrological processes. The evolution of electrically conductive fluids in porous media is imaged using time-lapse electrical resistance tomography. The related dynamic inversion problem is solved using Bayesian filtering techniques; that is, it is formulated as a sequential state estimation problem in which the target is an evolving posterior probability density of the system state. The dynamical inversion framework is based on the state space representation of the system which involves the construction of a stochastic evolution model and an observation model. The observation model that we use in this paper consists of the complete electrode model for ERT, with Archie's law relating saturations to electrical conductivity. The evolution model is an approximate model for simulating flow through partially saturated porous media. Unavoidable modeling and approximation errors in both the observation and evolution models are considered by computing approximate statistics for these errors. These models are then included in the construction of the posterior probability density of the estimated system state. This approximation error method allows the use of approximate, and therefore computationally efficient, observation and evolution models in the Bayesian filtering. We conside7r a synthetic example and show that the incorporation of an explicit model for the model uncertainties in the state space representation can yield better estimates than the frame-by-frame imaging approach.

Error handling strategies in multiphase inverse modeling

“…Second, the approximation error theory developed by Kaipio and Somersalo (2004) and demonstrated for a hydrogeophysical application by Lehikoinen et al (2010) will be implemented to address structural errors when highly simplified models are used during the inversion. Third, it is recognized (e.g., Moore and Doherty, 2006) that restricting the solution space by the zonation approach and other parameterization schemes may lead to calibrated models that have limited predictive power.…”

A truncated Levenberg–Marquardt algorithm for the calibration of highly parameterized nonlinear models

“…However, this interpretation assumes that the underlying conceptual model is error free, and that errors are independent. Since the residuals also include modeling errors, the validity of these assumptions is questionable, as demonstrated in Lehikoinen et al (2010) and Finsterle and Zhang (2010). From a practical perspective, the purpose of is to weight measurements of different quality, to scale observations of different types, to make the objective function dimensionless, and to weight the data fitting error relative to regularization terms.…”

Entropy-Bayesian Inversion of Time-Lapse Tomographic GPR data for Monitoring Dielectric Permittivity and Soil Moisture Variations

“…Chen, Rubin, Ma, & Baldocchi, 2008;Copty, Rubin, & Mavk, 1993;Dubreuil-Boisclair et al, 2011;Hou, Rubin, Hoversten, Vasco, & Chen, 2006;Hou & Rubin, 2005;Hubbard et al, 2001;M. B. Kowalsky et al, 2005;Michael B. Kowalsky et al, 2004;Lehikoinen, Huttunen, Finsterle, Kowalsky, & Kaipio, 2010) .…”