.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "examples/real/Perth_basin.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code. .. rst-class:: sphx-glr-example-title .. _sphx_glr_examples_real_Perth_basin.py: Perth Basin, Western Australia ================================= Regional stratigraphy and major normal faults of the northern Perth Basin This model covers the northern Perth Basin, Western Australia, including eight major normal faults (among them the Darling Fault) offsetting a Permian-to-Cretaceous sedimentary sequence. .. GENERATED FROM PYTHON SOURCE LINES 11-22 .. code-block:: Python import os # Importing GemPy import gempy as gp import gempy_viewer as gpv # Importing auxiliary libraries import matplotlib matplotlib.rcParams['figure.figsize'] = (20.0, 10.0) .. GENERATED FROM PYTHON SOURCE LINES 23-29 .. code-block:: Python cwd = os.getcwd() if 'examples' not in cwd: data_path = os.getcwd() + '/examples' else: data_path = cwd + '/../..' .. GENERATED FROM PYTHON SOURCE LINES 30-40 .. code-block:: Python geo_model = gp.create_geomodel( project_name='Perth_Basin', extent=[337000, 400000, 6640000, 6710000, -18000, 1000], refinement=6, importer_helper=gp.data.ImporterHelper( path_to_orientations=data_path + "/data/input_data/perth_basin/Paper_GU2F_sc_faults_topo_Foliations.csv", path_to_surface_points=data_path + "/data/input_data/perth_basin/Paper_GU2F_sc_faults_topo_Points.csv", ) ) .. GENERATED FROM PYTHON SOURCE LINES 41-43 .. code-block:: Python geo_model.structural_frame .. raw:: html
Structural Groups: StructuralGroup:
Name:default_formation
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:Abrolhos_Transfer

StructuralElement:
Name:Cadda

StructuralElement:
Name:Cattamarra

StructuralElement:
Name:Coomallo

StructuralElement:
Name:Cretaceous

StructuralElement:
Name:Darling

StructuralElement:
Name:Eneabba

StructuralElement:
Name:Eneabba_South

StructuralElement:
Name:Hypo_fault_E

StructuralElement:
Name:Hypo_fault_W

StructuralElement:
Name:Lesueur

StructuralElement:
Name:Permian

StructuralElement:
Name:Urella_North

StructuralElement:
Name:Urella_South

StructuralElement:
Name:Woodada_Kockatea

StructuralElement:
Name:Yarragadee
Fault Relations:
default_fo...
default_formation
True
False


.. GENERATED FROM PYTHON SOURCE LINES 44-50 .. code-block:: Python del_surfaces = ['Cadda', 'Woodada_Kockatea', 'Cattamarra'] for s in del_surfaces: gp.remove_element_by_name(geo_model, s) geo_model.structural_frame .. raw:: html
Structural Groups: StructuralGroup:
Name:default_formation
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:Abrolhos_Transfer

StructuralElement:
Name:Coomallo

StructuralElement:
Name:Cretaceous

StructuralElement:
Name:Darling

StructuralElement:
Name:Eneabba

StructuralElement:
Name:Eneabba_South

StructuralElement:
Name:Hypo_fault_E

StructuralElement:
Name:Hypo_fault_W

StructuralElement:
Name:Lesueur

StructuralElement:
Name:Permian

StructuralElement:
Name:Urella_North

StructuralElement:
Name:Urella_South

StructuralElement:
Name:Yarragadee
Fault Relations:
default_fo...
default_formation
True
False


