SEARCH · Engineering Papers

Results for “rock”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Influence of interfacial responses of Berea Sandstone in low-salinity waterflooding environments

Wetting conditions and fluid-fluid interfacial properties are essential to understanding oil recovery in the context of low-salinity and smart waterflooding. Many have deemed wettability alteration to be the dominant interaction, while some highlight the importance of fluid-fluid interactions. In some rock-fluid systems, benefits from the use of additives can arise to favor hydrocarbon production. Understanding the balance between the fluid-fluid and rock-fluid interactions when utilizing additives is crucial to evaluate benefits in improved efficiency. Ultimately, this understanding will shed light on recovery mechanisms. Regarding oil recovery, previous studies indicate that the injection of a naphthenic acid (NA) blend dissolved in an injection brine improves oil recovery in Berea Sandstone, particularly under low-salinity conditions. Interfacial rheology results for the brine-oil interface correlated well with the recovery response. This observation was speculated to serve as a potential improved recovery mechanism through oil-trapping suppression. However, the said studies did not evaluate potentially concurrent rock wettability alteration during recovery or the initial wetting condition. The purpose of this work is to determine the production behavior in systems with NAs, namely NA-10 (3-Cyclohexanepropinoic acid), NA-14 (Cyclopentane carboxylic acid), NA-18 (Cyclohexane carboxylic acid), and NA-mix (Naphthenic acid mixture), at different initial wetting conditions under spontaneous imbibition. In addition, we tested the NAs ability to alter wettability in Berea Sandstone and determine each “additive’s” capacity to improve recovery. Here, time-domain nuclear magnetic resonance (TD-NMR), contact angle, and Amott cell tests, i.e. spontaneous imbibition, were used to analyze the influence of the selected naphthenic acids on oil displacement under different initial wetting conditions. Our results indicate that individual naphthenic acids indeed have potential as additives in low-salinity waterflooding for oil-wet systems while negatively impacting recovery in water-wet systems. However, results suggest that the evolution of the interfacial properties and wettability conditions are equally important on the recovery behavior.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Unravelling the origins of shale nanoporosity using small-angle neutron scattering (SANS)

Hydrocarbon production from tight rocks is constrained by slow diffusion within the shale matrix, limited by small pore sizes and low permeability. The nanopore proportion and size distribution significantly influence matrix permeability, a key property for optimizing hydrocarbon recovery and supporting hydrogen production while minimizing environmental impacts. Small-angle neutron scattering (SANS) has been an important tool for exploring the characteristics and structure of shale nanopores. This study used SANS to analyze nanoporosity and pore size distribution (<100 nm) in tight rocks with varying compositions to determine the influence of rock heterogeneity on SANS measurements. Results showed that nanoporosity correlates with clay content, with the highest clay-rich shale (52.48 wt% clay) exhibiting 8.8 % nanoporosity. SANS also revealed more nanopores than traditional nitrogen adsorption measurements, affirming its ability to reflect bulk mineralogy and upholding the relevance of experimental findings using this technique to optimize field operational approaches.

03 NATURAL GAS↗

DOE OSTI2020

Deposition of sulfate aerosols with positive Δ 33 S in the Neoarchean

Anomalous sulfur isotope compositions present in Archean rocks have been intensely scrutinized over the last 20 years because they record key aspects of Earth's atmospheric composition prior to the appearance of free molecular oxygen ca. 2.3 billion years ago. These isotopic compositions can be described as mass anomalous fractionations (MAF) and are produced in the atmosphere as UV light interacts with SO 2 molecules. Most interpretations suggest that atmospheric processes generate a reduced S-phase with a positive ( 33 S-enriched) MAF signature, as measured in pyrites, and an oxidized S-phase with a negative anomaly, as measured in bedded barite deposits. However, recent data for carbonate-associated sulfate (CAS) — a direct proxy for the isotopic composition of sulfur from seawater sulfate — in Neoarchean rocks showed no such negative values, but rather the opposite. To understand if the positive MAF anomalies we measured in Neoarchean CAS reflect secondary processes (diagenetic, metamorphic, handling) instead of original signals of Archean seawater sulfate, we collected additional sample suites with various degrees of preservation and metamorphic alteration across the Campbellrand-Malmani platform in South Africa. Results illustrate that within this comprehensive suite, less-altered samples all contain positive MAF values while secondary processes tend to either remove CAS from the sample and/or decrease the 33 S-enrichment. This positive MAF signal in sulfate is therefore reasonably interpreted as a primary depositional origin, and implies that the assumption that sulfate always carries a negative MAF anomaly throughout the Archean rock record needs to be reconsidered. Our CAS observations suggest that future experiments and calculations should also consider atmospheric and/or sulfur cycling processes that can produce oxidized sulfur with a positive MAF signature.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Influence of microstructure on replacement and porosity generation during experimental dolomitization of limestones

