Eaton Method with Seismic Velocity Data

[2]:
from __future__ import print_function, division, unicode_literals
%matplotlib inline
import matplotlib.pyplot as plt

plt.style.use(['seaborn-paper', 'seaborn-whitegrid'])
plt.rcParams['font.sans-serif']=['SimHei']
plt.rcParams['axes.unicode_minus']=False

import numpy as np

import pygeopressure as ppp

Create survey CUG:

[3]:
# set to the directory on your computer
SURVEY_FOLDER = "M:/CUG_depth"

survey = ppp.Survey(SURVEY_FOLDER)

Retrieve well CUG1:

[4]:
well_cug1 = survey.wells['CUG1']

Get a, b from well CUG1:

[5]:
a = well_cug1.params['nct']["a"]
b = well_cug1.params['nct']["b"]

Get n from well CUG1:

[6]:
n = well_cug1.params['n']

Retrieve seismic data:

[7]:
vel_cube = survey.seismics['velocity']
obp_cube = survey.seismics['obp_new']

View velocity section:

[8]:
fig_vel, ax_vel = plt.subplots()

im = vel_cube.plot(
    ppp.InlineIndex(8000), ax_vel, kind='img', cm='gist_rainbow')
fig_vel.colorbar(im)
fig_vel.set(figwidth=8)
[8]:
[None]
../_images/tutorial_eaton_seis_14_1.png

Pressure Prediction with Eaton method:

[9]:
eaton_cube = ppp.eaton_seis(
    "eaton_new", obp_cube, vel_cube, n=3,
    upper=survey.horizons['T16'], lower=survey.horizons['T20'])

eaton_seis function will automatically optimize the coefficients of Normal Compaction Trend, a and b.

View calculated pressure:

[10]:
from pygeopressure.basic.vawt import opendtect_seismic_colormap

fig_pres, ax_pres = plt.subplots()

im = eaton_cube.plot(
    ppp.InlineIndex(8000), ax_pres,
    kind='img', cm=opendtect_seismic_colormap())

fig_pres.colorbar(im)
fig_pres.set(figwidth=8)
[10]:
[None]
../_images/tutorial_eaton_seis_19_1.png