Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reservoir characterization”

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.

608 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Subtask 1.5 – CO2 Injection Monitoring with an Optimized Scalable, Automated, Semipermanent Seismic Array

The scalable, automated, semipermanent seismic array (SASSA) method is a flexible and relatively cost-effective surface geophysical method for regular time-lapse monitoring of the movement of injected carbon dioxide (CO2) in a reservoir for CO2 enhanced oil recovery (EOR) or geologic CO2 storage operations. It has the advantages of a low-environmental-footprint while monitoring regions of a reservoir from the surface without the need for a regular grid distribution of receivers. Automated data collection is possible. As only time-lapse amplitude changes at the reservoir level due to CO2 movement within the reservoir are monitored, the turnaround time to deliver results from the SASSA method can be short, without the need for long, time-consuming data-processing workflows. As data is collected and processed, incremental information can be provided to the field operator. The Energy & Environmental Research Center (EERC) conducted a SASSA field test from September 2018 to November 2020 in a portion of the Bell Creek Field in Montana, which implemented new CO2 EOR field activities during the study period. Lessons learned from a proof-of-concept study were incorporated to improve the data quality of the SASSA method and demonstrate the viability of the technology. The EERC implemented several enhancements to improve data quality, including 1) an iterative survey design, which allowed placing the receivers in strategic locations where the movement of the CO2 in the reservoir could be tracked with minimum interference by the cultural noise in the study area; 2) the use of powerful seismic sources in the form of surface orbital vibrators, and 3) data acquisition during optimal periods. History-matched reservoir simulation was performed to predict gas saturation and pressure response induced by CO2 injection in the study area. The results were compared with the SASSA-measured responses to CO2 injection as a partial validation technique. A match between the two methods was observed for most of the SASSA points predicted to have intersected a CO2 saturation change. The validated results provide confidence that the SASSA method can be used independently as a CO2 saturation monitoring technique. As data are collected and processed, incremental information can be provided to the field operator. The critical components of the SASSA workflow for a successful application of the method are the following: Iterative survey design with information about CO2 injection activities from the oilfield operator. A detailed CO2 injection plan is the key driver to select the strategic monitoring location of the SASSA sensors. After this information is incorporated in the initial distribution of sources and receivers in the study area, high-resolution satellite images are used to identify ground locations not affected by cultural noise sources, such as power lines, pipelines/flow lines, or roadways. In the next iteration of the survey design, a scouting trip to the study area is needed to understand more details of the noise sources identified in the previous step and the intensity of the field activities that can also generate noise during the monitoring. Integrating the information from the scouting trip into the survey design to select the optimum source and receiver locations is the final step. Noise attenuation. The variety of noise types during seismic monitoring of an oil field is enormous. Tailored noise characterization and processing at a node-by-node level can enhance the performance and sensitivity of the SASSA technique. Future advancements that could improve the efficiency and application of the SASSA technology include: Gaining a better understanding of the noise field produced by the seismic source to aid the choice of receiver location. Surface noise from the source can overwhelm the small signal changes due to CO2 that the SASSA method measures. Improved data-processing workflow to automatically analyze and adapt to dynamic noise conditions associated with industrial settings. This subtask was funded through the EERC–DOE Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE- FE0024233.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Hydroelectric Power and Hydrogen Production Integration

