Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Assessing the fidelity of shallow-water carbonates as records of the Ni isotope composition of surface seawater

Nickel is a bioessential metal that is used in enzymes important to the C, N, and O cycles, and changes in its marine abundance and bioavailability may have affected the evolutionary trajectory of early life. Changes over time in the Ni isotope composition (δ 60 Ni) of surface seawater, which reflects biological demand for Ni, could allow for the reconstruction of the dynamics of Ni demand over Earth’s history, but this approach would require geologic records of surface seawater. Here, we investigate the fidelity of shallow-water carbonates as a record of the Ni isotope composition of surface seawater by determining how Ni is first partitioned into natural carbonates and then how post-depositional processes influence the Ni signal. Our samples come from the Great Bahama Bank, which is a well-studied, modern carbonate platform often used to study ancient platforms. We found that Ni is fractionated from seawater upon incorporation into carbonates capturing shallow (<18 cm), recent deposition (0.1 ‰–0.4 ‰ lighter than seawater). Variation among these [Ni] and δ 60 Ni values may be controlled by variation in mineral proportions. Meteoric diagenesis shifts δ 60 Ni to lower values, which we attribute to isotopically light meteoric fluids. In contrast, carbonates that experienced sediment-buffered marine diagenesis with respect to Ca isotopes and Sr/Ca ratios do not appear to differ in δ 60 Ni values from sediments generally representative of their initial deposition. The sensitivity of δ 60 Ni to diagenetic reset in these samples appears comparable to the sensitivities of Ca isotopes and Sr/Ca ratios, to first order. Thus, in general, carbonates that experienced sediment-buffered marine diagenesis with respect to these elements may hold the most promise as a record of the δ 60 Ni of coeval surface seawater. Additionally, we use our results to infer that the fraction of Ni removed from seawater into carbonates is less than 10 % of the total Ni output from the global oceans and incorporation of this Ni sink into global biogeochemical models will only have a minor impact on the modeled modern Ni budget.

Great Bahama Bank carbonates↗
DOE OSTI2024

More high-impact atmospheric river-induced extreme precipitation events under warming in a high-resolution model

Extreme precipitation events, as occurred in Europe 2021, or western North America 2023, with an intensity of 75 mm/day or 100 mm/day, respectively, exert a catastrophic impact. Lower-resolution (~100 km) climate models cannot simulate the intensity, nor their response to greenhouse warming. Using a high-resolution (~25 km) model capable of simulating such events, here we show that frequency-weighted area of such events over western Europe and the west coast of North America likely expands by more than 80% under an approximately 4°C of global warming from the historical levels. Along the west coasts of Europe and North America, area impacted by atmospheric rivers-induced extreme precipitation events is projected to double, driven by intensified landfalling atmospheric rivers. Thermodynamic processes drive the increase, whereas dynamic processes reduce the intensity over western Europe but enhance it along the west coast of North America. Here our findings provide policy-relevant information for climate adaptation strategies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Stress Dependency of Brittle Creep in Granite: Insights into Source Mechanisms and Parameters

Creep in rocks refers to the gradual deformation of rock material over time under the influence of constant stress. Characterizing these deformations is of great importance for engineering design, geotechnical assessment, mining operations, geological studies, and understanding natural hazards. While laboratory experiments and a variety of numerical approaches have offered explanations for microcrack interaction and damage accumulation under the three stages of creep (primary, secondary and tertiary) in conventional creep experiments, the micromechanisms of the fractures produced in brittle creep and its dependency on the applied stress have not been explored in detail. The present study focused on investigating the fracturing mechanisms that occur during creep-induced fracturing at different stress levels and estimation of the source parameters and energy budget components. A series of uniaxial compression creep experiments have been conducted at different stress level ratios (70%, 75%, 80% and 85%), to the unconfined compressive strength (UCS) of double-flawed Barre granite specimen. Creep measurements were complemented with the Acoustic Emission (AE) measurements. The creep-induced fractures were classified into double-couple (DC), compensated linear vector dipole (CLVD) and isotropic (ISO) components using the AE moment tensor decomposition technique. The results show that non-double-couple sources dominated during creep at all the specified stress levels; however, their proportions decreased as the stress level was increased. The source parameters estimation indicated a significant increase in the magnitude of the events and the radiated seismic energy with increasing levels of stress and a slight increase in the evaluated source radius and stress drop with increasing stress levels. Furthermore, this study contributes to the existing knowledge of creep-induced fracturing by providing insights into the fracturing mechanisms and the radiated seismic energy produced, which can be helpful for the development of improved models and strategies for rock engineering and geoscience applications.

