SEARCH · Engineering Papers

Results for “ROCKS”

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 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS2014

Alteration of the Carbon and Nitrogen Isotopic Composition in the Martian Surface Rocks Due to Cosmic Ray Exposure

C-13/C-12 and N-15/N-14 isotopic ratios are pivotal for our understanding of the Martian carbon cycle, history of the Martian atmospheric escape, and origin of the organic compounds on Mars. Here we demonstrate that the carbon and nitrogen isotopic composition of the surface rocks on Mars can be significantly altered by the continuous exposure of Martian surface to cosmic rays. Cosmic rays can effectively produce C-13 and N-15 isotopes via spallation nuclear reactions on oxygen atoms in various Martian rocks. We calculate that in the top meter of the Martian rocks, the rates of production of both C-13 and N-15 due to galactic cosmic rays (GCRs) exposure can vary within 1.5-6 atoms/cm3/s depending on rocks' depth and chemical composition. We also find that the average solar cosmic rays can produce carbon and nitrogen isotopes at a rate comparable to GCRs in the top 5-10 cm of the Martian rocks. We demonstrate that if the total carbon content in a surface Martian rock is <10 ppm, then the "light," potentially "biological" C-13/C-12 ratio would be effectively erased by cosmic rays over 3.5 billion years of exposure. We found that for the rocks with relatively short exposure ages (e.g., 100 million years), cosmogenic changes in N-15/N-14 ratio are still very significant. We also show that a short exposure to cosmic rays of Allan Hills 84001 while on Mars can explain its high-temperature heavy nitrogen isotopic composition (N-15/N-14). Applications to Martian meteorites and the current Mars Science Laboratory mission are discussed.

cosmic ray exposure↗

NASA NTRS2017

Constraining the Origin of Basaltic Volcanic Rocks Observed by Opportunity Along the Rim of Endeavour Crater

The Mars Exploration Rover (MER) Opportunity continues its exploration along the rim of Endeavour Crater. While the primary focus for investigation has been to seek evidence of aqueous alteration, Opportunity has observed a variety of rock types, including some that are hard and relatively unaltered. These rocks tend to occur most commonly as "float rocks" or "erratics" where the geologic setting does not clearly reveal their origin. Along the rim of Endeavour crater (Fig. 1), such rocks, commonly noted in Panoramic Camera (Pancam) left eye composites as "blue rocks", are abundant components of some of the Endeavour crater rim deposits, scree slopes, and colluvium deposits. In this abstract, we examine the similarity of several of these rocks analyzed using Opportunity's Alpha Particle X-Ray Spectrometer (APXS), images and color from the Pancam, and textures observed with the Microscopic Imager (MI. At issue is the blue rocks origin; are they impact melt or volcanic, what is their age relative to Endeavour crater, and how they are related to each other?

Bouchard, M. C.↗

NASA NTRS2021

Near-infrared Reflectance of Rocks at high Temperature: Preliminary Results and Implications for Near-infrared Emissivity of Venus's Surface

Light emitted from Venus’s surface can be viewed through spectral “windows” in its atmosphere, in the nearinfrared (NIR) around 1000 nanometer (nm) wavelengths. The NIR emissivity of Venus’s surface can constrain rock types and their weathering state; emissivities can be measured directly or calculated from reflectances. We measured the reflectances of igneous and sedimentary rocks at Venus’s surface temperature, 400 °C–500 °C at 850 and 950 nm; samples were heated in a box furnace in air, illuminated by light-emitting diodes (LEDs), and imaged with a modified charge-coupled device (CCD) camera. Reflectances were also measured at 25 °C from 350 to 1400 nm. Rock reflectances at 850 and 950 nm and 400 °C–500 °C are nearly identical to those at 25 °C, except for the effects of nanophase hematite forming on some surfaces. Fresh basalts have reflectances (high and low temperatures) near 7.5%; a leucogranite similarly has reflectances near 50%. Pigmentary hematite has nearly identical reflectances at high- and low-temperature at these wavelengths. Pigmentary hematite appears dark brown 400 °C–500 °C because its absorption edge has shifted to beyond the limit of human vision. These rock reflectances imply that basalts should have emissivities near 0.9, and granite (and similar felsic rocks) should have lower emissivities ∼0.5. Thus, basalt and felsic rock should be easily distinguished in NIR emissivity measurements of Venus’s surface, such as are baselined in recent Venus mission proposals. Other sedimentary rocks should have even lower emissivities: quartz sand at around ∼0.3, and anhydrite as low as 0.1.

Allan H Treiman↗

NASA NTRSDate not supplied

Pink Spinel in Apollo Impact Melt Rock 68815: Implications for Mg-Suite Magmatism

Introduction: Magnesian rocks from the lunar highlands are collectively termed the Mg-suite. Characterized by high (>60) Mg# (molar 100×Mg/[Mg+Fe]) in mafic phases and calcic plagioclase, these rocks are plutonic to hypabyssal in origin, and include a range of bulk mineralogies such as troctolites, dunites, norites, gabbronorites, and spinel troctolites [1]. These Mg-suite lithologies have distinct trace element concentrations and ratios that differentiate them from other lunar rock types. These rocks are ancient, generally dated to between 4.5 and 4.1 Ga, although it is unknown if this represents the full range of Mg-suite ages [1–2]. Among the Mg-suite lithologies, the spinel troctolites are relatively rare, to date only found in polymict breccias [3]. Spinel troctolites, as their name suggests, consist of calcic plagioclase and forsteritic olivine, with minor amounts of spinel (MgAl2O4), ± pyroxene and cordierite [1,4]. This form of spinel is often called ‘pink’ spinel because of its appearance in thin section under plane polarized light (PPL; Fig. 1), due to minor amounts of Cr. Spinel troctolites are generally plutonic or hypabyssal in origin (subsequently exhumed and incorporated into polymict breccias), or formed through impact processes (e.g., crystalline impact melt) [5]. A spinel-rich lithology has also been found in the Moscoviense region of the Moon via the Moon Mineralogy Mapper (M3) and lacks other mafic phases [6]. Finally, while the Mg-suite sampled thus far consists of plutonic (or hypabyssal) rocks, the question remains if such magmas could have erupted on the surface of the Moon [7]. These magmas have much lower density than mare basalts, but little sample or remote sensing evidence has been found to support the idea that extrusive Mg-suite volcanism occurred [7]. Here, we present a coordinated microanalytical study of spinel-bearing lithic and mineral clasts found in Apollo sample 68815. These data will be used to understand their petrogenesis (magmatic or impact) and modification histories, and to shed light on the existence of volcanic Mg-suite rocks. Sample Description: Apollo sample 68815 is a polymict impact melt breccia containing a variety of lithic and mineral fragments embedded in devitrified impact melt. This sample was chipped off the top of a boulder at Station 8 during the Apollo 16 mission and had an original weight of nearly 1.8 kg. In this study, we investigated two polished thin sections of 68815: 68815,17 and 68815,148, both containing spinel. Methods: The thin sections of 68815 were studied using optical light microscopy (PPL, cross-polarized light, and reflected light) with a Keyence VHX-7100 Digital Microscope. Each section was then X-ray mapped for 13–14 elements using a Cameca SX100 electron probe microanalyzer (EPMA) located in the Kuiper Materials Imaging and Characterization Facility (KMICF) at the University of Arizona. We have obtained geochemical information about the phases (olivine, plagioclase, spinel, pyroxene) in the thin sections also using the EPMA. In addition, we have used ThermoScientific Helios NanoLab 660 Focused-Ion-Beam Scanning-Electron Microscope (FIB-SEM) and a Hitachi S-4800 SEM (both in KMICF) to obtain backscattered electron (BSE) images and energy dispersive Xray spectrometry (EDS) maps of areas of interest. Using a JEOL 7900F SEM at the Astromaterials Research & Exploration Science (ARES) at NASA Johnson Space Center (JSC), we have obtained electron backscatter diffraction (EBSD) maps of the spinel-bearing portions of the thin sections. The EBSD data were collected under beam conditions of 20 kV, and ~90 μA, with step sizes varying from 0.05 to 2 μm. Following EBSD data collection, we processed the data using AZtecCrystal and MTEX, a free MATLAB toolbox. Results: We have found clasts with subophitic textures, that consist of primarily olivine and plagioclase, with minor amounts of pink Mg-Al spinel and pyroxene (Fig. 1). These clasts are up to ~1 mm in length and contain spinels up to 50 μm across. We have additionally identified pink Mg-Al spinels within the impact melt (i.e., not contained in lithic clasts) in both thin sections. In one instance, a single spinel grain is approximately 300 μm across (Fig. 1b, 2). The spinel fragments embedded in impact melt have varying compositions, typically distinct from the compositions of spinels in the lithic clasts. Spinel-Bearing Clasts: Ten lithic clasts with similar textures and mineral compositions were identified between 68815,17 (two clasts) and ,148 (eight clasts). These clasts fall into two groups. The first has skeletal olivine with intergranular plagioclase, with minor amounts of pyroxene and spinel (Fig. 1a, 2c, 2d). The spinel in these clasts are found amid the plagioclase. The second group have an intergranular texture of olivine and plagioclase, again with minor spinel and pyroxene. The second group may contain spinels surrounded by plagioclase, and spinels enclosed in olivine. Spinels located within both clast types range from no apparent Cr zoning, to reverse zoning (Cr-enrichment inward; Fig. 2d), to normal zoning (Cr-enrichment outward). In the clasts thus far investigated with EPMA, plagioclase compositions range from An# (molar 100×Ca/[Ca+Na+K]) 92–96. Olivine Mg# ranged from 77 to 94, while pyroxene had Mg# from 54–84. Spinel in the clasts have Cr# (molar 100×Cr/[Cr+Al]) 2–4 and Mg# 88–91, which is within the range of pristine and plutonic spinel troctolites [8]. Isolated Spinels: These crystals are generally euhedral to subhedral, and can exhibit reverse Cr zoning (Cr enrichment inward) or no apparent Cr zoning. The spinels thus far investigated via EPMA have Cr# 9–14 and Mg# 65–82. The Cr# for these spinels is within the range reported by [8], but have lower Mg#. Future Work: We will continue to process the EBSD data for these lithic and mineral clasts. We will also continue to characterize these clasts using EPMA and SEM. By thoroughly characterizing the various spinels and spinel-bearing clasts, we aim to constrain the petrogenesis of these minerals and rock fragments. Acknowledgments: We thank NASA for the loan of these thin sections. We thank Ken Domanik and Jerry Chang for their support with data collection. Work was supported by a University of Arizona RII Core Facilities Pilot Program grant and start-up funds to JJB. We acknowledge support from NASA’s Planetary Science Research program for analysis performed at JSC. References: [1] Shearer C. K. et al. (2015) Am. Min. 100, 294–325. [2] Borg L. E. et al. (2020) GCA 290, 312–332. [3] Warren P. H. (1993) Am. Min. 78, 360–376. [4] Dymek R. F. et al. (1976) LPS VII, 2335–2378. [5] Treiman et al. (2019) Am. Min. 104, 370–384. [6] Pieters et al. (2011) JGR: Plan. 116:E00G08. [7] Prissel et al. (2016) Icarus 277, 319–329. [8] Prissel et al. (2016) Am. Min. 101, 1624–1635.

