Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geoscience”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

High-precision measurement of U-Pu-Np-Am concentrations and isotope ratios in environmental reference materials by mass spectrometry

Herein, we report results of precise and sensitive mass spectrometric measurements of uranium, plutonium, neptunium, and americium concentrations and isotope ratios in a variety of environmental reference materials. Most of our work has been done on NIST SRM 4350b, River Sediment, but we also present results for NIST SRM 4354, Lake Sediment; NIST SRMs 4355 and 4355a, Peruvian Soil; NIST SRM 4357, Ocean Sediment; NIST SRM 1648a, Urban Particulate Material; NIST SRM 1649b, Urban Dust; IAEA CRM 385, Ocean Sediment; USGS BCR-2, Columbia River Basalt; and USGS BHVO-2, Hawaiian Volcanic Observatory Basalt. These materials reflect a wide range in long-lived actinide concentrations (e.g. 1E4 to 1E10 atoms 239 Pu/g) and isotope ratios. Measurements were performed in a clean laboratory by isotope dilution, multi-collector thermal ionization and multi-collector inductively coupled plasma mass spectrometry. In general, our results are in agreement with, but lower the uncertainty of, literature or certificate values for these reference materials. Our uranium results for the basalts also confirm previously reported high-precision mass spectrometric results from our laboratory. In many cases our measurements of U-Pu-Np-Am nuclides appear to be novel.

58 GEOSCIENCES↗
DOE OSTI2021

Iron-rich microstructure records of high temperature multi-component silicate melt behavior in nuclear fallout

Above-ground nuclear explosions that interact with the surface of the earth entrain materials from the surrounding environment, influencing the resulting physical and chemical evolution of the fireball, which can affect the final chemical phase and mobility of hazardous radionuclides that are dispersed in the environment as fallout particles. The interaction of iron with a nuclear explosion is of specific interest due to the potential for iron to act as a redox buffer and because of the likelihood of significant masses of metals to be present in urban environments. We investigated fallout from a historic surface interacting nuclear explosion conducted on a steel tower and report the discovery of widespread and diverse iron-rich micro-structures preserved within the samples, including crystalline dendrites and micron-scale iron-rich spheres with liquid immiscibility textures. We assert these micro-structures reflect local redox conditions and cooling rates and can inform interpretation of high temperature events, enabling new insights into fireball condensation physics and chemistry when metals from the local environment (i.e. structural steel) are vaporized or entrained. Finally, these observations also significantly expand the availability of silicate immiscibility datasets applicable to rapidly quenched systems such as meteorite impact melt glass.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Characterizing major and trace element compositions in fallout melt glass from a near-surface nuclear test

We report the chemical and isotopic compositions of fallout melt glasses from nuclear tests contain a range of information constraining the physical conditions within the fireball and the mechanisms of fallout formation but historic studies tended to exclude the behavior of stable major and trace elements. Here, we present a large study specifically focused on major and trace element relationships within a population of macroscale fallout samples from a single event. We interpret these data to better constrain how fallout melt glass formation in near surface environments is influenced by that environment and demonstrate how major and trace element abundances can provide useful insights into chemical processes within the fireball. Data confirm that the uranium in the fallout glass population derives from two isotopically distinct endmembers: isotopically enriched uranium (presumably from the weapon), and natural composition uranium that may be a combination of anthropogenic and environmental materials from within the blast zone. The similarity between major and trace element concentrations in fallout and corresponding local soils from the event site confirm the local soils as the most probable source of entrained material into the fireball and the source of carrier material into which the bomb vapor was incorporated. The lack of correlation between major and trace element abundances with size indicates that volatility driven processes, such as condensation from the fireball, do not control the composition of macroscale fallout melt glass. Although the fallout has major and trace element chemical characteristics broadly similar to those of the local, associated soils, some systematic differences are observed between the two populations. Fallout melt glass is depleted in volatile elements such as K, Na, Tl and Pb, consistent with heating to temperatures above ~1000 °C for 3-10 s. This is supported by the results of laser heating experiments performed on rhyolitic soil at temperatures (1600-2200 °C) and timescales (1-120 s) that are broadly relevant to fallout formation conditions. Relative enrichments of metals such as Cu and Co do not correlate with the abundance of uranium, suggesting that fallout also records input of near field anthropogenic materials. Our observations suggest that major chemical features can be related to processing in the fireball and used to inform the thermal-chemical evolution of the system. Ultimately, these data are consistent with a fallout formation mechanism that involves rapid melting of surface materials to form carrier material melts with minor incorporation of bomb vapor and a degree of volumetric volatile loss due to heating.

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

