Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Subsurface geologic modeling”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Interpretation of Hydraulic Fractures Based on Microseismic Response in the Marcellus Shale, Monongalia County, West Virginia, USA: Implications for Shale Gas Production

Summary Hydraulic fracturing is critical for extracting shale gas in the subsurface. The treatment technique of multistage hydraulic fracturing is widely used to maximize production. In multistage hydraulic fracturing, not all pumping stages make the same contribution for production, although the designed stimulation process is almost same in every pumping stage. In this study, we characterize the microseismic responses of hydraulic fracturing. Significant variations of microseismic characteristics are observed among the different pumping stages. Pre-existing natural fractures are examined along the horizontal well in selected stages. Combined with the Mohr’s circle analysis, results of the fracture study show that induced hydraulic fractures can be captured by pre-existing natural fractures. Induced hydraulic fractures are simulated by the unconventional fracture model (UFM), and the result reveals that the stress-shadow effect diverts the direction of hydraulic fractures. The diverted hydraulic fracture affects the development of the hydraulic fracture network, which has influence on production. The treatment of two-step pumping is investigated by application of hydraulic diffusivity. Hydraulic fracturing performance and production could be optimized by the treatment of two-step pumping in a single stage. The second pumping creates new fractures and fills the fractures with additional proppants to maintain production for a long duration.

Energy & Fuels↗
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 OSTI2025

Development of an ERT‐Based Framework for Bentonite Buffers Monitoring From Laboratory Tests: 1. Characterizing Thermal–Hydrological–Mechanical Processes

Abstract Bentonite clay is widely used in engineered barrier systems for the permanent disposal of high‐level radioactive waste due to its low permeability, high swelling capacity, and thermal stability. However, the complex thermal‐hydrological‐mechanical (THM) processes induced by heating from decaying radioactive waste and hydration from surrounding rock can lead to heterogeneous changes that are difficult to measure and predict. This study develops an Electrical Resistivity Tomography (ERT)‐based framework for monitoring THM processes, progressing from sample‐scale to bench‐scale tests, to inform field‐scale applications. Sample‐scale tests analyzed small bentonite samples under controlled variations in water content, temperature, and porosity to establish fundamental resistivity relationships. Bench‐scale tests involved larger bentonite columns subjected to heating (up to 200°C) and hydration under controlled pressure, simulating repository conditions. ERT measurements, complemented by X‐ray CT imaging, temperature monitoring, and tracing sensors, revealed coupled THM processes, such as hydration‐induced compression, swelling, and thermal gradients, leading to complex resistivity patterns. The results demonstrate the potential of ERT for capturing THM‐induced resistivity changes, though challenges remain in upscaling and quantitative analysis. This study evaluates laboratory test capabilities and proposes future improvements for understanding THM‐induced resistivity responses. A conceptual framework for ERT implementation in field‐scale monitoring is presented, synthesizing findings from both scales and exploring how ERT data can inform long‐term modeling and reduce prediction uncertainties. Overall, this ERT‐based framework offers a robust method for monitoring bentonite buffers, aiding in early issue detection and supporting the safe long‐term disposal of radioactive waste in geological repositories, while highlighting the need for future development. Plain Language Summary Bentonite clay is crucial in engineered barrier systems (EBS) for containing high‐level radioactive waste due to its ability to absorb water, swell, seal and remain stable under high temperatures. When bentonite absorbs water and heats up from radioactive decay, it experiences complex changes in its physical and mechanical properties. Understanding these changes is important for ensuring the long‐term safety and effectiveness of EBS. This study used Electrical Resistivity Tomography (ERT), a non‐invasive method that measures electrical conductivity to monitor these changes during laboratory experiments. The ERT data revealed significant variations in resistivity corresponding to changes in water content, temperature, and density, providing detailed spatial and temporal insights into the behavior of bentonite. These findings enhance our ability to predict the long‐term performance of bentonite barriers, ensuring the safe containment of radioactive waste. By improving our understanding of bentonite's behavior, this research supports the development of more reliable and effective barrier systems for radioactive waste disposal, protecting the environment and public health. Key Points ERT monitoring was employed to capture resistivity changes in bentonite during controlled heating and hydration experiments, providing insights into THM processes ERT data reveal significant resistivity changes correlated with water content, temperature, and mechanical effects, enhancing the understanding of THM dynamics in bentonite This study explores the potential of the framework for application in field‐scale EBS monitoring, emphasizing the need for integrating additional geophysical methods for comprehensive subsurface imaging

