Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sandstone”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Time‐And‐Space Averaging Applied to Intermittent Multiphase Flow Experiments

Abstract Various researchers have studied fluctuations in pore‐scale phase occupancy during multiphase flow in porous media using synchrotron‐based X‐ray microcomputed tomography (micro‐CT). However, the impact of these fluctuations on the concept of a representative volume is not yet fully understood. In this study, we performed spatial and temporal averaging of multiphase flow experiments visualized with synchrotron‐based micro‐CT, focusing on oil saturation as the key parameter to determine a representative time‐and‐space average. Our findings revealed that a saturation value representative of both time and space was achieved during fractional flow experiments in drainage mode with fractional flows of 0.8, 0.5, and 0.3. Furthermore, we computed a range of relative permeabilities on the basis of whether momentaneous saturation or time‐and‐space averaged saturation was utilized for direct simulation. Our results highlighted the importance of time‐and‐space averaging in determining a representative relative permeability and indicated that the temporal and spatial scales covered in a typical micro‐CT flow experiment were sufficient to obtain a representative saturation value for sandstone rock under intermittent flow conditions.

Environmental Sciences & Ecology↗
DOE OSTI2024

Measurement of the Static Nonlinear Third-Order Elastic Moduli of Rocks: Problems and Applicability

The third-order elastic (TOE) model has been used to describe the widely observed nonlinear mechanical behaviors of earth materials. In addition to linear elastic constants ( λ , μ ), three nonlinear elastic moduli ( A , B , C ) are required for isotropic rocks. Contrary to previous research on dynamic TOE moduli, this study followed the protocol to measure strain and stress under uniaxial and hydrostatic compressive tests statically, which were later used to invert for the full set of TOE moduli for four standard rock types with differing pore structures of a porous oolitic limestone, quartz-rich sandstones, and a dense crystalline basalt. The applicability of the TOE model to characterize nonlinearity depends on the fulfillment of path-independence and small-strain assumptions. Using the measured static TOE moduli, the finite element model demonstrates that the stress in the vicinity of the wellbore is more amplified than the stress in the linear elastic case, which leads to a wider zone of rock failure around the wellbore. Due to the long-standing discrepancy between static and dynamic moduli, the rarely reported full set of static TOE moduli is necessary and will benefit future research in understanding the effect of rock nonlinearity on geophysical and geomechanical applications, such as long-term safe storage of CO 2 and generating process of geohazards.

58 GEOSCIENCES↗
DOE OSTI2024

Accelerating Multiphase Simulations With Denoising Diffusion Model Driven Initializations

This study introduces a hybrid fluid simulation approach that integrates generative diffusion models with physics‐based simulations, aiming at reducing the computational costs of flow simulations while still honoring all the physical properties of interest. Pore‐scale simulations enhance our understanding of applications such as assessing hydrogen and storage efficiency in underground reservoirs. Nevertheless, they are computationally expensive and the presence of non‐unique solutions can require multiple simulations within a single geometry. To overcome the computational cost hurdle, we propose a method that couples generative diffusion models and physics‐based simulations. While training the data‐driven model, we simultaneously generate initial conditions and perform physics‐based simulations using these. This integrated approach enables us to receive real‐time feedback on a single compute node equipped with both CPUs and GPUs. By efficiently managing these processes within a single compute node, we can continuously monitor performance and halt training once the model meets the specified criteria. To test our model, we generate realizations in a real Berea sandstone fracture which shows that our technique is up to 4.4 times faster than commonly used flow simulation initializations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Permeability and Induced Polarization of Mudstones

Electrical measurements can be used to estimate hydraulic properties such as permeability ( k ) in sedimentary rocks. Previous work has focused on sandstones, siltstones, and carbonates, while investigations on mudstones have rarely been reported. In this study, we report on electrical geophysical measurements for 23 mudstone samples using an experimental approach designed to reliably saturate these low permeability mudstones. The modified Hagen-Poiseuille model linking permeability to the formation factor ( F ) and an effective pore radius ( r ) provides an excellent fit to the data set with a near-constant pore radius, indicating that the effective porosity (1/ F ) is the controlling factor on k . In these samples, the surface area normalized to pore volume ( S por ), frequently used in permeability estimation models, varies by 1–2 orders of magnitude and is thus not a reliable proxy of the inverse effective hydraulic radius. The formation factor also exerts the primary control on induced polarization (IP) parameters, whereas Spor shows no relation to the IP parameters. A strong linear relationship is found between IP parameters (imaginary conductivity and normalized chargeability) and surface conductivity, although the proportionality factor is significantly lower than those observed in more permeable rocks and sediments. Apparent relationships between the polarization strength-derived and time constant-derived geophysical length scales and the effective hydraulic radius appear to be driven by variations in the electrochemical parameters (i.e., specific polarizability and diffusion coefficient). Overall, these findings emphasize that predicting hydraulic properties from electrical measurements in fine-grained rocks remains challenging and requires further investigation into the electrochemical properties involved.

