Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geologic Modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 559 records · Page 31

Insights From Dayflow: A Historical Streamflow Reanalysis Dataset for the Conterminous United States

Abstract Reconstructed historical streamflow time series can supplement limited streamflow gauge observations. However, there are common challenges of typical modeling approaches: process‐based hydrologic models can be data/computation‐intensive, and statistics‐based models can be region/stream‐specific. Here we present a nationally scalable modeling framework integrating the simulated runoff from the Variable Infiltration Capacity (VIC) model with the Routing Application for Parallel computatIon of Discharge (RAPID) routing model leveraging high‐performance computing. We demonstrate an efficient method of assimilating streamflow at US Geological Survey (USGS) streamflow monitoring sites using a simple hierarchical approach in the VIC‐RAPID framework. The result is a reconstructed 36‐year (1980–2015) daily and monthly streamflow dataset (Dayflow) at ∼2.7 million NHDPlusV2 stream reaches in the conterminous US (CONUS). We perform a comprehensive evaluation at 7,526 USGS sites and characterize their error statistics. The results demonstrate that 49% of the USGS sites demonstrate Kling–Gupta Efficiency (KGE) > 0.5 and 58% of the sites show percentage bias within ±20% for the daily naturalized streamflow. Streamflow data assimilation across CONUS shows an overall improvement over naturalized streamflow, notably in the western semiarid‐to‐arid regions. Comparison to other national and global streamflow reanalysis datasets such as the National Water Model and Global Reach‐scale A priori Discharge Estimates for SWOT demonstrates improved KGE, reduced bias, and directions for Dayflow improvements. Investigations of error statistics with key hydrologic, hydroclimatic, and geomorphologic basin characteristics reveal region‐specific patterns which may help improve future framework applications. Overall, Dayflow may enable a better understanding of hydrologic conditions in a changing environment, especially in locations currently not represented by streamflow monitoring networks.

54 ENVIRONMENTAL SCIENCES↗

Modeling Cost of Offshore Carbon Storage in Saline Reservoirs

Poster presentation to the AAPG-SEG-SPE CCUS 2024 Conference. Offshore saline reservoirs provide a significant and accessible resource for geologic carbon storage (CS). The offshore environment requires distinct approaches to site selection, operations, monitoring, and risk that affect the technoeconomic assessment of offshore CS projects. The National Energy Technology Laboratory (NETL) has developed a CS cost model for offshore saline reservoirs known as CO2_S_COM_Offshore. Based on NETL’s widely used CO2_S_COM cost model for onshore saline CS, CO2_S_COM_Offshore enables technoeconomic analysis of CS in offshore areas. This model comprehensively incorporates multiple facets of offshore CS projects, from regional evaluation and site selection to permitting, transport, operations, monitoring, site closure, and decommissioning. Developed to model cost for offshore United States (US) Exclusive Economic Zones, aspects of this model can be adapted to international projects. Presented at the SPE/AAPG/SEG Carbon Capture Utilization and Storage Conference in Houston, TX, March 11-13, 2024.

Mark-Moser, Mackenzie K.↗

The Impact of Ethylenediaminetetraacetic Acid (EDTA) on the Sorption of Nd(III), Th{sup (IV)}, and U(VI) onto Dolomite in WIPP-Relevant Brines, GWB and ERDA-6 - 20356

Storage of legacy, transuranic waste at the Waste Isolation Pilot Plant (WIPP) requires accurate risk models for performance assessments. These models must predict the long-term fate and transport of contaminants in the WIPP and surrounding environment. Laboratory support is necessary for the development of an accurate model through better understanding of contaminant behavior in the WIPP. The WIPP is a deep geologic repository that is characterized by high ionic strength due to evaporite deposits in the Salado formation (porewaters up to 7.4 M). High concentrations of sodium chloride (NaCl), calcium chloride (CaCl{sub 2}), and magnesium chloride (MgCl{sub 2}) are particularly of concern due to their potential impact on solubility and ultimately transport of contaminants of concern. Additionally, ligands capable of forming strong complexes with metals are found in significant concentrations in the WIPP environment due to their disposal from waste processing and formation during degradation of repository components (e.g. cellulose degradation by calcium hydroxide in cement). Consequently, there is a need to investigate the fate and transport of contaminants in the WIPP environment with a focus on the effects of high ionic strength and strongly coordinating ligands. Contaminants of high concern include the actinides (U(VI), Pu(IV), and Am(III)). Laboratory support of the WIPP modeling efforts was provided by this research. A series of batch experiments were utilized to study the impact of ionic strength, contaminant concentration, and the presence of ethylenediaminetetraacetic acid (EDTA) on the sorption of Nd(III), Th(IV), and U(VI) onto dolomite, a carbonate mineral found within the Culebra formation. The Culebra formation is the most transmissive layer above the WIPP, making transport of contaminants the most likely route in previous risk assessments. Two brines specifically designed to imitate conditions within the WIPP environment were used in these experiments. The U.S. Energy Research and Development Administration Well 6 (ERDA-6, 95%) brine simulates a low Mg environment, while the generic weep brine (GWB, 95%) simulates the WIPP environment with a high Mg concentration. Mg is added to the repository in the form of MgO as an engineered barrier and is naturally present in some formations. It is expected to scavenge carbonate from solutions decreasing formation of the highly soluble and competitive uranium-carbonate complexes. EDTA is used during manipulations of nuclear materials due to its ability to form strong, stable complexes with the actinides and may be found in high concentrations (up to 0.08 mM) within the WIPP as waste canisters degrade. Further, it may provide additional avenues for unwanted transport of contaminants. Nd(III) and Th(IV) were used as stable analogues of Am(III) and Pu(IV), respectively, to represent the most common oxidation states of actinides in the WIPP environment in addition to U(VI). The potential for colloid formation was also investigated with low and high concentrations of contaminants (10 and 1000 μg/L) using three size-dependent particle separation steps. The results of this research provide a better understanding of potential behavior of actinides in the presence of EDTA and within WIPP-relevant brines. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NRAP-Open-IAM: Open-Source Integrated Assessment Model

NRAP Open IAM is the latest evolution in a series of quantitative models developed under the U.S. Department of Energy’s National Risk Assessment Partnership (NRAP) to address risk assessment, risk management and containment assurance questions of geological carbon storage (GCS). NRAP Open IAM allows a user to define a site specific GCS scenario, characterize important site features, events, and processes (considering uncertainty), couple together reduced order models of various system components to rapidly forecast response to CO 2 injection, and perform stochastic analysis to quantify containment effectiveness, leakage risk, and impact to receptors of concern.

Vasylkivska, Veronika S.↗

Modeling Aerial Gamma-Ray Backgrounds Using Non-negative Matrix Factorization

Airborne gamma-ray surveys are useful for many applications, ranging from geology and mining to public health and nuclear security. In all these contexts, the ability to decompose a measured spectrum into a linear combination of background source terms can provide useful insights into the data and lead to improvements in the techniques that use spectral energy windows. Multiple methods for the linear decomposition of spectra exist but are subject to various drawbacks, such as allowing negative photon fluxes or requiring detailed Monte Carlo modeling. In this work, we propose using non-negative matrix factorization (NMF) as a data-driven approach to spectral decomposition. Using aerial surveys that include flights over water, we demonstrate that the mathematical approach of NMF finds physically relevant structure in the aerial gamma-ray background, namely, that measured spectra can be expressed as the sum of nearby terrestrial emission, distant terrestrial emission, and radon and cosmic emission. Furthermore, these NMF background components are compared with the background components obtained by noise-adjusted singular value decomposition (NASVD), which contain negative photon fluxes and, thus, do not represent the emission spectra in as straightforward a way. Finally, we comment on the potential areas of research that are enabled by NMF decompositions, such as new approaches to spectral anomaly detection and data fusion.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Model Data Archive for Manuscript Titled "Evaluation of a Coupled Surface–Subsurface Hydrologic Model Using Dense Water‑Level Sensors in a Mixed Urban–Rural Watershed"

