Search NASA⌕ Search

SEARCH · Search NASA

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Hydrology Outweighs Temperature in Driving Production and Export of Dissolved Carbon in a Snowy Mountain Catchment

Abstract Terrestrial production and export of dissolved organic and inorganic carbon (DOC and DIC) to streams depends on water flow and biogeochemical processes in and beneath soils. Yet, understanding of these processes in a rapidly changing climate is limited. Using the watershed‐scale reactive‐transport model BioRT‐HBV and stream data from a snow‐dominated catchment in the Rockies, we show deeper groundwater flow averaged about 20% of annual discharge, rising to ∼35% in drier years. DOC and DIC production and export peaked during snowmelt and wet years, driven more by hydrology than temperature. DOC was primarily produced in shallow soils (1.94 ± 1.45 gC/m 2 /year), stored via sorption, and flushed out during snowmelt. Some DOC was recharged to and further consumed in the deeper subsurface via respiration (−0.27 ± 0.02 gC/m 2 /year), therefore reducing concentrations in deeper groundwater and stream DOC concentrations at low discharge. Consequently, DOC was primarily exported from the shallow zone (1.62 ± 0.96 gC/m 2 /year, compared to 0.12 ± 0.02 gC/m 2 /year from the deeper zone). DIC was produced in both zones but at higher rates in shallow soils (1.34 ± 1.00 gC/m 2 /year) than in the deep subsurface (0.36 ± 0.02 gC/m 2 /year). Deep respiration elevated DIC concentrations in the deep zone and stream DIC concentrations at low discharge. In other words, deep respiration is responsible for the commonly‐observed increasing DOC concentrations (flushing) and decreasing DIC concentrations (dilution) with increasing discharge. DIC export from the shallow zone was ~66% of annual export but can drop to ∼53% in drier years. Numerical experiments suggest lower carbon production and export in a warmer, drier future, and a higher proportion from deeper flow and respiration processes. These results underscore the often‐overlooked but growing importance of deeper processes in a warming climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Mechanistic understanding of carbon mineralization in fracture systems using microfluidics

Carbon mineralization in mafic and ultramafic rocks presents an opportunity for permanent carbon storage in the Earth's subsurface. However, due to their lower permeability, pre-existing fracture networks are key for mineralization to occur. Therefore, to fully develop this technology, a mechanistic understanding of the mineralization behavior in fractures with the consideration of hydrodynamic components is required. We use high-pressure microfluidics to investigate key mechanisms influencing dissolution–precipitation in a fracture network. The experiments were conducted in micromodels made of natural rocks with a comb-shaped flow channel to mimic a fracture network. This enabled studying the effect of injection rate on coupled dissolution–precipitation in advection and diffusion-dominated flow paths. We used gypsum carbonation as an analog reaction to allow for realistic experimental time frames due to its rapid reaction kinetics. The experimental work is coupled with high-fidelity numerical simulations to enhance our understanding of the parameters affecting the mineralization reaction. Our results demonstrate the importance of flow rate on the rate and nature of the gypsum carbonation reaction revealing that higher flow rates enable deeper penetration of the mineral precipitation front into the dead-end channels. This is an important finding since for sustained mineralization in a fracture network, precipitation in dead-ends while still allowing for flowing fractures is critical. Detailed characterization of the precipitates showed that lower flow rates led to porous and loose precipitates in the form of aragonite while higher flow rates mimicked supersaturation behavior leading to the formation of calcite. The reactive transport simulations further demonstrated the significance of flow velocity in advection-dominated channels to influence the efficiency of carbon mineralization in diffusion-dominated channels, potentially clogging of dead-end channels. These findings highlight the need for coupling chemical, mechanical, and hydrodynamic processes to evaluate the nature and extent of carbon mineralization in fractured media critical for permanent storage in mafic and ultramafic formations. This research further highlights the need for more investigation in potential subsurface fracture generation techniques to aid carbon mineralization.

25 ENERGY STORAGE↗
DOE OSTI2017

Measurement and simulation of heat exchange in fractured bedrock using inert and thermally degrading tracers

Multicomponent groundwater tracer tests were conducted in a well-characterized field site in Altona, NY using inert carbon-cored nanoparticles and a thermally degrading phenolic compound. Experiments were conducted in a mesoscale reservoir consisting of a single subhorizontal bedding plane fracture located 7.6 m below ground surface contained between two wells separated by 14.1 m. The reservoir rock, initially at 11.78C, was heated using 748C water. During the heating process, a series of tracer tests using thermally degrading tracers were used to characterize the progressive in situ heating of the fracture. Fiber-Optic Distributed Temperature Sensing (FODTS) was used to measure temperature rise orthogonal to the fracture surface at 10 locations. Recovery of the thermally degrading tracer’s product was increased as the reservoir was progressively heated indicating that the advancement of the thermal front was proportional to the mass fraction of the thermally degrading tracer recovered. Both GPR imaging and FODTS measurements reveal that flow was reduced to a narrow channel which directly connected the two wells and led to rapid thermal breakthrough. Computational modeling of inert tracer and heat transport in a two-dimensional discrete fracture demonstrate that subsurface characterization using inert tracers alone could not uniquely characterize the Altona field site. However, the inclusion of a thermally degrading tracer may permit accurate subsurface temperature monitoring. At the Altona field site, however, fluid-rock interactions appear to have increased reaction rates relative to laboratory-based measurements made in the absence of rock surfaces.

58 GEOSCIENCES↗
DOE OSTI2024

Utah FORGE 4-2492: Design and Implementation of Innovative Stimulation Treatments to Maximize Energy Recovery - 2024 Annual Workshop Presentation

This is a presentation on the Design and Implementation of Innovative Stimulation Treatments to Maximize Energy Recovery Efficiency by The University of Texas at Austin, presented by Mukul M. Sharma. This video slide presentation discusses the following objectives: (1) to place fractures uniformly in a horizontal well to ensure a uniform distribution of flow into the fractures; (2) to maximize the area of the created fracture; (3) to ensure connectivity of the created fracture network; and (4) to ensure fracture size is optimized. This presentation was featured in the Utah FORGE R&D Annual Workshop on August 14, 2024.

15 GEOTHERMAL ENERGY↗
DOE OSTI2014

Final report on "Modeling Diurnal Variations of California Land Biosphere CO2 Fluxes"

In Mediterranean climates, the season of water availability (winter) is out of phase with the season of light availability and atmospheric demand for moisture (summer). Multi-year half-hourly observations of sap flow velocities in 26 evergreen trees in a small watershed in Northern California show that different species of evergreen trees have different seasonalities of transpiration: Douglas-firs respond immediately to the first winter rain, while Pacific madrones have peak transpiration in the dry summer. Using these observations, we have derived species-specific parameterization of normalized sap flow velocities in terms of insolation, vapor pressure deficit and near-surface soil moisture. A simple 1-D boundary layer model showed that afternoon temperatures may be higher by 1 degree Celsius in an area with Douglas-firs than with Pacific madrones. The results point to the need to develop a new representation of subsurface moisture, in particular pools beneath the organic soil mantle and the vadose zone. Our ongoing and future work includes coupling our new parameterization of transpiration with new representation of sub-surface moisture in saprolite and weathered bedrock. The results will be implemented in a regional climate model to explore vegetation-climate feedbacks, especially in the dry season.

58 GEOSCIENCES↗
DOE OSTI2019

Numerical study of groundwater flow cycling controlled by seawater/freshwater interaction in Woodville Karst Plain

