Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “subsurface flow”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Dynamic Surface Incorporation of Pb 2+ Ions at the Actively Dissolving Calcite (104) Surface

The reaction of dissolved Pb 2+ with calcite surfaces at near-equilibrium conditions involves adsorption of Pb 2+ and precipitation of secondary heteroepitaxial Pb-carbonate minerals. A more complex behavior is observed under far-from-equilibrium conditions, including strong inhibition of calcite dissolution, development of microtopography, and near-surface incorporation of multiple monolayers (ML) of Pb 2+ without precipitation of secondary phases [where 1 ML ≡ 1 Ca/20.2 Å 2 , the crystallographic site density of the calcite (104) lattice plane]. However, the mechanistic controls governing far-from-equilibrium reactivity are not well understood. Here, in this study, we observe the interfacial incorporation of dissolved Pb 2+ during the dissolution of calcite (104) surfaces at pH ~3.7 in a flow-through reaction cell, revealing the formation of a ~1 nm thick Pb-rich calcite layer with a total Pb coverage of ~1.4 ML. These observations of the sorbed Pb distribution used resonant anomalous X-ray reflectivity, X-ray fluorescence, and nanoinfrared atomic force microscopy. We propose that this altered surface layer represents a novel sorption mode that is stabilized by conditions of sustained disequilibrium. This behavior may significantly impact the transport of dissolved metals during disequilibrium processes occurring in acid mine drainage and subsurface CO 2 injection and, if appropriately accounted for, could improve the predictive capability of geochemical reactive-transport models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Transient Pressure Interference During CO2 Injection in Saline Aquifers

Abstract CO2 injection in subsurface geological formations (e.g. deep saline aquifers) causes pressure perturbations over a large area surrounding the injection well. Observation wells are widely considered in geologic CO2 storage (GCS) projects where the pressure perturbation induced by CO2 injection is measured. In this work, we use analytical and numerical modeling tools along with field data to examine the pressure behavior in GCS projects before and after CO2 arrival at an observation well. Prior to CO2 arrival, a baseline pressure trend is established which corresponds to single-phase brine flow across the observation well (approximated by Theis solution). Therefore, analysis of early-time pressure data is straightforward, provides the single-phase flow characteristics (mobility and storativity), and helps establishing a baseline pressure change that can be extended beyond the single-phase flow period at the observation well. Upon CO2 arrival, a departure from this baseline trend is expected. For the pressure to detect the CO2 arrival at an observation well, the departure from baseline pressure behavior must be significant and well above the background noise levels. We use existing analytical models to determine the strength of the expected pressure departure signal from the baseline trend upon CO2 arrival. The strength of the expected pressure departure is found to be directly proportional to the mobility ratio. Accordingly, we establish a criterion to determine whether the pressure at an observation well can detect the CO2 arrival. We present an analysis approach through application to synthetic and field data and show the characteristic pressure behavior before and after CO2 arrival. We show that while generally the pressure can be either above or below the expected baseline pressure trend, it would be likely above the baseline upon CO2 arrival. This is because the mobility ratio becomes less than unity after CO2 arrival. We show that depending on the reservoir characteristics, changes in the pressure trend may or may not be sufficient to detect the CO2 arrival.

Engineering↗
NASA NTRS2010

Tectonics of the Outer Planet Satellites

Tectonic features on the satellites of the outer planets range from the familiar, such as clearly recognizable graben on many satellites, to the bizarre, such as the ubiquitous double ridges on Europa, the twisting sets of ridges on Triton, or the isolated giant mountains rising from Io's surface. All of the large and middle-sized outer planet satellites except Io are dominated by water ice near their surfaces. Though ice is a brittle material at the cold temperatures found in the outer solar system, the amount of energy it takes to bring it close to its melting point is lower than for a rocky body. Therefore, some unique features of icy satellite tectonics may be influenced by a near-surface ductile layer beneath the brittle surface material, and several of the icy satellites may possess subsurface oceans. Sources of stress to drive tectonism are commonly dominated by the tides that deform these satellites as they orbit their primary giant planets. On several satellites, the observed tectonic features may be the result of changes in their tidal figures, or motions of their solid surfaces with respect to their tidal figures. Other driving mechanisms for tectonics include volume changes due to ice or water phase changes in the interior, thermoelastic stress, deformation of the surface above rising diapirs of warm ice, and motion of subsurface material toward large impact basins as they fill in and relax. Most satellites exhibit evidence for extensional deformation, and some exhibit strike-slip faulting, whereas contractional tectonism appears to be rare. Io s surface is unique, exhibiting huge isolated mountains that may be blocks of crust tilting and foundering into the rapidly emptying interior as the surface is constantly buried by deposits from hyperactive volcanoes. Of the satellites, diminutive Enceladus is spectacularly active; its south polar terrain is a site of young tectonism, copious heat flow, and tall plumes venting into space. Europa's surface is pervasively tectonized, covered with a diverse array of exotic and incompletely understood tectonic features. The paucity of impact craters on Europa suggests that its tectonic activity is ongoing. Geysers on Triton show that some degree of current activity, while tectonic features that cross sparsely cratered terrain indicate that it may also be tectonically active. Ganymede and Miranda both exhibit ancient terrains that have been pulled apart by normal faulting. On Ganymede these faults form a global network, while they are confined to regional provinces on Miranda. Ariel, Dione, Tethys, Rhea, and Titania all have systems of normal faults cutting across their surfaces, though the rifting is less pronounced than it is on Ganymede and Miranda. Iapetus exhibits a giant equatorial ridge that has defied simple explanation. The rest of the large and middle-sized satellites show very little evidence for tectonic features on their surfaces, though the exploration of Titan's surface has just begun.

McKinnon, W. B.↗
DOE OSTI2021

Origin of water masses in Floridan Aquifer System revealed by 81 Kr

Groundwater in coastal aquifers serves as an essential resource in densely populated areas throughout the world. However, hydrological connection to the ocean leaves it at risk of salinization via changes in climate and hydrological cycles. Therefore, an accurate hydrological characterization of coastal aquifers is of the utmost importance. Although localized salinization of shallow aquifers has commonly been studied, geochemical constraints on the time scale of regional freshwater and seawater in coastal aquifers are limited especially in carbonate aquifers. Such information is not only crucial for water resource management, but also for understanding solute flux over the land-ocean boundary and for interpreting the geochemical signatures that each water mass carries. Toward this goal, we report tracer-based, subsurface residence times of groundwater in the Floridan Aquifer System in southern Florida, which is one of the most productive aquifers on Earth. The first application of Kr 81-an emerging tracer of old groundwater-to this aquifer has identified freshwater recharge during the last glacial period as anticipated from previous studies using other geochemical tools including radiocarbon, stable isotopes and noble gas concentrations. Moreover, freshwater of Holocene age was also detected at 100 km-distance from the recharge area, suggesting the possibility of an active flow system in the upgradient region. A contribution of fossil seawater that predates the last glacial maximum was locally identified, suggesting slow seawater circulation and limited but significant solute flux from the ocean into the aquifer for dolomitization. This also implies that the deep saline coastal groundwater can potentially serve as a paleo-seawater archive.

58 GEOSCIENCES↗
NASA NTRSDate not supplied

Imdr Regio: The Geology of A Possible Active Hot Spot on Venus and A Target for Future Exploration By the ESA Envision Mission.

With Venus being the focus of future new missions like ESA’s EnVision, we need to look for target sites for future investigations on Venus. The different instruments will help to study the planet in a unprecedented detail. In particular, the VenSAR instrument will map the 20% of the surface of Venus at a resolution of 5 to 10 meters/pixel, so it gets crucial to select what areas of Venus deserve to be targeted for such a high-resolution mapping. We propose here Imdr Regio as a good candidate location for combined high resolution SAR mapping and for study with the rest of the mission instruments (e.g., Subsurface Radar Sounder). Imdr Regio is a large topographic rise that extends approximately from 35ºS-50ºS and 195ºE-225ºE between the volcanic plains of Helen, Nsomeka and Wawalag Planitiae. Imdr Regio was classified as a volcano-dominated topographic rise with a minimum-maximum diameter of 1200-1400 km and a swell height of 1.6 km. The southeast of Imdr Regio is dominated by Idunn Mons (46.5ºS/214.5ºE). Studies using the 1 µm (derived) surface emissivity from the Venus Express mission of the volcanic flows surrounding Idunn Mons suggest that high emissivity values in the volcano flanks are related to relatively unweathered basaltic rock and therefore indicative of a recent or even ongoing volcanic activity. We have carried out geologic mapping in the area to constrain volcanic and tectonic structures, and the geologic history of the large topographic rise. This geologic mapping reveals that different styles of volcanism are present across Imdr Regio and that this volcanic activity takes place in all the hot spot in close relation with the formation of a rift throughout all time represented by in the area. The first stages on the evolution of this large volcano are characterized by the formation of a radial fracture system and contemporaneous large sheet flows. These large sheet flows are locally difficult to distinguish from regional plains and are also deformed by regional N-trending wrinkle ridges. After this initial phase, multiple overlapping digitate flow units form the flanks and summit of the volcano and are contemporaneous with NW-SE trending fractures and graben of Olapa Chasma, a rift system that cross the topographic rise. To the northwest another large volcano also presents a system of radial fractures but lacks clear large sheet flows. Numerous pit chains are associated to fractures and graben, suggesting that transport of magma under the surface is important. The geologic history based on the mapping suggests that activity in Imdr Regio started with the emplacement of a plume/diapir in the southeast, which resulted in the formation of a radial dyke system and the emplacement of large sheet flows in Idunn Mons. The presence of other large volcanoes to the northwest suggests the presence of another magmatic source (plume or diapir). After this initial stage, volcanism continues in Idunn Mons and activity in the topographic rise is strongly related to the formation of Olapa Chasma. Rift-related volcanic flows postdate these other volcanoes and are contemporaneous with the late activity in Idunn Mons. Monogenetic volcanism is present in the rift and in all the units that are in the area.

I. López↗
DOE OSTI2021

Injection data analysis using material balance time for CO 2 storage capacity estimation in deep closed saline aquifers

Estimating the ultimate storage capacity of deep saline aquifers is important to address the formation potential to store the envisioned large volumes of CO 2 . Injection data (i.e. injection rate, bottomhole pressure, and cumulative injected volume of CO 2 ) are routinely recorded during storage operations. These data contain valuable information on the subsurface (e.g. the reservoir pore volume and the formation storage capacity) that can be extracted. In this paper, we present a two-step graphical technique to infer the pore volume and the ultimate storage capacity of closed saline aquifers by analyzing the available injection data. First, the pore volume is inferred through adapting the concept of the material balance time. Material balance time is an approximate superposition time function developed to interpret production data from oil and gas wells operating at variable pressure/rate conditions during the boundary-dominated flow period. Using material balance techniques, the ultimate storage capacity is then estimated through linear extrapolation of the average pressure trend to the maximum allowable pressure the formation can withstand. The average pressure is not available in practice, but is can be obtained from the injection data. Two approaches are presented in this study to calculate the average pressure; namely the rigorous and the approximate approaches. Unlike the rigorous approach, the approximate approach does not require a prior knowledge of some reservoir properties (e.g. relative permeability, absolute permeability, formation porosity and thickness) to calculate the average pressure. To investigate its potential and reliability in analyzing CO 2 injection data, the proposed technique is applied to four synthetic cases representing different well operating conditions. Results indicate that the approximate approach consistently overestimates the actual (simulated) storage capacity as compared to the rigorous approach. The agreement - between the inferred and the simulated reservoir pore volume, and between the analytical and numerical estimates of storage capacity - validates the potential application of the technique to CO 2 storage in closed saline aquifers. The technique is further substantiated through application to a field data set utilized from a commercial-scale geological storage (CGS) project. Finally, field data interpretation shows that the proposed technique can be utilized to identify the degree of hydraulic continuity and reservoir compartmentalization within a target formation by interpreting the corresponding pressure and rate responses.

02 PETROLEUM↗
DOE OSTI2020

Creation of a Mixed-Mode Fracture Network at Mesoscale Through Hydraulic Fracturing and Shear Stimulation

Enhanced Geothermal Systems could provide a substantial contribution to the global energy demand if their implementation could overcome inherent challenges. Examples are insufficient created permeability, early thermal breakthrough, and unacceptable induced seismicity. In this study we report on the seismic response of a mesoscale hydraulic fracturing experiment performed at 1.5-km depth at the Sanford Underground Research Facility. We have measured the seismic activity by utilizing a 100-kHz, continuous seismic monitoring system deployed in six 60-m length monitoring boreholes surrounding the experimental domain in 3-D. The achieved location uncertainty was on the order of 1 m and limited by the signal-to-noise ratio of detected events. These uncertainties were corroborated by detections of fracture intersections at the monitoring boreholes. Three intervals of the dedicated injection borehole were hydraulically stimulated by water injection at pressures up to 33 MPa and flow rates up to 5 L/min. We located 1,933 seismic events during several injection periods. The recorded seismicity delineates a complex fracture network comprised of multistrand hydraulic fractures and shear-reactivated, preexisting planes of weakness that grew unilaterally from the point of initiation. We find that heterogeneity of stress dictates the seismic outcome of hydraulic stimulations, even when relying on theoretically well-behaved hydraulic fractures. Once hydraulic fractures intersected boreholes, the boreholes acted as a pressure relief and fracture propagation ceased. In order to create an efficient subsurface heat exchanger, production boreholes should not be drilled before the end of hydraulic stimulations.

58 GEOSCIENCES↗
DOE OSTI2026

Numerical Analysis of Offshore Wind-Farm-Induced Drag Effects on Coastal Upwelling Dynamics

Wind farms extract momentum from the atmospheric flow, generating wind-speed deficits both within the plant, and extending downstream. When located offshore, these deficits modulate air-sea coupling, potentially impacting coastal upwelling in sensitive regions. We investigate impacts of wind farms on coastal upwelling using kilometre-scale, three-way-coupled simulations with the coupled ocean-atmosphere-wave-sediment transport system for the US West Coast. Wind-farm effects are represented by a generalised turbine drag formulation, an idealised, height-dependent body force whose magnitude is systematically varied. This approach isolates the leading-order fluid-dynamical response in a realistic coastal configuration. The atmospheric adjustment exhibits an approximately linear relation between drag force and wind-speed deficit, with wakes that expand downstream and increase in magnitude as drag increases. An empirical orthogonal function analysis of sea-surface-temperature anomalies reveals the emergence of a canonical dipole pattern under strong drag forcing. Subsurface diagnostics show consistent shoaling of the mixed layer and suppressed upward velocities in areas near wind-farm region, accompanied by compensating enhancements of shoaling closer to the coast. These results identify turbine drag as a control parameter in assessing interactions between wind-farm wake and coastal upwelling and provide scaling relationships for understanding offshore wind-farm effects on the coastal circulation dynamics.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Pore-Scale Analysis of the Permeability Damage and Recovery during Cyclic Freshwater and Brine Injection in Porous Media Containing Non-Swelling Clays

Permeability damage in subsurface porous media caused by clay mobilization is encountered in many engineering applications, such as geothermal energy, water disposal, oil recovery, and underground CO2 storage. During the freshwater injection into rocks containing brine, the sudden decrease in salinity causes native clay fines to detach and clog pore throats, leading to a significant decline in permeability. The clay fines detach due to weakened net-attractive forces binding them to each other and the grain. Past experiments link this permeability damage on the immediate history of the salinity and the direction of flow. To better understand this phenomenon, we conducted pore-scale simulations of cyclic injection of freshwater and brine into sandstone containing Kaolinite clay. Our simulations establish a link between the clay-fine trajectory and the permeability trend observed by Khilar and Fogler (1983). For a uniform clay size of 3 microns, we observe a permeability decline by two orders of magnitude during freshwater injection with respect to brine injection. Increasing salinity and simultaneously reversing flow direction restores the permeability. The permeability restoration upon reversing the brine flow direction is attributed to the unblocking of pore throats in the reverse direction by the movement of the clay particles along the grain surfaces by the hydrodynamic force and the strong net-attractive force under high salinity.

