-------- Original Message -------- Subject: Dump file from Petrel Date: Tue, 3 Sep 2013 14:57:21 +0000 From: Zunsheng Jiao To: Philip H. Stauffer (stauffer@lanl.gov) Good morning Phil, Attached is file dumped from the Petrel. The number of the cells is 1,374,719, including I, j, and k number, X, y, and z coordinates, and zone id. If zone number is note defined (out of the model), it is assigned with -99.00. Could this file be used to generate FEHM files? Thanks, John Zunsheng Jiao, Chief Geologist Carbon Management Institute School of Energy Research, University of Wyoming Dept. 4302, 1000 University Avenue Laramie, WY 82071 (307) 766-6802 (O) (307)399-5745 (C) -rw-r--r-- 1 tam sftus 84020735 Sep 3 12:32 test_grid_Gslib.dat Regular grid 139x230x43 with null -99.00 x,y around irregular shape Spacing is 2500x2500 in x,y with variable z (spacing near 60-125 vertical). The outline boundary is irregular and can not be captured with a simple connect mesh method. With more work quad surfaces can be created from this point set and stack methods can be used. Note that top and bottom layers are not flat. Top layer has variable topology that pinches out some lower layers. The points still exist but are squished together close to top surface which causes neg ccoef in these areas. The extreme aspect ratio of xy to z combined with sloping interfaces causes difficulties with Delaunay connections such that nodes may be removed. The material interfaces are not guaranteed where slope is steep enough to "pop" elements up or down during connect. ==> test_grid_Gslib.dat <== PETREL: Properties 7 i_index unit1 scale1 j_index unit1 scale1 k_index unit1 scale1 x_coord unit1 scale1 y_coord unit1 scale1 z_coord unit1 scale1 Zone_ID unit1 scale1 1 1 43 -99.00 -99.00 -99.00 -99.00 ... 48 1 43 -99.00 -99.00 -99.00 -99.00 49 1 43 226050.00000000 3869700.00000000 -2362.99603271 1.000000 50 1 43 228550.00000000 3869700.00000000 -2339.46075439 1.000000 ... 41 2 43 -99.00 -99.00 -99.00 -99.00 42 2 43 208550.00000000 3872200.00000000 -2534.86547852 1.000000 43 2 43 211050.00000000 3872200.00000000 -2518.70272827 1.000000 ... 137 230 1 -99.00 -99.00 -99.00 -99.00 138 230 1 -99.00 -99.00 -99.00 -99.00 139 230 1 -99.00 -99.00 -99.00 -99.00 END readatt: reading 5 values from 1102042 lines. ATTRIBUTE NAME MIN MAX DIFFERENCE LENGTH xic 1.060500000E+05 4.510500000E+05 3.450000000E+05 1102042 yic 3.869700000E+06 4.442200000E+06 5.725000000E+05 1102042 zic -4.021226013E+03 1.698531097E+03 5.719757111E+03 1102042 xval 4.987690000E-01 1.600361180E+02 1.595373490E+02 1102042 nlayer 1 43 42 1102042 -------- Original Message -------- Subject: RE: mesh from gslib points Date: Wed, 4 Sep 2013 22:39:58 +0000 From: Zunsheng Jiao To: Stauffer, Philip Phil, We could start with a 50 km X 50 km domain. Four corners for this domain: SW x240000 y4043000 SE x290000 y4043000 NE x290000 y4093000 NW x240000 y4093000 Thanks, John From: Stauffer, Philip [mailto:stauffer@lanl.gov] Sent: Wednesday, September 04, 2013 1:56 PM To: Zunsheng Jiao Subject: Fwd: mesh from gslib points John, Terry is working on it. there are some issues. How much of the domain do you need ? If we trim some of the edges it will be much better computationally. Not sure how eclipse computes on such a bad mesh. Not voronoi at all in some places. P -------- Original Message -------- Subject: mesh from gslib points Date: Wed, 4 Sep 2013 13:51:30 -0600 From: Miller, Terry A To: Stauffer, Philip H Hi Phil, I have created awk and lagrit scripts to correct point values, subset, then connect into an FEHM mesh. This is a NS slice, 2 elements wide. 30x exaggeration in Z. I had to cut the north end off where layers go up to the surface and can not connect delaunay. This mesh has a few neg ccoef, but they are on the outside of the mesh so should be ok. I can make a big box shape subset of this in the morning. /scratch/fwo/tam/wy_gslib_v01/subset_NS