Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SYNTHETIC 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Applying 3D Geologic Modeling Workflows to the Argillite Reference Case (Rev. 1)

The objective of this short report is to document the application of our 3D geologic modeling workflow to an argillite (shale) host rock. Over the past four years, our team at Los Alamos National Laboratory has developed a geologic modeling workflow that can be applied to generic alluvial basins such as those found in the western United States. In “frontier” or “exploratory” basins where data are sparse, the first steps are to collect, evaluate and integrate available subsurface data into conceptual geologic models. Those models form the basis for constructing the geologic framework model, a 3D geocellular model ideally constrained by seismic and borehole data. To date we have constructed our models using “synthetic” well data derived from conceptual models, without the prospect of validating our workflow using “real” subsurface data. We were tasked to investigate whether our workflow designed for alluvial basin sediments could be applied to other potential repository host rocks. This task also provided the opportunity to work with high-quality subsurface data collected specifically for siting and evaluating a nuclear waste repository. Nagra, the Swiss governmental agency responsible for the disposal of the nation’s radioactive waste, generously provided us with data from two deep boreholes drilled through their argillaceous target formation. The aim of our proof-of-concept demonstration is to evaluate whether geostatistical methods offer a viable approach to property modeling in argillaceous rocks. Nagra provided us with the well data on the condition that we maintain confidentiality with all transferred information and results. Fortunately, Nagra posts numerous technical reports on its public website that describe the subsurface geology in great detail. All of the information and illustrations in this report related to the Swiss repository enterprise are taken from the Nagra public website.

58 GEOSCIENCES↗
DOE OSTI2024

Multivariate Bayesian Optimization of CoO Nanoparticles for CO 2 Hydrogenation Catalysis

The hydrogenation of CO 2 holds promise for transforming the production of renewable fuels and chemicals. However, the challenge lies in developing robust and selective catalysts for this process. Transition metal oxide catalysts, particularly cobalt oxide, have shown potential for CO 2 hydrogenation, with performance heavily reliant on crystal phase and morphology. Achieving precise control over these catalyst attributes through colloidal nanoparticle synthesis could pave the way for catalyst and process advancement. Yet, navigating the complexities of colloidal nanoparticle syntheses, governed by numerous input variables, poses a significant challenge in systematically controlling resultant catalyst features. We present a multivariate Bayesian optimization, coupled with a data-driven classifier, to map the synthetic design space for colloidal CoO nanoparticles and simultaneously optimize them for multiple catalytically relevant features within a target crystalline phase. The optimized experimental conditions yielded small, phase-pure rock salt CoO nanoparticles of uniform size and shape. These optimized nanoparticles were then supported on SiO 2 and assessed for thermocatalytic CO 2 hydrogenation against larger, polydisperse CoO nanoparticles on SiO 2 and a conventionally prepared catalyst. The optimized CoO/SiO 2 catalyst consistently exhibited higher activity and CH 4 selectivity (ca. 98%) across various pretreatment reduction temperatures as compared to the other catalysts. This remarkable performance was attributed to particle stability and consistent H* surface coverage, even after undergoing the highest temperature reduction, achieving a more stable catalytic species that resists sintering and carbon occlusion.

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

More Than 40 yr of Potentially Induced Seismicity Close to the San Andreas Fault in San Ardo, Central California

Evidence for fluid–injection–induced seismicity is rare in California hydrocarbon basins, despite widespread injection close to seismically active faults. We investigate a potential case of injection–induced earthquakes associated with San Ardo oilfield operations that began in the early 1950s. The largest potentially induced events occurred in 1955 (M L 5.2) and 1985 (M w 4.5) within ~6 km from the oilfield. We analyze Synthetic Aperture Radar interferometric images acquired by Sentinel–1A/B satellites between 2016 and 2020 and find surface deformation of up to 1.5 cm/yr, indicating pressure–imbalance in parts of the oilfield. Fluid injection in San Ardo is concentrated within highly permeable rocks directly above the granitic basement at a depth of ~800 m. Seismicity predominantly occurs along basement faults at 6–13 km depths. Seismicity and wastewater disposal wells are spatially correlated to the north of the oilfield. Temporal correlations are observed over more than 40 yr with correlation coefficients of up to 0.71 for seismicity within a 24 km distance from the oilfield. Such large distances have not previously been observed in California but are similar to the large spatial footprint of injection in Oklahoma. Furthermore, the San Ardo seismicity shows anomalous clustering with earthquakes consistently occurring at close spatial proximity but long interevent times. Similar clustering has previously been reported in California geothermal fields and may be indicative of seismicity driven by long–term, spatially persistent external forcing.

