Upgrade to Pro — share decks privately, control downloads, hide ads and more …

Open source software for simulations and inversions of airborne electromagnetic data

Open source software for simulations and inversions of airborne electromagnetic data

Lindsey J. Heagy, Seogi Kang, Rowan Cockett and Douglas W. Oldenburg

Inversions of airborne EM data are often an iterative process, not only requiring that the researcher be able to explore the impact of changing components such as the choice of regularization functional or model parameterization, but also often requiring that forward simulations be run and fields and fluxes visualized in order to build an understanding of the physical processes governing what we observe in the data. In the hope of facilitating this exploration and promoting reproducibility of geophysical simulations and inversions, we have developed the open source software package, SimPEG. The software has been designed to be modular and extensible with the goal of allowing researchers to interrogate all of the components and to facilitate the exploration of new inversion strategies. We present an overview of the software in its application to airborne EM and demonstrate its use for visualizing fields and fluxes in a forward simulation as well as its flexibility in formulating and solving the inverse problem. We invert a line of airborne TDEM data over a conductive vertical plate using a 1D voxel-inversion, a 2D voxel inversion and a parametric inversion, where all of the forward modelling is done on a 3D grid. The results in this paper can be reproduced by using the provided Jupyter notebooks. The Python software can also be modified to allow users to experiment with parameters and explore the physics of the electromagnetics and intricacies of inversion.

Presented at the 2018 conference on Airborne Electromagnetics (http://www.conferencemanager.dk/AEM2018). The abstract and notebooks used to generate the examples are available at: https://github.com/simpeg-research/heagy-2018-AEM/blob/master/Heagy_etal_2018_aem_workshop.pdf

Lindsey Heagy

June 19, 2018
Tweet

More Decks by Lindsey Heagy

Other Decks in Science

Transcript

  1. Open source software for simulations and inversions of airborne electromagnetic

    data Lindsey Heagy, Seogi Kang, Rowan Cockett, Doug Oldenburg & the SimPEG Team June 19, 2018
  2. building for researchers • Flexibility to experiment ◦ Prioritize organization

    and flexibility • Integration of information ◦ Geologic ◦ Multiple physics • Reproducible and transparent
  3. • Modular: building blocks ◦ organized in a framework ◦

    pieces available to manipulation • Declarative: express intent ◦ write what you mean ◦ looks like the math • Extensible: new research ◦ quantitative communication ◦ built in feedback loops • Open: for the future ◦ reproducible ◦ opportunities for collaboration building for researchers
  4. building for researchers • Modular: building blocks ◦ organized in

    a framework ◦ pieces available to manipulation • Declarative: express intent ◦ write what you mean ◦ looks like the math • Extensible: new research ◦ quantitative communication ◦ built in feedback loops • Open: for the future ◦ reproducible ◦ opportunities for collaboration
  5. a few questions • Understanding responses ◦ How does magnetic

    permeability or chargeability affect my data? • Survey design ◦ Which system should be used to excite / detect the target of interest? ◦ Would adding another receiver orientation improve detection / recovery of the target? • Inversion ◦ Where is a 3D inversion needed? ◦ Can a parametric inversion be used? 17
  6. what is your model? • Subsurface log conductivity • 1D,

    2D, 3D voxel model • Parametric model • … 25 • Need to map to forward simulation mesh • Keep track of derivatives for inversion
  7. 1D inversion • Model: 1D subsurface conductivity at each sounding

    location • Forward simulation mesh: 2D cylindrical mesh • Recent developments: SimPEG EM1D ◦ includes lateral and spatial constraints 26
  8. 2D voxel inversion • 2D model (tensor mesh) • Forward

    simulation mesh: 3D tensor mesh ◦ Map 2D section to 3D space 27
  9. 2D parametric inversion • 6 parameters describing model • Forward

    simulation mesh: 3D tensor mesh ◦ Map parameters to 2D mesh à fill 3D space 29
  10. building for researchers • Modular: building blocks ◦ organized in

    a framework ◦ pieces available to manipulation • Declarative: express intent ◦ write what you mean ◦ looks like the math • Extensible: new research ◦ quantitative communication ◦ built in feedback loops • Open: for the future ◦ reproducible ◦ opportunities for collaboration 1D inversion 2D inversion Parametric inversion
  11. framework • Electromagnetics: ◦ FDEM, TDEM, EM1D (analytic) • Vadose

    zone flow • Gravity • Magnetics • DC / IP • EMIP 31 In one consistent, modular framework
  12. coming soon… 33 • Large scale problems with OcTree meshes

    and domain decomposition • Efficiency improvements including parallelization • Joint and cooperative inversions with petrophysical information • …
  13. Thank you! 35 simpeg.xyz slack.simpeg.xyz github.com/simpeg-research/heagy2018-aem This work is licensed

    under a Creative Commons Attribution 4.0 International License > pip install SimPEG