Hydropower-based hydrogen production could introduce opportunities for new revenue streams for hydropower plants, including from energy storage and regeneration as well as from sale of the hydrogen product to external markets. Hydrogen-based energy storage and regeneration could also help support Idaho Power’s decarbonization goals by decreasing dependence on fossil-based peaking power plants. Additionally, integration of hydrogen production with hydropower generation could help address the issue of low dissolved oxygen river water conditions that commonly accompany hydropower plant operations by utilizing the oxygen byproduct from an electrolytic hydrogen production process as a resource for mitigation of low dissolved oxygen water conditions. Comprehensive techno-economic analysis of the hybrid hydroelectric and hydrogen energy storage system has revealed critical insights into the pathways and considerations for optimizing the economic value and environmental benefits of such systems. Among the three identified pathways of natural gas blending, regeneration, and direct sale of hydrogen, the direct sale of hydrogen emerges as the most profitable, particularly given the current pricing dynamics of hydrogen and electricity. However, as we anticipate a future grid characterized by higher renewable energy penetration, the landscape may evolve, featuring lower average electricity prices, heightened fluctuations, and more significant seasonal variations. Consequently, the attractiveness of electricity regeneration through hydrogen and the benefits of long-duration hydrogen storage are expected to increase substantially in such a dynamic energy scenario. The careful selection of component sizes within the hybrid system proves to be paramount for ensuring cost-effectiveness. Notably, the size of the hydrogen market has emerged as a critical determinant for the optimal size of the electrolyzer. Hydrogen storage should be sized to meet energy shifting requirements. The appropriate size of the fuel cell/microturbine generator hinges on the shape of electricity prices and the available revenue streams derived from participating in grid services. Striking the right balance among these components is essential for maximizing the overall efficiency and profitability of the hydrogen facility. Furthermore, the by-product of electrolysis, namely oxygen, introduces an additional dimension to the system's functionality. The oxygen generated can be effectively utilized for dissolved oxygen (DO) mitigation, particularly with larger electrolyzer sizes capable of satisfying the complete oxygen demand for this purpose. While the economic benefits derived from saved oxygen purchase costs may be relatively modest compared to other revenue streams, the environmental advantages of repurposing oxygen for DO mitigation could help Idaho Power meet their environmental obligations. In addition to the identified factors shaping the viability of hybrid hydrogen production and hydroelectric generation, it is noteworthy that the integration of hydrogen energy storage offers a unique advantage during unusually wet years. In such periods of increased water inflow, the hydrogen storage capacity serves as a valuable supplement to the reservoir. By utilizing hydrogen energy storage as a complementary reservoir, the system gains flexibility in reservoir management when dealing with fluctuations in water availability.

08 HYDROGEN↗
DOE OSTI2021

Core Analysis E-Reports Produced for the Southwest Regional Partnership on Carbon Sequestration for Wells 13-10A, 13-14, and 32-8 of the Farnsworth Unit, Texas

The submission is three electronic reports (e-reports) produced by Terra Tek, a former Schlumberger company, for and under the oversight of the Southwest Regional Partnership on Carbon Sequestration (SWP). The SWP designed and implemented a coring and core analysis program in conjunction with the previous field operator Chaparral Energy, L.L.C. The main SWP contacts who oversaw the coring and core analysis program, including geologic core descriptions and sample selection for testing by the Terra Tek, are listed below in the Authors section. The three reports are for Wells 13-10A, 13-14, and 32-8 of the Farnsworth Unit, TX, a site of the SWP's CO2 storage and enhanced oil recovery project. The e-reports include many data types that are not listed in full detail here, but may include: geomechanical, geochemical, and petrophysical measurements; white light and UV core photos; spectral gamma logs on the core; and Heterogeneous Rock Analysis logs. Fracture core review reports are included for Wells 13-14 and 32-8. The three e-reports each have many subfolders with the various types of data in spreadsheets, pdfs, image files, or other formats. The main e-reports for each well are listed by the well name. The sub-folder “E-Report” for each well has a Report.html file, which, when opened, will have hyperlinks to the various data sets and information. The files in the zipped folder will need to be extracted for the Report.html files to work properly. These authors thank Joseph Hall, a geologist formerly of Chaparral Energy, L.L.C., for assistance in designing the coring program, his geologic expertise of the Farnsworth Unit, and his assistance with characterizing the core. Christopher Gillespie, of Terra Tek at the time, managed the myriad geomechanical, geochemical, and petrophysical measurements provided by Terra Tek for this project. The core from which the samples were collected for this project are housed at the Subsurface Data and Core Libraries of the New Mexico Bureau of Geology and Mineral Resources, Socorro, New Mexico, USA. Funding for this project is provided by the U.S. Department of Energy's (DOE) National Energy Technology Laboratory (NETL) through the Southwest Regional Partnership on Carbon Sequestration (SWP) under Award No. DE-FC26-05NT42591. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. These e-reports describe objective technical results and analysis. Any subjective views or opinions that might be expressed in the reports do not necessarily represent the views of the US Department of Energy or the United States Government.

Caprock↗
DOE OSTI2023

Finding middle ground: Flow regimes designed for salmon and energy value

