Sandstone Greenstone Belt, Western Australia

An Archean greenstone belt model from the Yilgarn Craton, WA

This model covers the Sandstone Greenstone Belt in the Yilgarn Craton of Western Australia, combining an early granite intrusion, banded iron formation (BIF), and mafic units above a basement.

import os

# Importing gempy
import gempy as gp
import gempy_viewer as gpv

data_path = os.path.abspath('../../data/input_data/tut_SandStone')

# Importing the data from csv
geo_model = gp.create_geomodel(
    project_name='Greenstone',
    extent=[696000, 747000, 6863000, 6930000, -20000, 200],  # * Here we define the extent of the model
    refinement=6,
    # * 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 + "/SandStone_Foliations.csv",
        path_to_surface_points=data_path + "/SandStone_Points.csv",
        hash_surface_points=None,
        hash_orientations=None
    )
)

gpv.plot_2d(geo_model, direction=['z'])

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

gp.map_stack_to_surfaces(
    gempy_model=geo_model,
    mapping_object={
        "EarlyGranite_Series": 'EarlyGranite',
        "BIF_Series": ('SimpleMafic2', 'SimpleBIF'),
        "SimpleMafic_Series": 'SimpleMafic1', 'Basement': 'basement'
    }
)

Could not find element 'basement' in any group.

Structural Groups:	StructuralGroup: Name: EarlyGranite_Series Structural Relation: StackRelationType.ERODE Elements: StructuralElement: Name: EarlyGranite StructuralGroup: Name: BIF_Series Structural Relation: StackRelationType.ERODE Elements: StructuralElement: Name: SimpleMafic2 StructuralElement: Name: SimpleBIF StructuralGroup: Name: SimpleMafic_Series Structural Relation: StackRelationType.ERODE Elements: StructuralElement: Name: SimpleMafic1
Fault Relations:	EarlyGrani... BIF_Series SimpleMafi... EarlyGranite_Series BIF_Series SimpleMafic_Series
	True False

gp.compute_model(geo_model)

Setting Backend To: AvailableBackends.PYTORCH
GPU requested but unavailable; falling back to CPU (GEMPY_GPU_FALLBACK=True)
Setting Backend To: AvailableBackends.PYTORCH
Chunking done: 9 chunks
Chunking done: 10 chunks
Chunking done: 6 chunks
Chunking done: 44 chunks
Chunking done: 14 chunks
Chunking done: 24 chunks
Chunking done: 22 chunks
Chunking done: 7 chunks
Chunking done: 12 chunks

Solutions: 6 Octree Levels, 4 DualContouringMeshes

gpv.plot_2d(geo_model, cell_number=[-1], direction=['z'], show_data=False)

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

gpv.plot_2d(geo_model, cell_number=['mid'], direction='x')

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

sphinx_gallery_thumbnail_number = -1

gpv.plot_3d(geo_model, show_lith=True, show_boundaries=True, ve=None)

<gempy_viewer.modules.plot_3d.vista.GemPyToVista object at 0x7f99740692b0>