58 GEOSCIENCES↗
DOE OSTI2022

Experimental measurements of the viscosity and melt structure of alkali basalts at high pressure and temperature

Abstract Volcanic eruptions are shallow phenomena that represent the final stage of density- and viscosity- driven processes of melt migration from source rocks at upper mantle depths. In this experimental study, we investigated the effect of pressure (0.7–7.0 GPa) and temperature (1335–2000 °C) on the viscosity and the atomic melt structure of a synthetic anhydrous primitive alkaline basalt, an analogue of the pre-eruptive magma that likely feeds the Campi Flegrei Volcanic District at present day. Obtained viscosities (0.5–3.0 Pa s), mobility (0.1–0.4 g cm 3 Pa −1 s −1 ) and ascent velocity (1.5–6.0 m yr −1 ) are presented to support geochemical and geophysical observations of Campi Flegrei as a critical volcanic district currently undergoing gradual magma recharge at depth.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Impact of Concurrent Solubilization and Fines Migration on Fracture Aperture Growth in Shales during Acidized Brine Injection

Acidic hydraulic fracturing fluid chemically and physically alters shale rock fabric during injection and shut-in, creating a “reaction-altered zone” along the fracture faces. To better characterize the variable thickness and composition of this reaction-altered zone under advective flow, we take a coupled experimental and modeling approach. A fluidic cell, with six fiducial markers, is first fabricated to keep the rock sample in place during the core floods and to allow image alignment of acquired images. Then, we conduct a series of reactive core floods in a clay-rich siliceous Wolfcamp shale sample with 10 wt % carbonate, using a synthetic fracturing fluid under no confining stress and at room temperature. High-resolution computed tomography (CT) scans are periodically conducted to observe the spatial alteration of the fracture network. We then perform scanning electron microscopy (SEM-EDS) on the two orthogonal surfaces (fracture surface and freshly cut profile face) to generate high-resolution elemental maps that show the change in mineralogy, both with distance along a given flow path along a fracture surface and with depth from the fracture surface into the shale matrix. These results are contrasted against a two-dimensional (2D) advection-diffusion reaction model developed previously for batch reactions between shale and synthetic fracturing fluids. The model simulates the geochemical interaction occurring at the fracture/matrix interface and penetrating into the shale matrix during the reactive core flood. Both model and experimental results show that the acidic brine is neutralized during the core flood, corresponding to an increase in fracture aperture as a function of fluid volume injected with the greatest change near the inlet. SEM-EDS scans reveal significant dissolution of carbonates on the fracture surface without pyrite oxidation. The reactive transport model indicates that carbonate depletion into the shale interior should be observable, yet SEM-EDS shows no discernible loss of carbonate in the orthogonal profile face. The combination of these observations suggests an additional fracture evolution mechanism in the reactive system, i.e., fines migration. We show that fines migration enhances the access of fracturing fluid to the matrix resulting in a more pronounced fracture widening. Finally, we conclude that coupled mineral dissolution and fines migration govern fracture aperture growth during acidized brine injection. In this work, we effectively show the underlying risk of relying solely on models that do not include an important (transport) process that can alter the system significantly and propose a combined chemomechanical mechanism for fracture evolution appropriate for this shale mineralogy.

58 GEOSCIENCES↗
DOE OSTI2022

Geochemical Modeling of Celestite (SrSO 4 ) Precipitation and Reactive Transport in Shales

Celestite (SrSO 4 ) precipitation is a prevalent example of secondary sulfate mineral scaling issues in hydraulic fracturing systems, particularly in basins where large concentrations of naturally occurring strontium are present. Herein, we present a validated and flexible geochemical model capable of predicting celestite formation under such unconventional environments. Simulations were built using CrunchFlow and guided by experimental data derived from batch reactors. These data allowed the constraint of key kinetic and thermodynamic parameters for celestite precipitation under relevant synthetic hydraulic fracturing fluid conditions. Effects of ionic strength, saturation index, and the presence of additives were considered in the combined experimental and modeling construction. This geochemical model was then expanded into a more complex system where interactions between hydraulic fracturing fluids and shale rocks were allowed to occur subject to diffusive transport. We find that the carbonate content of a given shale and the presence of persulfate breaker in the system strongly impact the location and extent of celestite formation. The results of this study provide a novel multicomponent reactive transport model that may be used to guide future experimental design in the pursuit of celestite and other sulfate mineral scale mitigation under extreme conditions typical of hydraulic fracturing in shale formations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Seismic Imaging of the Mw 7.1 Ridgecrest Earthquake Rupture Zone From Data Recorded by Dense Linear Arrays

