SIMULATION¶

Defines the flow model to simulate. See below an example to define a GAS_WATER case running in isothermal mode:

SIMULATION
PROCESS_MODELS
SUBSURFACE_FLOW Flow
MODE GAS_WATER
OPTIONS
ISOTHERMAL
RESERVOIR_DEFAULTS
/                    ! end options
/                      ! end subsurface_flow
/                        ! end process_models
/                          ! end simulation


Here is an unlikely example showing every one of the keywords which will be described below:

SIMULATION
PREAMBLE
PROCESS_MODELS
SUBSURFACE_FLOW Flow
MODE GAS_WATER
OPTIONS
ISOTHERMAL
RESERVOIR_DEFAULTS
NOGASSOL
ANALYTICAL_JACOBIAN
HYSTERESIS
HYSTERESIS_PC 0.1
TL_OMEGA 0.5
WATER_VISCOSITY_NOSALT
/                        ! end options
/                          ! end subsurface_flow
/                            ! end process_models
CHECKPOINT
[checkpoint block]
/                            ! end checkpoint
restart test-23.0000y.h5
ERESTART base 3000 d
/                              ! end simulation


Below is a list of all options that can be specified in the SIMULATION block:

MODE¶

The MODE keyword selects the flow model for the simulation, and can take the following arguments: COMP, COMP3, COMP4, BLACK_OIL, GAS_WATER, TOIL_IMS, TODD_LONGSTAFF, SOLVENT_TL. MODE must be entered by the user as it cannot be defaulted.

For COMP mode, one must enter:

MODE COMP [Nc] EOS


where Nc is the total number of components (including water), for example:

MODE COMP 2 EOS


for two components, one of which is water.

COMP3 and COMP4 are two different modelling options of the same multigas component formulation. COMP3 models three components: GAS, SOLVENT and WATER. COMP4 models 4 components: OIL, GAS, SOLVENT and WATER. In both options the SOLVENT and GAS components share the gas phase, both may dissolve in the oil, and the solvent component may dissolve in the aqueous phase.

The degree of solution of the solvent component in the oil is defined by the table PVCOS, which follows the pattern of the PVCO table.

The PVT properties for the solvent may be specified explicitly using PVDS, or obtained from a property database (currently this is always CO2 properties generated using the Span-Wagner equation). In most respects the multi-gas model follows the existing solvent model – for example to inject the solvent component:

WELL_TYPE SOLVENT_INJECTOR
BHPL     400    Bar
TARG_SSV 2831684 m^3/day


ISOTHERMAL¶

Goes in the OPTIONS block of SIMULATION. If present, the model will run as isothermal, maintaining the same temperature as given in the initialization. All flow models run with the thermal option as default, however COMP3, COMP4 and COMP currently support only the isothermal mode.

NOSLVSOL¶

Goes in the OPTIONS block of SIMULATION. It is valid only for COMP3 and COMP4 modes: if present turns off the the solution of solvent in the aqueous phase.

NOGASSOL¶

Goes in the OPTIONS block of SIMULATION. See the following example:

SIMULATION
SIMULATION_TYPE SUBSURFACE
PROCESS_MODELS
SUBSURFACE_FLOW Flow
MODE GAS_WATER
OPTIONS
RESERVOIR_DEFAULTS
ISOTHERMAL
NOGASSOL
/                        ! end options
/                          ! end subsurface_flow
/                            ! end process_models
/                              ! end simulation


This is valid only for the GAS_WATER module: if present turns off the the solution of gas in the aqueous phase.

RESERVOIR_DEFAULTS¶

Goes in the OPTIONS block of SIMULATION.

If present causes PFLOTRAN to automatically set some values that are commonly useful for simulating reservoir models. It is possible to use RESERVOIR_DEFAULTS and still override any of the selected defaults by using the appropriate keyword in the corresponding card. A message will be output warning that a RESERVOIR_DEFAULTS default has been overridden in this way.

The values set by RESERVOIR_DEFAULTS are listed in reservoir defaults values.

ANALYTICAL_JACOBIAN¶

Goes in the OPTIONS block of SIMULATION.