Seawater intrusion due to sea level rise and climate change could significantly contaminate coastal groundwater resources, particularly in Florida, the flat low-land state in the United States. Based on the field investigation and hydrological measurements, a three-dimensional SEAWAT model is developed to evaluate the groundwater flow cycling and seawater intrusion to freshwater system in the Woodville Karst Plain (WKP), a typical karst groundwater system in the Floridan aquifer. The karst conduit network in the aquifer acts as fast flow pathway for groundwater flow and solute transport, so seawater could deeply intrude into the aquifer. Wakulla Spring, an inland spring 17 km from the coast and a coastal submarine spring, Spring Creek Spring Complex are connected through the conduit network. The flow direction between the two springs switches under various rainfall conditions in this region, thus the discharges at two karst springs are used to estimate the location of seawater/freshwater mixing interface. The SEAWAT modeling results indicate that the mixing interface, defined as 2 PSU (Practical Salinity Unit), intrudes 3 to 5 km through the subsurface karst conduit during the dry season and severely contaminates nearly 1 km width of groundwater around the conduit. The salinity distribution and the distance of seawater intrusion through the conduit system are very sensitive to precipitation variation and the sea level boundary condition. Furthermore, predictions are made for seawater intrusion to the aquifer under various sea level rise, precipitation scenarios and water pumping. The results show that the seawater intrusion could reach and contaminate inland freshwater systems if sea level rises 1.0 m or during a long-term no-precipitation season. Finally, this study provides insights for modeling and predicting the vulnerability of a coastal karst aquifer through the simulation of variable-density flow.

58 GEOSCIENCES↗
DOE OSTI2020

Geothermal Deep Direct Use for Turbine Inlet Cooling in East Texas

The National Renewable Energy Laboratory (NREL), the Southern Methodist University Geothermal Laboratory (SMU), Eastman Chemical (Longview, TX), and TAS (Houston, TX) evaluated the feasibility of using geothermal heat to improve the performance of a natural-gas power plant in East Texas. The area of interest is the Eastman Chemical plant in Longview, Texas, which is on the northwestern margin of a geologic region known as the Sabine Uplift. The feasibility study focused on determining the potential for accessing a subsurface hot-water geothermal resource within a 10-km radius of the site to provide thermal energy for absorption chillers. Wells within a 20-km radius are included for broader geological comparison to determine the heat flow, temperature-at-depth, field porosity and permeability. The lithologies of most interest are the Lower Cretaceous Trinity Group and Upper Jurassic Cotton Valley Group. The deeper Cotton Valley formations are hotter (averaging 117 to 130°C), yet permeability and porosity are low. The shallower Trinity Group contains more variability in permeability and porosity and lower temperatures averaging about 98 to 117°C. The shallower formations are considered despite the lower temperature because of increased ability to produce larger volumes of water and extract enough heat before reinjection. The complete SMU analysis is available in the National Geothermal Data System (NGDS). Tapping such deep geothermal sources for direct heating (as opposed to power generation) is known as geothermal deep direct use (DDU). Geothermal DDU has potential across a wide swath of the United States but is underutilized due to challenging project economics associated with developing a deep geothermal resource for what are typically small-scale, variable-demand projects. This project examines the feasibility of geothermal energy integration in a natural-gas combined cycle power station in East Texas. The DDU resource is tapped to drive absorption chillers (24/7) for production of chilled water at 5-10°C (41-50°F). This chilled water is stored until needed, which allows for continuous operation with a relatively small-capacity geothermal/absorption chiller system. When conditions are favorable, the chilled water is dispatched to cool the air entering the compressor stage of a gas combustion turbine. This process, known as turbine inlet cooling (TIC), boosts power production during periods of high temperature and high-power demand. Such systems can enhance grid reliability and reduce the cost for peak-demand power. A simulation model of the power plant was developed in IPSEpro software and validated against operational data from the plant. This model allowed the team to estimate the additional power that could be produced by applying TIC under different operating and ambient conditions. Absorption chiller performance was estimated from vendor sources to determine the production rate of chilled water from the geothermal resource. Geothermal drilling and development costs were estimated using NREL's GEOPHIRES 2.0. The expected lower drilling costs in this region led to an estimated cost of geothermal heat of about $4/MMBtu (1.4 cents/kWh t ). The estimated cost for the absorption chillers and TIC hardware were obtained from literature sources and project partners. Hourly data were obtained for weather, natural gas and electricity prices, and plant operating state for 2017, which served as a representative year. NREL estimated the capital cost, operating cost, and additional electricity production and revenue for different combinations of geothermal capacity, chiller capacity, and water storage-tank size. The analysis drove toward smaller geothermal and chiller systems to reduce equipment cost. A relatively low-cost water storage tank accumulated the near-continuous chilled water output for later use when TIC was most valued.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

FY2022 Progress on Imbibition Testing in Containment Science

Estimation of two-phase fluid flow properties is important to understand and predict water and gas movement through the vadose zone for agricultural, hydrogeological, and engineering applications, such as for vapor-phase contaminant transport and/or containment of noble gases in the subsurface. In this second progress report of FY22, we present two ongoing activities related to imbibition testing on volcanic rock samples. We present the development of a new analytical solution predicting the temperature response observed during imbibition into dry samples, as discussed in our previous first progress report for FY22. We also illustrate the use of a multi-modal capillary pressure distribution to simulate both early- and late-time imbibition data collected on tuff core that can exhibit multiple pore types. These FY22 imbibition tests were conducted for an extended period (i.e., far beyond the time required for the wetting front to reach the top of the sample), which is necessary for parameter estimation and characterization of two different pore types within the samples.

42 ENGINEERING↗
DOE OSTI2015

Using field observations to inform thermal hydrology models of permafrost dynamics with ATS (v0.83)

Climate change is profoundly transforming the carbon-rich Arctic tundra landscape, potentially moving it from a carbon sink to a carbon source by increasing the thickness of soil that thaws on a seasonal basis. Thus, the modeling capability and precise parameterizations of the physical characteristics needed to estimate projected active layer thickness (ALT) are limited in Earth system models (ESMs). In particular, discrepancies in spatial scale between field measurements and Earth system models challenge validation and parameterization of hydrothermal models. A recently developed surface–subsurface model for permafrost thermal hydrology, the Advanced Terrestrial Simulator (ATS), is used in combination with field measurements to achieve the goals of constructing a process-rich model based on plausible parameters and to identify fine-scale controls of ALT in ice-wedge polygon tundra in Barrow, Alaska. An iterative model refinement procedure that cycles between borehole temperature and snow cover measurements and simulations functions to evaluate and parameterize different model processes necessary to simulate freeze–thaw processes and ALT formation. After model refinement and calibration, reasonable matches between simulated and measured soil temperatures are obtained, with the largest errors occurring during early summer above ice wedges (e.g., troughs). The results suggest that properly constructed and calibrated one-dimensional thermal hydrology models have the potential to provide reasonable representation of the subsurface thermal response and can be used to infer model input parameters and process representations. The models for soil thermal conductivity and snow distribution were found to be the most sensitive process representations. However, information on lateral flow and snowpack evolution might be needed to constrain model representations of surface hydrology and snow depth.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2015

Using field observations to inform thermal hydrology models of permafrost dynamics with ATS (v0.83)

Climate change is profoundly transforming the carbon-rich Arctic tundra landscape, potentially moving it from a carbon sink to a carbon source by increasing the thickness of soil that thaws on a seasonal basis. However, the modeling capability and precise parameterizations of the physical characteristics needed to estimate projected active layer thickness (ALT) are limited in Earth System Models (ESMs). In particular, discrepancies in spatial scale between field measurements and Earth System Models challenge validation and parameterization of hydrothermal models. A recently developed surface/subsurface model for permafrost thermal hydrology, the Advanced Terrestrial Simulator (ATS), is used in combination with field measurements to calibrate and identify fine scale controls of ALT in ice wedge polygon tundra in Barrow, Alaska. An iterative model refinement procedure that cycles between borehole temperature and snow cover measurements and simulations functions to evaluate and parameterize different model processes necessary to simulate freeze/thaw processes and ALT formation. After model refinement and calibration, reasonable matches between simulated and measured soil temperatures are obtained, with the largest errors occurring during early summer above ice wedges (e.g. troughs). The results suggest that properly constructed and calibrated one-dimensional thermal hydrology models have the potential to provide reasonable representation of the subsurface thermal response and can be used to infer model input parameters and process representations. The models for soil thermal conductivity and snow distribution were found to be the most sensitive process representations. However, information on lateral flow and snowpack evolution might be needed to constrain model representations of surface hydrology and snow depth.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

GEOPHIRES Simulations for Deep Direct Use (DDU) Projects

This folder contains the GEOPHIRES codes and input files for running the base case scenarios for the six deep direct-use (DDU) projects. The six DDU projects took place during 2017-2020 and were funded by the U.S. Department of Energy Geothermal Technologies Office. They investigated the potential of geothermal deep direct-use at six locations across the country. The projects were conducted by Cornell University, West Virginia University (WVU), University of Illinois (U of IL), Sandia National Laboratory (SNL), Portland State University (PSU), and National Renewable Energy Laboratory (NREL). Four projects (Cornell, WVU, U of IL, SNL) investigated geothermal for direct heating of a local campus or community, the project by PSU considered seasonal subsurface storage of solar heating, and the NREL project investigated geothermal heating for turbine inlet cooling using absorption chillers. To allow comparison of techno-economic results across the six DDU projects, GEOPHIRES simulations were set up and conducted for each project. The GEOPHIRES code was modified for each project to simulate the local application and incorporate project-specific assumptions and results such as reservoir production temperature or financing conditions. The base case input file is included which simulates the base case conditions assumed by each project team. The levelized cost of heat (LCOH) is calculated and matches the base case LCOH reported by the project teams.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Model Development for Thermal-Hydrology Simulations of a Full-Scale Heater Experiment in Opalinus Clay

Disposal of commercial spent nuclear fuel in a geologic repository is studied. In situ heater experiments in underground research laboratories provide a realistic representation of subsurface behavior under disposal conditions. Here, this study describes process model development and modeling analysis for a full-scale heater experiment in opalinus clay host rock. The results of thermal-hydrology simulation, solving coupled nonisothermal multiphase flow, and comparison with experimental data are presented. The modeling results closely match the experimental data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2018

Numerical modeling and sensitivity analysis of seawater intrusion in a dual-permeability coastal karst aquifer with conduit networks

Long-distance seawater intrusion has been widely observed through the subsurface conduit system in coastal karst aquifers as a source of groundwater contaminant. In this study, seawater intrusion in a dual-permeability karst aquifer with conduit networks is studied by the two-dimensional density-dependent flow and transport SEAWAT model. Local and global sensitivity analyses are used to evaluate the impacts of boundary conditions and hydrological characteristics on modeling seawater intrusion in a karst aquifer, including hydraulic conductivity, effective porosity, specific storage, and dispersivity of the conduit network and of the porous medium. The local sensitivity analysis evaluates the parameters' sensitivities for modeling seawater intrusion, specifically in the Woodville Karst Plain (WKP). A more comprehensive interpretation of parameter sensitivities, including the nonlinear relationship between simulations and parameters, and/or parameter interactions, is addressed in the global sensitivity analysis. The conduit parameters and boundary conditions are important to the simulations in the porous medium because of the dynamical exchanges between the two systems. The sensitivity study indicates that salinity and head simulations in the karst features, such as the conduit system and submarine springs, are critical for understanding seawater intrusion in a coastal karst aquifer. The evaluation of hydraulic conductivity sensitivity in the continuum SEAWAT model may be biased since the conduit flow velocity is not accurately calculated by Darcy's equation as a function of head difference and hydraulic conductivity. In addition, dispersivity is no longer an important parameter in an advection-dominated karst aquifer with a conduit system, compared to the sensitivity results in a porous medium aquifer. In the end, the extents of seawater intrusion are quantitatively evaluated and measured under different scenarios with the variabilities of important parameters identified from sensitivity results, including salinity at the submarine spring with rainfall recharge, sea level rise, and a longer simulation time under an extended low rainfall period.

58 GEOSCIENCES↗
DOE OSTI2021

Influence of Agricultural Managed Aquifer Recharge (AgMAR) and Stratigraphic Heterogeneities on Nitrate Reduction in the Deep Subsurface