This archive provides scripts, input files, and datasets used for the implementation and evaluation of a fully coupled surface–subsurface hydrologic model in the Neches River Basin, southeast Texas. The study uses the Advanced Terrestrial Simulator (ATS) to simulate coupled surface–subsurface hydrologic processes over a mixed urban–rural watershed and evaluates model performance using a dense network of 136 in situ water-level sensors, nine U.S. Geological Survey (USGS) stream gauges, and SSEBop-derived evapotranspiration estimates during the period October 2014–June 2024. The workflow is implemented primarily in Python 3 using the Watershed Workflow package. The Jupyter notebooks can be executed using open-source software such as Anaconda JupyterLab or Visual Studio Code. Other data files include TXT, CSV, XML, SHP, TIF, NetCDF, HDF5, and ExodusII files, which can be processed using the provided Python scripts. ATS input files are provided in XML format and can be edited using any commonly used text editor. This archive contains: *Scripts and input files used to generate the ATS model setup, including watershed discretization, mesh generation, parameter mapping, and model configuration. *Jupyter notebooks used for preprocessing observational data, evaluating streamflow, water levels, and evapotranspiration, computing performance metrics, and generating the figures presented in the manuscript. *ATS simulation outputs and processed observational datasets, including OneRain and DD6 water-level sensors, USGS streamflow observations, GIS data, and supporting spatial datasets used throughout the study.

Dense water-level sensor network↗

Evaluation and Optimization of Well Completion Options for the Utah FORGE Site

Orientation and completion for well pairs that have been subjected to multi-zonal stimulation play a critical role in the long-term performance of an Enhanced Geothermal Reservoir. Enhanced geothermal systems often rely on preferential flow along fractures between well injection and production locations. Modeling this preferential flow using discrete fracture networks (DNF) relies on stochastic realizations of the DFN based on geological sampling. Here we present the development of a stochastic optimization methodology to determine well completion options in a discrete fracture network based on using parallel subset simulation. Stochastic optimization will provide insight into regions where placements of the injection and production wells are optimal. An example optimization of well-pair location optimization based on a deterministic-stochastic DFN model representing FORGE follows a discussion of the theory.

15 GEOTHERMAL ENERGY↗

Post-fire soil emissions of nitric oxide (NO) and nitrous oxide (N 2 O) across global ecosystems: a review

Wildfires may increase soil emissions of trace nitrogen (N) gases like nitric oxide (NO) and nitrous oxide (N 2 O) by changing soil physicochemical conditions and altering microbial processes like nitrification and denitrification. When 34 studies were synthesized, we found a significant increase in both NO and N 2 O emissions up to 1 year post-fire across studies spanning ecosystems globally. However, when fluxes were separated by ecosystem type, we found that individual ecosystem types responded uniquely to fire. Forest soils tended to emit more N 2 O after fire, but there was no significant effect on NO. Shrubland soils showed significant increases in both NO and N 2 O emissions after fires; often with extremely large but short-lived NO pulses occurring immediately after fire. Grassland NO emissions increased after fire, but the size of this effect was small relative to shrublands. N 2 O emissions from burned grasslands were highly variable with no significant effect. To better understand the variation in responses to fire across global ecosystems, more consistent measurements of variables recognized as important controls on soil fluxes of NO and N 2 O (e.g., N cycling rates, soil water content, pH, and substrate availability) are needed across studies. We also suggest that fire-specific elements like burn severity, microbial community succession, and the presence of char be considered by future studies. Our synthesis suggests that fires can exacerbate ecosystem N loss long after they burn, increasing soil emissions of NO and N 2 O with implications for ecosystem N loss, climate, and regional air quality as wildfires increase globally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Application of machine learning for modeling brønsted-guggenheim-scatchard specific ion interaction theory (SIT) coefficients

Machine learning methodologies can provide insight into Brønsted-Guggenheim-Scatchard specific ion interaction theory (SIT) parameter values where experimental data availability may be limited. This study develops and executes machine learning frameworks to model the SIT interaction coefficient, ε. Key findings include successful estimations of ε via artificial neural networks using clustering and value prediction approaches. Additionally, applicability to other chemical parameters is also assessed briefly. Models developed here provide support for a use-case of machine learning in geologic nuclear waste disposal research applications, namely in predictions of chemical behaviors of high ionic strength solutions (i.e., subsurface brines).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Toward more-robust, AI-enabled subsurface seismic imaging for geotechnical applications

