Mesh For Illinois Basin




Tet Mesh Stacked Layers

All files in a single tar/gz file:
all_files.tgz

Individual Files:
ill_01_lower.gmv
ill_01_lower.inp
ill_01_lower.table
ill_02_middle.gmv
ill_02_middle.inp
ill_02_middle.table
ill_03_upper.gmv
ill_03_upper.inp
ill_03_upper.table

The files *.gmv are GMV format files.
The files *.inp are AVS format finite element files.
The files *.table are the files we talked about on the phone that we created by just outputing the element attribute information from the AVS files.

*lower* are materials 1-5
*middle* is material 6
*upper is material 7

Material # vs Geologic Layer Names
7 Pennsylvanian-Mississippian Sandstones
6 New Albany Shale
5 Devonian Carbonates
4 Maqoketa Shale
3 Ordivician Carbonate
2 Eau Claire Shale
1 Mt Simon Sandstone
tesselate.lanl.gov:big_domain% more ill_01_lower.table
0 0 0 10 0
00010 1 1 1 1 1 1 1 1 1 1
evol, real
x_coord_elem, real
y_coord_elem, real
z_coord_elem, real
ie_index, integer
je_index, integer
dx, real
dy, real
av_thick, real
id_mat, integer
1 0.728560486400E+10 0.925000000000E+05 0.446250000000E+07 -0.347510772705E+03 1 108 0.500000000000E+03 0.500000000000E+03 0.291424194560E+05 1
2 0.484336435200E+10 0.925000000000E+05 0.446250000000E+07 -0.104931396484E+03 1 108 0.500000000000E+03 0.500000000000E+03 0.193734574080E+05 2

Where the columns are:

element number
evol = hex element volume m**3
x_coord_elem = x coord, centroid
y_coord_elem = y coord, centroid
z_coord_elem = z coord, centroid
ie_index = i,j indexing of column of elements
je_index = i,j indexing of column of elements
dx = hex cell dx in meters
dy = hex cell dy in meters
av_thick = cell average thickness computed from evol/(dx*dy)
id_mat = material ID number


Image Gallery Tet Mesh:

Click on image to enlarge.

Gallery written: Thu Jun 3 15:48:59 2010

Image Directory: /scratch/ymp/gable/grid_gen/www_mesh_catalog/nl_person/mt_simon01

12_surfs_01
LaGriT Finite Element Mesh
12_surfs_02
LaGriT Finite Element Mesh
12_surfs_03
LaGriT Finite Element Mesh
material_layers_02_100X_vert
LaGriT Finite Element Mesh
material_layers_100X_vert
LaGriT Finite Element Mesh
material_layers_100X_vert_xcut_01
LaGriT Finite Element Mesh
material_layers_100X_vert_xcut_02
LaGriT Finite Element Mesh
material_layers_100X_vert_xcut_03
LaGriT Finite Element Mesh
material_layers_100X_vert_ycut_01
LaGriT Finite Element Mesh
material_layers_100X_vert_ycut_02
LaGriT Finite Element Mesh
material_layers_100X_vert_ycut_03
LaGriT Finite Element Mesh
material_layers_grid_100X_vert
LaGriT Finite Element Mesh
tri1
LaGriT Finite Element Mesh
triangle_mesh
LaGriT Finite Element Mesh
tri_mesh_close
LaGriT Finite Element Mesh
voronoi_mesh_close
LaGriT Finite Element Mesh
voronoi_mesh
LaGriT Finite Element Mesh
voronoi_tri_mesh_close
LaGriT Finite Element Mesh




Hex Mesh Stacked Layers

stack_quads_sort.gmv
stack_quads_sort.inp
stack_quads_sort.table

LaGriT Control file and surfaces used to build 3D mesh

stack_clippedxy.lgi
xyclip_Mt_Simon_bottom.inp
xyclip_Top_ill01.inp
xyclip_Top_ill02.inp
xyclip_Top_ill03.inp
xyclip_Top_ill04.inp
xyclip_Top_ill05.inp
xyclip_Top_ill06a.inp
xyclip_Top_ill06b.inp
xyclip_Top_ill06.inp
xyclip_Top_ill07.inp
xyclip_Top_ill08.inp
xyclip_Top_ill09.inp
xyclip_Top_ill10.inp
xyclip_Top_ill11.inp