Noble gas adsorption to tuff

For this report a method was developed to measure trace noble gas element adsorption to the surfaces of geologic materials in the presence of a background gas that could potentially compete for surface adsorption sites. Adsorption of four noble gas elements (Ne, Ar, Kr, and Xe) at a concentration of 100 ppm in helium and nitrogen were measured on a sample of crushed tuff at 0, 15, 30, and 45°C. In addition, Ne, Ar, Kr, and Xe at 250 ppm and 500 ppm in nitrogen at 15°C were measured. Noble gas adsorption was found to increase with increasing atomic mass and decreasing temperature. It was also observed that the relative increase in noble gas element adsorption with decreasing temperature tends to increase with increasing atomic mass. Noble gas adsorption from nitrogen was approximately an order of magnitude greater than for the noble gases in helium. As the noble gas concentrations in nitrogen increased, adsorption increased in a slightly non-linear fashion which could be modeled using a Freundlich isotherm. For noble gas concentrations that were ≤100 ppm Henry's Law constant were calculated.

58 GEOSCIENCES↗
DOE OSTI2025

Preferential adsorption of noble gases in zeolitic tuff with variable saturation: A modeling study of counter-intuitive diffusive-adsorptive behavior

Noble gas transport through geologic media has important applications in the prediction and characterization of measured gas signatures related to underground nuclear explosions (UNEs). Retarding processes such as adsorption can cause significant species fractionation of radionuclide gases, which has implications for measured and predicted signatures used to distinguish radioxenon originating from civilian nuclear facilities or from UNEs. Accounting for the effects of variable water saturation in geologic media on tracer transport is one of the most challenging aspects of modeling gas transport because there is no unifying relationship for the associated tortuosity changes between different rock types, and reactive transport processes such as adsorption that are affected by the presence of water likewise behave differently between gas species. In this study, we perform numerical diffusive-adsorptive transport simulations to estimate gas transport parameters associated with bench-scale laboratory diffusion cell experiments measuring breakthrough in zeolitic and non-zeolitic rocks for a gaseous mixture of xenon, krypton, and SF 6 at varying degrees of water saturation (S w ). Counter-intuitive transport behavior was observed in the zeolitic rock experiments whereby breakthrough concentrations were significantly higher when the core was partially saturated (S w = 17 %) than under dry (S w = 0 %) conditions. Breakthrough of xenon was especially retarded in the dry core – likely due to comparatively high affinity of xenon for zeolitic adsorption sites – and estimated effective diffusion coefficients for all gases were approximately an order of magnitude lower than what is predicted by porosity-tortuosity models. We propose the counter-intuitive behavior observed is because water infiltration of zeolite nanopores reduces both the adsorptive capacity of the rock and the tortuosity of connected flow paths. We developed a two-site competitive kinetic Langmuir adsorption reaction for the porous media transport simulator in order to constrain transport parameters within zeolitic tuff, where differential adsorption to zeolite and non-zeolite pores was observed. We determined that liquid saturation-dependent diffusive-adsorptive transport is affected by subtle and at times competing processes that are specific to different gases, which have a significant overall influence on effective transport parameters.

58 GEOSCIENCES↗
DOE OSTI2025

A statistical approach to screening isotopic signatures in monitoring for underground nuclear explosions

The ability to differentiate between atmospheric radionuclide signatures from underground nuclear explosions (UNEs) and signals from other sources, such as medical isotope-production facilities and nuclear reactors, can be critical to the detection and monitoring of unannounced, low-yield nuclear events. Signatures having anomalously high amplitudes, compared to background levels, remain the best indicator in screening for a UNE. However, isotopic composition can further validate a suspected UNE signature, but separation from any atmospheric background composition is first necessary. To date, evaluating the challenges of performing this separation has typically involved comparing an observed background with a highly idealized deterministic model of radioxenon signature production by a UNE that does not consider the influence of post-detonation chemical/physical processes in the detonation cavity or the subsequent gas transport mechanisms that can also affect the isotopic composition of the detected gas signature. In addition, purely deterministic models, as previously employed, overlook the uncertainty inherent in estimating critical parameters characterizing the UNE and its detonation environment. In this paper, we create detailed, multi-parameter models of radionuclide evolution using the widely accepted England and Rider post-detonation radionuclide decay-chain network coupled to detailed models simulating physical production and transport processes affecting the gas signature. Because these models are governed by uncertain parameters including barometric fluctuations, realistic ranges of variation for each of the parameters influencing isotopic composition are then defined. A Latin-Hypercube sampling approach is used to obtain a random distribution of isotopic production and gas transport results associated with a given value of each parameter. We apply these results to background histories of two stations, one providing 4-isotope background measurements and the other providing two-isotope measurements associated with the 2013 DPRK announced UNE.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2024

