Two parallel noncomunicating zones in a reservoir

This one is only water flood (no low sal effect)

I'm injecting the seawater into a reservoir via a long horizontal well. The injectivity is better in some areas. I assign a larger injection rate to the 1D model for that zone, and a lower injection rate to the rest of the well, which is modelled again with a 1D model.



In [8]:

    
include("../FractionalFlow/FractionalFlow.jl")
using PyPlot, Dierckx, SetPyPlot
FF = FractionalFlow
setrcparam()









    



WARNING: replacing module FractionalFlow.



In [9]:

    
krw0_ww = 0.271636
krw0_ow = 0.226423
kro0_ww = 0.640774
kro0_ow = 0.53222
nw_ww = 3.74488
# nw_ow= 3.81109
# no_ww = 1.93454
nw_ow= 3.0
no_ww = 2.0
no_ow= 2.0844
sor_ww=0.333829
sor_ow=0.40325
sor_mw=0.1
swc_ww=0.212
swc_ow=0.212

u_inj_fast = 0.3e-5
u_inj_slow = 1.0e-5 
pv_inj_fast = 2.0
pv_inj_slow = u_inj_slow/u_inj_fast*pv_inj_fast
A_fast = 0.7 # normalized to a total injection area of 1.0
A_slow = 1-A_fast # normalized to a total injection area of 1.0
eq_const = 0.3









    Out[9]:





0.3



In [10]:

    
fluids_hs = FF.oil_water_fluids(mu_water=1.1e-3, mu_oil=2e-3)
fluids_ls = FF.oil_water_fluids(mu_water=1e-3, mu_oil=2e-3)
rel_perms_hs = FF.oil_water_rel_perms(krw0=krw0_ow, kro0=kro0_ow, 
        swc=swc_ow, sor=sor_ow, nw=nw_ow, no = no_ow)
rel_perms_ls = FF.oil_water_rel_perms(krw0=krw0_ww, kro0=kro0_ww, 
        swc=swc_ww, sor=sor_ww, nw=nw_ww, no = no_ww)
core_flood_fast = FF.core_flooding(u_inj=u_inj_fast, pv_inject=pv_inj_fast, p_back=1e5, sw_init=swc_ow, sw_inj=1.0, rel_perms=rel_perms_hs)
core_flood_slow = FF.core_flooding(u_inj=u_inj_slow, pv_inject=pv_inj_slow, p_back=1e5, sw_init=swc_ow, sw_inj=1.0, rel_perms=rel_perms_hs)

core_props = FF.core_properties()
ls_res_fast = FF.water_flood(core_props, fluids_ls, rel_perms_hs, core_flood_fast)
ls_res_slow = FF.water_flood(core_props, fluids_ls, rel_perms_hs, core_flood_slow)

FF.visualize(ls_res_fast)
FF.visualize(ls_res_slow)









    



┌ Info: difficulty finding the shock front saturation range!
└ @ Main.FractionalFlow /home/ali/projects/peteng/analytical/FractionalFlow/FractionalFlow.jl:401
┌ Info: difficulty finding the shock front saturation range!
└ @ Main.FractionalFlow /home/ali/projects/peteng/analytical/FractionalFlow/FractionalFlow.jl:401






    