.. GENERATED FROM PYTHON SOURCE LINES 51-66 .. code-block:: Python gp.map_stack_to_surfaces( gempy_model=geo_model, mapping_object={ "fault_Abrolhos_Transfer": ["Abrolhos_Transfer"], "fault_Coomallo": ["Coomallo"], "fault_Eneabba_South": ["Eneabba_South"], "fault_Hypo_fault_W": ["Hypo_fault_W"], "fault_Hypo_fault_E": ["Hypo_fault_E"], "fault_Urella_North": ["Urella_North"], "fault_Darling": ["Darling"], "fault_Urella_South": ["Urella_South"], "Sedimentary_Series": ['Cretaceous', 'Yarragadee', 'Eneabba', 'Lesueur', 'Permian'] } ) .. raw:: html
Structural Groups: StructuralGroup:
Name:fault_Abrolhos_Transfer
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:Abrolhos_Transfer

StructuralGroup:
Name:fault_Coomallo
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:Coomallo

StructuralGroup:
Name:fault_Eneabba_South
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:Eneabba_South

StructuralGroup:
Name:fault_Hypo_fault_W
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:Hypo_fault_W

StructuralGroup:
Name:fault_Hypo_fault_E
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:Hypo_fault_E

StructuralGroup:
Name:fault_Urella_North
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:Urella_North

StructuralGroup:
Name:fault_Darling
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:Darling

StructuralGroup:
Name:fault_Urella_South
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:Urella_South

StructuralGroup:
Name:Sedimentary_Series
Structural Relation:StackRelationType.ERODE
Elements:
StructuralElement:
Name:Cretaceous

StructuralElement:
Name:Yarragadee

StructuralElement:
Name:Eneabba

StructuralElement:
Name:Lesueur

StructuralElement:
Name:Permian
Fault Relations:
fault_Abro...fault_Coom...fault_Enea...fault_Hypo...fault_Hypo...fault_Urel...fault_Darl...fault_Urel...Sedimentar...
fault_Abrolhos_Transfer
fault_Coomallo
fault_Eneabba_South
fault_Hypo_fault_W
fault_Hypo_fault_E
fault_Urella_North
fault_Darling
fault_Urella_South
Sedimentary_Series
True
False


