Greenstone.

# Importing gempy
import gempy as gp

# Aux imports
import numpy as np
import matplotlib.pyplot as plt
import os

print(gp.__version__)

Out:

2.2.6
geo_model = gp.create_model('Greenstone')
data_path = 'https://raw.githubusercontent.com/cgre-aachen/gempy_data/master/'

# Importing the data from csv files and settign extent and resolution
geo_model = gp.init_data(geo_model, [696000, 747000, 6863000, 6930000, -20000, 200], [50, 50, 50],
                         path_o=data_path + "/data/input_data/tut_SandStone/SandStone_Foliations.csv",
                         path_i=data_path + "/data/input_data/tut_SandStone/SandStone_Points.csv")

Out:

Active grids: ['regular']
gp.plot_2d(geo_model, direction=['z'])
Cell Number: mid Direction: z

Out:

<gempy.plot.visualization_2d.Plot2D object at 0x7ff9774163d0>
gp.map_stack_to_surfaces(geo_model, {"EarlyGranite_Series": 'EarlyGranite',
                                     "BIF_Series": ('SimpleMafic2', 'SimpleBIF'),
                                     "SimpleMafic_Series": 'SimpleMafic1', 'Basement': 'basement'})
surface series order_surfaces color id
3 EarlyGranite EarlyGranite_Series 1 #728f02 1
0 SimpleMafic2 BIF_Series 1 #015482 2
1 SimpleBIF BIF_Series 2 #9f0052 3
2 SimpleMafic1 SimpleMafic_Series 1 #ffbe00 4
4 basement Basement 1 #443988 5


geo_model.add_surface_values([2.61, 2.92, 3.1, 2.92, 2.61])
surface series order_surfaces color id value_0
3 EarlyGranite EarlyGranite_Series 1 #728f02 1 2.61
0 SimpleMafic2 BIF_Series 1 #015482 2 2.92
1 SimpleBIF BIF_Series 2 #9f0052 3 3.10
2 SimpleMafic1 SimpleMafic_Series 1 #ffbe00 4 2.92
4 basement Basement 1 #443988 5 2.61


gp.set_interpolator(geo_model,
                    compile_theano=True,
                    theano_optimizer='fast_compile',
                    verbose=[])

Out:

Setting kriging parameters to their default values.
Compiling theano function...
Level of Optimization:  fast_compile
Device:  cpu
Precision:  float64
Number of faults:  0
Compilation Done!
Kriging values:
                        values
range                 8.7e+04
$C_o$                 1.8e+08
drift equations  [3, 3, 3, 3]

<gempy.core.interpolator.InterpolatorModel object at 0x7ff9c9afb9d0>
gp.compute_model(geo_model, set_solutions=True)

Out:

Lithology ids
  [5. 5. 5. ... 5. 5. 5.]
gp.plot_2d(geo_model, cell_number=[-1], direction=['z'], show_data=False)
Cell Number: -1 Direction: z

Out:

<gempy.plot.visualization_2d.Plot2D object at 0x7ff973cba710>
gp.plot_2d(geo_model, cell_number=[25], direction='x')
Cell Number: 25 Direction: x

Out:

<gempy.plot.visualization_2d.Plot2D object at 0x7ff973e79e50>

Out:

array([[2.61, 2.61, 2.61, ..., 2.61, 2.61, 2.61]])
p2d = gp.plot_2d(geo_model, cell_number=[25], block=geo_model.solutions.values_matrix,
           direction=['y'], show_data=True,
           kwargs_regular_grid={'cmap': 'viridis', 'norm':None})
Cell Number: 25 Direction: y

sphinx_gallery_thumbnail_number = 5

Greenstone

Out:

<gempy.plot.vista.GemPyToVista object at 0x7ff9fa5e5650>

Saving the model

gp.save_model(geo_model, path=os.pardir + ‘/data/gempy_models’)

gp.save_model(geo_model)

Out:

True

Total running time of the script: ( 0 minutes 19.272 seconds)

Gallery generated by Sphinx-Gallery