MATERIAL_PROPERTY with permeability and porosity¶
Note that if a mesh format different than GRDECL is used, the porosity and permeability values will be taken from MATERIAL_PROPERTY. Therefore, in addition to the data input described in MATERIAL_PROPERTY, the user must enter the PERMEABILITY and POROSITY block:
MATERIAL_PROPERTY <material_name>
ID <material_id>
POROSITY [<constant_porosity>, DATASET <porosity_dataset>]
<PERMEABILITY BLOCK>
CHARACTERISTIC_CURVES <ch_curve_name>
ROCK_DENSITY <rock_density>
SOIL_COMPRESSIBILITY <soil_compress_coeff>
SOIL_COMPRESSIBILITY_FUNCTION <compress_function_type>
SOIL_REFERENCE_PRESSURE [<ref_press> or INITIAL_PRESSURE]
SPECIFIC_HEAT <spec_heat> or HEAT_CAPACITY <heat_cap>
THERMAL_CONDUCTIVITY_DRY <th_cond_dry>
THERMAL_CONDUCTIVITY_WET <th_cond_wet>
END
PERMEABILITY¶
Note that permeabilities are specified in m2 , rather than milliDarcies, as used by Eclipse.
<PERMEABILITY BLOCK> depends on the type of permeability field to be used, as follows:
<PERMEABILITY BLOCK> for Isotropic and homogeneous permeability fields:
PERMEABILITY
ISOTROPIC: Turn off anisotropy.
PERM_ISO <iso_perm_value>: specifies a constant isotropic permeability value [m2]
END
<PERMEABILITY BLOCK> for anisotropy and homogeneous permeability fields:
PERMEABILITY
PERM_X <perm_x_value>: specifies a constant permeability value in x-direction [m2]
PERM_Y <perm_y_value>: specifies a constant permeability value in y-direction [m2]
PERM_Z <perm_z_value>: specifies a constant permeability value in z-direction [m2]
END
<PERMEABILITY BLOCK> for anisotropic and heterogeneous permeability fields:
In this case the permeabilities are taken from a ‘dataset’, a block of data which is normally loaded from a file:
PERMEABILITY
ANISOTROPIC or ISOTROPIC: specifies the type of permeability field
DATASET <dataset_name>: dataset where the permeability field is defined [m2]
END
In case of ANISOTROPIC permeability, three data sets must be present in the input deck, named <dataset_name>X, <dataset_name>Y and <dataset_name>Z respectively.
POROSITY¶
This can be defined in two alternative ways:
<constant_porosity>: assigning a constant value for all cells belonging to the material being defined.
DATASET <porosity_dataset>: assigning porosity values cell-by-cell, where the values are provided in a DATASET defined in the input deck.
The example below illustrates a heterogeneous anisotropic case. Note that in addition to this information, the required datasets must have been loaded - see DATASET for details of how to do this:
MATERIAL_PROPERTY formation
ID 1
POROSITY DATASET Porosity
ROCK_DENSITY 2.350d3
SPECIFIC_HEAT 1.0d3
THERMAL_CONDUCTIVITY_DRY 1.541d0
THERMAL_CONDUCTIVITY_WET 1.541d0
CHARACTERISTIC_CURVES ch1
PERMEABILITY
ANISOTROPIC
DATASET perm
/
/
DATASET¶
Specifies a data set to be associated with parameter sets in the model.
DATASET <dataset_name>
FILENAME <dataset_file_name>: specifies the external file name of the dataset
HDF5_DATASET_NAME <hdf5_dataset_name>: specifies the dataset name
END
See example below for a dataset defining the permeability in the x direction
DATASET permX
FILENAME perm_poro.h5
HDF5_DATASET_NAME PermeabilityX
END
Note such datasets must follow the conventions of an HDF5 file. Such datasets normally need to be prepared by external programs which link in the HDF5 libraries.