Model 2 - Anticline

A simple anticline structure. 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/'
path_to_data = data_path + "/data/input_data/jan_models/"
geo_data = gp.create_data('fold', extent=[0, 1000, 0, 1000, 0, 1000], resolution=[50, 50, 50],
                          path_o=path_to_data + "model2_orientations.csv",
                          path_i=path_to_data + "model2_surface_points.csv")

Out:

Active grids: ['regular']

geo_data.get_data().head()

		X	Y	Z	smooth	surface	G_x	G_y	G_z
surface_points	0	0	100	200	2.00e-06	rock1	NaN	NaN	NaN
	1	0	500	200	2.00e-06	rock1	NaN	NaN	NaN
	2	0	900	200	2.00e-06	rock1	NaN	NaN	NaN
	3	400	100	600	2.00e-06	rock1	NaN	NaN	NaN
	4	400	500	600	2.00e-06	rock1	NaN	NaN	NaN

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	rock1	Strat_Series	1	#015482	1
1	rock2	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 0x7fcb8b7774f0>

Calculating the model:

interp_data = gp.set_interpolator(geo_data, 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]

geo_data.orientations

	X	Y	Z	G_x	G_y	G_z	smooth	surface
1	500	500	620	1.00e-12	1.00e-12	1.0	0.01	rock1
0	500	500	820	1.00e-12	1.00e-12	1.0	0.01	rock2

sol = gp.compute_model(geo_data)

Displaying the result in y and x direction:

gp.plot_2d(geo_data, cell_number=15,
           direction='y', show_data=True)

Out:

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

sphinx_gallery_thumbnail_number = 2

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

gp.save_model(geo_data)

Out:

True