spinel↗

NASA NTRSDate not supplied

Using X-Ray Computed Tomography to Catalog Rock Fragments in Apollo Drive Tube 73002

Overview: The Apollo missions collected 382 kg of rock, regolith, and core samples from six locations on the nearside of the Moon. Approximately 84% by mass of the Apollo collection remains in pristine condition within the curation facility at Johnson Space Center (i.e., never allocated, continuously stored in dry-N2 purged cabinets, exposure history restricted to Teflon, stainless steel, and Al-metal). Although most Apollo samples have been well characterized, there are several types of samples that have remained wholly or largely unstudied since their return, and/or that have been cu-rated under special conditions, e.g., frozen samples, samples stored in He-purged environment, and previously unopened drive tubes. NASA solicited proposals for the Apollo Next Generation Sample Analysis Pro-gram (ANGSA), and 9 teams were selected to study a subset of the unopened and frozen samples [1]. The first sample opened as part of the ANGSA pro-gram was drive tube 73002. This was originally a ~30 cm long, 4 cm diameter drive tube collected on a land-slide deposit near Lara Crater at the Apollo 17 landing site. It was part of a ~60 cm long double drive tube collected, and the bottom half of the tube (73001) was sealed under vacuum on the Moon [2]. Prior to opening sample 73002, the sample was imaged with a high resolution X-ray Computed Tomography (XCT) scan of the entire tube [3], which provided invaluable information during the dissection process [4]. In addition to the pre-dissection XCT scans, individual >4 mm particles were separated from the 73002 regolith during processing and scanned at high resolution by XCT. Here we pre-sent the initial lithologic classification of 134 individual >4 mm rock fragments separated from the 73002 core during the dissection process. Methodology: Drive tube 73002 was manually dissected in 0.5 cm depth intervals in three passes (Fig. 1) [4,5]. Each interval from pass 1 and 2 was sieved to <1 mm and >1 mm size fractions, and >1 mm particles were further manually subdivided into 1-2, 2-4, 4-10, and >10 mm size fractions. Pass 3 was not sieved, but >10 mm clasts were separated manually. Each 4-10 mm and >10 mm fragment was individually weighed, triply bagged in Teflon, and scanned by XCT. There are 60 rock fragments in Pass 1, 64 rock fragments in Pass 2 (from the 4-10 mm and >10 mm size fractions), and 8 rock fragments from Pass 3 (>10 mm size fraction). Each individually bagged rock fragment was scanned using the 180 kV nano-focus transmission source on the Nikon XTH 320 XCT system at NASA Johnson Space Center [6]. Scanning conditions varied considerably for individual particles, in large part be-cause of the large variation in size (0.008-19.623 g). All scans fell within the following range of scan conditions: 2.8-20.6 um voxel size; 90-155 kV voltage; 18-39 uA current; 1891-3141 projections; and 902-2000 slices. Results and Discussion: The 132 rock fragments from sample 73002 fall into the following general cat-gories: agglutinates (n = 6); basalts (13); impact melts (5); impact melt breccias (IMB; 42); regolith breccias (62); and soil breccias (4); see Figure 2 for representative examples of each lithology. Within most of these broad lithologic groups are recognizable subgroups. For example, a significant portion of regolith breccia fragments contain some agglutinate-like glass (n = 9) or are dilithologic (7) because they contain a single large clast (~50% or more by volume). Subgroups can be based on similarities to previously identified lunar lithologies, such as high-Ti basalts (9) and VLT basalts (4; Fig. 3), or based on commonly seen features within the fragments, such as poikilitic ilmenite IMB (10), ilmenite-lath IMB (12), or vesicular IMB (11) “groups”. Particles in the same “group” are not necessarily intended to be genetically related, but rather identify particles that are similar and that follow-up studies can classify in more detail [7]. Conclusion: Identification of lithologies based on XCT is a powerful tool, but a more absolute classification will sometimes require additional textural information from thin sections (e.g., glassy-matrix regolith vs. impact-melt vs. granulitic breccia) or quantitative mineral compositions (e.g., basalt vs. monomict breccia).

