Engineering Papers⌕ Search

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.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Voucher Opportunity 5-15: Independent Assessment of Monitoring, Reporting, and Verification (MRV) Technologies and Practices for Enhanced Rock Weathering (CRADA 718) Abstract

Development of robust, transparent, and precise monitoring, reporting, and verification (MRV) technologies and practices is critical for carbon dioxide removal (CDR) project developers to comply with regulatory and permitting requirements, voluntary carbon market (VCM) protocols, and to ensure safety while reducing environmental impacts. Enhanced rock weathering (ERW)-based CDR technologies focus on removing atmospheric carbon through conversion into thermodynamically stable solid or aqueous carbonate forms for permanent storage (i.e., mineralization). This highly durable form of CDR enhances naturally occurring silicate rock weathering cycles by optimizing application of finely-ground silicate rock particles (i.e., from basalt) on terrestrial agricultural lands to accelerate natural silicate rock weathering and mineralization. Enhanced rock weathering may also provide improved crop yields and enhance soil health. A critical aspect for commercialization of these technologies is the development of MRV to quantify the net removal and durable storage of atmospheric CO 2 . For ERW systems, it is essential to accurately characterize the mineral feedstock selected for application to establish the baseline geochemical composition, mineral dissolution rates, and carbon removal potential of the feedstocks to estimate overall net removal. Given the difficulty with conducting MRV for ERW in diverse soil/environment types, over large application areas, and due to complex chemical reaction networks, this project will accelerate understanding towards consensus on best practices for MRV. The overall objectives of the proposed voucher project are to: 1) Characterize and analyze feedstock(s) intended for ERW field application by Lithos Carbon (“Voucher Recipient”/ “CRADA Participant”) to determine overall mineralization potential; 2) Facilitate knowledge transfer and documentation of experimental protocols, instrumentation, and other relevant best practices; and 3) Support the Voucher Recipient’s broader technology commercialization and ERW Research Facility development plans. This work will align with the Voucher Recipient’s MRV plans for field sites and build upon complementary efforts conducted by PNNL on mineralization MRV.

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

Geomechanical characterization of Rock Valley carbonates

In Rock Valley, Nevada, the USA, past earthquake sequences and shallow faulting in Paleozoic carbonates remain poorly understood. Carbonates are typically more ductile rocks that do not experience large stress drops or fracture coalescence, which contradicts previous observations in the region. As part of the Source Physics Experiment, an effort has been made to experimentally characterize the petrophysical and geomechanical properties of carbonates from Rock Valley. Well core and outcrop samples were used to determine the difference between shallow-buried Tertiary limestones and deeply buried Paleozoic limestones and dolostones. The carbonates possessed porosities between 3 and 9%, with V P and V S in the Tertiary carbonates around 2700 and 1800 m/s, respectively, and 6000 and 3100 m/s, respectively, in the Paleozoic carbonates. Microstructural characterization revealed that the Paleozoic carbonates contained significantly more pre-existing damage and alteration than the Tertiary carbonates, though the deformation varies from localized to diffuse. Unconfined Brazilian tests and triaxial tests with confining pressures between 0 and 50 MPa showed that, although samples all experienced failure, the dolostones typically failed at greater stresses and possessed greater E/ν ratios than the limestones. Furthermore, the Hoek–Brown failure criterion was used to construct failure envelopes with the compressive and tensile failure tests. Velocity, density, and porosity measurements were used to construct a hypothetical velocity-depth profile and compare the results with the theoretical petrophysical measurements. Using brittleness analysis, the likely failure conditions were determined to be low-porosity dolomitic rock for fault nucleation in the Paleozoic basement.

Brittleness↗
DOE OSTI2024

Validating the rhenium proxy for rock organic carbon oxidation using weathering profiles