This data package contains input files for TOUGHREACT for a modeling study examining the effects of managed aquifer recharge on agricultural lands on nitrate cycling and transport in the Central Valley of CA near Modesto. The files contain all the geochemical species, reactions, and hydrological parameters for the model. The files are text files used for the TOUGH family of code created by LBNL. To use the files a license is required. https://tough.lbl.gov/licensing-download/toughreact-licensing-download/Accompanying Paper Abstract: Agricultural managed aquifer recharge (AgMAR) is a proposed management strategy whereby surface water flows are used to intentionally flood croplands with the purpose of recharging underlying aquifers. However, legacy nitrate (NO3-) contamination in agriculturally-intensive regions poses a threat to groundwater resources under AgMAR. To address these concerns, we use a reactive transport modeling framework to better understand the effects of AgMAR management strategies (i.e., by varying the frequency, duration between flooding events, and amount of water) on N leaching to groundwater under different stratigraphic configurations and antecedent moisture conditions. In particular, we examine the potential of denitrification and nitrogen retention in deep vadose zone sediments (~15 m) using variable AgMAR application rates on two-dimensional representations of differently textured soils, soils with discontinuous bands/channels, and soils with preferential flow paths characteristic of typical agricultural field sites. Our results indicate that finer textured sediments, such as silt loams, alone or embedded within high flow regions, are important reducing zones providing conditions needed for denitrification. Simulation results further suggest that applying water all-at-once rather than in increments for a fixed volume of recharge transports higher concentrations of NO3- deeper into the profile, which has the potential to exacerbate groundwater quality. This transport into deeper depths can be aggravated by wetter antecedent soil moisture conditions. However, applying water all-at-once also increases denitrification within the vadose zone by promoting anoxic conditions. We conclude that AgMAR management strategies can be designed to enhance denitrification in the subsurface and reduce N leaching to groundwater, while specifically accounting for lithologic heterogeneity, antecedent soil moisture conditions, and depth to the water table. Our findings are potentially relevant to other systems that experience flooding inundation such as riparian corridors, floodplains, wetlands, and other managed landscapes like dedicated recharge basins.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Geochemical exports to river from the intrameander hyporheic zone under transient hydrologic conditions, Water Resources Research: Dataset

To understand how redox processes influence carbon, nitrogen, and iron cycling within the intrameander hyporheic zone, we developed a biotic and abiotic reaction network and incorporated it into the reactive transport simulator PFLOTRAN. Two dimensional reactive flow and transport simulations were performed (1) to evaluate how transient hydrological conditions control the lateral redox zonation within an intrameander region of the East River in Colorado and (2) to quantify the impact of a single meander on subsurface exports of carbon and other geochemical species to the river. This data package includes simulated hydrologic and geochemical data from 23 Oct 2015 to 22 Oct 2016.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Diagenesis is key to unlocking outcrop fracture data suitable for quantitative extrapolation to geothermal targets

Exceptionally large, well-exposed sandstone outcrops in New York provide insights into folds, deformation bands, and fractures that could influence permeability, heat exchange, and stimulation outcomes of geothermal reservoir targets. Cambrian Potsdam Sandstone with <5% porosity contains decimeter-scale open, angular-limbed monoclines <0.5 km apart with associated low-porosity mm-wide cataclastic deformation bands. Crossing and abutting relationships among sub-vertical opening-mode fractures show four chronological Sets A–D, striking NNW, NE, NW, and ENE, respectively. Fracture lengths and heights range from millimeters to tens of meters. Sets A and C macro-fractures, and possibly B and D, contain quartz deposits. All sets have abundant associated quartz cemented microfractures that also record set orientations and crosscutting relations. Quartz cement deposits—evidence of diagenesis—are the key to identifying attributes of outcrop fractures suitable for extrapolation to geothermal targets in sandstones because they show which fractures formed in the subsurface. Set A fluid inclusion homogenization temperatures (120°C–129°C) are compatible with fracture at >3 km depth. Fractures are stiff and those ≥0.05 mm (Set C) and ≥0.1 mm (Set A) are open and potentially conducive to flow. Sets A and D are abundant in outcrops with close fracture spacing—0.18 m and 0.68 m, respectively—and define a rectangular connectivity network dominated by crossing and abutting X and Y nodes. Set A aperture distributions follow a power law with slope –0.8 up to 0.15 mm; other sets have lognormal distributions. Set A and D microfractures are weakly clustered, while macro-fractures commonly have 1D anticlustered (regular or periodic) arrangements at shorter length scales (<0.2 m). Sub-horizontal fractures are barren and may have formed near the surface. Fracture heights, lengths, and spatial arrangements show good trace connectivity but low open connectivity. For geothermal applications, outcrop results predict low initial well-test permeabilities owing to quartz disconnecting open fractures, but stimulation of closely spaced microfractures and partly open macro-fractures could yield high surface area for heat exchange. Quantitative extrapolation of key fracture attributes like abundance, orientation, spatial arrangement, length, and open fracture connectivity is possible from outcrops to fractured reservoirs if differing thermal histories and diagenesis are accounted for.

