DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Scalable subsurface inverse modeling of huge data sets with an application to tracer concentration breakthrough data from magnetic resonance imaging

Abstract

When characterizing subsurface properties is crucial for reliable and cost-effective groundwater supply management and contaminant remediation. With recent advances in sensor technology, large volumes of hydro-geophysical and geochemical data can be obtained to achieve high-resolution images of subsurface properties. However, characterization with such a large amount of information requires prohibitive computational costs associated with “big data” processing and numerous large-scale numerical simulations. To tackle such difficulties, the Principal Component Geostatistical Approach (PCGA) has been proposed as a “Jacobian-free” inversion method that requires much smaller forward simulation runs for each iteration than the number of unknown parameters and measurements needed in the traditional inversion methods. PCGA can be conveniently linked to any multi-physics simulation software with independent parallel executions. In our paper, we extend PCGA to handle a large number of measurements (e.g. 106 or more) by constructing a fast preconditioner whose computational cost scales linearly with the data size. For illustration, we characterize the heterogeneous hydraulic conductivity (K) distribution in a laboratory-scale 3-D sand box using about 6 million transient tracer concentration measurements obtained using magnetic resonance imaging. Since each individual observation has little information on the K distribution, the data was compressed by the zero-th temporal moment of breakthroughmore » curves, which is equivalent to the mean travel time under the experimental setting. Moreover, only about 2,000 forward simulations in total were required to obtain the best estimate with corresponding estimation uncertainty, and the estimated K field captured key patterns of the original packing design, showing the efficiency and effectiveness of the proposed method. This article is protected by copyright. All rights reserved.« less

Authors:
 [1];  [2];  [1];  [3];  [4]
  1. Stanford Univ., Stanford, CA (United States)
  2. Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
  3. Univ. of Texas, Austin, TX (United States)
  4. Univ. of Illinois, Urbana-Champaign, IL (United States)
Publication Date:
Research Org.:
Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Energy Frontier Research Centers (EFRC) (United States). Center for Frontiers of Subsurface Energy Security (CFSES)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA); USDOE Office of Science (SC), Basic Energy Sciences (BES)
OSTI Identifier:
1257809
Report Number(s):
SAND2016-1277J
Journal ID: ISSN 0043-1397; 619422
Grant/Contract Number:  
AC04-94AL85000; SC0001114
Resource Type:
Accepted Manuscript
Journal Name:
Water Resources Research
Additional Journal Information:
Journal Name: Water Resources Research; Journal ID: ISSN 0043-1397
Publisher:
American Geophysical Union (AGU)
Country of Publication:
United States
Language:
English
Subject:
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND TECHNOLOGY; 58 GEOSCIENCES

Citation Formats

