Brick Version octree mesh
/scratch/fwo/tam/wy_2011/grid_brick
EV files
/scratch/fwo/wyo_ev_models/LaBarge1
FEHM Result files
/scratch/fwo/stauffer/Wyoming/Moxa
1 "Gallatin" Gallatin
2 "Bighorn" Bighorn
3 "Madison" Madison
4 "Amsden" Amsden
5 "Weber" Weber
6 "Phosphoria" Phosphoria
7 "Chugwater" Chugwater
8 "Ankareh" Ankareh
9 "Nugget" Nugget
10 "Morrison" Morrison
11 "Dakota" Dakota
12 "Frontier" Frontier
13 "Baxter" Baxter
14 "" Outside
First pass brick mesh with horizontal refinement at well location.
The point count is 97156.
The element count is 540932.
Background mesh with 256x256x32 meter spacing
Refined 3 levels horizontaly along x,y location of well
So resolution at well location is 32x32x32
MIN MAX DIFFERENCE NODES
METERS
xic 1.448530000E+05 1.497170000E+05 4.864000000E+03 97156
yic 1.426160000E+05 1.518320000E+05 9.216000000E+03 97156
zic -3.840000000E+03 -4.480000000E+02 3.392000000E+03 97156
FEET
xic 4.752395013E+05 4.911975066E+05 1.595800525E+04 97156
yic 4.679002625E+05 4.981364829E+05 3.023622047E+04 97156
zic -1.259842520E+04 -1.469816273E+03 1.112860892E+04 97156
Material Nodes
Baxter 13 has 21425 nodes. #nodes/nnodes is 0.220521628857
Frontier 12 has 6347 nodes. #nodes/nnodes is 0.653279274702E-01
Dakota 11 has 3804 nodes. #nodes/nnodes is 0.391535274684E-01
Morrison 10 has 9940 nodes. #nodes/nnodes is 0.102309688926
Nugget 9 has 12344 nodes. #nodes/nnodes is 0.127053394914
Ankareh 8 has 5329 nodes. #nodes/nnodes is 0.548499338329E-01
Chugwater 7 has 6997 nodes. #nodes/nnodes is 0.720181986690E-01
Phosphoria 6 has 4695 nodes. #nodes/nnodes is 0.483243428171E-01
Weber 5 has 3180 nodes. #nodes/nnodes is 0.327308662236E-01
Amsden 4 has 3051 nodes. #nodes/nnodes is 0.314031057060E-01
Madison 3 has 10764 nodes. #nodes/nnodes is 0.110790893435
Bighorn 2 has 4026 nodes. #nodes/nnodes is 0.414385125041E-01
Gallatin 1 has 5254 nodes. #nodes/nnodes is 0.540779791772E-01
THE PSET pxy at well location HAS 107 POINTS
THE PSET pinj at injection interval HAS 12 POINTS
tet_node_attributes.dat (where well xy is 50 and well interval is 51)
THE PSET pxy at well location HAS 107 POINTS
THE PSET pinj at injection interval HAS 12 POINTS
# node, integer
# imt1, integer
# vorvol, real
# iwell, integer
# xcoord, real
# ycoord, real
# zcoord, real
awk '{if ($4 > 50) print }' tet_node_attributes.dat
...
4650 2 0.327679999999E+05 50 0.147157000000E+06 0.145176000000E+06 -0.364800000000E+04
5390 2 0.327679999997E+05 50 0.147157000000E+06 0.145176000000E+06 -0.361600000000E+04
6130 2 0.327680000000E+05 51 0.147157000000E+06 0.145176000000E+06 -0.358400000000E+04
6870 3 0.327680000000E+05 51 0.147157000000E+06 0.145176000000E+06 -0.355200000000E+04
7610 3 0.327680000006E+05 51 0.147157000000E+06 0.145176000000E+06 -0.352000000000E+04
8350 3 0.327679999987E+05 51 0.147157000000E+06 0.145176000000E+06 -0.348800000000E+04
9090 3 0.327680000012E+05 51 0.147157000000E+06 0.145176000000E+06 -0.345600000000E+04
9830 3 0.327680000013E+05 51 0.147157000000E+06 0.145176000000E+06 -0.342400000000E+04
10570 3 0.327680000005E+05 51 0.147157000000E+06 0.145176000000E+06 -0.339200000000E+04
11310 3 0.327679999987E+05 51 0.147157000000E+06 0.145176000000E+06 -0.336000000000E+04
12050 3 0.327680000022E+05 51 0.147157000000E+06 0.145176000000E+06 -0.332800000000E+04
12790 3 0.327680000000E+05 51 0.147157000000E+06 0.145176000000E+06 -0.329600000000E+04
13530 3 0.327679999990E+05 51 0.147157000000E+06 0.145176000000E+06 -0.326400000000E+04
14270 3 0.327680000000E+05 51 0.147157000000E+06 0.145176000000E+06 -0.323200000000E+04
15010 3 0.327679999953E+05 50 0.147157000000E+06 0.145176000000E+06 -0.320000000000E+04
15750 4 0.327679999988E+05 50 0.147157000000E+06 0.145176000000E+06 -0.316800000000E+04
...
FEHM files:
-rw-r--r-- 1 tam pf 29076342 Feb 22 11:44 tet.fehmn
-rw-r--r-- 1 tam pf 18838609 Feb 22 11:45 tet.stor
-rw-r--r-- 1 tam pf 1069157 Feb 22 11:44 tet_material.zone
-rw-r--r-- 1 tam pf 1697350 Feb 22 11:44 tet_multi_mat.zone
-rw-r--r-- 1 tam pf 154411 Feb 22 11:44 tet_outside.zone
-rw-r--r-- 1 tam pf 890543 Feb 22 11:44 tet_outside_vor.area
Face top 1 has 908 nodes.
top Sum Voronoi Area_x Area_y Area_z
0.2949120E+06 0.1556480E+06 0.4482662E+08
Face bottom 2 has 908 nodes.
bottom Sum Voronoi Area_x Area_y Area_z
0.2949120E+06 0.1556480E+06 0.4482662E+08
Face left_w 3 has 3959 nodes.
left_w Sum Voronoi Area_x Area_y Area_z
0.3126067E+08 0.8683520E+06 0.2359296E+07
Face right_e 5 has 3959 nodes.
right_e Sum Voronoi Area_x Area_y Area_z
0.3126067E+08 0.8683520E+06 0.2359296E+07
Face back_n 6 has 2140 nodes.
back_n Sum Voronoi Area_x Area_y Area_z
0.8683520E+06 0.1649869E+08 0.1245184E+07
Face front_s 4 has 2140 nodes.
front_s Sum Voronoi Area_x Area_y Area_z
0.8683520E+06 0.1649869E+08 0.1245184E+07
*** Construct and Compress Sparse Matrix:3D ***
*** Compress Area Coefficient Values ***
AMatbld3d_stor: *****Zero Negative Coefficients ******
AMatbld3d_stor: Number of 'zero' (< 10e-8*max) coefs 0
AMatbld3d_stor: npoints = 97156 n connections = 681486
AMatbld3d_stor: Number of written coefs = 28
AMatbld3d_stor: Maximum num. connections to a node = 9
AMatbld3d_stor: Volume min = 1.6384000E+04
AMatbld3d_stor: Volume max = 2.0971520E+06
AMatbld3d_stor: Total Volume: 1.5205191E+11
AMatbld3d_stor: abs(Aij/xij) min = 0.0000000E+00
AMatbld3d_stor: abs(Aij/xij) max = 2.0480000E+03
AMatbld3d_stor: (Aij/xij) max = 0.0000000E+00
AMatbld3d_stor: (Aij/xij) min = -2.0480000E+03
Compression used for graph and coefficient values
The area coefficient values were written as scalar area/distance
ascii STOR file written with name tet.stor
* Project: Wyoming Seismic with John Jiao 2011
* Key: SequenceFile Data: seis3d.seq
* range FEET:
* 475200.000000 to 493200.000000
* 465500.000000 to 499900.000000
* -12600.000000 to 0.000000
* range METERS:
* 144840.96 to 150327.36
* 141884.4 to 152369.52
* -3840.48 to 0
* NOTE: There is a cut-out at min x,y corner that is aprox 450 meters
* Increase min y to avoid this area
*
* State Plane feet, zone: 4904
* The well location is at X 482800, Y 476300 ft (147157.44, 145176.24 meters)
*
*The mesh domain should be:
* X 475200 – 493200 ft
* Y 467100 – 499900 ft
* Z -1400 – -12600 ft elevation
* (-426.72 meters to -3840.48 meters)
* The target interval Weber Sandstone at the well for CO2 injection
* ranges from -10600 to -11800 ft (elevation) or depth 11300 to 18500 ft in depth.
* -3230.88 to -3596.64
*
* Target 50 m horizontal and 32 meter vertical near well
* Data collection points were every 32 meters vertical
* octree ratio length/sqrt(2) = height
* 50/sqrt(2) = 35.35
* 45/sqrt(2) = 31.81
* 64/sqrt(2) = 45.25
*
* select boundary coordinates within EV domain
* use integer values that divide whole
* use meters
* X length 17500 feet, 5334 meters
* 20*256=5120
* 19*256=4864
* increase x to cell center the well
* need to reduce by one cell length
define X0 144853.
define X1 149717.
* Y length 32000 feet, 9753.6 meters
* 38*256=9728
* increase y min 512m, 2 cells so 36*256=9216
* increase y to center well in cell
define Y0 142616.
define Y1 151832.
* Z length 11200 feet, 3413.76 meters
* -11200 to -1400 feet, 3840.48 to 426.72 meters
* 53*64=3392 106*32=3392
define Z0 -3840.
define Z1 -448.
define / NX / 20
define / NY / 37
define / NZ / 107
Begin forwarded message:
From: Zun Sheng Jiao <JJiao@uwyo.edu>
Date: February 7, 2011 10:06:58 AM MST
To: Phil Stauffer <pheeel.gud@gmail.com>, "Terry A. Miller" <tamiller@lanl.gov>
Subject: EV Model
Terry,
Attached are the figures from the Earthvision models.
The mesh domain should be:
X 475200 – 493200 ft
Y 467100 – 499900 ft
Z -1400 – -12600 ft elevation, (depth 8100 – 19300 ft)
The well location is at X 482800 ft and Y 476300 ft. The target interval Weber Sandstone at the well location for CO2 injection ranges from -10600 to -11800 ft (elevation) or depth 11300 to 18500 ft in depth.
For starting we may use 100 ft gridding size at the near injection location and then gradually coarsening away from the well. Let Phil to decide what interval we should be used.
Terry, could you please write down the procedures and scripts you used to create this mesh? I deeply appreciated your hard efforts and patience for me.
Best regards,
John Jiao
Wyoming State Geological Survey
PO Box 1347
Laramie WY 82073
(307)766-2286 Ext. 225
(307)766-2605 (fax)
Begin forwarded message:
From: Zun Sheng Jiao <JJiao@uwyo.edu>
Date: January 31, 2011 9:44:12 AM MST
To: Phil Stauffer <stauffer@lanl.gov>, "Terry A. Miller" <tamiller@lanl.gov>
Cc: Rodney Harold Debruin <RDebruin@uwyo.edu>, "Ronald C. Surdam" <RSurdam@uwyo.edu>
Subject: EV models
Phil, Terry,
I just upload the La Barge structural, porosity, and permeability models to LANL transfer site. The zip file Labarge_Seismic_PoroPerm_model.zip is included all files in EV folder. Please let me know they can be opened by your EV.
Thanks,
John Jiao
Wyoming State Geological Survey
PO Box 1347
Laramie WY 82073
(307)766-2286 Ext. 225
(307)766-2605 (fax)
/scratch/fwo/wyo_ev_models/LaBarge_Sesmic_PoroPerm_model
capillary.lanl.gov% cat README.txt
All horizons from Yuri *.txt file have datum 6700 ft
6700 ft - horizon picked = horizon in elevation.
int_vel.txt file: z is distance below to sea level (positive number),
Z*-1 to match with horizon in elevation
Good morning, Terry,
Three faces files you need to use, seis3d.faces for structure, por.por.faces for porosity, and perm_Vernik_perm.faces for permeability.
Thanks,
John Jiao
capillary.lanl.gov% ev_facedump seis3d.unsliced.faces
Version 6
Closed, Unsliced
Info
Key: SequenceFile Data: seis3d.seq
Key: WFPFile Data: H:Co2sequestrationUW-WSGS_CO23DSeismicPoroPerm_modelseis3d.wfp
Key: Title1 Data: Unknown
Key: Title2 Data: Unknown
Key: FacesFromType Data: ply
Key: Units Data: Unknown
Key: Accuracy Data: 0.250000
Space
range: 475200.000000 to 493200.000000, 465500.000000 to 499900.000000, -12600.000000 to 0.000000
x and y units: US feet
z units: feet
proj: State Plane, zone: 4904
ellipsoid: Clarke 1866/NAD27
projection parameters:
1) 6378206.400000 2) 6356583.800000 3) 0.000000
4) 0.000000 5) 0.000000 6) 0.000000
7) 0.000000 8) 0.000000 9) 0.000000
10) 0.000000 11) 0.000000 12) 0.000000
13) 0.000000 14) 0.000000 15) 0.000000
Total Space Volume: 7801920000000.000000
Fault Block Names:
[1] Block1
[2] Block2
[3] Block3
[4] Block4
[5] Block5
[6] Block6
Zone Names:
[1] Gallatin
[2] Bighorn
[3] Madison
[4] Amsden
[5] Weber
[6] Phosphoria
[7] Chugwater
[8] Ankareh
[9] Nugget
[10] Morrison
[11] Dakota
[12] Frontier
[13] Baxter
Volume[000] World /World/
Volume[001] Outside /World/Outside
Volume[002] Roi /World/Roi/
Volume[003] Zone /World/Roi/Block6/Baxter
Volume[004] FaultB /World/Roi/Block6/
Volume[005] Zone /World/Roi/Block6/Frontier
Volume[006] Zone /World/Roi/Block6/Gallatin
Volume[007] Zone /World/Roi/Block6/Bighorn
Volume[008] Zone /World/Roi/Block6/Madison
Volume[009] Zone /World/Roi/Block6/Amsden
Volume[010] Zone /World/Roi/Block6/Weber
Volume[011] Zone /World/Roi/Block6/Phosphoria
Volume[012] Zone /World/Roi/Block6/Chugwater
Volume[013] Zone /World/Roi/Block6/Ankareh
Volume[014] Zone /World/Roi/Block6/Nugget
Volume[015] Zone /World/Roi/Block6/Morrison
Volume[016] Zone /World/Roi/Block6/Dakota
Volume[017] Zone /World/Roi/Block1/Gallatin
Volume[018] Zone /World/Roi/Block2/Gallatin
Volume[019] Zone /World/Roi/Block5/Gallatin
Volume[020] Zone /World/Roi/Block3/Gallatin
Volume[021] Zone /World/Roi/Block4/Gallatin
Volume[022] Zone /World/Roi/Block4/Frontier
Volume[023] Zone /World/Roi/Block4/Dakota
Volume[024] Zone /World/Roi/Block5/Dakota
Volume[025] Zone /World/Roi/Block5/Frontier
Volume[026] Zone /World/Roi/Block5/Morrison
Volume[027] Zone /World/Roi/Block4/Bighorn
Volume[028] Zone /World/Roi/Block4/Madison
Volume[029] Zone /World/Roi/Block4/Morrison
Volume[030] Zone /World/Roi/Block4/Nugget
Volume[031] Zone /World/Roi/Block4/Baxter
Volume[032] Zone /World/Roi/Block5/Nugget
Volume[033] Zone /World/Roi/Block4/Ankareh
Volume[034] Zone /World/Roi/Block4/Chugwater
Volume[035] Zone /World/Roi/Block5/Chugwater
Volume[036] Zone /World/Roi/Block4/Phosphoria
Volume[037] Zone /World/Roi/Block5/Phosphoria
Volume[038] Zone /World/Roi/Block5/Baxter
Volume[039] Zone /World/Roi/Block4/Weber
Volume[040] Zone /World/Roi/Block5/Weber
Volume[041] Zone /World/Roi/Block4/Amsden
Volume[042] Zone /World/Roi/Block5/Amsden
Volume[043] Zone /World/Roi/Block5/Madison
Volume[044] Zone /World/Roi/Block5/Bighorn
Volume[045] Zone /World/Roi/Block5/Ankareh
Volume[046] Zone /World/Roi/Block2/Baxter
Volume[047] Zone /World/Roi/Block2/Frontier
Volume[048] Zone /World/Roi/Block2/Nugget
Volume[049] Zone /World/Roi/Block2/Phosphoria
Volume[050] Zone /World/Roi/Block2/Chugwater
Volume[051] Zone /World/Roi/Block2/Ankareh
Volume[052] Zone /World/Roi/Block2/Weber
Volume[053] Zone /World/Roi/Block2/Dakota
Volume[054] Zone /World/Roi/Block2/Amsden
Volume[055] Zone /World/Roi/Block2/Madison
Volume[056] Zone /World/Roi/Block2/Bighorn
Volume[057] Zone /World/Roi/Block2/Morrison
Volume[058] Zone /World/Roi/Block1/Bighorn
Volume[059] FaultB /World/Roi/Block1/
Volume[060] Zone /World/Roi/Block1/Madison
Volume[061] Zone /World/Roi/Block1/Amsden
Volume[062] Zone /World/Roi/Block1/Weber
Volume[063] Zone /World/Roi/Block1/Phosphoria
Volume[064] Zone /World/Roi/Block1/Chugwater
Volume[065] Zone /World/Roi/Block1/Ankareh
Volume[066] Zone /World/Roi/Block1/Nugget
Volume[067] Zone /World/Roi/Block1/Morrison
Volume[068] Zone /World/Roi/Block1/Dakota
Volume[069] Zone /World/Roi/Block1/Frontier
Volume[070] Zone /World/Roi/Block1/Baxter
Volume[071] FaultB /World/Roi/Block4/
Volume[072] Zone /World/Roi/Block3/Bighorn
Volume[073] FaultB /World/Roi/Block3/
Volume[074] Zone /World/Roi/Block3/Madison
Volume[075] Zone /World/Roi/Block3/Amsden
Volume[076] Zone /World/Roi/Block3/Weber
Volume[077] Zone /World/Roi/Block3/Phosphoria
Volume[078] Zone /World/Roi/Block3/Chugwater
Volume[079] Zone /World/Roi/Block3/Ankareh
Volume[080] Zone /World/Roi/Block3/Nugget
Volume[081] Zone /World/Roi/Block3/Morrison
Volume[082] Zone /World/Roi/Block3/Dakota
Volume[083] Zone /World/Roi/Block3/Frontier
Volume[084] Zone /World/Roi/Block3/Baxter
Volume[085] FaultB /World/Roi/Block2/
Volume[086] FaultB /World/Roi/Block5/