Chen, Hang↗
DOE OSTI2022

Plutonium reactive transport in fractured granite: Multi-species experiments and simulations

Plutonium (Pu) in the subsurface environment can transport in different oxidation states as an aqueous solute or as colloidal particles. The transport behavior of Pu is affected by the relative abundances of these species and can be difficult to predict when they simultaneously exist. This study investigates the concurrent transport of Pu intrinsic colloids, Pu(IV) (aq) and Pu(V-VI) (aq) through a combination of controlled experiments and semi-analytical dual-porosity transport modeling. Pu transport experiments were conducted in a fractured granite at high and low flow rates to elucidate sorption processes and their scaling behavior. In the experiments, Pu(IV) (aq) was the least mobile of the Pu species, Pu(V-VI) (aq) had intermediate mobility, and the colloidal Pu, which consisted mainly of precipitated and/or hydrolyzed Pu(IV), was the most mobile. The semi-analytical modeling revealed that the sorption of each Pu species was rate-limited, as the sorption could not be described by assuming local equilibrium in the experiments. The model was able to describe the sorption of the different Pu species that occurring either on fracture surfaces, in the pores of the rock matrix, or simultaneously in both locations. While equally good fits to the data could be achieved using any of these assumptions, a fracture-dominated process was considered to be the most plausible because it provided the most reasonable estimates of sorption rate constants. Importantly, a key result of this work is that the sorption rate constant of all Pu species tends to decrease with increasing time scales, which implies that Pu will tend to be more mobile at longer time scales than observations at shorter time scales suggest. This result has important implications for predicting the environmental impacts of Pu in the safety assessments of geologic repositories for radioactive waste disposal, and we explore potential mechanistic bases for upscaling the sorption rate constants to time and distance scales that cannot be practically evaluated in experiments.

58 GEOSCIENCES↗
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 OSTI2022

Discovering hidden geothermal signatures using non-negative matrix factorization with customized k-means clustering

Discovery of hidden geothermal resources is challenging. It requires the mining of large datasets with diverse data attributes representing subsurface hydrogeological and geothermal conditions. The commonly used play fairway analysis approach typically incorporates subject-matter expertise to analyze regional data to estimate geothermal characteristics and favorability. We demonstrate an alternative approach based on machine learning (ML) to process a geothermal dataset from southwest New Mexico (SWNM). The study region includes low- and medium-temperature hydrothermal systems. Several of these systems are not well characterized because of insufficient existing data and limited past explorative work. This study discovers hidden patterns and relations in the SWNM geothermal dataset to improve our understanding of the regional hydrothermal conditions and energy-production favorability. This understanding is obtained by applying an unsupervised ML algorithm based on non-negative matrix factorization coupled with customized k-means clustering (NMFk). NMFk can automatically identify (1) hidden signatures characterizing analyzed datasets, (2) the optimal number of these signatures, (3) the dominant data attributes associated with each signature, and (4) the spatial distribution of the extracted signatures. Here, in this study, NMFk is applied to analyze 18 geological, geophysical, hydrogeological, and geothermal attributes at 44 locations in SWNM. Using NMFk, we find data patterns and identify the spatial associations of hydrothermal signatures within two physiographic provinces (Colorado Plateau and Basin and Range) and two sub-regions of these provinces (the Mogollon-Datil volcanic field and the Rio Grande rift) in SWNM. The ML algorithm extracted five hydrothermal signatures in the SWNM datasets that differentiate between low (<90°C) and medium (90-150°C)-temperature hydrothermal systems. The algorithm also suggests that the Rio Grande rift and northern Mogollon-Datil volcanic field are the most favorable regions for future geothermal resource discovery. NMFk also identified critical attributes to identify medium-temperature hydrothermal systems in the study area. The resulting NMFk model can be applied to predict geothermal conditions and their uncertainties at new SWNM locations based on limited data from unexplored regions. The code to execute the performed analyses as well as the corresponding data can be found at https://github.com/SmartTensors/GeoThermalCloud.jl.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Detrital U-Pb zircon and 40 Ar/ 39 Ar muscovite geochronology from Middle Pennsylvanian strata in the Anadarko Basin, Texas Panhandle, USA