Dynamic modeling studies of basin-scale pressure interference and CO 2 plume evolution in multi-well geologic CO 2 storage

This research employs numerical modeling to analyze how CO 2 plumes and pressure fronts evolve when CO 2 is injected into a single storage formation from multiple projects located in proximity. Here, this analysis also evaluates the extent to which injection well spacing alleviates pressure buildup in the absence of active pressure management tactics. The simulation approach was based on a single, homogenous saline aquifer in which CO 2 injection occurs under a one-injector baseline case and several multi-well cases where well spacing varies. Analysis results show that the extent of pressure buildup is in the range of tens or a few hundreds of kilometers and contingent upon the defining pressure buildup demarcating the front edge. For the geological setting evaluated in this paper, our analysis suggests that without active basin pressure management strategies, commercial-scale projects would likely need to be sited far apart to avoid pressure interference from one another. Analysis results show the radius of CO 2 plume varies approximately from 2 to 3 km from injection wells (each injecting 1 Mt/year for 30 years) depending on cases and modeling parameters assumed. Given the pressure interference, this paper thus draws attention to the importance of greater coordination among storage operators and regulatory stakeholders. Because this analysis assumes a very specific geologic setting, this exploratory analysis bears further investigations across other geologic scenarios.

42 ENGINEERING↗
DOE OSTI2024

Mineralogy and reactive fluid chemistry evolution of hydraulically fractured Caney shale of Southern Oklahoma

Here, this study investigates geochemical rock-fluid interactions as a potential cause of rapid loss of permeability and productivity in hydraulically fractured shale reservoirs. It also interrogates the effects of these reactions in transforming depleted shale reservoirs into impermeable carbon storage units. The study employs batch reactor experiments where rock-powder samples are reacted with field fracturing fluid under reservoir temperature (95°C).

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

Spectral induced polarization (SIP) measurements across a PFAS-contaminated source zone

There is a pressing need for the development of field-scale, in situ screening technologies for assessing variations in aqueous film forming foam (AFFF) concentrations in soils at former fire training and storage sites. Field-scale Spectral Induced Polarization (SIP) geophysical measurements were acquired on a transect crossing an AFFF source zone. Soil samples were acquired from ten locations and used to determine variations in poly- and per-fluoroalkyl substances (PFAS) concentrations in soils and soil texture. These samples were also used to create triplicate soil columns for laboratory-grade SIP measurements. Field and laboratory observations provide evidence that SIP measurements are sensitive to the concentration of AFFF constituents associated with the pore surface in soils. The phase of the SIP measurements on the laboratory samples was linearly correlated with total soil-sorbed PFAS concentration. The phase from the field SIP measurements was highest over the location of maximum PFAS concentration measured on the laboratory samples, although a significant correlation between field-measured phase and laboratory-measured total PFAS concentration was not established. The sensitivity of the SIP response to the removal of soil PFAS using a methanol wash procedure (total PFAS concentration drop of 366 ppb) adds evidence for the case for SIP characterization of AFFF source zones. The results of these studies suggest that SIP might be developed into a field-scale technology for rapid, indirect assessment of AFFF source zones. Such a technology could improve the effectiveness of AFFF source zone characterization at reduced costs.

58 GEOSCIENCES↗
DOE OSTI2023

Potential Benefits in Remapping the Special Flood Hazard Area: Evidence from the U.S. Housing Market

