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.

This is the modeling data (input/output files of TOUGHREACT 4.10) used to simulate the reactive transport processes of the Aquifer Thermal Energy Storage (ATES) operations at Stockton University, NJ. Readme.txt lists all the files. TOUGHREACT 4.10 requires to reproduce the modeling output. The modeling data in this submission is related to the Aquifer Injection for Energy Storage purposes outlined in "Reactive Transport Modeling of Aquifer Thermal Energy Storage System at Stockton, NJ During Seasonal Operations".

This dataset contains input data, code, ReadMe files, output data, and figures that summarize the results of a stochastic analysis of geothermal reservoir production from two potential geothermal reservoirs that were evaluated for the Cornell University Deep Direct-Use project. These potential reservoirs are the Trenton-Black River (TBR) from 2.27-2.3 km depth, and basement rocks from 3.0-3.5 km depth and 3.5-4.0 km depth. Several utilization scenarios consisting of different injection fluid temperatures and flow rates were evaluated for each reservoir. Uncertainty in geologic properties, thermal properties, economic costs, and utilization efficiencies were evaluated using a Monte Carlo analysis of the reservoir simulations. Some reservoir simulations of the TBR were completed using the TOUGH2 software, as implemented in PetraSIM. The PetraSIM run files and associated data are provided with this submission. All other reservoir simulations were completed using the GEOPHIRES software, with some modifications to complete the uncertainty analyses. ReadMe files that describe additions to GEOPHIRES, the GEOPHIRES input data, and the output data are all provided, and references are provided to the code repository. Figures that summarize the reservoir heat production, temperature drawdown, and the probability of meeting targeted building heating demands with the produced heat and fluid temperatures are provided.

Preliminary geothermal reservoir simulations were performed using a homogeneous static model to evaluate and understand the effects of fluid and rock properties that could influence the delivery of thermal energy in a doublet system. A 5000 feet by 5100 feet by 500 feet homogeneous model having a constant porosity and permeability of 20% and 100 mD was used to perform preliminary geothermal reservoir simulations. The model was discretized in the x-, y-, and z-directions into 100, 101, and 100, gridblocks. Two wells were placed on the opposite ends of the central column of the discretized model. One of the wells was designated as a producer and the other an injector. Equal volumes of water was extracted and then injected into the reservoir via the production and injection wells. Water was extracted at a temperature of 109 deg F and re-injected at 50 deg F at the 1000 bbl/day. The files attached contains the input and output files of this simulation case. The input and some of the output files can be viewed in any text editor.

Preliminary geothermal reservoir simulations were performed using a homogeneous static model to evaluate and understand the effects of fluid and rock properties that could influence the delivery of thermal energy in a doublet system. A 5000 feet by 5100 feet by 500 feet homogeneous model having a constant porosity and permeability of 20% and 100 mD was used to perform preliminary geothermal reservoir simulations. The model was discretized in the x-, y-, and z-directions into 100, 101, and 100, gridblocks. Two wells were placed on the opposite ends of the central column of the discretized model. One of the wells was designated as a producer and the other an injector. Equal volumes of water was extracted and then injected into the reservoir via the production and injection wells. Water was extracted at a temperature of 109 deg F and re-injected at 50 deg F at the 1000 bbl/day. The files attached contains the input and output files of this simulation case. The input and some of the output files can be viewed in any text editor.

This is a link to the open access, published dataset and modeling that supports a feasibility study of Reservoir Thermal Energy Storage (RTES) in the Portland Basin, Oregon, USA. Citation: Burns, E.R., 2020, SUTRA model used to evaluate Saline or Brackish Aquifers as Reservoirs for Thermal Energy Storage in the Portland Basin, Oregon, USA: U.S. Geological Survey data release, https://doi.org/10.5066/P9A6D6XM.