The Late Mississippian-Permian Anadarko Basin formed in Texas and Oklahoma, USA as the result of inversion of Neoproterozic and Cambrian rift structures. Subsidence was driven by flexural loading of the Amarillo-Wichita Uplift, and this uplift may represent the easternmost element of the Ancestral Rocky Mountains system. The northwestern part of this basin has generally been interpreted to have been filled by sediment derived from the Ancestral Front Range Uplift, ~475 km to the northwest during the early stages of basin filling. We test this model using U-Pb detrital zircon and 40 Ar/ 39 Ar detrital muscovite results from three subsurface samples of the Morrow B sandstone in the northwestern part of the Anadarko Basin. We provide a new maximum depositional age of 310.9 ± 4.9 Ma that indicates the age of the Morrow B to be no older than late Atokan to early Desmoinesian Age (Moscovian), ~10 Myr younger than previously interpreted. In contrast to some previous interpretations, we propose that the most likely source for the sediment in the Morrow B is the Amarillo Uplift to the south. Detrital zircon and detrital muscovite data have age peaks at 900–1300 Ma, 1370 Ma and 1600–1800 Ma corresponding to derivation from Grenville, Granite-Rhyolite and Yavapai-Mazatzal basement provinces, respectively. A dominant detrital zircon peak at 1370 Ma suggests that Mesoproterozoic granites in the Amarillo Uplift were exposed by Middle Pennsylvanian time, and detritus eroded from the Amarillo Uplift dominated the lowstand sediment budget of the Texas Panhandle during this time; small volumes of sediment were likely sourced from the Ancestral Front Range. In conclusion, this study presents the first detrital geochronology data from the subsurface Anadarko Basin and the first detrital muscovite data from late Paleozoic southwestern Laurentia. The results presented here highlight the interpretive power of combined detrital zircon and muscovite datasets.

58 GEOSCIENCES↗
DOE OSTI2021

Discovering Hidden Geothermal Signatures using Unsupervised Machine Learning

Discovering hidden geothermal resources is a very challenging task. It requires the mining of large datasets, including various diverse data attributes representing subsurface hydrogeological and geothermal conditions. The commonly used Play Fairway Analysis (PFA) typically relies on subject-matter expertise to analyze site or regional data to estimate geothermal conditions and prospectivity. Here, we demonstrate an alternative approach based on machine learning (ML) to process a geothermal dataset of Southwest New Mexico (SWNM). The study region includes low- and medium-temperature hydrothermal systems. However, most of these systems are poorly characterized because of insufficient existing data and limited past explorative studies. This study aims to discover hidden patterns and relationships in the SWNM geothermal dataset to better understand regional hydrothermal conditions. This is achieved by applying an unsupervised machine learning algorithm based on non-negative matrix factorization coupled with customized k-means clustering (NMFk). NMFk can automatically identify (1) hidden (latent) signatures characterizing datasets, (2) the optimal number of these signatures, (3) dominant data attributes associated with each signature, and (4) spatial distribution of the extracted signatures. Here, NMFk is applied to analyze 18 geological, geophysical, hydrogeological, geothermal attributes at 44 locations in SWNM. NMFk successfully finds data patterns and identifies the spatial associations of hydrothermal signatures with the four physiographic provinces in SWNM (Colorado Plateau, Volcanic Field, Basin and Range, and the Rio Grande rift). The algorithm identified up to 5 hydrothermal signatures in the SWNM datasets that differentiate between low- and medium-temperature hydrothermal systems in different provinces. Also, the algorithm identifies two medium-temperature hydrothermal systems in SWNM that require further exploration for geothermal resource development. Based on our analyses, 12 of the attributes are important to identify medium-temperature hydrothermal systems, and the remaining six attributes are critical to characterize low-temperature hydrothermal systems. Based on the obtained results, we identify potential physiographic provinces for further exploration to characterize them as geothermal resources. The resulting NMFk model can be applied to predict geothermal conditions and their uncertainties at new SWNM locations based on limited data from unexplored areas.

58 GEOSCIENCES↗
DOE OSTI2025

On the viability of stimulated hydrogen generation from iron-rich formations