In regulated rivers, shaping seasonal flows to recover species at risk depends on understanding when to expect conflicts with competing water users and when their interests are aligned. Multi-objective optimization can be used to reveal such conflicts and commonalities. When species are involved, multi-objective optimization is challenged by the need to simulate complex species responses to flow regimes. Previously, we addressed that challenge by developing a simplified salmon model (Quantus) that defines cohorts of salmon by the river section and time in which they were spawned. Salmon in these space-time cohorts are tracked from the time redds (nests) are constructed until the cohort exits the tributary en route to the ocean. In this study, we modeled seasonal patterns in energy value and developed a Pareto-optimal frontier of seasonal flow patterns to maximize in-river salmon survival and hydropower value. Candidate flow regimes were characterized by two pulse flows varying in magnitude, timing, and duration and constrained by a total annual flow near the historical median. Our analysis revealed times when economic and salmon objectives were aligned and times when they differed. Pulse flows that favored higher energy value were timed to meet demand during extreme temperatures. Both salmon and hydropower objectives produced optimal flow regimes with pulse flows in early summer, but only solutions favoring hydropower value included high flows in mid-winter. Solutions favoring higher age-0 salmon survival provided an extended pulse flow in late winter/early spring, which suggests that access to productive floodplain habitat allowed faster growth and earlier out-migration and reduced the need for higher temperature-moderating flows later in spring. Minimum flows were also higher among solutions favoring salmon over energy. The tools used to produce these results can help to design simplified seasonal flow regimes by revealing compromise solutions that satisfy both fish and energy producers and highlighting when potential conflicts are likely.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Approaches for the Simulation of Coupled Processes in Evolving Fractured Porous Media Enabled by Exascale Computing

Models have historically represented fractured porous media with continuum descriptions that characterize the media using bulk parameters. The impact of small-scale features is not captured in these models, although they may be controlling the performance of subsurface applications. Pore-scale models can simulate processes in small-scale features by representing the pore space geometry explicitly but are computationally expensive for large domains. The alternative multiscale approach entails the combination of pore-scale and continuum-scale descriptions in a single framework. We use Chombo-Crunch, a computational capability that discretizes complex geometries with an adaptive, embedded boundary method to contrast these two approaches. Chombo-Crunch takes advantage of recent computational performance and memory bandwidth improvements resulting from the emergence of exascale computing resources. These combined improvements enable the efficient simulation of reactive transport in fractured media with a high degree of fidelity and the ability to capture the control small-scale processes exert on the overall medium evolution.

42 ENGINEERING↗
NASA NTRS1992

Venus steep-sided domes: Relationships between geological associations and possible petrogenetic models

Venus domes are characterized by steep sides, a circular shape, and a relatively flat summit area. In addition, they are orders of magnitude larger in volume and have a lower height/diameter ratio than terrestrial silicic lava domes. The morphology of the domes is consistent with formation by lava with a high apparent viscosity. Twenty percent of the domes are located in or near tessera (highly deformed highlands), while most other (62 percent) are located in and near coronae (circular deformational features thought to represent local mantle upwelling). These geological associations provide evidence for mechanisms of petrogenesis and several of these models are found to be plausible: remelting of basaltic or evolved crust, differentiation of basaltic melts, and volatile enhancement and eruption of basaltic foams. Hess and Head have shown that the full range of magma compositions existing on the Earth is plausible under various environmental conditions on Venus. Most of the Venera and Vego lander compostional data are consistent with tholeiitic basalt; however, evidence for evolved magmas was provided by Venera 8 data consistent with a quartz monzonite composition. Pieters et al. have examined the color of the Venus surface from Venera lander images and interpret the surface there to be oxidized. Preliminary modeling of dome growth has provided some interpretations of lava rheology. Viscosity values obtained from these models range from 10(exp 14) - 10(exp 17) pa*s, and the yield strength has been calculated to be between 10(exp 4) and 10(exp 6) Pa, consistent with terrestrial silicic rocks. The apparent high viscosity of the dome lavas suggests that the domes have a silicic composition or must augment their viscosity with increased visicularity or crystal content. Sixty-two percent of the Venus domes are associated with coronae, circular features that have been proposed as sites of mantle upwelling, and 20 percent of the domes are located near tessera, relatively high areas of complex deformed terrain. We have investigated several models that are consistent with these geologic associations. The first case involves the differentiation of basalt in a magma reservoir in the crust, perhaps produced by partial melting within a mantle plume. The second case is melting at the base of thickened basaltic crust, and the final case is volatile exsolution and enhancement within a basaltic magma reservoir. The association of domes with tessera might be explained by crustal remelting, while the association with coronae may be consistent with chemical differentiation of a magma reservoir or the exsolution and concentration of volatiles in the reservoir before eruption.