R A Zeigler↗

DOE OSTI2023

Fabric controls on fracture surface roughness of an architected rock material

Fluid flow through fractures is intimately linked to the fracture surfaces that define the void geometry through which fluids flow. Thus, an understanding of what controls fracture surface roughness is essential to the development of models for predicting fluid transport through fractured rock. The difficulty in predicting surface roughness arises from the complexity of rock which is inherently heterogeneous and nonuniform in composition, fabric, and structural components, even when samples are acquired from the same rock mass. Here, a benchmarked-simulation approach motivated from geo-architected 3D printed synthetic gypsum rocks is used to provide insight into the competing contributions from fabric and layering on fracture roughness formation. Simulation results from a discrete element model (Particle Flow Code, Itasca Consulting Group, Inc.) clearly indicate that the relative orientation between mineral layers and in-layer mineral fabric, and the variability in mineral bonding strengths determine whether anisotropic corrugated surfaces or isotropic surfaces are formed. Weak mineral layers oriented perpendicular to the applied load resulted in strong roughness anisotropy. Peak failure loads were found to vary up to 30% depending on the strength of the mineral fabric at the location of fracture initiation, which provides insight into the observed high variability in strength values of natural rock. The uniqueness of induced fracture roughness and peak failure load is intimately linked to layering, mineral fabric, and their distribution in the rock. These findings have important implications for any architected material fabricated through serial printing of layers with local compositional heterogeneity.

