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.

SourceResearch connectionFrom yearThrough year

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

At least 739 records · Page 41

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Continental Scale Hydrostratigraphy: Comparing Geologically Informed Data Products to Analytical Solutions

Abstract This study synthesizes two different methods for estimating hydraulic conductivity (K) at large scales. We derive analytical approaches that estimate K and apply them to the contiguous United States. We then compare these analytical approaches to three‐dimensional, national gridded K data products and three transmissivity (T) data products developed from publicly available sources. We evaluate these data products using multiple approaches: comparing their statistics qualitatively and quantitatively and with hydrologic model simulations. Some of these datasets were used as inputs for an integrated hydrologic model of the Upper Colorado River Basin and the comparison of the results with observations was used to further evaluate the K data products. Simulated average daily streamflow was compared to daily flow data from 10 USGS stream gages in the domain, and annually averaged simulated groundwater depths are compared to observations from nearly 2000 monitoring wells. We find streamflow predictions from analytically informed simulations to be similar in relative bias and Spearman's rho to the geologically informed simulations. R ‐squared values for groundwater depth predictions are close between the best performing analytically and geologically informed simulations at 0.68 and 0.70 respectively, with RMSE values under 10 m. We also show that the analytical approach derived by this study produces estimates of K that are similar in spatial distribution, standard deviation, mean value, and modeling performance to geologically‐informed estimates. The results of this work are used to inform a follow‐on study that tests additional data‐driven approaches in multiple basins within the contiguous United States.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

The age and potential causes of the giant Green Lake Landslide, Fiordland, New Zealand

Abstract Landslide deposits preserved in the geological record afford opportunities to better inform hillslope and seismic hazard and risk models, particularly in regions where observational records are short. In the Southern Alps of New Zealand, small coseismic landslides are frequent, but the geological record preserves several instances of more substantial (> 1 km 3 ) but infrequent mass failures. With an estimated volume of 27 km 3 , the giant Green Lake Landslide represents one of the largest subaerial landslides on Earth. Previous work has suggested this deep-seated mass movement was most likely triggered by high-intensity seismic shaking, but that local structural weakness and/or glacial debuttressing may help to explain the anomalously large failure volume. Resolving the potential contribution of the latter is important given predictions concerning alpine deglaciation in the coming decades to centuries. Key to resolution are secure chronologies of landslide emplacement and past glacier change. Here we present in situ cosmogenic 10 Be exposure ages from the Green Lake Landslide that suggest an emplacement age of 15.5 ± 0.7 ka. Recent work shows that glacial retreat in the region was underway by 19 ka, indicating that the Green Lake Landslide was emplaced 3–4 kyr after the onset of glacier retreat. Given the geometry of the former confining valley glacier, we expect that the deglaciation-landslide age gap is closer to the upper end of this estimate. If correct, this conclusion places greater weight on the roles of local geological structure and/or a great earthquake as factors contributing to the exceptionally large volume of this event.

58 GEOSCIENCES↗

DOE OSTI2020

2D Dam-Break Analysis of L Lake and PAR Pond Dams Using HEC-RAS

In 1991 a dam-break study was conducted for the high hazard dams located at L Lake and PAR Pond on the Savannah River Site. Two scenarios were considered, over topping from a Probable Maximum Flood (PMF), and a fair weather dam-break for either or both dams. Unfortunately, no inundation map was developed from the study. The purpose of this project was to redo the original dam-break study with improved data and methodology to generate Inundation maps to assist with emergency response and evacuation plans. The Hydrologic Engineering Center's River Analysis System (HEC-RAS) is a free to download river analysis modeling program developed by the US Army Corps of Engineers capable of 1D and 2D hydraulic calculations. Version 5.0.7 (released March 2019) was used for this project. Digital elevation models for the area were retrieved from the US Geological Survey database and converted to a .hdf file within the program. Then the 2D flood area was identified from the contours. Both L Lake and PAR Pond were inputted as 1D storage areas because DEM data does not contain elevation values under water bodies. An elevation vs volume curve was available for both storage areas. Initial elevations were set for both scenarios. Both dams are earthen dams. The Steel Creek dam at L Lake has 6 ft diameter conduit with an upper and lower sluice gate. PAR pond dam consist of a weir and sluice gate connected to an 8 x 8 ft channel. Both outlets were modeled with a pool elevation vs discharge curve. The steel creek dam sluice gates were assumed to be fully open in all cases. As in the previous study, the dams were set to breach when they were overtopped by 1.5 ft during PMF conditions (Figure 1). A fair weather breach was set to be due to a piping failure (Figure 1). In the dual dam break during fair weather conditions the PAR pond dam fails 3 hours after Steel Creek to achieve maximum flooding in the down stream reaches. Simulation was run 3 days for each case and with a 1 minute computational interval. Maximum flooding occurs under PMF conditions with the failure of both dams. PAR pond dam fails first 16 hours and 32 minutes after the start of the simulation with the Steel creek dam failing 6 minutes later. In all cases, the bridges and roads spanning Steel Creek and Lower Three Runs will be inundated and potentially washed away. The Burtons Ferry Highway south of the storage areas will be partially flooded during PMF failure, dual fair weather failure, and PAR pond failure under fair weather conditions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Comparing modelling approaches for a generic nuclear waste repository in salt