Gallery written: Fri Nov 19 12:24:30 2010

Image Directory: /scratch/ymp/gable/grid_gen/mt_simon/data_091211/big_domain/fehm_demo

fehm_mesh_01
LaGriT Finite Element Mesh
fehm_mesh_02
LaGriT Finite Element Mesh
fehm_mesh_03
LaGriT Finite Element Mesh
fehm_mesh_04
LaGriT Finite Element Mesh
fehm_mesh_composit
LaGriT Finite Element Mesh


Top of Stratigraphic units (in meters)

.inp file name

Unit Name
aquifer or confining unit

Permian = Land Surface

Top_ill01.inp





Pennsylvanian-Mississippian Sandstones??? aquifer

Pennsylvanian

Top _ill02.inp





??? ???

Upper Mississippian

Top _ill03.inp





??? ???

Lower Mississippian

Top _ill04.inp





??? ???

Upper Devonian

Top _ill05.inp





???
???
???
Top_ill06b.inp




???
???
???
Top_ill06a.inp




Devonian Carbonates??? ???

Middle Devonian

Top _ill06.inp





Devonian Carbonates??? ???

Upper SIlurian

Top _ill07.inp





??? ???

Lower Silurian

Top _ill08.inp





??? ???

Upper Ordovician

Top _ill09.inp





Ordivician Carbonate??? ???

Lower Ordovician

Top _ill10.inp





Ordivician Carbonate??? ???

Upper Cambrian (Mt. Simon top)

Top _ill11.inp





Mt Simon Sandstone aquifer

Mt. Simon bottom topography

Mt_Simon_bottom.inp





Quad and Tri surface meshes with isopach attributes

/scratch/ymp/gable/grid_gen/mt_simon/data_091211/big_domain

Compute elevation of base of aquifer and isopach of aquifer and confining units.

Starting with the following 8 surfaces as input:
/scratch/ymp/gable/grid_gen/mt_simon/data_091211/big_domain
xyclip_Top_ill01.inp
aquifer Pennsylvanian-Mississippian Sandstones
xyclip_Top_ill06b.inp
confining New Albany Shale
xyclip_Top_ill06a.inp
aquifer Devonian Carbonates
xyclip_Top_ill06.inp
confining Maqoketa Shale
xyclip_Top_ill07.inp
aquifer Ordivician Carbonate
xyclip_Top_ill10.inp
confining Eau Claire Shale
xyclip_Top_ill11.inp
aquifer Mt Simon Sandstone
xyclip_Mt_Simon_bottom.inp
Step 1 : Use top layer (xyclip_Top_ill01.inp) to truncate all the other layers. That is, if the elevation of a layer is greater than the elevation of the top layer, set the elevation to the elevation of the truncating layer.
Step 2 : Compute isopach from pairs of surface elevations.
Step 2a : Special case of the Eau Claire. If the Eau Claire isopach > 300, set it to 300. The reduction in thickness of the Eau Claire confining unit is added to the thickness of the overlaying Ordivician carbonate aquifer.
Step 3 : If all 4 nodes of a quad z_cu_isopach = 0.0
             AND
             If all 4 nodes of a quad z_aq_isopach = 0.0
             THEN Remove the element
Step 4 :
OUTPUT:
isopachs1_msss.inp
isopachs1_msss.gmv
isopachs2_ordcarb.inp
isopachs2_ordcarb.gmv
isopachs3_dsss.inp
isopachs3_dsss.gmv
isopachs4_pmss.inp
isopachs4_pmss.gmv

Each file has 3 node attributes:
z_aq_bot (elevation of base of aquifer, m)
z_aq_isopach (thickness of aquifer, m)
z_cu_isopach (thickness of overlaying confining unit, m)

