Article
 

Hybrid conjugate gradient-Occam algorithms for inversion of multifrequency and multitransmitter EM data

Public Deposited

Downloadable Content

Download PDF
https://ir.library.oregonstate.edu/concern/articles/t722hb49w

Descriptions

Attribute NameValues
Creator
Abstract
  • We describe novel hybrid algorithms for inversion of electromagnetic geophysical data, combining the computational and storage efficiency of a conjugate gradient approach with an Occam scheme for regularization and step-length control. The basic algorithm is based on the observation that iterative solution of the symmetric (Gauss-Newton) normal equations with conjugate gradients effectively generates a sequence of sensitivities for different linear combinations of the data, allowing construction of the Jacobian for a projection of the original full data space. The Occam scheme can then be applied to this projected problem, with the tradeoff parameter chosen by assessing fit to the full data set. For EM geophysical problems with multiple transmitters (either multiple frequencies or source geometries) an extension of the basic hybrid algorithm is possible. In this case multiple forward and adjoint solutions (one each for each transmitter) are required for each step in the iterative normal equation solver, and each corresponds to the sensitivity for a separate linear combination of data. From the perspective of the hybrid approach, with conjugate gradients generating an approximation to the full Jacobian, it is advantageous to save all of the component sensitivities, and use these to solve the projected problem in a larger subspace. We illustrate the algorithms on a simple problem, 2-D magnetotelluric inversion, using synthetic data. Both the basic and modified hybrid schemes produce essentially the same result as an Occam inversion based on a full calculation of the Jacobian, and the modified scheme requires significantly fewer steps (relative to the basic hybrid scheme) to converge to an adequate solution to the normal equations. The algorithms are expected to be useful primarily for 3-D inverse problems for which the computational burden is heavily dominated by solution to the forward and adjoint problems.
  • This is the publisher’s final pdf. The published article is copyrighted by the author and published by John Wiley & Sons, Inc. for the Royal Astronomical Society and can be found at: http://onlinelibrary.wiley.com/journal/10.1111/(ISSN)1365-246X/.
  • Keywords: Magnetotelluric, Inverse theory, Geomagnetic induction
  • Keywords: Magnetotelluric, Inverse theory, Geomagnetic induction
Resource Type
DOI
Date Available
Date Issued
Citation
  • Egbert, G. D. (2012), Hybrid conjugate gradient-Occam algorithms for inversion of multifrequency and multitransmitter EM data. Geophysical Journal International, 190: 255–266. doi: 10.1111/j.1365-246X.2012.05523.x
Journal Title
Journal Volume
  • 190
Journal Issue/Number
  • 1
Rights Statement
Funding Statement (additional comments about funding)
  • This work was supported by DOE grant DE-FG0302ER15318.
Publisher
Language
Replaces

Relationships

Parents:

This work has no parents.

Items

Thumbnail Title Date Uploaded Visibility Actions
file details EgbertGaryCEOASHybridConjugateGradient.pdf 2017-07-21 Public Download