,The HYDRUS program is a finite element model for simulating the one-dimensional movement of water, heat, and multiple solutes in variably saturated media. The program numerically solves the Richards' equation for saturated-unsaturated water flow and Fickian-based advection-dispersion equations for heat and solute transport. The flow equation incorporates a sink term to account for water uptake by plant roots. The heat transport equation considers conduction as well as convection with flowing water. The solute transport equations consider advective-dispersive transport in the liquid phase, and diffusion in the gaseous phase. The transport equations also include provisions for nonlinear and/or nonequilibrium reactions between the solid and liquid phases, linear equilibrium reactions between the liquid and gaseous phases, zero-order production, and two first-order degradation reactions: one which is independent of other solutes, and one which provides the coupling between solutes involved in sequential first-order decay reactions. The program may be used to analyze water and solute movement in unsaturated, partially saturated, or fully saturated porous media.,The flow region itself may be composed of nonuniform soils. Flow and transport can occur in the vertical, horizontal, or a generally inclined direction. The water flow part of the model can deal with (constant or time-varying) prescribed head and flux boundaries, boundaries controlled by atmospheric conditions, as well as free drainage boundary conditions. Soil surface boundary conditions may change during the simulation from prescribed flux to prescribed head type conditions (and vice-versa).,For solute transport the code supports both (constant and varying) prescribed concentration (Dirichlet or first-type) and concentration flux (Cauchy or third-type) boundary conditions. The dispersion coefficient includes terms reflecting the effects of molecular diffusion and tortuosity.,The unsaturated soil hydraulic properties are described using van Genuchten [1980], Brooks and Correy [1964] and modified van Genuchten type analytical functions. Modifications were made to improve the description of hydraulic properties near saturation. The HYDRUS code incorporates hysteresis by using the empirical model introduced by Scott et al. [1983] and Kool and Parker [1987]. This model assumes that drying scanning curves are scaled from the main drying curve, and wetting scanning curves from the main wetting curve. HYDRUS also implements a scaling procedure to approximate hydraulic variability in a given soil profile by means of a set of linear scaling transformations which relate the individual soil hydraulic characteristics to those of a reference soil.,Root growth is simulated by means of a logistic growth function. Water and salinity stress response functions can be defined according to functions proposed by Feddes et al. [1978] or van Genuchten [1987].,The governing flow and transport equations are solved numerically using Galerkin type linear finite element schemes. Integration in time is achieved using an implicit (backwards) finite difference scheme for both saturated and unsaturated conditions. Additional measures are taken to improve solution efficiency for transient problems, including automatic time step adjustment and adherence to preset ranges of the Courant and Peclet numbers. The water content term is evaluated using the mass-conservative method proposed by Celia et al. [1990]. Possible options for minimizing numerical oscillations in the transport solutions include upstream weighing, artificial dispersion, and/or performance indexing.,,

Low permeability geothermal reservoirs can be stimulated by hydraulic fracturing to create Enhanced (or Engineered) Geothermal Systems (EGS) with higher permeability and improved heat transfer to increase heat production. In this paper, we document our effort to develop a numerical simulator with explicit geomechanics-discrete flow network coupling by utilizing and further advancing the simulation capabilities of the Livermore Distinct Element Code (LDEC). The important modules of the simulator include an explicit finite element solid solver, a finite volume method flow solver, a joint model using the combined FEM-DEM capability of LDEC, and an adaptive remeshing module. The numerical implementation is verified against the classical KGD model. The interaction between two fractures with simple geometry and the stimulation of a relatively complex existing fracture network under different in-situ stress conditions are studied with the simulator.

This submission contains the source code of the Hydrothermal Finite Element Simulator used for the Treasure Island and UC Berkeley campus geothermal simulation. It contains a report that summarizes the development and validation of this Hydrothermal Finite Element Simulator, with a case study on Treasure Island site. It also contains a report that investigates the feasibility of upgrading the existing campus energy delivery system at UC Berkeley to a fifth-generation district heating and cooling system that includes geothermal heat/cold storage.

This submission contains the source code of the Hydrothermal Finite Element Simulator used for the Treasure Island and UC Berkeley campus geothermal simulation. It contains a report that summarizes the development and validation of this Hydrothermal Finite Element Simulator, with a case study on Treasure Island site. It also contains a report that investigates the feasibility of upgrading the existing campus energy delivery system at UC Berkeley to a fifth-generation district heating and cooling system that includes geothermal heat/cold storage.