Pavri, B.↗
DOE OSTI2022

Pore-Scale Microenvironments Control Anthropogenic Carbon Mineralization Outcomes in Basalt

Thin sections and hand samples from 50 sidewall cores from the Wallula Basalt Pilot Demonstration, a basaltic carbon sequestration demonstration, provided the opportunity for the in-depth analysis of carbon mineralization induced by the injection of supercritical CO 2 . In this study, we used optical petrography and scanning electron microscopy to characterize the physical and chemical characteristics of the basalt components influenced by carbon mineralization reactions from all available hand samples and thin sections within the three CO 2 injection zones and caprock flow interiors. We found extensive carbonate mineralization, mostly in the form of nodules that were shown to be chemically zoned: Ca-dominant in the core regions and Ca-bearing Fe-dominant in the outer regions. Carbonate mineralization also took the form of fracture-filling carbonate cement, and acicular aragonite was also observed. Overall, we clarified the structural and paragenetic relationships between newly formed minerals, identifying a new fibro-palagonite-like, poorly crystalline silicate phase that grew on the carbonate nodules and pore-lining zeolites. Here we observed Fe-dominant carbonate precipitates surrounding acicular aragonite and rhombohedral Ca-carbonate cores, whereas previous studies of these zoned nodules did not observe these structures. A comprehensive accounting of the carbon mineralization products is vital to understand and predict the behavior of supercritical CO 2 in the subsurface given both the diversity of the host rock between and within injection zones, especially considering that the morphology and chemistry of the diverse precipitates are influenced by the pore-scale microenvironments of the basalt.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2011

Spatially Resolved, In Situ Carbon Isotope Analysis of Archean Organic Matter

Spatiotemporal variability in the carbon isotope composition of sedimentary organic matter (OM) preserves information about the evolution of the biosphere and of the exogenic carbon cycle as a whole. Primary compositions, and imprints of the post-depositional processes that obscure them, exist at the scale of individual sedimentary grains (mm to micron). Secondary ion mass spectrometry (SIMS) (1) enables analysis at these scales and in petrographic context, (2) permits morphological and compositional characterization of the analyte and associated minerals prior to isotopic analysis, and (3) reveals patterns of variability homogenized by bulk techniques. Here we present new methods for in situ organic carbon isotope analysis with sub-permil precision and spatial resolution to 1 micron using SIMS, as well as new data acquired from a suite of Archean rocks. Three analytical protocols were developed for the CAMECA ims1280 at WiscSIMS to analyze domains of varying size and carbon concentration. Average reproducibility (at 2SD) using a 6 micron spot size with two Faraday cup detectors was 0.4 %, and 0.8 % for analyses using 1 micron and 3 micron spot sizes with a Faraday cup (for C-12) and an electron multiplier (for C-13). Eight coals, two ambers, a shungite, and a graphite were evaluated for micron-scale isotopic heterogeneity, and LCNN anthracite (delta C-13 = -23.56 +/- 0.1 %, 2SD) was chosen as the working standard. Correlation between instrumental bias and H/C was observed and calibrated for each analytical session using organic materials with H/C between 0.1 and 1.5 (atomic), allowing a correction based upon a C-13H/C-13 measurement included in every analysis. Matrix effects of variable C/SiO2 were evaluated by measuring mm to sub-micron graphite domains in quartzite from Bogala mine, Sri Lanka. Apparent instrumental bias and C-12 count rate are correlated in this case, but this may be related to a crystal orientation effect in graphite. Analyses of amorphous Archean OM suggest that instrumental bias is consistent for 12C count rates as low as 10% relative to anthracite. Samples from the ABDP-9 (n=3; Mount McRae Shale, approximately 2.5 Ga), RHDH2a (n=2; Carrawine Dolomite and Jeerinah Fm, approximately 2.6 Ga), WRL1 (n=3; Wittenoom Fm, Marra Mamba Iron Formation, and Jeerinah Fm, approximately 2.6 Ga), and SV1 (n=1; Tumbiana Fm, approximately 2.7 Ga) drill cores, each previously analyzed for bulk organic carbon isotope composition, yielded 100 new, in situ data from Neoarchean sedimentary OM. In these samples, delta C-13 varies between -53.1 and -28.3 % and offsets between in situ and bulk compositions range from -8.3 to 18.8%. In some cases, isotopic composition and mode of occurrence (e.g. morphology and mineral associations) are statistically correlated, enabling the identification of distinct reservoirs of OM. Our results support previous evidence for gradients of oxidation with depth in Neoarchean environments driven by photosynthesis and methane metabolism. The relevance of these findings to questions of bio- and syngenicity as well as the alteration history of previously reported Archean OM will be discussed.

