Model 1 - Horizontal stratigraphicΒΆ

This is the most simple model of horizontally stacked layers. We start by importing the necessary dependencies:

Importing GemPy

import gempy as gp

import pandas as pd
pd.set_option('precision', 2)

Creating the model by importing the input data and displaying it:

data_path = ''
geo_data = gp.create_data('horizontal', extent=[0, 1000, 0, 1000, 0, 1000], resolution=[50, 50, 50],
                          path_o=data_path + "/data/input_data/jan_models/model1_orientations.csv",
                          path_i=data_path + "/data/input_data/jan_models/model1_surface_points.csv")


Active grids: ['regular']

Setting and ordering the units and series:

gp.map_stack_to_surfaces(geo_data, {"Strat_Series": ('rock2', 'rock1'), "Basement_Series": ('basement')})
surface series order_surfaces color id
0 rock2 Strat_Series 1 #015482 1
1 rock1 Strat_Series 2 #9f0052 2
2 basement Basement_Series 1 #ffbe00 3

gp.plot_2d(geo_data, direction=['y'])
Cell Number: mid Direction: y


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

Calculating the model:

interp_data = gp.set_interpolator(geo_data, compile_theano=True,


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:
range            1.7e+03
$C_o$            7.1e+04
drift equations   [3, 3]

Displaying the result in x and y direction:

gp.plot_2d(geo_data, cell_number=[25],
           direction=['x'], show_data=True)
Cell Number: 25 Direction: x


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

sphinx_gallery_thumbnail_number = 2

gp.plot_2d(geo_data, cell_number=[25],
           direction=['y'], show_data=True)

Cell Number: 25 Direction: y



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

Gallery generated by Sphinx-Gallery