Hydrogen-based technologies present a promising solution for the global energy transition. In addition to electrolytic production, subsurface geological formations provide a potential natural source of hydrogen. Iron-rich ultramafic rocks, in particular, are favorable for hydrogen generation through natural processes such as serpentinization. Naturally occurring reactions and migration can be enhanced through various types of stimulation, including thermal, hydraulic, and chemical treatment. Through numerical simulations, we analyzed the complex interplay of factors influencing the production and migration within the subsurface, emphasizing the importance of different stimulation techniques, catalysts, and conditions. Our findings indicate that key parameters, such as damage zone permeability and width, significantly impact producible hydrogen mass. Our results indicate that a combination of large damage zone widths, high permeability, and a stimulated reaction rate of 1 × 10 -9 can yield economically viable production rates of up to 1 kg s -1 at the wellhead. Moreover, the availability of ferrous iron, rather than the serpentinization rate itself, has been identified as the primary limiting factor in achieving economically sustainable hydrogen production. In conclusion, while an unstimulated rock volume of 0.165 km 3 yields only 45t of hydrogen in two years, various stimulation techniques can increase production to 18500t.

08 - HYDROGEN↗
DOE OSTI2023

Recent advances in integrated hydrologic models: Integration of new domains

Over the past several decades, hydrologic models have advanced from independent models of the surface and subsurface to integrated models that can capture the terrestrial hydrologic cycle within one framework. In recent years, these coupled frameworks have seen the inclusion of biogeochemical processes, ecohydrology, sedimentation and erosion, cold region hydrology, anthropogenic activities, and atmospheric processes. This expansion is the result of increased computational, data, and modeling capabilities and capacities, as well as improved understanding of the processes that drive these integrated systems. Here, in this study, we review these recent advances to integrate new processes and systems into existing terrestrial hydrologic models and highlight the significant challenges and opportunities that remain. We identify that with so many models currently available and in development, selecting the most appropriate model is difficult, and we suggest a path for new or novice modelers to find the most appropriate code based on their needs. In addition, data required to parameterize and calibrate these models can often constrain their applicability and usefulness. However, advances in environmental sensors and measurement technology, in addition to data assimilation of non-traditional data (e.g. remote sensing, qualitative data) are providing new ways of addressing this issue. As we expand hydrologic models to integrate more processes and systems, our computational demands also increase. Recent and emerging advances in computational platforms, including cloud and quantum computing, in addition to the use of machine learning to capture some processes, will continue to support the use of increasingly larger and more complex, process-based models. Finally, we highlight that it is critical to develop state-of-the-science models that are accessible to all model users, not just those applied for research and development. We encourage continued development of diverse modeling platforms, considering the user needs, data availability, and computational resources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Methanogenesis and Acetogenesis in Hydrogenotrophy with Carbonate Minerals: Dependence on Mineral Surface Area, Biofilm Growth, and Microbial Community

The production, storage, and use of hydrogen are anticipated to grow substantially to achieve energy and climate goals. Consequently, microbial communities in many terrestrial and subsurface Earth environments could be exposed to elevated hydrogen concentrations. Hydrogen stimulates metabolic processes that reduce aqueous chemical species, such as bicarbonate or sulfate, that can exchange with solid mineral phases, but the controls on microbial hydrogenotrophy with mineral sources of electron acceptors are not fully understood. Herein, we applied laboratory experiments and biogeochemical modeling to study the response of a natural microbial community to an elevated partial pressure of hydrogen in the presence of carbonate minerals of varying composition, solubility, and size. Experimental incubations and simulation results showed that hydrogen consumption by microbial communities was initially dominated by sulfate reduction and, subsequently, transitioned to acetogenesis and methanogenesis. The rates of acetogenesis and methanogenesis were not correlated with the solubility of carbonate minerals. Instead, we observed strong linear correlations between the rates and surface area of carbonate minerals. Methane and acetate production slowed down in all incubations after about 2 weeks of incubation, although biogeochemical modeling predicted that the metabolic processes were not thermodynamically limited. Electron microscopy and infrared spectroscopy showed that biofilms with diverse microorganisms grew on the carbonates during this period. The methane δ13C value significantly increased, consistent with slower growth at elevated pH. This work highlights that microbial communities form biofilm on carbonate mineral surfaces as a response to hydrogen and that biofilm formation could pose a strong kinetic limitation to hydrogenotrophic metabolism utilizing carbonate minerals.