58 GEOSCIENCES↗
DOE OSTI2025

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗
DOE OSTI2024

Optimization of floodwater redistribution from Lake Nasser could recharge Egypt’s aquifers and mitigate its excessive floods

Extreme precipitation periods, possibly related to climate change, over the Nile River source areas caused flooding in Sudan and excess runoff reaching Lake Nasser in Egypt in 1998–2002 and 2019–2022. Excess water from the 1st event (25.5 × 10 9 m 3 ) was channeled to depressions within the plateau west of the Nile Valley, forming the Tushka Lakes, where it was left to evaporate, a fate that awaits the 53.5 × 10 9 m 3 from the 2nd event while the underlying fossil Nubian Sandstone Aquifer System is being depleted (–0.98 × 10 9 m 3 /yr). We simulated release scenarios of excess Lake Nasser waters (53.5 × 10 9 m 3 ) to proximal lowlands; preference was given to the scenario that recharged the aquifer through infiltration (74.3%) and minimized losses to evaporation (20.1%) and surface runoff (5.6%). Findings serve as an example of adaptations that replace catastrophic consequences of climate change with beneficial and sustainable development opportunities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

AnisONet: A deep neural operator-based anisotropic permeability upscaler from pore to Darcy scale

Directional permeability variations, which govern directional fluid flow in porous media with anisotropy, are important to accurately predict flow behavior, reactive transport, and fluid–solid interactions for various processes such as enhanced geothermal systems, energy storage devices, and biological systems. However, the intricate architecture of porous media makes it difficult to predict directional permeabilities. In this work, we present a novel machine learning (ML) framework, AnisONet, built upon an integration of a convolutional neural network, Swin transformer, and the deep operator network architecture, designed to predict anisotropic permeability and upscale predictions to larger spatial domains. First, AnisONet was evaluated with three classes of two-dimensional (2D) porous media, including synthetic circular and elliptical grains and natural sandstone grains from micro-computed tomography images. A lattice Boltzmann model (LBM) was used to calculate directional permeabilities at every 10° angle, producing 19 data points per image of porous media. AnisONet is then trained to predict permeability as a function of rotation angle. AnisONet showed strong predictive capability of directional permeability. Second, we tested our model for five upscaling cases with a large image size in the finite-element method (FEM) for 2D Darcy flow with various permeability tensor construction methods. Overall, upscaled permeability tensors in FEM simulations produce a reasonably good match with LBM results, highlighting the importance of selecting appropriate tensor formation strategies for accurate permeability upscaling. AnisONet, as a directional permeability estimator, could be further developed for more complex geometries, with the potential to develop a foundational ML model for various applications in porous media.

42 ENGINEERING↗
DOE OSTI2024

Slow dynamic elasticity at short times

It has been reported that slow dynamic nonlinear elastic relaxations, widely thought to proceed universally in proportion to the logarithm of time after cessation of mechanical conditioning, actually recover with a smaller slope at early times, with a time of transition that varies with the grain size of the material. Furthermore, this would constitute a heretofore unreported failure of the claimed universality, while suggesting application to nondestructive evaluation and structural health monitoring. Here we present new observations at short times, in the single bead system, in cement paste, mortar, concrete, sandstone and granite. Within the limits imposed by finite duration ring-down such that the effective instant of conditioning cessation is imprecise, and the corresponding ambiguity as to the time at which relaxation begins, we find no reliable sign of such a transition, even in samples of large grain size mortar and concrete similar to those described elsewhere as having clear and late cutoffs.

58 GEOSCIENCES↗
DOE OSTI2025

Wave ripples formed in ancient, ice-free lakes in Gale crater, Mars