.. GENERATED FROM PYTHON SOURCE LINES 67-70 Select which series are faults ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ .. GENERATED FROM PYTHON SOURCE LINES 70-87 .. code-block:: Python gp.set_is_fault( geo_model, [ "fault_Abrolhos_Transfer", "fault_Coomallo", "fault_Eneabba_South", "fault_Hypo_fault_W", "fault_Hypo_fault_E", "fault_Urella_North", "fault_Darling", "fault_Urella_South" ], ) print(geo_model.structural_frame.fault_relations) .. rst-class:: sphx-glr-script-out .. code-block:: none [[False False False False False False False False True] [False False False False False False False False True] [False False False False False False False False True] [False False False False False False False False True] [False False False False False False False False True] [False False False False False False False False True] [False False False False False False False False True] [False False False False False False False False True] [False False False False False False False False False]] .. GENERATED FROM PYTHON SOURCE LINES 88-90 .. code-block:: Python gpv.plot_2d(geo_model, direction=['z']) .. image-sg:: /examples/real/images/sphx_glr_Perth_basin_001.png :alt: Cell Number: mid Direction: z :srcset: /examples/real/images/sphx_glr_Perth_basin_001.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none .. GENERATED FROM PYTHON SOURCE LINES 91-93 .. code-block:: Python gp.set_topography_from_random(geo_model.grid) .. rst-class:: sphx-glr-script-out .. code-block:: none [-2800. 1000.] Active grids: GridTypes.OCTREE|TOPOGRAPHY|NONE Topography(_regular_grid=RegularGrid(resolution=array([192, 256, 64]), extent=array([ 3.37e+05, 4.00e+05, 6.64e+06, 6.71e+06, -1.80e+04, 1.00e+03]), values=array([[ 3.37164062e+05, 6.64013672e+06, -1.78515625e+04], [ 3.37164062e+05, 6.64013672e+06, -1.75546875e+04], [ 3.37164062e+05, 6.64013672e+06, -1.72578125e+04], ..., [ 3.99835938e+05, 6.70986328e+06, 2.57812500e+02], [ 3.99835938e+05, 6.70986328e+06, 5.54687500e+02], [ 3.99835938e+05, 6.70986328e+06, 8.51562500e+02]], shape=(3145728, 3)), mask_topo=array([], shape=(0, 3), dtype=bool), _transform=None, _base_resolution=array([6, 8, 2])), values_2d=array([[[ 3.37000000e+05, 6.64000000e+06, -8.39147067e+02], [ 3.37000000e+05, 6.64027451e+06, -7.87038174e+02], [ 3.37000000e+05, 6.64054902e+06, -7.43806087e+02], ..., [ 3.37000000e+05, 6.70945098e+06, -1.01702181e+03], [ 3.37000000e+05, 6.70972549e+06, -9.62665055e+02], [ 3.37000000e+05, 6.71000000e+06, -9.02652722e+02]], [[ 3.37329843e+05, 6.64000000e+06, -8.45509106e+02], [ 3.37329843e+05, 6.64027451e+06, -8.07154510e+02], [ 3.37329843e+05, 6.64054902e+06, -7.77055023e+02], ..., [ 3.37329843e+05, 6.70945098e+06, -9.91611095e+02], [ 3.37329843e+05, 6.70972549e+06, -9.45242811e+02], [ 3.37329843e+05, 6.71000000e+06, -8.97111509e+02]], [[ 3.37659686e+05, 6.64000000e+06, -8.25487379e+02], [ 3.37659686e+05, 6.64027451e+06, -7.88725276e+02], [ 3.37659686e+05, 6.64054902e+06, -7.43281635e+02], ..., [ 3.37659686e+05, 6.70945098e+06, -9.23235105e+02], [ 3.37659686e+05, 6.70972549e+06, -8.77676638e+02], [ 3.37659686e+05, 6.71000000e+06, -8.57097627e+02]], ..., [[ 3.99340314e+05, 6.64000000e+06, -1.26479236e+03], [ 3.99340314e+05, 6.64027451e+06, -1.18979385e+03], [ 3.99340314e+05, 6.64054902e+06, -1.13495191e+03], ..., [ 3.99340314e+05, 6.70945098e+06, -1.34370283e+03], [ 3.99340314e+05, 6.70972549e+06, -1.35875555e+03], [ 3.99340314e+05, 6.71000000e+06, -1.32592610e+03]], [[ 3.99670157e+05, 6.64000000e+06, -1.19395572e+03], [ 3.99670157e+05, 6.64027451e+06, -1.11955152e+03], [ 3.99670157e+05, 6.64054902e+06, -1.08116437e+03], ..., [ 3.99670157e+05, 6.70945098e+06, -1.30324540e+03], [ 3.99670157e+05, 6.70972549e+06, -1.32203653e+03], [ 3.99670157e+05, 6.71000000e+06, -1.27123499e+03]], [[ 4.00000000e+05, 6.64000000e+06, -1.13597243e+03], [ 4.00000000e+05, 6.64027451e+06, -1.06163021e+03], [ 4.00000000e+05, 6.64054902e+06, -1.01937065e+03], ..., [ 4.00000000e+05, 6.70945098e+06, -1.26587597e+03], [ 4.00000000e+05, 6.70972549e+06, -1.24971634e+03], [ 4.00000000e+05, 6.71000000e+06, -1.20352799e+03]]], shape=(192, 256, 3)), source=None, values=array([[ 3.37000000e+05, 6.64000000e+06, -8.39147067e+02], [ 3.37000000e+05, 6.64027451e+06, -7.87038174e+02], [ 3.37000000e+05, 6.64054902e+06, -7.43806087e+02], ..., [ 4.00000000e+05, 6.70945098e+06, -1.26587597e+03], [ 4.00000000e+05, 6.70972549e+06, -1.24971634e+03], [ 4.00000000e+05, 6.71000000e+06, -1.20352799e+03]], shape=(49152, 3)), resolution=(192, 256), raster_shape=()) .. GENERATED FROM PYTHON SOURCE LINES 94-96 .. code-block:: Python gpv.plot_3d(geo_model) .. image-sg:: /examples/real/images/sphx_glr_Perth_basin_002.png :alt: Perth basin :srcset: /examples/real/images/sphx_glr_Perth_basin_002.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none .. GENERATED FROM PYTHON SOURCE LINES 97-105 .. code-block:: Python gp.compute_model( gempy_model=geo_model, engine_config=gp.data.GemPyEngineConfig( backend=gp.data.AvailableBackends.PYTORCH, dtype="float64", ) ) .. rst-class:: sphx-glr-script-out .. code-block:: none Setting Backend To: AvailableBackends.PYTORCH GPU requested but unavailable; falling back to CPU (GEMPY_GPU_FALLBACK=True) Setting Backend To: AvailableBackends.PYTORCH Chunking done: 10 chunks Chunking done: 11 chunks Chunking done: 7 chunks Chunking done: 10 chunks Chunking done: 9 chunks Chunking done: 7 chunks Chunking done: 7 chunks Chunking done: 7 chunks Chunking done: 8 chunks Chunking done: 7 chunks Chunking done: 61 chunks Chunking done: 36 chunks Chunking done: 48 chunks Chunking done: 45 chunks Chunking done: 36 chunks Chunking done: 36 chunks Chunking done: 36 chunks Chunking done: 42 chunks Chunking done: 36 chunks Chunking done: 317 chunks Chunking done: 19 chunks Chunking done: 25 chunks Chunking done: 23 chunks Chunking done: 19 chunks Chunking done: 19 chunks Chunking done: 19 chunks Chunking done: 22 chunks Chunking done: 19 chunks Chunking done: 163 chunks Chunking done: 26 chunks .. raw:: html
Solutions: 6 Octree Levels, 13 DualContouringMeshes


