Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Solve oringinal GemPy strata order issue #2

Open
GeorgeLiang3 opened this issue Apr 11, 2023 · 0 comments
Open

Solve oringinal GemPy strata order issue #2

GeorgeLiang3 opened this issue Apr 11, 2023 · 0 comments

Comments

@GeorgeLiang3
Copy link
Owner

Problem: In Original GemPy, different grids configuration are used to compute the geological model and gravity. More specifically, the regular grid is used for geological models, and the centered grid is used to calculate gravity. The two grids are generated differently. The regular grid has the origin (0,0,0) at the lower corner, while the centered grid has the origin at the upper corner. This makes the order of the cell index in two directions. For regular grid, it's increasing upwards; for centered grid, it's increasing downward. This is not an issue of we deal with them separately, like what we have in the original GemPy. However, it becomes a problem if we would like to map the lithologies and the densities to the grid at the same time. It has the lithology and the densities in the wrong order. With the addition of the 'slope' technique, it is even harder to detect the correct value.
For example, a three layer model has lithology index (1,2,3) and densities values(4,5,6). In the original gempy, if mapping the two values to the same grid without care, it will be something like (1,2,3,6,-1,3)

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant