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 973 records · Page 54

CO 2 zonal injection rate allocation and plume extent evaluation through wellbore temperature analysis

Temperature analysis during a pause in injection operations, known as warmback analysis, has been used in the petroleum industry for evaluating the injection conformance and estimating the location of the flooded front in applications, such as waterflooding oil reservoirs. Here in this work, methods are introduced to extend the application of temperature warmback analysis to estimate the zonal CO 2 injection rate and zonal CO 2 plume extent during geologic CO 2 storage in a saline aquifer. First, novel analytical solutions are developed to model transient temperature in the aquifer during the injection and subsequent shut-in periods considering two-phase flow (gaseous CO 2 and aqueous brine) conditions in the aquifer. The solution involves a discretization of the aquifer into regions; the energy and mass conservation equations for the regions are solved simultaneously considering appropriate boundary conditions at the interfaces. Two solutions techniques are presented: multi-region and three-region solutions. Inverse models are developed accordingly to evaluate the injection profile and estimate the extent of the plume front in the reservoir during the injection period. The multi-region solution results in an inversion approach that requires regression analysis. However, the three-region formulation results in a simple graphical technique for inverse modeling. The analytical solutions are validated against a thermally coupled reservoir simulation tool using different synthetic cases for CO 2 injection in deep saline aquifers. The results of the developed solutions provide a good match with numerical results during forward and inverse modeling.

54 ENVIRONMENTAL SCIENCES↗

Cleaning up Hanford: A Deep Borehole Disposal Concept for the Cs/Sr Capsules - 16249

The Hanford site in Washington State is possibly the most radioactive-contaminated site in the USA presenting a major environmental problem requiring urgent clean-up. Sealed capsules containing halides of strontium (Sr) and cesium (Cs) account for over a third of the total radioactivity. Interim storage for a period of 10-half-lives is an unappealing option since many of the capsules contain Cs-135, which has a half-life of 2.3 million years, and the dose received from the shorter-lived isotopes would remain dangerously high for several generations. Deep geological disposal of the capsules would be a safer solution. We discuss a new and revised conceptual model for disposal of the entire Hanford capsule inventory in a single deep geological borehole with a diameter that could be as little as 31 cm. The concept involves placing thirty capsules inside a steel-disposal container. Holes up to this size are regularly drilled to depths much greater than 5 km by the hydrocarbon industry and are well within current drilling capabilities. Numerical modeling work demonstrates the feasibility of the conceptual model. In particular we examine the thermal envelope generated as a function of time enabling decisions to be made on possible fillings for the borehole annulus. Variations on the disposal concept include disposing different combinations of capsules to control the heat output and thus the temperature rise. The results lend support to a comment attributed to Secretary Moniz that the Hanford capsules could be the world's first HLW to be disposed of in a deep borehole. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Ice-Covered Lakes in Gale Crater Mars: The Cold and Wet Hypothesis

Recent geological discoveries from the Mars Science Laboratory provide evidence that Gale crater may have intermittently hosted a fluvio-lacustine environment during the Hesperian, with individual lakes lasting for a period of tens to hundreds of thousands of years. (Grotzinger et al., Science, 350 (6257), 2015). Estimates of the CO2 content of the atmosphere at the time the Gale sediments formed are far less than needed by any climate model to warm early Mars (Bristow et al., Geology, submitted), given the low solar energy input available at Mars 3.5 Gya. We have therefore explored the possibility that the lakes in Gale during the Hesperian were perennially covered with ice using the Antarctic Lakes as an analog. Using our best estimate for the annual mean surface temperature at Gale at this time (approx. 230K) we computed the thickness of an ice-covered lake. These thickness range from 10-30 meters depending on the ablation rate and ice transparency and would likely inhibit sediments from entering the lake. Thus, a first conclusion is that the ice must not be too cold. Raising the mean temperature to 245K is challenging, but not quite as hard as reaching 273K. We found that a mean annual temperature of 245K ice thicknesses range from 3-10 meters. These values are comparable to the range of those for the Antarctic lakes (3-6 m), and are not implausible. And they are not so thick that sediments cannot penetrate the ice. For the ice-covered lake hypothesis to work, however, a melt water source is needed. This could come from subaqueous melting of a glacial dam in contact with the lakes (as is the case for Lake Untersee) or from seasonal melt water from nearby glaciers (as is the case for the Dry Valley lakes). More work is needed to better assess these possibilities. However, the main advantage of the ice-covered lake model (and the main reason we pursued it) is that it relaxes the requirement for a long-lived active hydrological cycle involving rainfall and runoff, which no climate model is able to produce given known constraints on the early Mars environment.

