Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geologic characterization”

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 487 records · Page 27

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↗

Linking impact melt redox with crustal weathering regime

Abstract Impact melting and outgassing of impact melts were important processes shaping the crust and atmosphere of early Earth and Mars. The redox state of impact melts influences the gases they release, but the controls on impact melt redox are not well understood. Here, we present Fe micro–X-ray absorption near edge structure (μ-XANES) measurements of a suite of impact melt glasses from Lonar crater, India. As an impact into continental flood basalts, Lonar crater is an ideal analog for impacts into basaltic crust on early Earth and Mars. The Fe-μ-XANES technique permits robust characterization of the Fe3+/ΣFe ratio and therefore redox conditions in Lonar glasses. We found a range in Fe3+/ΣFe of 0.21–0.49 among Lonar glasses, including substantial intraclast variations. We conclude that a primary factor driving these variations was pre-impact oxidative weathering of target rocks. The coupling between crustal weathering regime and speciation of gases released from impact melts implies that impact melt outgassing could act as a feedback mechanism to reinforce prevailing surface redox conditions.

Geology↗

Deep Learning for Full Waveform Inversion of Elastic Active-Source Seismic Data to Estimate P-Wave Velocity Models

Seismic imaging methods are critical for Global Security and Energy & Homeland Security missions and activities that rely on subsurface characterization, but traditional methods remain computationally expensive and require significant labor hours and expertise to execute. Within the past few years, machine learning (ML), namely deep learning (DL), has been used to develop data-driven end-to-end full waveform inversion (FWI) methods to estimate 2D P-wave velocity (Vp) models in a fraction of the time as conventional FWI. These methods, however, are trained on simplistic acoustic wave seismic data and Vp models that are not realistic nor representative of real-world observations, leaving a large gap between the state-of-the-art and deployable, feasible, and practical DL FWI methods. Here, we generate a synthetic active-source, 3D, elastic wave seismic data set and a variety of Vp models with realistic geologic structure for training DL FWI methods. We evaluate six different methods that have performed well for acoustic DL FWI or medical imaging tasks using our more realistic dataset. We find that these six trained models do not match the performance of published acoustic end-to-end DL FWI methods, indicating more training data may be needed, physics may need to be incorporated to achieve good accuracy at the sacrifice of the end-to-end advantage, and/or novel methods need to be developed to enable end-to-end DL FWI methods to perform well for real-world seismic data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Borehole Disposal of Radioactive Waste in Israel - 20388

Israel is assessing borehole disposal of radioactive waste. With limited geological options for disposal, intermediate-depth borehole disposal is being considered in the arid Yamin Plain region of the northeastern Negev desert at depths of several hundred meters below ground surface in the vadose zone. Unlike deep borehole disposal of several kilometers, which relies on emplacement below the depth of recirculating groundwater, the safety case for intermediate-depth borehole disposal relies more on the aridity of the vadose zone, the robustness of the waste package, and the other engineered barriers (e.g., seals, backfill materials) in the disposal borehole. As a tradeoff, the shallower depth may allow the use of a larger diameter borehole than would be possible for a deeper borehole. To study the suitability of the Yamin Plain region for borehole disposal, a small-diameter characterization borehole is being planned to retrieve core samples and to better understand the vadose zone geomechanical and hydrogeochemical properties and percolation flux. The information from the characterization borehole will inform the safety case, and together with performance assessment analyses will help to identify key areas of uncertainty and guide future research and development activities aimed at demonstrating the feasibility of the intermediate-depth borehole disposal concept in Israel. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Probabilistic risk-based Area of Review (AoR) determination for a deep-saline carbon storage site