Replacement reactions commonly alter the multiscale pore structures of rocks during fluid-rock interactions. Analysis of these processes in various model fluid-rock systems during controlled laboratory experiments provides insights into the origins of microstructures found in natural materials. This study focused on understanding the effects of initial starting material permeability and resultant differences in transport pathways on porosity and mineralogical changes during limestone dolomitization.A series of replacement experiments (32–317 days in duration) have been conducted in which 1.59 cm (5/8 in.) diameter cores of two different limestones were reacted with saturated MgCl2 solutions at 200 °C. The Texas Cream (Austin Chalk) is a high-porosity, high-permeability limestone, whereas both the porosity and permeability of the Carthage Marble (Burlington Limestone) are relatively low. Altered limestones were imaged using scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM-EDX), Time-of-Flight Secondary Ion Mass Spectrometry (ToF–SIMS) and electron microprobe analysis (EMPA). A representative grain boundary of the low-porosity limestone was targeted for a focused ion beam (FIB) lift-out and characterized using transmission electron microscopy (TEM). These results were coupled with analyses of radial changes in the porosity distribution of the core derived from X-ray and neutron small- and ultra-small angle scattering ((U)SANS/(U)SAXS).The high-porosity/permeability limestone showed a four times faster bulk replacement rate than the lower-porosity/permeability material, and a different mechanism of porosity development. For the low-porosity limestone, a two-stage replacement occurred, with the reacted region of the core consisting of an inner rim in which the limestone was replaced by two calcite-dolomite solid solutions, and an outer rim in which the dolomite was replaced by magnesite. Elongated pores formed along grain boundaries at the initial limestone/dolomite reaction interface, and additional nanometer-scale porosity was formed at the secondary magnesite replacement rim. Grain boundaries were identified as preferential pathways for transport leading to dolomitization and a grain boundary diffusion rate was calculated based on microstructural characterization. In contrast, replacement in the high-porosity limestone was accompanied by porosity generation through replacement of individual grains by dolomite throughout the sample and, in longer runs, magnesite in outer parts of the core. These observations emphasize that both the mechanisms of the replacement reaction and the microstructure and chemistry of the replaced product are contingent on the initial structure of the starting material.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Oxidation rates and redox stabilization of ferrous iron in trioctahedral smectites

Iron(II)-bearing trioctahedral smectites (saponites) form during anoxic alteration of basaltic rock. They are predicted to have been widespread on the early Earth and are observed in the oceanic subsurface today. Smectite structures, including the occupancy of sites in the octahedral sheet, affect iron redox behavior but the rates and products of trioctahedral smectite oxidation have been largely unexplored to date. In this study we synthesized two Fe(II)-bearing trioctahedral smectites, one moderate (22 wt% Fe) and one high (27 wt% Fe) in iron content. We then examined the rate, extent, and products of their oxidation by dissolved oxygen, nitrite, and hydrogen peroxide. Dissolved oxygen caused partial oxidation of Fe(II) in the smectites with 14 to 43% of Fe(II) unoxidized after 20 to 30 days of exposure. The rate and extent of oxidation correlated with the dissolved oxygen concentration and the Fe(II) content of the clay. The incomplete oxidation in these experiments is consistent with the mixed-valent trioctahedral smectites observed in oxidized natural samples but contrasts with the complete reoxidation by oxygen shown by chemically- or microbially-reduced dioctahedral smectites. Oxidation of structural Fe(II) by 5 mmol L -1 nitrite was negligible for the moderate-iron smectite and yielded only ~17% oxidation after 54 days of reaction for the high-iron smectite. Hydrogen peroxide caused rapid and near-complete oxidation of both clays. Powder X-ray diffraction, variable-temperature Mössbauer spectroscopy, and extended X-ray absorption fine structure spectroscopy together detected no crystalline or short-range-ordered secondary phases and show that oxidized iron remained in the trioctahedral smectite structure. The recalcitrant Fe(II) pool in oxidized trioctahedral smectites exists in less distorted sites than Fe(II) in the initial clays. Its unreactive nature at prolonged reaction times indicates an elevated redox potential generated by the local coordination environment. Slower oxidation rates create a larger recalcitrant Fe(II) pool, suggesting kinetic competition between oxidation and a process involved in redox stabilization, such as electron exchange between octahedral iron sites or deprotonation of hydroxyl groups in the structure. The resistance to complete oxidation of trioctahedral ferrous smectites and their full retention of iron demonstrates that transitions from anoxic to oxic conditions generate mixed-valence smectites rather than a mixture of new phases. Identifying the diagenetic products of mixed-valent trioctahedral smectites may provide an indicator in the rock record of past redox cycling. In conclusion, substantial portions of structural Fe(II) in trioctahedral smectites display slow abiotic oxidation kinetics and represent potential electron donors for both microaerophilic iron oxidizing and nitrate-reducing, iron-oxidizing microorganisms in altered mafic rocks and related settings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Comparative Dissolution of Iron-Bearing Minerals by Catecholate and Hydroxamate Siderophores under Oxic and Anoxic Conditions

Siderophores play a crucial role in biological iron (Fe) acquisition and mobilization by promoting Fe mineral and rock weathering. While the effects of the hydroxamate siderophore desferrioxamine B (DFOB) have been extensively studied, the role of catecholates in the dissolution of Fe mineral and rock under varying redox conditions remain limited. Moreover, despite Fe being one of the most redox-active metals, the siderophore-mediated redox transformations of solid-phase Fe and their impact on mineral dissolution are not well understood. Herein, this study systematically investigated dissolution behavior of both Fe(II) and Fe(III)-bearing minerals and rocks (Fe(II)- bearing basalt and olivine, and Fe(III)- bearing nontronite and goethite), by two structurally distinct biological siderophores, catecholate protochelin and hydroxamate DFOB, under both oxic and anoxic conditions. Batch experiments quantified Fe and transition metals (Mn, Ni, Zn, Mo) released from the tested minerals in the presence of siderophores. Throughout the dissolution experiments, siderophore concentrations and Fe-siderophore complexation in solutions were measured using UV-vis spectrometry. Fe redox changes under oxic and anoxic conditions by siderophores were determined, and mineral surface alterations following siderophore treatments were characterized through scanning electron microscopy (SEM) and time-of-flight secondary ion mass spectroscopy (ToF-SIMS). Results revealed distinct interplays of dissolution mechanisms, including ligand-exchange promoted dissolution and reductive dissolution, along with Fe(III) reduction and Fe(II) oxidation, contributing to the Fe(II)- and Fe(III)-rich mineral weathering under varying redox conditions. Under oxic conditions, both protochelin and DFOB enhanced Fe release from Fe(II)-rich basalt and olivine more effectively than from Fe(III)-rich minerals. DFOB induced greater Fe(II)-mineral dissolution than protochelin. This difference was attributed to a higher level of Fe(II) oxidation by DFOB in contrast to protochelin, as well as the protochelin oxidation with the loss of Fe binding groups (catechols). Under anoxic conditions, both siderophores significantly reduced Fe(III) in nontronite and goethite, with protochelin demonstrating substantially stronger Fe(III) reduction capacity than DFOB. However, Fe(III) reduction negatively impacted Fe release from Fe(III)-rich minerals. Instead, Fe release from Fe(II)-rich minerals was enhanced under anoxic conditions due to the absence of Fe(II) oxidation and labile Fe(II). Variations in siderophore adsorption also contributed to dissolution differences. Beyond iron, the release of transition metals from tested minerals was influenced by elemental contents, mineral compositions, and siderophore types, revealing distinct metal- and siderophore-dependent patterns. This systematic investigation highlights the roles of siderophores in Fe redox cycling and metal mobilization, enhancing our understanding of different siderophore behaviors in siderophore-mediated microbial metal acquisition within redox-dynamic environments, with implications for bioleaching industries and applications in agriculture and climate change mitigation.

Guo, Dongyi↗

DOE OSTI2024

Modeling gas migration through clay-based buffer material using coupled multiphase fluid flow and geomechanics with stress-dependent gas permeability