Antarctic lakes↗

The interpretation of crustal dynamics data in terms of plate interactions and active tectonics of the Anatolian Plate and surrounding regions in the Middle East

The long-term objective of this project is to interpret NASA's Crustal Dynamics measurements (SLR) in the Eastern Mediterranean region in terms of relative plate movements and intraplate deformation. The approach is to combine realistic modeling studies with analysis of available geophysical and geological observations to provide a framework for interpreting NASA's measurements. This semi-annual report concentrates on recent results regarding the tectonics of Anatolia and surrounding regions from ground based observations. Also reported on briefly is progress in the use of the Global Positioning System to densify SLR observations in the Eastern Mediterranean. Reference is made to the previous annual report for a discussion of modeling results.

Toksoz, M. Nafi↗

Decoding Zeolite Crystallization and Stage III in Nuclear Waste Glasses by Coupled Modeling and Experiments

Under specific conditions of pH and temperature, nuclear waste immobilization borosilicate glasses may exhibit a sudden acceleration in their corrosion kinetics (stage III)—a behavior that has been associated with the formation of zeolite crystals. Such accelerated dissolution may compromise the integrity of nuclear wasteforms placed in geological depositories. However, thus far, none of the available models is able to predict the thermodynamic propensity and kinetics of zeolite precipitation as a function of the solution conditions due to (i) a lack of fundamental knowledge regarding the nucleation & growth mechanisms of zeolitic phases, (ii) uncertainty regarding the compositions (types) of zeolites that may form and the rate-limiting step in their precipitation as a function of the solution conditions, and (iii) the complexities that arise due to the vast parametric space (i.e., solution chemistry, temperature, number of secondary phases, etc.) that encompass these systems under conditions of environmental exposure. To resolve these challenges, this project aimed to unambiguously identify the thermodynamic propensity for zeolite precipitation and the kinetics thereof as a function of the solution conditions (composition, pH, and temperature). To achieve this goal: 1) We identified the solution conditions and zeolite phases relevant to nuclear glass dissolution. 2) We performed a series of ab initio molecular dynamics (AIMD) simulations to compute the thermodynamic properties of a group of characteristic zeolites that features a large range of compositions, various hydration levels, a wide range of framework structures, and partial atomic site occupancies. 3) We released a first-of-a-kind self-consistent thermodynamic database that can be used to assess the kinetics and the stability fields of zeolitic phases within a Gibbs energy minimization (GEM) framework. 4) We developed a robust geochemical modeling method allowing us to predict the stability of secondary phases (including zeolites, calcium–silicate–hydrate gels, and clays) upon the dissolution of nuclear waste immobilization glasses. 5) We introduced a model that predicts the dissolution kinetics of a series of borosilicate nuclear waste immobilization glasses in terms of the topology of their atomic network. 6) We investigated the roles of the solution composition on the crystallization kinetics of phillipsite zeolites and tobermorite silicate hydrates. Via PNNL’s collaboration and engagement, this project directly supports DOE’s nuclear waste immobilization activities by offering a technical, science-based foundation that will (i) facilitate predictions of the long-term corrosion rates and extents of existing nuclear waste immobilization glasses to help ensure safe and successful vitrification operations, (ii) inform the development of advanced glass formulations with enhanced durability, and, (iii) enable cost-savings that result from making more decisive and hence less conservative predictions while offering higher levels of nuclear waste embedment in smaller, more compact glass volumes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Coupled Hydromechanical Modeling of Induced Seismicity From CO 2 Injection in the Illinois Basin

