search

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

Define velocity property and export the property model #942

Open jcgeophysic opened 1 week ago

jcgeophysic commented 1 week ago

Hi! I hope you are doing well! I'm new user of Gempy and I tried to follow videos and tutorials for building a velocity model to use it later in Devito. Many things have changed, so I already built a model but I can not assing the velocity value to each layer in the model. It seems that GeophysicsInput module just work for gravity models. At the end, I like to export the velocity model, however, the way that is explained in the transform 21 video, is not working. I really appreciate your help.

Best regards

javoha commented 1 week ago

Hi, thanks for the question. Maybe I can offer a basic solution. When creating a gempy model we create a regular grid (at least as the default). The coordinates for this grid are stored in geo_data.grid.regular_grid.values. For each coordinate we have a corresponding lithology integer ID stored in geo_data.solutions.raw_arrays.lith_block. If you want to set a single specific velocity per rock unit (structural element) you can create an array of the same size and replace each lithology ID with a velocity value. Going more advanced you could also populate the volumes per rock with unit interpolated fields (e.g. using RBF, Kriging) if you have measurements or with Random fields (e.g. using Sequential Gaussian Simulation) if you dont. Does this help/answer your question? Cheers Jan

jcgeophysic commented 1 week ago

Hi Javoha, yes, it helped! Thank you very much.

Regards,

Jhonny

javoha commented 1 week ago

That's good to hear, let me know if you need more help.