Description
Mesh, properties, initial conditions, injection/withdrawal rates for modelling thermal, hydrological, and mechanical effects of fluid injection to and withdrawal from ground for Stockton University reservoir cooling system (aquifer storage cooling system), Galloway, New Jersey, for unscheduled two hour injection at 133 % designed capacity, on fine scale grid, with some results. Second simulation of J.T. Smith, E. Sonnenthal, P. Dobson, P. Nico, and M. Worthington, 2021. Thermal-hydrological-mechanical modeling of Stockton University reservoir cooling system, Proceedings of the 46th Workshop on Geothermal Reservoir Engineering, Stanford University, SGP-TR-218, from which Figures 6-9, pertain.
Resources
| Name | Format | Description | Link |
|---|---|---|---|
| displacement.03.tec.zip | ZIP | (Zipped) TReactMech mechanics displacement output file (after two hours injection). (a tecplot file, text file.) element corner x,y,z (m), x displacement (m), y displacement (m), z displacement (m). Output points are at element corners. Points are in x,y,z order, i.e., by horizontal planes, x varying fastest. | Open resource |
| File Overview.txt | Text | An over view of files in this submission. | Open resource |
| flowdata003.tec.zip | ZIP | (Zipped) TReactMech flow data output file (after two hours injection). (a tecplot file, text file.) Element center x,y,z, pressure, temperature, gas saturation, liquid saturation, other, other, assumed capillary pressure, gas density, liquid density, porosity, xx permeability, yy permeability, zz permeability, vertical permeability/ starting vertical permeability, assumed gas relative permability, assumed liquid relative permeability, specific heat capacity, thermal conductivity, bulk density, diffusivity. Units are m, m, m, Pa, deg. C, none, none, none, none, Pa, kg/m^3, kg/m^3, none, m^2, m^2, m^2, none, none, none, J/m/deg. C, kg/m^3, m^2 /s. Points are in x,y,z order, i.e., by horizontal planes, x varying fastest. | Open resource |
| flow.stockton.8layloc.inp | Unknown | Program TReactMech flow.inp input file (a text file). Global and local flow parameters. Local parameters organized by (MESH file) rock type. Formats for individual items described in TOUGH2 User's Guide, Version 2, appendix E. All quantities in S.I. (MKS) units. | Open resource |
| GENER.2well.stresstest.halfsrcs | Unknown | TReactMech GENER file (a text file). Fluid source and sink terms in TOUGH2 GENER file format. Sources and sinks are separate entries. For each source after line with element number; table of times (s), table of injection rates (kg/s), table of specific enthalpies (J/kg). For each sink after line with element number; table of times (s), table of (negative) injection rates (kg/s). Simulation injection/widthdrawl schedule time 0 at beginning of a November. Format described in TOUGH2 User's Guide, appendix E. Injection rates are 33.648 kg/s/well in two wells in the modelled domain (before reflections) and 33.648 kg/s in the remaining implied cold well (two positive plus one negative reflection, see paper Smith, et al., 2021), that is, 100.94 kg/s total. | Open resource |
| initialstress_sand0.5029_clay0.7003_1.zip | ZIP | (zipped) TReactMech initial stress file (a text file). Format described in file TReactMech_mech_3D_inp_format.txt. Units are Pa. | Open resource |
| mech_3D.stockton.7lay.notensclay.updated.inp | Unknown | Program TReactMech mech_3D.inp input file, (a text file). Global and local mechanics parameters. Local parameters organized by (MESH file) rock type. Structure described in file TReactMech_mech_3D_inp_format.txt All quantities in S.I. (MKS) units. | Open resource |
| MESH.stoctonloc.3.1.noflowymax.zip | ZIP | (Zipped) TReactMech MESH file. Grid information in TOUGH2 MESH file format. Elements are in x,z,y order; x varying fastest, z next, and y last; i.e., in successive vertical slices normal to the y direction. Lines after "ELEME" have; element name (characters 1-5), rock type (char. 16-20), volume (char 21-30), center x, y, z locations, (char 51-60,61-70,71-80). Lines after "CONNE" have connection information. Lines after "+++" have connection element number pairs (reduplicative). Formats described in TOUGH2 User's Guide, appendix E. Coordinates are approximately x = distance NE from cold well 1, y= distance NW from cold well 1, z = elevation (topography neglected). Volume stated for constant P T boundary condition cells is 10^50 times actual volume (e.g., see TOUGH2 Users Guide). | Open resource |
| Plasticity.tec | Unknown | TReactMech mechanical failure output file, at all times. (a text file.) At each time, a header, then X(m), Y(m), Z(m), ielem, N_shear, N_tensl, N_Surf, E_fail_xx, E_fail_yy, E_fail_zz, 2*E_fail_yz, 2*E_fail_xz, 2*E_fail_xy, E_fail_vol, at each mechanically failed (fractured) element. Output points are at element centers. Variables are; center x,y,z position (m), element number, maximum number of independent planes of shear failure, maximum number of independent planes of tensile failure, maximum number of independent failure planes, failure strain, volumetric failure strain. Failure strains are averaged over the 8 Gauss points of an element. Maximum numbers are over the 8 Gauss points. Only elements that have failed up to the reported time step are reported. | Open resource |
| SAVE.1.67e12s.zip | ZIP | (Zipped) TReactMech INCON file (a text file) (TOUGH2 compatible). Initial pressure and temperature by element name. File alternates lines of element name (with porosity and perm_xx, perm_yy, perm_zz), and pressure (Pa) and temperature (deg. C). | Open resource |
| stress_strain.03.tec.zip | ZIP | (Zipped) TReactMech mechanics stress and strain output file (after two hours injection). (a tecplot file, text file.) "X(m)", "Y(m)", "Z(m)", "Sigma_xx", "Sigma_yy", "Sigma_zz", "Sigma_yz", "Sigma_xz", "Sigma_xy", "Strain_xx", "Strain_yy", "Strain_zz", "Strain_yz", "Strain_xz", "Strain_xy", "E_fail_xx", "E_fail_yy", "E_fail_zz", "2*E_fail_yz","2*E_fail_xz","2*E_fail_xy","E_fail_vol", "Poro_elast", "Press". (center x,y,z position (m), total stress (Pa), total strain (non-dimen.), failure strain (non-dimen.), poroelastic porosity (non-dimen.), fluid pressure (Pa).) Output points are at element centers, stresses and strains are averaged over Gauss points. Points are in x,y,z order, i.e., by horizontal planes. | Open resource |
| TReactMech_mech_3D_inp_format.txt | Text | Program TReactMech mech_3D.inp input file keyword format descriptions. | Open resource |
| T-H-M Modeling of Stockton Univ. Reservoir Cooling System | Unknown | J.T. Smith, E. Sonnenthal, P. Dobson, P. Nico, and M. Worthington, 2021. Thermal-hydrological-mechanical modeling of Stockton University reservoir cooling system, Proceedings of the 46th Workshop on Geothermal Reservoir Engineering, Stanford University, SGP-TR-218. | Open resource |
| TOUGH2 User's Guide Version 2 | K. Pruess, C. Oldenburg, and G. Moridis, 2012. TOUGH2 User's Guide, Version 2, Lawrence Berkeley National Laboratory, LBNL-43134. | Open resource |