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.

SourceResearch connectionFrom yearThrough year

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

At least 505 records · Page 28

All sources NASA NTRS DOE OSTI

DOE OSTI2021

3-D Geologic Controls of Hydrothermal Fluid Flow at Brady geothermal field, Nevada, USA

In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a throughgoing network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density.

58 GEOSCIENCES↗

NASA NTRS1999

Mars Electromagnetic Sounding Experiment: MARSES

The MARSES is the sounding instrument developed of searching for water, water-ice or permafrost layers existing in some depth under the visible surface of Mars. There are many evidences that water once was abundant on Mars. There are stream lined islands formed by flowing water, flow patterns reminiscent of wadis in Earth deserts, and outflow channels thought to have been formed by sudden outrush of subterranean water. The secondary task is to measure the soil properties of the subsurface of Mars, which include porosity, electrical resistance of the liquid phase, thermal conductivity, temperature dependence. A main task of the MARSES system is to examine changes in subsurface properties of local areas regolith on the martian surface, and to relate them to optical images and other remote sensing data in order to understand the nature of different terrain forms. The dryed up regions of Martian frozen rocks is considered to have been developing during more than 3.5 bln years, so the upper layer boundary of permafrost can serve as an indicator reflecting the course of martian paleoclimate evolution. Additional information is contained in the original extended abstract.

Ozorovich, Y. R.↗

DOE OSTI2021

3-D Geologic Controls of Hydrothermal Fluid Flow at Brady Geothermal Field, Nevada using PCA

In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a through going network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density. This submission includes the published journal article detailing this work, the published 3-D geologic map of the Brady Geothermal Area used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity, 3-D well data, along which geologic data were sampled for PCA analyses, and associated metadata file. This work was done using existing R programs.

15 GEOTHERMAL ENERGY↗

DOE OSTI2023

Using Geophysical Information to Investigate Subsurface Structure within the High-Hydraulic Conductivity Analysis Zone

Within the 200 East Area of the Central Plateau and southeastward toward the Columbia River, a high-hydraulic conductivity zone (HCZ) has been interpreted to extend through the 200-PO-1 and 200-BP-5 operable units on the Hanford Site. The HCZ is a controlling hydraulic feature that impacts groundwater flow out of the 200 East Area and the fate of eastwardly migrating plumes from the 200 West Area. The lateral extent of the HCZ is highly uncertain, and despite strong evidence for the existence of the HCZ based on water-level data and contaminant plume tracking, there is still a limited understanding of how to define its boundaries. To provide additional information on the nature and extent of the HCZ, three surface geophysical methods – electrical resistivity tomography (ERT), time-domain electromagnetics (TEM), and seismic methods – were used to collect data south of 200 East. In addition, existing data from 200 East, consisting of surface seismic data, a borehole check shot survey in 699-37-47A, and borehole stratigraphic interpretations, were used to aid interpretations of newly collected seismic data south of 200 East. This work presumed that the contrast in subsurface geophysical properties would be a first-order aid identifying a transmissive zone(s) within the HCZ analysis area by imaging contrasts and/or anomalies in geophysical properties. While seismic, ERT, and TEM methods have sensitivity to overlapping physical properties (porosity, moisture content, lithology), the resolution and physics used to acquire each of these datasets are different, and therefore the information can also be different. Figure S.1 shows the locations of the geophysical data considered in this report.

58 GEOSCIENCES↗

DOE OSTI2020

Chemo‐mechanical Alterations Induced From CO 2 Injection in Carbonate‐Cemented Sandstone: An Experimental Study at 71 °C and 29 MPa

Abstract Carbon capture, utilization, and storage may lead to mechanical degradation of the subsurface reservoir from fluid‐rock interaction, which could lead to wellbore instability or reservoir compaction. To better understand potential relationship between mechanical degradation with various carbonate cement textures and compositions in sandstone reservoirs, six flow‐through experiments were conducted. Formation water (TDS = 5,390 mg/L) enriched with CO 2 flowed through two types of Pennsylvanian Morrow B Sandstone: an ankerite‐siderite‐cemented sandstone (disseminated cement texture) and a calcite‐cemented sandstone (poikilotopic cement texture). The experiments produced little change in permeability in the ankerite‐siderite‐cemented sandstone, but permeability increased up to more than 1 order of magnitude in the calcite‐cemented sandstone. Ultrasonic measurements and cylinder‐splitting tests (also known as Brazilian tests) suggested negligible mechanical degradation of the ankerite‐siderite‐cemented sandstone. Variable changes, with significant mechanical degradation in the static moduli, were observed in the calcite‐cemented sandstone. Thus, dissolution of the disseminated ankerite‐siderite cement (0.28–0.30%) had minimal impact on modifying the flow network and the mechanical integrity of the sandstone, whereas dissolution of the poikilotopic calcite cement (0.89–1.13%, quantified with fluid chemistry and visualized with X‐ray microcomputed tomography) impacted the mechanical strength of the sandstone by disconnecting framework grains. With the high water‐to‐rock mass ratios (7.3–8.2) and number of pore volumes (147–675) employed in these experiments, potential risks are most relevant to regions near injection wells. Ultimately, the chemo‐mechanical effects induced by CO 2 injection are strongly influenced by the cement texture and composition and the burial history of the reservoir rock.

Wu, Z.↗

DOE OSTI2026

Compound flood simulations indicate rising public exposure to sewage-contaminated waters in Waikīkī, Hawai‘i

Waikīkī, Hawai‘i serves as a valuable case study for flood vulnerability due to its low elevation, present-day exposure to multiple flood sources, and role as the state’s economic engine. This study applies the high-resolution, open-source WRF-Hydro-CUFA (Coastal Urban Flood Application) model to simulate compound flood interactions among pluvial, fluvial, coastal, storm drain-driven, and subsurface processes, and to evaluate how flood characteristics in Waikīkī evolve under projected sea-level rise (SLR). Simulated water levels in an estuarine drainage canal, serving as the primary outlet for gravity-flow drainage from Waikīkī, show strong agreement with observations (r ≥ 0.935, R 2 ≥ 0.874), accurately reproducing flooding observed during recent 5-year and 50-year storm events. Simulations incorporating SLR demonstrate a shift in dominant flood drivers as sea levels rise. In addition, results identify critical thresholds at which elevated estuarine tailwater conditions cause widespread drainage backflow, eventually occurring in the absence of rainfall. These findings suggest that tidally driven flooding will progressively dominate over rainfall-driven events, increasing public exposure to highly contaminated canal and groundwater sources. The results underscore the urgent need to integrate water quality issues and drainage failure into coastal flood management strategies for low-lying urban environments, particularly those reliant on polluted estuarine waterways for stormwater management.

Compound flooding↗

DOE OSTI2024

Parts-Per-Million Carbonate Mineral Quantification with Thermogravimetric Analysis–Mass Spectrometry

Mitigating the deleterious effects of climate change requires the development and implementation of carbon capture and storage technologies. To expand the monitoring, verification, and reporting (MRV) capabilities of geologic carbon mineralization projects, we developed a thermogravimetric analysis–mass spectrometry (TGA–MS) methodology to enable quantification of <100 ppm calcite (CaCO 3 ) in complex samples. We extended TGA–MS calcite calibration curves to enable a higher measurement resolution and lower limits of quantification for evolved CO 2 from a calcite–corundum mixture. We demonstrated <100 ppm carbonate mineral quantification with TGA–MS for the first time, an outcome applicable across earth, environmental, and materials science fields. We applied this carbonate quantification method to a suite of Columbia River Basalt Group (CRBG) well cuttings recovered in 2009 from Pacific Northwest National Laboratory’s Wallula #1 Well. Our execution of this new combined calcite and calcite–corundum calibration curve TGA–MS method on our CRBG sample suite indicated average carbonate contents of 0.050 wt % in flow interiors (caprocks) and 0.400 wt % in interflow zones (reservoirs) in the upper 1250 m of the Wallula #1 Well. Finally, by advancing our knowledge of continental flood basalt-hosted carbonates in the mafic subsurface and reaching new TGA–MS quantification limits for carbonate minerals, we expand MRV capabilities and support the commercial-scale deployment of carbon mineralization projects in the Pacific Northwest United States and beyond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Hydrologic Impacts of a Strike-Slip Fault Zone: Insights from Joint 3D Body-Wave Tomography of Rock Valley

