1887
  1. Home
  2. A-Z Publications
  3. Geophysical Prospecting
  4. Volume 63, Issue 2
  5. Article
Volume 63, Issue 2
  • E-ISSN: 1365-2478

Abstract

ABSTRACT

Over the past decade the typical size of airborne electromagnetic data sets has been growing rapidly, along with an emerging need for highly accurate modelling. One‐dimensional approximate inversions or data transform techniques have previously been employed for very large‐scale studies of quasi‐layered settings but these techniques fail to provide the consistent accuracy needed by many modern applications such as aquifer and geological mapping, uranium exploration, oil sands and integrated modelling. In these cases the use of more time‐consuming 1D forward and inverse modelling provide the only acceptable solution that is also computationally feasible.

When target structures are known to be quasi layered and spatially coherent it is beneficial to incorporate this assumption directly into the inversion. This implies inverting multiple soundings at a time in larger constrained problems, which allows for resolving geological layers that are undetectable using simple independent inversions. Ideally, entire surveys should be inverted at a time in huge constrained problems but poor scaling properties of the underlying algorithms typically make this challenging.

Here, we document how we optimized an inversion code for very large‐scale constrained airborne electromagnetic problems. Most importantly, we describe how we solve linear systems using an iterative method that scales linearly with the size of the data set in terms of both solution time and memory consumption. We also describe how we parallelized the core region of the code, in order to obtain almost ideal strong parallel scaling on current 4‐socket shared memory computers. We further show how model parameter uncertainty estimates can be efficiently obtained in linear time and we demonstrate the capabilities of the full implementation by inverting a 3327 line km SkyTEM survey overnight. Performance and scaling properties are discussed based on the timings of the field example and we describe the criteria that must be fulfilled in order to adapt our methodology for similar type problems.

Copyright © 2015 European Association of Geoscientists & Engineers © 2014 European Association of Geoscientists & Engineers
Loading

Article metrics loading...

/content/journals/10.1111/1365-2478.12200
2014-11-14
2024-04-26

  • From This Site

    /content/journals/10.1111/1365-2478.12200
    dcterms_title,dcterms_subject,pub_keyword
    -contentType:Journal -contentType:Contributor -contentType:Concept -contentType:Institution
    10
    5
Loading full text...

Full text loading...

