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 »
- Authors:
-
- Stanford Univ., Stanford, CA (United States)
- Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
- Univ. of Texas, Austin, TX (United States)
- 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}
}
Web of Science
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
A modified Levenberg–Marquardt algorithm for quasi-linear geostatistical inversing
journal, July 2004
- Nowak, Wolfgang; Cirpka, Olaf A.
- Advances in Water Resources, Vol. 27, Issue 7
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
Effect of radial flow on deviations from local equilibrium during sorbing solute transport through homogeneous soils
journal, November 1986
- Valocchi, Albert J.
- Water Resources Research, Vol. 22, Issue 12
Quasi-Linear Geostatistical Theory for Inversing
journal, October 1995
- Kitanidis, Peter K.
- Water Resources Research, Vol. 31, Issue 10
Principal Component Geostatistical Approach for large-dimensional inverse problems
journal, July 2014
- Kitanidis, P. K.; Lee, J.
- Water Resources Research, Vol. 50, Issue 7
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
Bayesian and Geostatistical Approaches to Inverse Problems
book, October 2010
- Kitanidis, P. K.
- Large‐Scale Inverse Problems and Quantification of Uncertainty
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
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
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
A Reassessment of the Groundwater Inverse Problem
journal, May 1996
- McLaughlin, Dennis; Townley, Lloyd R.
- Water Resources Research, Vol. 32, Issue 5
New Trends in Fiber-Optic Chemical and Biological Sensors
journal, October 2008
- Orellana, Guillermo; Haigh, David
- Current Analytical Chemistry, Vol. 4, Issue 4
Highly parameterized inverse estimation of hydraulic conductivity and porosity in a three-dimensional, heterogeneous transport experiment: HIGHLY PARAMETERIZED INVERSE ESTIMATION IN 3-D
journal, October 2012
- Yoon, Hongkyu; McKenna, Sean A.
- Water Resources Research, Vol. 48, Issue 10
Coupled simulation of DNAPL infiltration and dissolution in three-dimensional heterogeneous domains: Process model validation: 3-D Process Modeling of Dnapl Infiltration and Dissolution
journal, October 2013
- Kokkinaki, A.; O'Carroll, D. M.; Werth, C. J.
- Water Resources Research, Vol. 49, Issue 10
Preconditioned all-at-once methods for large, sparse parameter estimation problems
journal, November 2001
- Haber, E.; Ascher, U. M.
- Inverse Problems, Vol. 17, Issue 6
Robust inversion, dimensionality reduction, and randomized sampling
journal, June 2012
- Aravkin, Aleksandr; Friedlander, Michael P.; Herrmann, Felix J.
- Mathematical Programming, Vol. 134, Issue 1
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
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
Sensitivity of temporal moments calculated by the adjoint-state method and joint inversing of head and tracer data
journal, October 2000
- Cirpka, Olaf A.; Kitanidis, Peter K.
- Advances in Water Resources, Vol. 24, Issue 1
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
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
Numerical solution of saddle point problems
journal, April 2005
- Benzi, Michele; Golub, Gene H.; Liesen, Jörg
- Acta Numerica, Vol. 14
Imaging geochemical heterogeneities using inverse reactive transport modeling: An example relevant for characterizing arsenic mobilization and distribution
journal, February 2016
- Fakhreddine, Sarah; Lee, Jonghyun; Kitanidis, Peter K.
- Advances in Water Resources, Vol. 88
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
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
Fast full-wavefield seismic inversion using encoded sources
journal, November 2009
- Krebs, Jerome R.; Anderson, John E.; Hinkley, David
- GEOPHYSICS, Vol. 74, Issue 6
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
On Integrating Groundwater Transport Models with Wireless Sensor Networks
journal, August 2010
- Barnhart, Kevin; Urteaga, Iñigo; Han, Qi
- Ground Water, Vol. 48, Issue 5
Large-scale hydraulic tomography and joint inversion of head and tracer data using the Principal Component Geostatistical Approach (PCGA)
journal, July 2014
- Lee, J.; Kitanidis, P. K.
- Water Resources Research, Vol. 50, Issue 7
Numerical simulation of water flow in three dimensional heterogeneous porous media observed in a magnetic resonance imaging experiment: NUMERICAL SIMULATION OF WATER FLOW IN 3D
journal, June 2008
- Yoon, Hongkyu; Zhang, Changyong; Werth, Charles J.
- Water Resources Research, Vol. 44, Issue 6
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
Likelihood-informed dimension reduction for nonlinear inverse problems
journal, October 2014
- Cui, T.; Martin, J.; Marzouk, Y. M.
- Inverse Problems, Vol. 30, Issue 11
Orthonormal residuals in geostatistics: Model criticism and parameter estimation
journal, June 1991
- Kitanidis, Peter K.
- Mathematical Geology, Vol. 23, Issue 5
Acquisition using simultaneous sources
journal, July 2008
- Hampson, Gary; Stefani, Joe; Herkenhoff, Fred
- The Leading Edge, Vol. 27, Issue 7
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
Acquisition Using Simultaneous Sources
conference, January 2007
- Stefani, J.; Hampson, G.; F. Herkenhoff, E.
- 69th EAGE Conference and Exhibition incorporating SPE EUROPEC 2007
A Newton-CG method for large-scale three-dimensional elastic full-waveform seismic inversion
journal, May 2008
- Epanomeritakis, I.; Akçelik, V.; Ghattas, O.
- Inverse Problems, Vol. 24, Issue 3
Fast computation of uncertainty quantification measures in the geostatistical approach to solve inverse problems
journal, August 2015
- Saibaba, Arvind K.; Kitanidis, Peter K.
- Advances in Water Resources, Vol. 82
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
Geosensing for Developing Sustainable Responses to Environmental Hazards Underground
conference, March 2014
- Pamukcu, Sibel; Ghazanfari, Ehsan
- Geo-Congress 2014 Keynote Lectures
Numerical methods for large eigenvalue problems
journal, January 2002
- Sorensen, Danny C.
- Acta Numerica, Vol. 11
Acquisition using simultaneous sources
conference, December 2008
- Hampson, Gary; Stefani, Joe; Herkenhoff, Fred
- SEG Technical Program Expanded Abstracts 2008
Fast computation of uncertainty quantification measures in the geostatistical approach to solve inverse problems
text, January 2014
- Saibaba, Arvind K.; Kitanidis, Peter K.
- arXiv
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
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
Randomized Truncated SVD Levenberg-Marquardt Approach to Geothermal Natural State and History Matching
text, January 2017
- Bjarkason, Elvar K.; Maclaren, Oliver J.; O'Sullivan, John P.
- arXiv