def test_create(setup_spec):
spec = setup_spec
vel = "TEST.segy"
create_segy_file(vel, spec)
grid = EclGridGenerator.createRectangular(dims=(2, 2, 2), dV=(100, 100, 100))
res_surface = OTSResSurface(grid=grid, above=0)
OTSVelSurface(res_surface, vel)
def test_z3d(setup_spec):
spec = setup_spec
vel = "TEST.segy"
int_val = [50, 150]
create_segy_file(vel, spec, xl=int_val, il=int_val, cdp_x=int_val, cdp_y=int_val)
grid = EclGridGenerator.createRectangular(dims=(2, 2, 2), dV=(100, 100, 100))
res_surface = OTSResSurface(grid=grid, above=0)
ots_s = OTSVelSurface(res_surface, vel)
assert (4, 2) == ots_s.z3d.shape
def __init__(
self,
eclbase,
mapaxes,
seabed,
youngs,
poisson,
rfactor,
convention,
above,
velocity_model,
):
"""
The OTS class manages the information required to calculate
overburden timeshift.
The constructor will look for the Eclipse files INIT, EGRID
and UNRST based on the input case, if some of the files are
missing an exception will be raised. It will then instantiate
a EclSubsidence object will be used to manage the rest of the
overburden timeshift calculations.
"""
case = os.path.splitext(eclbase)[0]
self._init_file = EclFile("%s.INIT" % case)
self._rst_file = EclFile("%s.UNRST" % case)
self._grid = EclGrid("%s.EGRID" % case, apply_mapaxes=mapaxes)
self.subsidence = EclSubsidence(self._grid, self._init_file)
self._seabed = seabed
self._youngs_modulus = youngs * 1e9
self._poisson_ratio = poisson
self._r_factor = rfactor
self._convention = convention
self._surface = res_surface = OTSResSurface(grid=self._grid,
above=above)
if velocity_model is not None:
self._surface = OTSVelSurface(res_surface=res_surface,
vcube=velocity_model)
self._restart_views = {}