Modeling electromagnetics on cylindrical meshes with applications to steel-cased wells

Overview

This repository contains the notebooks used to generate the examples shown in "Modeling electromagnetics on cylindrical meshes with applications to steel-cased wells" by Lindsey J. Heagy and Douglas W. Oldenburg and Chapter 3 of the thesis Electromagnetic imaging of subsurface injections. There are two directories of notebooks

notebooks: contains the notebooks used to generate figures in the paper and thesis
- 1_TDEM_Comparison_w_Commer:
  - Figures 3 and 4 in Heagy & Oldenburg (2018)
  - Figures 3.3 and 3.4 in the thesis
- 2_DC_Kaufman_infinite_well:
  - Figure 5 in Heagy & Oldenburg (2018)
  - Figure 3.5 in the thesis
- 3_DC_Kaufman_finite_well:
  - Figure 6 in Heagy & Oldenburg (2018)
  - Figure 3.6 in the thesis
- 4_TDEM_current_plots:
  - Figure 7 in Heagy & Oldenburg (2018)
- 5_FDEM_following_Augustin_Fig3:
  - Figures 8 - 10 in Heagy & Oldenburg (2018)
  - Figures 3.7 - 3.9 in the thesis
- 6_FDEM_Permeability:
  - Figures 11 and 12 in Heagy & Oldenburg (2018)
  - Figures 3.10 and 3.11 in the thesis
- 7_TDEM_Permeability:
  - Figures 13 and 14 in Heagy & Oldenburg (2018)
  - Figures 3.12 and 3.12 in the thesis
appendix_notebooks: supplementary notebooks that contain notebooks examining the impact of formulation choice for the discretization of Maxwell's equations and additional tests for the 2D and 3D cyl codes
- compare_2D_and_3D_with_target: compares results of 2D and 3D problems for a cylindrically symmetric source and target
- compare_3D_FDEM_horizontal_electric_dipole: compares the 3D cyl code with an analytic for a horizontal electric dipole (FDEM) that is offset from the axis of symmetry
- compare_3D_cyl_and_tensor_sources_downhole: this example compares results from a 3D tensor simulation and a 3D cyl simulation in a grounded-source FDEM experiment
- compare_TDEM_horizontal_electric_dipole: compares the 3D cyl code with an analytic for a horizontal electric dipole (TDEM) that is offset from the axis of symmetry
- cyl3D_DC_horizontal_electric_dipole: compares the 3D cyl code with an analytic for a horizontal electric dipole (DC) that is offset from the axis of symmetry
- cyl3D_DC_vertical_electric_dipole: compares the 2D and 3D cyl codes with an analytic for a vertical electric dipole (DC) that is offset from the axis of symmetry
- electric_3D_EB_HJ: compares the two different discretizations of Maxwell's equations (EB, HJ) for an electric dipole positioned on the axis of symmetry
- extension_5_FDEM_following_Augustin_Fig3_pipe_radius: extension of notebook 5 above to examine the impact of pipe radius on the extent of the static shielding region

Launching the notebooks

Online

The notebooks can be run online through mybinder or azure notebooks.

Locally

To run them locally, you will need to have python installed, preferably through anaconda.

You can then clone this repository. From a command line, run

git clone https://github.com/simpeg-research/heagy-2018-emcyl.git

Then cd into the heagy-2018-emcyl

cd heagy-2018-emcyl

To setup your software environment, we recommend you use the provided conda environment

conda env create -f environment.yml
source activate emcyl-environment

alternatively, you can install dependencies through pypi

pip install -r requirements.txt

You can then launch Jupyter

jupyter notebook

Jupyter will then launch in your web-browser.

Running the notebooks

Each cell of code can be run with shift + enter or you can run the entire notebook by selecting cell, Run All in the toolbar.

For more information on running Jupyter notebooks, see the Jupyter Documentation

Citation

Heagy, L. J., & Oldenburg, D. W. (2018). Modeling electromagnetics on cylindrical meshes with applications to steel-cased wells. http://arxiv.org/abs/1804.07991

@article{Heagy2018a,
    archivePrefix = {arXiv},
    arxivId = {1804.07991},
    author = {Heagy, Lindsey J. and Oldenburg, Douglas W.},
    eprint = {1804.07991},
    month = {apr},
    title = {{Modeling electromagnetics on cylindrical meshes with applications to steel-cased wells}},
    url = {http://arxiv.org/abs/1804.07991},
    year = {2018}
}

Issues

If you run into problems or bugs, please let us know by creating an issue in this repository.

License

These notebooks are licensed under the MIT License which allows academic and commercial re-use and adaptation of this work.

Version

Version: 0.0.2

simpeg-research / heagy-2018-emcyl

readme