Abstract We analyze seismograms recorded by four arrays (B1–B4) with 100 m station spacing and apertures of 4–8 km that cross the surface rupture of the 2019 Mw 7.1 Ridgecrest earthquake. The arrays extend from B1 in the northwest to B4 in the southeast of the surface rupture. Delay times between P wave arrivals associated with ∼1,200 local earthquakes and four teleseismic events are used to estimate local velocity variations beneath the arrays. Both teleseismic and local P waves travel faster on the northeast than the southwest side of the fault beneath arrays B1 and B4, but the velocity contrast is less reliably resolved at arrays B2 and B3. We identify several 1–2 km wide low‐velocity zones with much slower inner cores that amplify S waveforms, inferred as damage zones, beneath each array. The damage zones at arrays B2 and B4 also generate fault‐zone head and trapped waves. An automated detector, based on peak ground velocities and durations of high‐amplitude waves, identifies candidate fault‐zone trapped waves (FZTWs) in a localized zone for ∼600 earthquakes at array B4. Synthetic waveform modeling of averaged FZTWs, generated by ∼30 events with high‐quality signals, indicates that the trapping structure at array B4 has a width of ∼300 m, depth of 3–5 km, S wave velocity reduction of ∼20% with respect to the surrounding rock, Q ‐value of ∼30, and S wave velocity contrast of ∼4% across the fault (faster on the northeast side). The results show complex fault‐zone internal structures (velocity contrasts and low‐velocity zones) that vary along fault strike.

Qiu, Hongrui↗
DOE OSTI2020

The transformation of diamond to graphite: Experiments reveal the presence of an intermediate linear carbon phase

Natural diamonds that have been partially replaced by graphite have been observed to occur in natural rocks. While the graphite-to-diamond phase transition has been extensively studied the opposite of this (diamond to graphite) remains poorly understood. We performed high-pressure and temperature hydrous and anhydrous experiments up to 1.0 GPa and 1300 °C using Amplex premium virgin synthetic diamonds (20–40 μm size) as the starting material mixed with Mg(OH) 2 as a source of H 2 O for the hydrous experiments. The experiments revealed that the diamond-to-graphite transformation at P = 1.0 GPa and T = 1300 °C was triggered by the presence of H 2 O and was accomplished through a three-stage process. Stage 1: diamond reacts with a supercritical H 2 O producing an intermediate 200–500 nm size “globular carbon” phase. This phase is a linear carbon chain; i.e. a polyyne or carbyne. Stage 2: the linear carbon chains are unstable and highly reactive, and they decompose by zigzagging and cross-linking to form sp 2 -hybridized structures. Stage 3: normal, disordered, and onion-like graphite is produced by the decomposition of the sp-hybridized carbon chains which are re-organized into sp 2 bonds. Our experiments show that there is no direct transformation from sp 3 C-bonds into sp 2 C-bonds. Our hydrous high-pressure and high-temperature experiments show that the diamond-to-graphite transformation requires an intermediate metastable phase of a linear hydrocarbon. This process also provides a simple mechanism for the substitution of other elements into the graphite structure (e.g. H, S, O).

36 MATERIALS SCIENCE↗
DOE OSTI2022

Deep Learning At Depth: Estimating subsurface parameters from geophysical monitoring data

Geophysical imaging techniques are a non-invasive way to image the subsurface and understand both subsurface solid (rock/soil) and fluid property distributions and their evolution in time. Inversions of the geophysical data, such as Electrical Resistance Tomography (ERT) data, are solved to estimate the subsurface property distributions, such as conductivity, and many inversion techniques smooth out sharp gradients in rock or fluid property distributions. Sharp gradients in subsurface properties tend to be present in situations with complex subsurface structures, which are common in many subsurface applications. We have successfully demonstrated that it is possible to inform, or constrain, inversions with neural networks trained on synthetic data with complex subsurface structures. Initial results suggest this process may be optimizable to yield property distributions that better represent the true property distributions than the same inversion process without the neural network constraint. Future work would optimize the neural network performance for this application and then apply the synthetic-data trained neural network to real data to understand the utility and performance of this technique for real data sets.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

RAPID UNSEPARATED RARE EARTH ELEMENT ANALYSES BY ISOTOPE DILUTION MULTICOLLECTOR INDUCTIVELY COUPLED MASS SPECTROMETRY (ID-MC-ICP-MS)

We present a method for analyzing rare earth elements (REEs) by isotope dilution without separation of individual elements. Utilizing a Nu Plasma 2 multicollector inductively coupled plasma mass spectrometer (MC-ICP-MS) equipped with an Aridus II desolvation nebulizer, sample introduction tuning parameters (torch position, gas flows) were investigated and used to minimize oxide formation (e.g., NdO/Nd <0.05%) and remove the need to analyze the REEs as separated “cuts.” Ratios for spike/sample isotopes in mixed synthetic standards containing varying amounts of Ce, Nd, Sm, Eu, Gd, Dy, Er, and Yb were measured accurately within 1% of IUPAC values. La and Lu in the standards were more than 2% off from expected values due to more demanding interference corrections. A USGS BCR-2 rock standard was also processed with flux fusion, iron coprecipitation and a single column pass bulk separation followed by analysis. The concentrations for Ce, Nd, Sm, Eu, Dy, Er, and Yb fell within 2% of consensus values. A United States Geological Survey (USGS) BIR-1A rock standard was also processed in triplicate by hydrofluoric and nitric acid digestion, yielding rare earth concentrations within 3% of consensus values with 0.3% relative standard deviation (n=3). This demonstrates the ability to produce high quality REE patterns of comparable quality to traditional methods involving labor and time intensive REE separation, with sample preparation and analytical times that are comparable to analyses against isotope dilution by single-collector ICP-MS.

Shen, Steve D.↗
DOE OSTI2023

Compaction of crushed salt for safe containment – a summary of the KOMPASS projects

Abstract. For the underground disposal of high-level nuclear waste in rock salt formations, the safety concept includes the backfilling of open cavities with crushed salt. For the prognosis of the sealing function of the backfill for the safe containment of the radioactive waste, it is crucial to have a comprehensive process understanding of the crushed-salt compaction behavior. The crushed-salt compaction process is influenced by internal properties (e.g., grain size, mineralogy, and moisture content) and boundary conditions (e.g., temperature, stress state, and compaction rate) and, therefore, involves several coupled thermal–hydro–mechanical (THM) processes (Hansen et al., 2014; Kröhn et al., 2017). With the paradigm shift from the limited release of radionuclides to safe containment due to the German Repository Site Selection Act passed in 2017, the importance of crushed salt as geotechnical barrier has increased, with a focus on the evolution of its hydraulic properties. Based on the knowledge gaps in the current process understanding, the “Compaction of crushed salt for safe containment” (KOMPASS) projects were initiated to improve the scientific basis behind using crushed salt for the long-term isolation of high-level nuclear waste within rock salt repositories. The efforts to improve the prediction of crushed-salt compaction begun during the first phase of the KOMPASS projects (Czaikowski et al., 2020) and were followed up in a second phase ending in June 2023. The primary achievements of the projects are as follows (Czaikowski et al., 2020; Friedenberg et al., 2022): specification of the KOMPASS reference material, an easily available and reproducible synthetic crushed-salt material, for generic investigations; development of pre-compaction methods and successful production of samples in the short term and under in situ loading conditions; formulation of an extended laboratory program addressing the isolated investigation of known relevant factors influencing the compaction behavior of crushed salt (Düsterloh et al., 2022); execution of long-term compaction tests addressing isotropic and deviatoric load changes, temperature, and compaction state; construction of a backfill body using the KOMPASS reference material in the Sondershausen mine through collaboration with the SAVER (Entwicklung eines salzgrusbasierten Versatzkonzepts unter der Option Rückholbarkeit) project (Schaarschmidt and Friedenberg, 2022); advancement of the tools for microstructure investigation methods (Svensson and Laurich, 2022); generation (first stages) of a microphysical process list combining literature research with our own findings; benchmarking of long-term compaction test for model development and optimization of various existing models as well as the development of new models; application of a virtual demonstrator (2D model representing a backfilled drift in rock salt) for the visualization of developments and the quantification of the models (Rabbel, 2022). In summary, the KOMPASS projects contributed to the reduction of uncertainties and the strengthening of the safety case for using crushed salt within rock salt repositories.