References

  1. AkinJ.E.1982. Application and Implementation of Finite Element Methods . First Edition. Academic Press.
  2. AmestoyP.R., DavisT.A. and DuffI.S.1996. An approximate minimum degree ordering algorithm. SIAM Journal on Matrix Analysis and Applications17, 886–905.
     [Google Scholar]
  3. AmestoyP.R., DuffI.S., L'ExcellentJ.Y. and KosterJ.2002. A fully asynchronous multifrontal solver using distributed dynamic scheduling. SIAM Journal on Matrix Analysis and Applications23, 15–41.
     [Google Scholar]
  4. AukenE., ChristiansenA.V., JacobsenB.H., FogedN. and SørensenK.I.2005. Piecewise 1D Laterally Constrained Inversion of resistivity data. Geophysical Prospecting53, 497–506.
     [Google Scholar]
  5. AukenE., ChristiansenA.V., JacobsenL. and SørensenK.I.2008. A resolution study of buried valleys using laterally constrained inversion of TEM data. Journal of Applied Geophysics65, 10–20.
     [Google Scholar]
  6. AukenE., ChristiansenA.V., KirkegaardC., FiandacaG., SchamperC., BehroozmandA.A., BinleyA., NielsenE., EffersøF., ChristensenN.B., SørensenK.I., FogedN. and VignoliG.2014. An overview of a highly versatile forward and stable inverse algorithm for airborne, ground‐based and borehole electromagnetic and electric data. Exploration Geophysics , 1–13, DOI: 10.1071/EG13097.
  7. BalayS., GroppW.D., McInnesL.C. and SmithB.F.1997. Efficient management of parallelism in object oriented numerical software libraries. Modern Software Tools in Scientific Computing, 163–202.
     [Google Scholar]
  8. BhowmickS., TothB. and RaghavanP.2009. Towards low‐cost, high‐accuracy classifiers for linear solver selection. Computional Science – ICCS 2009. pp. 463–472. Springer.
     [Google Scholar]
  9. BrodieR.2010. Holistic Inversion of Airborne Electromagnetic Data. PhD thesis, The Australian National University, Canberra, Australia.
     [Google Scholar]
  10. BrodieR. and SambridgeM.2006. A holistic approach to inversion of frequency‐domain airborne EM data. Geophysics71, G301–G312.
     [Google Scholar]
  11. ChapmanB., GabrieleJ. and van der PasR.2007. Using OpenMP: Portable Shared Memory Parallel Programming . MIT press.
  12. ChenJ. and RaicheA.1998. Inverting AEM data using a damped eigenparameter method. Exploration Geophysics29, 128–132.
     [Google Scholar]
  13. ChristensenN.B.2002. A generic 1D imaging method for transient electromagnetic data. Geophysics67, 438–447.
     [Google Scholar]
  14. ChristensenN.B., ReidJ.E. and HalkjaerM.2009. Fast approximate inversion of SkyTem airborne electromagnetic data. 1–7.
  15. ChristiansenA.V., AukenE. and ViezzoliA.2011. Quantification of modeling errors in airborne TEM caused by inaccurate system description. Geophysics76, F43–F52.
     [Google Scholar]
  16. CostelloeM.T., WhitakerA., BrodieR., FisherA. and SorensenC.2007. Paterson airborne electromagnetic survey, onshore energy and minerals, Geoscience Australia. ASEG Extended Abstracts 2007, 1.
     [Google Scholar]
  17. CoxL.H., WilsonG.A. and ZhdanovM.S.2010. 3D inversion of airborne electromagnetic data using a moving footprint. Exploration Geophysics41, 250–259.
     [Google Scholar]
  18. CuthillE. and MckeeJ.1969. Reducing the bandwidth of sparse symmetric matrices. ACM ‘69 Proceedings of the 1969 24th National Conference, pp. 157–172.
  19. ErnA., GiovangigliV., KeyesD.E. and SmookeM.D.1994. Towards polyalgorithmic linear system solvers for nonlinear elliptic problems. SIAM Journal on Scientific Computing15, 681–703.
     [Google Scholar]
  20. FarquharsonC.G., OldenburgD.W. and LiY.1999. An approximate inversion algorithm for time‐domain electromagnetic surveys. Journal of Applied Geophysics42, 71–80.
     [Google Scholar]
  21. FarquharsonC.G., OldenburgD.W. and RouthP.S.2003. Simultaneous 1D inversion of loop‐loop electromagnetic data for magnetic susceptibility and electrical conductivity. Geophysics68, 1857–1869.
     [Google Scholar]
  22. FletcherR.1976. Conjugate gradient methods for indefinite systems. Lecture Notes in Mathematics. pp. 73–89. Springer.
     [Google Scholar]
  23. GeorgeT., GuptaA. and SarinV.2008. A recommendation system for preconditioned iterative solvers.ICDM ‘08 Proceedings of the 2008 Eighth IEEE International Conference on Data Mining, pp. 803–808.
  24. GouldN.I.M., ScottJ.A. and HuY.2007. A numerical evaluation of sparse direct solvers for the solution of large sparse symmetric linear systems of equations. ACM Transactions on Mathematical Software33.
     [Google Scholar]
  25. GroppW.D., LuskE. and SkjellumA.1999. Using MPI: Portable Parallel Programming with the Message‐Passing Interface 2nd Edition . MIT press.
  26. GuillemoteauJ., SailhacP. and BehaegelM.2012. Fast approximate 2D inversion of airborne TEM data: Born approximation and empirical approach. Geophysics77, WB89–WB97.
     [Google Scholar]
  27. HestenesM.R. and StiefelE.1952. Methods of conjugate gradients for solving linear systems. Journal of Research of the National Bureau of Standards49, 409–436.
     [Google Scholar]
  28. JørgensenF., ScheerW., ThomsenS., SonnenborgT.O., HinsbyK., WiederholdH. et al. 2012. Transboundary geophysical mapping of geological elements and salinity distribution critical for the assessment of future sea water intrusion in response to sea level rise. Hydrology and Earth System Sciences16, 1845–1962.
     [Google Scholar]
  29. KalvigP.2008. A conceptual geological model of the Voltaian basin based on GEOTEM data. The Voltaian Basin, Ghana.Workshop and Excursion, March 10–17, 2008, Abstract Volume 25–29.
  30. KirkegaardC., SchamperC., ChristiansenA.V., VignoliG. and AukenE.2013. An efficient hybrid inversion scheme combining approximate and full forward solutions for AEM. SAGA AEM, South Africa, 1–2.
     [Google Scholar]
  31. LawrieK.C., TanP., HalasL., AppsH., CullenK., BrodieC.et al. 2010. Using AEM Data as Part of an Integrated Assessment of the Salinity Hazard and Risk to Gunbower State Forest and the River Murray Floodplain in the Gunbower Island‐Barr Creek Reach of the Murray River, SE Australia.ASEG Extended Abstracts 2010, 1.
     [Google Scholar]
  32. LiX.S. and DemmelJ.W.2003. SuperLU_DIST: A scalable distributed‐memory sparse direct solver for unsymmetric linear systems. ACM Transactions on Mathematical Software29, 110–140.
     [Google Scholar]
  33. MacnaeJ., KingA., StolzN., OsmakoffA. and BlahaA.1998. Fast AEM data processing and inversion. Exploration Geophysics29, 163–169.
     [Google Scholar]
  34. MarquartD.1963. An Algorithm for Least Squares Estimation of Nonlinear Parameters. SIAM, Journal of Applied Mathematics11, 431–441.
     [Google Scholar]
  35. MenkeW.1989. Geophysical Data Analysis Discrete Inverse Theory . Academic Press.
  36. Mueller‐PetkeM. and YaramanciU.2010. QT inversion – Comprehensive use of the complete surface NMR data set. Geophysics75, WA199–WA209.
     [Google Scholar]
  37. NachtigalN.M., ReddyS.C. and TrefethenL.N.1992. How fast are nonsymmetric matrix iterations. SIAM Journal on Matrix Analysis and Applications13, 778–795.
     [Google Scholar]
  38. PodgorskiJ.E., KgothangL., NgwisanyiT., PlougC., AukenE. and GreenA.G.2010. Introducing the Okavango Delta, Botswana, Airborne TEM Survey. EAGE Near Surface 2010.
  39. PressW.H., TeukolskyS.A., VetterlingW.T. and FlanneryB.P.2007. Numerical Recipes3rd Edition. Cambridge University Press.
     [Google Scholar]
  40. ReidJ.E. and FullagarP.K.1998. Conductivity‐depth transformation of slingram transient electromagnetic data. Exploration Geophysics, 570–576.
     [Google Scholar]
  41. RückerC., GüntherT. and SpitzerK.2006. Three‐dimensional modelling and inversion of dc resistivity data incorporating topography – I. Modelling. Geophysical Journal International166, 495–505.
     [Google Scholar]
  42. SaadY.1990. SPARSKIT: A basic tool kit for sparse matrix computations. Technical Report RIACS‐90–20, Research Institute for Advanced Computer Science, NASA Ames Research Center.
  43. SaadY.1994. ILUT: A dual threshold incomplete LU factorization. Numerical Linear Algebra with Applications1, 387–402.
     [Google Scholar]
  44. SaadY.2003. Iterative Methods for Sparse Linear Systems (2nd edition). SIAM.
     [Google Scholar]
  45. SaadY. and SchultzM.H.1986. GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM Journal on Scientific and Statistical Computing7, 856–869.
     [Google Scholar]
  46. SaadY. and Van der VorstH.A.2000. Iterative solution of linear systems in the 20th century. Journal of Computational and Applied Mathematics123, 1–33.
     [Google Scholar]
  47. SantosF.A.M.2004. 1D laterally constrained inversion of EM34 profiling data. Journal of Applied Geophysics56, 123–134.
     [Google Scholar]
  48. SattelD.2005. Inverting airborne electromagnetic (AEM) data with Zohdy's method. Geophysics70, G77–G85.
     [Google Scholar]
  49. SchenkO. and GärtnerK.2004. Solving unsymmetric sparse systems of linear equations with PARDISO. Future Generation Computer Systems20, 475–487.
     [Google Scholar]
  50. SørensenK.I. and AukenE.2004. SkyTEM – A new high‐resolution helicopter transient electromagnetic system. Exploration Geophysics35, 191–199.
     [Google Scholar]
  51. StreichR., BeckenM. and RitterO.2011. 2.5D controlled‐source EM modeling with general 3D source geometries. Geophysics76, F387–F393.
     [Google Scholar]
  52. TarantolaA. and ValetteB.1982. Generalized nonlinear inverse problems solved using a least squares criterion. Reviews of Geophysics and Space Physics20, 219–232.
     [Google Scholar]
  53. TartarasE. and BeamishD.2005. Laterally constrained inversion of fixed‐wing frequency‐domain AEM data. 12th European Meeting of Environmental and Near Surface Geophysics.
  54. ToledoS., ChenD. and RothkinV.2001. Taucs – A library of sparse linear solvers: http://www.tau.ac.il/~stoledo/taucs/.
  55. ValléeM.A. and SmithR.S.2009. Inversion of airborne time‐domain electromagnetic data to a 1D structure using lateral constraints. Near Surface Geophysics7, 63–71.
     [Google Scholar]
  56. Van der VorstH.A.1992. Bi‐CGSTAB: A fast and smoothly converging variant of Bi‐CG for the solution of nonsymmetric linear systems. SIAM Journal on Scientific and Statistical Computing13, 631–644.
     [Google Scholar]
  57. ViezzoliA., ChristiansenA.V., AukenE. and SørensenK.I.2008. Quasi‐3D modeling of airborne TEM data by Spatially Constrained Inversion. Geophysics73, F105–F113.
     [Google Scholar]
  58. WilsonG.A., RaicheA. and SugengF.2006. 2.5D inversion of airborne electromagnetic data. Exploration Geophysics, 363–371.
     [Google Scholar]
  59. WolfgramP., SattelD. and ChristensenN.B.2003. Approximate 2D inversion of AEM data. Exploration Geophysics34, 29–33.
     [Google Scholar]
  60. WuX.2003. A 3D finite‐element algorithm for DC resistivity modelling using the shifted incomplete Cholesky conjugate gradient method. Geophysical Journal International154, 947–956.
     [Google Scholar]
  61. YangD. and OldenburgD.W.2012. Three‐dimensional inversion of airborne time‐domain electromagnetic data with applications to a porphyry deposit. Geophysics77, B23–B34.
     [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.1111/1365-2478.12200
Loading
A parallel, scalable and memory efficient inversion code for very large‐scale airborne electromagnetics surveys
Geophysical Prospecting 63, 495 (2015); https://doi.org/10.1111/1365-2478.12200
/content/journals/10.1111/1365-2478.12200
/content/journals/10.1111/1365-2478.12200
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): Inversion code; Modelling

Most Cited This Month Most Cited RSS feed

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