Acetobacterium↗
DOE OSTI2023

Integration of seismic-pressure-petrophysics inversion of continuous active-seismic monitoring data for monitoring and quantifying CO 2 plume (Final Report)

The overall objective of this project is to develop and validate an integrated package of joint seismic-pressure-petrophysics inversion (jSPPI) of continuous active-source seismic monitoring dataset capable of providing real-time monitoring of CO 2 plume during geologic carbon sequestration (GCS). The three specific developments include: (a) the methodologies for fast seismic full waveform inversion of continuous active source seismic monitoring, (CASSM) datasets for simultaneously estimating velocity and attenuation, and with data assimilation; (b) joint Bayesian petrophysical inversion of seismic models and pressure data for providing and updating CO 2 saturation models; (c) the methods using multiple datasets including (Crainfield and Frio-II borehole) synthetic, laboratory, and field CASSM datasets. The outcomes of jSPPI include (a) a workflow for processing CASSM data, (b) Bayesian inversion algorithms using CASSM data and pressure response data, and (c) integration with data assimilation algorithms for continuously updating site-specific models used for prediction and reservoir management. The validation of joint FWI will be conducted using synthetic models based on the Cranfield and Frio experiments as well as field CASSM datasets collected as part of the Frio-II pilot injection. To quantify and map the mass and distribution of CO 2 (saturation), we will jointly invert velocity and attenuation measurements from the FWI with a Bayesian approach using a rock physics model for attenuation (e.g., White’s attenuation model with two selected patch sizes (White, 1976; Dutta and Seriff, 1979)). The Bayesian inversion will be applied to each time step in the CASSM survey in an updating scheme, which integrates with an ensemble of reservoir simulations at each step. A more complete experimental validation dataset will be collected as part of a mesoscale (2-3 m) gas-CO 2 injection experiment utilizing a higher frequency version of the CASSM system developed for laboratory studies; the integrated inversion will be demonstrated using this dataset which will provide both a dense geometry as well as more precise secondary confirmation measurements (e.g. saturation) typically not available in the field. The resulting real-time map of CO 2 saturation is able to provide a deeper scientific understanding of the complex, time-varying dynamics of subsurface fluid flow migration path as well as the rapid detection of CO 2 leakage hazards.

25 ENERGY STORAGE↗
NASA NTRS2010

Plasma IMS Composition Measurements for Europa and the Other Galilean Moons

NASA and ESA are planning the joint Europa Jupiter System Mission (EJSM) to the Jupiter system with specific emphasis to Europa and Ganymede, respectively. The Japanese Space Agency is also planning an orbiter mission to explore Jupiter's magnetosphere and the Galilean satellites. For NASA's Jupiter Europa Orbiter (JEO) we are developing the 3D Ion Mass Spectrometer (IMS) with two main goals which can also be applied to the other Galilean moons, 1) measure the plasma interaction between Europa and Jupiter's magnetosphere and 2) infer the 4 pi surface composition to trace elemental and significant isotopic levels. The first goal supports the magnetometer (MAG) measurements, primarily directed at detection of Europa's sub-surface ocean, while the second gives information about transfer of material between the Galilean moons, and between the moon surfaces and subsurface layers putatively including oceans. The measurement of the interactions for all the Galilean moons can be used to trace the in situ ion measurements of pickup ions back to either Europa's or Ganymede's surface from the respectively orbiting spacecraft. The IMS instrument, being developed under NASA's Astrobiology Instrument Development Program, would maximally achieve plasma measurement requirements for JEO and EJSM while moving forward our knowledge of Jupiter system composition and source processes to far higher levels than previously envisaged. The composition of the global surfaces of Europa and Ganymede can be inferred from the measurement of ejected neutrals and pick-up ions using at minimum an in situ payload including MAG and IMS also fully capable of meeting Level 1 mission requirements for ocean detection and survey. Elemental and isotopic analysis of potentially extruded oceanic materials at the moon surfaces would further support the ocean objectives. These measurements should be made from a polar orbiting spacecraft about Europa or Ganymede at height 100 km. The ejecta produced by sputtering of the surfaces of Europa and Ganymede has been shown to be representative of the surface composition. Level 2 science on surface geology and composition can then be further enhanced by addition of the following: 3D Ion Neutral Mass Spectrometer (INNS), 3D plasma electron spectrometer (ELS), and hot plasma energetic particle instrument. The measurement approach is to alternate between times measuring pickup ions and times measuring plasma and magnetic field parameters along the spacecraft trajectory. By measuring the pickup ion energy, arrival direction and mass-per-charge, the ion can be traced back along the ejection trajectory to the approximate area of origin if the 3-D electric field and magnetic field are known. In situ observations of plasma flow velocities and vector magnetic fields can be used to determine the local convective electric field (E = -VXB) along the spacecraft trajectory. By combining this information with models of the magnetospheric interaction with Europa, one can generate 3D maps of the electric and magnetic field and compute the trajectories of the pickup ions back to the surface or exospheric points of origin. In the case of Ganymede there is the additional complexity of its own internal dipole magnetic field, while Io's volcanic activity introduces the complexity of a highly structured denser atmosphere. Callisto with its less globally extended exosphere will have a simpler interaction than for Europa (i.e., more like our moon). We will discuss these differences in light of the above proposed technique. Finally, the INNS observations and neutral exosphere models are needed to estimate production rates of pickup ions. The hot plasma measurements are needed to correct for sputtering rates which can be time dependent and electron plasma observations for electron impact ionization rates. Instrument characteristics, field-of-view requirements, modes of operation and effects of radiation on instrument functionality will be discussed.