Friedenberg, Larissa↗
DOE OSTI2022

Solution-Grown Ternary Semiconductors: Nanostructuring and Stereoelectronic Lone Pair Distortions in I–V–VI 2 Materials

Alkali pnictogen dichalcogenides–I–V–VI 2 or APnCh 2 –have been identified as promising semiconducting materials for energy conversion devices. However, the controlled nanoscale synthesis and our understanding of the effects of cation ordering and stereochemically active lone pairs on the structures of these ternary compounds remain underdeveloped. Here, we use solution-phase chemistry to synthesize a family of APnCh 2 materials, including LiSbSe 2 , NaSbS 2 , NaSbSe 2 , NaBiS 2 , and NaBiSe 2 . Our approach utilizes alkali metal hydrides (AH) or carboxylates, A(O 2 CR), PnPh 3 , and elemental chalcogens as synthetic precursors and oleylamine or 1-octadecene as solvents. Synthetic manipulation via fine-tuning of reaction temperature enables control over the degree of ordering caused by the Sb 5s 2 lone pair-induced distortions in NaSbS 2 . Pair distribution function analysis demonstrates that the structure of the Sb-containing phases deviates much more from a disordered rock salt structure than that of the Bi-containing phases. This local distortion, induced by the Sb lone pair, leads to a previously unreported noncentrosymmetric NaSbS 2 crystal structure, which is additionally supported by second-harmonic generation measurements. Infrared and multinuclear solid-state NMR spectroscopies show that oleylamine or chelating carboxylates and, in some cases, unreacted precursors (LiH and PnPh 3 ) remain bound to the nanocrystalline surfaces. Further, a deeper understanding of the local atomic environment, long-range ordering, surface chemistry, and optoelectronic properties of these materials may speed up their fundamental study and application.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Synthetic Scientific Image Generation with VAE, GAN, and Diffusion Model Architectures

Generative AI (genAI) has emerged as a powerful tool for synthesizing diverse and complex image data, offering new possibilities for scientific imaging applications. This review presents a comprehensive comparative analysis of leading generative architectures, ranging from Variational Autoencoders (VAEs) to Generative Adversarial Networks (GANs) on through to Diffusion Models, in the context of scientific image synthesis. We examine each model's foundational principles, recent architectural advancements, and practical trade-offs. Our evaluation, conducted on domain-specific datasets including microCT scans of rocks and composite fibers, as well as high-resolution images of plant roots, integrates both quantitative metrics (SSIM, LPIPS, FID, CLIPScore) and expert-driven qualitative assessments. Results show that GANs, particularly StyleGAN, produce images with high perceptual quality and structural coherence. Diffusion-based models for inpainting and image variation, such as DALL-E 2, delivered high realism and semantic alignment but generally struggled in balancing visual fidelity with scientific accuracy. Importantly, our findings reveal limitations of standard quantitative metrics in capturing scientific relevance, underscoring the need for domain-expert validation. We conclude by discussing key challenges such as model interpretability, computational cost, and verification protocols, and discuss future directions where generative AI can drive innovation in data augmentation, simulation, and hypothesis generation in scientific research.

Generative Adversarial Networks↗
DOE OSTI2026

The solubility of ErPO 4 and Er speciation in hydrothermal fluids at varying pH and salinity between 350 and 450 °C