.. GENERATED FROM PYTHON SOURCE LINES 106-108 .. code-block:: Python gpv.plot_2d(geo_model, cell_number="mid") .. image-sg:: /examples/real/images/sphx_glr_Perth_basin_003.png :alt: Cell Number: mid Direction: y :srcset: /examples/real/images/sphx_glr_Perth_basin_003.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none .. GENERATED FROM PYTHON SOURCE LINES 109-111 .. code-block:: Python gpv.plot_2d(geo_model, cell_number="mid", series_n=-1, show_scalar=True) .. image-sg:: /examples/real/images/sphx_glr_Perth_basin_004.png :alt: Cell Number: mid Direction: y :srcset: /examples/real/images/sphx_glr_Perth_basin_004.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none .. GENERATED FROM PYTHON SOURCE LINES 112-114 .. code-block:: Python gpv.plot_2d(geo_model, cell_number=[12], direction=["y"], show_data=True, show_topography=True) .. image-sg:: /examples/real/images/sphx_glr_Perth_basin_005.png :alt: Cell Number: 12 Direction: y :srcset: /examples/real/images/sphx_glr_Perth_basin_005.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none .. GENERATED FROM PYTHON SOURCE LINES 115-116 sphinx_gallery_thumbnail_number = 6 .. GENERATED FROM PYTHON SOURCE LINES 116-116 .. code-block:: Python gpv.plot_3d(geo_model, show_lith=True, show_boundaries=True, show_topography=True, ve=None) .. image-sg:: /examples/real/images/sphx_glr_Perth_basin_006.png :alt: Perth basin :srcset: /examples/real/images/sphx_glr_Perth_basin_006.png :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none .. rst-class:: sphx-glr-timing **Total running time of the script:** (1 minutes 40.010 seconds) .. _sphx_glr_download_examples_real_Perth_basin.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: Perth_basin.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: Perth_basin.py ` .. container:: sphx-glr-download sphx-glr-download-zip :download:`Download zipped: Perth_basin.zip ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_