Chemical weathering over geological timescales acts as a source or sink of atmospheric carbon dioxide (CO 2 ), while influencing long-term redox cycling and atmospheric oxygen (O 2 ) at Earth's surface. There is a growing recognition that the oxidative weathering of rock organic carbon (OC petro ) can release more CO 2 than is locally drawn down by silicate weathering, and may vary due to changes in erosion and climate. The element rhenium (Re) has emerged as a proxy to track the oxidative weathering of OC petro , yet uncertainties in its application remain namely that we lack a systematic assessment of the comparative mobility of Re and OC petro during sedimentary rock weathering. Here we measure Re and OC petro loss across gradients in rock weathering at 9 global sites, spanning a range of initial OC petro values from ~0.2 % to >10 %. We use titanium to account for volume changes during weathering and assess Re and OC petro loss alongside major elements that reflect silicate (Na, Mg), carbonate (Ca, Mg) and sulfide (S) weathering. Across the dataset, Re loss is correlated with OC petro loss but not with loss of any other major element. Further, across the weathering profiles, the average molar ratio of OC petro to Re loss was 0.84 ± 0.15, with 8 out of 9 sites having a ratio >0.74. At one site (Marcellus Shale), the average ratio was lower at 0.58 ± 0.11. The excess loss of Re matches expectations that, typically, between ~0 and 20 % of the Re liberated by sedimentary rock weathering derives from silicate or sulfide phases, while some OC petro may be physically or chemically protected from weathering. Overall, our measurements provide validation for the Re proxy of OC petro oxidation and allow future work to further improve our knowledge of regional and global-scale rates of this important source of CO 2 in the geochemical carbon cycle.

58 GEOSCIENCES↗
DOE OSTI2026

Stimulated geologic hydrogen: from mechanistic control to engineered rock transformation

Geologic hydrogen (GeoH 2 ) generated from subsurface iron-rich rock–water reaction (i.e., serpentinization) is emerging as a promising candidate for the next primary energy source. Yet, accelerating GeoH 2 production from geological to human timescales via enhanced serpentinization remains a formidable scientific and technical challenge. Here, in this Review, we decipher the mechanistic control and explore strategies for accelerating in situ, engineered iron-rich rock transformation into carbon-free GeoH 2 by orders of magnitude. Serpentinization rate is hindered by low porosity and permeability of source rocks, suboptimal temperatures, unfavorable water chemistry, inefficient Fe 2+ -to-Fe 3+ conversion, thermodynamic constraints, and low reactive surface area. While closed-system experiments provide valuable mechanistic insights, open-system conditions with fluid circulation are more crucial for economically viable GeoH 2 production. We assess stimulation techniques from enhanced hydrocarbon and geothermal recovery as tools to be adopted or adapted for increasing reactive surface areas for stimulated GeoH2 production. We estimate that 7.40 × 10 5 to 1.73 × 10 6 million metric tons (Mt) of hydrogen could be engineered over 20 to 50 years from about 10% iron-rich rocks within 10 km depth of continental crust. Enabling GeoH 2 as a viable energy source requires not only advancing scientific frontiers but also forming a global GeoH 2 research network and innovation ecosystem to address the critical scientific, technical, societal, economic, and policy challenges.

08 HYDROGEN↗
DOE OSTI2025

Helium Release During Fracture and Granular Fragmentation of Rocks

Geogenic Helium-4 ( 4 He) in-situ increases locally in regions of large deformation generated naturally or anthropogenically. This gas release by deformation is a potential geochemical precursor signal for subsurface deformation. To evaluate the applicability of 4 He degassing for correlating deformation in different lithologies, we conducted high force crush tests, up to 97,800 N axial load, to assess the total 4 He released during fragmentation of the rocks. We observed that the highest 4 He released occurred in the sedimentary rocks and that release correlated strongly with lithologic age and U/Th content. Microstructural changes of the pre- and post-test rocks indicate that the degree of grain size reduction relates directly to the total 4 He released during crushing. The range of in-place 4 He was calculated based XRF measurements of uranium and thorium in each lithology, with the results indicating that the majority of the trapped 4 He was not released. However, the 4 He released by deformation depended upon how the each rock deformed during deformation and the degree of grain size reduction. We postulate that 4 He precursor signals can be used to understand subsurface deformation only if geomechanical and geochemical conditions for 4 He enrichment in a lithology are met.