This paper contains a comparison of five modelling approaches for a simplified nuclear waste repository in a domal salt formation. It is the result of a four-year collaboration between five international teams on Task F of the DECOVALEX-2023 project on performance assessment modelling. The primary objectives of Task F are to build confidence in the models, methods, and software used for performance assessment (PA) of deep geologic nuclear waste repositories, and/or to bring to the fore additional research and development needed to improve PA methodologies. This work demonstrates how these objectives are accomplished through staged development and comparison of the models and methods used by participating teams in their PA frameworks. Participating teams made a wide range of model assumptions, ranging from compartmentalized networks to full 3D models of the salt formation and repository. Despite differences in the modelling strategies, all models indicate that salt compaction and diffusion of radionuclides in brine are key processes in the repository. For the isothermal spent nuclear fuel and vitrified waste scenario with multiple early failures considered, all models indicate little of the disposed radionuclides will migrate beyond the repository seal over the 100,000-year simulations. In general, the model output quantities have the largest differences over the short term and near the waste. Disparities between the models are believed to be due to differing simplifications from the conceptual model.

DECOVALEX↗

DOE OSTI2023

DECOVALEX-2023: Task F2-Salt Final Report

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

DECOVALEX-2023: Task F2 Salt Final Report

The subject of Task F of DECOVALEX-2023 concerns performance assessment modelling of radioactive waste disposal in deep mined repositories. The primary objectives of Task F are to build confidence in the models, methods, and software used for performance assessment (PA) of deep geologic nuclear waste repositories, and/or to bring to the fore additional research and development needed to improve PA methodologies. In Task F2- (salt), these objectives have been accomplished through staged development and comparison of the models and methods used by participating teams in their PA frameworks. Coupled-process submodels and deterministic simulations of the entire PA model for a reference scenario for waste disposal in domal salt have been conducted. The task specification has been updated continuously since the initiation of the project to reflect the staged development of the conceptual repository model and performance metrics. Thermal, hydrological, mechanical, and chemical properties of individual components of the engineered and natural system were chosen for relevance by participating teams. The salt reference case system was characterized using data and measurements collected at relevant underground research laboratories (URLs), field sites, and simulation results from teams with specialized modelling capability. Participating teams made a wide range of model assumptions from compartmentalized networks to full 3D models of the salt formation. No single contributed model includes full-fidelity representation of all the features, events, and processes (FEPs) detailed in the task specification, but almost all features and processes are represented in at least one model. Despite differences in the modelling strategies developed by participating teams, all models indicate that salt compaction and radionuclide diffusion are key processes in the repository, and for the FEPs and model scenario considered, little of the disposed radionuclides will migrate beyond the repository seal over the 100,000 year simulations. In general, the model output quantities have the largest differences over the short term and near the waste. The models tend to be more similar further from waste and at later time. Disparities between the models are believed to be due to differing simplifications from the task specification, some of which are chosen simplifications to reduce complexity, and some are restrictions imposed by the modelling tools. A second round of this task has been accepted for DECOVALEX-2027 in conjunction with Task F1 on crystalline PA modelling. The future round includes waste package heating, improved modelling of salt creep closure, additional comparisons of coupled-process sub-models, and the impact of repository engineering design on radionuclide migration in the repository. Participants will also propose and finalize a set of uncertain inputs for the reference case simulations, propagate these uncertainties in a set of realizations, and conduct sensitivity analyses on the simulation results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS2016