Bhuvankar, Pramod (ORCID:000000027739973X)↗
NASA NTRS1991

Large area QNDE inspection for airframe integrity

Quantitative Nondestructive Evaluation (QNDE) technology is being developed to provide new options for cost effective inspection of airframes. An R&D effort based on five NDE technologies is addressing questions of structural bonding assessment, corrosion detection, multisite damage detection, and fatigue characterization. The research/applications are being conducted by prioritized focussing and staging of the following technologies: (1) thermal NDE; (2) ultrasonic NDE; (3) coherent optical NDE; (4) magnetic imaging NDE; and (5) radiographic NDE. The focus here is on the most recent applications of thermal NDE technology to large area inspection of lap-joint and stiffener bonds. The approach is based on pulsed radiant heating of the airframe and measurement of the surface temperature of the structure with an infrared imager. Several advantages of the technique are that it is noncontacting, inspects one square meter area in a period of less than 2 minutes and has no difficulty inspecting typical curvatures of the fuselage. Numerical models of heat flow in these geometries are used to determine appropriate techniques for reduction of the infrared images, thereby delineating regions of disbonds. These models are also used to determine the optimum heating and measurement times for maximizing the contrast between bonded and unbonded structures. Good agreement is found between these results and experimental measurements, and a comparison of the two are presented. Also presented are results of measurements on samples with fabricated defects which show the technique is able to clearly indicate regions of disbonds. Measurements on an airframe also clearly image subsurface structure.

Winfree, William P.↗
DOE OSTI2024

Lagrangian Decomposition of the Atlantic Ocean Heat Transport at 26.5°N

Abstract The Atlantic Meridional Overturning Circulation (AMOC) plays a critical role in the global climate system through the redistribution of heat, freshwater and carbon. At 26.5°N, the meridional heat transport has traditionally been partitioned geometrically into vertical and horizontal circulation cells; however, attributing these components to the AMOC and Subtropical Gyre (STG) flow structures remains widely debated. Using water parcel trajectories evaluated within an eddy‐rich ocean hindcast, we present the first Lagrangian decomposition of the meridional heat transport at 26.5°N. We find that water parcels recirculating within the STG account for 37% (0.36 PW) of the total heat transport across 26.5°N, more than twice that of the classical horizontal gyre component (15%). Our findings indicate that STG heat transport cannot be meaningfully distinguished from that of the basin‐scale overturning since water parcels cooled within the gyre subsequently feed the northward, subsurface limb of the AMOC.

Tooth, Oliver J.↗
DOE OSTI2024

Modeled estimates of permafrost hydrology and related fields for pan-Arctic region over the period 1980-2100

This dataset contains gridded estimates of total runoff (mm/day); subsurface runoff (mm/day); snow water equivalent (SWE, mm/day); vertical upward water flux of combined surface evaporation, evapotranspiration, and snow sublimation (mm/day); and soil temperature (C) across the pan-Arctic basin from simulations of the Permafrost Water Balance Model (PWBM v4). These data were generated to allow researchers to advance understanding of the changes unfolding in the Arctic's terrestrial water cycle under the influences of climate warming, hydrological cycle intensification, and permafrost thaw. The dataset will enable researchers to address important questions spanning multiple disciplines. For example, how will rising river flows exported to the coast impact critical zone ecosystems, near-shore ice dynamics, and ocean biogeochemical processes. Spatial resolution is 25x25 km with N=39926 grid cells spanning the domain. Three forcing datasets were used to drive simulations at an implicit daily time step over the period 1980-2100. Daily air temperature, precipitation, and wind speed data from the WATCH Forcing Data methodology applied to ERA5 reanalysis (W5E5) were used for meteorological forcings over the period 1980-2019. Simulations for the period 1980-2100 used meteorological data from two global climate models; IPSL-CM6A-LR and MPI-ESM1-2-HR. The PWBM hydrological outputs are organized by grid and month in gzipped ASCII text files within each archived tar file. Soil temperature is organized by grid cell and day of year from DOY 205 to DOY 260. No special software is required to work with these data. The package contains spatial domain and river basin information files. These data are most useful for analyses of the spatial and temporal variations in terrestrial hydrology and soil thermal dynamics. A manuscript describing the model simulations and associated analysis has been accepted for publication in The Cryosphere (Rawlins and Karmalkar, 2024).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Mechanisms for Microseismicity Occurrence Due to CO 2 Injection at Decatur, Illinois: A Coupled Multiphase Flow and Geomechanics Perspective

