Search NASA⌕ Search

SEARCH · Search NASA

Results for “Geological 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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Modeling gas migration through clay-based buffer material using coupled multiphase fluid flow and geomechanics with stress-dependent gas permeability

A model for gas migration through clay-based buffer material is developed for modeling gas generation and migration associated with deep geologic nuclear waste disposal. The model is based on a multiphase fluid flow and geomechanics simulator that is adapted to consider enhanced gas flow when gas pressure is high enough to approach the confining stress magnitude. A key feature in the model is a direct coupling between gas permeability and stress, through a non-linear stress-dependent permeability function. The model was first tested and calibrated by modelling two different laboratory gas migration tests on Wyoming (MX-80) bentonite samples. The calibrated model was then applied to model gas migration through a bentonite buffer of a large-scale gas injection test (Lasgit) conducted at the Äspö Hard Rock Laboratory in Sweden. Observed preferential gas migration along interfaces (between compacted blocks and along the canister surface) required explicit representation of such interfaces in the model. The model with the stress-dependent gas permeability accurately captured observed experimental responses in terms of gas breakthrough time, peak gas pressure, and cumulative gas flow rates. The calibrated model was finally applied to simulate migration of hydrogen gas generated within a breached nuclear waste canister over 10,000 years, involving migration of much larger gas volumes. For the considered gas generation rate and host rock properties, the generated gas could migrate through the bentonite buffer and released into the surrounding host rock at a maximum gas pressure somewhat higher than the initial total stress, though a significant amount of hydrogen remained within the buffer. This modelling sets the stage for further detailed analysis of the impact of hydrogen gas generation on the long-term performance of nuclear waste repositories.

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

Modeling the distribution of the endangered Jemez Mountains salamander (Plethodon neomexicanus) in relation to geology, topography, and climate

The Jemez Mountains salamander (Plethodon neomexicanus; hereafter JMS) is an endangered salamander restricted to the Jemez Mountains in north-central New Mexico, United States. This strictly terrestrial and lungless species requires moist surface conditions for activities such as mating and foraging. Threats to its current habitat include fire suppression and ensuing severe fires, changes in forest composition, habitat fragmentation, and climate change. Forest composition changes resulting from reduced fire frequency and increased tree density suggest that its current aboveground habitat does not mirror its historically successful habitat regime. However, because of its limited habitat area and underground behavior, we hypothesized that geology and topography might play a significant role in the current distribution of the salamander. We modeled the distribution of the JMS using a machine learning algorithm to assess how geology, topography, and climate variables influence its distribution. The best habitat suitability model indicates that geology type and maximum winter temperature (November to March) were most important in predicting the distribution of the salamander (23.5% and 50.3% permutation importance, respectively). Minimum winter temperature was also an important variable (21.4%), suggesting this also plays a role in salamander habitat. Our habitat suitability map reveals low uncertainty in model predictions, and we found slight discrepancies between the designated critical habitat and the most suitable areas for the JMS. Because geological features are important to its distribution, we recommend that geological and topographical data are considered, both during survey design and in the description of localities of JMS records once detected.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2014

Petrology and Geochemistry of Unbrecciated Harzburgitic Diogenite MIL 07001: A Window Into Vestan Geological Evolution

There is a strong case that asteroid 4 Vesta is the parent of the howardite, eucrite and diogenite (HED) meteorites. Models developed for the geological evolution of Vesta can satisfy the compositions of basaltic eucrites that dominate in the upper crust. The bulk compositional characteristics of diogenites - cumulate harzburgites and orthopyroxenites from the lower crust - do not fit into global magma ocean models that can describe the compositions of basaltic and cumulate eucrites. Recent more detailed formation models do make provision for a more complicated origin for diogenites, but this model has yet to be completely vetted. Compositional studies of bulk samples has led to the hypothesis that many diogenites were formed late by interaction of their parent melts with a eucritic crust, but those observations may alternatively be explained by subsolidus equilibration of trace elements between orthopyroxene and plagioclase and Ca-phosphate in the rocks. Differences in radiogenic Mg-26 content between diogenites and eucrites favors early formation of the former, not later formation. Understanding the origin of diogenites is crucial for understanding the petrologic evolution of Vesta. We have been doing coordinated studies of a suite of diogenites including petrologic investigations, bulk rock major and trace element studies, and in situ trace element analyses of orthopyroxene. Here we will focus on an especially unusual, and potentially key, diogenite, MIL 07001.

