1887
ASEG2013 - 23rd Geophysical Conference
  • ISSN: 2202-0586
  • E-ISSN:

Abstract

The inversion of large-scale magnetic data sets has historically been successfully achieved through integral transforms of the large, dense sensitivity matrix. Two well-known transforms, the discrete Fourier and multi- dimensional wavelet, reduce the required storage and ultimately speed of the inversion by storing only the necessary transform coefficients without losing accuracy. The main drawback of the approaches is the required calculation of the entire dense sensitivity matrix prior to the transform. This process can be much more costly than the inversion itself. We solve the magnetostatic Maxwell’s equation using a finite volume technique on an ocTree-based mesh. The ocTree mesh greatly reduces the time required for the inversion process. When working in the differential equation domain it is not necessary to explicitly form the sensitivity matrix; this decreases the storage requirement of the problem and increases the overall speed of the inversion. The principal mesh is broken up into sub-domain ocTree grids to further enable parallelization of the forward problem. These grids extend the entire domain of the principal mesh to include large regional features that may influence the data. We present the discretization of the equations and verify the accuracy of the modelling both with the principal mesh and with multiple sub-domains. We show a synthetic example and a large field example consisting of over 4 million data and 5 million model cells that was inverted on a desktop computer.

Loading

Article metrics loading...

/content/journals/10.1071/ASEG2013ab209
2013-12-01
2026-01-14
Loading full text...

Full text loading...

References

  1. Ascher, U. M. and Haber, E., 2001, Grid refinement and scaling for distributed parameter estimation problems: Inverse Problems, 17, 517-590.
  2. Bhattacharyya, B. K., 1964, Magnetic anomalies due to prism- shaped bodies with arbitrary magnetizations: Geophysics, 29, 517-531.
  3. Cordell, L., 1992, A scattered equivalent-source method for interpolation and gridding of potential field data in three dimensions: Geophysics, 57, 629-636.
  4. Haber, E., Heldmann, S., and Ascher, U., 2007, Adaptive finite volume method for distributed non-smooth parameter identification: Inverse problems, 23, 1659-1676.
  5. Lelièvre, P. G. and Oldenburg, D. W., 2006, Magnetic forward modelling and inversion for high susceptibility: Geophysical Journal International, 166, 76-90.
  6. Li, Y., and Oldenburg, D. W., 2000, Joint inversion of surface and three-component borehole magnetic data: Geophysics, 65, 540-552.
  7. Li, Y., and Oldenburg, D. W., 2003, Fast inversion of large- scale magnetic data using wavelet transforms and a logarithmic barrier method: Geophysics Journal International, 152, 251-265.
  8. Li, Y., and Oldenburg, D. W., 2010, Rapid construction of equivalent sources using wavelets: Geophysics, 75, L51-L59.
  9. Mareschal, J. C., 1985, Inversion of potential field data in Fourier tranform domain: Geophysics, 50, 685-691.
  10. Phillips, N., Nguyen, T.N.H., and Thomson, V., 2009, QUEST Project: 3D inversion modelling, integration, and visualization of airborne gravity, magnetic, and electromagnetic data, BC, Canada; Geoscience BC, Report 2009-15, 79 p.
  11. Pilkington, M., 1997, 3-D magnetic imaging using conjugate gradients: Geophysics, 62, 1132-1142.
  12. Vogel, C. R., 2002, Computation Methods for Inverse Problems (Frontiers in Applied Mathematics): Society for Industrial Mathematics.
/content/journals/10.1071/ASEG2013ab209
Loading
  • Article Type: Research Article
Keyword(s): Inversion; large-scale; magnetostatic Maxwell’s equations; potential fields
This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error