Sittler, Edward↗
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 OSTI2025

Stratigraphy‐Induced Localization of Microseismicity During CO 2 Injection in Illinois Basin

Abstract Subsurface fluid injection stimulates complex hydromechanical interaction, necessitating the integration of geomechanical data across spatial and temporal scales to consider the sophisticated behavior. Induced seismic response is usually associated with the complex reservoir architecture and pre‐existing features that are three‐dimensional, such as local stratigraphy, fractures, faults, and other discontinuities. This study encompasses laboratory characterization of the coupled hydromechanical response of cores extracted from rock formations in Illinois Basin: reservoir ‐ Mt. Simon sandstone, basal seal ‐ Argenta sandstone, and crystalline basement ‐ Precambrian rhyolite. High‐resolution numerical modeling allows considering the three‐dimensional complexity of the Illinois Basin Decatur Project with spatial resolution comparable to one of the active seismic surveys. A detailed reconstruction of the evolving state of stress in formations lacking direct stress measurements is achieved by numerical modeling that integrated laboratory‐derived hydromechanical properties, a porosity‐permeability relationship, active seismic data, and an inverted three‐dimensional porosity distribution. It appears that the microseismic clusters, mainly observed in the crystalline basement during the injection, are linked to zones experiencing more critically stressed conditions prior to injection. These zones have a potential for reactivation during the injection and are attributed to the specific local stratigraphy of the injection site, as well as transfer of triggering perturbations during the injection.

