Search NASA⌕ Search

SEARCH · Search NASA

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

CO 2 -Soluble Nonionic Surfactants for Enhanced CO 2 Storage via In Situ Foam Generation

Geologic carbon storage (GCS) is a rapidly evolving technology with the potential to reduce the environmental impact of fos-sil fuels. Saline aquifers, which comprise a sandstone matrix with brine in the pores, make up much of the pore space availa-ble for CO 2 storage in the United States. When CO 2 is injected in saline aquifers, however, capillary fingering occurs and only a small percentage of the pore space is filled with CO 2 . This is due to the low viscosity of CO 2 , which is roughly ten times less viscous than brine. To address this problem, we tested the ability of inexpensive, commercially available nonionic surfactants to be dissolved in the injected CO 2 and increase the apparent viscosity of CO 2 by generating CO 2 -in-water foams in situ.

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

Porous Colloidal Nanoparticles as Injectable Multimodal Contrast Agents for Enhanced Geophysical Sensing

Injecting fluids into underground geologic structures is crucial for the development of long-term strategies for managing captured carbon and facilitating sustainable energy extraction operations. Here, we have previously reported that the injection of metal–organic frameworks (MOFs) into the subsurface can enhance seismic monitoring tools to track fluids and map complex structures, reduce risk, and verify containment in carbon storage reservoirs because of their absorption capacity of low-frequency seismic waves. Here, we demonstrate that water-based Cr/Zn/Zr MOF colloidal suspensions (nanofluids) are multimodal geophysical contrast agents that enhance near-wellbore logging tools. Based on experimental fluid-only measurements, MIL-101(Cr), ZIF-8, and UiO-66 nanofluids have distinct complex conductivity and/or low-field nuclear magnetic resonance (NMR) signatures that are relevant to field-deployed technologies, implying the potential to enhance near-wellbore monitoring of CO 2 injection and associated processes with downhole logging tools. Small- and wide-angle X-ray scattering characterization of ~0.5 wt % MIL-101(Cr) suspensions confirmed phase stability and provided insight into the fractal nature of colloidal nanoparticles. Finally, low-field (2 MHz) NMR measurements of MIL-101(Cr) nanofluid injection into a prototypical Berea sandstone demonstrate how paramagnetic high-surface area MOFs may dominate the relaxation times of hydrogen-bearing fluids in porous geologic matrices, enhancing the mapping of near-surface and near-wellbore transport pathways and advancing sustainable subsurface energy technologies.

NMR↗
DOE OSTI2020

Influence of Clay Wettability Alteration on Relative Permeability

Abstract Understanding the wettability of porous materials is important to model fluid flow in the subsurface. One of the critical factors that influences wetting in real reservoirs is the composition of geologic materials. The wetting properties for clay minerals can have a particularly strong impact on flow and transport. In this work, we analyze the chemical composition of a Mt. Simon sandstone core to resolve the microscopic structure of clay regions and assess how alterations to the local wetting properties influence multiphase transport based on core flooding experiments and relative permeability simulations. We show that whichever fluid has greater affinity toward clay minerals will tend to accumulate within these high surface area regions, leading to dramatic shifts in the relative permeability. This work establishes the essential importance of the mineral composition and associated wetting properties in the modeling of flow and transport in reservoir‐scale systems.

Fan, Ming↗
DOE OSTI2021

Porosity-Permeability Evolution During Simultaneous Mineral Dissolution and Precipitation

