Step06 Obtain Facies Grid for Forward sim initial state

 

In the present workflow step the Spatial_grid containing the facies distribution for the initial forward simulation state (including formations 01 to 07) will be obtained. Such spatial grid should not be confused with the ones output from restoration, which contain the facies distribution for each layer at depositional state / location.

 

The facies grid will be obtained by reading the restoration results of the corresponding plot file (using Read_plt_file) and performing a spatial grid output.

 

The main data can be found in Geol_003\Case02\Step06_Obtain_Facies_Grid\Data and comprises:

 

1.Geol_003_Case02_Step06_gridout.dat with data defined to read the target plot file and output a spatial grid with the facies variable.

2.Geol_003_Case01_Step04_rest_004.plt which corresponds to the plot file output at the end of the restoration stage in which formation08 was translated to the restoration surface. Note that the file Geol_003_Case01_Step04_rest_004.xmf is not required but is provided so that the user can visualize the restoration results in ParaView.

 

 

 

 

Geol_003_Case02_Step06_gridout.dat

 

The data in Geol_003_Case02_Step06_gridout.dat file comprises:

 

1.Read_plt_file to read the target restoration result plot file.

2.Spatial_grid_output to output the facies spatial grid.

 

 

Below a detailed description of the data structures is provided.

 

 

Read_plt_file

 

Data File


 

* Read_plt_file

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

 File_name    "Geol_003_Case01_Step04_rest_004.plt"

 

The name of the plot file being read is specified.

It should be noted that the relative location of any data with the Read_plt_file should be the same as it would be with Control_data (i.e. as Spatial_grid_output would be defined after a Control_data, it must also be defined after the Read_plt_file data).

 

 

 

Spatial_grid_output

 

Data File


 

* Spatial_grid_output

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

 Name                       "Fac"

 Stage_name_flag                0

 Grid_name_flag                 1

 File_format              "ascii"  

 Mesh_topo_format               0  

 Element_variables  IDM=1

  "Facies"

 Groups  IDM=7

  "formation01"

  "formation02"

      ...

  "formation07"

 

The Name of the grid to be output is defined as "Fac"

Stage_name_flag set to 0 indicates that the simulation stage name should not be appended to the output .spat file.

Grid_name_flag set to 1 indicates that the name of the grid (in this case "Fac") must be appended to the output .spat file.

The File_format is set to ascii (standard ParaGeo data structures and keywords) so that a .spat file is output.

Mesh_topo_format indicates the type of grid being output. In this case is set to 0 which corresponds to a Spatial_grid with the current mesh topology with data output on a group by group basis (each formation would have its corresponding Spatial_grid_group data).

Element_variables is used to specify that Facies element variable must be output.

Groups keyword is used to define the name of the groups for which the grid is going to be output. In the present case it is compulsory to define it as otherwise the grid is output for all the groups by default and the restoration plot contains the top surface membrane group (which would lead to a fatal error upon attempting the output of the Facies variable from it).

 

 

 

 

 

Results

 

The results for the present workflow steps are provided in Geol_003\Case2\Step06_Obtain_Facies_Grid\Results. The main result consists on the Geol_003_Case02_Step06_gridout_Fac.spat file which contains the spatial grid with the facies distribution on a group by group basis, i.e.:

 

 

Data File


 

* Spatial_grid                    NUM=1

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

 Name                       "Fac"

 Type                     "Group"

 Groups              IDM=7

    "formation01"

    "formation02"

        ...

    "formation07"

 Element_variables   IDM=1

    "Facies"

 

 

* Spatial_grid_group               NUM=2

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

 Name                     "formation01"

 Group_assignment         "formation01"

 Grid_name                "Fac"

 Element_type             "TET4"

 Configuration_type       "Current"

 Element_variables                IDM=1

    "Facies"

 Topology  IDM=4   JDM=13101

    11146     1046    11160    30778

    11185    11244    11178    29914

             ...

    30512    11875    30810    30869

 Node_numbers  IDM=3288

    29    30   ...     30898

 Coordinates  IDM=3  JDM=3288

   50951.6        28000.0        1470.69    

   56000.0        28000.0        1753.00

              ...  

   53983.5        30760.3        1372.76    

 Element_values  IDM=1  JDM=13101

   4.00000    

   4.00000    

     ...

   4.00000

   

The output .spat file contains a Spatial_grid of Type "Group" with values for the element variable "Facies" for formations 01 to 07 being defined. Such type of grid requires a Spatial_grid_group data for each of the formations listed in Groups keyword.

 

The .spat file contains also the Spatial_grid_group for each of the formations comprising the following data:

Name of the spatial grid group (which coincides with the Group_assignment)

The "parent" Spatial_grid name (in this case "Fac")

The Element_type defining the topology of the grid

Configuration_type set to "Current" indicating that the grid contains the data in the configuration at the output time

The Element_variables with values defined in the present grid

The Topology of the grid with the node numbers for each element/cell in a separate row (which corresponds to the mesh of the restoration at the output time)

The Node_numbers defining the topology of the grid

The Coordinates for each of the nodes defining the grid topology

The Element_values defining the values for the Facies property at each grid element/cell.