[0.0, 8.57888e-8, 1.83728e-7, 2.8803e-7, 3.97402e-7, 5.11244e-7, 6.29219e-7, 7.51116e-7, 8.76802e-7, 1.00619e-6, 1.13921e-6, 1.27583e-6, 1.41604e-6, 1.55981e-6, 1.70714e-6, 1.85806e-6, 2.01255e-6, 2.17064e-6, 2.33236e-6, 2.49771e-6, 2.66674e-6, 2.83946e-6, 3.0159e-6, 3.19611e-6, 3.3801e-6, 3.56791e-6, 3.75959e-6, 3.95516e-6, 4.15466e-6, 4.35813e-6, 4.56561e-6, 4.77713e-6, 4.99274e-6, 5.21247e-6, 5.43636e-6, 5.66446e-6, 5.89679e-6, 6.13341e-6, 6.37434e-6, 6.61964e-6, 6.86933e-6, 7.12346e-6, 7.38206e-6, 7.64517e-6, 7.91283e-6, 8.18507e-6, 8.46194e-6, 8.74346e-6, 9.02967e-6, 9.32061e-6, 9.6163e-6, 9.91678e-6, 1.02221e-5, 1.05322e-5, 1.08473e-5, 1.11672e-5, 1.14921e-5, 1.18219e-5, 1.21567e-5, 1.24965e-5, 1.28413e-5, 1.31911e-5, 1.3546e-5, 1.3906e-5, 1.4271e-5, 1.46411e-5, 1.50164e-5, 1.53967e-5, 1.57821e-5, 1.61727e-5, 1.65683e-5, 1.69691e-5, 1.73749e-5, 1.77859e-5, 1.82019e-5, 1.86229e-5, 1.9049e-5, 1.94801e-5, 1.99162e-5, 2.03573e-5, 2.08032e-5, 2.12541e-5, 2.17099e-5, 2.21704e-5, 2.26357e-5, 2.31057e-5, 2.35804e-5, 2.40597e-5, 2.45435e-5, 2.50318e-5, 2.55245e-5, 2.60215e-5, 2.65227e-5, 2.70281e-5, 2.75376e-5, 2.8051e-5, 2.85683e-5, 2.90893e-5, 2.9614e-5, 3.01422e-5, 3.01422e-5, 6.02845e-5]
[0.59675, 0.595891, 0.595032, 0.594172, 0.593313, 0.592454, 0.591595, 0.590735, 0.589876, 0.589017, 0.588158, 0.587298, 0.586439, 0.58558, 0.584721, 0.583862, 0.583002, 0.582143, 0.581284, 0.580425, 0.579565, 0.578706, 0.577847, 0.576988, 0.576128, 0.575269, 0.57441, 0.573551, 0.572692, 0.571832, 0.570973, 0.570114, 0.569255, 0.568395, 0.567536, 0.566677, 0.565818, 0.564958, 0.564099, 0.56324, 0.562381, 0.561522, 0.560662, 0.559803, 0.558944, 0.558085, 0.557225, 0.556366, 0.555507, 0.554648, 0.553788, 0.552929, 0.55207, 0.551211, 0.550352, 0.549492, 0.548633, 0.547774, 0.546915, 0.546055, 0.545196, 0.544337, 0.543478, 0.542618, 0.541759, 0.5409, 0.540041, 0.539182, 0.538322, 0.537463, 0.536604, 0.535745, 0.534885, 0.534026, 0.533167, 0.532308, 0.531448, 0.530589, 0.52973, 0.528871, 0.528012, 0.527152, 0.526293, 0.525434, 0.524575, 0.523715, 0.522856, 0.521997, 0.521138, 0.520278, 0.519419, 0.51856, 0.517701, 0.516842, 0.515982, 0.515123, 0.514264, 0.513405, 0.512545, 0.511686, 0.212, 0.212]
[0.0, 5.0e-6, 1.0e-5, 1.5e-5, 2.0e-5, 2.5e-5, 3.0e-5, 3.5e-5, 4.0e-5, 4.5e-5, 5.0e-5, 5.5e-5, 6.0e-5, 6.5e-5, 7.0e-5, 7.5e-5, 8.0e-5, 8.5e-5, 9.0e-5, 9.5e-5, 0.0001, 0.000105, 0.00011, 0.000115, 0.00012, 0.000125, 0.00013, 0.000135, 0.00014, 0.000145, 0.00015, 0.000155, 0.00016, 0.000165, 0.00017, 0.000175, 0.00018, 0.000185, 0.00019, 0.000195, 0.0002, 0.000205, 0.00021, 0.000215, 0.00022, 0.000225, 0.00023, 0.000235, 0.00024, 0.000245, 0.00025, 0.000255, 0.00026, 0.000265, 0.00027, 0.000275, 0.00028, 0.000285, 0.00029, 0.000295, 0.0003, 0.000305, 0.00031, 0.000315, 0.00032, 0.000325, 0.00033, 0.000335, 0.00034, 0.000345, 0.00035, 0.000355, 0.00036, 0.000365, 0.00037, 0.000375, 0.00038, 0.000385, 0.00039, 0.000395, 0.0004, 0.000405, 0.00041, 0.000415, 0.00042, 0.000425, 0.00043, 0.000435, 0.00044, 0.000445, 0.00045, 0.000455, 0.00046, 0.000465, 0.00047, 0.000475, 0.00048, 0.000485, 0.00049, 0.000495, 0.0005, 0.000505, 0.00051, 0.000515, 0.00052, 0.000525, 0.00053, 0.000535, 0.00054, 0.000545, 0.00055, 0.000555, 0.00056, 0.000565, 0.00057, 0.000575, 0.00058, 0.000585, 0.00059, 0.000595, 0.0006, 0.000605, 0.00061, 0.000615, 0.00062, 0.000625, 0.00063, 0.000635, 0.00064, 0.000645, 0.00065, 0.000655, 0.00066, 0.000665, 0.00067, 0.000675, 0.00068, 0.000685, 0.00069, 0.000695, 0.0007, 0.000705, 0.00071, 0.000715, 0.00072, 0.000725, 0.00073, 0.000735, 0.00074, 0.000745, 0.00075, 0.000755, 0.00076, 0.000765, 0.00077, 0.000775, 0.00078, 0.000785, 0.00079, 0.000795, 0.0008, 0.000805, 0.00081, 0.000815, 0.00082, 0.000825, 0.00083, 0.000835, 0.00084, 0.000845, 0.00085, 0.000855, 0.00086, 0.000865, 0.00087, 0.000875, 0.00088, 0.000885, 0.00089, 0.000895, 0.0009, 0.000905, 0.00091, 0.000915, 0.00092, 0.000925, 0.00093, 0.000935, 0.00094, 0.000945, 0.00095, 0.000955, 0.00096, 0.000965, 0.00097, 0.000975, 0.00098, 0.000985, 0.00099, 0.000995]
[0.59675, 0.569226, 0.552695, 0.540078, 0.529567, 0.520334, 0.511918, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212]






    












    












    












    












    












    












    












    












    












    












    