Injection of CO 2 for geologic carbon sequestration into deep sedimentary formations involves fluid pressure increases that engage hydromechanical processes that can cause seismicity by activation of existing faults. In this work, we use a coupled multiphase fluid flow and geomechanical simulator to model spatiotemporal fluid pressure and stress changes in order to study the poroelastic effect of CO 2 injection on faults in crystalline basement rock below the injection zone. The seismicity rate along features interpreted to be basement faults is modeled using Dieterich's rate-and-state earthquake nucleation model. The methodology is applied to microseismicity detected during CO 2 injection into the Mount Simon formation during the Illinois Basin—Decatur Project. The modeling accurately captures an observed reduction in seismicity rate when the injection in the second well was into a slightly shallower zone above the base of the Mount Simon formation. Moreover, the modeling shows that it is important to consider poroelastic stress changes, in addition to fluid pressure changes for accurately modeling of the observed seismicity rate.

58 GEOSCIENCES↗

Quantification of the Impact of Acidified Brine on Fracture-Matrix Transport in a Naturally Fractured Shale Using in Situ Imaging and Modeling

Understanding flow, transport, chemical reactions, and hydromechanical processes in fractured geologic materials is key for optimizing a range of subsurface processes including carbon dioxide and hydrogen storage, unconventional energy resource extraction, and geothermal energy recovery. Flow and transport processes in naturally fractured shale rocks have been challenging to characterize due to experimental complexity and the multiscale nature of quantifying continuum scale descriptions of mass exchange between micrometer-scale fractures and nanometer-scale pores. In this study, we use positron emission tomography (PET) to image the transport of a conservative tracer in a naturally fractured Wolfcamp shale core before and after the core was exposed to low pH brine conditions. Image-based experimental observations are interpreted by fitting an analytical transport model to fracture-containing voxels in the core. Results of this analysis indicate subtle increases in matrix diffusivity and a slightly more uniform fracture velocity distribution following exposure to low pH conditions. These observations are compared with a multicomponent one-dimensional reactive transport model that indicates the capacity for a 10% increase in porosity at the fracture-matrix interface as a result of the low pH brine exposure. This porosity change is the result of the dissolution of carbonate minerals in the shale matrix to low pH conditions. Further, this image-based workflow represents a new approach for quantifying spatially resolved fracture-matrix transport processes and provides a foundation for future work to better understand the role of coupled transport, reaction, and mechanical processes in naturally fractured rocks.

58 GEOSCIENCES↗

Enhancing Regional Seismic Velocity Models With Higher‐Resolution Local Results Using Sparse Dictionary Learning

Abstract We use sparse dictionary learning to develop transformations between seismic velocity models of different resolution and spatial extent. Starting with data in the common region of both models, the method can enhance a regional lower‐resolution model to match the style and resolution of local higher‐resolution results while preserving its regional coverage. The method is demonstrated by applying it to two‐dimensional V S and three‐dimensional V P and V S regional and local velocity models in southern California. The enhanced reconstructed regional results exhibit clear visual improvements, especially in the reconstructed V P / V S ratios, and better correlations with geological features. Moreover, the reconstructed regional V P , V S models outperform the original ones in comparison of simulated earthquake waveforms to observations. The improved fitting to observed waveforms extends beyond the domain of the overlapping region. The developed dictionary learning approach provides physically interpretable results and offers a powerful tool for additional applications of data enhancement in earth sciences.

Geochemistry & Geophysics↗

Hydrothermal hydration of Martian crust: illustration via geochemical model calculations