A typical U.S. homebuyer's understanding of whether a property faces flood risk is based on whether the property is located inside the National Flood Insurance Program's (NFIP) Special Flood Hazard Area (SFHA). The SFHA boundary, however, may bias homebuyers’ perceptions of flood risk relative to unobserved true risk because the SFHA is an incomplete, and sometimes inaccurate, representation of flood hazards. Using a national dataset of property transactions, flood hazard data from the First Street Foundation, state-level measures of flood disclosure laws, and a spatially restricted triple-difference design, we distinguish capitalization effects of policy-driven (SFHA) versus quasi-objective (First Street) indicators of flood hazard. We identify these effects by assessing thousands of local spatial interactions between property SFHA designations, measures of quasi-objective flood hazard, and being in a state that mandates flood history disclosures. Being inside the SFHA generates a risk discount, but the signal is muted relative to underlying hazard exposure. Further, the SFHA signal can result in inefficient discounts for properties erroneously mapped in the SFHA. Disclosure requirements about flood history accentuate price effects for hazardous properties and result in more adequate risk internalization. In the absence of disclosures, however, we find either no or a weak risk signal for houses outside the SFHA that face flood hazard. Our results highlight potential benefits of updating SFHA boundaries to include all houses that may experience flooding, and argue in favor of requiring flood-related disclosures from sellers to improve the market's ability to internalize flood risk.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A new soil mixing layer model for simulating conservative solute loss from initially saturated soil to surface runoff

Surface water pollution due to solute loss from soil to surface runoff is a serious threat to the environment, and it is of great significance to predict the solute loss. While soil mixing layer theory has been widely used to predict the solute loss, it has been observed that the conventional soil mixing layer models with constant mixing layer depth underestimate solute concentration in surface runoff, and thus underestimate the solute loss after surface runoff starts. This study introduced a new soil mixing layer model with a new concept that the mixing layer depth increases downward with time after surface runoff occurs. This new concept was implemented by adding a mixing front between the soil mixing layer and its underlying soil layer, which is conceptually similar to the wetting front of soil infiltration given by the Green-Ampt equation. This is equivalent to treating the soil mixing layer as a deforming control volume. A mechanistic explanation for the new model was presented. This study also developed a closed form solution of the new model. The new model was evaluated by simulating five laboratory experiments conducted by two groups of researchers. The first two experiments have ponding water, and consider different drainage and runoff conditions at both early and late runoff times. The other three experiments do not have ponding water, and consider different infiltration conditions at the bottom of soil boxes. Simulation results of the five experiments show that the new model resolves the underestimation problem of conventional soil mixing layer models. The new model is subject to a total of six limits such as that the new model only considers conservative solute. Three limits are related to the assumptions used to derive the closed form solution, and they may be alleviated in a future research as briefly discussed at the end of the paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Characterization of flow and transport in a fracture network at the EGS Collab field experiment through stochastic modeling of tracer recovery

Energy extraction from subsurface reservoirs is important for addressing the increasing energy demand and environmental concerns such as global warming. However, the characterization of subsurface reservoirs, particularly reservoirs dominated by fracture networks remains a challenge due to the lack of means to directly observe subsurface processes. This study explores the feasibility and efficacy of characterizing fracture flow and transport processes in an enhanced geothermal system (EGS) testbed through stochastic tracer modeling. There are two enabling factors that allow application of stochastic modeling to characterize a subsurface reservoir. First, an abundance of geological and geophysical measurements enables the development of a high-fidelity and well-constrained fracture network model. Second, high-performance computing (HPC) allows running massive realizations efficiently. Six conservative tracer tests were stochastically modeled and produced satisfactory realizations that successfully reproduce field tracer recovery data from each tracer test. The evolution of flow and transport processes in the fracture network was then analyzed from these satisfactory realizations. This work demonstrates that stochastic tracer modeling on a high-fidelity fracture network model is feasible and can provide important insights regarding flow and transport characteristics in subsurface fractured reservoirs.

58 GEOSCIENCES↗
DOE OSTI2021

Can artificial intelligence and data-driven machine learning models match or even replace process-driven hydrologic models for streamflow simulation?: A case study of four watersheds with different hydro-climatic regions across the CONUS

