Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SANDSTONES”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Calculating Trajectories Associated With Solute Transport in a Heterogeneous Medium

Here, we present a trajectory-based technique for calculating solute transport in a porous medium that has several advantages over existing methods. Unlike streamlines, the extended trajectories are influenced by each of the important parameters governing transport, including molecular diffusion and transverse dispersion. The approach is complete and does not require any additional techniques, such as operator splitting or particle tracking, in order to account for the full dispersion tensor. The semianalytic expressions make it clear how the flow field, the concentration distribution, and the dispersion tensor contribute to the velocity field of an injected solute. The equations are valid for an arbitrary porous medium, including those with rapid spatial variations in properties, overcoming limitations faced by previous approaches based upon asymptotic techniques. A test on a layered model with sharp boundaries indicates that the extended trajectories are compatible with the results of a numerical simulator and differ from streamlines. We also describe a new form of the dispersion tensor that incorporates a known asymmetry. The trajectories indicate that the modifications of the dispersion tensor lead to more focused transport within regions of high conductivity. The trajectories are used to define a semianalytic relationship between solute travel times and variations in solute velocities along a path that may be used for tomographic imaging. In an application to the injection of a radioactive tracer into a Berea sandstone core, monitored using micropositron emission tomographic (micro-PET) observations, the sensitivities are used to map the spatial variations of permeability within the core.

42 ENGINEERING↗
DOE OSTI2019

High-Speed Quantification of Pore-Scale Multiphase Flow of Water and Supercritical CO 2 in 2-D Heterogeneous Porous Micromodels: Flow Regimes and Interface Dynamics

The pore-scale flow of CO 2 and water in 2-D heterogeneous porous micromodels over a Ca range of nearly three orders of magnitude was explored experimentally. The porous geometry is a close reprint of real sandstone, and the experiments were performed under reservoir-relevant conditions (i.e., 8 MPa and 21°C), thus ensuring relevance to practical CO 2 operations. High-speed fluorescent microscopy and image processing were employed to achieve temporally and spatially resolved data, providing a unique view of the dynamics underlying this multiphase flow scenario. Under conditions relevant to CO 2 sequestration, final CO 2 saturation was found to decrease and increase logarithmically with Ca within the capillary and viscous-fingering regimes, respectively, with a minimum occurring during regime crossover. Specific interfacial length generally scales linearly with CO 2 saturation, with higher slopes noted at high Ca due to stronger viscous and inertial forces, as supported by direct pore-scale observations. Additionally, statistical analysis of the interfacial movements revealed that pore-scale events are controlled by their intrinsic dynamics at low Ca, but overrun by the bulk flow at high Ca. During postfront flow, while permeability is typically correlated with total CO 2 saturation in the porous domain (regardless of its mobility), the saturation of active CO 2 pathways in the current study correlated very well with permeability. This alternate approach to characterize relative permeability could serve to mitigate hysteresis in relative permeability curves. Taken together, these results provide unique insights that address inconsistent observations in the literature and previously unanswered questions about the underlying flow dynamics of this important multiphase flow scenario.

42 ENGINEERING↗
DOE OSTI2019

Creep Deformation in Vaca Muerta Shale From Nanoindentation to Triaxial Experiments

In this work, we tested samples cored from the Vaca Muerta shale reservoir using nanoindentation (2 min) and triaxial (12 hr) creep experiments in which confining pressure and differential stresses were <40 MPa. In all cases, we observed transient creep wherein strain increased as the logarithm of time. Creep was always compactional, led to increased moduli, was triggered by changes in either hydrostatic or deviatoric stress, and occurred under loads well below the failure stress. Our results are consistent with yield cap models proposed to describe shear-enhanced compaction of sandstones and carbonates, assuming that the yield surface depends on strain rate. We compared our results to earlier studies that observed the transition from transient creep to approximately constant strain rate behavior. If short-term creep can be quantified by a creep modulus, $C$, and long-term creep, after a transition time, $t_c$, by a linear viscosity, $η$, then $η = Ct_c$. Owing to heterogeneity, the local values of Young's and creep moduli, E and C, measured during nanoindentation, varied by several orders of magnitude. Simple averaging of the indentation results overestimated E and C, as compared to their triaxial counterparts. The kinetics of log time creep, which is seen in various materials and loading circumstances, has been represented by a simple conceptual model incorporating the interplay of viscoelastic elements with widely distributed characteristic times. Thus, we argue that log time creep is an emerging phenomenon, independent of the identities of the underlying physical mechanisms.

58 GEOSCIENCES↗
DOE OSTI2019

Stabilizing Effect of High Pore Fluid Pressure on Slip Behaviors of Gouge-bearing Faults

We conducted experiments to investigate the influence of pore fluid pressure on the frictional strength and slip behavior of gouge bearing faults. Saw cut porous sandstone samples with a layer of gouge powders placed between the precut surfaces were deformed in the conventional triaxial loading configuration. A series of velocity-step tests were performed to measure the response of the friction coefficient to variations in sliding velocity. Pore volume changes were monitored during shearing of the gouge. Our results demonstrate that under constant effective pressure, increasing pore pressure stabilizes the frictional slip of faults with all four gouge materials including antigorite, olivine, quartz, and chrysotile. The stabilizing effect is the strongest in antigorite gouge, which shows an evolution of friction parameters from velocity-weakening toward velocity-strengthening behavior with increasing pore pressure. Experiments with controlled pore volume show that the pore volume reduction diminishes under high pore fluid pressures, implying an increasing dilation component at these conditions. The dilatant hardening mechanism can explain the observed strengthening. These results provide a possible explanation to the observed spatial correlation between slow slip events and high pore pressure in many subduction zones.

58 GEOSCIENCES↗
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 OSTI2018

Differential depth distribution of microbial function and putative symbionts through sediment-hosted aquifers in the deep terrestrial subsurface

An enormous diversity of previously unknown bacteria and archaea has been discovered recently, yet their functional capacities and distributions in the terrestrial subsurface remain uncertain. Here, we continually sampled a CO 2 -driven geyser (Colorado Plateau, Utah, USA) over its 5-day eruption cycle to test the hypothesis that stratified, sandstone-hosted aquifers sampled over three phases of the eruption cycle have microbial communities that differ both in membership and function. Genome-resolved metagenomics, single-cell genomics and geochemical analyses confirmed this hypothesis and linked microorganisms to groundwater compositions from different depths. Autotrophic Candidatus “Altiarchaeum sp.” and phylogenetically deep-branching nanoarchaea dominate the deepest groundwater. A nanoarchaeon with limited metabolic capacity is inferred to be a potential symbiont of the Ca. “Altiarchaeum”. Candidate Phyla Radiation bacteria are also present in the deepest groundwater and they are relatively abundant in water from intermediate depths. During the recovery phase of the geyser, microaerophilic Fe- and S-oxidizers have high in situ genome replication rates. Autotrophic Sulfurimonas sustained by aerobic sulfide oxidation and with the capacity for N 2 fixation dominate the shallow aquifer. Overall, 104 different phylum-level lineages are present in water from these subsurface environments, with uncultivated archaea and bacteria partitioned to the deeper subsurface.

59 BASIC BIOLOGICAL SCIENCES↗
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↗