Regulatory oversight of a geologic carbon sequestration (GCS) project relies on iterative estimations, throughout the project lifetime, of the area where increased risks to underground sources of drinking water (USDWs) may occur due to injection of CO2. This area, referred to as Area of Review (AoR), is typically delineated by predicting the migration of fluid between the reservoir and the lowermost USDW via an open wellbore using predictions from physics-based reservoir simulators. The inherent uncertainty in input parameters used in reservoir modeling therefore affects the accuracy of determining the AoR for a project. Furthermore, the standard analytical approaches for calculating a critical pressure to delineate the risk area yield an infinite AoR for cases where the injection reservoir is overpressured relative to the USDW. A methodology is presented in this paper to better characterize the risk to USDWs while accounting for the uncertainty in reservoir modeling, with an application to a permitted GCS project with an overpressured injection formation, FutureGen 2.0. The methodology is demonstrated using the National Risk Assessment Partnership’s open-source integrated assessment model (NRAP-Open-IAM) to develop a probabilistic estimate of impact risk to USDW quality. CO2 and pressure predictions from the reservoir modeling conducted using the STOMP-CO2 simulator for the FutureGen 2.0 site are used in a NRAP-Open-IAM model with reservoir, wellbore, and aquifer components to: (1) assess the extent of potential leakage into the USDW for the predicted reservoir pressure conditions; (2) evaluate the extent of potential impact using “no-net-degradation” thresholds; and (3) account for uncertainty in reservoir permeabilities.

03 NATURAL GAS↗

Single-Particle Isotopic Analysis Using the Liquid Sampling-Atmospheric Pressure Glow Discharge Microplasma Coupled to an Orbitrap Mass Spectrometer

Isotope ratio (IR) determinations on individual particles can provide valuable information relative to the sample, including processing and dating, which could be valuable to geological and nuclear material analysis communities. In comparison to “bulk” measurements, in which all particle information is lost and homogenized within a sample, “particle” measurements allow for fingerprinting and can provide unique insights related to the sample’s nano- and micro- compositions. Furthermore, the complex nature of particles, with potential isobaric interferences from either the matrix or other elements within the same particle, poses a significant analytical challenge for conventional mass spectrometry platforms employed for elemental analysis due to their limited mass resolution. Presented here is the first demonstration of an ultrahigh resolution method for single particle (SP) isotopic analysis using the liquid sampling-atmospheric pressure glow discharge (LS-APGD) microplasma ionization source coupled to an Orbitrap mass spectrometer and FTMS Booster X2T acquisition/processing system. For proof of concept, well-characterized CeO 2 microparticles with a nominal diameter of 1.0 μm were analyzed at a mass resolution of ∼330,000. Important instrumental parameters were optimized to enable the detection of single particles. The 142 Ce/ 140 Ce ratio of individual particles and the population average were determined and compared to an SP inductively coupled plasma time-of-flight mass spectrometry (ICP-TOF-MS) analysis. The isotopic accuracy and precision were determined by two methods and found to be in good agreement with the SP-ICP-TOF-MS method, with the linear regression slope method providing a more accurate ratio, showing a ∼1.4% relative difference from the ratio determined from a particle digest. The encouraging performance of the method was further supported by a determined detection limit of 2.9 fg ( 142 Ce). The developed method is anticipated to overcome many of the challenges posed by isobaric interferences encountered in conventional mass spectrometry techniques when analyzing real-world samples for particle isotopic composition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unsupervised Clustering of Microseismic Events and Focal Mechanism Analysis at the CO 2 Injection Site in Decatur, Illinois

Characterization of induced microseismicity at a carbon dioxide (CO 2 ) storage site is critical for preserving reservoir integrity and mitigating seismic hazards. We apply a multilevel machine learning (ML) approach that combines the nonnegative matrix factorization and hidden Markov model to extract spectral representations of microseismic events and cluster them to identify seismic patterns at the Illinois Basin-Decatur Project. Unlike traditional waveform correlation methods, this approach leverages spectral characteristics of first arrivals to improve event classification and detect previously undetected planes of weakness. By integrating ML-based clustering with focal mechanism analysis, we resolve small-scale fault structures that are below the detection limits of conventional seismic imaging. Our findings reveal temporal bursts of microseismicity associated with brittle failure, providing insights into the spatio-temporal evolution of fault reactivation during CO 2 injection. This approach enhances seismic monitoring capabilities at CO 2 injection sites by improving fault characterization beyond the resolution of standard geophysical surveys.

