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
941 stars 232 forks source link

Modeling sedimentary stratas or alluvial aquifer #849

Closed Alfaruuki closed 7 months ago

Alfaruuki commented 7 months ago

Please, is it possible to use gempy to create a 3d geological model of alluvial layers such that the heterogeneity is displayed in each layers of the model.

Specifically, from the attached model output, you will see each geologic unit spread allover the surface which is not true for sedimentary stratas.

image

I am expecting heterogeneity in the layers. for example something like the bellow.

image

Am I missing something? please help......

Alfaruuki commented 7 months ago

Please, share link to videos that can help me. Thank you

javoha commented 7 months ago

Hello, the models you are showing in your example have a high variation on a rather small scale. This is not necesssarily the best use case for gempys univeral co-kriging approach. I am no expert but I think there are several publications with other approaches to better represent such setttings. What you can do: Use gempy to model bigger struccutral units and thern perform something like SGS or SIS within these units to create models similar to your example. It all depends on the amount of data you have and the purpose of your model.

Alfaruuki commented 7 months ago

Thank you for your response Javoha, It really meant alot to me. But what do you mean by SGS or SIS. To be more explicit about the purpose of my model. I would like to couple the resulting geological model with Flopy. Hence the variability within each layers is very important to me. Find attach for another sample. Please help.

image

javoha commented 7 months ago

I think you should look a little more into the mathematical foundation behind gempy to be able to evaluate if it works for your use-case. SGS is Sequential Gaussian Simulation, SIS is Sequential Indicator Simulation. But whatever approach you choose really depends on your available data, purpose, etc.

Alfaruuki commented 7 months ago

Okay. Thank you for your response