The rare earth elements (REE) are important for the green-energy transition and can be incorporated into the REE phosphates, such as xenotime-(Y), which also hosts heavy REE (Tb– Lu). Xenotime-(Y) is a common accessory mineral in metamorphic rocks and a range of mineral deposits where it controls the mobility of heavy REE, however, the impact of high temperature aqueous fluids on the behavior of heavy REE is largely unknown. Thermodynamic modeling can be utilized as a tool to predict the mobility of REE in hydrothermal aqueous fluids, but must be supported by accurate experimental data. Here, we measured the solubility of endmember synthetic xenotime-structured ErPO4 in NaCl-HCl-NaOH-bearing aqueous solutions at 350 °C and water vapor saturation pressure, at 400 and at 450 °C and 500 bar using batch-type Inconel reactors. Erbium speciation was investigated as a function of pH from 2.8 to 8, where Er chloride species are predominant at acidic conditions (pH <3) and Er hydroxyl complexes are predominant at near- neutral to alkaline conditions (pH >3). At pH 7–9, the measured ErPO 4 solubility (-9.8 to -7.5 log m Er ) is up to 2.5 orders of magnitude lower than thermodynamic predictions (-9.4 to -6.7 log m Er ) using existing thermodynamic databases. At pH 2–3, the predicted ErPO 4 solubility is ~0.5 orders of magnitude higher at 350 °C and ~1 order of magnitude lower at 450 °C compared to experimentally measured Er concentrations. The thermodynamic properties of aqueous Er species were therefore revised in this study. The partial molal Gibbs energy of formation (Δ f G 0 T,P ) for aqueous Er hydroxyl and chloride species are optimized using GEMSFITS and the logarithmic formation constants (logβ n ( Cl,OH) ) were derived at each experimental temperature and pressure. The updated thermodynamic properties for Er hydroxyl species (Er(OH) +2 , Er(OH) 2 + , and Er(OH) 3 0 ) show that their stability shifts to more acidic conditions at and below 400 °C. The Er chloride species (ErCl +2 and ErCl 2 + ) show increased stability compared to Er hydroxyl species at temperatures of 450 °C and 0.01 mol/kg NaCl. The updated thermodynamic properties are implemented into the GEM-Selektor modeling package to investigate the mobility of Er in saline hydrothermal fluids in equilibrium with alkaline rocks. Importantly, the updated properties for Er hydroxyl species result in low Er solubility at rock equilibrated pH conditions due to an expanded hydroxyl predominance zone, but lower aqueous complex stability overall, whereas previous models suggest greater stability for aqueous Er species. Furthermore, ErPO 4 solubility increases with decreasing temperature due to the deprotonation of HCl, which increases the acidity of hydrothermal fluids and the availability of Cl - to complex with the REE. These simulations highlight how fluid-rock reaction and temperature affect the mobility of REE in hydrothermal ore-forming systems.

58 GEOSCIENCES↗
DOE OSTI2023

Integration of seismic-pressure-petrophysics inversion of continuous active-seismic monitoring data for monitoring and quantifying CO 2 plume (Final Report)

The overall objective of this project is to develop and validate an integrated package of joint seismic-pressure-petrophysics inversion (jSPPI) of continuous active-source seismic monitoring dataset capable of providing real-time monitoring of CO 2 plume during geologic carbon sequestration (GCS). The three specific developments include: (a) the methodologies for fast seismic full waveform inversion of continuous active source seismic monitoring, (CASSM) datasets for simultaneously estimating velocity and attenuation, and with data assimilation; (b) joint Bayesian petrophysical inversion of seismic models and pressure data for providing and updating CO 2 saturation models; (c) the methods using multiple datasets including (Crainfield and Frio-II borehole) synthetic, laboratory, and field CASSM datasets. The outcomes of jSPPI include (a) a workflow for processing CASSM data, (b) Bayesian inversion algorithms using CASSM data and pressure response data, and (c) integration with data assimilation algorithms for continuously updating site-specific models used for prediction and reservoir management. The validation of joint FWI will be conducted using synthetic models based on the Cranfield and Frio experiments as well as field CASSM datasets collected as part of the Frio-II pilot injection. To quantify and map the mass and distribution of CO 2 (saturation), we will jointly invert velocity and attenuation measurements from the FWI with a Bayesian approach using a rock physics model for attenuation (e.g., White’s attenuation model with two selected patch sizes (White, 1976; Dutta and Seriff, 1979)). The Bayesian inversion will be applied to each time step in the CASSM survey in an updating scheme, which integrates with an ensemble of reservoir simulations at each step. A more complete experimental validation dataset will be collected as part of a mesoscale (2-3 m) gas-CO 2 injection experiment utilizing a higher frequency version of the CASSM system developed for laboratory studies; the integrated inversion will be demonstrated using this dataset which will provide both a dense geometry as well as more precise secondary confirmation measurements (e.g. saturation) typically not available in the field. The resulting real-time map of CO 2 saturation is able to provide a deeper scientific understanding of the complex, time-varying dynamics of subsurface fluid flow migration path as well as the rapid detection of CO 2 leakage hazards.