If present, causes the derivatives in the Jacobian matrix to be computed analytically instead of approximated. Leads to better convergence of the nonlinear solver in most cases.

Note that this option is automatically selected as part of RESERVOIR_DEFAULTS.

HYSTERESIS¶

Goes in the OPTIONS block of SIMULATION. See the following example:

SIMULATION
SIMULATION_TYPE SUBSURFACE
PROCESS_MODELS
SUBSURFACE_FLOW Flow
MODE GAS_WATER
OPTIONS
RESERVOIR_DEFAULTS
ISOTHERMAL
HYSTERESIS
/                      ! end options
/                        ! end subsurface_flow
/                          ! end process_models
END                          ! end simulation


Activate the relative permeability hysteresis modelling in the simulation using the Carlson model [C+81]. The user must supply the indices of the required CHARACTERISTIC_CURVES sections to be used for the imbibition process using the IMBNUM grid keyword.

HYSTERESIS is available for the GAS_WATER, COMP3, COMP4 and COMP modes. If data is supplied, GAS_WATER can model relative permeability hysteresis of both wetting and non-wetting phases. Instead COMP3, COMP4 and COMP can only model relative permeability hysteresis of the gas phase.

HYSTERESIS_PC, EHYST¶

Goes in the OPTIONS block of SIMULATION.

HYSTERESIS_PC specifies the curvature parameter for capillary pressure hysteresis, [Kil76]. It takes a number as argument, which must be greater than zero. See example below:

HYSTERESIS_PC 0.1


When HYSTERESIS_PC is entered both the capillary pressure and relative permeability hysteresis are active. The user must supply the indices of the required CHARACTERISTIC_CURVES sections to be used for the imbibition process using the IMBNUM grid keyword.

This option is only available for the GAS_WATER mode.

RESTART¶

Goes in the SIMULATION block of SIMULATION.

RESTART specifies the file that contains a previous state of the simulation in PFLOTRAN format, from which the simulation can be restarted. It takes a file name as argument. For more details see RESTART. An example is given below:

SIMULATION
[process models etc blocks]
RESTART restart_file.h5
/


If not specified, the simulation starts from its initial state.

ERESTART¶

Goes in the SIMULATION block of SIMULATION.

ERESTART specifies the file that contains a previous state of the simulation in Eclipse format restart file, from which the simulation can be restarted. Takes as an argument a file name. For more detail see ERESTART. An example is given below:

SIMULATION
[process models etc blocks]
ERESTART base 3000 d
/                                      ! end simulation


If not specified, the simulation starts from its initial state.

CHECKPOINT¶

Goes in the SIMULATION block of SIMULATION.

Define instructions to save the simulation states for a given set of times. For more details see CHECKPOINT. An example is given below:

SIMULATION
[process models etc blocks]
CHECKPOINT
FORMAT HDF5
PERIODIC TIMESTEP 20
/                          ! end checkpoint
/                            ! end simulation


TL_OMEGA¶

Goes in the OPTIONS block of SIMULATION.

Only required for the TODD_LONGSTAFF and SOLVENT_TL modes.

Takes one argument, which specifies the Todd-Longstaff model omega value, which may lie between 0 and 1. The default value is 0. An example is given below:

SIMULATION
SUBSURFACE_FLOW Flow
MODE  TODD_LONGSTAFF
OPTIONS
ISOTHERMAL
RESERVOIR_DEFAULTS
TL_OMEGA 0.5
/                      ! end options
/                        ! end subsurface_flow
/                          ! end simulation


Goes in the SIMULATION block of SIMULATION.

For advanced users or debugging purposes only.

Causes additional information about PFLOTRAN_OGS to be written to the .out file, under the heading PROVENANCE.

The time and date of compilation and the current PFLOTRAN_OGS Git commit can always be found in this section. When the keyword ABOUT is included, much more detailed information about the compilation and version control is included.

PREAMBLE¶

Goes in the SIMULATION block of SIMULATION.

Causes some additional information about the run setup to be output to the screen before the run starts. This information can also be found in the run’s .out file by default.