Some statistics:
Mount Simon Sandstone (msss)
 z_aq_isopach       3.725516012E+01  6.193169592E+02 5.820617991E+02
 z_cu_isopach       6.685144404E+01  3.000000000E+02 2.331485560E+02

Ord. Carbonates (orcarb)
 z_aq_isopach       7.876838593E+00  1.052749671E+03 1.044872833E+03
 z_cu_isopach       2.690042130E-02  2.379591396E+02 2.379322392E+02

Dev-Suliurian SS (dsss)
 z_aq_isopach       6.357165042E-01  2.500000000E+02 2.493642835E+02
 z_cu_isopach       0.000000000E+00  7.500000000E+01 7.500000000E+01

Penn-Miss Ss (pmss)
 z_aq_isopach       0.000000000E+00  1.359867474E+03 1.359867474E+03
 z_cu_isopach       0.000000000E+00  0.000000000E+00 0.000000000E+00

Similarly, create surfaces and attributes for triangulated surfaces:
/scratch/ymp/gable/grid_gen/mt_simon/data_091211/big_domain

lagrit < interpolate_onto_tri.lgi
   interp.mlgi

OUTPUT:
isopachs1_tri_msss.inp
isopachs1_tri_msss.gmv
isopachs1_tri_ordcarb.inp
isopachs1_tri_ordcarb.gmv
isopachs1_tri_dsss.inp
isopachs1_tri_dsss.gmv
isopachs1_tri_pmss.inp
isopachs1_tri_pmss.gmv

DOWNLOADS:
The following file:
all_isopachs.tgz

Contains these files:
isopachs1_msss.gmv
isopachs1_msss.inp
isopachs1_tri_msss.gmv
isopachs1_tri_msss.inp
isopachs2_ordcarb.gmv
isopachs2_ordcarb.inp
isopachs2_tri_ordcarb.gmv
isopachs2_tri_ordcarb.inp
isopachs3_dsss.gmv
isopachs3_dsss.inp
isopachs3_tri_dsss.gmv
isopachs3_tri_dsss.inp
isopachs4_pmss.gmv
isopachs4_pmss.inp
isopachs4_tri_pmss.gmv
isopachs4_tri_pmss.inp

LaGriT Control files:
Quad surfaces:
lagrit < compute_isopachs_modify_ec_v02.lgi
   truncate_layers.mlgi
   add_att_copy_att.mlgi
   delete_elem_mod_att.mlgi

Tri surfaces:
lagrit < interpolate_onto_tri.lgi
   interp.mlgi


Image Gallery:

Click on image to enlarge.

Gallery written: Tue Nov 30 13:54:48 2010

Image Directory: /scratch/ymp/gable/grid_gen/mt_simon/data_091211/big_domain/images_isopachs

dsss_aq_bottom
LaGriT Finite Element Mesh
dsss_aq_isopach
LaGriT Finite Element Mesh
dsss_cu_isopach
LaGriT Finite Element Mesh
msss_aq_bottom
LaGriT Finite Element Mesh
msss_aq_isopach
LaGriT Finite Element Mesh
msss_cu_isopach
LaGriT Finite Element Mesh
ordcarb_aq_bottom
LaGriT Finite Element Mesh
ordcarb_aq_isopach
LaGriT Finite Element Mesh
ordcarb_cu_isopach
LaGriT Finite Element Mesh
pmss_aq_bottom
LaGriT Finite Element Mesh
pmss_aq_isopach
LaGriT Finite Element Mesh
pmss_cu_isopach
LaGriT Finite Element Mesh
quad_tri_mesh_01
LaGriT Finite Element Mesh
quad_tri_mesh_02
LaGriT Finite Element Mesh
quad_tri_mesh_03
LaGriT Finite Element Mesh
stacked_modified_eau_claire
LaGriT Finite Element Mesh



Cross Sections Extracted From 3D Mesh

INPUT:
stack_fill.gmv

LaGriT Control File:
extract_x_sections.lgi

OUTPUT:
merge_y_xs_tri.inp
merge_x_xs_tri.inp

Cross Section X

Cross Section Y