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 Date: February 7, 2011 10:06:58 AM MST To: Phil Stauffer , "Terry A. Miller" 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 Date: January 31, 2011 9:44:12 AM MST To: Phil Stauffer , "Terry A. Miller" Cc: Rodney Harold Debruin , "Ronald C. Surdam" 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/