With recent developments in computational techniques, Data-driven Machine Learning Models (DMLs) have shown great potential in simulating streamflow and capturing the rainfall-runoff relationship in given watersheds, which are traditionally fulfilled by Process-based Hydrologic Models (PHMs). There are debates on whether the DMLs can outperform and possibly replace the classical PHMs for streamflow simulation and river forecasting, but no clear conclusions have been made. This study aims to investigate whether the newer DMLs have any potential in further improving the simulation accuracy of classical PHMs, and vice versa. To do this, we compared a few popular PHMs and DMLs over four watersheds across the Continental US (CONUS) that are associated with different input, climate, and regional conditions. A total of five hydrologic models were chosen, including (1) two classical lumped models, i.e., the Sacramento Soil Moisture Accounting (SAC-SMA) and Xinanjiang (XAJ); (2) one modern distributed model, termed Coupled Routing and Excess Storage (CREST); (3) and two DMLs including an Artificial Neural Networks (ANN) and a deep learning model, termed Long Short Term Memory (LSTM). Our results demonstrated that the DMLs still significantly biased when using the baseline input scenario with the PHMs. However, the DMLs fed with delayed input scenarios had great potential and can reach high simulation accuracy. The DMLs, especially the ANN, outperformed other employed models under the rainfall-runoff relationship in which rainfall dominantly drives. Furthermore, the DMLs also showed better performance in the high-flow regime, while the PHMs had a better performance for the low-flow regime, implying both PHMs and DMLs have their own merits and are worthy of joint development. In general, our study indicated a great potential of using DMLs to simulate streamflow, but further studies are still needed to verify the transferability and scalability of DMLs in large-scale experiments, such as the Distributed Model Intercomparison Projects 1&2 conducted by National Weather Services but to compare modern DMLs and PHMs.

58 GEOSCIENCES↗
DOE OSTI2021

Geostatistical interpolation of streambed hydrologic attributes with addition of left censored data and anisotropy

Spatial geostatistical interpolation of point measurements of streambed attributes in the hyporheic zone may be constrained by the streambed anisotropy, and data density and spatial distribution may significantly impact the results. Spatial clustering and low spatial data density can be caused by bedrock outcropping at the streambed limiting installation of in-stream piezometers. This study examines parameter error variability of the geostatistical interpolation using anisotropic interpolation methods and increasing the data density by adding left censored values (i.e., data below measurement limit) to locations where measurements were limited by exposed bedrock lining the streambed. The reduction in relative standard error of the interpolation was determined for the spatial distributions of streambed attributes including hydraulic conductivity, seepage flux, and mercury solute flux measured in two different years along a study reach in East Fork Poplar Creek, Tennessee, USA. Here, two methods to impute the left censored values were compared including the conventional half the detection limit substitution method, and the Stochastic Approximation of Expectation-Maximization (SAEM) algorithm, which both had comparable results. Imputing left censored data increased the data density to recommended ranges, reduced data clustering, increased the spatial dependence for some attributes, and reduced the standard error for each of the three attributes. For the reach considered herein, addition of the left censored values resulted in a larger error reduction than the consideration of anisotropy within the interpolation, which confirms the benefit of data addition to increase data density within data-limited river corridors.

58 GEOSCIENCES↗
DOE OSTI2021

Insights into storm direction effect on flood response

In this study, we investigate the directional influence of storm movement on catchment flood peak response using the synthetic circular basin. Due to the complexity in defining storm movements that require meteorological modeling, we adopt a novel approach of combining the basin rotation method (BRM) with a circular basin construct. A systematic basin rotation approach provides a proxy for the storm direction relative to the river network topology to study its hydrologic response. Using Stage-IV rainfall with 4-km by 1-hour resolution for 16 years period from 2003 to 2018 in Iowa, U.S., we systematically analyze consequent flood peak response due to storm directions (basin rotations) using a distributed hydrologic model called the Hillslope-Link Model (HLM). HLM has demonstrated the ability to reproduce observed streamflow for real river basins studied in Iowa. The hydrologic simulation results show that the BRM can quantify the relationship between the rainstorm direction and catchment hydrologic response. Also, the maximum flood peak response is strongly dependent on the interaction of the peak of the rainstorm with the peak of the width function. The results indicate significant differences in runoff volume and runoff peak and their interannual variability due to rainstorm direction. The study has important implications for flood frequency predictions and developing flood resilience and watershed management strategies, particularly due to changing storm tracks under a warming climate.

42 ENGINEERING↗
DOE OSTI2022

Multiscale petrophysical characterization and flow unit classification of the Minnelusa eolian sandstones