The Pliocene Model Intercomparison Project (PlioMIP) Phase 2: Scientific Objectives and Experimental Design

The Pliocene Model Intercomparison Project (PlioMIP) is a co-ordinated international climate modelling initiative to study and understand climate and environments of the Late Pliocene, as well as their potential relevance in the context of future climate change. PlioMIP examines the consistency of model predictions in simulating Pliocene climate and their ability to reproduce climate signals preserved by geological climate archives. Here we provide a description of the aim and objectives of the next phase of the model intercomparison project (PlioMIP Phase 2), and we present the experimental design and boundary conditions that will be utilized for climate model experiments in Phase 2. Following on from PlioMIP Phase 1, Phase 2 will continue to be a mechanism for sampling structural uncertainty within climate models. However, Phase 1 demonstrated the requirement to better understand boundary condition uncertainties as well as uncertainty in the methodologies used for data-model comparison. Therefore, our strategy for Phase 2 is to utilize state-of-the-art boundary conditions that have emerged over the last 5 years. These include a new palaeogeographic reconstruction, detailing ocean bathymetry and land-ice surface topography. The ice surface topography is built upon the lessons learned from offline ice sheet modelling studies. Land surface cover has been enhanced by recent additions of Pliocene soils and lakes. Atmospheric reconstructions of palaeo-CO2 are emerging on orbital timescales, and these are also incorporated into PlioMIP Phase 2. New records of surface and sea surface temperature change are being produced that will be more temporally consistent with the boundary conditions and forcings used within models. Finally we have designed a suite of prioritized experiments that tackle issues surrounding the basic understanding of the Pliocene and its relevance in the context of future climate change in a discrete way.

climate models↗

DOE OSTI2022

Using Cosmic Ray Muons to Assess Geological Characteristics in the Subsurface

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

Gadey, Harish Reddy↗

DOE OSTI2024

Modelling of failure and fracture development of the Callovo-Oxfordian claystone during an in-situ heating experiment associated with geological disposal of high-level radioactive waste

Not Available

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1990

Tectonic evolution of Eastern Ishtar Terra, Venus

Tectonic deformation is cited as the cause for the complex morphology and topography of Eastern Ishtar Terra. A number of models proposed to account for this deformation are reviewed, and their predictions for the geology of Fortuna Tessera are examined. The detailed geology, topography, and structural fabric of Eastern Ishtar are reviewed, and these observations are compared to the predictions of the existing models for the formation and evolution of Fortuna Tessera. On the basis of the analysis, a number of conclusions are drawn regarding the evolution of Eastern Ishtar Terra. These include the recognition of primary morphological and topographical patterns which allow Eastern Ishtar Terra to be divided into four major physiographic regions named Maxwell Montes, Eastern Fortuna Tessera, Central Fortuna Tessera, and Northern Fortuna Tessera. A model is suggested that accounts for the morphology and high topography of Eastern Ishtar Terra.

Vorder Bruegge, R. W.↗

NASA NTRS2011

The Dawn Topography Investigation

The objective of the Dawn topography investigation is to derive the detailed shapes of 4 Vesta and 1 Ceres in order to create orthorectified image mosaics for geologic interpretation, as well as to study the asteroids' landforms, interior structure, and the processes that have modified their surfaces over geologic time. In this paper we describe our approaches for producing shape models, plans for acquiring the needed image data for Vesta, and the results of a numerical simulation of the Vesta mapping campaign that quantify the expected accuracy of our results. Multi-angle images obtained by Dawn's framing camera will be used to create topographic models with 100 m/pixel horizontal resolution and 10 m height accuracy at Vesta, and 200 m/pixel horizontal resolution and 20 m height accuracy at Ceres. Two different techniques, stereophotogrammetry and stereophotoclinometry, are employed to model the shape; these models will be merged with the asteroidal gravity fields obtained by Dawn to produce geodetically controlled topographic models for each body. The resulting digital topography models, together with the gravity data, will reveal the tectonic, volcanic and impact history of Vesta, and enable co-registration of data sets to determine Vesta's geologic history. At Ceres, the topography will likely reveal much about processes of surface modification as well as the internal structure and evolution of this dwarf planet.

dawn↗