A model for gas migration through clay-based buffer material is developed for modeling gas generation and migration associated with deep geologic nuclear waste disposal. The model is based on a multiphase fluid flow and geomechanics simulator that is adapted to consider enhanced gas flow when gas pressure is high enough to approach the confining stress magnitude. A key feature in the model is a direct coupling between gas permeability and stress, through a non-linear stress-dependent permeability function. The model was first tested and calibrated by modelling two different laboratory gas migration tests on Wyoming (MX-80) bentonite samples. The calibrated model was then applied to model gas migration through a bentonite buffer of a large-scale gas injection test (Lasgit) conducted at the Äspö Hard Rock Laboratory in Sweden. Observed preferential gas migration along interfaces (between compacted blocks and along the canister surface) required explicit representation of such interfaces in the model. The model with the stress-dependent gas permeability accurately captured observed experimental responses in terms of gas breakthrough time, peak gas pressure, and cumulative gas flow rates. The calibrated model was finally applied to simulate migration of hydrogen gas generated within a breached nuclear waste canister over 10,000 years, involving migration of much larger gas volumes. For the considered gas generation rate and host rock properties, the generated gas could migrate through the bentonite buffer and released into the surrounding host rock at a maximum gas pressure somewhat higher than the initial total stress, though a significant amount of hydrogen remained within the buffer. This modelling sets the stage for further detailed analysis of the impact of hydrogen gas generation on the long-term performance of nuclear waste repositories.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Prediction of the thermal-hydraulic-mechanical response of a geological repository at large scale and sensitivity analyses

To better understand the Thermal-Hydraulic-Mechanical (THM) response of the rock caused by the high-level radioactive waste (HLW) released heat in the Callovo-Oxfordian formation (COx) in the near-field and far-field areas, a series of coupled THM modelling is performed. This study presents a case study based on the French geological repository for HLW (Cigéo project) to better assess modelling of deep geological repositories (DGR) within the DECOVALEX-2019 framework. In this study, the proposed coupled THM model, implemented in the platform of COMSOL, is validated by comparing the modelled THM results for a point heat source in an infinite rock mass with an analytical solution. It is also validated by comparing its thermal component for the Base Case with accurate results calculated using another numerical method. To shorten the calculation time, the DGR is modelled using six cuboid blocks representing 162 HLW cells and only six central HLW cells of interest modelled in detail. This simplification is validated. This model can easily be used to model not only the near-field THM response of the rock but also the far-field THM response. The influence of the boundary conditions applied on the gallery wall, HLW cell walls, and external surfaces of the model and the influence of different vertical dimensions of the geometry are studied. The sensitivity analyses of the THM parameters of COx on the THM response at different locations are performed. Finally, the influence of using a 2-dimensional (2D) model to represent a 3-dimensional (3D) repository is also investigated.

58 GEOSCIENCES↗

DOE OSTI2024

The role of annealing and grain boundary controls on the mechanical properties of limestones and marbles

Chemical and mechanical processes are coupled in many geological and geochemical environments. Reactive processes anneal defects and restructure grain boundaries, modifying their elastic properties, levels of internal friction, wave propagation rates and fracture behaviors. The nature of these changes is, however, contingent on the initial state of the rock. Here, in this study, impulse excitation was used to measure changes in mechanical properties as a function of dry and steam heating time at 300 °C for three carbonate rocks: Carrara marble, Carthage marble (Burlington Limestone), and Texas Cream limestone (Edwards limestone), with initial porosities of about 1 %, 3 %, and 27 %, respectively. Frequency-dependent phenomena along with mineral recrystallization were observed. This was coupled with small-angle X-ray and neutron scattering analysis to determine the relationship between changes in bulk mechanical and microstructural properties. An observed decrease in both the Young's and shear moduli in the experimental limestones as a function of heating time reflects and quantifies a reduction in the stiffness of the rock due to annealing. The internal friction of the samples first increases then decreases with reaction time, reflecting defect annealing, but this was balanced against grain boundary dissolution apparently driven by condensation of steam in the confined grain-boundary environment. This suggests an increase in the boiling point under confinement leading to dissolution, increased porosity and widening of the grain boundaries. In addition, comparison of small-angle X-ray and neutron scattering results suggests that it is inappropriate to assume that pores are empty for quantitative analysis of typical small-angle scattering samples.

58 GEOSCIENCES↗

DOE OSTI2022

A gradient-based deep neural network model for simulating multiphase flow in porous media

