def test_cell_containment(self):
grid_location = "local/ECLIPSE/faarikaal/faarikaal%d.EGRID"
well_location = "local/ECLIPSE/faarikaal/faarikaal%d.txt"
for i in range(1, 8):
grid_file = self.createTestPath(grid_location % i)
well_file = self.createTestPath(well_location % i)
grid = EclGrid(grid_file)
# Load well data
with open(well_file, "r") as f:
lines = [line.split() for line in f.readlines()]
points = [map(float, line[:3:]) for line in lines]
exp_cells = [tuple(map(int, line[3::])) for line in lines]
msg = "Expected point %s to be in cell %s, was in %s."
for point, exp_cell in zip(points, exp_cells):
reported_cell = grid.find_cell(*point)
self.assertEqual(
exp_cell,
reported_cell,
msg % (str(point), str(exp_cell), str(reported_cell)),
)
def test_cell_containment(self):
grid_location = "local/ECLIPSE/faarikaal/faarikaal%d.EGRID"
well_location = "local/ECLIPSE/faarikaal/faarikaal%d.txt"
for i in range(1, 8):
grid_file = self.createTestPath(grid_location % i)
well_file = self.createTestPath(well_location % i)
grid = EclGrid(grid_file)
# Load well data
with open(well_file, "r") as f:
lines = [line.split() for line in f.readlines()]
points = [map(float, line[:3:]) for line in lines]
exp_cells = [tuple(map(int, line[3::])) for line in lines]
msg = "Expected point %s to be in cell %s, was in %s."
for point, exp_cell in zip(points, exp_cells):
reported_cell = grid.find_cell(*point)
self.assertEqual(
exp_cell,
reported_cell,
msg % (str(point), str(exp_cell), str(reported_cell))
)
def test_update_simdata_from_rft_norne(norne_data):
"""Similar test as the reek version, but the Norne RFT file
does not contain saturations, and in libecl terms contains EclPLTCell
as opposed to EclRFTCell as in Reek"""
grid = EclGrid(ECL_BASE_NORNE + ".EGRID")
rft = EclRFTFile(ECL_BASE_NORNE + ".RFT")
rft_well_date = rft.get("C-3H", datetime.date(1999, 5, 4))
# A trajectory point for an active cell in Norne, picked
# from a line in gendata_rft_input_files/C-3H.txt
point = TrajectoryPoint(455752.59771598293, 7321015.949386452,
2785.78173828125, 2785.78173828125)
assert point.grid_ijk is None
assert point.pressure is None
assert point.swat is None
assert point.sgas is None
assert point.soil is None
point.set_ijk(
grid.find_cell(point.utm_x, point.utm_y, point.true_vertical_depth))
assert point.grid_ijk == (8, 12, 20) # Zero-indexed integers.
point.update_simdata_from_rft(rft_well_date)
# There is no saturation data in the Norne binary output, then these
# should be None
assert point.swat is None
assert point.sgas is None
assert point.soil is None
# Construct a Trajectory from the point
traj = Trajectory([])
traj.trajectory_points = [point] # (can't initialize from list of points)
dframe = traj.to_dataframe()
assert {"i", "j", "k", "pressure"}.issubset(set(dframe))
assert "swat" not in dframe
def test_update_simdata_from_rft_reek(reek_data):
"""Test data extraction from the binary Eclipse files
for a single well point using TrajectoryPoint.update_simdata_from_rft()"""
grid = EclGrid(ECL_BASE_REEK + ".EGRID")
rft = EclRFTFile(ECL_BASE_REEK + ".RFT")
rft_well_date = rft.get("OP_1", datetime.date(2000, 2, 1))
# A trajectory point for an active cell in reek:
point = TrajectoryPoint(462608.57, 5934210.96, 1624.38, 1624.38)
assert point.grid_ijk is None
assert point.pressure is None
assert point.swat is None
assert point.sgas is None
assert point.soil is None
point.set_ijk(
grid.find_cell(point.utm_x, point.utm_y, point.true_vertical_depth))
assert point.grid_ijk == (28, 27, 7)
point.update_simdata_from_rft(rft_well_date)
assert np.isclose(point.pressure, 304.37)
assert np.isclose(point.swat, 0.151044)
assert np.isclose(point.soil, 1 - 0.151044)
assert np.isclose(point.sgas, 0.0)
# Construct a Trajectory from the point
traj = Trajectory([])
traj.trajectory_points = [point] # (can't initialize from list of points)
dframe = traj.to_dataframe()
assert {"i", "j", "k", "pressure", "soil", "sgas",
"swat"}.issubset(set(dframe))
def test_update_simdata_outside_grid(tmpdir):
grid = EclGrid(ECL_BASE_REEK + ".EGRID")
rft = EclRFTFile(ECL_BASE_REEK + ".RFT")
rft_well_date = rft.get("OP_1", datetime.date(2000, 2, 1))
# A point outside the grid:
point = TrajectoryPoint(45000, 60000000, 1, 1)
point.set_ijk(
grid.find_cell(point.utm_x, point.utm_y, point.true_vertical_depth))
assert point.grid_ijk is None # There is no Exception raised by set_ijk()
point.update_simdata_from_rft(rft_well_date)
assert point.pressure is None # Since we are outside the grid.
# Construct a Trajectory from the point
traj = Trajectory([])
traj.trajectory_points = [point] # (can't initialize from list of points)
dframe = traj.to_dataframe()
assert not set(dframe).intersection({"i", "j", "k", "pressure", "swat"})
def test_update_simdata_outside_well(tmpdir):
grid = EclGrid(ECL_BASE_REEK + ".EGRID")
rft = EclRFTFile(ECL_BASE_REEK + ".RFT")
rft_well_date = rft.get("OP_1", datetime.date(2000, 2, 1))
# A point in the grid, but not related to the well
point = TrajectoryPoint(462825.55, 5934025.52, 1623.19, 1623.19)
point.set_ijk(
grid.find_cell(point.utm_x, point.utm_y, point.true_vertical_depth))
# NB: grid_ijk ints start at zero, ResInsight and ecl2df report this as (29, 29, 7)
assert point.grid_ijk == (28, 28, 6)
point.update_simdata_from_rft(rft_well_date)
assert point.pressure is None
assert point.swat is None
assert point.sgas is None
assert point.soil is None
# Construct a Trajectory from the point
traj = Trajectory([])
traj.trajectory_points = [point] # (can't initialize from list of points)
dframe = traj.to_dataframe()
assert {"i", "j", "k"}.issubset(set(dframe))
assert not set(dframe).intersection({"pressure", "swat", "soil", "sgas"})