If hydrothermal Systems existed on Mars, hydration of crustal rocks may have had the potential to affect the water budget of the planet. We have conducted geochemical model calculations to investigate the relative roles of host rock composition, temperature, water-to-rock ratio, and initial fluid oxygen fugacity on the mineralogy of hydrothermal alteration assemblages, as well as the effectiveness of alteration to store water in the crust as hydrous minerals. In order to place calculations for Mars in perspective, models of hydrothermal alteration of three genetically related Icelandic volcanics (a basalt, andesite, and rhyolite) are presented, together with results for compositions based on SNC meteorite samples (Shergotty and Chassigny). Temperatures from 150 degrees C to 250 degrees C, water-to-rock ratios from 0.1 to 1000, and two initial fluid oxygen fugacities are considered in the models. Model results for water-to-rock ratios less than 10 are emphasized because they are likely to be more applicable to Mars. In accord with studies of low-grade alteration of terrestrial rocks, we find that the major controls on hydrous mineral production are host rock composition and temperature. Over the range of conditions considered, the alteration of Shergotty shows the greatest potential for storing water as hydrous minerals, and the alteration of Icelandic rhyolite has the lowest potential.

NASA Discipline Exobiology↗

Carbon Limitation Leads to Thermodynamic Regulation of Aerobic Metabolism

Organic matter (OM) metabolism in freshwater ecosystems is a critical source of uncertainty in global biogeochemical cycles, yet the processes regulating aerobic respiration in aquatic environments remain poorly understood. Aerobic respiration is typically predicted through kinetic controls such as organic carbon and oxygen concentrations, while assuming thermodynamic regulation has no influence. The role of OM thermodynamics has been mainly explored under anaerobic conditions; however, recent observations have suggested that thermodynamic properties of OM may influence aerobic respiration rates. Here, we explicitly test OM thermodynamics as a key regulator of aerobic respiration. We pair controlled microcosm experiments with ultrahigh-resolution OM characterization to show that aerobic respiration increases with OM thermodynamic favorability under carbon limitation. We also demonstrate a shift in the regulation of aerobic respiration from OM thermodynamics to nitrogen content when carbon is in excess, highlighting a central role for OM thermodynamics particularly in carbon-limited ecosystems. This illuminates a structural gap in aquatic biogeochemical models by demonstrating that both kinetic and thermodynamic constraints can modulate aerobic respiration. Our work therefore proposes a new conceptual model in which thermodynamic and nutrient limitations dually control aerobic respiration.

54 ENVIRONMENTAL SCIENCES↗

GDSA Framework Development and Process Model Integration FY2023

The Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Spent Fuel & Waste Disposition (SFWD) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and high-level nuclear waste (HLW). A high priority for SFWST disposal R&D is disposal system modeling (Sassani et al. 2021). The SFWST Geologic Disposal Safety Assessment (GDSA) work package is charged with developing a disposal system modeling and analysis capability for evaluating generic disposal system performance for nuclear waste in geologic media.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The chemical reactivity of the Martian soil and implications for future missions

Possible interpretations of the results of the Viking Biology Experiments suggest that greater than 1 ppm of a thermally labile oxidant, perhaps H2O2, and about 10 ppm of a thermally stable oxidant are present in the martian soil. We reexamine these results and discuss implications for future missions, the search for organics on Mars, and the possible health and engineering effects for human exploration. We conclude that further characterization of the reactivity of the martian regolith materials is warrented-although if our present understanding is correct the oxidant does not pose a hazard to humans. There are difficulties in explaining the reactivity of the Martian soil by oxidants. Most bulk phase compounds that are capable of oxidizing H2O to O2 per the Gas Exchange Experiment (GEx) are thermally labile or unstable against reduction by atmospheric CO2. Models invoking trapped O2 or peroxynitrates (NOO2(-)) require an unlikely geologic history for the Viking Lander 2 site. Most suggested oxidants, including H2O2, are expected to decompose rapidly under martian UV. Nonetheless, we conclude that the best model for the martian soil contains oxidants produced by heterogeneous chemical reactions with a photochemically produced atmospheric oxidant. The GEx results may be due to catalytic decomposition of an unstable oxidizing material by H2O. We show that interfacial reaction sites covering less than 1% of the available soil surfaces could explain the Viking Biology Experiments results.