[0.0, 2.85963e-7, 6.12427e-7, 9.601e-7, 1.32467e-6, 1.70415e-6, 2.0974e-6, 2.50372e-6, 2.92267e-6, 3.35395e-6, 3.79736e-6, 4.25277e-6, 4.72012e-6, 5.19936e-6, 5.69048e-6, 6.19352e-6, 6.7085e-6, 7.23548e-6, 7.77453e-6, 8.32571e-6, 8.88913e-6, 9.46486e-6, 1.0053e-5, 1.06537e-5, 1.1267e-5, 1.1893e-5, 1.2532e-5, 1.31839e-5, 1.38489e-5, 1.45271e-5, 1.52187e-5, 1.59238e-5, 1.66425e-5, 1.73749e-5, 1.81212e-5, 1.88815e-5, 1.9656e-5, 2.04447e-5, 2.12478e-5, 2.20655e-5, 2.28978e-5, 2.37449e-5, 2.46069e-5, 2.54839e-5, 2.63761e-5, 2.72836e-5, 2.82065e-5, 2.91449e-5, 3.00989e-5, 3.10687e-5, 3.20543e-5, 3.30559e-5, 3.40736e-5, 3.51074e-5, 3.61575e-5, 3.7224e-5, 3.83068e-5, 3.94062e-5, 4.05222e-5, 4.16548e-5, 4.28042e-5, 4.39704e-5, 4.51534e-5, 4.63532e-5, 4.757e-5, 4.88038e-5, 5.00546e-5, 5.13223e-5, 5.26071e-5, 5.39089e-5, 5.52278e-5, 5.65636e-5, 5.79164e-5, 5.92862e-5, 6.06728e-5, 6.20764e-5, 6.34967e-5, 6.49337e-5, 6.63873e-5, 6.78575e-5, 6.93442e-5, 7.08471e-5, 7.23662e-5, 7.39014e-5, 7.54524e-5, 7.70191e-5, 7.86014e-5, 8.0199e-5, 8.18117e-5, 8.34394e-5, 8.50816e-5, 8.67383e-5, 8.84091e-5, 9.00937e-5, 9.17919e-5, 9.35033e-5, 9.52275e-5, 9.69644e-5, 9.87133e-5, 0.000100474, 0.000100474, 0.000200948]
[0.59675, 0.595891, 0.595032, 0.594172, 0.593313, 0.592454, 0.591595, 0.590735, 0.589876, 0.589017, 0.588158, 0.587298, 0.586439, 0.58558, 0.584721, 0.583862, 0.583002, 0.582143, 0.581284, 0.580425, 0.579565, 0.578706, 0.577847, 0.576988, 0.576128, 0.575269, 0.57441, 0.573551, 0.572692, 0.571832, 0.570973, 0.570114, 0.569255, 0.568395, 0.567536, 0.566677, 0.565818, 0.564958, 0.564099, 0.56324, 0.562381, 0.561522, 0.560662, 0.559803, 0.558944, 0.558085, 0.557225, 0.556366, 0.555507, 0.554648, 0.553788, 0.552929, 0.55207, 0.551211, 0.550352, 0.549492, 0.548633, 0.547774, 0.546915, 0.546055, 0.545196, 0.544337, 0.543478, 0.542618, 0.541759, 0.5409, 0.540041, 0.539182, 0.538322, 0.537463, 0.536604, 0.535745, 0.534885, 0.534026, 0.533167, 0.532308, 0.531448, 0.530589, 0.52973, 0.528871, 0.528012, 0.527152, 0.526293, 0.525434, 0.524575, 0.523715, 0.522856, 0.521997, 0.521138, 0.520278, 0.519419, 0.51856, 0.517701, 0.516842, 0.515982, 0.515123, 0.514264, 0.513405, 0.512545, 0.511686, 0.212, 0.212]
[0.0, 5.0e-6, 1.0e-5, 1.5e-5, 2.0e-5, 2.5e-5, 3.0e-5, 3.5e-5, 4.0e-5, 4.5e-5, 5.0e-5, 5.5e-5, 6.0e-5, 6.5e-5, 7.0e-5, 7.5e-5, 8.0e-5, 8.5e-5, 9.0e-5, 9.5e-5, 0.0001, 0.000105, 0.00011, 0.000115, 0.00012, 0.000125, 0.00013, 0.000135, 0.00014, 0.000145, 0.00015, 0.000155, 0.00016, 0.000165, 0.00017, 0.000175, 0.00018, 0.000185, 0.00019, 0.000195, 0.0002, 0.000205, 0.00021, 0.000215, 0.00022, 0.000225, 0.00023, 0.000235, 0.00024, 0.000245, 0.00025, 0.000255, 0.00026, 0.000265, 0.00027, 0.000275, 0.00028, 0.000285, 0.00029, 0.000295, 0.0003, 0.000305, 0.00031, 0.000315, 0.00032, 0.000325, 0.00033, 0.000335, 0.00034, 0.000345, 0.00035, 0.000355, 0.00036, 0.000365, 0.00037, 0.000375, 0.00038, 0.000385, 0.00039, 0.000395, 0.0004, 0.000405, 0.00041, 0.000415, 0.00042, 0.000425, 0.00043, 0.000435, 0.00044, 0.000445, 0.00045, 0.000455, 0.00046, 0.000465, 0.00047, 0.000475, 0.00048, 0.000485, 0.00049, 0.000495, 0.0005, 0.000505, 0.00051, 0.000515, 0.00052, 0.000525, 0.00053, 0.000535, 0.00054, 0.000545, 0.00055, 0.000555, 0.00056, 0.000565, 0.00057, 0.000575, 0.00058, 0.000585, 0.00059, 0.000595, 0.0006, 0.000605, 0.00061, 0.000615, 0.00062, 0.000625, 0.00063, 0.000635, 0.00064, 0.000645, 0.00065, 0.000655, 0.00066, 0.000665, 0.00067, 0.000675, 0.00068, 0.000685, 0.00069, 0.000695, 0.0007, 0.000705, 0.00071, 0.000715, 0.00072, 0.000725, 0.00073, 0.000735, 0.00074, 0.000745, 0.00075, 0.000755, 0.00076, 0.000765, 0.00077, 0.000775, 0.00078, 0.000785, 0.00079, 0.000795, 0.0008, 0.000805, 0.00081, 0.000815, 0.00082, 0.000825, 0.00083, 0.000835, 0.00084, 0.000845, 0.00085, 0.000855, 0.00086, 0.000865, 0.00087, 0.000875, 0.00088, 0.000885, 0.00089, 0.000895, 0.0009, 0.000905, 0.00091, 0.000915, 0.00092, 0.000925, 0.00093, 0.000935, 0.00094, 0.000945, 0.00095, 0.000955, 0.00096, 0.000965, 0.00097, 0.000975, 0.00098, 0.000985, 0.00099, 0.000995]
[0.59675, 0.585937, 0.577924, 0.571245, 0.565443, 0.560277, 0.555596, 0.5513, 0.547317, 0.543589, 0.540078, 0.536752, 0.533584, 0.53055, 0.527637, 0.524825, 0.522104, 0.519462, 0.516889, 0.514377, 0.511918, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212, 0.212]






    