58 GEOSCIENCES↗
DOE OSTI2022

Dynamic Earth Energy Storage: Terawatt-year, Grid-scale Energy Storage Using Planet Earth as a Thermal Battery (GeoTES): Phase I Project (Final Report)

Grid-scale energy storage has been identified by the U.S. Department of Energy’s (DOE) Energy Storage Grand Challenge as a necessary technology to support the continued build-out of intermittent renewable energy resources required to attain a carbon-free energy future. To meet this goal, the 2018 Department of Energy Research and Innovation Act mandated the creation of a comprehensive program to accelerate the development and commercialization of next-generation energy storage technologies. One of numerous energy storage technology options is the storage of excess energy as heated geothermal brine in suitable geologic formations. This concept, known as reservoir thermal energy storage (RTES), geologic thermal energy storage (GeoTES), aquifer thermal energy storage (ATES), etc., relies on the storage of thermal energy in geologic formations for recovery and use in large-scale direct use geothermal (e.g., district heating, industrial processes, etc.) and electrical power generation applications. This thermal energy is derived from excess or waste heat from any high-temperature heat source, such as concentrated solar or from conventional thermal/nuclear generation. As such, RTES can potentially play a significant role in meeting the energy storage shortfall in the coming decades. RTES can provide energy arbitrage through both the storage and production of thermal energy stored in geologic formations for direct use applications and can serve as a source of hot fluids that can be used to generate electricity to support peak demand ramping, thus easing stress on transmission and distribution. This energy storage option has geographic benefits in that energy can be stored locally or regionally depending on the various needs/loads. RTES can also be located across an enormous geographic area, without the need for a traditional hydrothermal resource but where thermal gradients and hydrogeology allow economic exploitation of subsurface heat. The work conducted for this project includes (1) a review of lessons learned from past high-temperature RTES international projects; (2) geochemical experimental investigation and numerical simulations of potential domestic sedimentary reservoirs and (3) development of a thermo-hydrological-mechanical (THM) numerical simulation tool for optimizing formation properties and design parameters to maximize thermal energy storage performance.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Evidence for ground-state electron capture of K 40

Potassium-40 is a widespread, naturally occurring isotope whose radioactivity impacts estimated geological ages spanning billions of years, nuclear structure theory, and subatomic rare-event searches—including those for dark matter and neutrinoless double-beta decay. The decays of this long-lived isotope must be precisely known for its use as a geochronometer, and to account for its presence in low-background experiments. There are several known decay modes for potassium-40, but a predicted electron-capture decay directly to the ground state of argon-40 has never been observed. The existence of this decay mode impacts several fields, while theoretical predictions span an order of magnitude. Here we report on the first, successful observation of this rare decay mode, obtained by the KDK (potassium decay) Collaboration using a novel combination of a low-threshold x-ray detector surrounded by a tonne-scale, high-efficiency γ -ray tagger at Oak Ridge National Laboratory. A blinded analysis reveals a distinctly nonzero ratio of intensities of ground-state electron-captures ( I EC 0 ) over excited-state ones ( I EC * ) of I EC 0 / I EC * = 0.0095 ± stat 0.0022 ± sys 0.0010 (68% CL), with the null hypothesis rejected at 4 σ [Stukel et al. , Phys. Rev. Lett. 131 , 052503 (2023) ]. In terms of branching ratio, this unambiguous signal yields I EC 0 = 0.098 % ± stat 0.023 % ± sys 0.010 % , roughly half of the commonly used prediction. This first observation of a third-forbidden unique electron capture improves our understanding of low-energy backgrounds in dark-matter searches and has implications for nuclear-structure calculations. For example, a shell-model based theoretical estimate for the neutrinoless double-beta decay half-life of calcium-48 is increased by a factor of 7 - 2 + 3 . Our nonzero measurement shifts geochronological ages by up to a percent; implications are illustrated for Earth and solar system chronologies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

The Changing Influence of Precipitation on Soil Moisture Drought With Warming in the Mediterranean and Western North America