25 ENERGY STORAGE↗
DOE OSTI2021

Unnamed Pt(Cu 0.67 Sn 0.33 ) from the Bolshoy Khailyk River, Western Sayans, Russia, and a Review of Related Compounds and Solid Solutions

We describe a potentially new species of a platinum cupride–stannide mineral (PCSM) of composition Pt(Cu 0.67 Sn 0.33 ). It occurs in a placer deposit in the River Bolshoy Khailyk, southern Krasnoyarskiy kray, Russia. A synthetic equivalent of PCSM was obtained and characterized. The PCSM occurs as anhedral or subhedral grains up to 15 μm × 30 μm in association with various platinum-group minerals, Rh–Co-rich pentlandite and magnetite, all hosted by a placer grain of Cu–Au–Pt alloy. Synchrotron micro-Laue diffraction studies indicate that the PCSM mineral is tetragonal and belongs to the inferred space-group P4/mmm (#123). Its unit-cell parameters are a = 2.838 (3) Å, c = 3.650 (4) Å, and V = 29.40 (10) Å 3 , and Z = 1. The c:a ratio calculated from the unit-cell parameters is 1.286. These characteristics are in good agreement with those obtained for specimens of synthetic Pt(Cu 0.67 Sn 0.33 ). A review on related minerals and unnamed phases is provided to outline compositional variations and extents of solid solutions in the relevant systems PtNi–PtFe–PtCu, PdCu–PdHg–PdAu, PdHg–PtHg, and AuCu–PtCu. The PCSM-bearing mineralization appears to be related genetically with an ophiolitic source-rock of the Aktovrakskiy complex of the western Sayans. The unnamed phase likely crystallized from microvolumes of a highly fractionated melt rich in Cu and Sn.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Unveiling stimulation fluid-driven alterations in shale pore architecture through combined interpretation of TD-NMR and multi-component gas adsorption

The hydraulic fracturing stimulation process can alter petrophysical and flow properties of the shale matrix in unconventional reservoirs. These alterations can result from interactions with native and exogenous fluids, and can be detrimental to reservoir productivity. Different experimental procedures have been proposed to study reactive processes and thereby understand the aforementioned interactions. Most published works do not consider the rock initial saturation with formation fluid, or the effect of organic matter (OM). To gain insights into these effects, two reactive experimental sets, differing in experimental conditions, were conducted. In the first set, at 125°C and 45 MPa, a stimulation fluid surrogate was put in contact with Baxter Shale samples, previously saturated and equilibrated with synthetic formation fluid. A baseline consisted of exposure to formation fluid exclusively. In the second set, at 0.08 MPa and 25°C, contact with stimulation fluid took place as in the first set. Two additional assays in this second set either skipped the initial formation fluid saturation or removed the OM, respectively. Petrophysical properties analyses, focused on the pore architecture, relied on the combined interpretation of gas adsorption and Time-Domain Nuclear Magnetic Resonance (TD-NMR). Results show that the presence of the formation fluid significantly restricts pore space accessibility to the stimulation fluid and mitigates the pore-architecture alteration to the macro-pore region. Mineral dissolution and precipitation dominate over surface or wettability alteration. Next, results suggest that the OM restricts accessibility to some pores. Ultimately, the equilibration process with synthetic formation fluid drives mainly surface alteration, unhindered by OM.

58 GEOSCIENCES↗
DOE OSTI2024

Assessing the Value of Seismic Amplitude Versus Offset (AVO) Attributes for CO2 Storage Project Using a Bayesian Network Model for Decision Support

Attributes versus offset (AVO) are a set of measurements to analyze how the characteristics of reflected seismic waves change as a function of the offset. It can be useful for monitoring CO2 storage sites because the presence of leaked CO2 into the overlying aquifer can change the properties of the rocks and pore fluids composition that can alter the way seismic waves reflect and their amplitudes. The time-lapse changes in AVO attributes derived from repeat seismic surveys can help identify anomalies or shifts in the subsurface that could potentially be used as an indicator for CO2 leak detection. Our study leverages multiple seismic attributes derived from synthetic seismic data and Bayesian network model to quantify the probability of leak detection in the overlying aquifer above the storage reservoir. It helps to quantify the value of individual seismic attributes at multiple monitoring periods based upon their sensitivities.

Kumar, Abhash↗