def test_oil_simulator(self):
sim = OilSimulator()
sim.addWell("OP1", seed=1)
sim.addBlock("6,6,6", seed=2)
for report_step in range(10):
sim.step(scale=1.0 / 10.0)
values = [
sim.fopr(),
sim.fopt(),
sim.fgpr(),
sim.fgpt(),
sim.fwpr(),
sim.fwpt(),
sim.fwct(),
sim.fgor(),
sim.opr("OP1"),
sim.gpr("OP1"),
sim.wpr("OP1"),
sim.gor("OP1"),
sim.wct("OP1"),
sim.bpr("6,6,6"),
]
self.assertAlmostEqual(values[0], values[8]) # fopr = opr:op1
self.assertAlmostEqual(values[2], values[9]) # fgpr = gpr:op1
self.assertAlmostEqual(values[4], values[10]) # fwpr = wpr:op1
self.assertAlmostEqual(sim.foip(), sim.ooip - sim.fopt())
self.assertAlmostEqual(sim.fgip(), sim.goip - sim.fgpt())
self.assertAlmostEqual(sim.fwip(), sim.woip - sim.fwpt())
self.assertAlmostEqualList(values, EXPECTED_VALUES[report_step])
def test_oil_simulator(self):
sim = OilSimulator()
sim.addWell("OP1", seed=1)
sim.addBlock("6,6,6", seed=2)
expected_values = [
[
0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0,
0.0, 1.0
],
[
0.3359771423145687, 0.3359771423145687, 0.25672494192349865,
0.25672494192349865, 0.010039005455891323,
0.010039005455891323, 0.029013112597713192, 0.7641143089523995,
0.3359771423145687, 0.25672494192349865, 0.010039005455891323,
0.7641143089523995, 0.029013112597713192, 0.8462347957619747
],
[
0.7252470407619624, 1.0612241830765312, 0.5175173529594699,
0.7742422948829686, 0.017973831236885583, 0.028012836692776905,
0.02418370085402209, 0.7135738912023045, 0.7252470407619624,
0.5175173529594699, 0.017973831236885583, 0.7135738912023045,
0.02418370085402209, 0.6888364145396828
],
[
0.7723163496234255, 1.8335405326999568, 0.5742386073607806,
1.3484809022437492, 0.11041673583737899, 0.1384295725301559,
0.12508507685507134, 0.7435277106858791, 0.7723163496234255,
0.5742386073607806, 0.11041673583737899, 0.7435277106858791,
0.12508507685507134, 0.6403046565762696
],
[
0.6038799675664164, 2.437420500266373, 0.6888868738548185,
2.037367776098568, 0.267892132439122, 0.4063217049692779,
0.3072960610203287, 1.140767885762087, 0.6038799675664164,
0.6888868738548185, 0.267892132439122, 1.140767885762087,
0.3072960610203287, 0.5205364945011657
],
[
0.23016535126253962, 2.6675858515289126, 0.721655666522216,
2.7590234426207836, 0.35552466124555465, 0.7618463662148325,
0.6070184801736589, 3.135379250454838, 0.23016535126253962,
0.721655666522216, 0.35552466124555465, 3.135379250454838,
0.6070184801736589, 0.4800677649914682
],
[
0.026293782934652718, 2.693879634463565, 0.7131990780527108,
3.4722225206734945, 0.6392372725163122, 1.4010836387311447,
0.8647254356377257, 7.131990780527107, 0.026293782934652718,
0.7131990780527108, 0.6392372725163122, 7.131990780527107,
0.8647254356377257, 0.3872974839053025
],
[
0.0, 2.693879634463565, 0.8676997908824122, 4.339922311555907,
0.8580356376693129, 2.2591192764004577, 0.8956197493411856,
8.676997908824122, 0.0, 0.8676997908824122, 0.8580356376693129,
8.676997908824122, 0.8956197493411856, 0.22557165737149715
],
[
0.10560669451549878, 2.799486328979064, 0.869082212788759,
5.209004524344666, 0.8903674796589355, 3.1494867560593933,
0.8939664328113363, 8.229423492288294, 0.10560669451549878,
0.869082212788759, 0.8903674796589355, 8.229423492288294,
0.8939664328113363, 0.1340241573819292
],
[
0.08615885630000791, 2.885645185279072, 0.44074890315982446,
5.64975342750449, 0.9425699260811738, 4.0920566821405675,
0.9040831722665535, 4.407489031598244, 0.08615885630000791,
0.44074890315982446, 0.9425699260811738, 4.407489031598244,
0.9040831722665535, 0.13404047971467026
]
]
for report_step in range(10):
sim.step(scale=1.0 / 10.0)
values = [
sim.fopr(),
sim.fopt(),
sim.fgpr(),
sim.fgpt(),
sim.fwpr(),
sim.fwpt(),
sim.fwct(),
sim.fgor(),
sim.opr("OP1"),
sim.gpr("OP1"),
sim.wpr("OP1"),
sim.gor("OP1"),
sim.wct("OP1"),
sim.bpr("6,6,6")
]
self.assertAlmostEqual(values[0], values[8]) # fopr = opr:op1
self.assertAlmostEqual(values[2], values[9]) # fgpr = gpr:op1
self.assertAlmostEqual(values[4], values[10]) # fwpr = wpr:op1
self.assertAlmostEqual(sim.foip(), sim.ooip - sim.fopt())
self.assertAlmostEqual(sim.fgip(), sim.goip - sim.fgpt())
self.assertAlmostEqual(sim.fwip(), sim.woip - sim.fwpt())
self.assertAlmostEqualList(values, expected_values[report_step])