3D printed rock↗

DOE OSTI2016

Fully coupled hydro-mechanical numerical manifold modeling of porous rock with dominant fractures

Coupled hydro-mechanical (HM) processes are significant in geological engineering such as oil and gas extraction, geothermal energy, nuclear waste disposal and for the safety assessment of dam foundations and rock slopes, where the geological media usually consist of fractured rock masses. We developed a model for the analysis of coupled hydro-mechanical processes in porous rock containing dominant fractures, by using the numerical manifold method (NMM). In the current model, the fractures are regarded as different material domains from surrounding rock, i.e., finite-thickness fracture zones as porous media. Compared with the rock matrix, these fractured porous media are characterized with nonlinear behavior of hydraulic and mechanical properties, involving not only direct (poroelastic) coupling but also indirect (property change) coupling. By combining the potential energy associated with mechanical responses, fluid flow and solid–fluid interactions, a new formulation for direct HM coupling in porous media is established. For indirect coupling associated with fracture opening/closure, we developed a new approach implicitly considering the nonlinear properties by directly assembling the corresponding strain energy. Compared with traditional methods with approximation of the nonlinear constitutive equations, this new formulation achieves a more accurate representation of the nonlinear behavior. We implemented the new model for coupled HM analysis in NMM, which has fixed mathematical grid and accurate integration, and developed a new computer code. We tested the code for direct coupling on two classical poroelastic problems with coarse mesh and compared the results with the analytical solutions, achieving excellent agreement, respectively. We tested for indirect coupling on models with a single dominant fracture and obtained reasonable results. The current poroelastic NNM model with a continuous finite-thickness fracture zone will be further developed considering thin fractures in a discontinuous approach for a comprehensive model for HM analysis in fractured porous rock masses.

58 GEOSCIENCES↗

DOE OSTI2026

Geochemical evolution in Cacapon member: Fluid-rock interaction experiments and model insights for Appalachian Basin geothermal development

Here, this study combines recirculated flow-through experimental results conducted for 17 days at 90C and 200 PSI with reactive transport modeling to estimate fluid-rock interactions occurring in a sandy mudstone using an interbedded sandstone-shale sample from the Cacapon Member of the upper Tuscarora Sandstone/lower Rose Hill Formation for the purpose of geothermal exploration. Results suggest that the fluid and rock are likely to be in or near partial-equilibrium after approximately one year. In addition, after >400 h of continuous injection at 0.05 mL/min (5*10 −8 m 3 /min) the reactive front is restricted to the first ∼13% of the 4 cm*1.6 cm experimental rock length, whereas after >9100 h (∼1 year), the reactive front extends to 30% total length. The rate of changes in dissolution or precipitation are however, very minimal, with all major rock-forming minerals having rates <10 −11 mol/L porous media/s. Reservoir rock in the presence of dilute brine as may be the case during operation of an enhanced geothermal system would experience little alteration during the shut-in phase, and possibly up to one year. These results have utility in geothermal exploration for reservoirs at similar temperatures as well as general fluid-sandy mudstone rock interaction geochemistry.

Appalachian Basin↗

DOE OSTI2019

A visco-plastic framework for interface processes in sedimentary reservoir rocks at HPHT conditions

As energy operations face the challenge of reservoirs at ever-increasing depths, modelling the response of reservoir rocks at high pressure and high temperature (HPHT) conditions is a crucial step for successfully unlocking new resources. At these conditions, rocks can experience thermal and pressure sensitivity, as well as admit internal phase changes at the pore-solid interfaces that determine their macroscopic response to external loading. It is therefore expected that constitutive laws of Geomechanics should be enriched to accommodate such effects. In this work, a multi-physics constitutive theory for sedimentary rocks is proposed within the general framework of viscoplasticity. The viscosity of the material is assumed to be a function of the temperature, pore-pressure and energy required to alter the inter-granular interfaces. This energy is expressed through the chemical potentials of the phases involved during processes like chemical dissolution/precipitation of the interfaces or mechanical debonding. The resulting flow law and corresponding stress equilibrium are coupled to the energy and mass conservation laws, constituting a closed system of equations, which is solved using the Finite Element simulator REDBACK. A series of numerical tests for performance and calibration is then successfully performed against different types of reservoir rocks, confining pressures and temperatures (from room temperature to over 800K). Here, we show that the mechanical response of sedimentary porous rocks at strains usually achieved in laboratory testing can be explained by accounting for the interface processes taking place at the cementitious material bonding the grains, a process that can also determine the brittle-to-ductile transition of sedimentary rocks.

