def evaluate(self, wat_sat):
wat_index = 0
krw_index = 1
Table = TableInterpolation()
if (wat_sat < self.SWOF[0][0]
or wat_sat > self.SWOF[len(self.SWOF) - 1][0]):
krw = Table.SCALExtraP(self.SWOF, wat_sat, wat_index, krw_index)
else:
krw = Table.LinearInterP(self.SWOF, wat_sat, wat_index, krw_index)
return krw
def evaluate(self, wat_sat):
wat_index = 0
pc_index = 3
Table = TableInterpolation()
if (wat_sat < self.SWOF[0][0]
or wat_sat > self.SWOF[len(self.SWOF) - 1][0]):
pc = Table.SCALExtraP(self.SWOF, wat_sat, wat_index, pc_index)
else:
pc = Table.LinearInterP(self.SWOF, wat_sat, wat_index, pc_index)
return pc
def evaluate(self, state):
pres = state[0]
pres_index = 0 # first column in the table - pressure
muo_index = 2 # third column in the table - viscosity
Table = TableInterpolation()
if (pres < self.table[0][0]
or pres > self.table[len(self.table) - 1][0]):
muo = Table.LinearExtraP(self.table, pres, pres_index, muo_index)
else:
muo = Table.LinearInterP(self.table, pres, pres_index, muo_index)
return muo
def evaluate(self, state):
pres = state[0]
pres_index = 0 # first column in the table - pressure
Bo_index = 1 # second column in the table - volume factor
Table = TableInterpolation()
if (pres < self.table[0][0]
or pres > self.table[len(self.table) - 1][0]):
Bo = Table.LinearExtraP(self.table, pres, pres_index, Bo_index)
else:
Bo = Table.LinearInterP(self.table, pres, pres_index, Bo_index)
return self.oil_dens / Bo
def evaluate(self, state):
water_saturation = dead_oil_water_saturation_evaluator(
self.pvdo, self.pvtw, self.surface_oil_dens, self.surface_wat_dens)
wat_sat = water_saturation.evaluate(state)
wat_index = 0
pc_index = 3
Table = TableInterpolation()
if (wat_sat < self.table[0][0]
or wat_sat > self.table[len(self.table) - 1][0]):
pc = Table.SCALExtraP(self.table, wat_sat, wat_index, pc_index)
else:
pc = Table.LinearInterP(self.table, wat_sat, wat_index, pc_index)
return pc