Deformation signals↗
DOE OSTI2025

Temporal relationships among lunar crustal rocks

Temporal relationships among the three most common suites of lunar crustal rocks have been investigated by obtaining new high precision ages on Felsic/Alkali-suite Quartz monzodiorite Clast B from breccia 15405 and Magnesian-suite norite 78235/6/8/55/56 and comparing them to previously dated ferroan anorthosite sample 60025. The weighted average age of 4337.19 ± 0.49 Ma of 15405 Clast B is defined by zircon U-Pb and Pb-Pb ages as well as mineral isochron Sm-Nd and Nd-Nd ages. It is identical to the weighted average age for Apollo 17 norite 78235/6/8/55/56 of 4334.1 ± 3.5 Ma which is defined by Pb-Pb ages measured on baddeleyites in this investigation and less precise Pb-Pb and Sm-Nd ages reported in the literature. Both ages are ∼ 25 Ma younger than the weighted average of Sm-Nd and Pb-Pb ages reported in the literature on ferroan anorthosite 60025 of 4359.3 ± 2.3 Ma. The fact that ages of all three samples are defined by multiple U-Pb, Pb-Pb, Sm-Nd, and 142 Nd- 143 Nd chronometers provide confidence that they record the igneous crystallization history of the samples and do not represent disturbances or mixing lines with no temporal significance. Here, the extent to which these three ages represent broader scale magmatism is difficult to evaluate. Nevertheless, the age defined for 15405 Clast B, 78235/6/8/55/56, and 60025 are contemporaneous with the peak of ages observed in detrital zircons from the Apollo 12, 14, 15, and 17 landing sites (4340 ± 20 Ma), a Mg-suite Sm-Nd whole rock isochron defined by samples from Apollo 14, 15, 16, and 17 landing sites (4348 ± 25 Ma), and a Ferroan Anorthosite-suite Sm-Nd whole rock isochron defined by samples from the Apollo 15 and 16 landing sites (4354 ± 29 Ma). This implies that Ferroan Anorthosite-suite magmatism is temporally distinct and earlier than magmatism associated with the Mg-suite and the Felsic/Alkali-suite, as predicted by the lunar magma ocean model of lunar differentiation. The short 35 ± 10 Ma interval between primary ferroan anorthosite magmatism and secondary magmatism suggests that the lunar crust formed over a limited period of time. Although heat from decay of long-lived isotopes, large impacts, tidal heating associated with interactions between the Earth and Moon, and density driven overturn of the magma ocean have all been invoked to explain production of ancient secondary crustal magmatism, only tidal heating and cumulate overturn are consistent with the apparent short duration of secondary crustal magmatism and the great depth of crystallization implied for some Mg-suite samples. The initial ε 143 Nd values derived from the 15405 Clast B and 78238 Mg-suite norite isochrons, as well as a Mg-suite whole rock isochron are −0.23 ± 0.11, −0.27 ± 0.74, and −0.25 ± 0.09, respectively. They are identical within uncertainty indicating that Mg-suite and Felsic/Alkali-suite magmas were derived from materials that had the same time averaged Sm/Nd ratios since the formation of the solar system. This, combined with the contemporaneous nature of 15405 Clast B and 78235/6/8/55/56 Mg-suite norite, is consistent with evolution of both samples, and likely both magma suites, from a common source through closed system fractional crystallization or partial melting processes.

Chronology↗
DOE OSTI2025

Apatite geochemistry as a tool for understanding the petrogenesis of layered mafic-ultramafic rocks in the Bushveld Complex, South Africa