Here, we numerically investigate the mechanisms that resulted in induced seismicity occurrence associated with CO 2 injection at the Illinois Basin–Decatur Project (IBDP). We build a geologically consistent model that honors key stratigraphic horizons and 3D fault surfaces interpreted using surface seismic data and microseismicity locations. We populate our model with reservoir and geomechanical properties estimated using well-log and core data. We then performed coupled multiphase flow and geomechanics modeling to investigate the impact of CO 2 injection on fault stability using the Coulomb failure criteria. We calibrate our flow model using measured reservoir pressure during the CO 2 injection phase. Our model results show that pore-pressure diffusion along faults connecting the injection interval to the basement is essential to explain the destabilization of the regions where microseismicity occurred, and that poroelastic stresses alone would result in stabilization of those regions. Slip tendency analysis indicates that, due to their orientations with respect to the maximum horizontal stress direction, the faults where the microseismicity occurred were very close to failure prior to injection. These model results highlight the importance of accurate subsurface fault characterization for CO 2 sequestration operations.

58 GEOSCIENCES↗
DOE OSTI2024

Creating microfluidic channels functionalized with micro- and nano-scale features via femtosecond laser surface processing

We demonstrate the ability to create functionalized microfluidic channels using femtosecond laser surface processing (FLSP). FLSP is an emerging advanced manufacturing technology used to modify the surface properties of materials directly and permanently by producing self-organized quasi-periodic micro- and nano-scale surface features along with surface and subsurface chemical and grain structure changes. We demonstrate on Hastelloy X that by controlling the laser fluence and pulse count, the depth of the microchannels and height of the FLSP microstructures within the microchannels can be controlled independently.

Ultrashort pulse laser applications, micro- and na↗
DOE OSTI2025

Corrosion of 316 L stainless steel under the natural circulation of molten NaCl-MgCl 2 salt

The corrosion behavior of 316 L stainless steel (SS) was studied via the natural circulation of molten eutectic NaCl-MgCl 2 salt through a microloop. The post-corrosion tested 316 L SS microloop sections were characterized with microscopy techniques to determine the microstructural and microchemical changes that occurred at the alloy/salt interface. It was found that 316 L SS showed heterogeneous dissolution at the hot-leg, whereas deposition of corrosion products occurred at the cold-leg. For the first time, experimentally obtained molten salt flow-induced corrosion of 316 L SS results were combined with computational thermodynamic-kinetic models to validate the dissolution and deposition in terms of elemental compositional changes at the alloy/salt interface. The thermodynamic-kinetic modeling predicted that the heterogeneous dissolution of Cr and Fe from the hot-leg section of 316 L SS persisted throughout the salt circulation. The model also estimated that, despite Cr deposition starting earlier than Fe, the total redeposition of Fe is expected to be significantly greater than that of Cr over the circulation of salt. Furthermore, the modeling accurately predicted the subsurface enrichment of Mo which is attributed to the reduced Cr activity and the relatively higher diffusion rate of Mo within the alloy matrix. Here, the agreement between modeling and experimental results confirms that Fe chlorides dissolve at the hot-leg and subsequently deposit at the cold-leg due to activity changes driven by the thermal gradient. By contrast, Cr was not detected in the cold-leg deposits, which is attributed to its weaker temperature dependence on activity, limiting its redeposition under these conditions.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University: Preprint

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ("CUBO") was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗
DOE OSTI2025

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ('CUBO') was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