NASA NTRS1987

Vector constraints on western U.S. deformation from space geodesy, neotectonics, and plate motions

The rate-of-slip vector on the San Andreas fault in central California estimated from geodetic and Holocene geological data is inconsistent with the prediction of rigid plate models such as RM2. This well-known 'San Andreas discrepancy' is diagnostic of plate deformation distributed both east of the fault in the Basin and Range and west of the fault along the California continental margin. Constraints on the integrated deformation rates across these two regions consistent with: (1) the kinematical boundary conditions imposed by the rigid plate model; (2) neotectonic and paleoseismic estimates of deformation rates; (3) ground-based geodetic measurements; and (4) rates of change observed by very long baseline interferometry along seven baselines to western U.S. sites are constructed. The space-geodetic data on Basin and Range extension taken over a 4-year interval are compatible with geological observations averaged over the Holocene; the best estimate of its integrated deformation rate, provided by the joint inversion of both data types, is 9.7 + or - 2.1 mm/yr, N 56 deg W + or - 10 deg, too small and in the wrong direction to account entirely for the San Andreas discrepancy. The integral of this deformation, estimated by subtracting the Basin and Range contribution from the discrepancy vector, requires significant right-lateral shear parallel to the San Andreas (13 + or - 5 mm/yr) and some compression perpendicular to it (9 + or - 3 mm/yr).

Minster, J. Bernard↗

DOE OSTI2022

HFTS-1 Natural Joint Data and Engineering Summary

This report provides an analysis of engineering and geologic data collected from the Hydraulic Fracturing Test Site #1 (HTSF-1) project in the southern Midland Basin, Reagan County, Texas. The site is being studied as part of the Science-informed Machine Learning for Accelerating Real-Time Decisions in Subsurface Applications (SMART) Initiative at the National Energy Technology Laboratory (NETL). The data collected is intended to provide a basis of understanding of the site and to construct simulation models. This report provides a summary analysis of the engineering and geologic data collected from the project to provide a basis of understanding of the site and to construct simulation models. The report examines various possible correlations in engineering properties based on natural fracture data from the program, which consisted of 11 horizontal wells, one vertical well, and one slant well. The horizontal wells are in two horizons: 1) Upper and 2) Middle Wolfcamp formations. Program data include: fracture frequency and fracture orientation data from four core runs in a slant well; fracture spacing and orientation from a vertical pilot well; laboratory triaxial testing and mineralogical determinations; and porosity results from magnetic resonance analyses from various wells, together with observations based on the data collection. In addition, available references were reviewed on the site for additional insights. As the focus of the report is on the natural system, hydraulic fracture data from the site were not examined in detail in this report. Data variability is the chief observation in examination of the database. Fracture frequency in the slant well can range from sections with values as high as five fractures per ft to sections up to 100+ ft in length with no natural observed fractures. Fracture spacing across is typically less than 10 ft, but can range up to hundreds of feet. Apparent fracturing shows the trends in two predominate orientations, E-SW and WNW-ESE, but minor variations exist. The rock units vary across the site from siliceous mudstones to calcareous mudstones, showing a general layering with depth. The laboratory properties such as strength and modulus show no apparent trend with depth, but appear to correlate with rock mineralogy with high strength and modulus values where calcium content is high. In addition, an attempt to examine variability and mineralogy on a larger scale was made using a color-coded system based on gamma ray measurements. A staged colored approach was adopted, presuming that lower gamma ray values indicate higher value of calcium content (blue scale) and that higher gamma ray values indicate higher clay mineral content (orange scale). As provided in report appendices, the system correlated well with visual examination of the slant core and the petrofabric analyses of the vertical pilot well. The results showed large variability in mineral content along the horizontal plane across the site. The change in mineral content was also rapid, on a scale less than that of the average hydraulic fracture stage length of about 180 ft.

42 ENGINEERING↗

NASA NTRS2009

Proceedings of the 40th Lunar and Planetary Science Conference