Non-invasive seismic imaging has the potential to cost-effectively evaluate large volumes of subsurface material to inform geotechnical site investigation. However, seismic imaging using full waveform inversion (FWI) requires significant computational time and is dependent on an initial starting model. As a result, FWI has not yet been widely adopted into geotechnical practice. Previous efforts, on relatively simple two-layered models, indicate that data-driven artificial intelligence (AI) models may be as effective as FWI at predicting 2D images of shear wave velocity (V s ). Furthermore, the AI model predictions can be made almost instantaneously after data acquisition and do not require an initial starting model. We examine the generality of these findings by developing a new AI model for subsurface seismic imaging, whereby we make several notable contributions. First, we architect a multimodal AI model that combines time- and frequency-domain representations of the seismic wavefield to predict a 50 m by 20 m subsurface image of V s . Second, we developed a new diverse dataset of 100,000 images with their corresponding seismic wavefields to train the AI model. Third, we propose four physics-informed data augmentations for data-driven seismic imaging. Fourth, we develop two prediction consistency tests to evaluate the model’s performance when the true subsurface is unknown. Our final model, which has been made publicly available, is capable of predicting a subsurface V s image from a single seismic wavefield with an average, mean absolute percent error (MAPE) of 24 %. The predictive model is applied to a field dataset and shown to be consistent with local geology and shear-wave refraction measurements from the same location.

Artificial intelligence↗

Patterns and controls on island‐wide aboveground biomass accumulation in second‐growth forests of Puerto Rico

Abstract Understanding the heterogeneity of biomass accumulation in second‐growth tropical forests following land use abandonment is important for informing ecosystem carbon models and forest restoration efforts. There is an urgent need for a broad sample of second‐growth forests to enhance our knowledge of carbon accumulation in human‐dominated landscapes, especially for older forests. Puerto Rico has predominantly second‐growth forests, ranging in age from approximately 25 to more than 80 years. We used an island‐wide sample of airborne lidar from the NASA Goddard Lidar, Hyperspectral, and Thermal (G‐LiHT) Airborne Imager collected on March 2017, forest inventory data, and data on forest age, precipitation, soils, and land use to estimate aboveground biomass stocks in moist and wet, second‐growth tropical forests. Biomass accumulation rates in Puerto Rico were lower, on average, than in other Neotropical forests. Median biomass across >16,700 ha of older second‐growth forests was 105 Mg ha −1 , and sampled biomass rarely surpassed 250 Mg ha −1 . Differences in biomass by age were large and persistent across different substrates and land uses, with a plateau in the pattern of island‐wide biomass accumulation after about 33 years. A spatial regression model showed that multiple factors were related to biomass accumulation, including time since abandonment, geologic substrate, past land use as coffee or pasture, precipitation, topographic wetness index, and slope. Our findings have important consequences for the total carbon storage and expected climate mitigation benefits of large‐scale reforestation efforts, and highlight the value of airborne lidar for quantifying biomass variability in complex tropical landscapes. Abstract in Spanish is available with online material.

54 ENVIRONMENTAL SCIENCES↗

Shallow EGS Regional Resource Potential and Map - Snake River Plain

SMU Geothermal Lab developed a methodology to estimate shallow (1 km to 4 km) Enhanced Geothermal Systems (EGS) resource potential using an approach that utilizes recent geology and geophysical research along with new well data to improve the thermal conductivity model, mitigate impacts from groundwater flow in the thermal model, and examine radioactivity data variations. By incorporating the results of the most recent projects with the SMU shallow methodology, we developed a more accurate, updated resource estimate for the Snake River Plain (SRP). The resulting maps and resource estimates can be used by the National Renewable Energy Lab (NREL), Bureau of Land Management (BLM), and the public to determine how best to move forward with future project development.

15 GEOTHERMAL ENERGY↗

WHOLESCALE Catalog of Rock Samples at San Emidio Nevada collected in January 2021 Version 2.0

This submission includes an update to the WHOLESCALESamples2021January.xls file where the strikes and dips initially reported have been corrected to comply with the right hand rule. The updated excel file and a link to the original submission are included in this report. The original submission "WHOLESCALE Catalog of Rock Samples at San Emidio Nevada collected in January 2021" contains information on thirty-six rock samples collected from San Emidio, Nevada during January, 2021 for Subtask 2.3 of the WHOLESCALE project. The following resources include a .zip of rock sample photos taken in the field, a .zip of rock sample photos taken in the laboratory at UW-Madison, and an excel catalog of rock samples with information on sample name, rock type, coordinates of sample location, structural measurements, field notes, observations for plug preparation (e.g., weathering, ability to be cut and cored), and rock descriptions. It should be noted that not every sample was photographed in the field. Names and descriptions of rock formation units are taken from Rhodes et al. (2011). The README.txt file is a description of this submission.

