Aarhus University Seal

An overview of a highly versatile forward and stable inverse algorithm for airborne, ground-based and borehole electromagnetic and electric data

Research output: Contribution to journal/Conference contribution in journal/Contribution to newspaperJournal articleResearchpeer-review


  • Esben Auken
  • Anders Vest Christiansen
  • Casper Kirkegaard
  • ,
  • Gianluca Fiandaca
  • ,
  • Cyril Schamper, Sorbonne Universités, UPMC Univ Paris 06
  • ,
  • Ahmad Ali Behroozmand, Stanford University
  • ,
  • Andrew Binley, University Lancaster
  • ,
  • Emil Nielsen, Technical University of Denmark
  • ,
  • Flemming Effersø, SkyTEM Surveys ApS
  • ,
  • Niels Bøie Christensen
  • Kurt Sørensen, SkyTEM Surveys ApS
  • ,
  • Nikolaj Foged
  • Giulio Vignoli, Geological Survey of Denmark and Greenland

We present an overview of a mature, robust and general algorithm providing a single framework for the inversion of most electromagnetic and electrical data types and instrument geometries. The implementation mainly uses a 1D earth formulation for electromagnetics and magnetic resonance sounding (MRS) responses, while the geoelectric responses are both 1D and 2D and the sheet's response models a 3D conductive sheet in a conductive host with an overburden of varying thickness and resistivity. In all cases, the focus is placed on delivering full system forward modelling across all supported types of data. Our implementation is modular, meaning that the bulk of the algorithm is independent of data type, making it easy to add support for new types. Having implemented forward response routines and file I/O for a given data type provides access to a robust and general inversion engine. This engine includes support for mixed data types, arbitrary model parameter constraints, integration of prior information and calculation of both model parameter sensitivity analysis and depth of investigation. We present a review of our implementation and methodology and show four different examples illustrating the versatility of the algorithm. The first example is a laterally constrained joint inversion (LCI) of surface time domain induced polarisation (TDIP) data and borehole TDIP data. The second example shows a spatially constrained inversion (SCI) of airborne transient electromagnetic (AEM) data. The third example is an inversion and sensitivity analysis of MRS data, where the electrical structure is constrained with AEM data. The fourth example is an inversion of AEM data, where the model is described by a 3D sheet in a layered conductive host.

Original languageEnglish
JournalExploration Geophysics
Pages (from-to)223-235
Number of pages13
Publication statusPublished - 2015

See relations at Aarhus University Citationformats

ID: 91006960