Bondarenko, N. [University of Illinois Urbana‐Cham↗
DOE OSTI2020

Preliminary Study of Potential Utilization of Triassic Rift Basins for Long-Term Carbon Sequestration

Triassic rift basins of the eastern United States present a potential option for long–term carbon dioxide (CO 2 ) sequestration. Exposed and buried basins are located near large point CO 2 sources. Because of their similar origins, comparable fill successions indicate a level of reproducibility that can be conveyed between basins. The thick rock successions of the exposed Culpeper and Gettysburg basins served as proxies for understanding of the basin sequences. Their study allowed recognition of five mappable assemblages of rock types, herein termed lithofacies associations. These lithofacies associations were formed by alluvial fan, braided and meandering streams and marginal and distal lake depositional processes.Stratigraphic architecture of the lithofacies associations within exposed basins suggests a vertical succession that consists of an initial coarse-grained fluvial succession that is progressively replaced upward and basinward by finer grained lacustrine deposits. Along the faulted margins of the basins alluvial fans transitioning to fluvial delta deposits reflect lateral progradation and contemporaneous lateral components of the basin filling.To test the model devised in exposed basins, facies associations were developed for the buried Taylorsville Basin. This basin preserves over 8,000 feet of Triassic rocks that are concealed beneath more than 2,000 feet of Cretaceous and Tertiary Coastal Plain sediments. Composition and stratigraphic architecture of these rocks are similar to those observed in exposed basins. Infilling of the basin was the result of vertical aggradation and lateral progradation of coarse-grained to fine-grained facies. Based upon study of outcrops in the Culpeper and Gettysburg basins, groups of recurring lithologic facies were identified as the fundamental constructs of the basin’s sediment infilling. These groups of facies, termed lithofacies associations, represent an amalgamation of lithologic components from broadly similar depositional systems. These lithologic associations were then extrapolated into the exposed and buried portions of the Taylorsville basins. This effort determined that Triassic rift basins provide several avenues for potential study of geologic sequestration of carbon. The characteristic rift basin succession in the basins presents possible reservoir targets within the marginal fluvial deposits. Furthermore, intrusive and extrusive mafic bodies present a potential source of fracture porosity that could serve as CO 2 reservoirs. Porosity values from thin section analysis of exposed basin samples are higher in alluvial fan and braided fluvial lithofacies. Porosity is highest in samples that show less compaction, due to a lack of ductile lithic fragments and/or higher stratigraphic position. Log data from the Taylorsville Basin indicates that porosity and permeability values in basin marginal fluvial strata are elevated. Near the center of the basin, porosity and permeability values are greatly reduced, primarily owing to the fine-grained character of the lacustrine deposits. Some Triassic basins also contain thick intervals of concordant extrusive and intrusive mafic igneous rocks. These igneous bodies are potential CO 2 reservoirs for several reasons. Firstly, extrusive lava flows provide potential storage in the layers of primary porosity that occur in the vesicules formed at the top of lava flows. Secondly, both lava flows and subsurface igneous sills exhibit extensive fracture porosity produced by the rapid cooling of the flows and intrusions. Thirdly, these igneous rocks are mafic in composition, and studies have shown that iron- and magnesium-rich mafic rocks provide sequestration opportunities through carbonate remineralization. Lastly, extrusive and intrusive igneous rocks are invariably preserved within fine-grained lake deposits. These lacustrine sediments can serve as a fine-grained confining layer that encases the igneous rocks both above and below.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Using Cosmic Ray Muons to Assess Geological Characteristics in the Subsurface

Cosmic rays are energetic nuclei and elementary particles that originate from stars and intergalactic events. The interaction of these particles with the upper atmosphere produces a wide range of secondary particles that reach the surface of the earth, of which muons are the most prominent. With enough energy, muons can travel up to a few kilometers beneath the surface of the earth before being stopped completely. The terrestrial muon flux profile and associated zenith angle can be utilized to determine geological characteristics of a location (e.g., rock overburden and density) without having to use conventional methods such as boreholes. This work uses a low-power plastic scintillator-based muon detection system as a prototype for this non-destructive geological assay methodology. Four custom designed 102 cm x 51 cm x 5 cm plastic scintillation panels are used to realize two orthogonal detection planes. Optical photons from each scintillation panel are read using OnSemi J-Series 4x4 silicon photomultiplier (SiPM) arrays in conjunction with preamplifiers. Simultaneous triggers between detectors from two planes indicate a coincidence event which is recorded using the QuarkNet data acquisition system (DAQ) from Fermi National Accelerator Laboratory. A custom detector holder was designed to securely mount the detection system and rotate the panels along the zenith to collect data at variable angles. In order to quantify the systematic uncertainties associated with the detector, such as energy depositions and angular resolution of the detector design, a Monte Carlo (MC) simulation using Geant4 is being developed. Cosmic ray flux prediction will be included in the project by adding the CORSIKA MC code to the simulation toolchain. Simulated and experimental data will drive the development and validation of a reconstruction algorithm that, upon completion, is expected to predict average overburden and rock density. Extended detector exposure to muons can be used as a means to understand changes in the surrounding environment like rock porosity. On the experimental front, muons will initially be measured at the surface, establishing the baseline flux. This is followed by recording the muon flux at variable depths and zenith angles, where the data will be used by the reconstruction algorithm to predict the overburden. The result will be benchmarked against geological surveys. The measured flux data will also be used to benchmark independent and established models. Successful proof-of-concept demonstration of this technology can open doors for long term non-invasive geological monitoring. The detector design, experimental methodology, and the benchmarking efforts are detailed in this work.

Gadey, Harish Reddy↗