Abstract Anthropogenic climate change has already affected drought severity and risk across many regions, and climate models project additional increases in drought risk with future warming. Historically, droughts are typically caused by periods of below‐normal precipitation and terminated by average or above‐normal precipitation. In many regions, however, soil moisture is projected to decrease primarily through warming‐driven increases in evaporative demand, potentially affecting the ability of negative precipitation anomalies to cause drought and positive precipitation anomalies to terminate drought. Here, we use climate model simulations from Phase Six of the Coupled Model Intercomparison Project (CMIP6) to investigate how different levels of warming (1, 2, and 3°C) affect the influence of precipitation on soil moisture drought in the Mediterranean and Western North America regions. We demonstrate that the same monthly precipitation deficits (25th percentile relative to a preindustrial baseline) at a global warming level of 2°C increase the probability of both surface and rootzone soil moisture drought by 29% in the Mediterranean and 32% and 6% in Western North America compared to the preindustrial baseline. Furthermore, the probability of a dry (25th percentile relative to a preindustrial baseline) surface soil moisture month given a high (75th percentile relative to a preindustrial baseline) precipitation month is 6 (Mediterranean) and 3 (Western North America) times more likely in a 2°C world compared to the preindustrial baseline. For these regions, warming will likely increase the risk of soil moisture drought during low precipitation periods while simultaneously reducing the efficacy of high precipitation periods to terminate droughts.

Environmental Sciences & Ecology↗
DOE OSTI2023

The role of international collaboration in the United States geologic disposal research program

Abstract. More than a decade ago, the United States disposal program discontinued all research activities focused on the unsaturated fractured tuff formation at Yucca Mountain as the geologic disposal site for spent fuel and high-level radioactive waste. A new research and development (R&D) program was initiated to provide a sound technical basis for alternative disposal options across clay, crystalline, and salt rocks. The goals of this broad program were (and still are) to (1) increase confidence in the robustness of generic disposal concepts, (2) develop the science and engineering tools needed to support disposal concept implementation, and (3) conduct R&D on the direct disposal of existing dual-purpose (storage and transportation) canisters. Recognizing the benefits of international collaboration toward the common goal of safely and efficiently managing the back end of the nuclear fuel cycle, the program emphasized international cooperation as an effective strategy for sharing information and knowledge. In a multi-laboratory effort coordinated by Lawrence Berkeley National Laboratory, the United States Department of Energy (DOE) program established formal and informal cooperation partnerships with several international initiatives and institutions and developed a number of collaborative R&D activities in important research areas, such as engineered barrier integrity, near-field perturbations, radionuclide transport, performance assessment, and methods for characterization and monitoring of engineered and natural barriers. This presentation gives an overview of these R&D activities, with a specific focus on activities that improve our current understanding of the coupled thermal–hydrological–mechanical–chemical (THMC) processes occurring in engineered and natural barriers. We start with a brief review of selected international collaboration initiatives and then describe a few specific collaboration projects. We focus specifically on such studies that use experimental data sets provided by international research cooperation for joint modeling work to increase confidence in performance-relevant predictions of coupled processes. Overall, the focus on international collaboration has allowed deep engagement of US researchers with the international waste management R&D community in terms of best practices, new scientific advances, state-of-the-art simulation tools, new monitoring and performance confirmation approaches, and lessons learned. The joint R&D with international researchers, worldwide sharing of knowledge and experience, and access to relevant data and experiments from a variety of host rocks have helped our researchers to significantly improve their understanding of the current technical basis for disposal in a range of potential host rock environments. International collaboration also provides ample opportunity for training and educating junior staff that are well suited to move the United States disposal research program forward into the next decades, a promising avenue for developing a next-generation workforce of disposal scientists.

Birkholzer, Jens T.↗
DOE OSTI2020

A synthesis of methane dynamics in thermokarst lake environments

