perm_powers_truths = np.log10(np.array([
1.00000000e-15, 1.00000000e-15, 1.00000000e-15,
2.00000000e-15, 2.00000000e-15, 8.00000000e-15,
5.00000000e-16, 5.00000000e-16, 5.00000000e-16,
1.00000000e-16, 1.00000000e-16, 1.00000000e-16,
5.00000000e-14, 5.00000000e-14, 5.00000000e-16,
5.00000000e-14, 5.00000000e-14, 5.00000000e-16,
5.00000000e-14, 1.00000000e-14, 3.00000000e-14,
5.00000000e-15, 5.00000000e-15, 5.00000000e-15,
5.00000000e-14, 5.00000000e-14, 5.00000000e-16,
5.00000000e-14, 5.00000000e-14, 5.00000000e-16]))
#---coarse process model
process_model_coarse = GC.GeoModel(name='test_process_model_kerinci',
datfile_name='input-files/kerinci/coarse-model/Keriv0_027',
incon_name='input-files/kerinci/coarse-model/Keriv0_027',
geom_name='input-files/kerinci/coarse-model/gKerinci_v0')
list_of_obs_wells = ['LP002','LP001','KRDB1']
process_model_coarse.rename_wells_as_obs(list_of_obs_wells)
process_model_coarse.set_rock_permeabilities(perm_powers=perm_powers_truths)
# ---Layered
# perm_powers_truths = np.log10(np.array([
# 1.00000000e-15, 1.00000000e-15,
# 2.00000000e-15, 8.00000000e-15,
# 5.00000000e-16, 5.00000000e-16,
# 1.00000000e-16, 1.00000000e-16,
# 5.00000000e-14, 5.00000000e-16,
# 5.00000000e-14, 5.00000000e-16,
5.00000000e-14, 1.00000000e-14, 1.00000000e-16, 1.00000000e-16,
1.00000000e-15, 2.00000000e-14, 2.50000000e-14, 1.00000000e-14,
5.00000000e-16, 5.00000000e-16, 5.00000000e-16, 1.00000000e-14
]))
#---coarse process model
#process_model_coarse = GC.GeoModel(name='test_process_model_coarse',
# datfile_name='input-files/elvar-new/coarse-model/2DC002',
# incon_name='input-files/elvar-new/coarse-model/2DC002_IC',
# geom_name = 'input-files/elvar-new/coarse-model/g2coarse',
# islayered=True)
#---medium process model
process_model_medium = GC.GeoModel(
name='test_process_model_medium',
datfile_name='input-files/elvar-new-tighter-cap/medium-model/2DM002',
incon_name='input-files/elvar-new-tighter-cap/medium-model/2DM002_IC',
geom_name='input-files/elvar-new-tighter-cap/medium-model/g2medium',
islayered=True)
#---fine process model
process_model_fine = GC.GeoModel(
name='test_process_model_fine',
datfile_name='input-files/elvar-new-tighter-cap/fine-model/2DF002',
incon_name='input-files/elvar-new-tighter-cap/fine-model/2DF002_IC',
geom_name='input-files/elvar-new-tighter-cap/fine-model/g2fine',
islayered=True)
#---fine synthetic model
synthetic_model_fine = GC.GeoModel(
name='test_synthetic_model_fine',
datfile_name='input-files/elvar-new-tighter-cap/fine-model/2DF002',
perm_powers_truths = np.log10(np.array([
1.00000000e-15, 1.00000000e-15, 1.00000000e-15,
2.00000000e-15, 2.00000000e-15, 8.00000000e-15,
5.00000000e-16, 5.00000000e-16, 5.00000000e-16,
1.00000000e-16, 1.00000000e-16, 1.00000000e-16,
5.00000000e-14, 5.00000000e-14, 5.00000000e-16,
5.00000000e-14, 5.00000000e-14, 5.00000000e-16,
5.00000000e-14, 1.00000000e-14, 3.00000000e-14,
5.00000000e-15, 5.00000000e-15, 5.00000000e-15,
5.00000000e-14, 5.00000000e-14, 5.00000000e-16,
5.00000000e-14, 5.00000000e-14, 5.00000000e-16]))
#---coarse process model
process_model_coarse = GC.GeoModel(name='test_process_model_kerinci',
datfile_name='input-files/kerinci/coarse-model/Keriv0_027',
incon_name='input-files/kerinci/coarse-model/Keriv0_027',
geom_name='input-files/kerinci/coarse-model/gKerinci_v0')
list_of_obs_wells = ['LP002','LP001','KRDB1']
process_model_coarse.rename_wells_as_obs(list_of_obs_wells)
process_model_coarse.set_rock_permeabilities(perm_powers=perm_powers_truths)
#list_of_obs_wells = ['LP002'] #['LP002','LP001','KRDB1']'KRDB1']
#---fine process model
process_model_fine = GC.GeoModel(name='test_process_model_kerinci_fine',
datfile_name='input-files/kerinci/fine-model/Keriv1_027',
incon_name='input-files/kerinci/fine-model/Keriv1_027',
geom_name='input-files/kerinci/fine-model/gKerinci_v1')