┌ Info: difficulty finding the shock front saturation range!
└ @ Main.FractionalFlow /home/ali/projects/peteng/analytical/FractionalFlow/FractionalFlow.jl:401
┌ Info: difficulty finding the shock front saturation range!
└ @ Main.FractionalFlow /home/ali/projects/peteng/analytical/FractionalFlow/FractionalFlow.jl:401






    Out[10]:





PyObject <matplotlib.legend.Legend object at 0x7f4cd7a267f0>



In [11]:

    
plot(ls_res_fast.water_cut_time[:,1], ls_res_fast.water_cut_time[:,2], 
    ls_res_slow.water_cut_time[:,1], ls_res_slow.water_cut_time[:,2])
legend(["Fast", "Slow"])









    












    Out[11]:





PyObject <matplotlib.legend.Legend object at 0x7f4cd7550b70>

The above figure shows the water cut in each zone. I need to add them together based on the total flow rate in each zone and plot the final value versus time. I need to know the total injection rate, the total surface area of each zone, and the water cut versus time. Then I calculate the new water cut versus time by calculating the total flow of oil and normalizing it.



In [12]:

    
t_fast = ls_res_fast.water_cut_time[:,1]
WC_fast = ls_res_fast.water_cut_time[:,2]
t_slow = ls_res_slow.water_cut_time[:,1]
WC_slow = ls_res_slow.water_cut_time[:,2]
WC_f =  Spline1D(t_fast, WC_fast, k=1, bc="nearest")
WC_s =  Spline1D(t_slow, WC_slow, k=1, bc="nearest")









    Out[12]:





