Airborne electromagnetic modelling options and their consequences in target definition

Given the range of geological conditions under which airborne EM surveys are conducted, there is an expectation that the 2D and 3D methods used to extract models that are geologically meaningful would be favoured over 1D inversion and transforms. We do after all deal with an Earth that constantly undergoes, faulting, intrusions, and erosive processes that yield a subsurface morphology, which is, for most parts, dissimilar to a horizontal layered earth. We analyse data from a survey collected in the Musgrave province, South Australia. It is of particular interest since it has been used for mineral prospecting and for a regional hydro-geological assessment. The survey comprises abrupt lateral variations, more-subtle lateral continuous sedimentary sequences and filled palaeovalleys. As consequence, we deal with several geophysical targets of contrasting conductivities, varying geometries and at different depths. We invert the observations by using several algorithms characterised by the different dimensionality of the forward operator. Inversion of airborne EM data is known to be an ill-posed problem. We can generate a variety of models that numerically adequately fit the measured data, which makes the solution non-unique. The application of different deterministic inversion codes or transforms to the same dataset can give dissimilar results, as shown in this paper. This ambiguity suggests the choice of processes and algorithms used to interpret AEM data cannot be resolved as a matter of personal choice and preference. The degree to which models generated by a 1D algorithm replicate/or not measured data, can be an indicator of the data’s dimensionality, which perse does not imply that data that can be fitted with a 1D model cannot be multidimensional. On the other hand, it is crucial that codes that can generate 2D and 3D models do reproduce the measured data in order for them to be considered as a plausible solution. In the absence of ancillary information, it could be argued that the simplest model with the simplest physics might be preferred. Given the range of geological conditions under which airborne EM surveys are conducted, there is an expectation that 2D and 3D methods used to extract models of geological significance would be favoured over 1D inversion and transforms. We analyse data from the Musgrave province, South Australia, used for mineral and for hydro-geological prospecting.

[1]  Robert G. Ellis,et al.  Inversion of airborne electromagnetic data , 1998 .

[2]  Esben Auken,et al.  Sharp spatially constrained inversion with applications to transient electromagnetic data , 2015 .

[3]  Michael S. Zhdanov,et al.  Focusing geophysical inversion images , 1999 .

[4]  Andrew Binley,et al.  An overview of a highly versatile forward and stable inverse algorithm for airborne, ground-based and borehole electromagnetic and electric data , 2015 .

[5]  A. Christiansen,et al.  Quasi-3D modeling of airborne TEM data by spatially constrained inversion , 2008 .

[6]  James Macnae,et al.  Fast AEM data processing and inversion , 1998 .

[7]  Michael S. Zhdanov,et al.  3D inversion of airborne electromagnetic data , 2012 .

[8]  Estimating Noise Levels in AEM Data , 2003 .

[9]  J. Guillemoteau,et al.  Regularization strategy for the layered inversion of airborne transient electromagnetic data: application to in‐loop data acquired over the basin of Franceville (Gabon) , 2011, 1203.6539.

[10]  D. Sattel,et al.  An Overview of Helicopter Time-Domain EM Systems , 2009 .

[11]  Esben Auken,et al.  Layered and laterally constrained 2D inversion of resistivity data , 2004 .

[12]  Bill Russell-Cargill Vessel based solutions for the acquisition of geotechnical data in subsea site investigations , 2010 .

[13]  Niels B. Christensen,et al.  A generic 1-D imaging method for transient electromagnetic data , 2002 .

[14]  D. Oldenburg,et al.  Estimating depth of investigation in DC resistivity and IP surveys , 1999 .

[15]  Breaks in lithology: Interpretation problems when handling 2D structures with a 1D approximation , 2010 .

[16]  Michael S. Zhdanov,et al.  Geophysical Electromagnetic Theory and Methods , 2009 .

[17]  Esben Auken,et al.  A Global Measure for Depth of Investigation , 2011 .

[18]  Giorgio Cassiani,et al.  Statistical multioffset phase analysis for surface-wave processing in laterally varying media , 2011 .

[19]  A. Peter Annan,et al.  A comparison of data from airborne, semi‐airborne, and ground electromagnetic systems , 2001 .

[20]  G. Ross The Geology of South Australia: Volume 1. The Precambrian , 1997 .

[21]  Niels B. Christensen,et al.  Resolution analyses for selecting an appropriate airborne electromagnetic (AEM) system , 2012 .

[22]  I. Roach The Frome airborne electromagnetic survey, South Australia , 2012 .

[23]  J. Macnae,et al.  Detection and correction of SPM effects in airborne EM surveys , 2013 .

[24]  A. Christiansen,et al.  An integrated processing scheme for high-resolution airborne electromagnetic surveys, the SkyTEM system , 2009 .

[25]  Michael S. Zhdanov,et al.  3D inversion of airborne electromagnetic data using a moving footprint , 2010 .

[26]  Pascal Sailhac,et al.  Fast approximate 2D inversion of airborne TEM data: Born approximation and empirical approach , 2012 .

[27]  Michael S. Zhdanov,et al.  Sharp boundary inversion in crosswell travel-time tomography , 2006 .

[28]  H. Johnson,et al.  A comparison of 'traditional' and multimedia information systems development practices , 2003, Inf. Softw. Technol..

[29]  J. G. The Geology of South Australia , 1919, Nature.

[30]  K. Kubik,et al.  Compact gravity inversion , 1983 .