The sources of the magmas that formed the Rustenburg Layered Suite of the Bushveld Complex in South Africa remain debated, despite decades of research. Vertical and lateral variation in bulk rock and mineral separate Sr-Nd isotopic compositions, which generally indicate enriched sources, demonstrate that the layered sequence was formed by the emplacement of multiple batches of magma, crucially resulting in episodes of PGE-Cr-V mineralisation. The Lu-Hf isotope compositions of zircon are, however, at odds with the bulk rock Sr-Nd isotopic heterogeneity as they show near homogeneous compositions throughout the layered sequence (εHf (2.06 Ga) =−8). This lack of variation in Hf isotope composition has been attributed to deep, continental lithospheric mantle-related and/or crustal contamination of plume-derived Bushveld magmas. In this study, we analysed the major, trace element and Sr-Nd isotope geochemistry of apatite in the Rustenburg Layered Suite. Apatite occurs as an intercumulus mineral in the lowermost regions and a cumulus mineral in the uppermost regions of the layered sequence and can therefore be used to test existing models for the isotopic disequilibrium between bulk rock Sr-Nd and zircon Hf isotopic compositions. Apatite is largely chlorapatite in the lowermost regions and fluorapatite in the uppermost regions of the layered sequence. The Merensky Reef is unusual in that it comprises both chlorapatite and fluorapatite. Apatite throughout the layered sequence is generally unzoned and shows no evidence of late-stage alteration. Trace element data show that apatite is enriched in L/HREE, with common negative Eu-Sr anomalies. These trace element signatures are consistent with a magmatic origin for the apatite grains, with prior, or concurrent, plagioclase crystallization from the same melt. Variability in in situ Sr and Nd isotope compositions of apatite is recorded throughout the layered sequence with εNd (2.06 Ga) compositions varying between −2.5 and − 10.2 and initial 87 Sr/ 86 Sr compositions varying between 0.7079 and 0.7103 (for the Marikana dikes only). The variability in Sr-Nd isotope compositions of apatite is consistent with the bulk rock (and mineral separate) variation in Sr-Nd isotope compositions, suggesting apatite preserves primary magmatic compositions in the Rustenburg Layered Suite.

Apatite↗
DOE OSTI2024

Modeling nuclear waste disposal in crystalline rocks at the Forsmark and Olkiluoto repository sites – Evaluation of potential thermal–mechanical damage to repository excavations

We conduct coupled thermo-hydro-mechanical modeling of a KBS-3V repository design in crystalline rocks, using data and conditions from the Forsmark in Olkiluoto repository sites in Sweden and Finland. The study focuses on repository performance related to the impact of thermal and hydraulic evolution on the potential for thermal–mechanical damage to underground repository excavations. For the designs and conditions considered at the Forsmark and Olkiluoto repository sites, the simulations show a peak temperature well under the adopted performance target of a 100°C maximum temperature, whereas there is still a high potential for thermal–mechanical damage to the KBS-3V waste deposition holes. The thermal–mechanical damage is much more likely if rock permeability is so low that it delays saturation and swelling of bentonite-clay-based backfill beyond the time for the thermal–mechanical peak, which occurs 50 to 100 years after nuclear waste deposition. We also found that sidewalls of the KBS-3V emplacement tunnels are vulnerable to tensile fracturing due to the combined effect of thermal stressing and backfill swelling. The study highlights a strong interaction between bentonite-based backfill and host rock through capillary suction along with induced rock desaturation. A careful design and selection of the bentonite-clay-based backfill materials for KBS-3V tunnels and deposition holes can facilitate a timely saturation and backfill swelling that in turn can minimize thermal–mechanical damage.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Measurement of the Static Nonlinear Third-Order Elastic Moduli of Rocks: Problems and Applicability

