Tutorial: Geophysical modeling & inversion with pyGIMLi

Florian Wagner, Carsten Rücker, Thomas Günther, Andrea Balza

Tutorial Info:
- Tutorial Slack-channel: #t21-mon-pygimli

Many tasks in applied geosciences cannot be solved by a single measurement method and require the integration of geophysical, geological, geotechnical and hydrological methods. Such data integration tasks need flexible open-source software. PyGIMLi is a Python library for modeling and inversion, which provides finite element and finite volume solvers on structured and unstructured 2D/3D meshes, various ready-to-use geophysical forward operators and method managers, as well as a generalized deterministic Gauss-Newton inversion algorithm with flexible means of incorporating prior knowledge (e.g., through region-specific regularization). In this tutorial, we will provide an overwiew of pyGIMLi‘s functionality and applications, solve stationary heat conduction on a 2D mesh in a Jupyter Notebook with pyGIMLi‘s modeling level and create and invert synthetic geophysical data on a user-specified subsurface geometry. The tutorial will conclude with further development plans, integration with other software packages and guidelines to contribute.

= = Content = =

00:00:00 - Introduction
00:01:09 - Main features, conda installer, API doc
00:09:30 - 2D meshtools demonstration
00:31:10 - Equation level: 2D heat equation
00:48:19 - Crosshole traveltime forward modeling
01:04:42 - Method Manager: Traveltime inversion
01:24:47 - Inverting electrical resistivity field data
01:41:18 - Inversion with own forward operator
01:49:54 - Homepage with examples, papers, contribution guide

[chhunkeanthai]

Thanks for this helpful tutorial, Sir. May you please make another tutorials relate to cross-hole electrical resistivity modeling & inversion with pyGIMLi?

[juanacosta]

Does pyGIMLi have inversions for gravity and magnetic data?

[turtle]

thanks for you work with pygimli ! I am a beginner and try to use it but I do not know how to use the data created thanks to the solver to make calculations (as a normal numpy matrix for instance)? For instance with the 2D heat equation I would like to export the Temperature as a function of x, y (as a numpy matrix perhaps ?) . Is there a way to do that ? Thanks for your response

[appliedofexplorationgeophy]

Dear Sir
I have MATLAB code as (seg2read), and I need to modify it to (seg2 output); thus, I am asking for your help.

[abdullahsahil]

Thanks for the good information can you make a video on 3D VP model and VS model by using seismic refraction i will be thankful to you

[G.L.Eng.]

Hi, using pyGIMLi is it possible to do a seismic tomography?