Large‐scale inverse modeling with an application in hydraulic tomography

[1] Inverse modeling has been widely used in subsurface problems, where direct measurements of parameters are expensive and sometimes impossible. Subsurface media are inherently heterogeneous in complex ways, which implies that the number of unknowns is usually large. Furthermore, technologies such as hydraulic tomography and electric resistivity tomography allow the collection of more indirect measurements, and at the same time, there is an increased appreciation of the value of detailed characterization of the subsurface media in, for example, remediation projects. Hence, we need efficient inverse methods that can assimilate a large volume of measurements to estimate even larger numbers of parameters, i.e., large-scale inverse modeling. In this paper, we present a Bayesian method that employs a sparse formulation, and we applied this method to a laboratory hydraulic tomography problem, where we successfully estimated half a million unknowns that represent the hydraulic conductivity field of the sandbox at a fine scale. The inversion took about 2 h with a single core.

[1]  Jet-Chau Wen,et al.  A simultaneous successive linear estimator and a guide for hydraulic tomography analysis , 2009 .

[2]  A. Sahuquillo,et al.  Stochastic simulation of transmissivity fields conditional to both transmissivity and piezometric data—I. Theory , 1997 .

[3]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[4]  Minghui Jin,et al.  AN ITERATIVE STOCHASTIC INVERSE METHOD: CONDITIONAL EFFECTIVE TRANSMISSIVITY AND HYDRAULIC HEAD FIELDS , 1995 .

[5]  P. Kitanidis,et al.  Geostatistical inversing for large-contrast transmissivity fields , 2009 .

[6]  A Revil,et al.  A Potential‐Based Inversion of Unconfined Steady‐State Hydraulic Tomography , 2009, Ground water.

[7]  Walter A Illman,et al.  Practical Issues in Imaging Hydraulic Conductivity through Hydraulic Tomography , 2008, Ground water.

[8]  Geir Nævdal,et al.  Reservoir Monitoring and Continuous Model Updating Using Ensemble Kalman Filter , 2005 .

[9]  Philip E. Gill,et al.  Practical optimization , 1981 .

[10]  Reply [to “Comment on ‘A reassessment of the groundwater inverse problem’ by D. McLaughlin and L. R. Townley”] , 1997 .

[11]  E. G. Vomvoris,et al.  A geostatistical approach to the inverse problem in groundwater modeling (steady state) and one‐dimensional simulations , 1983 .

[12]  P. Kitanidis,et al.  Estimation of historical groundwater contaminant distribution using the adjoint state method applied to geostatistical inverse modeling , 2004 .

[13]  Y. Rubin,et al.  Analysis of the Cape Cod tracer data , 1997 .

[14]  J. Scales,et al.  Resolution of seismic waveform inversion: Bayes versus Occam , 1997 .

[15]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[16]  Junfeng Zhu,et al.  Laboratory sandbox validation of transient hydraulic tomography , 2007 .

[17]  P. Kitanidis,et al.  A geostatistical approach to contaminant source identification , 1997 .

[18]  D. A. Zimmerman,et al.  A comparison of seven geostatistically based inverse approaches to estimate transmissivities for modeling advective transport by groundwater flow , 1998 .

[19]  Peter K. Kitanidis,et al.  Generalized covariance functions associated with the Laplace Equation and Their use in interpolation and inverse problems , 1999 .

[20]  Andres Alcolea,et al.  Inverse problem in hydrogeology , 2005 .

[21]  Hiromitsu Saegusa,et al.  Hydraulic tomography in fractured granite: Mizunami Underground Research site, Japan , 2009 .

[22]  Yoram Rubin,et al.  Prediction of tracer plume migration in disordered porous media by the method of conditional probabilities , 1991 .

[23]  Gene H. Golub,et al.  Calculating the singular values and pseudo-inverse of a matrix , 2007, Milestones in Matrix Computation.

[24]  R. Parker,et al.  Occam's inversion; a practical algorithm for generating smooth models from electromagnetic sounding data , 1987 .

[25]  W. Nowak,et al.  Geostatistical inference of hydraulic conductivity and dispersivities from hydraulic heads and tracer data , 2006 .

[26]  P. Kitanidis,et al.  An Application of the Geostatistical Approach to the Inverse Problem in Two-Dimensional Groundwater Modeling , 1984 .

[27]  G. Dagan.  Stochastic Modeling of Groundwater Flow by Unconditional and Conditional Probabilities: The Inverse Problem , 1985 .

[28]  Walter A. Illman,et al.  Steady-state hydraulic tomography in a laboratory aquifer with deterministic heterogeneity: Multi-method and multiscale validation of hydraulic conductivity tomograms , 2007 .

[29]  W. Gragg,et al.  The generalized Schur algorithm for the superfast solution of Toeplitz systems , 1987 .

[30]  F. Hoog.  A new algorithm for solving Toeplitz systems of equations , 1987 .

[31]  P. Kitanidis.  Quasi‐Linear Geostatistical Theory for Inversing , 1995 .

[32]  W. Nowak,et al.  Geostatistical inverse modeling of transient pumping tests using temporal moments of drawdown , 2005 .

[33]  Wolfgang Nowak,et al.  Efficient Computation of Linearized Cross-Covariance and Auto-Covariance Matrices of Interdependent Quantities , 2003 .

[34]  W. Yeh.  Review of Parameter Identification Procedures in Groundwater Hydrology: The Inverse Problem , 1986 .

[35]  C. Paige.  Bidiagonalization of Matrices and Solution of Linear Equations , 1974 .

[36]  P. Kitanidis.  Introduction to Geostatistics: Applications in Hydrogeology , 1997 .