The Rock Valley fault zone (RVFZ), an intraplate strike-slip fault zone in the southern Nevada National Security Site (NNSS), hosted a series of very shallow (<3 km) earthquakes in 1993. The RVFZ may also have hydrological significance within the NNSS, potentially playing a role in regional groundwater flow, but there is a lack of local hydrological data. In the Spring of 2021, we collected active-source accelerated weight drop seismic data over part of the RVFZ to better characterize the shallow subsurface. Here, we manually picked ~17,000 P-wave travel times and over 14,000 S-wave travel times, which were inverted for P-wave velocity (V P ), S-wave velocity (V S ), and V P /V S ratio in a 3D joint tomographic inversion scheme. Seismic velocities are imaged as deep as ~700 m in areas and generally align with geologic and structural expectations. V P and V S are relatively reduced near mapped and inferred faults, with the most prominent lower V P and V S zone around the densest collection of faults. We image V P /V S ratios ranging from ~1.5 to ~2.4, the extremes of which occur at a depth of ~100 m and are juxtaposed across a fault. One possible interpretation of the imaged seismic velocities is enhanced fault damage near the densest collection of faults with relatively higher porosity and/or crack density at ~100 m depth, with patches of semiperched groundwater present in the sedimentary rock in higher V P /V S areas and drier rock in lower V P /V S areas. A relatively higher V P /V S area beneath the densest faults persists at depth, which suggests percolation of groundwater via the fault damage zone to the regionally connected lower carbonate aquifer. Potentially, the presence and movement of groundwater may have played a role in the 1993 earthquake aftershocks.

58 GEOSCIENCES↗

NASA NTRS1981

Ground ice on Mars - Inventory, distribution, and resulting landforms

A model is presented for the distribution of ground ice on Mars, and model predictions are compared with Viking Orbiter photographs of landforms possibly related to ground ice. Estimates of the amount of water originally outgassed on Mars and the amount of water presently on Mars are presented which show that approximately 90% of the estimated outgassed volume may be stored in the Martian subsurface as ground ice. The extent of the Martian cryosphere, the zone in which the temperature is below the freezing point of water and in which ground ice can exist, is examined, and it is shown that, in the presence of a protective surface layer approximately 10 m thick, ground ice may occur beneath the entire Martian surface. Observed features on the Martian surface considered most likely to reflect the presence of ground ice are discussed, including thermokarst-like pits and debris flows, and possibly polygonally patterned ground, curvilinear features and pingo-like mounds, and the geographic distributions of these features as seen in Viking photographs are examined. It is found that the possibly ice-related features are concentrated in the northern midlatitudes, the equatorial zone near Olympus Mons, and the Southern Hemisphere near the edge of the plains, indicating that subsurface ice may be present over the entire planet. The origin of the ice-related landforms is explained by a combination of volcanic heating and variations in insolation.

Rossbacher, L. A.↗

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↗

NASA NTRS1989

Spaceborne radar observations: A guide for Magellan radar-image analysis

Geologic analyses of spaceborne radar images of Earth are reviewed and summarized with respect to detecting, mapping, and interpreting impact craters, volcanic landforms, eolian and subsurface features, and tectonic landforms. Interpretations are illustrated mostly with Seasat synthetic aperture radar and shuttle-imaging-radar images. Analogies are drawn for the potential interpretation of radar images of Venus, with emphasis on the effects of variation in Magellan look angle with Venusian latitude. In each landform category, differences in feature perception and interpretive capability are related to variations in imaging geometry, spatial resolution, and wavelength of the imaging radar systems. Impact craters and other radially symmetrical features may show apparent bilateral symmetry parallel to the illumination vector at low look angles. The styles of eruption and the emplacement of major and minor volcanic constructs can be interpreted from morphological features observed in images. Radar responses that are governed by small-scale surface roughness may serve to distinguish flow types, but do not provide unambiguous information. Imaging of sand dunes is rigorously constrained by specific angular relations between the illumination vector and the orientation and angle of repose of the dune faces, but is independent of radar wavelength. With a single look angle, conditions that enable shallow subsurface imaging to occur do not provide the information necessary to determine whether the radar has recorded surface or subsurface features. The topographic linearity of many tectonic landforms is enhanced on images at regional and local scales, but the detection of structural detail is a strong function of illumination direction. Nontopographic tectonic lineaments may appear in response to contrasts in small-surface roughness or dielectric constant. The breakpoint for rough surfaces will vary by about 25 percent through the Magellan viewing geometries from low to high Venusian latitudes. Examples of anomalies and system artifacts that can affect image interpretation are described.

Ford, J. P.↗

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 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↗

NASA NTRS1985

Volcanic Debris Flows in the Elysium Region of Mars

Photogeologic studies of the Elysium volcanic province appear to provide a specific example of the importance of volcanic-ice interaction to produce the channels of Hrad and Granicus Valles. In addition, these studies shows that the channels lie on the surface of a large sedimentary deposit which is interpreted as an accumulation of volcanic debris flows or lahars. In spite of some similarities with Martian outflow channels, this latter difference may distinguish the Elysium channels from other types of Martian channels. Geologic relations are described which demonstrate that the debris flows formed amidst other volcanic activity in the Elysium region thereby suggesting that the magmatism was important to the generation of the mobilizing liquid. The lahars resulted from the melting of ground ice and liquefaction of subsurface materials. The intersection of this fluid reservoir with the regional fracture system lead to the rapid expulsion of a muddy slurry down the steep western slope of the province.