We report simulation of multiphase flow in porous media is crucial for the effective management of subsurface energy and environment-related activities. The numerical simulators used for modeling such processes rely on spatial and temporal discretization of the governing mass and energy balance partial-differential equations (PDEs) into algebraic systems via finite-difference/volume/element methods. These simulators usually require dedicated software development and maintenance, and suffer low efficiency from a runtime and memory standpoint for problems with multi-scale heterogeneity, coupled-physics processes or fluids with complex phase behavior. Therefore, developing cost-effective, data-driven models can become a practical choice, and in this work, we choose deep learning approaches as they can handle high dimensional data and accurately predict state variables with strong nonlinearity. In this paper, we describe a gradient-based deep neural network (GDNN) constrained by the physics related to multiphase flow in porous media. We tackle the nonlinearity of flow in porous media induced by rock heterogeneity, fluid properties, and fluid-rock interactions by decomposing the nonlinear PDEs into a dictionary of elementary differential operators. We use a combination of operators to handle rock spatial heterogeneity and fluid flow by advection. Since the augmented differential operators are inherently related to the physics of fluid flow, we treat them as first principles prior knowledge to regularize the GDNN training. We use the example of pressure management at geologic CO 2 storage sites, where CO 2 is injected in saline aquifers and brine is produced, and apply GDNN to construct a predictive model that is trained with physics-based simulation data and emulates the physics process. We demonstrate that GDNN can effectively predict the nonlinear patterns of subsurface responses, including the temporal and spatial evolution of the pressure and saturation plumes. We also successfully extend the GDNN to convolutional neural network (CNN), namely gradient-based CNN (GCNN), and validate its capability to improve the prediction accuracy. GDNN has great potential to tackle challenging problems that are governed by highly nonlinear physics and enable the development of data-driven models with higher fidelity.

42 ENGINEERING↗

DOE OSTI2022

Gas diffusion through variably-water-saturated zeolitic tuff: Implications for transport following a subsurface nuclear event

Noble gas transport through geologic media has important applications in the characterization of underground nuclear explosions (UNEs). Without accurate transport models, it is nearly impossible to distinguish between xenon signatures originating from civilian nuclear facilities and UNEs. Understanding xenon transport time through the earth is a key parameter for interpreting measured xenon isotopic ratios. One of the most challenging aspects of modeling gas transport time is accounting for the effect of variable water saturation of geological media. In this study, we utilize bench-scale laboratory experiments to characterize the diffusion of krypton, xenon, and sulfur hexafluoride (SF6) through intact zeolitic tuff under different saturations. Here, we demonstrate that the water in rock cores with low partial saturation dramatically affects xenon transport time compared to that of krypton and SF6 by blocking sites in zeolitic tuff that preferentially adsorb xenon. This leads to breakthrough trends that are strongly influenced by the degree of the rock saturation. Xenon is especially susceptible to this phenomenon, a finding that is crucial to incorporate in subsurface gas transport models used for nuclear event identification. We also find that the breakthrough of SF6 diverges significantly from that of noble gases within our system. When developing field scale models, it is important to understand how the behavior of xenon deviates from chemical tracers used in the field, such as SF6 (Carrigan et al., 1996). These new insights demonstrate the critical need to consider the interplay between rock saturation and fission product sorption during transport modeling, and the importance of evaluating specific interactions between geomedia and gases of interest, which may differ from geomedia interactions with chemical tracers.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

DOE OSTI2020

Thermo-hydrologic processes in maar eruptions: The role of vapor transport and condensation

Phreatomagmatic eruptions occur when magma interacts directly with water or slurries, commonly from groundwater in the shallow crust, resulting in rapid conversion of thermal energy to mechanical energy in a set of processes known as molten fuel-coolant interaction (MFCI). The efficiency of MFCI is controlled by properties of the magma and coolant, contact surface geometry, and system conditions, such as pressure and temperature; phreatomagmatic behavior can range from passive thermal granulation to violent thermohydraulic explosions (Zimanowski et al., 2015). Maar-diatreme eruptions, characterized by repetitive phreatomagmatic explosions, can form in a wide range of near-surface environments, from soft sediment substrate (e.g., Tecuitlapa, Mexico, Ort and Carrasco-Núñez, 2009; Hopi Buttes, AZ, Lefebvre et al., 2013) to hard, fractured country rock (e.g., West Eifel volcanic field, Germany and Massif Central, France, Lorenz, 2003). Country-rock structure, including faults and joints, rock type, and shallow crustal hydrologic properties, can influence an eruption’s behavior by controlling how much and at what rate water can be supplied to the magma (hydraulic flow rate), and determining the sediment: water ratio in the case of soft-sediment substrate eruptions. In a “soft” substrate environment, volcanic tremor can lead to liquefaction of saturated sediment, producing a slurry that may then interact with magma in unique ways (White, 1996; Auer, et al., 2006).