Williford, Kenneth H.↗
DOE OSTI2022

The weathered bedrock vadose zone: A hidden control on water availability in the western United States (Final Technical Report)

This report summarizes the research completed during the project entitled "The weathered bedrock vadose zone: A hidden control on water availability in the western United States" (Award: DE-SC0018039). The report covers work completed between 9/15/17 and 9/14/20. The objective of the exploratory project work was to evaluate the overarching hypothesis that the subsurface underlying actively eroding landscapes is systematically structured as a result of a balance between weathering and erosion and that this structure leads to dynamic moisture storage within weathered rock. To accomplish project objectives, we characterized hillslope weathering profile structure and conducted subsurface moisture monitoring across four study sites with comparable bedrock types but differing erosion rate and hydroclimate: The Angelo Coast Range Reserve (ACRR), Sagehorn-Russell Ranch (SRR), Antelope Valley Ranch (AVR), and the LBNL Watershed Function Scientific Focus Area (SFA). New characterization and monitoring infrastructure at the AVR and SFA sites was established through this project and enables long-term monitoring of hydrologic and carbon fluxes at the sites for comparison with long-term monitoring in other sites. Deep drilling, downhole logging, and characterization of chemistry and porosity of rock recovered from drilling reveals two weathering fronts across the sites. A shallower front is marked by the extent of pervasive fracturing and oxidation of pyrite and organic carbon. A deeper weathering front is marked by the appearance of open fractures and discoloration. At the AVR site, multiple ridge tops of different hillslope length were investigated to establish the relationship between the weathering front and topography. Drilling reveals that the fraction of the hillslope relief that is weathered scales with hillslope length, consistent with the theoretical framework that guided the project relating groundwater drainage with landscape evolution. With the establishment of new boreholes, we have documented significant water storage in weathered bedrock using low-field borehole nuclear magnetic resonance (NMR) across all four sites. In the upper 4 m of bedrock, the average water contents at ACRR, AVR, SRR, and SFA are approximately 25%, 16%, 25%, 23% respectively. At the ACRR site, where long-term monitoring via neutron probe surveys and time-domain transmission sensors are available, we have demonstrated that NMR reliably records moisture content and importantly, moisture content dynamics in weathered bedrock. Across all sites, our NMR monitoring has also revealed that a considerable fraction of the water storage occurs in fractures or pores that are significantly larger than the fine-grained bedrock matrix of the parent rock. To identify the proportion of water storage occurring in fractures and large pores, we analyzed the NMR signal by evaluating both the sum of echoes and an inversion of the NMR signal for a distribution of T2 relaxation times. Both methods reveal that at least 20% of the water storage occurs in fractures and larger pores across all sites. To constrain the interpretation of the in-situ field-scale NMR measurements, we compare the NMR response in the field to independent pore-structure information obtained via helium pycnometry. Results support the inference that dynamic seasonal water storage is dominantly restricted to the fracture network, which can reach 7% of the total volume. This dynamic range in water content is comparable to that of some soils, underscoring the need to understand how this region of the Earth System is structured, and functions as a hydrologic reservoir.

58 GEOSCIENCES↗
DOE OSTI2023