Mittlefehldt, D. W.↗
DOE OSTI2024

NORTH DAKOTA CARBONSAFE PHASE III: SITE CHARACTERIZATION AND PERMITTING OF GEOLOGIC STORAGE OF CARBON DIOXIDE

The Energy & Environmental Research Center (EERC), in partnership with Minnkota Power Cooperative Inc. (Minnkota), SLB, and Computer Modelling Group Ltd. (CMG), supported wide-scale deployment of carbon capture and storage (CCS) as part of the U.S. Department of Energy (DOE) National Energy Technology Laboratory Carbon Storage Assurance Facility Enterprise (CarbonSAFE) Initiative Phase III. This phase included the acquisition, analysis, and development of information to fully characterize two storage complexes to demonstrate viable storage resources for commercial volumes of CO2 (defined by DOE as a minimum of 50 million tonnes [MMt] of CO2 within a 30-year period) (National Energy Technology Laboratory, 2024). Phase III also involved the preparation, submission, and approval of North Dakota underground injection control (UIC) Class VI storage facility permits (SFPs)—required precursors to applications for Class VI injection well permits. The presumed viability of commercial-scale CCS, situated adjacent to Minnkota’s Milton R. Young Station (MRYS), is validated by Minnkota’s continued pursuit of Project Tundra—an initiative to build the world’s largest lignite-based CCS project in central North Dakota (www.projecttundrand.com). Project Tundra comprises two scopes of work, Tundra Capture (installation of postcombustion CO2 capture at MRYS) and Tundra SGS (secure geologic storage). The efforts of North Dakota CarbonSAFE Phase III, Site Characterization and Permitting, supported Tundra SGS. Extensive site-specific characterization activities included a successful multimeasurement geophysical approach and drilling a stratigraphic test well (J-ROC 1, subsequently renamed Liberty-1) adjacent to MRYS. Core collection and analyses, downhole testing and fluid sampling, and geophysical logging were performed on J-ROC 1 and on a nearby stratigraphic test well (known as J-LOC 1), which was drilled, cored, and tested under a complementary project funded by the North Dakota Lignite Research Program. The injection tests performed on J-LOC 1 positively impacted the CarbonSAFE project, resulting in fewer proposed injection wells and significant construction, operations, and monitoring cost savings. The characterization data collected and analyses performed were integrated into geologic models, and successive numerical simulations were run to determine CO2 plume extent and subsurface pressure buildup associated with the planned CO2 injection rate of nearly 4 MMt per year. The latter doubles the CarbonSAFE Initiative goal with an estimated 100 MMt of CO2 stored in 20 years. Application of the U.S. Environmental Protection Agency’s (EPA’s) method for estimating the Class VI Rule area of review (AOR) to the overpressurized Broom Creek Formation inspired an alternative method of calculation, called risk-based AOR delineation. This peer-reviewed method was applied for the first time during the storage facility-permitting process. The two SFP applications submitted in 2021 successfully resulted in North Dakota Industrial Commission (NDIC) orders in 2022 authorizing the creation of the storage facility areas and amalgamation of pore space as well as establishing financial responsibility requirements. After approval of the SFPs, Minnkota filed in 2022 applications for permits to reenter the J-ROC 1 well and to drill two new wells—all with the intended purpose to become Class VI injection wells. To establish eligibility under the Internal Revenue Code for Section 45Q tax incentives, a monitoring, reporting, and verification (MRV) plan was prepared and submitted by Minnkota to EPA in 2022, resulting in the first such plan approved in North Dakota. Also in 2022, under the National Environmental Policy Act (NEPA), Minnkota prepared and submitted an environmental information volume (EIV) describing the proposed CCS project and associated potential environmental impacts. Based on the EIV, DOE determined that the proposed construction project required an environmental assessment, and Minnkota submitted the first draft in 2023 and a revised draft in 2024. Both submissions were followed by a public comment period. Subsequently, DOE issued a finding of no significant impact (FONSI) on September 13, 2024. A successful outreach program, strongly based in the production, presentation, and dissemination of informational material, fostered an environment to aid stakeholders in making informed decisions regarding the planned project. Opportunities for public input were provided at various steps along the way, including at county planning and zoning meetings, before and during the SFP administrative hearing, and during environmental assessment public comment periods. In addition, land/pore space owners and mineral owners had various points of contact, including granting access rights, securing pore space leasing, and mineral owner notifications. Based upon the successful storage facility permitting issued by NDIC, approval of the MRV plan by EPA, and receipt of a FONSI under the NEPA, Minnkota is continuing its pursuit of Project Tundra. In December 2023, the Office of Clean Energy Demonstrations under its Carbon Capture Demonstrations Projects Program announced funding for the capture system (Office of Clean Energy Demonstrations, 2023) and a proposal for CarbonSAFE Phase IV: Construction funding was submitted in March 2024 for the storage project. A go/no-go decision to proceed with construction and operations in the Broom Creek Formation is anticipated in 2024. References National Energy Technology Laboratory, CarbonSafe Initiative, https://netl.doe.gov/carbon-management/carbon-storage/carbonsafe (accessed August 2024). Office of Clean Energy Demonstrations, 2023, OCED selects three projects in CA, ND, and TX to reduce harmful carbon pollution, create new economic opportunities, and advance carbon reducing technologies, December, www.energy.gov/oced/articles/oced-selects-three-projects-ca-nd-and-tx-reduce-harmful-carbon-pollution-create-new (accessed August 2024).