58 GEOSCIENCES↗

DOE OSTI2021

Factors controlling pore network development of thermally mature Early Palaeozoic mudstones from the Baltic Basin (N Poland)

Understanding the formation of pore space, especially in low porosity shales (as source rocks and as unconventional resources), is critical to the oil and gas industry, since pores control the space available for hydrocarbon and participate in hydrocarbon transport. We examined 87 Ordovician and Silurian mudstone samples collected from four wells located in the Pomeranian part of the Baltic Basin (northern Poland), one of the primary Polish targets for hydrocarbon exploration. These samples represent the Pelplin, the Paslek, the Jantar, the Prabuty, and the Sasino Formations, which still requires more detailed porosity studies. Our study aimed to identify factors controlling porosity development, by applying bulk techniques (organic petrology and TOC analyses, quantitative mineralogy, and porosimetry) as well as nano- to microscale techniques (thin section petrography, electron microscopy). The studied samples are mainly argillaceous mudstones. The results of porosimetry measurements, combined with image analysis, indicate that the pores of all studied rocks are dominated by micropores (pores < 2 nm in diameter), mesopores (2-50 nm in diameter) and small macropores. The SEM images showed three main pore types: a) voids related to clay mineral aggregates, b) pores inside organic matter particles, and c) pores between other mineral grains. In the Jantar and Sasino mudstones, the organic matter content and its thermal maturity control porosity. The occurrence of solid bitumen in the rocks from these formations reduces samples’ mesoporosity because of the pore-clogging effect. In contrast, in the Paslek and Prabuty Formations, there is low organic matter content and specific surface area and the volume of mesopores increase with clay minerals content. In the Pelplin mudstones, there are no prevailing factors controlling porosity. As a result, we suggest that a combination of SEM image analysis and dual liquid porosity (DLP) measurements is a powerful method to assess porosity available for petroleum flow in mudstones.

03 NATURAL GAS↗

DOE OSTI2020

Macro- and micro-compression testing of shales

Here, macro- and micro-geomechanical properties of shales were investigated at millimeter- and micrometer-scale using an example of the Woodford Shale. Geomechanical properties, such as uniaxial (unconfined) compressive strength and Young's modulus, were quantified at the millimeter-scale and compared with results from micro-compression testing experiments of FIB-SEM-nanofabricated micro-pillars. Size-scale and compositional/structural heterogeneity effect on elasto-plastic deformation and failure behavior of shales were investigated. Also, relationship between elemental/mineral composition and unconfined compressive strength and Young's modulus, at a micrometer-scale, was discussed. It was shown that geomechanical properties of shales are scale-dependent and are strongly affected by compositional/structural anisotropy. In particular, micro-compression testing experiments showed non-uniformly distributed compressive strength and Young's modulus within the investigated Woodford Shale rock samples. It was demonstrated, that geomechanical properties of shales, investigated at the micrometer-scale, tend to be affected by different minerals and/or pores, while the same properties, investigated at the millimeter-scale, are governed by different micro-facies (and/or micro-fractures) present within the rock. This shows that geomechanical properties of shales can be dramatically different depending on the scale of investigation and compositional/structural heterogeneity of these rocks, and therefore are not easily transferable across the scales.

58 GEOSCIENCES↗

DOE OSTI2020

Estimating diffusion coefficients of shale oil, gas, and condensate with nano-confinement effect

Molecular diffusion plays an important role in oil and gas migration in tight shale formations. However, there are insufficient reference data in the literature to specify the diffusion coefficients within porous media. This study aims at calculating diffusion coefficients of shale gas, shale condensate, and shale oil at reservoir conditions with CO 2 injection for EOR/EGR. The large nano-confinement effects including large gas-oil capillary pressure and critical property shifts could alter the phase behaviors. This paper estimates the diffusivities of shale fluids in nanometer-scale shale rock from two perspectives: 1) examining the shift of diffusivity caused by nanopore confinement effects from phase change (phase composition and fluid property) perspective, and 2) calculating the effective diffusion coefficient in porous media by incorporating rock intrinsic properties (porosity and tortuosity factor). The tortuosity is obtained by using tortuosity-porosity relations as well as the measured tortuosity of shale from 3D imaging techniques. The results indicated that nano-confinement effects could affect the diffusion coefficient through altering the phase properties, such as phase compositions and densities. Compared to bulk phase diffusivity, the effective diffusion coefficient in porous shale rock is reduced by 10 2 to 10 4 times as porosity decreases from 0.1 to 0.03.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2023