Spline1D(knots=[0.0,150.754 … 29849.2,30000.0] (200 elements), k=1, extrapolation="nearest", residual=0.0)



In [13]:

    
t = range(0, 20000, length=20000)
WC = (u_inj_fast*A_fast.*WC_f(t)+u_inj_slow*A_slow.*WC_s(t))./(u_inj_fast*A_fast+u_inj_slow*A_slow)









    Out[13]:





20000-element Array{Float64,1}:
 0.0               
 0.0               
 0.0               
 0.0               
 0.0               
 0.0               
 0.0               
 0.0               
 0.0               
 0.0               
 0.0               
 0.0               
 0.0               
 ⋮                 
 0.9943416003180039
 0.9943420593422865
 0.994342518366569 
 0.9943429773908519
 0.9943434364151345
 0.9943438954394171
 0.9943443544637   
 0.9943448134879825
 0.9943452725122651
 0.9943457315365479
 0.9943461905608305
 0.9943466495851132



In [14]:

    
total_pv = pv_inj_fast*A_fast+pv_inj_slow*(1-A_fast)



plot(ls_res_fast.water_cut_time[:,1]./t_fast[end].*total_pv, ls_res_fast.water_cut_time[:,2], 
    ls_res_slow.water_cut_time[:,1]./t_fast[end].*total_pv, ls_res_slow.water_cut_time[:,2])

plot(t./t_fast[end].*total_pv, WC, linewidth = 3, "--")

legend(["Late breakthrough", "Early breakthrough", "Overall"])

xlabel("Injected PV [-]")
ylabel("Water cut [-]")

savefig("results/water_cut_WF_zones.png")