The Theory Behind GemPy
Scientific Background
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:
Figure from de la Varga et al. (2019).
References
-
de la Varga, M., Schaaf, A., and Wellmann, F. (2019). GemPy 1.0: open-source stochastic geological modeling and inversion, Geosci. Model Dev., 12, 1-32.
-
Wellmann, F., & Caumon, G. (2018). 3-D Structural geological models: Concepts, methods, and uncertainties. In Advances in Geophysics (Vol. 59, pp. 1-121). Elsevier.
-
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.