Environmental impact of solution pH on the formation and migration of iron colloids in deep subsurface energy systems

Deep subsurface stimulation processes often promote fluid-rock interactions that can lead to the formation of small colloidal particles that are suspected to migrate through the rock matrix, partially or fully clog pores and microfractures, and promote the mobilization of contaminants. Thus, the goal of this work is to understand the geochemical changes of the host rock in response to reservoir stimulation that promote the formation and migration of colloids. Two different carbonate-rich shales were exposed to different solution pHs (pH = 2 and 7). In this study, iron and other mineral transformations at the shale-fluid interface were first characterized by synchrotron-based XRF mapping. Then, colloids that were able to migrate from the shale into the bulk fluid were characterized by synchrotron-based extended X-ray absorption structure (EXAFS), scanning electron microscopy (SEM), and single-particle inductively coupled plasma time-of-flight mass spectrometry (sp-icpTOF-MS). When exposed to the pH = 2 solution, extensive mineral dissolution and secondary precipitation was observed; iron-(oxyhydr)oxide colloids colocated with silicates were observed by SEM at the fluid-shale interfaces, and the mobilization of chromium and nickel with these iron colloids into the bulk fluid was detected by sp-icpTOF-MS. Iron EXAFS spectra of the solution at the shale-fluid interface suggests the rapid (within minutes) formation of ferrihydrite-like nanoparticles. Thus, we demonstrate that the pH neutralization promotes the mobilization of existing silicate minerals and the rapid formation of new iron colloids. These Fe colloids have the potential to migrate through the shale matrix and mobilize other heavy metals (such as Cr and Ni, in this study) and impacting groundwater quality, as well produced waters from these hydraulic fracturing operations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Geology and elevation shape bacterial assembly in Antarctic endolithic communities

Ice free areas of continental Antarctica are among the coldest and driest environments on Earth, and yet, they support surprisingly diverse and highly adapted microbial communities. Endolithic growth is one of the key adaptations to such extreme environments and often represents the dominant life-form. Despite growing scientific interest, little is known of the mechanisms that influence the assembly of endolithic microbiomes across these harsh environments. Here, we used metagenomics to examine the diversity and assembly of endolithic bacterial communities across Antarctica within different rock types and over a large elevation range. While granite supported richer and more heterogeneous communities than sandstone, elevation had no apparent effect on taxonomic richness, regardless of rock type. Conversely, elevation was clearly associated with turnover in community composition, with the deterministic process of variable selection driving microbial assembly along the elevation gradient. The turnover associated with elevation was modulated by geology, whereby for a given elevation difference, turnover was consistently larger between communities inhabiting different rock types. Overall, selection imposed by elevation and geology appeared stronger than turnover related to other spatially-structured environmental drivers. Our findings indicate that at the cold-arid limit of life on Earth, geology and elevation are key determinants of endolithic bacterial heterogeneity. This also suggests that warming temperatures may threaten the persistence of such extreme-adapted organisms.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Intensity measure adequacy assessment for nonlinear site response using Information Theory

Nonlinear site response modeling is a crucial aspect of Probabilistic Seismic Hazard Analysis. Site amplification models routinely rely on a rock intensity measure to characterize the strength of the bedrock motion. However, the adequacy of such intensity measures towards predicting amplifications across the oscillator period range has not been investigated in the literature. This paper analyzes the adequacy of rock intensity measures using state of the art criterion established in Performance-Based Earthquake Engineering and techniques from Information Theory. The efficiency and the sufficiency of several rock intensity measure are assessed. It was found that spectral accelerations at low periods usually are adequate for predicting amplifications across the period range. This supports the current practice of using Peak Ground Acceleration in Ground Motion Models. However, for extremely soft sites, which demonstrate nonlinear effects well into the period range, it is better practice to ensure that amplification and spectral acceleration share the same oscillator period. Lastly, for predicting the peak shear strain (an important parameter that controls nonlinearity of site response), Peak Ground Velocity is generally adequate, and this conclusion is in line with the definition of proxy shear strain.

42 ENGINEERING↗