GemPy
Modelling an onlap surface across more than one series #784
Replies: 5 comments 6 replies
-
|
Hi, |
Beta Was this translation helpful? Give feedback.
-
|
Hi Jan! Thanks for responding - please see the data I'm using below! Troian_2020_gempy_input_data.csv Cheers, Luciano PS: Edited the csv file to add a few points and consider for terrain elevation in the subsurface points. |
Beta Was this translation helpful? Give feedback.
-
|
Is that a mix of surface points and orientations? Looks like all orientations on first glance. |
Beta Was this translation helpful? Give feedback.
-
|
Hi Jan, I have the surfaces and orientations on the same file (same number of orientations and surfaces). I do not know if this is not a good practice, but I considered all my horizontal stratigraphy points (units 1, 2, 3, 4 and pelotas_lower) to have orientations (so azimuth = 0, dip = 0 and polarity = 1 for each), and my onlap surface between those and the basement was a trial and error to get to the point seen in the image I sent earlier in the discussion. The onlap surfaces named "dummy_point" are the outcropping basement, which I defined based on topography in QGIS. Hope that helps, otherwise please let me know and I will share more info! Cheers, Luciano |
Beta Was this translation helpful? Give feedback.
-
|
Hey @lucianomarquetto, I removed the "pelotas lower" points, as those looked like they were in the dataframe twice. You will achieve the correct outcropping (geological map) by adding a topography. I renamed the series a little for my convenience and set the bottom relation of the horizontal strata series to "Onlap". gp.map_stack_to_surfaces(geo_model, {"Horizontal_Series": ('unit4', 'unit3','unit2', 'unit1'),
"Onlap_Series": ('onlap_surface5'),
"Basement_Series":('basement')})
geo_model.set_bottom_relation("Horizontal_Series", "Onlap")Is this what you were looking for? |
Beta Was this translation helpful? Give feedback.
-
|
Hi @javoha! This is pretty much what I am looking for at the present, thanks! Based on your comment I was able to reproduce the model in a similar way (took the pelotas_lower out and it worked). I think my hassle came from a confusion I made about onlap surface vs regular surface - I added the "pelotas_lower" to represent the undifferentiated strata below my four units, and the "onlap surface" to define the boundary between the horizontal strata (pelotas_lower and units 1-4) and the basement (so the "pelotas_lower" had fewer points than the "onlap surface"). I understand by reading your comment that having only the onlap surface is enough for this. I went further, based on your code, and defined my four units as one series each, which also worked now. In the future I plan to detail them, so better to already try how it would work. Anyway, one thing that caught my attention now is that my 3D plots are odd. At first I thought I was running something different from you, because in your model I could see the basement (red) in the northwest and the onlap surface (purple) in the southeast, and all surfaces seemed to be in the right depth comparing to the translucent blocks. In my plots (both defining the four units as one series only or one series each), my southeast section becomes red as if basement, and surfaces seem off compared to the translucent block limits. But then I noticed my top surface, which should be blue, is appearing with the color of the second surface:
I believe this is then a pyvista problem? My pyvista is 0.37.0, and I tried plotting as basic and background running Anyway, I will try finding a solution to this color and surface issue, and if you happen to know what I can do please share. So far I can only thank you and the Gempy community for all help and support so far, it has been an enriching experience to work with the package! |
Beta Was this translation helpful? Give feedback.
-
|
Ok just noticed this: #754 Probably my pyvista problems are related to my version then! |
Beta Was this translation helpful? Give feedback.
-
|
Yeah I get the same problem with the colors - for the image I accidentally used one of my older environments where everything was still working properly. Hopefully most of these dependency related issues will be fixed in the next rellease, but of course this will always be part of open-source development. If its super urgent for you to fix this, you could try to export the gempy surfaces and work directly with pyvista to create "your own" plots, but this will probably take some time. Edit: I will actually reopen this, so people can find it easier when looking for help, maybe you can mark it as answered. |
Beta Was this translation helpful? Give feedback.
-
|
Hey Jan, I am not in a hurry to have it right, but will definitely try exporting the model and working directly on pyvista if needed, thanks! Further investigating I noticed it is not pyvista (at least not alone) that is giving me those off surfaces and color problems. I tried creating a new virtual environment with python 3.8.1 and getting the gempy required packages from a requirements.txt file running Would you mind sharing your 'older environment' yml file so I can try reproducing it here? Now that I know it is just the 3D visualization and not the model that is odd I can keep working with no problems, but still would like to try and have a better visualization if I can, until the next gempy release. Cheers! |
Beta Was this translation helpful? Give feedback.
-
|
Hi, name: gempy_env
|
Beta Was this translation helpful? Give feedback.
-
|
Hi Jan, Thanks for sharing your environment - I will try to reproduce it or at least find the specific packages related to 3D visualization and compare to my versions. It sure will help me go further in this matter! When (if) I reach any conclusion will share here! Cheers, Luciano. |
Beta Was this translation helpful? Give feedback.
-
Hi all,
I am currently working on a coastal area with 5 series of horizontal stratigraphy, but with the basement outcropping in the northwest. I am trying to model an onlap surface to represent this basement coming from depth and outcropping, but can't get it to work (the horizontal layers do not compute properly).
This is my stratigraphy:
I was able to model the basement satisfactorily:
And the relation to the horizontal stratigraphy I'm trying to reach is something like below:
What I have already tried:
I am able to get variable but not satisfactory results. Any suggestions? Is there a code function that allows an onlap surface to cross all series above it instead of only the series right above?
I am using gempy 2.2.11 in Ubuntu 22.04.
Thanks!
Beta Was this translation helpful? Give feedback.
All reactions