Christiansen, E. H.↗

NASA NTRS2022

On the Stability and Phase Behavior of Titan's Subsurface Liquid Columns

On Titan, liquid hydrocarbon may stay in the subsurface porous permeable crust known as the alkanofer, analogous to water in Earth’s aquifer. In addition to pressure gradient, the subsurface liquid in alkanofers is subject to vertical compositional grading due to the gravity and temperature gradient. The common wisdom is that the liquid would normally stay underground in a stability established by the pressure that increases with depth as observed in aquifers on Earth. However, Titan’s liquids consist of nitrogen and hydrocarbons, mainly methane and ethane, the behavior of which is very sensitive to temperature and pressure. Consequently, the liquid density does not always increase with depth, thus may introduce a reverse density profile that leads to vertical convective instability of the liquid column. If reverse density profiles are present, capillary pressures arising from liquid trapped within small pores in the crust can help with the column stability. The liquid held in the capillaries can seal the space below it thus helping with the stability, unless the overpressure built from underneath becomes larger than the capillary pressure, which causes leakage to allow the liquids to seep upward from the deep. This situation is analogous to hydrocarbon seeps on Earth, where oil and natural gas escape the reservoir and flow slowly through network of cracks to the surface. An algorithm based on an extended Gibbs equation commonly used in petroleum reservoir engineering is employed in this work to produce pressure, density, and compositional profiles for evaluating the stability and phase behavior of Titan’s subsurface fluids.

Sugata Tan↗

NASA NTRS2013

The Face that Launched a Thousand Slips

Helene, (approximately 17.6 kilometers mean radius) is an L4 Trojan co-orbital of Saturn's moon Dione. Its hemisphere features an unusual morphology consisting of broad depressions and a generally smooth surface patterned with streaks and grooves. The streaks appear to be oriented down-gradient, as are the grooves. This pattern suggests intensive mass-wasting as a dominant process on the leading hemisphere. Kilometer-scale impact craters are very sparse on the leading hemisphere other than the degraded kilometer-scale basins defining the overall satellite shape, and many small craters have a diffuse appearance suggesting ongoing mass wasting. Thus mass wasting must dominate surface-modifying processes at present. In fact, the mass wasting appears to have been sufficient in magnitude to narrow the divides between adjacent basins to narrow septa, similar, but in lower relief, to the honeycomb pattern of Hyperion. The prominent groves occur primarily near topographic divides and appear have cut into a broad, slightly lower albedo surface largely conforming to the present topography but elevated a few meters above the smooth surfaces undergoing mass wasting flow. Low ridges and albedo markings on the surface suggest surface flow of materials traveling up to several kilometers. Diffusive mass wasting produces smooth surfaces - such a pattern characterizes most of the low-lying surfaces. The grooves, however, imply that the transport process is advective at those locations where they occur, that is, erosion tends to concentrate along linear pathways separated by divides. In fact, in many places grooves have a fairly regular spacing of 125-160 meters, defining a characteristic erosional scale. Several questions are prompted by the unusual morphology of Helene: 1) What is the nature of the surface materials? 2) Are the transport processes gradual or catastrophic motion from one or a few events? 3) What mechanisms drive mass wasting and groove development? 4) Have the formative processes been active in the recent past? 5) Finally, is the surface accreting or eroding? The smooth character of the leading edge hemisphere of Helene and the dominance of mass wasting suggest that the surface is composed of fine-grained debris, probably dominated by dust-size to small gravel particles. The Lagrangian points of saturnian satellites are locations where planetesimals might have accreted to form co-orbital satellites such as and they may capture ejecta from their master moon. Published models suggest that Helene is a site of net accretion, but we find no extant explanation for the dominance of fine grain sizes for the surface (and probable subsurface) composition of Helene and the other Lagrangian satellites. Observation of the mass wasting tracks on Helene suggests the presence of well-defined streams of debris with low bordering levees that may be depositional features or remnants of the dissected higher surface. Some flows in grazing illumination appear to have a convex cross-section. This mass-flow morphology might be consistent with the occurrence of large-scale flow events, but which might have occurred through rapid emplacement or slow glacier-like creep. On the other hand, small craters appear to have been "softened" by creep-like processes, indicating ongoing mass wasting.

honeycomb pattern of Hyperion↗

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↗