Permeability in reactive porous media evolves in complex ways that are hard to predict. Macroscopic empirical equations are often used to estimate permeability evolution but fail to reflect the impact of pore scale reactions on permeability. This work aims to understand the evolution of permeability in systems where mineral dissolution and precipitation occur simultaneously. Pore network models are used for permeability simulation of a Paluxy sandstone using pore and pore-throat size distributions from the analysis of X-ray CT images. Pore and pore-throat radii are increased and decreased to reflect the effects of dissolution and precipitation reactions. Varying spatial distributions of reactions are simulated and the resulting porosity and permeability values compared with commonly used macroscopic equations. It is observed that when dissolution occurs at the inlet and precipitation at the outlet, porosity increases while permeability decreases. Similar results are also observed in the respective opposite scenario where precipitation occurs at the inlet and dissolution at the outlet. When dissolution and precipitation reactions are randomly distributed throughout the network, porosity increases with little change in permeability that is not captured by the empirical equations. When reactions are controlled by pore and pore-throat sizes, such as dissolution occurring in pores and pore-throats of small size and precipitation in pores and pore-throats of larger size, porosity and permeability decrease. When precipitation occurs in pores and pore-throats of smaller size and dissolution in pores and pore-throats of larger size, porosity and permeability increase and is described relatively well by the Verma-Pruess porosity-permeability equation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Type–Curve Approach for Evaluating Aquifer Properties by Interpreting Shallow Strain Measured During Well Tests

Strains occur at shallow depths in response to pressure changes during well tests in an underlying aquifer, and recent developments in instrumentation have made it feasible to measure essentially the full strain tensor. Simulations using poroelastic analyses indicate that shallow normal strains are approximately proportional to the logarithm of time when a well is injecting into or pumping from a deep aquifer or reservoir. The drawdown is also a linear function of log time, as shown by the classic Cooper-Jacob type-curve analysis. The time when the semilog straight line intercepts the zero-strain axis is similar to the time determined from the Cooper-Jacob pressure analysis, and it can be used to estimate hydraulic diffusivity, suggesting that horizontal strain data can be used directly to estimate aquifer properties. This approach was validated using measurements from shallow (30-m deep) borehole strainmeters during an injection test at a 530-m-deep sandstone aquifer/reservoir in Oklahoma. The results show intercept times for the shallow normal strain data are essentially the same as for deep pressure data from an equivalent radial distance. The slopes of the semilog plots of the pressure and the strain increase at the same time, suggesting that they both respond to a lateral aquifer boundary. Significantly, though, strain was measured at shallow depths while the pressure data were measured at 530-m depth. This suggests that strain data from shallow depths could be an effective way to improve the characterization of an underlying aquifer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Distributed Fiber Optics Measurements of Rock Deformation and Failure in Triaxial Tests

The implementation of distributed strain measurement methods in triaxial and uniaxial tests have demonstrated the development of strain localization, even at early stages. This implies that single point measurement methods are location dependent. The use of distributed methods is required not only to improve the interpreted constitutive parameters obtained from triaxial tests, but also to understand the implications of strain localization in the failure process. Here, this work uses optical fibers in triaxial tests. The developed distributed measurement method was implemented on granite, gabbro, and sandstone samples and tested under different confining pressures, reaching 200 MPa on the granite sample. Using a temporal resolution of 0.25 Hz and 5 mm of spatial resolution, the strain evolution at over 300 locations at the sample surface was measured during testing. When compared to point sensing methods, the use of optical fiber greatly increases the number of measurements at the surface of the sample, and their interpolation provides the entire deformation of the sample surface. In all the tests performed, strain localizations were revealed before failure. A three-dimensional interpretation of a test, combining an optical scan of the sample and the distributed measurements, show good correlation between the fractures and the strain localization.

58 GEOSCIENCES↗
DOE OSTI2023

Using the Shallow Strain Tensor to Characterize Deep Geologic Reservoirs

Abstract Storing and recovering water, carbon, and heat from geologic reservoirs is central to managing resources in a changing climate. We tested the hypothesis that the strain tensor caused by injecting or producing fluids can be measured at shallow depths and interpreted to advance understanding of underlying deep aquifers or reservoirs. Geodetic‐grade strainmeters were deployed at 30 m depth overlying the Bartlesville Formation, a 500‐m‐deep sandstone near Tulsa, OK. The strainmeters are 220 m east of injection well 9A completed in a permeable lens at the base of the Bartlesville Formation. Water was injected into well 9A at approximately 1.0 L/s during four tests that ranged in duration from a few hours to a few weeks. The horizontal strain increased (tension) and the circumferential strain was a few times larger than the radial strain. The vertical strain decreased (compression) during injection. Strain rates were approximately 100 nε/day during the first few hours, but the rates decreased and were approximately 10 nε/day during most of the tests. Four independent methods of poroelastic simulation and inversion predict reservoir properties and geometries that are similar to each other and consistent with independent information about the reservoir. All strain interpretations predict that a boundary to the permeable lens occurs beneath the vicinity of the strainmeters, which is consistent with core data from the site. The boundary of the permeable lens is located by matching the vertical, radial and circumferential strains, which demonstrates the value of measuring the strain tensor.

Murdoch, Lawrence C.↗
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 OSTI2022

Assessment of chemo-mechanical impacts of CO 2 sequestration on the caprock formation in Farnsworth oil field, Texas

This study evaluates the chemo-mechanical influence of injected CO 2 on the Morrow B sandstone reservoir and the upper Morrow shale caprock utilizing data from the inverted 5-spot pattern centered on Well 13-10A within the Farnsworth unit (FWU). This study also seeks to evaluate the integrity of the caprock and the long-term CO 2 storage capability of the FWU. The inverted 5-spot pattern was extracted from the field-scale model and tuned with the available field observed data before the modeling work. Two coupled numerical simulation models were utilized to continue the study. First, a coupled hydro-geochemical model was constructed to simulate the dissolution and precipitation of formation minerals by modeling three intra-aqueous and six mineral reactions. In addition, a coupled hydro-geomechanical model was constructed and employed to study the effects of stress changes on the caprock’s porosity, permeability, and ground displacement. The Mohr–Coulomb circle and failure envelope were used to determine caprock failure. In this work, the CO 2 -WAG injection is followed by the historical field-observed strategy. During the forecasting period, a Water Alternating Gas (WAG) injection ratio of 1:3 was utilized with a baseline bottom-hole pressure constraint of 5500 psi for 20 years. A post-injection period of 1000 years was simulated to monitor the CO 2 plume and its effects on the CO 2 storage reservoir and caprock integrity. The simulation results indicated that the impacts of the geochemical reactions on the porosity of the caprock were insignificant as it experienced a decrease of about 0.0003% at the end of the 1000-year post-injection monitoring. On the other hand, the maximum stress-induced porosity change was about a 1.4% increase, resulting in about 4% in permeability change. It was estimated that about 3.3% of the sequestered CO 2 in the formation interacted with the caprock. Despite these petrophysical property alterations and CO 2 interactions in the caprock, the caprock still maintained its elastic properties and was determined to be far from its failure.

04 OIL SHALES AND TAR SANDS↗
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 OSTI2020

A Study of Temperature Distribution and Thermal Stresses in a Hot Rock Due to Rapid Cooling

Thermal stresses may be induced in a hot dry rock when a cold fluid is injected in the well. To study this problem, we look at the thermoelastic response of a hot rock that is suddenly cooled. The cooling is assumed to be either at a constant temperature or at a constant heat flux per unit depth. Our approach is to nondimensionalize the equations and perform a parametric study and look at the temperature distribution and the induced-thermal stresses. The results indicate that depending on the extent of cooling and the cooling time, thermal stresses can be induced. Numerical simulations on sandstone, with an initial uniform temperature of 473 K, are also carried out. The results show that if the cooling is due to the surface temperature maintained at 463 K (10 °C lower than the initial temperature of the hot rock), thermal stresses that are larger than the rock tensile strength could be induced. When the cooling is due to a constant surface heat flux, this temperature can be reached after about 777 days of cooling with a minimum value of a heat flux of -20 W/m.

Tran, Phuoc↗