A computationally efficient parallel Levenberg‐Marquardt algorithm for highly parameterized inverse model analyses

Inverse modeling seeks model parameters given a set of observations. However, for practical problems because the number of measurements is often large and the model parameters are also numerous, conventional methods for inverse modeling can be computationally expensive. We have developed a new, comp...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:Water resources research Jg. 52; H. 9; S. 6948 - 6977
Hauptverfasser: Lin, Youzuo, O'Malley, Daniel, Vesselinov, Velimir V.
Format: Journal Article
Sprache:Englisch
Veröffentlicht: Washington John Wiley & Sons, Inc 01.09.2016
American Geophysical Union (AGU)
Schlagworte:
Algorithms
Computational efficiency
Computer applications
Damping
dimensionality reduction
Earth Sciences
Efficiency
Environment models
Fields
Frameworks
Heterogeneity
Hydraulic conductivity
hydraulic inverse modeling
Krylov subspace approximation
Levenberg‐Marquardt method
LSQR iterative method
Mathematical models
MATHEMATICS AND COMPUTING
Methods
Modelling
Parameterization
Parameters
Problems
Software
subspace recycling
Three dimensional models
Transmissivity
ISSN:0043-1397, 1944-7973
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:Inverse modeling seeks model parameters given a set of observations. However, for practical problems because the number of measurements is often large and the model parameters are also numerous, conventional methods for inverse modeling can be computationally expensive. We have developed a new, computationally efficient parallel Levenberg‐Marquardt method for solving inverse modeling problems with a highly parameterized model space. Levenberg‐Marquardt methods require the solution of a linear system of equations which can be prohibitively expensive to compute for moderate to large‐scale problems. Our novel method projects the original linear problem down to a Krylov subspace such that the dimensionality of the problem can be significantly reduced. Furthermore, we store the Krylov subspace computed when using the first damping parameter and recycle the subspace for the subsequent damping parameters. The efficiency of our new inverse modeling algorithm is significantly improved using these computational techniques. We apply this new inverse modeling method to invert for random transmissivity fields in 2‐D and a random hydraulic conductivity field in 3‐D. Our algorithm is fast enough to solve for the distributed model parameters (transmissivity) in the model domain. The algorithm is coded in Julia and implemented in the MADS computational framework (http://mads.lanl.gov). By comparing with Levenberg‐Marquardt methods using standard linear inversion techniques such as QR or SVD methods, our Levenberg‐Marquardt method yields a speed‐up ratio on the order of ∼101 to ∼102 in a multicore computational environment. Therefore, our new inverse modeling method is a powerful tool for characterizing subsurface heterogeneity for moderate to large‐scale problems. Key Points We generate a Krylov subspace and obtain a much smaller approximated problem by projecting the original problem down to the subspace We employ both coarse and fine‐grained parallelism to our LM method to parallelize the implementation of the LM algorithm in two levels We employ a subspace recycling technique to take the advantage of the solution space previously generated
Bibliographie:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
content type line 23
USDOE
AC52-06NA25396
LANL EP Program
LA-UR-16-22377
ISSN:0043-1397
1944-7973
DOI:10.1002/2016WR019028