PDE-based geophysical modelling using finite elements: examples from 3D resistivity and 2D magnetotellurics

We present a general finite-element solver, escript, tailored to solve geophysical forward and inverse modeling problems in terms of partial differential equations (PDEs) with suitable boundary conditions. Escript's abstract interface allows geoscientists to focus on solving the actual problem without being experts in numerical modeling. General-purpose finite element solvers have found wide use especially in engineering fields and find increasing application in the geophysical disciplines as these offer a single interface to tackle different geophysical problems. These solvers are useful for data interpretation and for research, but can also be a useful tool in educational settings. This paper serves as an introduction into PDE-based modeling with escript where we demonstrate in detail how escript is used to solve two different forward modeling problems from applied geophysics (3D DC resistivity and 2D magnetotellurics). Based on these two different cases, other geophysical modeling work can easily be realized. The escript package is implemented as a Python library and allows the solution of coupled, linear or non-linear, time-dependent PDEs. Parallel execution for both shared and distributed memory architectures is supported and can be used without modifications to the scripts.

[1]  H. Maurer,et al.  Advances in three-dimensional geoelectric forward solver techniques , 2009 .

[2]  Peter N. Shive,et al.  Singularity removal: A refinement of resistivity modeling techniques , 1989 .

[3]  S. Dyksterhuis,et al.  Finite-element modelling of contemporary and palaeo-intraplate stress using ABAQUS™ , 2005, Comput. Geosci..

[4]  Thomas Poulet,et al.  escript-RT: Reactive transport simulation in Python using escript , 2012, Comput. Geosci..

[5]  Matthew J. Yedlin,et al.  Some refinements on the finite-difference method for 3-D dc resistivity modeling , 1996 .

[6]  David Andrs,et al.  A parallel, fully coupled, fully implicit solution to reactive transport in porous media using the preconditioned Jacobian-Free Newton-Krylov Method , 2013 .

[7]  Brian R. Spies,et al.  The use and misuse of apparent resistivity in electromagnetic methods , 1986 .

[8]  Christopher C. Pain,et al.  Effective multidimensional resistivity inversion using finite-element techniques , 2002 .

[9]  Lutz Gross,et al.  Large Scale Joint Inversion of Geophysical Data using the Finite Element Method in escript , 2013 .

[10]  P. Silvester,et al.  A finite-element program package for magnetotelluric modelling , 1984 .

[11]  A. Dey,et al.  Resistivity modeling for arbitrarily shaped three-dimensional structures , 1979 .

[12]  Yuguo Li,et al.  Three‐dimensional DC resistivity forward modelling using finite elements in comparison with finite‐difference solutions , 2002 .

[13]  Qinzhong Ma.  The boundary element method for 3-D dc resistivity modeling in layered earth , 2002 .

[14]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[15]  Coupled hydromechanical modeling to study the integrity and safety of geological storage of CO2 , 2009 .

[16]  Michael A. Mooney,et al.  Saline pulse test monitoring with the self‐potential method to nonintrusively determine the velocity of the pore water in leaking areas of earth dams and embankments , 2012 .

[17]  Sam L. Butler,et al.  Forward modeling of applied geophysics methods using Comsol and comparison with analytical and laboratory analog models , 2012, Comput. Geosci..

[18]  Lutz Gross,et al.  Interface modeling in incompressible media using level sets in Escript , 2007 .

[19]  Chester J. Weiss,et al.  Adaptive finite-element modeling using unstructured grids: The 2D magnetotelluric example , 2005 .

[20]  Michael S. Zhdanov,et al.  Methods for modelling electromagnetic fields Results from COMMEMI—the international project on the comparison of modelling methods for electromagnetic induction , 1997 .

[21]  Vladimir I. Dmitriev,et al.  On two-dimensional interpretation of magnetotelluric soundings , 1998 .

[22]  H. Johansen,et al.  FAST HANKEL TRANSFORMS , 1979 .

[23]  T. Günther,et al.  Three‐dimensional modelling and inversion of dc resistivity data incorporating topography – II. Inversion , 2006 .

[24]  L. Berselli,et al.  Disperse Two-Phase Flows, with Applications to Geophysical Problems , 2014, Pure and Applied Geophysics.

[25]  Alan G. Jones,et al.  The magnetotelluric method : theory and practice , 2012 .

[26]  Vipin Kumar,et al.  Parallel static and dynamic multi‐constraint graph partitioning , 2002, Concurr. Comput. Pract. Exp..

[27]  Z. Ren,et al.  A goal-oriented adaptive finite-element approach for multi-electrode resistivity system , 2014 .

[28]  Sand erosion as an internal boundary value problem , 2015 .

[29]  Timothy A. Davis,et al.  Algorithm 832: UMFPACK V4.3---an unsymmetric-pattern multifrontal method , 2004, TOMS.

[30]  G. Keller,et al.  The Geoelectrical Methods in Geophysical Exploration , 1994 .

[31]  H. Muhlhaus,et al.  Shear Band Formation in Numerical Simulations Applying a Continuum Damage Rheology Model , 2012, Pure and Applied Geophysics.

[32]  J. Wait.  Propagation of radio waves over a stratified ground , 1953 .

[33]  Weerachai Siripunvaraporn,et al.  GJI Geomagnetism, rock magnetism and palaeomagnetism An efficient modified hierarchical domain decomposition for two-dimensional magnetotelluric forward modelling , 2010 .

[34]  Philip E. Wannamaker,et al.  A stable finite element solution for two-dimensional magnetotelluric modelling , 1987 .

[35]  Hans Petter Langtangen,et al.  A FEniCS-Based Programming Framework for Modeling Turbulent Flow by the Reynolds-Averaged Navier-Stokes Equations , 2011, ArXiv.

[36]  G. W. Hohmann,et al.  An investigation of finite-element modeling for electrical and electromagnetic data in three dimensions , 1981 .

[37]  Heather A. Sheldon,et al.  Conforming Finite-Element Methods for Modeling Convection in an Incompressible Rock Matrix , 2013, Transport in Porous Media.

[38]  Induced seismicity of a normal blind undetected reservoir-bounding fault influenced by dissymmetric fractured damage zones , 2014 .

[39]  H. Chauris,et al.  Resistivity Modelling with Topography , 2012 .

[40]  Jeffrey S. Ovall,et al.  A parallel goal-oriented adaptive finite element method for 2.5-D electromagnetic modelling , 2011 .

[41]  Anders Logg,et al.  Automated Solution of Differential Equations by the Finite Element Method: The FEniCS Book , 2012 .

[42]  Klaus Spitzer,et al.  Adaptive unstructured grid finite element simulation of two-dimensional magnetotelluric fields for arbitrary surface and seafloor topography , 2007 .