In [1]:
import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline
import welly
welly.__version__
Out[1]:
In [2]:
from welly import Well
In [3]:
w = Well.from_las('P-129_out.LAS')
The make_synthetic() method assumes that the DT log is called 'DT' and has units of µs/m, and the RHOB log is called 'RHOB' and has units of kg/m3. Our logs have the right names, but the wrong units, so let's adjust them:
In [4]:
w.data['DT'] = w.data['DT'].extrapolate() / 0.3048
w.data['RHOB'] = 1000 * w.data['RHOB']
In [5]:
w.make_synthetic(srd=100)
In [6]:
s = w.data['Synthetic']
s.plot()
In [ ]: