SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
For presentation at the Geological Society of America Annual Meeting, Pittsburgh, PA, October 15-18, 2023.
For presentation at the Geological Society of America Annual Meeting, Pittsburgh, PA, October 15-18, 2023.
This study implemented a data-intensive supervised machine learning approach, through a quasi-experimental framework, with the objective of quantifying the impact of oil and gas operator-specific proprietary data on ML-based predictive model performance as compared to using oil and gas datasets that are publicly available.
This work was focused on developing novel methods to facilitate rapid data to knowledge transformation for subsurface applications, particularly those related to unconventional oil and gas.
This is the four-page proceedings paper for the ANS IHLRWM conference planned in Nov. 2025.
Fifty-eight rock chips from fifteen samples of sedimentary rocks from the Ramah Group (approximately 1.9 Ga) in northeastern Labrador, Canada, were analyzed for major and minor elements, including C and S, to elucidate weathering processes on the Earth's surface about 1.9 Ga ago. The samples come from the Rowsell Harbour, Reddick Bight, and Nullataktok Formations. Two rock series, graywackes-gray shales of the Rowsell Harbour, Reddick Bight and Nullataktok Formations, and black shales of the Nullataktok Formation, are distinguishable on the basis of lithology, mineralogy, and major and trace element chemistry. The black shales show lower concentrations than the graywackes-gray shales in TiO2 (0.3-0.7 wt% vs. 0.7-1.8 wt%), Al2O3 (9.5-20.1 wt% vs. 13.0-25.0 wt%), and sigma Fe (<1 wt% vs. 3.8-13.9 wt% as FeO). Contents of Zr, Th, U, Nb, Ce, Y, Rb, Y, Co, and Ni are also lower in the black shales. The source rocks for the Ramah Group sediments were probably Archean gneisses with compositions similar to those in Labrador and western Greenland. The major element chemistry of source rocks for the Ramah Group sedimentary rocks was estimated from the Al2O3/TiO2 ratios of the sedimentary rocks and the relationship between the major element contents (e.g., SiO2 wt%) and Al2O3/TiO2 ratios of the Archean gneisses. This approach is justified, because the Al/Ti ratios of shales generally retain their source rock values; however, the Zr/Al, Zr/Ti, and Cr/Ni ratios fractionate during the transport of sediments. The measured SiO2 contents of shales in the Ramah Group are generally higher than the estimated SiO2 contents of source rocks by approximately 5 wt%. This correction may also have to be applied when estimating average crustal compositions from shales. Two provenances were recognized for the Ramah Group sediments. Provenance I was comprised mostly of rocks of bimodal compositions, one with SiO2 contents approximately 45 wt% and the other approximately 65 wt%, and was the source for most sedimentary rocks of the Ramah Group, except for black shales of the Nullataktok Formation. The black shales were apparently derived from Provenance II that was comprised mostly of felsic rocks with SiO2 contents approximately 65 wt%. Comparing the compositions of the Ramah Group sedimentary rocks and their source rocks, we have recognized that several major elements, especially Ca and Mg, were lost almost entirely from the source rocks during weathering and sedimentation. Sodium and potassium were also leached almost entirely during the weathering of the source rocks. However, significant amounts of Na were added to the black shales and K to all the rock types during diagenesis and/or regional metamorphism. The intensity of weathering of source rocks for the Ramah Group sediments was much higher than that of typical Phanerozoic sediments, possibly because of a higher PCO2 in the Proterozoic atmosphere. Compared to the source rock values, the Fe3+/Ti ratios of many of the graywackes and gray shales of the Ramah Group are higher, the Fe2+/Ti ratios are lower, and the sigma Fe/Ti ratios are the same. Such characteristics of the Fe geochemistry indicate that these sedimentary rocks are comprised of soils formed by weathering of source rocks under an oxygen-rich atmosphere. The atmosphere about 1.9 Ga was, therefore, oxygen rich. Typical black shales of Phanerozoic age exhibit positive correlations between the organic C contents and the concentrations of S, U, and Mo, because these elements are enriched in oxygenated seawater and are removed from seawater by organic matter in sediments. However, such correlations are not found in the Ramah Group sediments. Black shales of the Ramah Group contain 1.7-2.8 wt% organic C, but are extremely depleted in sigma Fe (<1 wt% as FeO), S (<0.3 wt%), U (approximately l ppm), Mo (<5 ppm), Ni (<2 ppm), and Co (approximately 0 ppm). This lack of correlation, however, does not imply that the approximately 1.9 Ga atmosphere-ocean system was anoxic. Depletion of these elements from the Ramah Group sediments may have occurred during diagenesis.
We report mineral precipitation due to reactions with injected fluids during unconventional fracture stimulation is a well-recognized problem. The goal of this study is to evaluate secondary mineral precipitation and permeability attenuation under chemical injection scenarios specific to the Delaware basin. Whole cylindrical cores (2.54 cm diameter and 2.54 cm height) and ground shale (150–250 μm) from the carbonate-rich Bone Spring Formation, Delaware Basin TX (Leonardian), were reacted at 80 °C and 85 bar using a hydraulic fracturing fluid (HFF) recipe and an injection sequence typical of the Delaware Basin. The reacted shales and solutions were analyzed using a variety of laboratory- and synchrotron-based techniques to characterize both the chemical and spatial distributions of secondary mineral precipitation and identify changes in permeability and mineralogy. This carbonate-rich shale (>84% calcite) rapidly neutralized the acidic HFF. Synchrotron-based X-ray fluorescence mapping coupled with X-ray absorption spectroscopy (both bulk and micro) showed that most of the iron was in an oxidized form prior to exposure to HFF and that almost all iron(II) became fully oxidized after the reaction. Scanning electron microscopy images of the ground shale samples primarily identified iron(oxyhydr)oxide microcrystals on grain surfaces. A few small isolated iron-rich areas also contained sulfur, suggesting that some pyrite was preserved when isolated within a calcite crystal but that most was oxidized. The rapid neutralization of the acid spearhead in these carbonate-rich samples demonstrates that the acid spearhead is useful for initiating fractures in extremely calcite-rich rocks but does little to enhance rock permeability. This suggests that the impact of the acid spearhead is significantly smaller for carbonate-rich shales compared to clay-rich shales, which has broad implications for acidizing in carbonate-rich shale formations and iron transformations within these shales.
Montana State University’s (MSU) Energy Research Institute (ERI), in collaboration with the Center for Biofilm Engineering (CBE) and the Department of Civil Engineering (CE), has conducted a long‐term research program aimed at developing a novel cementing agent to address wellbore integrity and reduce the unwanted upward migration of fluids and greenhouse gases from the subsurface. The primary technology developed through this research program is known as ureolysis‐induced calcite precipitation (UICP), which harnesses bio‐chemical processes to precipitate calcium carbonate (CaCO 3 ). The same general process can also be called microbially-induced calcium carbonate precipitation (MICP) when microbes provide the process-catalyzing urease enzyme. Both terms are used in this report. Results have conclusively demonstrated that, if properly controlled, UICP can successfully seal fractures, high permeability zones, and compromised cement in the vicinity of wellbores and in nearby caprock. This technology has been successfully deployed to mitigate annular leakage in two test wells and over sixty commercial wells with a 100% success rate. This success in downhole deployment generates consideration of other subsurface applications where UICP could provide benefit to the energy sector, such as shale property modification for unconventional oil and gas recovery. The focus of this research project was to investigate fundamental material and mechanical properties of select shale cores and analyze how these properties change due to engineered mineral precipitation with the intent to control these properties to achieve a range of engineering objectives. Ultimately, the project aim was to identify valuable new areas where application of UICP might contribute to national energy security and environmental protection. The research workplan coupled UICP treatment of core samples, nuclear magnetic resonance (NMR) characterization, and mechanical strength testing at MSU with advanced X‐Ray micro-computed tomography (μCT) imaging and numerical modeling performed by collaborators at two national laboratories, the National Energy Technology Laboratory (NETL) and Lawrence Berkeley National Laboratory (LBNL). Experimental results are useful to inform geo-mechanical models which could be applied to predict mineralized rock formation behavior at field scale. Our findings suggest that NMR and μCT methods to detect and quantify biomineral formation in shale fractures are complementary and consistent with each other. Either could be used to estimate the volume of new mineral formed by UICP in shale fractures. The use of surfactants and guar gum to enhance biomineral precipitation in shale fractures merits further research. UICP can, under some conditions, increase the tensile strength of sealed shale fractures beyond that of the intact shale. These findings demonstrate that continued research in this area may be valuable to understanding and improving shale resource recovery techniques.
Throughout the Proterozoic Era, sedimentary organic carbon burial helped set the pace of global oxygenation and acted as a major modulator of atmospheric CO 2 and climate. Although Proterozoic rocks generally contain low concentrations of organic matter (OM), there are key exceptions to this rule, including the relatively OM-rich Arctic Bay shales from Baffin Island, Canada (Bylot Supergroup, Borden Basin, ~1.05 Ga). The mechanisms driving elevated OM concentrations in these and other Proterozoic shales remain poorly understood. In the Mesozoic and Cenozoic, organic matter sulfurization can be a major driver of enhanced OM burial across a range of redox conditions comparable to those inferred for many Proterozoic environments. Therefore, in this study, we evaluate the role of sulfurization in driving OM preservation in the Mesoproterozoic Borden Basin and discuss its relevance to Proterozoic systems in general. We present the first evidence for syngenetic-to-early-diagenetic OM sulfurization in a Proterozoic basin, which begins to fill a several-billion-year gap in our record of organic S across Earth history. We find that OM sulfurization was particularly extensive in shales from a relatively shallow-water section (Alpha River) but less extensive in shales deposited in deeper water (Shale Valley), which is consistent with models that infer sulfidic ‘wedges’ or O 2 -minimum-zone-type structures on shelf margins at least intermittently at this time. At the shallower site, organic S and pyrite are similarly 34 S-depleted and thus likely formed at roughly the same time near the sediment–water interface under conditions previously interpreted to have been ferruginous to intermittently sulfidic. In contrast, at the deeper-water site, large S-isotope differences between pyrite and organic S along with low apparent OM sulfurization intensities indicate that pyrite formation was favored over OM sulfurization during early sedimentation under variable but primarily ferruginous conditions. Although Mesoproterozoic biomass can be substantially sulfurized, indicators of sulfurization intensity are not correlated with OM concentrations, and therefore sulfurization does not appear to have been the primary driver of enhanced OM concentrations in Arctic Bay Formation shales. The link between sulfurization and total OM preservation may have been modulated during the deposition of Arctic Bay Formation shales by differences in iron availability, nutrient cycling, and particle dynamics in the Mesoproterozoic.
Background: Measuring the physical properties of shale is critical for optimizing engineering activities such as geothermal energy generation and hydraulic fracturing. Shale is a transversely isotropic material. Furthermore, this material can also include micro and macro cracks at different locations and orientations that cause it to behave differently under tensile or compressive loading. Objective: In this work, a combined experimental–numerical approach is proposed to evaluate the bi-modulus elastic properties of isotropic and transversely isotropic materials. Methods: Full-field strain measurements for a circular disk under diametral compression are integrated with a regression analysis technique to evaluate the elastic properties of bi-modulus materials subjected to tensile and compressive loads using two loading configurations on the same specimen. Digital Image Correlation (DIC) is used to measure the full-field strains. Subsequently, in the case of an isotropic material, a linear least-squares approach is utilized to process the experimentally determined strains in conjunction with analytical expressions of the stress fields (in terms of far-field loading) to determine the elastic modulus E, the shear modulus G, and the Poisson’s ratio $v$. In the case of a transversely isotopic material, such as shale, a finite element model is implemented to determine the stress fields (again in terms of far-field loading), which is followed by repeating the previous regression analysis in an iterative process to estimate the elastic parameters. Results: The results show that the proposed technique successfully provides a complete set of elastic properties as a function of both the loading condition and the principal material directions. The technique is validated by measurements on a known isotropic material and then applied to determine the properties of shale. Conclusion: In this work, the proposed approach is successfully used to calculate the bi-modulus elastic response of poly(methyl meth- acrylate) (PMMA) and shale. As expected, PMMA exhibits an isotropic response with no bi-modulus effect, however, shale exhibits both transverse isotropy and a bi-modulus effect. Therefore, this approach holds promise for investigating the elastic properties of materials like rocks and fiber-reinforced composite laminates as functions of the principal material directions and the loading conditions.
Deep subsurface stimulation processes often promote fluid-rock interactions that can lead to the formation of small colloidal particles that are suspected to migrate through the rock matrix, partially or fully clog pores and microfractures, and promote the mobilization of contaminants. Thus, the goal of this work is to understand the geochemical changes of the host rock in response to reservoir stimulation that promote the formation and migration of colloids. Two different carbonate-rich shales were exposed to different solution pHs (pH = 2 and 7). In this study, iron and other mineral transformations at the shale-fluid interface were first characterized by synchrotron-based XRF mapping. Then, colloids that were able to migrate from the shale into the bulk fluid were characterized by synchrotron-based extended X-ray absorption structure (EXAFS), scanning electron microscopy (SEM), and single-particle inductively coupled plasma time-of-flight mass spectrometry (sp-icpTOF-MS). When exposed to the pH = 2 solution, extensive mineral dissolution and secondary precipitation was observed; iron-(oxyhydr)oxide colloids colocated with silicates were observed by SEM at the fluid-shale interfaces, and the mobilization of chromium and nickel with these iron colloids into the bulk fluid was detected by sp-icpTOF-MS. Iron EXAFS spectra of the solution at the shale-fluid interface suggests the rapid (within minutes) formation of ferrihydrite-like nanoparticles. Thus, we demonstrate that the pH neutralization promotes the mobilization of existing silicate minerals and the rapid formation of new iron colloids. These Fe colloids have the potential to migrate through the shale matrix and mobilize other heavy metals (such as Cr and Ni, in this study) and impacting groundwater quality, as well produced waters from these hydraulic fracturing operations.
Catalytic pyrolysis of methane (CH4) is a promising approach to generate hydrogen (H2). The apparent activation energy of this process has a significant influence on the efficiency and required temperature for H2 generation. Preliminary experiments indicate that minerals present in shales have catalytic effects during in-situ H2 generation from shale reservoirs under electromagnetic (EM) heating. However, the quantitative role of such natural catalysts on apparent activation energy is not explored yet. This research evaluated the role of reservoir rock (shale) in EM heating and on the apparent activation energy of methane pyrolysis (MP) for in-situ H2 generation. Experiments are conducted in a customized EM reactor with frequency 2.45 GHz, and reaction temperature and generated gases are measured by a real-time IR pyrometer and gas analyzer, respectively. It is found that shale samples experienced thermal runway (TR) at 420 ºC under 0.15 kW EM power without any artificial heating promoter. In the presence of spent shale, CH4 conversion started approximately at 600 ºC with negligible amount of carbon dioxide (CO2) generated during this process. We further calculated apparent reaction order and apparent activation energy for MP process in the presence of shale under EM heating, which are 0.4357 and 98.12 kJ/mol, respectively. This research paves a way for leveraging the role of minerals as natural catalysts for enhancing H2 generation under EM heating in petroleum reservoirs.
This report provides an analysis of the clad barrier function associated with the direct disposal of dual purpose canisters (DPCs) under hypothetical conditions in a shale repository and in an alluvial repository, including the effect of a postulated criticality event inside a disposed DPC. Should a postulated criticality event occur in a hypothetical shale repository, cladding will primarily degrade by general corrosion. Stress corrosion cracking, hydride cracking, creep failure, pitting and crevice corrosion, rod pressurization, and clad unzipping are calculated to have little impact on cladding persistence. At the higher temperature expected during a postulated criticality event in a saturated shale repository, general corrosion of cladding would be rapid - on the order of 0.034 microns/yr. A few hundred years after onset of a postulated criticality event in a shale repository complete general corrosion of fuel assembly grid spacer walls and guide tubes will likely result in settling of fuel rods upon each other. This rod consolidation should displace the water moderator and possibly terminate a postulated criticality. The primary potential degradation pathway for cladding in a hypothetical alluvial repository is localized corrosion by fluoride, which cannot occur in a shale repository. Fluoride-enhanced corrosion of cladding would be accelerated under the slightly higher (< 100°C) temperatures associated with a postulated criticality event. The impact of criticality in both cases (shale and alluvial) would be to increase the amount of failed cladding. But it would require very specialized transport pathways.
A regional pattern of three closely spaced erosional events at and above the K/T boundary was determined from six Cretaceous/Tertiary boundary sections in the Negev of Israel. The sections were collected from locations throughout the central and northern Negev. All sections are lithologically similar. The Maastrichtian consists of a sequence of limestone beds intercalated with thin marly beds. In some sections, the last limestone bed is followed by 1 to 2 m of calcareous marls grading upwards into several meters of grey shale. In other sections the limestone bed is followed directly by grey shale with the contact containing particles of limestone and marl. A 5 to 20 cm thick dark grey organic-rich clay layer is present about 1.5 to 2.5 m above the base of the grey shale. The grey shale grades upwards into increasingly carbonate rich marls. No unconformities are apparent in field outcrops. During field collection the dark grey clay layer was believed to represent the K/T boundary clay. Microfossil analysis however identified the boundary at the base of the grey shale. The black shale represents a low productivity anoxic event similar to, but younger than, the K/T boundary clay in other K/T boundary sections. High resolution planktic foraminiferal and carbonate analysis of these sections (at 5 to 10 cm intervals) yield surprising results. The K/T boundary is marked by an erosional event which removed part or all of the uppermost Maastrichtian marls above the last limestone bed. Percent carbonate data for four Negev sections are illustrated and show the regional similarities in carbonate sedimentation. Faunal and carbonate data from the Negev sections thus show three closely spaced short erosional events at the K/T boundary and within the first 50,000 to 100,000 years of the Danian. These K/T boundary erosional events may represent global climatic or paleoceanographic events.
CO 2 injection is a promising method for enhanced oil recovery (EOR) in unconventional shale reservoirs. Here, in this work, we postulate that CO 2 EOR may be improved by the dissolution of surfactants into CO 2 . Although CO 2 is a relatively good solvent for oil, we show that CO 2 and Eagle Ford oil are immiscible at compositions above 70 wt % CO 2 , even at pressures as high as 62 MPa. The presence of a CO 2 –oil interface at reservoir conditions indicates that the addition of a surfactant has the potential to improve oil recovery–via wettability alteration from oil-wet to CO 2 -wet, CO 2 –oil interfacial tension (IFT) reduction, or both. Three nonionic surfactants (branched tridecyl ethoxylate Indorama SURFONIC TDA-9, branched nonylphenol ethoxylate Indorama SURFONIC N-100, and linear dodecyl ethoxylate Indorama SURFONIC L12-6) were evaluated for CO 2 -solubility, shale wettability alteration, effect on CO 2 –oil IFT, ability to generate CO 2 –oil foams, and ability to increase oil extraction from Eagle Ford, Mancos, and Bakken shale cores. Each surfactant dissolved in CO 2 up to 1 wt % at pressures and temperatures commensurate with CO 2 EOR. CO 2 -dissolved surfactants did not significantly affect CO 2 –oil IFT or generate CO2–oil foams, but they did induce a dramatic change in the contact angle of an oil droplet on an oil-aged shale chip in CO 2 from strongly oil-wet (11°) toward intermediate CO 2 –oil wettability (82°) (at 80 °C, 27.6 MPa). The branched tridecyl ethoxylated surfactant, SURFONIC TDA-9, afforded the highest oil recovery in core soaking experiments–75%, compared to 71% by pure CO 2 . Analysis of oil extracts by gas chromatography revealed that heavier oil components were produced when the surfactant was added to CO 2 . These results indicate that CO 2 -dissolved surfactants may increase oil recovery from shale by wettability alteration from oil-wet toward CO 2 -wet.
Mineral precipitation within hydraulically fractured shale may affect fluid flow pathways and impact longterm hydrocarbon production. The ability to predict geochemical reactions that lead to problematic mineral precipitation will lead to active reservoir management strategies to improve overall production. Using the Marcellus Shale as a case study, a combination of laboratory experiments and reaction path modeling was applied to determine which reactions are likely to occur upon introduction of hydraulic fracturing fluid into the shale reservoir. Experimental results indicate that contact between fracturing fluid and shale will result in dissolution of primary minerals (quartz, feldspars, kaolinite, chlorite, pyrite) and secondary mineral precipitation over time periods of less than one week. Precipitation of barite, Fe-oxides, feldspars, amorphous silica and clay is likely to occur within the reservoir during shut in and early flowback due to mixing between fracturing fluid and reservoir brine as based on modeling saturation indices using experimental fluid data. Reaction path modeling of the experimental scenarios corroborates the specific dissolution and precipitation reactions observed experimentally. Comparison of the results to injected and produced waters from a Marcellus Shale well pad in Greene County, PA, USA, shows that the mineral reactions occur during the hydraulic fracturing, shut in, and early flowback periods. The results presented in this paper demonstrate the value in applying experimental approaches to identify mineral precipitation/dissolution reactions that may significantly impact reservoir performance. The good agreement between geochemical models and experimental results provides confidence that numerical models can be applied to screen the potential fluid-mineral and fluid-mixing reactions in unconventional reservoirs that result in undesired mineral scale precipitation.
Image-based evaluation methods are a valuable tool for source rock characterization. The time and resources needed to obtain images has spurred development of machine-learning generative models to create synthetic images of pore structure and rock fabric from limited image data. While generative models have shown success, existing methods for generating 3D volumes from 2D training images are restricted to binary images and grayscale volume generation requires 3D training data. Shale characterization relies on 2D imaging techniques such as scanning electron microscopy (SEM), and grayscale values carry important information about porosity, kerogen content, and mineral composition of the shale. Here, we introduce RockFlow, a method based on generative flow models that creates grayscale volumes from 2D training data. We apply RockFlow to baseline binary micro-CT image volumes and compare performance to a previously proposed model. We also show the extension of our model to 2D grayscale data by generating grayscale image volumes from 2D SEM and dual modality nanoscale shale images. The results show that our method underestimates the porosity and surface area on the binary baseline datasets but is able to generate realistic grayscale image volumes for shales. With improved binary data preprocessing, we believe that our model is capable of generating synthetic porous media volumes for a very broad class of rocks from shale to carbonates to sandstone.
The Arctic Alaska petroleum province is geologically and geochemically complex. Mixed hydrocarbon charge from multiple source rocks and/or levels of thermal maturity is common within an individual oil pool. Biomarker and chemometric statistical analyses were used to correlate twenty-nine oils to five oil families derived from: (1) Triassic Shublik Formation (calcareous organofacies), (2) Triassic Shublik Formation (shaly organofacies), (3) Jurassic Kingak Shale, (4) Cretaceous shale (pebble shale unit and Hue Shale), and (5) Paleogene shale (Canning Formation). Age-diagnostic and source-related oil biomarker parameters establish clear genetic relationships between the normal oil-window components and their putative source designations. However, application of diamondoid analyses reveals mixed-oil accumulations with postmature charge contributions (diamondoid-rich and biomarker-poor) in many oils. Most sampled reservoirs contain a predominant charge derived from a single oil-window source plus a minor contribution from one or more higher maturity source(s). Variations in source organofacies also are recognized in the Shublik, Kingak, and Cretaceous oil families. In some cases, oils from multiple pools within a single field display relatively homogeneous geochemical profiles, suggesting a common source and migration pathway. For example, oil from the significant Pikka discovery is inferred to originate mainly from the calcareous Shublik Formation. In other cases, variability among oils from multiple pools within a single field (e.g., Milne Point, Colville River, and Northstar) likely indicates a more complex source, migration, and charge history. The results of this study may be useful for anticipating the composition of oil charge in stratigraphic traps with low-permeability sandstone reservoirs, where oil gravity and other chemical parameters may influence economic viability.