The third-order elastic (TOE) model has been used to describe the widely observed nonlinear mechanical behaviors of earth materials. In addition to linear elastic constants ( λ , μ ), three nonlinear elastic moduli ( A , B , C ) are required for isotropic rocks. Contrary to previous research on dynamic TOE moduli, this study followed the protocol to measure strain and stress under uniaxial and hydrostatic compressive tests statically, which were later used to invert for the full set of TOE moduli for four standard rock types with differing pore structures of a porous oolitic limestone, quartz-rich sandstones, and a dense crystalline basalt. The applicability of the TOE model to characterize nonlinearity depends on the fulfillment of path-independence and small-strain assumptions. Using the measured static TOE moduli, the finite element model demonstrates that the stress in the vicinity of the wellbore is more amplified than the stress in the linear elastic case, which leads to a wider zone of rock failure around the wellbore. Due to the long-standing discrepancy between static and dynamic moduli, the rarely reported full set of static TOE moduli is necessary and will benefit future research in understanding the effect of rock nonlinearity on geophysical and geomechanical applications, such as long-term safe storage of CO 2 and generating process of geohazards.

58 GEOSCIENCES↗
DOE OSTI2026

Cosmogenic 3 He exposure dating in mafic rocks by “Virtual mineral separation” of pyroxene

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.

Geosciences↗
DOE OSTI2026

Seismic velocity modeling and earthquake relocations for the southern Nevada National Security Site (NNSS), with a focus on Rock Valley

The Rock Valley Direct Comparison (RV/DC) is the third phase of the Source Physics Experiment (SPE; Snelson et al., 2014), a project aimed at improving seismic discrimination between explosive sources and natural earthquakes. Phases I and II of SPE focused on studying the generation of shear waves in endmember geologies (granite and alluvium, respectively), while Phase III aims to detonate a chemical explosion co-located with the source region of the 1993 Rock Valley earthquake sequence, located in the Rock Valley Fault Zone (RVFZ), in the southern portion of the Nevada Test Site (now the Nevada National Security Site, or NNSS).

58 GEOSCIENCES↗
DOE OSTI2025

Nucleation-promoting and growth-limiting synthesis of disordered rock-salt Li-ion cathode materials

Disordered rock-salt oxides and oxyfluorides are promising positive electrode materials for high-performance lithium-ion batteries free of nickel and cobalt. However, conventional synthesis methods rely on post-synthesis pulverization to achieve cycling-appropriate particle sizes, offering limited control over particle microstructure and crystallinity. This accelerates degradation and complicates secondary particle processing. Here we present a synthesis strategy that enhances nucleation while suppressing particle growth and agglomeration across various disordered rock-salt compositions, including lithium–manganese–titanium oxide, lithium–manganese–niobium oxide, and lithium–nickel–titanium oxide systems. Applied to Li 1.2 Mn 0.4 Ti 0.4 O 2 , this method yields highly crystalline, well-dispersed sub-200 nm particles that form homogeneous electrode films with stable cycling behavior. Tested in cells with lithium metal as the counter electrode, these electrodes deliver ~200 mAh/g with 85% capacity retention relative to the first cycle after 100 cycles (20 mA/g, 1.5–4.8 V), and an average discharge voltage loss of 4.8 mV per cycle, compared to 38.6% retention and 7.5 mV loss per cycle for electrodes derived from pulverized solid-state particles. This approach suggests a route to enhance the performance and durability of disordered rock-salt electrodes for sustainable lithium-ion batteries.

Batteries↗
DOE OSTI2024

An integrated experimental–modeling approach to identify key processes for carbon mineralization in fractured mafic and ultramafic rocks

Abstract Controlling atmospheric warming requires immediate reduction of carbon dioxide (CO2) emissions, as well as the active removal and sequestration of CO2 from current point sources. One promising proposed strategy to reduce atmospheric CO2 levels is geologic carbon sequestration (GCS), where CO2 is injected into the subsurface and reacts with the formation to precipitate carbonate minerals. Rapid mineralization has recently been reported for field tests in mafic and ultramafic rocks. However, unlike saline aquifers and depleted oil and gas reservoirs historically considered for GCS, these formations can have extremely low porosities and permeabilities, limiting storage volumes and reactive mineral surfaces to the preexisting fracture network. As a result, coupling between geochemical interactions and the fracture network evolution is a critical component of long-term, sustainable carbon storage. In this paper, we summarize recent advances in integrating experimental and modeling approaches to determine the first-order processes for carbon mineralization in a fractured mafic/ultramafic rock system. We observe the critical role of fracture aperture, flow, and surface characteristics in controlling the quantity, identity, and morphology of secondary precipitates and present where the influence of these factors can be reflected in newly developed thermo-hydro-mechanical–chemical models. Our findings provide a roadmap for future work on carbon mineralization, as we present the most important system components and key challenges that we are overcoming to enable GCS in mafic and ultramafic rocks.

