Complete annotated TOML configuration#

The following blocks reproduce the annotated input.toml configuration used throughout the original documentation. They are presented in file order. Copy them into one TOML file, preserving the order where convenient.

For explanations, units, model relationships, and figures, use Integrated models and geometry, Boundary conditions and projection, Rock, saturation, hysteresis, and salinity, and Wells, schedules, and TUNING.

View input.toml | View raw input.toml

 1# Set mpirun, the full path to the flow executable, and simulator flags (except --output-dir)
 2flow = "flow --relaxed-max-pv-fraction=0 --enable-opm-rst-file=true --newton-min-iterations=1"
 3
 4# Set the model parameters
 5regional_dims = [45000,15000,81] # Regional aquifer length, width, and depth [m]
 6regional_x_n = [3,5,5] # Variable array of x-refinement (Regional)
 7regional_y_n = [5,5,5] # Variable array of y-refinement (Regional)
 8regional_z_n = [1,1,1] # Variable array of z-refinement (Regional)
 9reference_x_n = [1,1,1,5,5,15,15,15,15,15,5,5,5,5,5] # Variable array of x-refinement (Reference)
10reference_y_n = [25,25,25] # Variable array of y-refinement (Reference)
11reference_z_n = [6,5,5,5,3,3,3,3,3] # Variable array of z-refinement (Reference)
12site_location = [18000,5000,0,27000,10000,81] # Site xi, yi, zi, xf, yf, and zf box positions [m]
13fault_regional = [10000,11000,0.01,10,22.5] # Regional fault x, and y positions [m], x and y multipliers for the trans, and height of the fault jump [m]
14fault_site = [[21583,5710],[24081,8233],[0,0]] # Site fault x, and y positions [m] (initial and final) and x and y multipliers for the trans
15thickness = [9,9,9,9,9,9,9,9,9] # Thickness of the layers [m]
16pressure = 2E2  # Pressure on the reservoir top [Bar]
17temperature = [60,50] # Top and bottom temperatures [C]
18rock_comp = 6.11423e-5 # Rock compressibility [1/Bar]
19sensor_coords = [20000,8000,0] # Sensor position x, y, and z to assess the error over time w.r.t the reference solution [m]
20z_xy = "(20-20*mt.sin((2*mt.pi*(x+y)/10000)))" # The function for the reservoir surface
21hysteresis = ["Killough", "Both"] # Activate hysteresis: 1) Killough or Carlson and 2) Both, Pc, or Kr (HYSTOPT in EHYSTR keyword, see the OPM manual)
22salinity = 2.92 # Add salinity (value [1e-3 kg-M/kg])
23
24# Set the boundary conditions
25# Use open/closed/porv for the Regional aquifer (if porv, enter the bottom, right, top, and left values (e.g, porv 1e8 1e7 1e6 1e5))
26regional_bctype = ["open"]
27# Use open/closed/porv/porvproj/flux/pres/pres2p/wells for the BC site (if porv; bottom, right, top, and left values (e.g, porv 1e4 1e3 1e2 1e1))
28# For pres/flux BC's, add 'interp' to use linear interpolation in time
29site_bctype = ["open"]
30
31# Set the saturation functions
32krw = "krw * ((sw - swi) / (1.0 - sni -swi)) ** nkrw"        # Wetting rel perm saturation function [-]
33krn = "krn * ((1.0 - sw - sni) / (1.0 - sni - swi)) ** nkrn" # Non-wetting rel perm saturation function [-]
34pcap = "pen * ((sw - swi) / (1.0 - swi)) ** (-(1.0 / npen))" # Capillary pressure saturation function [Bar]
35
36# Properties sat functions: 1) swi [-], 2) sni [-], 3) krw [-], 4) krn [-], 5) pen [Bar], 6) nkrw [-], 7) nkrn [-],
37# 8) npen [-], 9) threshold cP evaluation, and 10) npoints [-] (entry per layer, if hysteresis, additional entries per layer)
38safu = [[0.12,0.10,1,1,3060e-5,2,2,2,1e-4,1000],
39[0.12,0.10,1.0,1.0,3870.63e-5,2,2,2,1e-4,1000],
40[0.12,0.10,1.0,1.0,3060.00e-5,2,2,2,1e-4,1000],
41[0.12,0.10,1.0,1.0,3870.63e-5,2,2,2,1e-4,1000],
42[0.32,0.10,1.0,1.0, 193531e-5,2,2,2,1e-4,1000],
43[0.14,0.10,1.0,1.0,8654.99e-5,2,2,2,1e-4,1000],
44[0.12,0.10,1.0,1.0,6120.00e-5,2,2,2,1e-4,1000],
45[0.14,0.10,1.0,1.0,8654.99e-5,2,2,2,1e-4,1000],
46[0.12,0.10,1.0,1.0,6120.00e-5,2,2,2,1e-4,1000],
47[0.22,0.30,1.0,1.0,3060.00e-5,2,4,2,1e-4,1000],
48[0.12,0.30,1.0,1.0,3870.63e-5,2,4,2,1e-4,1000],
49[0.12,0.30,1.0,1.0,3060.00e-5,2,4,2,1e-4,1000],
50[0.12,0.30,1.0,1.0,3870.63e-5,2,4,2,1e-4,1000],
51[0.32,0.30,1.0,1.0, 193531e-5,2,4,2,1e-4,1000],
52[0.14,0.30,1.0,1.0,8654.99e-5,2,4,2,1e-4,1000],
53[0.12,0.30,1.0,1.0,6120.00e-5,2,4,2,1e-4,1000],
54[0.14,0.30,1.0,1.0,8654.99e-5,2,4,2,1e-4,1000],
55[0.12,0.30,1.0,1.0,6120.00e-5,2,4,2,1e-4,1000]]
56
57# Properties rock: 1) Kxy [mD], 2) Kz [mD], and 3) phi [-] (entry per layer)
58rock = [[1013.25,101.325,0.25],
59[506.625,50.6625,0.20],
60[1013.25,101.325,0.25],
61[506.625,50.6625,0.20],
62[0.10132,0.01013,0.10],
63[101.324,10.1324,0.20],
64[202.650,20.2650,0.20],
65[101.324,10.1324,0.20],
66[202.650,20.2650,0.20]]
67
68# Wells position: 1) x, 2) y, 3) zi, and 4) zf positions [m] (entry per well)
69well_coords = [[21180,7068,0,81],[24200,7800,15,65],[21718,7122,45,81],[14518,11377,0,50],[31679,8883,0,30],[28477,2732,0,81]]
70
71# Define the injection values (entry per change in the schedule):
72# 1) injection time [d], 2) time step size to write results regional [d], 3) time step size to write results site/reference [d]
73# 1) fluid (0 wetting, 1 non-wetting) well 0, 2) injection rates [kg/day] well 0, 3) fluid ... well n, injection, ...well n (as many as num of wells)
74# if 'wells' for site_bctype, then 1) bottom, 2) right, 3) top, and 4) left values (0(prod)/1(inj), pressure [Bar]))
75# if --enable-tuning=1, then 1) for TUNING values as described in the OPM manual
76inj = [[[365,73,73],[1,3e5,1,3e5,1,3e5,1,5e6,1,5e6,0,1e7]],[[365,73,73],[1,3e5,1,3e5,1,3e5,1,5e6,1,0,0,1e7]]]