The 40th Lunar and Planetary Science Conference included sessions on: Phoenix: Exploration of the Martian Arctic; Origin and Early Evolution of the Moon; Comet Wild 2: Mineralogy and More; Astrobiology: Meteorites, Microbes, Hydrous Habitats, and Irradiated Ices; Phoenix: Soil, Chemistry, and Habitability; Planetary Differentiation; Presolar Grains: Structures and Origins; SPECIAL SESSION: Venus Atmosphere: Venus Express and Future Missions; Mars Polar Caps: Past and Present; SPECIAL SESSION: Lunar Missions: Results from Kaguya, Chang'e-1, and Chandrayaan-1, Part I; 5 Early Nebula Processes and Models; SPECIAL SESSION: Icy Satellites of Jupiter and Saturn: Cosmic Gymnasts; Mars: Ground Ice and Climate Change; SPECIAL SESSION: Lunar Missions: Results from Kaguya, Chang'e-1, and Chandrayaan-1, Part II; Chondrite Parent-Body Processes; SPECIAL SESSION: Icy Satellites of Jupiter and Saturn: Salubrious Surfaces; SNC Meteorites; Ancient Martian Crust: Primary Mineralogy and Aqueous Alteration; SPECIAL SESSION: Messenger at Mercury: A Global Perspective on the Innermost Planet; CAIs and Chondrules: Records of Early Solar System Processes; Small Bodies: Shapes of Things to Come; Sulfur on Mars: Rocks, Soils, and Cycling Processes; Mercury: Evolution and Tectonics; Venus Geology, Volcanism, Tectonics, and Resurfacing; Asteroid-Meteorite Connections; Impacts I: Models and Experiments; Solar Wind and Genesis: Measurements and Interpretation; Mars: Aqueous Processes; Magmatic Volatiles and Eruptive Conditions of Lunar Basalts; Comparative Planetology; Interstellar Matter: Origins and Relationships; Impacts II: Craters and Ejecta Mars: Tectonics and Dynamics; Mars Analogs I: Geological; Exploring the Diversity of Lunar Lithologies with Sample Analyses and Remote Sensing; Chondrite Accretion and Early History; Science Instruments for the Mars Science Lander; . Martian Gullies: Morphology and Origins; Mars: Dunes, Dust, and Wind; Mars: Volcanism; Early Solar System Chronology; Seek Out and Explore: Upcoming and Future Missions; Mars: Early History and Impact Processes; Mars Analogs II: Chemical and Spectral; Achondrites and their Parent Bodies; and Planning for Future Exploration of the Moon The poster sessions were: Lunar Missions: Results from Kaguya, Chang'e-1, and Chandrayaan-1; LRO and LCROSS; Geophysical Analysis of the Lunar Surface and Interior; Remote Observation and Geologic Mapping of the Lunar Surface; Lunar Spectroscopy; Venus Geology, Geophysics, Mapping, and Sampling; Planetary Differentiation; Bunburra and Buzzard Coulee: Recent Meteorite Falls; Meteorites: Terrestrial History; CAIs and Chondrules: Records of Early Solar System Processes; Volatile and Organic Compounds in Chondrites; Crashing Chondrites: Impact, Shock, and Melting; Ureilite Studies; Petrology and Mineralogy of the SNC Meteorites; Martian Meteorites; Phoenix Landing Site: Perchlorate and Other Tasty Treats; Mars Polar Atmospheres and Climate Modeling; Mars Polar Investigations; Mars Near-Surface Ice; Mars: A Volatile-Rich Planet; Mars: Geochemistry and Alteration Processes; Martian Phyllosilicates: Identification, Formation, and Alteration; Astrobiology; Instrument Concepts, Systems, and Probes for Investigating Rocks and Regolith; Seeing is Believing: UV, VIS, IR, X- and Gamma-Ray Camera and Spectrometer Instruments; Up Close and Personal: In Situ Analysis with Laser-Induced Breakdown Spectroscopy and Mass Spectrometry; Jupiter and Inscrutable Io; Tantalizing Titan; Enigmatic Enceladus and Intriguing Iapetus; Icy Satellites: Cryptic Craters; Icy Satellites: Gelid Geology/Geophysics; Icy Satellites: Cool Chemistry and Spectacular Spectroscopy; Asteroids and Comets; Comet Wild 2: Mineralogy and More; Hypervelocity Impacts: Stardust Models, LDEF, and ISPE; Presolar Grains; Early Nebular Processes: Models and Isotopes; Solar Wind and Genesis: Measurements and Interpretation; Education and Public Outreach; Mercury; Pursuing Lunar Exploration; Sources and Eruptionf Lunar Basalts; Chemical and Physical Properties of the Lunar Regolith; Lunar Dust and Transient Surface Phenomena; Lunar Databases and Data Restoration; Meteoritic Samples of the Moon; Chondrites, Their Clasts, and Alteration; Achondrites: Primitive and Not So Primitive; Iron Meteorites; Meteorite Methodology; Antarctic Micrometeorites; HEDs and Vesta; Dust Formation and Transformation; Interstellar Organic Matter; Early Solar System Chronology; Comparative Planetology; Impacts I: Models and Experiments; Impacts II: Craters and Ejecta; Mars: Volcanism; Mars: Tectonics and Dynamics; Martian Stratigraphy: Understanding the Geologic History of Mars Through the Sedimentary Rock Record; Mars: Valleys and Valley Networks; Mars: Aqueous Processes in Valles Marineris and the Southern Highlands; Mars: Aqueous Geomorphology; Martian Gullies: Morphology and Origins; Mars: Dunes, Dust, and Wind; Mars: Remote Sensing; Mars: Geologic Mapping, Photogrammetry, and Cratering; Martian Mineralogy: Constraints from Missions and Laboratory Investigations; Mars Analogs: Chemical and Physical; Mars Analogs: Sulfates and Sulfides; Missions: Approaches, Architectures, Analogs, and Actualities; Not Just Skin Deep: Electron Microscopy, Heat Flow, Radar, and Seismology Instruments and Planetary Data Systems, Techniques, and Interpretation.