Peck, Wesley↗
NASA NTRS1981

Ascent and eruption of basaltic magma on the earth and moon

The ascent and emplacement of basaltic magma on the earth and moon is modeled by the application of geological and physical observations and constraints. Relatively simple mathematical models of the motion of gas/liquid mixtures are shown to be adequate in the treatment of basaltic eruptions, provided that allowance is made for the coalescence of gas bubbles and that realistic geological and petrochemical constraints are applied to the numerical values of variables. Because gas exsolution from magmas on the earth and moon commonly occur at depths of less than 2 km, it is generally convenient to consider separately the rise of bubble-free magmatic liquid at depth in a planetary crust and the more complex motions occurring near the surface with gas exsolution.

Wilson, L.↗
DOE OSTI2023

Implementation and verification of a user-defined element (UEL) for coupled thermal-hydraulic-mechanical-chemical (THMC) processes in saturated geological media

Efficient and accurate modeling of the coupled thermal-hydraulic-mechanical-chemical (THMC) processes in various rock formations is indispensable for designing energy geo-structures such as underground repositories for high-level nuclear wastes. Here, this work focuses on developing and verifying an implicit finite element solver for generic coupled THMC problems in geological settings. Starting from the mass, momentum, and energy balance laws, a specialized set of governing equations and a thermoporoelastic constitutive model is derived. This system is then solved by an implicit finite element (FE) scheme. Specifically, the residuals and the Jacobians are scripted in a user-defined element (UEL) subroutine which is then combined with the general-purpose FE software Abaqus Standard to solve initial-boundary value problems. Considering the complexity of the system, the UEL development follows a stepwise manner by first solving the coupled hydraulic-mechanical (HM) and thermal-hydraulic-mechanical (THM) equations before moving on to the full THMC problem. Each implementation step consists of at least one verification test by comparing computed results with closed-form analytical solutions to ensure that the various coupling effects are correctly realized. To demonstrate the robustness of the algorithm and to validate the UEL, a three-dimensional case study is performed with reference to the in-situ heating test of ATLAS at Belgium in 1980s. A hypothetical radionuclide leakage event is then simulated by activating the chemical-concentration degree of freedom and prescribing a constant high concentration at the heater's surface. The model predicts a limited contaminated regime after six years considering both diffusion and advection effects on species transport.

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

Coupled THM modeling of bentonite heating and hydration in tank tests with a new temperature-dependent water retention model

