Val_005d Kerogen thickness

Problem Description

 

The user is assumed to have undertaken the example Val_005c beforehand.

 

The influence of the specification of kerogen layer thickness is demonstrated here.  The present example is a variation of the oil and gas expulsion example for kerogen B. The data files are identical to Val_005c for kerogen B with the exception of the data required to include kerogen layer thickness.  Note that this includes definition of stratigraphy data (described in detail later). When thickness is defined, the change in porosity and pore pressure are both scaled by "thickness/formation thickness". Therefore, for this single element example with dimensions of 10m x 10m (i.e. formation thickness = 10m), if thickness of e.g. 5m is defined, then the change in both porosity and pore pressure is reduced by a factor of 0.5.  Conversely, as in Val_005c, when thickness is not defined, then the computed element volume with thickness = formation thickness would be used.

 

Two examples are simulated here:

1Thickness set to 10m to validate against the result with no thickness definition - results should be exactly the same for changes in porosity and pore pressure.

2Thickness set to 5m to compare against thickness=10m results - results should be reduced by a factor of 0.5.

 

 

Data File Description

 

The data files for the oil and gas expulsion process for kerogen B with kerogen thickness definition are in ParaGeo Examples\Validation\Val_005\Val_005d\Data. Only the key data structures required for the kerogen thickness specification are described here, these include:

Kerogen_kinetics - include definition of kerogen thickness.

Geometry_set - include geometry definition for base and top of example model.

Stratigraphy data is required for the ParaGeo code to identify the base and top of the model to get the thickness:

oStratigraphy_definition - defines the unit name and basal horizon.

oStratigraphy_horizon - defines the top horizon data.

oStratigraphy_surface_load - define data to de-active the default stratigraphy surface loading.

 

Kerogen_kinetics

 

Kerogen_kinetics defines the assignment and processing data associated with kerogen kinetics property name "kerogen1", e.g. kerogen organofacies type, total organic carbon content (TOC), transformation index (TI), hydrocarbon index (HI), gas concentration ratio (G), weight fraction carbon (W), Universal gas constant, kerogen density and kerogen layer thickness.

 

Data File

 

 

* Kerogen_kinetics               NUM=1

! ---------------------------------

 Name    "kerogen1"

 Kerogen_facies_type       "B"

 TOC                           8.0

 Transformation_index        0.018

 Hydrocarbon_index           0.592

 Gas_concentration            0.17

 Weight_fraction_carbon       0.75

 Universal_gas_constant    0.00831446261815324

 Kerogen_density            1200.0

 Thickness                    10.0

 

1Thickness defines the kerogen layer thickness as 10m.

 

 

 

 

Geometry_set data

 

Geometry_set data defines the geometry entities for the base and top of the model.

 

Data File

 

 

* Geometry_set               NUM=1  

! -------------------------------

 Name   "Base"

 Lines  IDM=1

   1

 

* Geometry_set               NUM=12  

! -------------------------------

 Name   "Top"

 Lines  IDM=1

   3

 

1Two geometry sets "Base" and "Top" are defined to identify the geometry lines for the base and top of the model, respectively.

 

 

 

Stratigraphy data

 

Stratigraphy data is required for the ParaGeo code to identify the base and top of the model to get the thickness:

oStratigraphy_definition - defines the unit name and basal horizon.  For simplicity of stratigraphy data definition, the stratigraphy unit name and stratigraphy horizon name are defined with identical name to the Group_data name of "Default-1".  When this same naming convention is used, Horizon_numbers and Group_numbers are automatically set and need not be defined.

oStratigraphy_horizon - defines the top horizon geometry for the stratigraphy unit.

oStratigraphy_surface_load - define data to de-active the default stratigraphy surface loading.

 

 

Data File

 

 

* Stratigraphy_definition

! ---------------------------------

 Units  IDM=1

   "Default-1"

 Basal_horizon     Base

 

 

* Stratigraphy_horizon NUM=1

! ---------------------------------

 Name         "Default-1"

 Geometry_set      "Top"

 

 

* Stratigraphy_surface_load

! ---------------------------------

 Active_flag      -1

 

1Stratigraphy definition data:

a.Stratigraphy unit name defined as "Default-1" - same as the Group_data name.

b.Basal horizon defined as geometry set named "Base".

2Stratigraphy horizon data for stratigraphy unit "Default-1":

a.Stratigraphy horizon name also defined as "Default-1" for consistency.

b.Top surface of unit is defined by geometry set named "Top".

3Stratigraphy surface load defined with Active_flag = 1 to de-activate any default stratigraphy surface loading.

 

 

 

 

 

 

Results for porosity change and pore pressure generation with different kerogen thickness definition

 

The results files for the oil and gas expulsion for kerogen organofacies type B with kerogen thickness definitions are in ParaGeo Examples\Validation\Val_005\Val_005d\Results. In this directory, the spreadsheet mat_kerogenB_expulsion_thick.xlsx contains:

history results output from the simulations.

comparison plots of porosity change and pore pressure generation for Val_005c (no thickness definition) vs thickness=10m vs thickness=5m (Figure 1).

 

The figure below shows the ParaGeo results for (a) porosity change and (b) pore pressure generation with different kerogen thickness definition for kerogen B.  Observed are:

identical results for the case with no thickness definition (Val_005c) and the case with thickness=10m.  This is as expected because the kerogen thickness definition of 10m is the same as the formation thickness definition (i.e. 10m) which is used by default if no kerogen thickness is defined as in Val_005c.

results with thickness=5m are scaled down by a factor of 0.5 as expected due to kerogen thickness being halved.

 

Note that  the amount of pore pressure generated here is artificially high due to the pseudo-artificial conditions of the single element model.  A more suitable example to observe overpressure due to hydrocarbon generation would be a uniaxial column with sediment compaction and increasing temperature with burial Kin_001.

 

Influence of Kerogen Thickness (ParaGeo Results)

(a) Porosity change

(b) Pore pressure generation

Val_005d_001a

Val_005d_001b

Figure 1: ParaGeo results for (a) porosity change and (b) pore pressure generation resulting from the oil and gas expulsion for organofacies type B with different kerogen thickness definition