Source record↗

DOE OSTI2020

A Perspective on the Successes of the NNSS Underground Test Area (UGTA) Activity - 20221

Between 1951 and 1992, 828 underground nuclear tests were performed at the Nevada National Security Site (NNSS). Underground testing was done in five major testing areas, which included (1) Frenchman Flat, (2) Rainier Mesa/Shoshone Mountain, (3) Yucca Flat/Climax Mine (4) Central Pahute Mesa, and (5) Western Pahute Mesa. About one third of the underground tests were detonated near, at or below the water table, and thus radioactive contamination was introduced to the groundwater system. The U.S. Department of Energy's Underground Test Area (UGTA) Activity was established in the late 1990's to address the fate and movement of residual radioactivity in groundwater, and characterize the risk that it may pose to NNSS workers and the offsite public. It has accomplished this goal through a process of iterative groundwater sampling and numerical groundwater flow and transport modeling, and long-term monitoring. The UGTA Activity draws on the expertise of scientists in the fields of geology, hydrology, radiochemistry, and risk assessment from the U.S. DOE staff, the lead contractor (currently Navarro) and many other organizations, in cooperation with the governing regulatory body, the State of Nevada's Division of Environmental Protection (NDEP). Los Alamos National Laboratory (LANL) is one of the participating research organizations involved with the UGTA Activity. Under the direction of DOE and the lead contractor, LANL's role in UGTA has evolved over time from an initial focus on geologic characterization, to conducting field and laboratory experiments, and finally to its current role of providing modeling and geochemistry expertise to characterize the rates and directions of groundwater and radionuclide movement. In its current role, LANL has either developed or provided the numerical tools for developing flow and transport models in each of the four major testing areas. In addition, LANL has provided an independent assessment of future groundwater flow and transport behavior through the analysis of naturally-occurring geochemical and isotopic tracers in groundwater. The five major testing areas are now in different stages of investigation or closure: (1) Frenchman Flat is in its fifth year of post closure modeling; (2) The Rainier Mesa/Shoshone Mountain Closure Report has been submitted to NDEP and is awaiting approval; (3) The Yucca Flat/Climax Mine Closure Report is under development; and (4) Central and (5) Western Pahute Mesa is still undergoing investigation and flow and transport model development. The success of the UGTA Activity in reaching or approaching closure in three of the five major testing areas is primarily due to the DoE's focus on regulatory closure requirements, but also is in part attributable to a rigorous internal review process and the involvement of NDEP and Nye County water professionals as participants and observers in the reviews. County and state involvement in the internal review process has promoted trust that the U.S. DOE is deeply committed to ensuring the present and future safely of NNSS workers and the surrounding communities. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Simulation insights into wetting properties of hydrogen-brine-clay for hydrogen geo-storage