Zent, Aaron P.↗

A Model for Formation of Dust, Soil and Rock Coatings on Mars: Physical and Chemical Processes on the Martian Surface

This model is one of many possible scenarios to explain the generation of the current surface material on Mars using chemical, magnetic and spectroscopic data From Mars and geologic analogs from terrestrial sites. One basic premise of this model is that the dust/soil units are not derived exclusively from local rocks, but are rather a product of global, and possibly remote, weathering processes. Another assumption in this model is that there are physical and chemical interactions of the atmospheric dust particles and that these two processes create distinctly different results on the surface. Physical processes distribute dust particles on rocks and drift units, forming physically-aggregated layers; these are reversible processes. Chemical reactions of the dust/soil particles create alteration rinds on rock surfaces and cohesive, crusted surface units between rocks, both of which are relatively permanent materials. According to this model the dominant components of the dust/soil particles are derived from alteration of volcanic ash and tephra, and contain primarily nanophase and poorly crystalline ferric oxides/oxyhydroxide phases as well as silicates. These phases are the alteration products that formed in a low moisture environment. These dust/soil particles also contain a smaller amount of material that was exposed to more water and contains crystalline ferric oxides/oxyhydroxides, sulfates and clay silicates. These components could have formed through hydrothermal alteration at steam vents or fumeroles, thermal fluids, or through evaporite deposits. Wet/dry cycling experiments are presented here on mixtures containing poorly crystalline and crystalline ferric oxides/oxyhydroxides, sulfates and silicates that range in size from nanophase to 1-2 pm diameter particles. Cemented products of these soil mixtures are formed in these experiments and variation in the surface texture was observed for samples containing smectites, non-hydrated silicates or sulfates. Reflectance spectra were measured of the initial particulate mixtures, the cemented products and ground versions of the cemented material. The spectral contrast in the visible/near-infrared and mid-infrared regions is significantly reduced for the cemented material compared to the initial soil, and somewhat reduced for the ground, cemented soil compared to the initial soil. The results of this study suggest that diurnal and seasonal cycling on Mars will have a profound effect on the texture and spectral properties of the dust/soil particles on the surface. The model developed in this study provides an explanation for the generation of cemented or crusted soil units and rock coatings on Mars and may explain albedo variations on the surface observed near large rocks or crater rims.

Bishop, Janice↗

The faint young Sun paradox: an observational test of an alternative solar model

We report the results of deep observations at radio (3.6 cm) wavelengths of the nearby solar-type star pi 01 Ursa Majoris with the Very Large Array (VLA) intended to test an alternative theory of solar luminosity evolution. The standard model predicts a solar luminosity only 75% of the present value and surface temperatures below freezing on Earth and Mars at 4 Ga, seemingly in conflict with geologic evidence for liquid water on these planets. An alternative model invokes a compensatory mass loss through a declining solar wind that results in a more consistent early luminosity. The free-free emission from an enhanced wind around nearby young Sun-like stars should be detectable at microwave frequencies. Our observations of pi 01 UMa, a 300 million year-old solar-mass star, place an upper limit on the mass loss rate of 4-5 x 10(-11) M(solar) yr-1. Total mass loss from such a star over 4 Gyr would be less than 6%. If this star is indeed an analog of the early Sun, it casts doubt on the alternative model as a solution to the faint young Sun paradox, particularly for Mars.

NASA Discipline Exobiology↗

Physically rigorous reduced-order flow models of fractured subsurface environments without explosive computational cost