Greenhouse gas emissions from physical permafrost thaw disturbance and subsidence, including the formation and expansion of thermokarst (thaw) lakes, may double the magnitude of the permafrost carbon feedback this century. These processes are not accounted for in current global climate models. Thermokarst lakes, in particular, have been shown to be hotspots for emissions of methane (CH4), a potent greenhouse gas with 32 times more global warming potential than carbon dioxide (CO2) over a 100-year timescale. Here, we synthesize several studies examining CH4 dynamics in a representative first-generation thermokarst lake (Vault Lake, informal name) to show that CH4 production and oxidation potentials vary with depth in thawed sediments beneath the lake. This variation leads to depth-dependent differences in both in situ dissolved CO2:CH4 ratios and net CH4 production responses to additional warming. Comparing CH4 production, oxidation, and flux values from studies at Vault Lake suggests up to 99% of produced CH4 is oxidized and/or periodically entrapped before entering the atmosphere. We summarize these findings in the context of CH4 literature from thermokarst lakes and identify future research directions for incorporating thermokarst lake CH4 dynamics into estimates of the permafrost carbon feedback.

58 GEOSCIENCES↗
DOE OSTI2023

Non-Markovian Hole Excess Noise in Avalanche Amorphous Selenium Thin Films

Enhancing the signal-to-noise ratio in avalanche photodiodes by utilizing impact ionization gain requires materials exhibiting low excess noise factors. Amorphous selenium (a-Se) as a wide bandgap at ~2.1 eV, a solid-state avalanche layer, demonstrates single-carrier hole impact ionization gain and manifests ultralow thermal generation rates. A comprehensive study of the history dependent and non-Markovian nature of hot hole transport in a-Se was modeled using a Monte Carlo (MC) random walk of single hole free flights, interrupted by instantaneous phonon, disorder, hole–dipole, and impact-ionization scattering interactions. The hole excess noise factors were simulated for 0.1–15 μm a-Se thin-films as a function of mean avalanche gain. The hole excess noise factors in a-Se decreases with an increase in electric field, impact ionization gain, and device thickness. The history dependent nature of branching of holes is explained using a Gaussian avalanche threshold distance distribution and the dead space distance, which increases determinism in the stochastic impact ionization process. An ultralow non-Markovian excess noise factor of ~1 was simulated for 100 nm a-Se thin films corresponding to avalanche gains of 1000. Future detector designs can utilize the nonlocal/non-Markovian nature of the hole avalanche in a-Se, to enable a true solid-state photomultiplier with noiseless gain.

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

The damage Mechanics challenge Results: Participant predictions compared with experiment

In this article, We present results from a recent exercise where participating organizations were asked to provide model-based blind predictions of damage evolution in 3D-printed geomaterial analogue test articles. Participants were provided with a range of data characterizing both the undamaged state (e.g., ultrasonic measurements) and damage evolution (e.g., 3-point bending, unconfined compression, and Brazilian testing) of the material. In this paper, we focus on comparisons between the participants’ predictions and the previously secret challenge problem experimental observations. We present valuable lessons learned for the application of numerical methods to deformation and failure in brittle-ductile materials. The exercise also enables us to identify which specific types of calibration data were of most utility to the participants in developing their predictions. Further, we identify additional data that would have been useful for participants to improve the confidence of their predictions. Consequently, this work improves our understanding of how to better characterize a material to enable more accurate prediction of damage and failure propagation in natural and engineered brittle-ductile materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Integrated Hydro-terrestrial Modeling 2.0: Progress and Path Forward on Building a National Capability

It is the role of the U.S. federal government and its supporting agencies, including academia and future scientists, to ensure that its people have sustained and equitable freshwater services, as well as the critical knowledge necessary to make decisions about the future as it relates to freshwater services. Clear and consistent information and guidance from federal agencies is critical. Integrated Hydro-Terrestrial Modeling (IHTM), as a United States (U.S.) national capability, focuses on understanding, quantifying, and managing the replenishment of water supply through hydrologic cycle processes and their governing forces. To provide that information, we need enhanced IHTM capabilities that capitalize on the strengths of each U.S. governmental agency and its core mission. The first IHTM workshop was held in 2019, and its subsequent report was published in 2020. The U.S. Global Change Research Program (USGCRP) and member agencies held a second IHTM workshop (IHTM 2.0) from October 31 to November 2, 2023 in Reston, Virginia. The IHTM 2.0 workshop focused on the need to support a multiscale framework to accelerate research insights, better integrate operational and planning perspectives, and bridge national-to-regional capabilities to address major interdependent societal water challenges. The workshop was organized according to a “WHAT” and “HOW” framework, with the common underlying “WHY” being the integrated water resource challenges and the “WHO” defined through interagency and cooperating academic partners. The report provides a summary of plenary presentations and breakout discussions, and a road map that focuses on near-term activities.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