Hydrogen geo-storage is attracting substantial interdisciplinary interest as a cost-effective and sustainable option for medium- and long-term storage. Hydrogen can be stored underground in diverse formations, including aquifers, salt caverns, and depleted oil and gas reservoirs. The wetting dynamics of the hydrogen-brine-rock system are critical for assessing both structural and residual storage capacities, and ensuring containment safety. Through molecular dynamics simulations, we explore how varying concentrations of cushion gases (CO 2 or CH 4 ) influence the wetting properties of hydrogen-brine-clay systems under geological conditions (15 MPa and 333 K). We employed models of talc and the hydroxylated basal face of kaolinite (kaoOH) as clay substrates. Our findings reveal that the effect of cushion gases on hydrogen-brine-clay wettability is strongly dependent on the clay-brine interactions. Notably, CO 2 and CH 4 reduce the water wettability of talc in hydrogen-brine-talc systems, while exerting no influence on the wettability of hydrogen-brine-kaoOH systems. Detailed analysis of free energy of cavity formation near clay surfaces, clay-brine interfacial tensions, and the Willard-Chandler surface for gas-brine interfaces elucidate the molecular mechanisms underlying wettability changes. Our simulations identify empirical correlations between wetting properties and the average free energy required to perturb a flat interface when clay-brine interactions are less dominant. Here, our thorough thermodynamic analysis of rock-fluid and fluid-fluid interactions, aligning with key experimental observations, underscores the utility of simulated interfacial properties in refining contact angle measurements and predicting experimentally relevant properties. These insights significantly enhance the assessment of gas geo-storage potential. Prospectively, the approaches and findings obtained from this study could form a basis for more advanced multiscale simulations that consider a range of geological and operational variables, potentially guiding the development and improvement of geo-storage systems in general, with a particular focus on hydrogen storage.

25 ENERGY STORAGE↗

DOE OSTI2025

Seismic H2: Version 1.0

Seismic-H2 is an integrated software package for geological hydrogen reservoir simulation, optimization, and leakage monitoring. The package includes multiple components: (1) code used for modeling seismic wave propagation in 3D heterogeneous elastic media based on finite-difference method to support detection of geological hydrogen storage reservoir leakage; (2) 3D reservoir simulations of leaks from an underground reservoir and 3D simulations of saline aquifers and depleted gas reservoirs; (3) seismic monitoring costs of passive and active seismic monitoring required for UHS; (4) rock physics calculations and interpolations for converting the reservoir simulations from part (2) into the elastic media models in part (1); (5) pre-processing seismic data; and lastly (6), a GUI interface that combines these different components.

Creasy, Neala↗

DOE OSTI2021

Stepwise Dynamic Calibration of a Hydromechanical Simulation Using Time-Lapse Vertical Seismic Profile

This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse Vertical Seismic Profiling (VSP) integration. The selected study site is the 13-10A injector well within the ongoing CO2-EOR operation of the Farnsworth Field Unit. The Time-lapse VSP dataset carries the combined effects of fluid substitution and mean effective stress changes, thereby providing a dataset amenable for the calibration of production and injection-induced stress changes. The concept is similar to calibrating a reservoir simulation model in that the process honor real field data to set up an inverse problem. The solution optimizes the independent and impactful geomechanical parameters that replicate the observed time-lapse seismic velocity changes. This stress calibration is enabled by 4D geomechanical modeling and the VSP Integration workflow. This calibration benefits from extensive geological, geophysical and geomechanical characterization through 3D seismic data, geophysical well logs, and core assessed as part of the 1D MEM conducted on the 13-10A subject well. These data are used to develop a site-specific rock physics model. The Biot Gassmann workflow combines rock physics and reservoir simulation outputs to determine the fluid substitution contribution to seismic velocity change. Additionally, modeled seismic velocity attributed to mean effective stress are determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from the ultrasonic seismic velocity measurements on the extracted Morrow B core. A penalty function is then formed between the modeled seismic velocities and the observed time-lapse VSP dataset. Four independent and impactful geomechanical parameters have been determined. These are the bulk modulus and shear modulus for zero porosity and the shear and compressional seismic velocity to mean effective stress derivatives. The dataset of numerous coupled hydromechanical- geomechanical simulation realizations is built by combining variations of the four stated geomechanical parameters. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer are used to converge on the optimal geomechanical parameters. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling.

02 PETROLEUM↗