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']
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'])
Cell Number: mid Direction: y

Out:

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

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


Displaying the result in y and x direction:

gp.plot_2d(geo_data, cell_number=15,
           direction='y', show_data=True)
Cell Number: 15 Direction: y

Out:

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

sphinx_gallery_thumbnail_number = 2

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

gp.save_model(geo_data)
Cell Number: 25 Direction: x

Out:

True

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

Gallery generated by Sphinx-Gallery