1.3b: 2-D sections

Importing

import gempy as gp
import gempy_viewer as gpv
import numpy as np

np.random.seed(1234)

Setup the model

Importing the data from CSV-files and setting extent and resolution

data_path = 'https://raw.githubusercontent.com/cgre-aachen/gempy_data/master/'

geo_model: gp.data.GeoModel = gp.create_geomodel(
    project_name='Tutorial_ch1_1_Basics',
    extent=[0, 2000, 0, 2000, 0, 750],
    resolution=[20, 20, 20],  # * Here we define the resolution of the voxels
    refinement=4,  # * Here we define the number of octree levels. If octree levels are defined, the resolution is ignored.
    importer_helper=gp.data.ImporterHelper(
        path_to_orientations=data_path + "/data/input_data/getting_started/simple_fault_model_orientations.csv",
        path_to_surface_points=data_path + "/data/input_data/getting_started/simple_fault_model_points.csv",
        hash_surface_points="4cdd54cd510cf345a583610585f2206a2936a05faaae05595b61febfc0191563",
        hash_orientations="7ba1de060fc8df668d411d0207a326bc94a6cdca9f5fe2ed511fd4db6b3f3526"
    )
)

gp.map_stack_to_surfaces(
    gempy_model=geo_model,
    mapping_object=  # TODO: This mapping I do not like it too much. We should be able to do it passing the data objects directly
    {
        "Fault_Series": 'Main_Fault',
        "Strat_Series": ('Sandstone_2', 'Siltstone', 'Shale', 'Sandstone_1')
    }
)

gp.set_is_fault(
    frame=geo_model.structural_frame,
    fault_groups=['Fault_Series']
)

Surface points hash:  4cdd54cd510cf345a583610585f2206a2936a05faaae05595b61febfc0191563
Orientations hash:  7ba1de060fc8df668d411d0207a326bc94a6cdca9f5fe2ed511fd4db6b3f3526

Structural Groups:	StructuralGroup: Name: Fault_Series Structural Relation: StackRelationType.FAULT Elements: StructuralElement: Name: Main_Fault StructuralGroup: Name: Strat_Series Structural Relation: StackRelationType.ERODE Elements: StructuralElement: Name: Sandstone_2 StructuralElement: Name: Siltstone StructuralElement: Name: Shale StructuralElement: Name: Sandstone_1
Fault Relations:	Fault_Seri... Strat_Seri... Fault_Series Strat_Series
	True False

Add sections

pass section dictionary with startpoint, endpoint and resolution for every section:

gp.set_section_grid(
    grid=geo_model.grid,
    section_dict={
        'section1': ([0, 0], [2000, 2000], [100, 80]),
        'section2': ([800, 0], [800, 2000], [150, 100]),
        'section3': ([0, 200], [1500, 500], [200, 150])
    }  # p1,p2,resolution
)

Active grids: GridTypes.NONE|SECTIONS|DENSE

	start	stop	resolution	dist
section1	[0, 0]	[2000, 2000]	[100, 80]	2828.427125
section2	[800, 0]	[800, 2000]	[150, 100]	2000.000000
section3	[0, 200]	[1500, 500]	[200, 150]	1529.705854

Add topography

gp.set_topography_from_random(
    grid=geo_model.grid,
    fractal_dimension=1.2,
    d_z=np.array([300, 750]),
    topography_resolution=np.array([50, 50])
)

Active grids: GridTypes.NONE|SECTIONS|TOPOGRAPHY|DENSE

<gempy.core.data.grid_modules.topography.Topography object at 0x7fdad93306d0>

Active grids:

geo_model.grid.active_grids

<GridTypes.NONE|SECTIONS|TOPOGRAPHY|DENSE: 1050>

gpv.plot_section_traces(geo_model)

<function plot_section_traces at 0x7fdaecff2830>

geo_model.interpolation_options.mesh_extraction = False
sol = gp.compute_model(geo_model)

Setting Backend To: AvailableBackends.numpy
Chunking done: 8 chunks

gpv.plot_2d(geo_model, section_names=['topography'])

/Users/flow/opt/anaconda3/envs/gempy3_pre/lib/python3.10/site-packages/gempy_viewer/API/_plot_2d_sections_api.py:105: UserWarning: Section contacts not implemented yet. We need to pass scalar field for the sections grid
  warnings.warn(

<gempy_viewer.modules.plot_2d.visualization_2d.Plot2D object at 0x7fda98ba58d0>

gpv.plot_2d(geo_model, section_names=['section1'])

/Users/flow/opt/anaconda3/envs/gempy3_pre/lib/python3.10/site-packages/gempy_viewer/API/_plot_2d_sections_api.py:105: UserWarning: Section contacts not implemented yet. We need to pass scalar field for the sections grid
  warnings.warn(

<gempy_viewer.modules.plot_2d.visualization_2d.Plot2D object at 0x7fdaed888fd0>

gpv.plot_2d(
    model=geo_model,
    section_names=['section1', 'section2', 'section3', 'topography'],
    show_topography=True
)


# sphinx_gallery_thumbnail_number = 4

/Users/flow/opt/anaconda3/envs/gempy3_pre/lib/python3.10/site-packages/gempy_viewer/API/_plot_2d_sections_api.py:105: UserWarning: Section contacts not implemented yet. We need to pass scalar field for the sections grid
  warnings.warn(

<gempy_viewer.modules.plot_2d.visualization_2d.Plot2D object at 0x7fda98ba52a0>