58 GEOSCIENCES↗
DOE OSTI2024

Determining Stress Orientation in Rock Valley, Nevada, Using Ambient Seismic Noise

The stress field and the mechanical properties of rocks are important to consider for nuclear explosion monitoring due to their effect on seismic wave radiation from earthquakes and explosions. At the Rock Valley Direct Comparison site, the regional orientation of the maximum horizontal compressive stress (SH max ) is well constrained, but it is unknown whether there are local heterogeneities. Here, I show that stress-induced anisotropy in nonlinear elasticity can be used to estimate the orientation of SH max . Rocks have compliant internal contacts, such as fractures and mineral grain boundaries, that respond to applied strains more strongly than individual mineral crystals. This strain response is asymmetric between compressive and dilatational strains and is affected by anisotropy in the ambient stress field. Traditional seismic velocity measurements are of strain-averaged velocity, which is less sensitive than nonlinear elasticity to the behavior of compliant internal contacts and therefore the stress field and fracture behavior. My results show that the orientation of SH max measured using ambient seismic noise aligns with regional estimates, with some heterogeneity that can also be explained by the limitations of the seismic array. I demonstrate the potential of a passive technique for monitoring the stress field in places that currently lack local measurements.

58 GEOSCIENCES↗
DOE OSTI2024

Life Cycle Analysis of Enhanced Rock Weathering

This study aims to evaluate the life cycle impacts of weathering agricultural fields in the U.S. Midwest with igneous rock (dunite/olivine) and industrial waste minerals using literature and estimated data from a techno-economic analysis (TEA). The GWP impact with igneous rock and industrial waste materials ranges 34–263 and 23–180 kg CO2-eq/tonne CO2 sequestered, respectively, with comminution and transportation as the main environmental hotspots for igneous rock, and transportation for industrial waste.

Izar-Tenorio, Jorge↗
DOE OSTI2026

DEMO-FTES: Development, Monitoring, and Control of Fracture Thermal Energy Storage in Crystalline Rock Formations (CRADA Final Report)

The DEMO-FTES project investigated the feasibility of Fracture Thermal Energy Storage (FTES) as a seasonal energy storage solution in crystalline rock formations. FTES leverages hydraulically induced fractures to exchange heat between circulating fluids and the surrounding rock mass, enabling long-term thermal energy retention due to the high specific heat and low thermal conductivity of rock. This approach has the potential to reduce heating and cooling energy demands and enhance building energy resilience. The project combined dimensional analysis, numerical modeling, laboratory experiments, and meso-scale field tests to evaluate FTES performance and advance its technology readiness level from 3 to 5. Scaling analysis identified key dimensionless parameters governing heat transfer and fluid flow, ensuring laboratory and field tests were representative of larger-scale systems. Numerical simulations using TOUGH and iTOUGH2 frameworks supported experiment design and interpretation, modeling fracture geometry, thermal-hydraulic behavior, and thermo-mechanical coupling. Laboratory tests at EPFL involved creating single and multiple fractures in 25 cm cubic samples of Gabbro and Granite under true triaxial stress.

25 ENERGY STORAGE↗
DOE OSTI2026

Modeling supercritical CO 2 flow and mineralization in reactive host rocks with PFLOTRAN v7.0