02 PETROLEUM↗
DOE OSTI2025

Thermal and Well Flow Performance of Closed-Loop Geothermal in the Wattenberg Area

Closed-loop geothermal systems provide an alternative to resource-constrained hydrothermal systems and stimulation-intensive enhanced geothermal systems. In this work, we apply the slender-body theory (SBT) model, to simulate the well flow and heat transfer performance of U-loop well designs in the Wattenberg area of the Denver-Julesburg Basin. Three U-loop well patterns are investigated, including a single-, double-, and multi-lateral design. The subsurface within the area of interest is characterized by deep, hot (> 200 degrees C) igneous/metamorphic basement rock underlying multiple sedimentary formations. The lateral section(s) of the U-loop lie(s) within a target depth of 6 km, where temperatures are estimated to approach 300 degrees C. As a base case, conduction-only heat transfer is investigated through simulations with the SBT model within U-loops with open-hole laterals that exchange heat directly with the hot, dry rock using water as a working fluid. The utilization of supercritical CO2 as a heat transfer fluid is also considered. For each scenario, the system performance in terms of annual heat production and temperature profile over a 20-year project lifetime are assessed. Also, the levelized costs of heat and electricity (LCOH and LCOE) are determined using a top-down techno-economic analysis model. The results show that the performance- and cost-optimized U-loop design is one having an injection-production well spacing of 1,000 meters with ten 50-meter-spaced laterals that traverse a subsurface system with a temperature gradient of 60 degrees C/km. By injecting 20 degrees C-water at a rate of 60 kg/s through this loop, an average heat production of 19 MWth (i.e., 2.2 MWe net plant output) can be achieved, resulting in an LCOE and LCOH of $136/MWhe and $1.53/GJ, respectively, over a 20-year project life.

closed-loop geothermal↗
DOE OSTI2020

Feasibility of using reactive silicate particles with temperature-responsive coatings to enhance the security of geologic carbon storage

The large-scale deployment of geologic carbon storage will require minimizing the risk of CO 2 leakage. Here we report on a novel strategy to seal leakage pathways by injecting functionalized reactive calcium silicate (CaSiO 3 ) particles. The particles consist of a mineral silicate core within a temperature-sensitive Poly(N,N-dimethylaminoethyl methacrylate) coating that limits contact between the core and the surrounding CO 2 (aq) above a switching temperature making them relatively unreactive deep in the subsurface. If the CO 2 leaks to shallower depths and cooler temperatures, the coatings would extend and enable the silicate core to react with CO 2(aq) to form solid carbonate precipitates that reduce the permeability of the fracture or flow path. A combined laboratory- and simulation-based investigation was conducted to explore key early-stage research questions related to the feasibility of this leakage management strategy. Sand columns were injected with either coated or uncoated CaSiO 3 particles and the permeability was measured before and after exposure to CO 2 . Experiments were conducted at 35 °C and 70 °C, a range that spans the switching temperature for the coatings. Results show that at the cooler temperature the coating effectively enabled reaction and permeability reduction. At higher temperatures, the coated particles limited mineral carbonation and the permeability decrease was minimal compared to columns with uncoated particles. Here, the morphology and location of the precipitates provides insight about permeability evolution. Model simulations of 2D radial particle transport from an injection point provided insight into the effects of particle size, coating properties, and concentrations of particles in several representative rock formations.

54 ENVIRONMENTAL SCIENCES↗