DEEPEN 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. Part of the DEEPEN project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology is based on the method proposed by Poux et al. (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., 2022. For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. More information on this process is described in the following sections. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty. There are two different versions of the Leapfrog model and associated favorability models: - v1.0: The first release in June 2023 - v2.1: The second release, with improvements made to the earthquake catalog (included additional identified events, removed duplicate events), to the temperature model (fixed a deep BHT), and to the index models (updated the seismicity-heat source index models for supercritical and EGS, and the resistivity-insulation index models for all three play types). Also uses the jet color map rather than the magma color map for improved interpretability. - v2.1.1: Updated to include v2.0 uncertainty results (see below for uncertainty model versions) There are two different versions of the associated uncertainty models: - v1.0: The first release in June 2023 - v2.0: The se...

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Ordovician-Cambrian Units: Hierarchical Evaluation of Geologic Carbon Storage Resource Estimates

The Indiana Geological and Water Survey (IGWS) led subtask 1.1 to assess the regional distribution and estimate the storage capacity of Ordovician-Cambrian stratigraphic units located within the partnership region. A comprehensive data set of wireline logs and petrophysical information was used to generate these interpretations. These data include core analysis for porosity and permeability, mercury injection capillary pressure (MICP), and existing well data including location and stratigraphic information. This report includes storage resource estimates (SREs) for three potential storage reservoirs(limestone and dolostone from the Upper Ordovician Trenton Limestone/Black River Group and equivalent units, the Middle Ordovician St. Peter Sandstone, and primary target reservoir rocks of the Lower Ordovician and Upper Cambrian Knox Supergroup and equivalent units) calculated using six methodologies: (1) a fixed value of porosity of 10 percent in all units evaluated; (2) a unique average porosity (per well) from wireline-derived porosity (neutron, sonic, and/or density porosity for each unit); (3) porosity values from core analysis; (4) a depth-dependent porosity model (Knox Supergroup only); (5) porosity based on a model based on petrophysical facies; and (6) SREs using National Energy Technology Laboratory’s CO2 Storage prospeCtive Resource Estimation Excel aNalysis (CO2-SCREEN beta V2). All methods used the same values for thickness for each unit. However, the areal extent of each assessment was limited by the data available for each method. Estimated volumes were calculated in 1-by-1 kilometer grid cells and summarized as county and total stratigraphic unit volumes. The resultant SREs mass are displayed using boxplots, which allow for comparing data statistics (mean values and variability) between methods. Differences observed in SRE results from the six methods are mainly attributable to differences in the data and conceptual models used to interpret or estimate porosity in each method. Based on this systematic variability between methods, it is inferred that methods 1, 4, and 6 are best used for regional-scale reconnaissance estimates of storage capacity while methods 2, 3, and 5 are more appropriate for local scales where more data is required. All estimates are data-density dependent and different methods require different amounts of data for reasonable assessments. ArcMap 10.5.1 software was used to portray SREs to help visualize spatial variance of estimates for each methodology, and more importantly, to highlight those areas having the greatest total storage potential estimates.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Structurally Frustrated Remeika Phases: Tolerance Factor and Non-Kramers’ Ion Effects

Quantum materials that are found on the verge of structural, magnetic, and electronic instabilities are deep reservoirs for exotic phenomena. Our criteria to study the reported Ln 5 Ru 6 Sn 18 (Ln = Gd, Tb) are guided by the tolerance factor for Remeika phases─a class of intermetallic compounds with the general formula A 3 M 4 X 13 (where A = rare-earth element, M = transition metal, and X = tetrel), which can be viewed as structural analogues of pseudoperovskites (4 × ABX 3 ), with one of the tetrel atoms occupying the A-site. We have grown single crystals, up to 0.5 cm on the largest facet, of Gd 5 Ru 6 Sn 18 (a = 13.8052(12) Å) and Tb 5 Ru 6 Sn 18 (a = 13.7825(9) Å), which we expect to be on the verge of a formation instability. Bulk magnetic, heat capacity, and electrical transport measurements reveal unusual behavior that originates from the trivalent Gd and Tb ions, with deviations from typical metallic behavior and magnetic ordering that is short-range or disrupted in some other way. Furthermore, to better understand these phenomena, single-crystal neutron diffraction is investigated, which reveals short-range correlations among the partially occupied Tb and Sn sublattices, with 3D-ΔPDF analysis supporting direction-dependent local order. Taken together, these measurements show the utility of using tolerance factors in identifying materials that are likely to exhibit complex phenomena and uncover environments where a formation instability may strongly impact emergent bulk magnetic and electronic phenomena.

36 MATERIALS SCIENCE↗
NASA NTRS2016

Magnesium and Titanium Isotopic Compositions of an Unusual Hibonite-Perovskite Refractory Inclusion from Allende: It Is Fun

Introduction: Hibonite-rich refractory inclusions are among the first solids that formed in the solar nebula, and thus provide constraints on the earliest environment in the Solar System. An unusual hibonite-perovskite inclusion from Allende, SHAL, consists of a large (approximately 500 by 200 microns) single hibonite crystal and coexisting blocky perovskite (approximately 200 microns in size). The hibonite is characterized by chemical and oxygen isotopic compositions similar to those in the FUN (Fractionated and Unknown Nuclear anomalies) inclusion HAL. However, the rare earth element (REE) patterns measured at different spots of SHAL hibonite are highly variable, ranging from Group II-like (light REEs enriched relative to heavy REEs) to Group III-like (relatively flat with slight Eu depletions), but overall contrast largely with that of HAL, especially in the Ce and Yb abundances. This implies that SHAL hibonite formed and underwent distillation processes under more reducing conditions. Interestingly, the accompanying perovskite has uniform, unfractionated oxygen isotopic compositions (averaging delta (sup 17) O equals delta (sup 18) O equals −7 per mille) and REE abundances that are completely different from those of SHAL hibonite. This has been interpreted that perovskite and hibonite may not be co-genetic. Here we performed Al-Mg and Ti isotopic measurements of SHAL hibonite and perovskite to determine if the FUN characteristics are observed in these two isotope systems, and to further constrain the origin and evolution of SHAL. Results: Isotopic measurements of Al-Mg and Ti in SHAL were performed on the UCLA CAMECA ims-1290 ion microprobe by following the analytical protocols described in [1]. The Al-Mg and Ti data obtained in both terrestrial standards and SHAL hibonite and perovskite are shown below. Both SHAL hibonite and perovskite, despite very high (sup 27) Al to (sup 24) Mg ratios, are devoid of (sup 26) Mg excesses that can be attributed to the decay of (sup 26) Al. Delta (sup 25) Mg (mass-dependent fractionation) in hibonite is approximately −5 per mille per atomic mass unit relative to Madagascar hibonite, but is not well constrained for perovskite due to very large uncertainties owing to extremely low Mg contents. Similar to Mg isotopes, SHAL hibonite and perovskite show essentially the same Ti isotopic compositions, with anomalies in (sup 50) Ti of approximately 14 per mille, but the former shows greater Ti isotope fractionation than the latter (2.5 per mille per atomic mass unit versus 0 per mille). Discussion and Conclusions: The Al-Mg and Ti isotopic compositions of SHAL hibonite are consistent with those of HAL, suggesting that SHAL hibonite is a FUN inclusion and likely formed prior to homogenization of (sup 26) Al and Ti isotope variations in the solar nebula. However, the formation mechanisms for SHAL and HAL differ, given the differences in the REE patterns and degrees of oxygen mass-dependent fractionation. The Group-II to Group-III like REE patterns, the Yb depletions, and negative delta (sup 25) Mg observed in SHAL hibonite are all consistent with condensation of the hibonite precursor in a reducing environment.. The lack of large Ce depletions in SHAL hibonite implies that distillation processes that fractionated hibonite's oxygen isotopes must have taken place under a reducing condition, but the extent to which SHAL hibonite was distilled appears to be less than HAL because of the smaller degree of oxygen mass-dependent fractionation. The perovskite shares essentially the same Ti and Mg isotopic compositions as hibonite and probably formed in the same reservoir.. The ultrarefractory REE pattern seen in perovskite likely resulted from gas-solid fractionation which depleted HREEs in this reservoir. This process also explains why SHAL hibonite is generally depleted in HREEs relative to LREEs.

Liu, M.-C.↗
NASA NTRSDate not supplied

Oxygen Isotopic Composition of Refractory Inclusions from the Miller Range (MIL) 090019 Carbonaceous Chondrite

Introduction: Primitive meteorites contain Calcium-Aluminum-rich Inclusions (CAIs)that preserve the records of the earliest times in the Solar System. CAIs record O isotopic variations with time and/or location from solar (16O-rich) to planetary (16O-poor) compositions [e.g., 1, 2].However, CAIs were also affected by parent body processing to varying degrees. The goal of this work was to identify the processes that resulted in the fine-scale spatial variations in oxygen isotopic compositions within the most ancient Solar System solids. We report oxygen isotopic composition of CAIs from MIL 090019 CO3.1 carbonaceous chondrite. MIL 090019 contains a high abundance of CAIs showing mineralogical as well as textural variations [3]. Methods: We analyzed 10 CAIs with different mineralogical assemblages, classified as corundum-bearing, grossite-bearing, hibonite-bearing, melilite-bearing, anorthite-bearing inclusions and amoeboid olivine aggregates(AOA’s). The mineralogical and petrological characterization of CAIs was performed using the JOEL Hyperprobe 8530 electron microprobe at NASA JSC. Oxygen isotopic imaging of the CAIs was done using the Cameca NanoSIMS 50L ion microprobe at NASA JSC. We followed the analytical protocol described in [4]. O-isotopic maps of the CAI were acquired by rastering a ~3 pA primary Cs+beam at 16 keVover an area of 20 ×20μm2over a period of ~7 hours. Negative secondary ions of 16O−, 17O−, 18O−, 28Si−, 24Mg16O−, 27Al16O−, and 40Ca16O−were simultaneously acquired using electron multiplier detectors at a mass resolving power sufficient to resolve the 16OH−interference from the 17O−peak, where the contribution of 16OH−was <0.1%.An electron flood gun was used to mitigate sample charging during the analyses. We used San Carlos olivine and Madagascar hibonite as isotopic standards to correct for the instrumental mass fractionation. The O-isotopic ratios were corrected for the quasi simultaneous arrival (QSA) effect and the detector dead time. All reported errors are 1 sigma. Results: The MIL 090019 carbonaceous chondrite hosts refractory inclusions varying in their mineralogies and textures[2, 3]. We found the oxygen isotopic compositions of mineral phases in MIL 090019 CAIs show large variations with Δ17O values varying from -27 to 0‰. A corundum-bearing inclusion records isotopic variations with corundum, melilite, perovskite, and anorthite showing Δ17Ovalues of -20 ±2.5, -10.2 ± 6.5, -10.7 ± 1.3,and -1.2 ± 3.2‰ respectively. A grossite-hibonite-bearing inclusion also records O-isotopic heterogeneity, where grossite, hibonite, Mg-rich spinel, and melilite show -3.7 ± 2.7, -22.9 ± 2.5, -21.4 ± 2.6, and -17.3 ± 2.8‰ respectively, whereas perovskite grains record a heterogenous isotopic composition ranging from -22.6 to -13.9‰. A hibonite-bearing inclusion shows average Δ17Ο values of -17.6 ± 5.3 ‰ and -15.1 ± 4.8 ‰ for hibonite and spinel respectively. The two perovskite grains also show relatively 16O-poor composition (Δ17Ο = -11.8 ± 6.4 ‰ and -8.6 ± 5.2 ‰). Anorthite has ~0.1 wt% Na2O and ~0.52 wt% FeO and a Δ17Ο value of -7.2 ± 5.0 ‰, whereas diopside records a Δ17Ο value of -3.5 ± 5.0 ‰. A hibonite-spinel-rich inclusion contains spinel and hibonite that are16O-rich (Δ17Οvalues -21.3 ± 2.8‰ and -19.4 ± 3.0‰, respectively). Melilite shows heterogenous oxygen isotopic composition, with an inner 16O-rich region and an outer 16O-poor region separated by a sharp boundary, with a bulk Δ17Ο= -16.2 ± 2.1‰, whereas two perovskite grains are relatively 16O-poor (Δ17Ο= -13.2 ± 3.7‰ and -7.1 ± 3.4‰).A spinel-rich inclusion contains perovskite grains with Δ17Οranging from -27.0 to -22.8 ‰and spinel with average Δ17Οvalue of -21.4 ± 2.8 ‰. In a melilite and spinel-rich CAI, perovskite, spinel, diopside, melilite, and Al-Ti-rich pyroxene record Δ17Ο values of-26.9 ± 4.6‰, -22.4 ± 4.3‰,-23.4 ± 4.9, -19.3 ± 3.4, and -19.0 ± 4.3 ‰ respectively. Olivine and spinel components of 2 AOA’s record a 16O-rich composition with Δ17O~-23‰. Discussion: The coexistence of 16O-rich and 16O-poor minerals within CAIs has been attributed to differing degrees of O isotopic exchange, reflecting a wide range of O diffusion rates[5]. However, the O isotopic heterogeneity observed in these CAIs cannot be explained by oxygen diffusion in nebular or parent body settings alone. The isotopic imaging of CAIs from the MIL 090019 meteorite shows variations in the O isotopic composition recorded in their mineral components, suggesting that these CAIs record the oxygen isotopic heterogeneity in the nebular gas from which they condensed. The coexistence of 16O-rich and 16O-poor mineral components in the CAI-forming region suggests that the gas in this region was not well-mixed and consisted of distinct 16O-rich and 16O-poor gaseous reservoirs.

P Mane↗