Lee, Jonghyun, Yoon, Hongkyu, Kitanidis, Peter K., Werth, Charles J., and Valocchi, Albert J. Scalable subsurface inverse modeling of huge data sets with an application to tracer concentration breakthrough data from magnetic resonance imaging. United States: N. p., 2016. Web. doi:10.1002/2015wr018483.
Lee, Jonghyun, Yoon, Hongkyu, Kitanidis, Peter K., Werth, Charles J., & Valocchi, Albert J. Scalable subsurface inverse modeling of huge data sets with an application to tracer concentration breakthrough data from magnetic resonance imaging. United States. https://doi.org/10.1002/2015wr018483
Lee, Jonghyun, Yoon, Hongkyu, Kitanidis, Peter K., Werth, Charles J., and Valocchi, Albert J. Thu . "Scalable subsurface inverse modeling of huge data sets with an application to tracer concentration breakthrough data from magnetic resonance imaging". United States. https://doi.org/10.1002/2015wr018483. https://www.osti.gov/servlets/purl/1257809.
@article{osti_1257809,
title = {Scalable subsurface inverse modeling of huge data sets with an application to tracer concentration breakthrough data from magnetic resonance imaging},
author = {Lee, Jonghyun and Yoon, Hongkyu and Kitanidis, Peter K. and Werth, Charles J. and Valocchi, Albert J.},
abstractNote = {When characterizing subsurface properties is crucial for reliable and cost-effective groundwater supply management and contaminant remediation. With recent advances in sensor technology, large volumes of hydro-geophysical and geochemical data can be obtained to achieve high-resolution images of subsurface properties. However, characterization with such a large amount of information requires prohibitive computational costs associated with “big data” processing and numerous large-scale numerical simulations. To tackle such difficulties, the Principal Component Geostatistical Approach (PCGA) has been proposed as a “Jacobian-free” inversion method that requires much smaller forward simulation runs for each iteration than the number of unknown parameters and measurements needed in the traditional inversion methods. PCGA can be conveniently linked to any multi-physics simulation software with independent parallel executions. In our paper, we extend PCGA to handle a large number of measurements (e.g. 106 or more) by constructing a fast preconditioner whose computational cost scales linearly with the data size. For illustration, we characterize the heterogeneous hydraulic conductivity (K) distribution in a laboratory-scale 3-D sand box using about 6 million transient tracer concentration measurements obtained using magnetic resonance imaging. Since each individual observation has little information on the K distribution, the data was compressed by the zero-th temporal moment of breakthrough curves, which is equivalent to the mean travel time under the experimental setting. Moreover, only about 2,000 forward simulations in total were required to obtain the best estimate with corresponding estimation uncertainty, and the estimated K field captured key patterns of the original packing design, showing the efficiency and effectiveness of the proposed method. This article is protected by copyright. All rights reserved.},
doi = {10.1002/2015wr018483},
journal = {Water Resources Research},
number = ,
volume = ,
place = {United States},
year = {Thu Jun 09 00:00:00 EDT 2016},
month = {Thu Jun 09 00:00:00 EDT 2016}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 26 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Inverse problem in hydrogeology
journal, February 2005

  • Carrera, Jes�s; Alcolea, Andr�s; Medina, Agust�n
  • Hydrogeology Journal, Vol. 13, Issue 1
  • DOI: 10.1007/s10040-004-0404-7

A modified Levenberg–Marquardt algorithm for quasi-linear geostatistical inversing
journal, July 2004


Hybrid Krylov Methods for Nonlinear Systems of Equations
journal, May 1990

  • Brown, Peter N.; Saad, Youcef
  • SIAM Journal on Scientific and Statistical Computing, Vol. 11, Issue 3
  • DOI: 10.1137/0911026

Effect of radial flow on deviations from local equilibrium during sorbing solute transport through homogeneous soils
journal, November 1986


Quasi-Linear Geostatistical Theory for Inversing
journal, October 1995

  • Kitanidis, Peter K.
  • Water Resources Research, Vol. 31, Issue 10
  • DOI: 10.1029/95WR01945

Principal Component Geostatistical Approach for large-dimensional inverse problems
journal, July 2014


Finite elements for elliptic problems with stochastic coefficients
journal, February 2005

  • Frauenfelder, Philipp; Schwab, Christoph; Todor, Radu Alexandru
  • Computer Methods in Applied Mechanics and Engineering, Vol. 194, Issue 2-5
  • DOI: 10.1016/j.cma.2004.04.008

Bayesian and Geostatistical Approaches to Inverse Problems
book, October 2010


Randomized algorithms for generalized Hermitian eigenvalue problems with application to computing Karhunen-Loève expansion: RANDOMIZED ALGORITHMS FOR GHEP
journal, November 2015

  • Saibaba, Arvind K.; Lee, Jonghyun; Kitanidis, Peter K.
  • Numerical Linear Algebra with Applications, Vol. 23, Issue 2
  • DOI: 10.1002/nla.2026

GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems
journal, July 1986

  • Saad, Youcef; Schultz, Martin H.
  • SIAM Journal on Scientific and Statistical Computing, Vol. 7, Issue 3
  • DOI: 10.1137/0907058

Analysis of hydraulic tomography using temporal moments of drawdown recovery data: HYDRAULIC TOMOGRAPHY
journal, February 2006

  • Zhu, Junfeng; Yeh, Tian-Chyi J.
  • Water Resources Research, Vol. 42, Issue 2
  • DOI: 10.1029/2005WR004309

A Reassessment of the Groundwater Inverse Problem
journal, May 1996

  • McLaughlin, Dennis; Townley, Lloyd R.
  • Water Resources Research, Vol. 32, Issue 5
  • DOI: 10.1029/96WR00160

New Trends in Fiber-Optic Chemical and Biological Sensors
journal, October 2008


Preconditioned all-at-once methods for large, sparse parameter estimation problems
journal, November 2001


Robust inversion, dimensionality reduction, and randomized sampling
journal, June 2012

  • Aravkin, Aleksandr; Friedlander, Michael P.; Herrmann, Felix J.
  • Mathematical Programming, Vol. 134, Issue 1
  • DOI: 10.1007/s10107-012-0571-6

Temporal Moment-Generating Equations: Modeling Transport and Mass Transfer in Heterogeneous Aquifers
journal, August 1995

  • Harvey, Charles F.; Gorelick, Steven M.
  • Water Resources Research, Vol. 31, Issue 8
  • DOI: 10.1029/95WR01231

Solution of Sparse Indefinite Systems of Linear Equations
journal, September 1975

  • Paige, C. C.; Saunders, M. A.
  • SIAM Journal on Numerical Analysis, Vol. 12, Issue 4
  • DOI: 10.1137/0712047

Sensitivity of temporal moments calculated by the adjoint-state method and joint inversing of head and tracer data
journal, October 2000


Hydraulic tomography using temporal moments of drawdown recovery data: A laboratory sandbox study: LAB VALIDATION OF HYDRAULIC TOMOGRAPHY
journal, January 2009

  • Yin, Danting; Illman, Walter A.
  • Water Resources Research, Vol. 45, Issue 1
  • DOI: 10.1029/2007WR006623

Characterization of NAPL Source Zone Architecture and Dissolution Kinetics in Heterogeneous Porous Media Using Magnetic Resonance Imaging
journal, May 2007

  • Zhang, Changyong; Werth, Charles J.; Webb, Andrew G.
  • Environmental Science & Technology, Vol. 41, Issue 10
  • DOI: 10.1021/es061675q

Numerical solution of saddle point problems
journal, April 2005


Imaging geochemical heterogeneities using inverse reactive transport modeling: An example relevant for characterizing arsenic mobilization and distribution
journal, February 2016


Large-scale sandbox experiment on longitudinal effective dispersion in heterogeneous porous media: LONGITUDINAL EFFECTIVE DISPERSION
journal, December 2004

  • Jose, Surabhin C.; Rahman, M. Arifur; Cirpka, Olaf A.
  • Water Resources Research, Vol. 40, Issue 12
  • DOI: 10.1029/2004WR003363

Fast Algorithms for Bayesian Uncertainty Quantification in Large-Scale Linear Inverse Problems Based on Low-Rank Partial Hessian Approximations
journal, January 2011

  • Flath, H. P.; Wilcox, L. C.; Akçelik, V.
  • SIAM Journal on Scientific Computing, Vol. 33, Issue 1
  • DOI: 10.1137/090780717

Fast full-wavefield seismic inversion using encoded sources
journal, November 2009

  • Krebs, Jerome R.; Anderson, John E.; Hinkley, David
  • GEOPHYSICS, Vol. 74, Issue 6
  • DOI: 10.1190/1.3230502

Geostatistical inference of hydraulic conductivity and dispersivities from hydraulic heads and tracer data: GEOSTATISTICAL INFERENCE
journal, August 2006

  • Nowak, Wolfgang; Cirpka, Olaf A.
  • Water Resources Research, Vol. 42, Issue 8
  • DOI: 10.1029/2005WR004832

On Integrating Groundwater Transport Models with Wireless Sensor Networks
journal, August 2010


An Effective Method for Parameter Estimation with PDE Constraints with Multiple Right-Hand Sides
journal, January 2012

  • Haber, Eldad; Chung, Matthias; Herrmann, Felix
  • SIAM Journal on Optimization, Vol. 22, Issue 3
  • DOI: 10.1137/11081126X

Inverse problems: A Bayesian perspective
journal, May 2010


Likelihood-informed dimension reduction for nonlinear inverse problems
journal, October 2014


Orthonormal residuals in geostatistics: Model criticism and parameter estimation
journal, June 1991


Acquisition using simultaneous sources
journal, July 2008

  • Hampson, Gary; Stefani, Joe; Herkenhoff, Fred
  • The Leading Edge, Vol. 27, Issue 7
  • DOI: 10.1190/1.2954034

Hydraulic Tomography: Continuity and Discontinuity of High- K and Low- K Zones : D. L. Hochstetler
journal, June 2015

  • Hochstetler, David L.; Barrash, Warren; Leven, Carsten
  • Groundwater, Vol. 54, Issue 2
  • DOI: 10.1111/gwat.12344

Acquisition Using Simultaneous Sources
conference, January 2007

  • Stefani, J.; Hampson, G.; F. Herkenhoff, E.
  • 69th EAGE Conference and Exhibition incorporating SPE EUROPEC 2007
  • DOI: 10.3997/2214-4609.201401462

A Newton-CG method for large-scale three-dimensional elastic full-waveform seismic inversion
journal, May 2008


Fast computation of uncertainty quantification measures in the geostatistical approach to solve inverse problems
journal, August 2015


Statistical inverse problems: Discretization, model reduction and inverse crimes
journal, January 2007

  • Kaipio, Jari; Somersalo, Erkki
  • Journal of Computational and Applied Mathematics, Vol. 198, Issue 2
  • DOI: 10.1016/j.cam.2005.09.027

Geosensing for Developing Sustainable Responses to Environmental Hazards Underground
conference, March 2014


Numerical methods for large eigenvalue problems
journal, January 2002


Acquisition using simultaneous sources
conference, December 2008

  • Hampson, Gary; Stefani, Joe; Herkenhoff, Fred
  • SEG Technical Program Expanded Abstracts 2008
  • DOI: 10.1190/1.3063930

Works referencing / citing this record:

Randomized Truncated SVD Levenberg‐Marquardt Approach to Geothermal Natural State and History Matching
journal, March 2018

  • Bjarkason, Elvar K.; Maclaren, Oliver J.; O'Sullivan, John P.
  • Water Resources Research, Vol. 54, Issue 3
  • DOI: 10.1002/2017wr021870

A review of geochemical–mechanical impacts in geological carbon storage reservoirs
journal, May 2019

  • Akono, Ange‐Therese; Druhan, Jennifer L.; Dávila, Gabriela
  • Greenhouse Gases: Science and Technology, Vol. 9, Issue 3
  • DOI: 10.1002/ghg.1870

Randomized Truncated SVD Levenberg-Marquardt Approach to Geothermal Natural State and History Matching
text, January 2017