58 GEOSCIENCES↗

DOE OSTI2015

A study of changes in deep fractured rock permeability due to coupled hydro-mechanical effects

This paper presents a numerical study of the hydro-mechanical behaviour of a fractured rock domain at 1000. m depth below the land surface as a function of different levels of fluid pore pressure. A 2D fractured rock domain is adopted based on data obtained from outcrop mapping, displaying multiple fracture sets, fracture intersections, dead-end and curved fractures. A continuum based numerical model is used to evaluate the effects of compressive boundary stresses, cracking by tension failure in the intact rock and fractures and shear displacement along fractures on its equivalent permeability. Two in situ stress boundary conditions are considered: an isotropic case SR1 with the two horizontal boundary compressive stresses having the same magnitude, and an anisotropic case SR2 with the ratio between these compressive stress components set to be 2. In the SR2 case, changes in the local stress and stress ratio distributions due to different fluid pore pressure levels are anisotropic and more significant than in the SR1 case, because of tension failures in the intact rock forming bridges between fractures. These failure regions opened new flow connections between fractures and thereby caused important anisotropic changes in the flow paths, and significant decrease in local gradients of fluid pore pressure. The equivalent permeability increases sharply when the fluid pore pressure is approximately 90% of the magnitude of the minimum stress at the boundaries of the fractured rock domain. Results show that the equivalent permeability of the fractured rock domain is most sensitive to the fractures normal stiffness, the permeability of the tension failure regions and the power-law exponent for permeability change.

58 GEOSCIENCES↗

DOE OSTI2014

Fractured rock stress‐permeability relationships from in situ data and effects of temperature and chemical‐mechanical couplings

Abstract The purpose of this paper is to (i) review field data on stress‐induced permeability changes in fractured rock; (ii) describe estimation of fractured rock stress‐permeability relationships through model calibration against such field data; and (iii) discuss observations of temperature and chemically mediated fracture closure and its effect on fractured rock permeability. The field data that are reviewed include in situ block experiments, excavation‐induced changes in permeability around tunnels, borehole injection experiments, depth (and stress) dependent permeability, and permeability changes associated with a large‐scale rock‐mass heating experiment. Data show how the stress‐permeability relationship of fractured rock very much depends on local in situ conditions, such as fracture shear offset and fracture infilling by mineral precipitation. Field and laboratory experiments involving temperature have shown significant temperature‐driven fracture closure even under constant stress. Such temperature‐driven fracture closure has been described as thermal overclosure and relates to better fitting of opposing fracture surfaces at high temperatures, or is attributed to chemically mediated fracture closure related to pressure solution (and compaction) of stressed fracture surface asperities. Back‐calculated stress‐permeability relationships from field data may implicitly account for such effects, but the relative contribution of purely thermal‐mechanical and chemically mediated changes is difficult to isolate. Therefore, it is concluded that further laboratory and in situ experiments are needed to increase the knowledge of the true mechanisms behind thermally driven fracture closure, and to further assess the importance of chemical‐mechanical coupling for the long‐term evolution of fractured rock permeability.

58 GEOSCIENCES↗

DOE OSTI2019

A scale-consistent method for imaging porosity and micrite in dual-porosity carbonate rocks