This study presents a coupled thermo-hydro-mechanical (THM) model for simulating the heating and hydration behavior of bentonite, a buffer material in deep geological repositories (DGRs). The model incorporates a new temperature-dependent soil water retention curve which captures the thermal-induced shift in water retention behavior. It also distinguishes between liquid and gas permeability, modeling intrinsic gas permeability as a function of accessible porosity to improve vapor transport and desaturation predictions. The model was validated against two large-scale tank tests, demonstrating good agreement with measured temperature, relative humidity, and water inflow data. It revealed a complex porosity evolution driven by thermal expansion, vapor movement, vapor condensation, and hydration-induced swelling during heating and hydration processes. The simulation results also suggest that the permeability of the hydration layer plays a critical role in controlling water intake. Clogging of this layer can significantly reduce the volume of water inflow during the hydration phase. Furthermore, while the model effectively captures key THM behavior, further development of the mechanical constitutive law is required to account for possible thermo-elasto-plastic volume changes and microstructural effects. Overall, the model provides a robust tool for evaluating the evolution of bentonite-based barrier material in DGRs.

Guo, Guanlong [Lawrence Berkeley National Laborato↗
NASA NTRS1982

Evaluation of thermal data for geologic applications

Sensitivity studies using thermal models indicated sources of errors in the determination of thermal inertia from HCMM data. Apparent thermal inertia, with only simple atmospheric radiance corrections to the measured surface temperature, would be sufficient for most operational requirements for surface thermal inertia. Thermal data does have additional information about the nature of surface material that is not available in visible and near infrared reflectance data. Color composites of daytime temperature, nighttime temperature, and albedo were often more useful than thermal inertia images alone for discrimination of lithologic boundaries. A modeling study, using the annual heating cycle, indicated the feasibility of looking for geologic features buried under as much as a meter of alluvial material. The spatial resolution of HCMM data is a major limiting factor in the usefulness of the data for geologic applications. Future thermal infrared satellite sensors should provide spatial resolution comparable to that of the LANDSAT data.

Kahle, A. B.↗
DOE OSTI2022

Uncertainty propagation in pore water chemical composition calculation using surrogate models

Performance assessment in deep geological nuclear waste repository systems necessitates an extended knowledge of the pore water chemical conditions prevailing in host-rock formations. In the last two decades, important progress has been made in the experimental characterization and thermodynamic modeling of pore water speciation, but the influence of experimental artifacts and uncertainties of thermodynamic input parameters are seldom evaluated. In this respect, we conducted an uncertainty propagation study in a reference geochemical model describing the pore water chemistry of the Callovian-Oxfordian clay formation. Nineteen model input parameters were perturbed, including those associated to experimental characterization (leached anions, exchanged cations, cation exchange selectivity coefficients) and those associated to generic thermodynamic databases (solubilities). A set of 13 quantities of interest were studied by the use of polynomial chaos expansions built non-intrusively with a least-squares forward stepwise regression approach. Training and validation sets of simulations were carried out using the geochemical speciation code PHREEQC. The statistical results explored the marginal distribution of each quantity of interest, their bivariate correlations as well as their global sensitivity indices. The influence of the assumed distributions for input parameters uncertainties was evaluated by considering two parametric domain sizes.

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

Pore‐Scale Modeling of Carbon Dioxide and Hydrogen Transport During Geologic Gas Storage

Abstract Geologic storage of CO 2 and H 2 are climate‐positive techniques for meeting the energy transition. While similar formations could be considered for both gases, the flow dynamics could differ due to differences in their thermophysical properties. We conduct a rigorous pore‐scale study of water/CO 2 and water/H 2 systems at relevant reservoir conditions in a Bentheimer rock sample using the lattice Boltzmann method to quantify the effects of capillary, viscous, inertial, and wetting forces during gas invasion. At similar conditions, H 2 invasion is weaker compared to CO 2 due to unfavorable viscosity ratios. Increasing flow rate, however, increases the breakthrough saturation for both gas systems in the range of capillary numbers studied. At isolated conditions of flow rate, viscosity ratio, and wettability, local inertial effects are found to be critical and show consistent increase in the invaded gas saturation. The effect of inertial forces persits for both gases across all field conditions tested.

58 GEOSCIENCES↗
NASA NTRS1997

Bounds on Lithospheric Thickness on Venus from Magellan Gravity and Topography Data

The primary objective of the work executed under NAGW-4784 is to provide constraints on the thermal and tectonic evolution of Venus. Establishing thermal and tectonic evolution models requires not only geological, but geophysical constraints, in particular the nature of temporal and spatial variations in crustal and lithospheric thickness. The major topics of study completed under NAGW-4784 (described more fully below) are: (1) detailed analyses of the resolution of Magellan Line-Of-Site (LOS) Doppler data to establish the minimum resolvable wavelength in the gravity data; (2) calculations of the global strain field in the venusian lithosphere and comparisons with global strain patterns from geological mapping; (3) study of the geological history of coronae at E. Eistla Regio; (4) estimation of crustal and lithospheric thickness by modeling of topography at asymmetric and symmetric rift-like chasmata; (5) preliminary investigations of spatial versus temporal variations in lithospheric thickness. Both the PI and Co-I have presented papers based on these topics at national and international meetings (American Geophysical Union Meetings, Lunar and Planetary Science Conferences, Chapman Conference on the Geodynamics of Venus).

Johnson, Catherine L.↗
DOE OSTI2024

Generalization of Deep-Learning Models for Classification of Local Distance Earthquakes and Explosions across Various Geologic Settings

Although accurately classifying signals from earthquakes and explosions at local distance (<250 km) remains an important task for seismic network operations, the growing volume of available seismic data presents a challenge for analysts using traditional source discrimination techniques. In recent years, deep-learning models have proven effective at discriminating between low-magnitude earthquakes and explosions measured at local distances, but it is not clear how well these models are capable of generalizing across different geological settings. To address the issue of generalization between regions, we train deep-learning models (convolutional neural networks [CNNs]) on time–frequency representations (scalograms) of three-component earthquake and explosion signals from eight different regions in the continental United States. We explore scenarios where models are trained on data from all regions, individual regions, or all but one region. We find that although CNN models trained on individual regions do not necessarily generalize well across different settings, models trained on multiple regions that include diverse path coverage generalize to new regions, with station-level accuracy of up to 90% or more for data sets from unseen regions. In general, CNN-based discrimination models significantly outperform models based on uncorrected P/S ratio (measured in the 10–18 Hz frequency band), even when CNN models are tested on data from entirely unseen regions.

58 GEOSCIENCES↗
DOE OSTI2024

Sensitivity of modelled passive margin stratigraphy to variations in sea level, sediment supply and subsidence

Abstract We produced a 10 Myr synthetic stratigraphic section using a forward stratigraphic model that generates marine deltaic stratigraphy over geological timescales. We recursively fit the model using a Bayesian inversion algorithm to test: (1) if it could be accurately reconstructed; (2) if the parameters used to create it could be recovered; and (3) the sensitivity of the model output to given model parameters and the attendant physical processes. The original synthetic stratigraphic section was produced with cyclical sea‐level variations of 40 and 30 m with 2.4 and 10 Myr periods respectively. Sediment was also supplied cyclically, in 2.4 and 10 Myr cycles with amplitudes of 30 and 80 tons/100 kyr, respectively, varying from a mean of 232 tons/100 kyr. Parameter values were sampled to fit the model using a Markov chain Monte Carlo algorithm, resulting in a ±5 m (1σ) variation between the experimental output and the original. Sea level varied by ±7 m (1σ) within the posterior distribution of parameters. As a result, both the 10 Myr and 2.4 Myr sea‐level cycles could be extracted from the original output. The variation in sediment supply was approximately ±38 tons/100 kyr (1σ) and, as a result, only the larger long‐term supply variations could be accurately recovered in refitting the model. The variation in thermal, flexural and total subsidence across those parameter sets is less than ±10 m (1σ). The original section experienced 150 m of total subsidence at the depocentre. Our results demonstrate the distinct and interpretable imprint of sea level and subsidence on continental margin stratigraphy can be quantified. Moreover, we conclude that sea‐level change produces a defined effect on the geometries of stratigraphic architecture, and that techniques applied for the purpose of delineating sea‐level variation from continental margin strata have a well‐founded conceptual basis.

Schmelz, William J.↗
DOE OSTI2016

LETS v. 3.0

Program LETS calculates the electric current distribution (in space and time) along an electrically energized steel-cased geologic borehole situated within the subsurface earth. The borehole is modeled as an electrical transmission line that “leaks” current into the surrounding geology. Parameters pertinent to the transmission line current calculation (i.e., series resistance and inductance, shunt capacitance and conductance) are obtained by sampling the electromagnetic (EM) properties of a three-dimensional (3D) geologic earth model along a (possibly deviated) well track.

Aldridge, David F.↗