Journal
PHYSICS OF THE EARTH AND PLANETARY INTERIORS
Volume 150, Issue 1-3, Pages 3-14Publisher
ELSEVIER SCIENCE BV
DOI: 10.1016/j.pepi.2004.08.023
Keywords
magnetotellurics; data-space method; 3D inversion; Occam's inversion
Categories
Ask authors/readers for more resources
A three-dimensional magnetotelluric (MT) minimum structure inversion algorithm has been developed based on a data-space variant of the Occam approach. Computational costs associated with construction and inversion of model-space matrices make a model-space Occam approach to 3D NIT inversion impractical. These difficulties are overcome with a data-space approach, where matrix dimensions depend on the size of the data set, rather than the number of model parameters. With the transformation to data space it becomes feasible to invert modest 3D NIT data sets on a PC. To reduce computational time, a relaxed convergence criterion is used for the iterative forward modeling code used to compute the sensitivity matrix. This allows reduction in computational time by more than 70%, without affecting the inversion results. Numerical experiments with synthetic data show that reasonable fits can be obtained within a small number of iterations, with a few additional iterations required to eliminate unnecessary structure and find the model with minimum norm. (c) 2004 Elsevier B.V. All rights reserved.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available