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 = 'https://raw.githubusercontent.com/cgre-aachen/gempy_data/master/'
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")

Out:

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'])

Out:

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

Calculating the model:

interp_data = gp.set_interpolator(geo_data, compile_theano=True,
                                  theano_optimizer='fast_compile')

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             1732.05
$C_o$            71428.57
drift equations    [3, 3]

sol = gp.compute_model(geo_data)

Displaying the result in x and y direction:

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

Out:

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

sphinx_gallery_thumbnail_number = 2

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

gp.save_model(geo_data)

Out:

True