Integration of petrophysical and geological information is critical to simulation of subsurface carbon storage (GCS). In this sense, two depositional facies were identified from the core description and well-log interpretation, namely massive (MS) and cross-bedded (CB) facies groups. Additionally, pore-scale characteristics were studied by a combination of techniques, e.g. Nuclear Magnetic Resonance (NMR) and mercury intrusion capillary pressure (MICP). Scanning electron microscope (SEM) and petrographic analyses show that the pore structure is dominantly controlled by the depositional environment and dolomite cementation. NMR-T 2 distributions of MS and CB facies show triple and quadruple modes, respectively. In addition, MICP of high- and low-permeability MS facies samples, and their CB facies group mixtures were collected. The MS sample pore-throat size distribution is uni-modal, while the triple-modal characteristic of the mixtures indicates heterogeneous pore structures at the sub-core scale for CB facies. The reliably estimates of porosity and permeability for both facies groups via NMR techniques and the MLR (Multiple Linear Regression) approach demonstrate the applicability of these techniques to eolian sandstone. Moreover, irreducible water saturation via the T 2 -cutoff method correlates strongly with T 2LM instead of porosity. Finally, the rock quality index and flow zone indicator were calculated based on Combinable Magnetic Resonance (CMR) log interpretations. This provides direct connection to properties measured in the well. Four flow units were classified for both facies groups. Finally, results show that better reservoir quality with significant heterogeneities is observed in the CB facies. This study highlights the importance integrating multiscale petrophysical properties including facies, pore architecture and diagenesis analysis with core- to log-scale property characterization. The results herein validate our reservoir characterization and flow unit classification in eolian reservoirs.

58 GEOSCIENCES↗
DOE OSTI2022

A robust deep learning workflow to predict multiphase flow behavior during geological C O 2 sequestration injection and Post-Injection periods

Simulation of multiphase flow in porous media is essential to manage the geologic CO 2 sequestration (GCS) process, and physics-based simulation approaches usually take prohibitively high computational cost due to the nonlinearity of the coupled physics. This paper contributes to the development and evaluation of a deep learning workflow that accurately and efficiently predicts the temporal-spatial evolution of pressure and CO 2 plumes during injection and post-injection periods of GCS operations. Based on a Fourier Neural Operator, the deep learning workflow takes input variables or features including rock properties, well operational controls and time steps, and predicts the state variables of pressure and CO 2 saturation. To further improve the predictive fidelity, separate deep learning models are trained for CO 2 injection and post-injection periods due to the difference in primary driving force of fluid flow and transport during these two phases. We also explore different combinations of features to predict the state variables. We use a realistic example of CO 2 injection and storage in a 3D heterogeneous saline aquifer, and apply the deep learning workflow that is trained from physics-based simulation data and emulate the physics process. Through this numerical experiment, we demonstrate that using two separate deep learning models to distinguish post-injection from injection period generates the most accurate prediction of pressure, and a single deep learning model of the whole GCS process including the cumulative injection volume of CO 2 as a deep learning feature, leads to the most accurate prediction of CO 2 saturation. For the post-injection period, it is key to use cumulative CO 2 injection volume to inform the deep learning models about the total carbon storage when predicting either pressure or saturation. The deep learning workflow not only provides high predictive fidelity across temporal and spatial scales, but also offers a speedup of 250 times compared to full physics reservoir simulation, and thus will be a significant predictive tool for engineers to manage the long-term process of GCS.

58 GEOSCIENCES↗
DOE OSTI2023

Pore connectivity influences mass transport in natural rocks: Pore structure, gas diffusion and batch sorption studies

For this work, six rocks (one granodiorite, one limestone, two chalks, one mudstone, and one dolostone) with different extents of heterogeneity at six different particle sizes (from 75 to 8000 μm) were studied to describe the effects of pore connectivity on mass transport. The methods applied were (i) porosity measurement of granular rocks, (ii) analyses of gas-phase diffusive transport in a bed of packed particles, along with a solid quartz method at these six particle sizes being developed to identify the contribution of intraparticle diffusion, and (iii) batch sorption tests of multiple ions (anions and cations) with subsequent analyses of inductively coupled plasma-mass spectrometry. Granular porosity measurement results reveal that with decreasing particle sizes, the effective porosities for the “heterogenous” group of rocks (Grimsel granodiorite and Edwards limestone) increase, whereas the porosities of another “homogeneous” group (two Israel chalk samples, Japan mudstone, and Wyoming dolostone) remain constant. Gas diffusion results show that the intraparticle gas diffusion coefficient among these two sample groups, varying in the magnitude of 10 -8 to 10 -6 m 2 /s, are not directly correlated to the porosity differences. Moreover, the batch sorption work displays a different affinity of rocks for various tracers. For Grimsel granodiorite, Japan mudstone, and Wyoming dolostone, the adsorption capacity of Sm 3+ and Eu 3+ increases as the particle size decreases. In general, this integrated research of grain size distribution, granular rock porosity, intraparticle diffusivity, and ionic sorption capacity gives insights into the pore connectivity effect on both physical and chemical transport behaviors for different lithologies and/or different particle sizes.

58 GEOSCIENCES↗