15 GEOTHERMAL ENERGY↗

Prestack Waveform Inversion at GC 955: Trials and sensitivity of PWI to high-resolution seismic data (Data Report)

Understanding the spatial extent and variability of hydrate deposits has important implications for potential economic production, climate change, and assessing natural hazards risks. Seismic reflection techniques are capable of determining the extent of gas hydrate deposits and can shed light on hydrate concentrations and properties. Using high resolution prestack time migrated seismic data collected during a 2013 US Geological Survey expedition and prestack waveform inversion (PWI) techniques, we produce highly resolved velocity models and compare them to co-located well logs. Coupling our PWI results with velocity-porosity relationships and nearby well control, we compare hydrate properties at site 955 in the Green Canyon and site 313 at Walker Ridge.

03 NATURAL GAS↗

Decovalex-2019 (Task A Final Report)

The DECOVALEX Project is an on-going international research collaboration, established in 1992, to advance the understanding and modeling of coupled Thermal (T), Hydrological (H), Mechanical (M) and Chemical (C) processes in geological in geological systems. This document is the final report of Task A which was proposed by the DECOVALEX Project and coordinated by the British Geological Survey (BGS), presenting the technical definitions of the problems studied, approaches applied, achievements made and outstanding issues for future research. The purpose of Task A under DECOVALEX-2019 is to better understand the processes governing the advective movement of gas in both low-permeability argillaceous repository host rocks and clay-based engineered barriers. Special attention is given to the mechanisms controlling gas entry, flow and pathway sealing and their impact on the performance of the engineered clay barrier. Previous work suggests gas flow may be accompanied by the creation of dilatant pathways whose properties change temporally and spatially within the medium. Thus, new numerical representations for the quantitative prediction of gas migration fluxes through argillaceous rock formations have been developed. These provide an invaluable tool with which to assess the impact of gas flow on repository layout and therefore design of any future facility. In addition, experience gained through this task is of direct relevance to other clay-based engineering issues where immiscible gas flow is a consideration including shale gas, hydrocarbon migration, carbon capture and storage and landfill design.

58 GEOSCIENCES↗

DECOVALEX-2019 (Task G Final Report)

The DECOVALEX Project is an on-going international research collaboration, established in 1992, to advance the understanding and modeling of coupled Thermal (T), Hydrological (H), Mechanical (M) and Chemical (C) processes in geological in geological systems. This document is the final report of Task G which was proposed by the Swedish Radiation Safety Authority (SSM) and coordinated by geomecon GmbH, presenting the technical definitions of the problems studied, approaches applied, achievements made and outstanding issues for future research. Task G in DECOVALEX-2019 is primarily concerned with the evolution of transmissivity throughout the lifetime of a repository for radioactive waste, and in particular of spent nuclear fuel, in sparsely fractured and competent crystalline rock masses, and the representation of these processes via numerical simulation. The numerical simulations were to be evaluated for their suitability, and to be validated. Over the course of the work presented here, not only were transmissivity changes and the best approaches to simulate those change examined, but also strategies were to be developed on how to monitor transmissivity change. Guidelines on what a control program and a monitoring system for a repository should include and be able to measure, and how this can be done in practice, was a key objective of the research project.

58 GEOSCIENCES↗

DECOVALEX-2019 (Task D: INBEB Final Report)

The DECOVALEX Project is an on-going international research collaboration, established in 1992, to advance the understanding and modeling of coupled Thermal (T), Hydrological (H), Mechanical (M) and Chemical (C) processes in geological in geological systems. This document is the final report of Task D which was proposed by the DECOVALEX Project and coordinated by CIMNE_UPC, presenting the technical definitions of the problems studied, approaches applied, achievements made and outstanding issues for future research. Task D of the DECOVALEX 2019 project is devoted to the study of the hydro-mechanical (HM) and thermo-hydro-mechanical (THM) Interactions in Bentonite Engineered Barriers (INBEB). Special attention is paid to the evolution of barrier heterogeneity under transient conditions. A proper understanding of the processes occurring during the transient period requires the appropriate simulation of the behavior of the engineered barrier by means of coupled numerical analyses. In this context, the main objective of the Task is to assess the capabilities of numerical formulations and codes to simulate and interpret the HM and THM behavior of the bentonite barrier, with a particular focus on the final state reached on saturation.

58 GEOSCIENCES↗