Learn about the theory behind GemPy

The core algorithm of GemPy is based on a universal cokriging interpolation method devised by Lajaunie et al. (1997) and extended by Calcagno et al. (2008). Its implicit nature allows the user to automatically generate complex 3D structural geological models through the interpolation of: (1) surface contact points (3D coordinates of points marking the boundaries between different features, such as layer interfaces, fault planes, unconformities), and (2) orientation measurements (orientations of the poles perpendicular to the dipping of surfaces at any point in the 3D space).

The theory behind GemPy has been described in details in the publication "GemPy 1.0: open-source stochastic geological modeling and inversion" by de la Varga et al. 2019. You can view and download the scientific paper for free:

  • Calcagno, P., Chilès, J. P., Courrioux, G., & Guillen, A. (2008). Geological modelling from field data and geological knowledge: Part I. Modelling method coupling 3D potential-field interpolation and geological rules. Physics of the Earth and Planetary Interiors, 171(1-4), 147-157.

  • Lajaunie, C., Courrioux, G., & Manuel, L. (1997). Foliation fields and 3D cartography in geology: principles of a method based on potential interpolation. Mathematical Geology, 29(4), 571-584.

Read the code documentation

Extensive documentation of GemPy's code can be found on Read the Docs:


© 2019 by the GemPy Team.