Fractured media models comprise discontinuities of multiple lengths (e.g. fracture lengths and apertures, wellbore area) that fall into the relatively insignificant length scales spanning millimeter-scale fractures to centimeter-scale wellbores in comparison to the extensions of the field of interest, and challenge the conventional discretization methods imposing highly-fine meshing and formidably large numerical cost. By utilizing the recent developments in the finite element analysis of electromagnetics that allow to represent material properties on a hierarchical geometry, this project develops computational capabilities to model fluid flow, heat conduction, transport and induced polarization in large-scale geologic environments that possess geometrically-complex fractures and man-made infrastructures without explosive computational cost. The computational efficiency and robustness of this multi-physics modeling tool are demonstrated by considering various highly-realistic complex geologic environments that are common in many energy and national security related engineering problems.

42 ENGINEERING↗

Gas diffusion through variably-water-saturated zeolitic tuff: Implications for transport following a subsurface nuclear event

Noble gas transport through geologic media has important applications in the characterization of underground nuclear explosions (UNEs). Without accurate transport models, it is nearly impossible to distinguish between xenon signatures originating from civilian nuclear facilities and UNEs. Understanding xenon transport time through the earth is a key parameter for interpreting measured xenon isotopic ratios. One of the most challenging aspects of modeling gas transport time is accounting for the effect of variable water saturation of geological media. In this study, we utilize bench-scale laboratory experiments to characterize the diffusion of krypton, xenon, and sulfur hexafluoride (SF6) through intact zeolitic tuff under different saturations. Here, we demonstrate that the water in rock cores with low partial saturation dramatically affects xenon transport time compared to that of krypton and SF6 by blocking sites in zeolitic tuff that preferentially adsorb xenon. This leads to breakthrough trends that are strongly influenced by the degree of the rock saturation. Xenon is especially susceptible to this phenomenon, a finding that is crucial to incorporate in subsurface gas transport models used for nuclear event identification. We also find that the breakthrough of SF6 diverges significantly from that of noble gases within our system. When developing field scale models, it is important to understand how the behavior of xenon deviates from chemical tracers used in the field, such as SF6 (Carrigan et al., 1996). These new insights demonstrate the critical need to consider the interplay between rock saturation and fission product sorption during transport modeling, and the importance of evaluating specific interactions between geomedia and gases of interest, which may differ from geomedia interactions with chemical tracers.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Unconfined Compression Results on Core from Boreholes U-15n#12 and U-15n#13, NNSS in support of SPE

The Nevada National Security Site (NNSS) serves as the geologic setting for a Source Physics Experiment (SPE) program. The SPE provides ground truth data to create and improve strong ground motion and seismic S-wave generation and propagation models. The NNSS was chosen as the test bed because it provides a variety of geologic settings ranging from relatively simple to very complex. Each series of SPE testing will comprise the setting and firing of explosive charges (source) placed in a central borehole at varying depths and recording ground motions in instrumented boreholes located in two rings around the source, positioned at different radii. Modeling using advanced simulation codes will be performed both before and after each test to predict ground response and to improve models based on acquired field data, respectively. A key component in the predictive capability and ultimate validation of the models is the full understanding of the intervening geology between the source and the instrumented boreholes including the geomechanical behavior of the site's rock/structural features. This report summarizes unconfined compression testing (UCS) from coreholes U-15n#12 and U-15n#13 and compares those datasets to UCS results from coreholes U-15n and U-15n#10. U-15n#12 corehole was drilled at -60° to the horizontal and U-15n#13 was drilled vertically in granitic rock (quartz monzonite) after the third SPE shot. Figure 1 illustrates at the surface, U 15n#12 and U-15n#13 coreholes were approximately 30 meters and 10 meters from the central SPE borehole (U-15n) respectively. Corehole U-15n#12 intersects the central SPE borehole (U 15n) at a core depth of 174 feet (approximately 150 feet vertical depth). The location of U 15n#12 and U-15n#13 is the site of the first, second and third SPE's, in Area 15 of the NNSS.

58 GEOSCIENCES↗