Model 4 - PinchoutΒΆ

This script demonstrates how to create a geological model with a layer of varying thickness (pinchout) using GemPy, a Python-based, open-source library for implicit geological modeling.

# Import necessary libraries
import gempy as gp
import gempy_viewer as gpv

# sphinx_gallery_thumbnail_number = 3

Generate the model Define the path to data

data_path = ''
path_to_data = data_path + "/data/input_data/jan_models/"
# Create a GeoModel instance
data = gp.create_geomodel(
    extent=[0, 1000, 0, 1000, 0, 1000],
    resolution=[50, 50, 50],
        path_to_orientations=path_to_data + "model4_orientations.csv",
        path_to_surface_points=path_to_data + "model4_surface_points.csv"
# Map geological series to surfaces
    mapping_object={"Strat_Series": ('rock2', 'rock1')}
# Compute the geological model
geo_data = data
Surface points hash:  76e41015e8bc57232f56f90ff5f9fc56d0f5d55b67a5030f9bc5c95edfed0674
Orientations hash:  575806bb62944f0604e8fecccdc9d53bfabd758bee293b1ad8f3575f173429f1
Setting Backend To: AvailableBackends.PYTORCH

Plot the initial geological model in the y direction

gpv.plot_2d(geo_data, direction=['y'], show_results=False)

# Plot the result of the model in the x and y direction with data
gpv.plot_2d(geo_data, direction=['x'], show_data=True)
gpv.plot_2d(geo_data, direction=['y'], show_data=True)
  • Cell Number: mid Direction: y
  • Cell Number: mid Direction: x
  • Cell Number: mid Direction: y
<gempy_viewer.modules.plot_2d.visualization_2d.Plot2D object at 0x7f00cb4856f0>
gpv.plot_3d(geo_data, show_surfaces=True)
d04 pinchout
< object at 0x7f00cb60bd90>

Total running time of the script: (0 minutes 1.003 seconds)

Gallery generated by Sphinx-Gallery