Symmetrical wave ripples identified with NASA’s Curiosity rover in ancient lake deposits at Gale crater provide a key paleoclimate constraint for early Mars: At the time of ripple formation, climate conditions must have supported ice-free liquid water on the surface of Mars. These features are the most definitive examples of wave ripples on another planet. The ripples occur in two stratigraphic intervals within the orbitally defined Layered Sulfate Unit: a thin but laterally extensive unit at the base of the Amapari member of the Mirador formation, and a sandstone lens within the Contigo member of the Mirador formation. In both locations, the ripples have an average wavelength of ~4.5 centimeters. Internal laminae and ripple morphology show an architecture common in wave-influenced environments where wind-generated surface gravity waves mobilize bottom sediment in oscillatory flows. Their presence suggests formation in a shallow-water (<2 meters) setting that was open to the atmosphere, which requires atmospheric conditions that allow stable surface water.

Science & Technology - Other Topics↗
DOE OSTI2026

AmeriFlux FLUXNET-1F US-CLF Cole Farm

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-CLF Cole Farm. This is the FLUXNET version of the carbon flux data for the site US-CLF Cole Farm produced by applying the standard ONEFlux (1F) software. Site Description - The Cole Farm catchment (0.65km2) is located ~ 4 km southwest of the Shale Hills site, draining orthogonally to a syncline axis of the Wills Creek Formation, a calcareous shale containing interbedded siltstone, sandstone, shaly limestone, and dolomite. Even though the farm adopted no-till practices in the 1970s, the axial channel of Cole Farm flows over a thick (>2.5 m) package of sediment in the valley floor. Soils range in texture from silty clay at the ridge top to sandy loam in the valley floor. Data was collected and funded by the Critical Zone Observatory Network.