Global rules for translating land-use change (LUH2) to land-cover change for CMIP6 using GLM2

Anthropogenic land-use and land-cover change activities play a critical role in Earth system dynamics through significant alterations to biogeophysical and biogeochemical properties at local to global scales. To quantify the magnitude of these impacts, climate models need consistent land-cover change time series at a global scale, based on land-use information from observations or dedicated land-use change models. However, a specific land-use change cannot be unambiguously mapped to a specific land-cover change. Here, nine translation rules are evaluated based on assumptions about the way land-use change could potentially impact land cover. Utilizing the Global Land-use Model 2 (GLM2), the model underlying the latest Land-Use Harmonization dataset (LUH2), the land-cover dynamics resulting from land-use change were simulated based on multiple alternative translation rules from 850 to 2015 globally. For each rule, the resulting forest cover, carbon density and carbon emissions were compared with independent estimates from remote sensing observations, U.N. Food and Agricultural Organization reports, and other studies. The translation rule previously suggested by the authors of the HYDE 3.2 dataset, that underlies LUH2, is consistent with the results of our examinations at global, country and grid scales. This rule recommends that for CMIP6 simulations, models should (1) completely clear vegetation in land-use changes from primary and secondary land (including both forested and non-forested) to cropland, urban land and managed pasture; (2) completely clear vegetation in land-use changes from primary forest and/or secondary forest to rangeland; (3) keep vegetation inland-use changes from primary non-forest and/or secondary non-forest to rangeland. Our analysis shows that this rule is one of three (out of nine)rules that produce comparable estimates of forest cover, vegetation carbon and emissions to independent estimates and also mitigate the anomalously high carbon emissions from land-use change observed in previous studies in the 1950s. According to the three translation rules, contemporary global forest area is estimated to be 37.42×106 km2, within the range derived from remote sensing products. Likewise, the estimated carbon stock is in close agreement with reference biomass datasets, particularly over regions with more than 50 % forest cover

58 GEOSCIENCES↗
DOE OSTI2024

Climate, Hydrology, and Nutrients Control the Seasonality of Si Concentrations in Rivers

Abstract The seasonal behavior of fluvial dissolved silica (DSi) concentrations, termedDSi regime, mediates the timing of DSi delivery to downstream waters and thus governs river biogeochemical function and aquatic community condition. Previous work identified five distinct DSi regimes across rivers spanning the Northern Hemisphere, with many rivers exhibiting multiple DSi regimes over time. Several potential drivers of DSi regime behavior have been identified at small scales, including climate, land cover, and lithology, and yet the large‐scale spatiotemporal controls on DSi regimes have not been identified. We evaluate the role of environmental variables on the behavior of DSi regimes in nearly 200 rivers across the Northern Hemisphere using random forest models. Our models aim to elucidate the controls that give rise to (a) average DSi regime behavior, (b) interannual variability in DSi regime behavior (i.e., Annual DSi regime), and (c) controls on DSi regime shape (i.e., minimum and maximum DSi concentrations). Average DSi regime behavior across the period of record was classified accurately 59% of the time, whereas Annual DSi regime behavior was classified accurately 80% of the time. Climate and primary productivity variables were important in predicting Average DSi regime behavior, whereas climate and hydrologic variables were important in predicting Annual DSi regime behavior. Median nitrogen and phosphorus concentrations were important drivers of minimum and maximum DSi concentrations, indicating that these macronutrients may be important for seasonal DSi drawdown and rebound. Our findings demonstrate that fluctuations in climate, hydrology, and nutrient availability of rivers shape the temporal availability of fluvial DSi.

Environmental Sciences & Ecology↗
DOE OSTI2021

Experimental study of three-dimensional CO 2 -water drainage and fracture-matrix interactions in fractured porous media

