Creating a 3D velocity model in GemPy

[adleram]

Hi All,

Started recently to work with GemPy, targeting automatic generation of realistic 3D velocity models. So far, I have not found a simple way to perform this, and wondered if anyone has accomplished this before, or saw sample code that performs this?

Best,

Amir

[javoha]

Hello @adleram,

thanks for bringing this to github. Just for everyone to follow: This is not implemented directly in Gempy at the moment - But as you have the lithology for each grid location it should be straightforward to calculate velocity by combining lithology and some form of depth relation.

Just to clarify access, after computing the model:

sol.lith_block - Gives you the lithology values for each grid point.
sol.grid.values - Gives you the corresponding coordinates of the grid points

Again, I hope this is helpfull and feel free to ask more questions or contribute your function to the repository if you code a nice solution.

Cheers, Jan

[AlexanderJuestel]

This should easily be doable I think when attributing velocities to each lithology. In that case, we would assume constant velocities for the same lithology at depth. But it could be a first approximation. We could provide such functionality in GemGIS as post-processing option if wished @javoha :)

[adleram]

Thanks Alexaner, we are very interested in coding such capability. We have tried to accomplish this, however, could not find a way to access the full grid information of the model, in order to assign a certain velocity to a corresponding lithology. Would you be bale to provide a baseline python script that points the required main line of code? we thank you in a advance for your help, and we will share our complete solution with the entire Gempy community.

[AlexanderJuestel]

Hello,

I will be moving that for now to GemGIS. I have no update on that yet but this way it will not get lost here in GemPy :)

[AlexanderJuestel]

By combining this with density, even a simple synthetic seismic section or cube could be calculated.

[AlexanderJuestel]

Closing because of duplicate of #64