gempy-project / gempy

GemPy is an open-source, Python-based 3-D structural geological modeling software, which allows the implicit (i.e. automatic) creation of complex geological models from interface and orientation data. It also offers support for stochastic modeling to address parameter and model uncertainties.
https://gempy.org
European Union Public License 1.2
980 stars 233 forks source link

How layers influence each other? #314

Closed domist07 closed 4 years ago

domist07 commented 4 years ago

How do the layers within a series influence each other during the calculation in GemPy? Regarding my attached graphic I think that currently variant 1 is the case. But in principle, in my opinion, it would be possible to implement variant 2 in GemPy. But for this, the definition of the layers would have to be adapted in a way that it is possible to define for each layer individually by which other layers it is influenced and one would have to be able to define in which order the calculation of the layers is executed. SurfaceInfluence

javoha commented 4 years ago

Hi domist07, the interpolation in GemPy is based on Universal Cokriging. Thus all surfaces (meaning the input data related to these surfaces) that belong to one series influence each other. This is also desired as we consider these geolgoic units to have sedimented under roughly the same geometric setting, forcing them to be more or less parallel. In case of unconformities new series allow for a separate interpolation for a different set of surfaces. For more backgorund maybe check out the original literature:

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.

Can you maybe give an example of a geologic setting where you want one surface to influence the location of another but not vice versa?

Edit: I understand that you want to go from top to bottom, but again, can you elaborate on why?

domist07 commented 4 years ago

Thanks for the explanation. I have no example where this is absolutely necessary. The question came to my mind because I have seen that in Leapfrog Works the layers are calculated one after the other and you can specify which ones affect them. I was asking myself, how this makes a difference in the resulting model.