Comparison of Mineralogy Variability and Major, Minor, Trace Element Compositions in Dunite and Peridotite Rock Samples with Use of Chemical, Optical and Spectroscopic Methods
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Here, we propose an innovative containment strategy based on time-dependent soft solids, namely 2 wt% Laponite suspensions, to reinforce natural subsurface seals and engineer flow barriers, with an eye toward making H 2 subsurface storage scalable and geographically agnostic. This suspension can be injected at its initial low viscosity and elasticity into a porous medium, allowing for easy pumping and targeted delivery. Once inside the target zone, it matures into a soft solid with much higher viscosity and elasticity, acting as a potential flow barrier. We discuss the rheological properties of the suspensions and demonstrate that hydrogen does not adversely affect their microstructure, but rather increases the suspensions’ viscosity and elasticity. Moreover, the suspensions enhance the rock samples’ compressive strength, while hydrogen exposure increases their stiffness and ductility. The ability of rock samples saturated with the suspensions to contain higher injected gas pressures is enhanced by aging at higher temperatures.
This study explores the geochemical reactions that can cause permeability loss in hydraulically fractured reservoirs. The experiments involved the reaction of powdered-rock samples with produced brines in batch reactor system at temperature of 95 °C and atmospheric pressure for 7-days and 30-days respectively. Results show changes in mineralogy and chemistry of rock and fluid samples respectively, therefore confirming chemical reactions between the two during the experiments. The mineralogical changes of the rock included decreases of pyrite and feldspar content, whilst carbonate and illite content showed an initial stability and increase respectively before decreasing. Results from analyses of post-reaction fluids generally corroborate the results obtained from mineralogical analyses. Integrating the results obtained from both rocks and fluids reveal a complex trend of reactions between rock and fluid samples which is summarized as follows. Dissolution of pyrite by oxygenated fluid causes transient and localized acidity which triggers the dissolution of feldspar, carbonates, and other minerals susceptible to dissolution under acidic conditions. The dissolution of minerals releases high concentrations of ions, some of which subsequently precipitate secondary minerals. On the field scale, the formation of secondary minerals in the pores and flow paths of hydrocarbons can cause significant reduction in the permeability of the reservoir, which will culminate in rapid productivity decline. This study provides an understanding of the geochemical rock–fluid reactions that impact long term permeability of shale reservoirs.
Measurements of lead isotopic compositions in detrital K-feldspar have been increasingly used as a tool to assess sediment provenance. We compiled a database of previously published Pb isotope data from 700 bedrock K-feldspar samples and 1,423 age-corrected bedrock whole rock samples from western North American igneous and metamorphic bodies. Additionally, we report 66 new K-feldspar Pb isotope data for plutons throughout the Pacific Northwest region of the United States and British Columbia. Results show that the Pb isotope values of plutonic K-feldspar depend on the isotopically juvenile or evolved nature of underlying crust. Samples obtained from the mid Cretaceous – mid Eocene Coast Plutonic Complex, North Cascades, and Intermontane superterrane that occur west of the initial 87 Sr/ 86 Sr (Sri) = 0.706 isopleth exhibit a highly restricted 207 Pb/ 206 Pb and 208 Pb/ 206 Pb values centred upon 0.83 and 2.03, respectively. Conversely, rocks overlying older continental crust further east such as the Middle Jurassic – Late Cretaceous Omineca crystalline belt, Idaho batholith, and Boulder Batholith exhibit far greater variation of Pb isotope values that parallel the 100 Ma isochron calculated from a two-stage Pb evolution model. We demonstrate that Pb isotopic results from the Idaho and Boulder Batholith region can be used to define distinctive subregions for Pb isotopic provenance analysis, and compare these signatures to the Mojave-Salinian batholith of southern California and western Arizona, as these two areas have previously been proposed as source regions for extraregional sediment that was deposited within the Nanaimo Basin during the Campanian – Maastrichtian. Future Pb isotopic analysis of detrital K-feldspar from the Nanaimo Basin of southwestern British Columbia may effectively distinguish between potential extraregional sources separated by thousands of kilometres.
Geological samples are inherently multi-scale. Understanding their bulk physical and chemical properties requires characterization down to the nano-scale. A powerful technique to study the three-dimensional microstructure is X-ray tomography, but it lacks information about the chemistry of samples. To develop a methodology for measuring the multi-scale 3D microstructure of geological samples, correlative X-ray micro- and nanotomography were performed on two rocks followed by scanning electron microscopy with energy-dispersive spectroscopy (SEM-EDS) analysis. The study was performed in five steps: (i) micro X-ray tomography was performed on rock sample cores, (ii) samples for nanotomography were prepared using laser milling, (iii) nanotomography was performed on the milled sub-samples, (iv) samples were mounted and polished for SEM analysis and (v) SEM imaging and compositional mapping was performed on micro and nanotomography samples for complimentary information. Correlative study performed on samples of serpentine and basalt revealed multiscale 3D structures involving both solid mineral phases and pore networks. Significant differences in the volume fraction of pores and mineral phases were also observed dependent on the imaging spatial resolution employed. This highlights the necessity for the application of such a multiscale approach for the characterization of complex aggregates such as rocks. Information acquired from the chemical mapping of different phases was also helpful in segmentation of phases that did not exhibit significant contrast in X-ray imaging. Adoption of the protocol used in this study can be broadly applied to 3D imaging studies being performed at the Advanced Light Source and other user facilities.
Underground Hydrogen Storage (UHS) in depleted oil and gas reservoirs could provide a cost-effective solution to balance seasonal fluctuations in renewable energy generation. However, data and knowledge on UHS at subsurface conditions are limited so it is difficult to estimate how effective this type of storage could be. In this study, we perform high pressure experiment to measure the effectiveness of cyclic hydrogen (H 2 ) storage in a specimen of Temblor sandstone retrieved from the San Joaquin Basin of California. Our experiment mimics reservoir pressure conditions to measure H 2 -brine relative permeability and fluid-rock interactions over the course of ten charging and discharging cycles. Initial gas breakthrough occurred at 15 % to 25 % H2 saturation in the specimen with 3 % NaCl brine as the resident fluid. Continuing injecting to 4 pore volumes (PV) of H 2 yielded an asymptotic H 2 saturation of 38 % to 41 %, a level often referred to as the irreducible gas saturation based on two-phase flow. The boundary condition in this study mimics the near wellbore region, which experiences bi-directional H 2 flow. This bi-directional flow led to evaporative drying of the specimen resulting in 94 % H 2 saturation at the end of 10th cycle. This indicates that cyclic flow and evaporative drying can lead to more efficient reservoir storage where a larger fraction of the reservoir porosity is usable to store H 2 . The produced gas stream consisted of H 2 mixed with 8 % to 22 % H 2 O, indicating formation dry-out by evaporation. Meanwhile, produced water chemistry indicated calcite and silicate dissolution, with calcite sourced from fossil fragments. This led to a loss of cementation and weakened the rock sample. Combined, our results indicate dry-out, compaction, increased H 2 saturation, rock weakening, and permeability loss during cyclic UHS. Overall, we anticipate that the combined effects should lead to higher than anticipated UHS storage efficiency per volume of sandstone reservoir rock.
We describe a “virtual mineral separation” method for measuring the cosmogenic 3 He concentration in pyroxene in mafic rocks that consist mainly of plagioclase and pyroxene, without physically separating the minerals. This approach is significantly faster and more cost-effective than the conventional method, which requires physical separation and purification of pyroxene grains by time-consuming and labor-intensive crushing, acid cleaning, magnetic separation, HF etching, and handpicking under a microscope. The premise of the method is that helium diffusivity is much higher in plagioclase than in pyroxene, so controlled preheating of a mixed whole-rock sample can degas 3 He from plagioclase while retaining all 3 He in pyroxene. A second heating step releases all 3 He from pyroxene for measurement. To then obtain a 3 He concentration in pyroxene rather than the whole rock, we determine the pyroxene weight fraction in the sample using X-ray computed tomography (CT). A comparison of 3 He concentrations in pyroxene measured using virtual mineral separation with those measured in the same samples by physical mineral separation in previous work shows no evidence of systematic bias between the methods. Virtual mineral separation greatly simplifies the workflow for 3 He exposure-dating of mafic rocks, reduces time, effort, and cost, and permits measurements on very small samples. This enables new emerging applications of exposure dating, such as quantifying stochastic surface processes, ecosystem studies, and potential subglacial bedrock exposure dating.
In this paper, we explore the application of convolutional neural networks (CNNs) for predicting the chemical composition of complex geologic samples in a simulated Martian atmospheric environment. Specifically, we aim to characterize oxide weight percentages (wt.%) of rock samples analyzed by remote Laser-Induced Breakdown Spectroscopy (LIBS), framing the problem as a multi-target regression task . Neural networks trained on LIBS spectra are prone to overfitting due to high spectral complexity, limited labeled data, and measurement noise. While regularization is critical for improving generalization, common methods (e.g., ℓ 2 regularization) impose constraints not directly tied to data distribution properties. We propose a novel regularization method based on a specific ƒ-divergence induced by a graph-based estimator, designed to constrain the distributional discrepancy between predictions and targets. This regularizer serves a dual purpose: (a) mitigating overfitting by enforcing a constraint on the distributional difference between predictions and noisy targets, and (b) acting as an auxiliary loss that penalizes large divergences. To enable backpropagation, we develop a differentiable approximation of this particular ƒ-divergence, making the method feasible for neural networks. Experiments on ChemCam and SuperCam LIBS calibration spectra show that mathematical equation-divergence regularization outperforms or matches standard regularization methods (ℓ 1 , ℓ 2 , dropout) and the classical baseline, partial least squares (PLS). Combining ƒ-divergence regularization with standard regularization yields further performance gains, indicating that distributional regularization is useful in this context giving a promising direction for robust model training in planetary science applications. Source code is publicly available at Klein and Li (2025), https://doi.org/10.11578/dc.20250530.7.
Abstract. In situ cosmogenic 36Cl analysis by accelerator mass spectrometry (AMS) is routinely employed to date Quaternary surfaces and assess rates of landscape evolution. However, standard laboratory preparation procedures for 36Cl dating require the addition of large amounts of isotopically enriched chlorine spike solution; these solutions are expensive and increasingly difficult to acquire from commercial sources. In addition, the typical workflow for 36Cl dating involves measuring both 35Cl/37Cl and 36Cl/Cl concurrently on the high-energy (post-accelerator) end of the AMS system, but 35Cl/37Cl determinations using this technique can be complicated by isotope fractionation and system memory during measurement. The traditional workflow also does not provide 36Cl extraction laboratories with the data needed to calculate native Cl concentrations in advance of 36Cl/Cl measurements. In light of these concerns, we present an improved workflow for extracting and measuring chlorine in geologic materials. Our initial step is to characterize 35Cl/37Cl on sample aliquots of up to ∼1 g prepared in Ag(Cl, Br) matrices, which greatly reduces the amount of isotopically enriched spike solution required to measure native Cl content in each sample. To avoid potential issues with isotope fractionation through the accelerator, 35Cl/37Cl is measured on the low-energy, pre-accelerator end of the AMS line. Then, for 36Cl/Cl measurements, we extract Cl as AgCl or Ag(Cl, Br) in analytical batches with a consistent total Cl load across all samples; this step is intended to minimize source memory effects during 36Cl/Cl measurements and allows the preparation of AMS standards that are customized to match known Cl contents in the samples. To assess the efficacy of this extraction and measurement workflow, we compare chlorine isotope ratio measurements on seven geologic samples prepared using standard procedures and the updated workflow. Measurements of 35Cl/37Cl and 36Cl/Cl are consistent between the two workflows, and 35Cl/37Cl values measured using our methods have considerably higher precision than those measured following standard protocols. The chemical preparation and measurement workflow presented here (1) reduces the amount of isotopically enriched chlorine spike used per rock sample by up to 95 %; (2) identifies rocks with high native Cl concentrations, which may be lower priority for 36Cl surface exposure dating, at an early stage of analysis; and (3) allows laboratory users to maintain control over the total chlorine content within and across analytical batches. These methods can be incorporated into existing laboratory and AMS protocols for 36Cl analyses and will increase the accessibility of 36Cl dating for geologic applications.
Ex-situ carbon mineralization is a process by which CO2 is reacted with alkaline silicate minerals and rocks to produce stable carbonate materials, which can be used for other industrial processes. Arizona, U.S.A., hosts abundant surficial mafic rocks in three young volcanic fields, Geronimo-San Bernardino, San Francisco, and Springerville, and other distributed locations throughout the state. We created a Mafic Rock Resource Inventory (MRRI) that categorizes geochemical, physical, and textural characteristics of a diverse suite of surficial mafic rock samples and provide a benchmark reaction dataset parameterizing the temperature, pressure, and pH conditions best suited ex-situ mineralization in different rock types. MRRI data is publicly available online via a map-viewer. We establish two reaction condition sets, varied in temperature, pressure, and pH, where crystal-rich and glassy rocks reach maximum reaction extent and different carbonate phases are formed. Systematic ex-situ mineralization experiments on 21 diverse rock types show trends in geochemical, mineralogical, and reactivity behavior and establish maximum effective capture capacity. From this, scoria cones in three Arizona volcanic fields have a ~62 Gt effective CO₂ storage capacity with one of the fields having a ~42 Gt storage capacity in lava flows. Reactivity results have applications to alkaline mafic rock resources exposed globally, including producing additional effective storage capacity estimates and scaled commercialization of mafic rock ex-situ mineralization, should reaction extents be improved through advances in mineralization techniques. MRRI data were used to create a Direct Air Capture to Mineralization (DACM) systems model, technoeconomic analysis, and life-cycle assessment. These documents are presented as three appendices.
Electrofracturing deeply buried shale formations could be used to increase reservoir permeability and improve reservoir production without requiring large volumes of freshwater. This paper describes a novel experimental system and initial test results to electrofracture shale under high confining pressures. Core-scale laboratory testing was performed on twelve rock samples recovered from a shale gas reservoir. Each sample was subjected to confining pressures of 20.7 MPa (3000 psi) or 58.6 MPa (8000 psi), representative of overburden pressures at depth. Samples were then subjected to application of high voltage until specimen fracture. The experiments produced deformed samples with multiple fracture types, both parallel and oblique to bedding planes. Electrofracturing increased permeabilities by up to nine orders of magnitude for extended time periods. Rock fracture and throughgoing fractures were demonstrated. Computed tomography images revealed the creation of fractures and tube/tunnel flow channels, which resisted closure under hydrostatic pressures up to 58.6 MPa. The breakdown energy and permeability changes in the sample were independent of applied confining pressure. The cumulative energy input required for fracture depended on applied confining pressure and sample length. The energy required to fracture samples up to 9 cm in length is generally more than 0.5 kJ/cm, but no greater than 1 kJ/cm. Our results show that electrofracture of shales under confining pressure is possible and could be a possible water-free mechanism for reservoir stimulation.
Accurate permeability predictions remain pivotal for understanding fluid flow in porous media, influencing crucial operations across petroleum engineering, hydrogeology, and related fields. Traditional approaches, while robust, often grapple with the inherent heterogeneity of reservoir rocks. With the advent of deep learning, convolutional neural networks (CNNs) have emerged as potent tools in image-based permeability estimation, capitalizing on micro-CT scans and digital rock imagery. This paper introduces a novel paradigm, employing vision transformers (ViTs)—a recent advancement in computer vision—for this crucial task. ViTs, which segment images into fixed-sized patches and process them through transformer architectures, present a promising alternative to CNNs. We present a methodology for implementing ViTs for permeability prediction, its results on diverse rock samples, and a comparison against conventional CNNs. The prediction results suggest that, with adequate training data, ViTs can match or surpass the predictive accuracy of CNNs, especially in rocks exhibiting significant heterogeneity. This study underscores the potential of ViTs as an innovative tool in permeability prediction, paving the way for further research and integration into mainstream reservoir characterization workflows.
Wave propagation and diffusive transport phenomena in a geological rock sample undergoing crack evolution process are expected to interact with the mechanical discontinuities in the medium. The measurements of the signals associated with these phenomena can be used to assess and monitor the crack-driven micromechanical alterations in the rock. Different wave/diffusion phenomena, such as sonic propagation, pressure diffusion, and acoustic emission (AE), are sensitive to different elements of the mechanical discontinuities generated during the evolution of the crack clusters from initiation to coalescence. Sonic propagation, AE, and pressure diffusion monitoring have the potential to map the crack evolution because the transmitter-receiver arrays can be designed, arranged and tuned to (1) achieve maximum recovery of the scattered waveforms and travel times, (2) capture the later arrivals and multiple reflections, and (3) illuminate large rock volume. However, the structural/topological complexities of the mechanical discontinuities, complex distribution of the stress fields, complex mechanical alterations in media, and fluid redistribution in the crack system pose serious challenges for the detection and modeling of the crack evolution process (from here on, we will use the term ‘crack evolution process’ to mean that the crack evolution occurred under shallow crustal conditions). For purposes of accurately accounting such complexities and heterogeneities in the absence of reliable physical laws, simulation methods, and signal processing techniques, my early-career research proposal will develop and apply novel data-driven machine learning methods to: (1) extract signal components relevant to the various phases of crack evolution and (2) generate a 2D visual map of the crack evolution process.
Microwave-assisted hydrogen generation from hydrocarbon-bearing reservoir rocks is strongly influenced by mineralogy, methane activation, carbonate reactions, and thermal runaway behavior. This study investigates a new approach in which carbonate phases are generated in-situ through the reaction of internally produced CO2 with Ca(OH)2 under microwave heating conditions. The objective is to evaluate how rock mineralogy, methane injection, and Ca(OH)2 addition influence hydrogen generation, carbon redistribution, and stage-dependent reaction pathways during microwave exposure. Microwave heating experiments were conducted using Permian Basin reservoir rocks under three experimental conditions: rock-only experiments under Ar atmosphere, CH4–Ar experiments without additive, and CH4–Ar experiments containing 5 wt% Ca(OH)2. Methane-assisted experiments were performed under continuous injection of 30 standard cubic centimeters per minute (sccm) CH4 and 30 sccm Ar. Based on thermal runaway behavior, each experiment was divided into three operational stages: Before Thermal Runaway (BR), After Thermal Runaway–Decrease in Microwave Power (ARD), and After Thermal Runaway– Increase in Microwave Power (ARI). Temperature and gas composition were continuously monitored throughout the experiments. The rock-only experiments demonstrated that hydrogen generation can occur intrinsically from hydrocarbon-bearing rocks under microwave heating, even without externally injected methane. However, hydrogen production did not correlate solely with kerogen content, indicating that mineralogy strongly influences hydrogen-generation pathways. Correlation analyses suggested that kerogen decomposition initially generated CH4, CO, and CO2, followed by secondary hydrocarbon reactions associated with H2 and C2 hydrocarbon formation. Methane-assisted experiments substantially increased hydrogen production; however, identical methane injection rates produced significantly different hydrogen yields among the rock samples, confirming that mineralogical composition controls methaneconversion behavior under microwave heating conditions. The addition of Ca(OH)2 significantly altered carbon evolution behavior in a stage-dependent manner. During the BR stage, Ca(OH)2 reduced gas-phase CO2 production, particularly in carbonate-rich rocks, indicating favorable conditions for in-situ carbonation and carbonate deposition prior to extensive thermal decomposition. The suppression of CO2 during BR became more pronounced with increasing carbonate content of the rock system. After thermal runaway, carbonate-containing systems exhibited enhanced hydrogen generation behavior, suggesting that carbonatederived mineral transformations and carbonate-mediated reactions contribute to high-temperature hydrogen-generation pathways. In carbonate-poor rocks, Ca(OH)2 addition enabled simultaneous URTeC 4493775 2 enhancement of hydrogen production and partial suppression of CO2 release during the post-runaway stages. Overall, the results demonstrate that microwave-assisted hydrogen generation is governed by dynamically evolving interactions among kerogen decomposition, methane activation, mineral transformations, carbonate formation/decomposition, and thermal runaway behavior. This work introduces in-situ carbonate engineering with Ca(OH)2 as a strategy for coupling hydrogen generation with partial insitu carbon management under microwave heating conditions.
In this paper we test new approaches for predicting the amount of element oxides in rock samples from the ChemCam instrument suite onboard the NASA Curiosity rover by focusing on K 2 O. Using the expanded dataset compiled by Gasda et al. (2021) with and without the Earth to Mars (E2M and NoE2M) transformation discussed in Clegg et al. (2017) we trained blended submodels using the “double blending” technique and compared these to ensemble methods (Random Forest, ExtraTrees, and Gradient Boosting Regression). We found that ensemble methods performed similar to blended submodels when looking at RMSE-P on the laboratory spectra and provided significant advantages when looking at spectra coming from Mars. For the full model, blended submodels achieved an RMSE-P of 0.62 and 0.60 (E2M and NoE2M respectively) while Gradient Boosting Regression resulted in a slightly improved RMSE-P of 0.59 and 0.60. More importantly, by employing a local RMSE-P estimation technique where model performance is evaluated based on nearby test samples we found that using ensemble methods can lower the quantification limit for K 2 O from the current value of ≈0.6 wt% to ≈0.08 wt% using Extra Trees and Random Forest. This would allow for a much larger range of K 2 O values to be quantified on Mars with greater certainty given that most targets seen on Mars tend to have <1 wt% K2O. Finally, we used both Mean Decrease in Impurity (MDI) and permutation importance techniques to investigate the wavelengths used by the ensemble methods and found that they correspond to known potassium emission lines. This suggests that ensemble methods can provide an easier to train and improved alternative to blended submodels for predicting potassium compositions from Laser Induced Breakdown Spectroscopy (LIBS) data.
Abstract Geologic storage of CO 2 and H 2 are climate‐positive techniques for meeting the energy transition. While similar formations could be considered for both gases, the flow dynamics could differ due to differences in their thermophysical properties. We conduct a rigorous pore‐scale study of water/CO 2 and water/H 2 systems at relevant reservoir conditions in a Bentheimer rock sample using the lattice Boltzmann method to quantify the effects of capillary, viscous, inertial, and wetting forces during gas invasion. At similar conditions, H 2 invasion is weaker compared to CO 2 due to unfavorable viscosity ratios. Increasing flow rate, however, increases the breakthrough saturation for both gas systems in the range of capillary numbers studied. At isolated conditions of flow rate, viscosity ratio, and wettability, local inertial effects are found to be critical and show consistent increase in the invaded gas saturation. The effect of inertial forces persits for both gases across all field conditions tested.
This repository contains experimental data from a series of fluid injection-induced shearing tests conducted on Utah FORGE granitoid. The experiments were performed using an aluminum triaxial pressure vessel (TEMCO) apparatus at Pennsylvania State University. The primary aim was to investigate the impact of temperature on fault seismicity and to explore how different pre-stress ratios influence shearing behavior. The data includes results from experiments conducted under varying conditions of temperature, pre-stress ratios, and pore pressure increments. The rock samples used in these experiments were 60-grit granitoid with a single inclined fracture oriented at 60 degrees with respect to the horizontal cross-section. The dataset also includes measurements of normal stiffness at different temperatures and images of the experimental apparatus. The included README file details each experiential setup and outlines which data files represent which experimental conditions.
Young’s modulus is an important parameter for characterizing the strength of, and wave propagation through, a given material. This study improves the estimation of Young’s modulus using the impulse excitation technique based on an experimental analysis of 19 borosilicate glass bars. Analysis of the frequency equations relating Young’s modulus to the out-of-plane and in-plane flexural resonant frequencies of rectangular prisms has been conducted for both the fundamental frequency and its overtones at higher orders of vibration. Furthermore, the Young’s moduli of three novaculite rocks with various porosities were then measured up to the seventh order of vibration to validate the optimum shear constant equation for estimating Young’s modulus. Young’s modulus was found to be nearly frequency independent for these rock samples.