Unlike many other clastic rocks, relating velocity and permeability to porosity for micrite-bearing carbonate rocks has been largely unsuccessful. Recent studies have shown that additional parameters, most notably the distribution and/or proportion of micrite, can be used to parameterize the velocity and permeability behavior. However, there is currently no scale-consistent, 3D methodology for differentiating macroporosity and microporosity from the total porosity measured on bench-top laboratory equipment. Previous studies estimated microporosity and micrite content by combining total porosity measurements conducted on whole 50 mm cores with measurements of phase volumes on 1 mm digital rocks (i.e., scale-inconsistent). As a step forward from those, we imaged dual-porosity carbonate rocks using X-ray microcomputed tomography and then leveraged a recently developed, optimization-based technique, called data-constrained modeling, to map the macroporosity and microporosity distribution of our samples. We evaluate the volumetric proportions of macropores, micropores, and coarse-grained calcite as a function of micrite content — with their respective uncertainties — all measured on the same digital rock and with the same method. Finally, we determine how measurements of the volumetric phase proportions could be extended using standard effective medium models to predict reservoir physical properties. Here, the sensitivity of these models to the proportion of micrite and microporosity within the micrite is evidence that the nonuniqueness among permeability, velocity, and porosity that is commonly observed of micrite-bearing carbonate rocks can be explained by a variation of micrite content and microporosity at a similar porosity.

58 GEOSCIENCES↗

DOE OSTI2024

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

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

58 GEOSCIENCES↗

DOE OSTI2022

Experimental study of compaction localization in carbonate rock and constitutive modeling of mechanical anisotropy

Abstract Sedimentary rocks are inherently anisotropic and prone to strain localization. While the influence of rock anisotropy on the brittle/dilative regime has been studied extensively, its influence on the ductile/compactive regime is much less explored. This paper discusses the anisotropic behavior of a high‐porosity carbonate rock from central Europe (the Maastricht Tuffeau). A set of triaxial tests with concurrent x‐ray tomography has been performed at different confining pressures. The anisotropic characteristics of this rock have been investigated by testing samples cored at different inclinations of the bedding, thus revealing non‐negligible effects of the coring direction on yielding and compaction behavior. Specifically, samples cored perpendicular to bedding display higher strength and longer stages of post‐yielding deformation before manifesting re‐hardening. Despite such alterations of the inelastic response, Digital Image Correlation has revealed that the strain localization mode is independent of the coring direction, thus being primarily affected by the confinement level. To capture the observed interaction between material anisotropy and compaction behavior at the continuum‐scale, an elastoplastic constitutive law has been proposed. For this purpose, a set of tensorial bases has been introduced to replicate how the oriented rock fabric modulates the yielding and plastic flow characteristics of the material. The analyses show that the impact of the coring direction on yield function and plastic flow rule is fundamentally different, thus requiring the use of distinct projection strategies (a strategy here defined heterotopic mapping). The performance of the model, studied through parametric analyses and by calibrating the experimental results, illustrates the improved capability of the proposed constitutive approach when applied to strongly anisotropic porous rocks.

Shahin, Ghassan↗

DOE OSTI2017

Characterization and differentiation of rock varnish types from different environments by microanalytical techniques

We investigated rock varnishes collected from several locations and environments worldwide by a wide range of microanalytical techniques. These techniques were selected to address the challenges posed by the chemical and structural complexity within the micrometer- to nanometer-sized structures in these geological materials. Femtosecond laser ablation-inductively coupled plasma-mass spectrometry (fs LA-ICP-MS), scanning transmission X-ray microscopy-near edge X-ray adsorption fine structure spectroscopy (STXM-NEXAFS) in combination with scanning electron microscopy (SEM) of focused ion beam (FIB) ultra-thin (100–200 nm) sections, conventional and polarization microscopy, as well as electron paramagnetic resonance (EPR) measurements were used to obtain information about these rock varnishes. Rock varnishes from different environments, which cannot readily be distinguished based on their macroscopic appearance, differ significantly in their constituent elemental mass fractions, e.g., of Mn, Fe, Ni, Co, Ba, and Pb, and their rare earth element (REE) patterns. Structural characteristics such as the particle sizes of embedded dust grains, internal structures such as layers of Mn-, Fe-, and Ca -rich material, and structures such as cavities varied between varnishes from different environments and regions in the world. The EPR spectra were consistent with aged biogenic Mn oxides in all samples, but showed subtle differences between samples of different origin. Our observations allow us to separate rock varnishes into different types, with differences that might be indicators of distinct geneses. Five different types of rock varnish could be distinguished, Type I–V, of which only Type I might be used as potential paleoclimate archive. Each varnish type has specific characteristics in terms of their elemental composition, element distribution, and structures. The combination of element ratios (Mn/Ba, Al/Ni, Mn/REY, Mn/Ce, Mn/Pb, La N /Yb N , and Ce/Ce*), total REE contents, and structures can be used to separate the different types of rock varnish from each other.

