.. 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. ============ .. GENERATED FROM PYTHON SOURCE LINES 5-17 .. 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) os.environ["aesara_FLAGS"] = "mode=FAST_RUN,device=cuda" .. GENERATED FROM PYTHON SOURCE LINES 18-24 .. 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 25-36 .. code-block:: Python geo_model: gp.data.GeoModel = gp.create_geomodel( project_name='Perth_Basin', extent=[337000, 400000, 6640000, 6710000, -18000, 1000], resolution=[100, 100, 100], refinement=4, 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 37-39 .. 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 40-46 .. 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 47-62 .. 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 63-66 Select which series are faults ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ .. GENERATED FROM PYTHON SOURCE LINES 66-82 .. code-block:: Python gp.set_is_fault( geo_model, fault_groups=[ "fault_Abrolhos_Transfer", "fault_Coomallo", "fault_Eneabba_South", "fault_Hypo_fault_W", "fault_Hypo_fault_E", "fault_Urella_North", "fault_Darling", "fault_Urella_South" ], ) .. raw:: html
Structural Groups: StructuralGroup:
Name:fault_Abrolhos_Transfer
Structural Relation:StackRelationType.FAULT
Elements:
StructuralElement:
Name:Abrolhos_Transfer

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

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

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

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

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

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

StructuralGroup:
Name:fault_Urella_South
Structural Relation:StackRelationType.FAULT
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 83-84 gp.set_fault_relation(geo_model, fr) .. GENERATED FROM PYTHON SOURCE LINES 84-87 .. code-block:: Python 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-89 %matplotlib inline .. GENERATED FROM PYTHON SOURCE LINES 89-91 .. 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 92-94 .. code-block:: Python gp.set_topography_from_random(geo_model.grid) .. rst-class:: sphx-glr-script-out .. code-block:: none [-2800. 1000.] Active grids: ['regular' 'topography'] .. GENERATED FROM PYTHON SOURCE LINES 95-97 .. 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 98-106 .. 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 Chunking done: 21 chunks Chunking done: 21 chunks Chunking done: 21 chunks Chunking done: 21 chunks Chunking done: 21 chunks Chunking done: 21 chunks Chunking done: 21 chunks Chunking done: 21 chunks Chunking done: 21 chunks .. raw:: html
Solutions: 4 Octree Levels, 13 DualContouringMeshes


.. GENERATED FROM PYTHON SOURCE LINES 107-109 .. 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 110-112 .. 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 113-115 .. 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 116-117 sphinx_gallery_thumbnail_number = 6 .. GENERATED FROM PYTHON SOURCE LINES 117-117 .. code-block:: Python gpv.plot_3d(geo_model, show_topography=True) .. 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:** (0 minutes 38.530 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 ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_