Estimation of CO 2 storage capacity in fractured porous reservoirs requires a better understanding of the CO 2 -water flow fundamentals in fracture-matrix systems. There are few core-flood experiments on three-dimensional CO 2 -water drainage in a fracture-matrix system, yet none have examined the impacts of the capillary continuity of fractures and fracture-matrix interactions. In this study, twelve drainage experiments were conducted in four fracture-matrix columns, each comprising a vertical stack of a cylindrical rock core, a filter paper serving as an analogue to a horizontal fracture, and a ceramic plate. Here, the core sample was surrounded by open space at the top and circumferential sides to model fractures, allowing for 3-D CO 2 -water drainage in the rock core and displaced water draining across the horizontal fracture. X-ray computed tomography was conducted to visualize the dynamic invasion/drainage processes in four rock core samples showing contrasts in anisotropy, permeability, and heterogeneity. Experimental results show (1) the equilibrium CO 2 saturations in the rock matrix vary from 0.10 to 0.60 at controlled capillary pressures up to 200 kPa, (2) the CO 2 saturations in the matrix increase with water saturation in the horizontal fracture resulting in better capillary continuity for water to drain across the fracture, and (3) the fracture water saturation can be enhanced by the non-uniform fracture capillary pressure and countercurrent flow of CO 2 and water across the fracture-matrix interface. In the case of high fracture water saturation, matrix CO 2 saturation largely depends on matrix anisotropy and heterogeneity. The core-scale experimental results contribute to understand the fracture-matrix interactions and CO 2 storage efficiency in fractured porous media.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Multi-year assessment of the impact of ship-borne radiosonde observations on polar WRF forecasts in the Arctic

Abstract To compensate for the lack of conventional observations over the Arctic Ocean, ship-borne radiosonde observations have been regularly carried out during summer Arctic expeditions and the observed data have been broadcast via the global telecommunication system since 2017. With these data obtained over the data-sparse Arctic Ocean, observing system experiments were carried out using a polar-optimized version of the Weather Research and Forecasting (WRF) model and the WRF Data Assimilation (WRFDA) system to investigate their effects on analyses and forecasts over the Arctic. The results of verification against reanalysis data reveal: (1) DA effects on analyses and forecasts; (2) the reason for the year-to-year variability of DA effects; and (3) the possible role of upper-level potential vorticity in delayed DA effects. The overall assimilation effects of the extra data on the analyses and forecasts over the Arctic are positive. Initially, the DA effects are the most apparent in the temperature variables in the middle/lower troposphere, which spread to the wind variables in the upper troposphere. The effects decrease with time but reappear after approximately 120 h, even in the 240-h forecasts. The effects on forecasts vary depending on the proximity of the radiosonde observation locations to the high synoptic variability. The upper-level potential vorticity is known to play an important role in the development of Arctic cyclones, and it is suggested as a possible explanation for the delayed DA effects after about 120 h.

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

Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations

One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.

1D assumption↗
DOE OSTI2023

Low-Temperature Geothermal Geospatial Datasets: An Example from Alaska

This project is a component of a broader effort focused on geothermal heating and cooling (GHC) with the aim of illustrating the numerous benefits of incorporating GHC and geothermal heat exchange (GHX) into community energy planning and national decarbonization strategies. To better assist private sector investment, it is currently necessary to define and assess the potential of low-temperature geothermal resources. For shallow GHC/GHX fields, there is no formal compilation of subsurface characteristics shared among industry practitioners that can improve system design and operations. Alaska is specifically noted in this work, because heretofore, it has not received a similar focus in geothermal potential evaluations as the contiguous United States. The methodology consists of leveraging relevant data to generate a baseline geospatial dataset of low-temperature resources (less than 150 degrees C) to compare and analyze information accessible to anyone trying to understand the potential of GHC/GHX and small-scale low-temperature geothermal power in Alaska (e.g., energy modelers, communities, planners, and policymakers). Importantly, this project identifies data related to (1) the evaluation of GHC/GHX in the shallow subsurface, and (2) the evaluation of low-temperature geothermal resource availability. Additionally, data is being compiled to assess repurposing of oil and gas wells to contribute co-produced fluids toward the geothermal direct use and heating and cooling resource potential. In this work we identified new data from three different datasets of isolated geothermal systems in Alaska and bottom-hole temperature data from oil and gas wells that can be leveraged for evaluation of low-temperature geothermal resource potential. The goal of this project is to facilitate future deployment of GHC/GHX analysis and community-led programs and update the low-temperature geothermal resources assessment of Alaska. A better understanding of shallow potential for GHX will improve design and operations of highly efficient GHC systems. The deployment and impact that can be achieved for low-temperature geothermal resources will contribute to decarbonization goals and facilitate widespread electrification by shaving and shifting grid loads.

15 GEOTHERMAL ENERGY↗