Understanding the flow and reactivity of CO 2 injected into geological reservoirs is important for many subsurface applications including secure geologic carbon storage (GCS), critical mineral extraction, enhanced geothermal systems (EGS), and enhanced oil recovery (EOR). Traditionally, subsurface CO 2 injection for GCS applications has focused on geologic formations with favorable subsurface configurations for CO 2 migration and trapping through non-reactive mechanisms such as structural, solubility, and petrophysical trapping. Recently, CO 2 -reactive rocks such as mafic and ultramafic basalts have been investigated for their potential to react with injected CO 2 in situ to simultaneously dissolve host rock minerals and mineralize CO 2 as carbonates. Engineering rapid CO 2 mineralization in the subsurface is attractive because of the increased density of stored CO 2 , the additional safety factors associated with solidification, and the potential to extract valuable critical minerals. Here we present recent developments in the parallel flow and reactive transport simulator PFLOTRAN to model coupled CO 2 -brine flow and reactive transport for a wide range of injection and production applications involving reactive CO 2 -brine systems. These developments are based on the well established and trusted CO 2 flow capabilities in the STOMP-CO 2 simulator. New capabilities added to PFLOTRAN include new CO 2 -brine equations of state with optional thermal coupling, several new constitutive relationships like capillary pressure smoothing and scanning path hysteresis, a fully implicit well model, and native linkage with PFLOTRAN's well-established reactive transport libraries. A series of benchmarks between PFLOTRAN and STOMP-CO 2 verify the newly developed CO 2 -brine flow capabilities. Demonstrations of coupled CO 2 -brine flow modeling and reactive transport show how CO 2 mineralization can be engineered in reactive host rocks. Finally, an example use case involving copper leaching by CO 2 and critical mineral extraction is presented to showcase the strengths of this new implementation. Several limitations still remain, including limited availability of field data to parameterize models. Future work should constrain the evolution of mineral surface area during mineralization and the temperature and/or pH dependence of geochemical reactions for specific systems of interest.

Critical Minerals↗
DOE OSTI2024

Effects of Interactions Between Produced Formation Fluid and Rock Matrix on Pore Structure of Caney Shale, Southern Oklahoma

ABSTRACT: Rock-fluid interactions change properties of shales during exploitation. To investigate effects of rock-fluid interactions on pore structure of shales matrix after hydraulic fracturing, powder samples from two late Mississippian Caney Shale cores in the Ardmore Basin, southern Oklahoma, were used to react with formation produced fluid from the field in the batch reactor analysis. X-ray diffraction for mineralogy and Low-pressure nitrogen adsorption isotherms for pore structure were measured for original, after-7days, and after-30days samples. Results show that the samples consist mainly of quartz, followed by clay minerals, carbonates, and feldspar. The pore sizes of micropore (<2 nm) and mesopore (2-50 nm) increase 14%-233% due to dissolution of pyrite, feldspar, and carbonates after 7 days. Due to the transformation from smectite to illite and the increase of pore size, the specific surface area (SSA) decreases after 7-days interactions. After 30-days interactions, the micropore volume slightly increases and the mesopore and macropore volume decreases. Due to the decrease of pore size, the SSA of 30-days reacted samples increases correspondingly and is lower (for the clay-rich sample) or higher (for the calcareous sample) than that of the unreacted samples. Findings improve our understanding of dynamic alteration of shale properties during production. 1. INTRODUCTION Energy demand will continuously grow owing to the increasing global population as well as energy consumption (EIA, 2023). On the other hand, shale gas and oil reshaped the energy market in the United States, enabling the United States to become a net-export of natural gas country in 2017 (EIA, 2023). However, shale reservoirs are challenging tight formations that are still poorly understood in the extraction and production of hydrocarbons (Ross and Bustin, 2009; Curtis et al., 2012; Xiong et al., 2015, 2021a; Li Y. et al., 2016; Gong et al., 2019a; Benge et al., 2021; Awejori et al., 2022; Huang et al., 2022). One of the most challenging topics is the rock-fluid interactions post hydraulic fracturing and its subsequent impacts on the pore structures of fractured formation matrix.

Xiong, Fengyang↗