Willis, Rachel Marie [Sandia National Laboratories↗

Defining Infrastructure Feasibility for Hub-Scale Offshore Atlantic Carbon Storage in the Northeastern United States

In the Northeast U.S., deep rock formations along the Atlantic outer continental shelf may have the potential to sequester 150–1136 million metric tons of CO2. However, the design and infrastructure necessary to develop offshore carbon storage in this region is not well defined because there has been little oil and gas exploration and no commercial production. Consequently, an infrastructure feasibility design was completed for a hub-scale offshore CO2 storage system along the Northeast U.S. Atlantic. The design included development of a detailed, site-specific geological model for a location near the Great Stone Dome geological structure in the Baltimore Canyon Trough off the coast of Delaware, Maryland, and New Jersey. A field injection system topology design was completed to portray a design with eight wells in two clusters connected by central manifolds. Reservoir simulations were completed for the injection system that showed the hub may be able to inject 17 million metric tons (MMT) of CO2 per year for thirty years, but injection rates varied substantially across the eight wells. A CO2 pipeline design determined feasible routes from the east coast shoreline to the injection field. Finally, the CO2 injection system design included subsea injection trees, manifolds, and power umbilicals. This is the first study to define large-scale carbon storage design and infrastructure options for the offshore Atlantic, which can help to progress this region towards field characterization and early-mover deployment for future decarbonization in the region.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Surrogate models for development of unconventional shale reservoirs by an integrated numerical approach of hydraulic fracturing, flow and geomechanics, and machine learning

We develop well-completion surrogate models by taking an integrated workflow of hydraulic fracturing, flow, geomechanics, and machine learning simulation. There are three steps in the proposed workflow. First, history-matching processes are conducted with the field data including pumping and production data for characterization. Second, full-physics simulation is performed with various parameters of the field development (e.g., cluster spacing, clusters per stage, pumping rates and times, amount of proppant, and well spacing) to generate multiple simulation results by changing the parameters of the completion design with well-known hydraulic fracturing, reservoir, geomechanics simulators to calculate fracture geometry, reservoir depressurization, induced stress changes. The workflow is demonstrated over a field in the Southern Midland Basin. Here, we take two completion scenarios: a single well case followed by a multi-well case. Finally, a Long Short-Term Memory (LSTM) machine learning algorithm is employed to create surrogate models that can replicate the full-physics simulation results. Furthermore, results show that the trained models applied in the single well and multi-well cases for a particular geological system can provide good accuracy close to those provided by full-physics simulations. Specifically, the site-specific surrogate models can predict fracture parameters (length, height, and surface area) and cumulative production accurately with computational efficiency, suggesting our proposed workflow can be used as a pragmatic tool for expediting the well completion optimization process.

Geomechanics↗

Augmentation of WRF-Hydro to simulate overland-flow- and streamflow-generated debris flow susceptibility in burn scars

In steep wildfire-burned terrains, intense rainfall can produce large runoff that can trigger highly destructive debris flows. However, the ability to accurately characterize and forecast debris flow susceptibility in burned terrains using physics-based tools remains limited. Here, we augment the Weather Research and Forecasting Hydrological modeling system (WRF-Hydro) to simulate both overland and channelized flows and assess postfire debris flow susceptibility over a regional domain. We perform hindcast simulations using high-resolution weather-radar-derived precipitation and reanalysis data to drive non-burned baseline and burn scar sensitivity experiments. Our simulations focus on January 2021 when an atmospheric river triggered numerous debris flows within a wildfire burn scar in Big Sur – one of which destroyed California's famous Highway 1. Compared to the baseline, our burn scar simulation yields dramatic increases in total and peak discharge and shorter lags between rainfall onset and peak discharge, consistent with streamflow observations at nearby US Geological Survey (USGS) streamflow gage sites. For the 404 catchments located in the simulated burn scar area, median catchment-area-normalized peak discharge increases by ~ 450 % compared to the baseline. Catchments with anomalously high catchment-area-normalized peak discharge correspond well with post-event field-based and remotely sensed debris flow observations. We suggest that our regional postfire debris flow susceptibility analysis demonstrates WRF-Hydro as a compelling new physics-based tool whose utility could be further extended via coupling to sediment erosion and transport models and/or ensemble-based operational weather forecasts. Given the high-fidelity performance of our augmented version of WRF-Hydro, as well as its potential usage in probabilistic hazard forecasts, we argue for its continued development and application in postfire hydrologic and natural hazard assessments.

54 ENVIRONMENTAL SCIENCES↗

Geomechanical Characterization of the Mount Simon Sandstone and Eau Claire Formation of Northern Illinois Basin

The purpose of this paper is to analyze publicly available geomechanical data from two carbon storage sites in the Illinois Basin: the Illinois Basin-Decatur Project (IBDP) in Macon County, IL and the FutureGen2.0 project in Morgan County, IL, cancelled before injection began. This study estimates the magnitudes and directions of the least compressive principal stress gradient using hydraulic fracture-based stress measurement tests, borehole breakouts, and sonic logs in the Eau Claire Formation (primary confining zone) and the Mount Simon Sandstone (targeted reservoir). The range of stresses expected in the underlying Precambrian basement is also investigated. The evaluation of the state of stress is carried out in a probabilistic manner using the State of Stress Analysis Tool (SOSAT) developed under the National Risk Assessment Partnership (NRAP) program. Finally, the paper also provides a discussion on distinctions between fracture initiation pressure, propagation pressure, and fracture closure pressure, given that hydraulic fracture testing results show a significant difference between these three values. The discussion also addresses why the fracture closure pressure should be preferred in injection pressure determinations, a critical point in geological carbon storage projects.

de Toledo Camargo, Julia↗

On the Interplay Between Distributed Bulk Plasticity and Local Fault Slip in Evolving Fault Zone Complexity

Abstract We numerically investigate the role of plastic strain accumulation on the mechanical response of a planar strike‐slip fault. Our models show that fault‐zone strength significantly impacts the ensuing sequence of earthquakes. Weaker fault zones accumulating more plastic strain promote more complexity in the seismicity pattern through aperiodic earthquake occurrences and intermittent episodes of rupture and arrest. However, if the fault zone strength is high enough, the overall earthquake sequence is characterized by periodic fault‐spanning events. We find that both the fault normal stress and the fault geometric profile evolve throughout the earthquake sequence, suggesting a self‐roughening mechanism. Despite the significant impact of plasticity on the fault response, the width of the plastically deforming region in the fault zone is small compared to the fault length. Our results suggest a rich behavior in dynamically evolving fault zones and support the need for further high‐resolution studies of the highly non‐linear near‐fault region.

Geology↗

National‐Scale Assessment of Total Gaseous Mercury Isotopes Across the United States

Abstract With the 2011 promulgation of the Mercury and Air Toxics Standards by the U.S. Environmental Protection Agency, and the successful negotiation by the United Nations Environment Program of the Minamata Convention, emissions of mercury (Hg) have declined in the United States. While the declines in atmospheric Hg concentrations in North America are encouraging, linking the declines to changing domestic and global source portfolios remains challenging. To address these research gaps, the U.S. Geological Survey initiated the first national‐scale effort to establish a baseline of total gaseous mercury stable isotope values at 31 sites distributed across the United States. Results indicated that unique Hg sources, such as Hg evasion from an elemental Hg contaminated site or free tropospheric intrusions in high altitude sites, were distinguishable from background atmospheric values. Minor gradients were observed across the nation, with regions of heavy industrial activity demonstrating lower δ 202 Hg, but no consistent changes in other isotopes such as Δ 199 Hg and Δ 200 Hg were observed. Furthermore, δ 202 Hg was impacted by foliar uptake and senescence but trends varied between forested regions in the northeastern and midwestern United States. These data demonstrate regional emission sources and other environmental variables can impact total gaseous Hg (TGM) isotope values, highlighting the need to characterize atmospheric Hg isotopes over larger geographical areas to evaluate changes related to national and international Hg regulations.

Tate, Michael T.↗

Williston Basin Carbon Ore, Rare Earth, and Critical Minerals (CORE-CM) Program

The U.S. Department of Energy’s Office of Fossil Energy awarded 13 CORE-CM programs as part of the CORE-CM Initiative, designed to develop the technology and upstream and midstream supply chains to extract rare-earth elements (REEs) and critical minerals (CMs) from the nation’s coal supplies. The intent is to catalyze regional economic growth and job creation while strengthening the use of domestic resources. The Williston Basin CORE-CM Program aims to drive the expansion and transformation of coal and coal-based resource usage within the Williston Basin to produce REEs, CMs, and nonfuel carbon-based products (CBPs). The work constitutes Phase 1 of a long-term program. The objectives of Phase 1 are to identify the existing knowledge base and gaps and to develop a series of assessments/plans, including an initial basinal assessment; a characterization and data acquisition plan; a waste stream reuse plan; a basinal strategies assessment for infrastructure, industries, and business; a technology assessment, development, and field testing plan; a technology innovation center plan(s); and a stakeholder outreach and education plan. The goal of this project is to initiate the development of a new industry in the resource-rich Williston Basin that will catalyze economic growth and job creation through a coalition team of private industry, university, and state, local, and federal government entities. A coalition team of nearly 30 private industry, university, and state, local, and federal government partners was formed for this program. Partners include the University of North Dakota Energy & Environmental Research Center, Institute for Energy Studies, and Nistler College; North Dakota State University; Montana Tech University; Pacific Northwest National Laboratory; Critical Materials Institute; North Dakota Geological Survey; South Dakota Geological Survey; U.S. Geological Survey; North American Coal Corporation; BNI Energy; Basin Electric Power Cooperative; Minnkota Power Cooperative; and many more.

Kay, John P.↗

Development of a Conditioned System-Level Groundwater Model to Evaluate Long-Term Groundwater Impacts within a Performance Assessment - 20115

A preliminary performance assessment (PA) of single-shell tank Waste Management Area (WMA) A-AX located at the U.S. Department of Energy (DOE)'s Hanford Site in southeastern Washington is being conducted to satisfy the requirements of DOE Order 435.1 [1] as it relates to closure of the single-shell radioactive waste tanks in the WMA A-AX tank farms. A PA assesses the fate, transport, and impacts of radionuclides within a low-level radioactive waste disposal facility in its assumed closure configuration and the subsequent potential doses to humans over a 1,000-year compliance period and a 10,000-year performance evaluation period. The WMA A-AX preliminary PA evaluation is structured around the complementary use of process-level and system-level models to calculate the facility performance against established DOE Order 435.1 [1] performance objectives. Process-level models are those that represent a detailed phenomenological representation of processes of concern in the PA. Process models typically only represent one or a few of the components of the PA, such as groundwater flow and transport, and must be integrated with other modeling elements to perform PA calculations. System-level models are those that are abstracted from the process models, retaining the essential features of the process model, while allowing integration of all aspects of the PA in a single modeling framework. System-level models are often characterized by coarser numerical discretization, lower dimensionality, or other similar simplifications compared to the process-level model. Traditionally, a three-dimensional (3-D) process-level model is utilized primarily to evaluate the long-term impact on groundwater and the potential doses to individuals who consume contaminated groundwater. System-level models are also utilized to evaluate the groundwater pathway in PAs. These models typically have reduced dimensionality (1-D) and are conditioned utilizing flow fields (Darcy fluxes) and moisture content distributions that are abstracted from the process level models. The abstraction approach assures that the flow field in both models is consistent for a specific set of input parameters for flow, differing only in the discretization and dimensionality of the two models. The preliminary WMA A-AX PA incorporates a detailed representation of the geological system and hydraulic properties within the 3-D model STOMP{sup C} numerical code so that the effects of relevant features and processes on water flow and radionuclide transport in the subsurface can be evaluated. The complementary system-level model is developed utilizing the GoldSim{sup C} code to implement a simplified, 1-D equivalent model to represent the groundwater pathway. Contaminant transport through the vadose zone and unconfined aquifer for Tc-99 and I-129 were evaluated in each of the models. Adjustments to the saturated portion of the GoldSim{sup C} 1-D model were required to mimic the dispersion effect captured with the 3-D STOMP{sup C} model. Once this conditioning was conducted, highly similar results for the transport of Tc-99 and I-129 were achieved at the point of calculation, located 100 m downgradient from the WMA A-AX fenceline. These radionuclides represent elements that are regarded as primary dose drivers in the PA groundwater pathway analysis. The high degree of conformance between the two models suggests that the use of the equivalent 1-D system model is suitable for evaluating the full suite of radionuclides that will be released from the tank sources within WMA A-AX over the 1,000-year compliance period and over the 10,000-year evaluation period. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Utah FORGE: Powder X-ray Diffraction Data from Well 16A(78)-32 Core

This dataset from Lawrence Livermore National Laboratory (LLNL) consists of four raw X-ray diffraction (XRD) scans and preliminary results of quantitative XRD analysis. The scanned samples were prepared from four subcores, which came from various depths of the FORGE well 16A(78)-32 core. Desired core lengths were selected from available core photos (on GDR), provided by FORGE personnel, and subcored at LLNL. The XRD scans were collected in May 2023 at LLNL as pre-experimental characterization data for these subcores, which will be used in core-flooding experiments at LLNL and in triaxial direct shear experiments at Los Alamos National Laboratory as part of DOE Project 5-2428. XRD scans are in RAW file format (e.g., FORGE-5477-full.raw) and are suitable for viewing and analysis using open-source quantitative XRD software (e.g., Profex; www.profex-xrd.org ) and/or other proprietary instrument software.

15 GEOTHERMAL ENERGY↗

A Play-Based Exploration of CO 2 Storage in the Illinois Basin: Community Baseline Assessment for Carbon Storage Engagement Planning in the Illinois Basin

This report presents a stakeholder mapping analysis conducted as part of the Play-Based Exploration of CO 2 Storage in the Illinois Basin project (DOE Cooperative Agreement DE-FE0032366). The analysis was undertaken to systematically identify, categorize, and characterize stakeholder groups and community conditions relevant to the planning and deployment of carbon capture and storage (CCS) infrastructure across the Illinois Basin. Its purpose is to establish a contextual baseline that informs targeted engagement strategies, supports equitable project planning, and strengthens the social dimensions of CCS site screening.

01 COAL, LIGNITE, AND PEAT↗

Conditional Pseudo-Reversible Normalizing Flow for Surrogate Modeling in Quantifying Uncertainty Propagation

We introduce a conditional pseudo-reversible normalizing flow (PR-NF) that directly learns conditional probability distributions from noisy physical models to efficiently quantify both forward and inverse uncertainty propagation. Traditional surrogate modeling approaches approximate only the deterministic component of physical models, requiring separate noise characterization and computationally expensive sampling methods for inverse problems. Here, in this work, we develop the conditional PR-NF model to directly learn and efficiently generate samples from the conditional probability density functions (PDFs). The training process utilizes dataset consisting of input-output pairs without requiring prior knowledge about the noise and the function. Once trained, our model efficiently generates samples from conditional PDFs for any input within the training domain. Moreover, the pseudo-reversibility feature allows for the use of fully connected neural network architectures, which simplifies the implementation and enables theoretical analysis. We provide a rigorous convergence analysis of the conditional PR-NF model, showing its ability to converge to the target conditional PDF using the Kullback−Leibler divergence. To demonstrate the effectiveness of our method, we apply it to several benchmark tests and a real-world geologic carbon storage problem.

97 MATHEMATICS AND COMPUTING↗