Python for HPC geophysical electromagnetic applications: experiences and perspectives

Nowadays, the electromagnetic modelling are a fun-damental tool in geophysics due to their wide field of application: hydrocarbon and mineral exploration, reservoir monitoring, CO��� storage characterization, geothermal reservoir imaging and many others. In particular, the 3D CSEM forward modelling (FM) is an established tool in the oil & gas industry because of the hope that the application of such methods would eventually lead to the direct detection of hydrocarbons through their insulating properties. Although 3D CSEM FM is nowadays a well-known geophysical prospecting tool and his fundamental mathematical theory is well-established, the state-of-art shows a relative scarsity of robust, flexible, modular and open-source codes to simulate these problems on HPC platforms, which is crucial in the future goal of solving inverse problems. In this talk we describe our experience and perspectives in the development of an HPC python code for the 3D CSEM FM, namely, PETGEM. We focus on three points: 1) 3D CSEM FM theory from a practical point of view, 2) PETGEM features and Python potential for HPC applications, and 3) Modelling results of real-life 3D CSEM FM cases. These points depict that PETGEM could be an attractive and competitive HPC tool to simulate real-scale of 3D CSEM FM in geophysics.

[1]  Art Raiche,et al.  Efficient solution of full domain 3D electromagnetic modelling problems , 2000 .

[2]  José María Cela,et al.  Improving edge finite element assembly for geophysical electromagnetic modelling on shared-memory architectures , 2016, 2016 IEEE 7th Annual Ubiquitous Computing, Electronics & Mobile Communication Conference (UEMCON).

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

[4]  Gregory A. Newman,et al.  Three-dimensional induction logging problems, Part 2: A finite-difference solution , 2002 .

[5]  J. Shadid,et al.  Three‐dimensional wideband electromagnetic modeling on massively parallel computers , 1996 .

[6]  J. T. Smith,et al.  Three-dimensional electromagnetic modeling using finite difference equations: The magnetotelluric example , 1994 .

[7]  Shikun Dai,et al.  Finite element modelling of marine controlled-source electromagnetic responses in two-dimensional dipping anisotropic conductivity structures , 2011 .

[8]  Adaptive finite-element modeling using unstructured grids : The 2 D magnetotelluric example , 2018 .

[9]  José María Cela,et al.  Edge-based parallel framework for the simulation of 3d CSEM surveys , 2016 .

[10]  Karsiti,et al.  1 D Modeling of Controlled-Source Electromagnetic ( CSEM ) Data using Finite Element Method for Hydrocarbon Detection in Shallow Water , 2011 .

[11]  Jelena Koldan Numerical solution of 3-D electromagnetic problems in exploration geophysics and its implementation on massively parallel computers , 2013 .