skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Large-scale inverse model analyses employing fast randomized data reduction: INVERSE MODELING WITH DATA REDUCTION

Journal Article · · Water Resources Research
DOI:https://doi.org/10.1002/2016WR020299· OSTI ID:1507327

When the number of observations is large, it is computationally challenging to apply classical inverse modeling techniques. We have developed a new computationally efficient technique for solving inverse problems with a large number of observations (e.g., on the order of 107 or greater). Our method, which we call the randomized geostatistical approach (RGA), is built upon the principal component geostatistical approach (PCGA). We employ a data reduction technique combined with the PCGA to improve the computational efficiency and reduce the memory usage. Specifically, we employ a randomized numerical linear algebra technique based on a so-called “sketching” matrix to effectively reduce the dimension of the observations without losing the information content needed for the inverse analysis. In this way, the computational and memory costs for RGA scale with the information content rather than the size of the calibration data. Our algorithm is coded in Julia and implemented in the MADS open-source high-performance computational framework (http://mads.lanl.gov). We apply our new inverse modeling method to invert for a synthetic transmissivity field. Compared to a standard geostatistical approach (GA), our method is more efficient when the number of observations is large. Moreover, our method is capable of solving larger inverse problems than the standard GA and PCGA approaches. Thus, our new model inversion method is a powerful tool for solving large-scale inverse problems. The method can be applied in any field and is not limited to hydrogeological applications such as the characterization of aquifer heterogeneity.

Research Organization:
Los Alamos National Laboratory (LANL), Los Alamos, NM (United States)
Sponsoring Organization:
USDOE Office of Science (SC)
Grant/Contract Number:
89233218CNA000001; SC11145687; 11145687
OSTI ID:
1507327
Alternate ID(s):
OSTI ID: 1374898
Report Number(s):
LA-UR-17-20239
Journal Information:
Water Resources Research, Vol. 53, Issue 8; ISSN 0043-1397
Publisher:
American Geophysical Union (AGU)Copyright Statement
Country of Publication:
United States
Language:
English
Citation Metrics:
Cited by: 19 works
Citation information provided by
Web of Science

References (53)

Inverse problem in hydrogeology journal February 2005
A modified Levenberg–Marquardt algorithm for quasi-linear geostatistical inversing journal July 2004
LSRN: A Parallel Iterative Solver for Strongly Over- or Underdetermined Systems journal January 2014
Extensions of Lipschitz mappings into a Hilbert space book January 1984
Scalable subsurface inverse modeling of huge data sets with an application to tracer concentration breakthrough data from magnetic resonance imaging: SCALABLE INVERSE MODELING WITH A HUGE MRI DATA SET journal July 2016
Large-scale stochastic linear inversion using hierarchical matrices: Illustrated with an application to crosswell tomography in seismic imaging journal July 2013
Fast iterative implementation of large-scale nonlinear geostatistical inverse modeling: FAST ITERATIVE NONLINEAR INVERSE MODELING journal January 2014
Quasi-Linear Geostatistical Theory for Inversing journal October 1995
Principal Component Geostatistical Approach for large-dimensional inverse problems journal July 2014
A fast randomized algorithm for overdetermined linear least-squares regression journal September 2008
Parameter estimation by ensemble Kalman filters with transformed data: Approach and application to hydraulic tomography: PARAMETER ESTIMATION BY tEnKFs journal April 2012
Active Subspace Methods in Theory and Practice: Applications to Kriging Surfaces journal January 2014
A computationally efficient parallel Levenberg-Marquardt algorithm for highly parameterized inverse model analyses: PARALLEL LEVENBERG-MARQUARDT FOR INVERSE MODELING journal September 2016
GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems journal July 1986
Analysis of hydraulic tomography using temporal moments of drawdown recovery data: HYDRAULIC TOMOGRAPHY journal February 2006
A statistical approach to the inverse problem of aquifer hydrology: 1. Theory journal August 1979
Should hydraulic tomography data be interpreted using geostatistical inverse modeling? A laboratory sandbox investigation: Geostatistics for hydraulic tomography? journal May 2015
A data-scalable randomized misfit approach for solving large-scale PDE-constrained inverse problems journal May 2017
Waveform inversion with source encoding for breast sound speed reconstruction in ultrasound computed tomography journal March 2015
UPRE method for total variation parameter selection journal August 2010
Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 1. Maximum Likelihood Method Incorporating Prior Information journal February 1986
Three-dimensional numerical inversion of pneumatic cross-hole tests in unsaturated fractured tuff: 2. Equivalent parameters, high-resolution stochastic imaging and scale effects journal December 2001
Faster least squares approximation journal October 2010
Three-dimensional numerical inversion of pneumatic cross-hole tests in unsaturated fractured tuff: 1. Methodology and borehole effects journal December 2001
An iterative geostatistical inverse method for steady flow in the vadose zone journal January 1997
Faster dimension reduction journal February 2010
Solution of Sparse Indefinite Systems of Linear Equations journal September 1975
Sensitivity of temporal moments calculated by the adjoint-state method and joint inversing of head and tracer data journal October 2000
A statistical approach to the inverse problem of aquifer hydrology: 2. Case study journal February 1980
Hydraulic tomography using temporal moments of drawdown recovery data: A laboratory sandbox study: LAB VALIDATION OF HYDRAULIC TOMOGRAPHY journal January 2009
Stochastic Fusion of Information for Characterizing and Monitoring the Vadose Zone journal January 2002
Implementing Randomized Matrix Algorithms in Parallel and Distributed Environments journal January 2016
Introduction to hierarchical matrices with applications journal May 2003
Efficient methods for large-scale linear inversion using a geostatistical approach: EFFICIENT GEOSTATISTICAL APPROACH journal May 2012
Fast full-wavefield seismic inversion using encoded sources journal November 2009
Geostatistical inference of hydraulic conductivity and dispersivities from hydraulic heads and tracer data: GEOSTATISTICAL INFERENCE journal August 2006
Blendenpik: Supercharging LAPACK's Least-Squares Solver journal January 2010
RandNLA: randomized numerical linear algebra journal May 2016
Are Models Too Simple? Arguments for Increased Parameterization journal May 2007
Sparser Johnson-Lindenstrauss Transforms journal January 2014
Large-scale hydraulic tomography and joint inversion of head and tracer data using the Principal Component Geostatistical Approach (PCGA) journal July 2014
Blind source separation for groundwater pressure analysis based on nonnegative matrix factorization journal September 2014
The minimum structure solution to the inverse problem journal October 1997
Geostatistical reduced-order models in underdetermined inverse problems: Groms in Underdetermined Inverse Problems journal October 2013
A hybrid regularized inversion methodology for highly parameterized environmental models: HYBRID REGULARIZATION METHODOLOGY journal October 2005
A new optimization approach for source-encoding full-waveform inversion journal May 2013
Numerical linear algebra in the streaming model conference January 2009
A sparse Johnson: Lindenstrauss transform conference January 2010
Breast ultrasound tomography with two parallel transducer arrays conference March 2016
Randomized Sketching for Large-Scale Sparse Ridge Regression Problems conference November 2016
Improved Approximation Algorithms for Large Matrices via Random Projections
  • Sarlos, Tamas
  • 2006 47th Annual IEEE Conference on Foundations of Computer Science, 2006 47th Annual IEEE Symposium on Foundations of Computer Science (FOCS'06) https://doi.org/10.1109/FOCS.2006.37
conference October 2006
Computational Methods for Inverse Problems book January 2002
Efficient implementation of ultrasound waveform tomography using source encoding conference June 2012

Cited By (5)

Randomized Truncated SVD Levenberg‐Marquardt Approach to Geothermal Natural State and History Matching journal March 2018
Robust Predictive Design of Field Measurements for Evapotranspiration Barriers Using Universal Multiple linear Regression journal November 2019
A Reduced‐Order Successive Linear Estimator for Geostatistical Inversion and its Application in Hydraulic Tomography journal March 2018
An approach to quantum-computational hydrologic inverse analysis journal May 2018
Randomized Truncated SVD Levenberg-Marquardt Approach to Geothermal Natural State and History Matching text January 2017