58 GEOSCIENCES↗

DOE OSTI2015

Finite element modeling of rock cutting and its fragmentation process

Rock cutting is a challenging problem from a modeling perspective. The challenges come from the complexity of the physics from the tool-rock interaction to the fracture process and propagation of the quasibrittle rocks. This study was aimed at developing a finite element procedure that was capable of providing reasonable estimates of cutting forces and, at the same time, capturing the essential characteristics of the fragmentation process. Published laboratory rock scratch tests were used as modeling targets since these tests encompass all essential characteristics of rock. Both shallow cuts and deep cuts from a rectangular cutter were analyzed first, and followed by modeling of shallow cuts from a disc cutter. It was concluded that rock cutting could be reasonably modeled by using a plasticity model, an element erosion scheme that remove an element when its energy release equals fracture energy, together with a proper selection of modeling parameters.

02 PETROLEUM↗

DOE OSTI2014

Use of nanotomographic images for structure analysis of carbonate rocks

Carbonate rocks store more than 50% of world's petroleum. These rocks' structures are highly complex and vary depending on many factors regarding their formation, e.g., lithification and diagenesis. In order to perform an effective extraction of petroleum it is necessary to know petrophysical parameters, such as total porosity, pore size and permeability of the reservoir rocks. Carbonate rocks usually have a range of pore sizes that goes from nanometers to meters or even dozen of meters. The nanopores and micropores might play an important role in the pores connectivity of carbonate rocks. X-ray computed tomography (CT) has been widely used to analyze petrophysical parameters in recent years. This technique has the capability to generate 2D images of the samples' inner structure and also allows the 3D reconstruction of the actual analyzed volume. CT is a powerful technique, but its results depend on the spatial resolution of the generated image. Spatial resolution is a measurement parameter that indicates the smallest object that can be detected. There are great difficulties to generate images with nanoscale resolution (nanotomographic images). In this work three carbonate rocks, one dolomite and two limestones (that will be called limestone A and limestone B) were analyzed by nanotomography. The measurements were performed with the SkyScan2011 nanotomograph, operated at 60 kV and 200 μA to measure the dolomite sample and 40 kV and 200 μA to measure the limestone samples. Each sample was measured with a given spatial resolution (270 nm for the dolomite sample, 360 nm for limestone A and 450 nm for limestone B). The achieved results for total porosity were: 3.09 % for dolomite, 0.65% for limestone A and 3.74% for limestone B. This paper reports the difficulties to acquire nanotomographic images and further analysis about the samples' pore sizes.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Validation of earthquake ground-motion models in southern California, USA, using precariously balanced rocks

Accurate estimates of earthquake ground shaking rely on uncertain ground-motion models derived from limited instrumental recordings of historical earthquakes. A critical issue is that there is currently no method to empirically validate the resultant ground-motion estimates of these models at the timescale of rare, large earthquakes; this lack of validation causes great uncertainty in ground-motion estimates. Here, we address this issue and validate ground-motion estimates for southern California utilizing the unexceeded ground motions recorded by 20 precariously balanced rocks. We used cosmogenic 10Be exposure dating to model the age of the precariously balanced rocks, which ranged from ca. 1 ka to ca. 50 ka, and calculated their probability of toppling at different ground-motion levels. With this rock data, we then validated the earthquake ground motions estimated by the Uniform California Earthquake Rupture Forecast, Version 3 (UCERF3) seismic-source characterization and the Next Generation Attenuation (NGA)-West2 ground-motion models. We found that no ground-motion model estimated levels of earthquake ground shaking consistent with the observed continued existence of all 20 precariously balanced rocks. The ground-motion model I14 estimated ground-motion levels that were inconsistent with the most rocks; therefore, I14 was invalidated and removed. At a 2475 year mean return period, the removal of this invalid ground-motion model resulted in a 2−7% reduction in the mean and a 10−36% reduction in the 5th−95th fractile uncertainty of the ground-motion estimates. Our findings demonstrate the value of empirical data from precariously balanced rocks as a validation tool for removing invalid ground-motion models and, in turn, reducing the uncertainty in earthquake ground-motion estimates.

Geology↗