Davis, Kenneth J. [Department of Meteorology, Eart↗
DOE OSTI2024

Geologic Characterization of the South Georgia Rift Basin for Source Proximal CO2 Storage

The project Geologic Characterization of the South Georgia Rift Basin for Source Proximal CO2 Storage is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The South Carolina Research Foundation and partners evaluated the feasibility of CCS in the Jurassic/ Triassic (J / TR) saline formations of the buried Mesozoic South Georgia Rift (SGR) Basin that extends from South Carolina into Georgia. The J / TR sequence, based on preliminary assessment of limited geologic and geophysical data, appears to have both the appropriate areal extent and multiple horizons to permanently and safely store CO2 The presence of several igneous rock layers within the sequence may potentially provide adequate seals to prevent upward CO2 migration into the Coastal Plain aquifer systems. Approximately 81 kilometers of 2-D seismic reflection data were collected by Bay Geophysical, Inc. to explore a portion of the SGR located in southern Georgia. The 81 kilometers were divided into two lines approximately 40.5 kilometers each, with Line 1 intersecting Georgia well GGS 3457. Line 2 intersects Line 1 at the southern portion of Line 1 to maximize the extent of coverage away from GGS-3457 (a deep well drilled in the 1980s for oil and gas exploration). This well had a set of usable logs, including gamma and neutron logs that provided promising results related to CO2 storage. Results showed sandstone with porosity values greater than 10 percent and a thickness of 120 meters. The design of the seismic shot was to extrapolate information away from the well and to better define the extent of the SGR and the necessary reservoir and caprock for a successful CO2 injection. A numerical simulation model of CO2 Injection and migration was developed based on the geology log for the GGS-3457 well. The simulation model was used to investigate the feasibility of injecting 30 million metric tons of CO2 into SGR J / TA sediments and integrity of the diabase layers as seals to prevent CO2 migration.

2-D seismic↗
DOE OSTI2024

Comparative Assessment of U-Net-Based Deep Learning Models for Segmenting Microfractures and Pore Spaces in Digital Rocks

Segmentation of high-resolution X-ray microcomputed tomography (µCT) images is crucial in digital rock physics (DRP), affecting the characterization and analysis of microscale phenomena in the porous media. The complexity of geological structures and nonideal scanning conditions pose significant challenges to conventional image segmentation approaches. Motivated by the recent increasing popularity of deep learning (DL) techniques in image processing, this work undertakes a comparative study of DL models, specifically U-Net and its variants, for segmenting multiple targets with distinguished features in digital rocks, including discrete fracture networks (DFNs), pore spaces, and solid rock. Particularly, DFNs have a smaller volumetric fraction over others, bringing in a substantial challenge of imbalanced segmentation. The primary focus is to evaluate the architecture and feature enhancement strategies of various DL models, including U-Net, attention U-Net, residual U-Net, U-Net++, and residual U-Net++. The models were designed as 2.5D, utilizing a central 2D image and its two adjacent upper and lower 2D images as input to provide a pseudo-3D context. In addition, because the ground truth of segmentation was unknown for real-world digital rocks, we created a benchmark data set following the inverse operations of segmentation. The data synthesis started from the label images (i.e., solid rock, pore spaces, and DFNs), followed by simulating partial volume blurring, adding random background noise, and introducing ring artifacts to mimic real raw X-ray µCT images. The data set, which included various rock types (i.e., sandstone and artificial data), scanning resolution, and magnitudes of noise and artifacts, was divided into training and testing data sets with a 90% and 10% ratio, respectively. Moreover, in addition to the conventional pixel-wise evaluation metrics, the physics-based metric of the lattice-Boltzmann method (LBM) simulated permeability provided more comprehensive assessments. The results demonstrated that the residual connections, nested architectures, and redesigned skip connections contribute to the model performance and give the residual U-Net++ the highest accuracy. The improvements were mainly on the boundaries and small targets, especially the DFNs, which dominate the interconnectivity and therefore affect the permeability greatly. This study also rigorously evaluated the efficiency and generalization of each model, demonstrating that the sophisticated architectures achieved excellent practicability and maintained robust performance on completely unseen data, ensuring their suitability for diverse and challenging DRP applications.

58 GEOSCIENCES↗
DOE OSTI2024

The Forgotten Marshallese of Wotho, Rongelap, and Utirik Atolls

Two of the Laboratory’s major roads are named for Bikini and Enewetak Atolls, where sixty-six nuclear tests were conducted between 1946 and 1958. These signs honor the sacrifices made by the natives of both atolls, who left their ancestral islands to make room for Operations Crossroads, Sandstone, Greenhouse, Ivy, Castle, Redwing, and Hardtack I. Lost in the mist of time, however, are similar sacrifices made by the islanders of Wotho, Rongelap, and Utirik Atolls.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Well-Log Derived Geomechanical Analysis of Microseismicity in the Mt. Simon Saline Aquifers (Illinois Basin - Decatur Project)

The Illinois Basin Decatur Project (IBDP) successfully demonstrated the safe geologic storage of carbon dioxide at a commercial scale. Within the IBDP project three deep wells (injection (CCS1), monitoring (VW1), geophysical (GM1)) were competed and geophysical logs were recorded. During injection and post-injection periods microseismic monitoring was conducted to create a miscoseismic catalog. The correlations between microseimic attributes and geomechanical well logs define major geomechanical drivers of microseismic expression to understand a reservoir response to CO2 injection in geological context. Utilizing standard sonic and density well logs, the dynamic elastic moduli were calculated and employed to correlate with microseismic pseudo-logs. A multi-dimensional Mu-rho and Lambda-rho (MRLR) hyperdimensional plots display of meaningful data and uncovered subtle relationships between elastic properties of sandstones and the seismological attributes of recorded microseismicity.

Myshakin, Evgeniy↗
DOE OSTI2024

Computed Tomography Scanning and Geophysical Measurements of the Integrated Mid-Continent Stacked Carbon Storage Hub Sleepy Hollow Reagan Unit 86A Well

The Computed Tomography (CT) facilities, the Multi-Sensor Core Logger (MSCL), and the Geologic Storage Core Flow laboratory at the National Energy Technology Laboratory (NETL) in Morgantown, West Virginia and Pittsburgh, Pennsylvania were used to characterize a core through Upper Pennsylvanian strata (limestones, mudstones, and sandstones) from the Sleepy Hollow Reagan Unit (SHRU) 86A well in southwest-central Nebraska. The Integrated Mid-Continent Stack Carbon Storage Hub (IMSCS-HUB) core from the vertical well was obtained as part of the Department of Energy’s (DOE) effort to assess the feasibility of stacked storage complexes in Nebraska and Kansas to support a commercial-scale CO2 storage hub. The Sleepy Hollow Field (SHF) is one of three sites within the IMSCS-HUB corridor. Bulk scans of core were obtained from the IMSCS-HUB SHRU 86A well. This report, and the associated scans, include detailed datasets not typically made available to the public. The data sets presented in this report can be accessed from NETL's Energy Data eXchange (EDX) online system.

58 GEOSCIENCES↗
DOE OSTI2024

Geophysical Signatures of Crack Network Coalescence in Rocks at Multiple Length Scales

The main goal of the research project was to identify the geophysical signatures of fracture growth in natural rocks by utilizing novel geophysical techniques. The research objectives were to (a) investigate the potential for geophysical methods to determine when cracks initiate, the types and locations of propagated cracks, and the coalescence of networks of cracks in natural rocks at multiple scales, (b) determine how damage at the microscale evolved into damage at the macroscale and then link the microscopic and macroscopic observations, (c) quantify crack coalescence in rocks under realistic stress conditions using coupled mechanical-geophysical-optical visualization, and (d) identify the precursors in geophysical signals to crack coalescence. The following research thrusts were explored to achieve the research objectives: (1) uniaxial compression testing of rock specimens with and without a set of pre‐existing flaws and (2) triaxial compression testing of natural rock specimens. These thrusts allowed for exploring fracturing in rocks under realistic in situ environments and at multiple scales. This project provided educational opportunities for nine graduate and undergraduate students and resulted in 27 peer-reviewed publications. This first research thrust focused on investigating the micromechanics of fractures in rocks through uniaxial compression testing combined with advanced geophysical and imaging techniques, specifically acoustic emission (AE) monitoring, ultrasonic imaging, and 2-dimensional Digital Image Correlation (2D-DIC). By examining damage processes under time-independent and time-dependent loading conditions, insights into damage localization, crack initiation, and fracturing mechanisms were gained. It was observed that the AE signals and the strain-based measurements directly reflect the state of damage in the rock specimen and could be used to identify the cracking levels, such as the crack initiation (CI) and crack damage (CD), and the mode of deformation. A novel calibration apparatus was developed to enhance the accuracy of AE sensors, allowing for the estimation of key parameters such as magnitude, source dimension, stress drop, and radiated seismic energy associated with the fractures. The findings highlighted significant variations in the temporal evolution of AE source parameters during the primary, secondary, and tertiary stages of creep, identifying tensile cracking as the primary deformation mode. The second research thrust focused on enhancing the understanding of fracturing processes in natural rocks through triaxial compression testing, real-time AE monitoring, and ultrasonic monitoring. We investigated the impact of various factors such as fracture propagation regimes, injection parameters, rock types, and pre-existing conditions on the hydraulic fracture (HF) behavior using scaled true-triaxially loaded specimens of Barre granite and Lyons sandstone. Custom sensor housing facilitated concurrent active and passive monitoring to analyze hydro-mechanical responses and microseismicity associated with different HF scenarios. A coupled investigation of passive microseismicity and active signal attributes permitted a detailed comprehension of the various HF processes (aseismic deformation, fracture initiation and propagation, fluid permeation, and leak-off) and their dependence on the specific rock type. The findings of this research demonstrated the effectiveness of AE monitoring techniques in providing valuable insights into the impact of various factors on the behavior and dynamics of HF processes. The advancements in monitoring techniques, offering a more thorough and precise approach, represent a significant step towards optimizing HF practices and ensuring sustainable resource extraction.

58 GEOSCIENCES↗
DOE OSTI2024

Upscaling Methods Applied to a Fine-Scale Reservoir Model

This study was conducted as part of the Southwest Regional Partnership on Carbon Sequestration (SWP) project to evaluate how upscaling fine-scale simulation models to coarse-scale simulation models impacted the results. The focus was on the Farnsworth Unit (FWU) and its Morrow' B' Sandstone reservoir, specifically the west half of the field. Due to data limitations and the geologic characteristics of the surrounding area, the upscaling was limited to the west half of the FWU rather than a broader basinscale model. The primary aim was to explore how upscaling impacts numerical simulation models, particularly regarding CO 2 -enhanced oil recovery (EOR) and storage capacity predictions. Upscaling was necessary to reduce computational demands when transitioning from high-resolution geological models to coarser grids, as large-scale simulations with finer grids can be computationally prohibitive. This study expands on previous work by the SWP to understand how additional upscaling, applied to already fine-scale numerical simulation models, affects reservoir performance simulations (Ampomah, Balch, & Grigg, 2015). This is key to understanding how loss of resolution can affect coarsescale model results that may be used for large sensitivity analyses, uncertainty quantifications, and training data for machine learning applications.

02 PETROLEUM↗