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 0x7f3e1439eed0>

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 enabled. Using device: cuda
GPU device count: 1
Current GPU device: 0
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 0x7f3e3d4da1d0>

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

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

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 0x7f3eea0f1d30>