Case 5 Porosity Initialization via Spatial Grids

 

In this exercise the initialization of porosity will be performed via a Spatial_grid located within the file porosity.spat. Note that if porosity is initialized via a grid and the constitutive model is elasto-plastic, the volumetric plastic strain corresponding to each porosity value should also be provided in the grid in order to initialize at the correct pre-consolidation pressure. The spatial grid file (porosity.spat) contains a spatial_grid with:

19 cells in X direction (1000 m each)

22 cells in Y direction (1000 m each)

38 cells in Z direction (500 m each)

4X, Y and Z coordinates for the grid points

5Porosity and volumetric plastic strain at grid points

 

The porosity distribution defined by the grid is as follows:

1In X=0m the porosity vs depth trend follows the NCT (same porosity distribution as in Case 1).

2In X=9000m (outside from the domain) the porosity at a given depth is 0.8 · Porosity in NCT. The lateral variation of porosity is defined to be lineal (see the provided exel file in Geost_002\Exercices\Case5\Data\000 - Porosity Grid Definition.xlsx).

 

To perform this exercise copy the datafile from Case4 as a template. Then:

1Delete the Prorosity_spatial keyword from the Geostatic_data. If desired Spatial_variation_definition and Spatial_variation_values can be deleted as well as will not be used.

2Include the Spatial_grid in file porosity.spat using the Include data structure.

3Assign the values for porosity and volumetric plastic strain included in the grid using the Spatial_state_set data structure.

 

 

Data File

 

 

* Geostatic_data            NUM=1

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

 Groups                   IDM=4

  "formation0"

 "formation1"

 "formation2"

 "formation3"

 Pore_pressure_distribution   "Hydrostatic"

 Time_curve                   100

 

 

* Include

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

 Filename "porosity.spat"

 

 

* Spatial_state_set         NUM=1

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

 Spatial_grid  "Porosity_Spatial"

 Element_variables    IDM=1

  "Porosity"

 State_variables      IDM=1

  "P_strnv"        

1Porosity_spatial keyword is removed from Geostatic_data.

2The spatial grid within the file named porosity.spat is included in the analysis.

3The values for porosity and volumetric plastic strain are assigned as initial values by using the Spatial_state_set data structure.

 

Geost_002_Porosity_Spatial

Porosity initialized via a Spatial Grid (note the porosity variation in X direction)

 

 

 

The data  and result files for the project are in directory: geost_002\Exercises\Case5