Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “drilling analysis”

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.

485 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Astrobiology Survey of a Lava Cave at Lava Beds National Monument by a Rover Carrying a Remote Sensing Instrument Payload

We report here on a survey of a lava tube cave by a rover that is instrumented for astrobiology missions. The NASA Ames testbed rover, CaveR, was deployed in Valentine Cave in Lava Beds National Monument (N. CA, USA) during August of 2018. The rover carried an instrument package consisting of Near Infrared and Visible Spectrometer System (NIRVSS) a point spectrometer operating in 1590-3400 nm range, sensitive to H2O and -OH bearing minerals, pyroxenes, and carbonates (Roush, et al 2018); the bore sighted Drill Operations Camera (DOC), a monochrome imager illuminated by LEDs at 410, 540, 640, 740, 905 and 940 nm; a Realsense ™ depth sensor system for 3D model generation; and a high resolution DSLR stereo camera. The payload was mounted on a tiltable instrument platform attached to the left side of the rover. The rover was driven manually in the cave by field operators, following instructions from a remote science operations team, and simulating a mission concept with science-guided autonomy. A simulated mission took place for 3 days with a team of 3 scientists selecting targets and interpreting data from the payload. To begin the mission, the rover drove along one wall of the cave imaging continuously with the Realsense in 20 m cave segments, three total. At the start of each day, the images from a 20m segment and a panorama stitched from them were provided to the science team to examine. The science team used these data to prioritize specific points along the cave wall for the collection of NIRVSS, DOC, and DSLR data. The objective of the data collection was to identify and study putative biological and mineralogical features in the cave. The data were delivered in xGDS, a customized mapping, planning, and data base management software developed at NASA Ames (Lee, et al 2013). Once the targets for further observations were selected, a plan for collecting the observations (positions in the cave and pointing for each requested observation) was constructed using xGDS and delivered to a rover team to execute the science data collection plan. Acquired data were delivered back to the science team for analysis. Preliminary results from the experiment illustrate the utility of the system (rover plus payload) to study the cave geology and mineralogy and its potential for identifying biomineral features.

Wong, Uland↗
DOE OSTI2021

Microstructure and water absorption of ancient concrete from Pompeii: An integrated synchrotron microtomography and neutron radiography characterization

There is renewed interest in using advanced techniques to characterize ancient Roman concrete due to its exceptional durability and low-carbon footprint. In the present work, samples were drilled from the “Hospitium” in Pompeii and were analyzed by synchrotron microtomography (μCT) and neutron radiography to study how the microstructure, including the presence of induced cracks, affects their water adsorption. The water distribution and absorptivity were quantified by neutron radiography. The 3D crack propagation, pore size distribution and orientation, tortuosity, and connectivity were analyzed from μCT results using advanced imaging methods. Porosity was also measured by mercury intrusion porosimetry (MIP) as a reference. Ductile fracture patterns were observed once cracks were introduced. Compared to Portland cement mortar/concrete, the Pompeii samples had relatively high porosity, low connectivity, and a similar coefficient of capillary penetration. In addition, permeability was predicted from models based on percolation theory and pore structure data to evaluate the fluid transport properties. Understanding the microstructure of ancient Pompeii concrete is important because it could inspire the development of modern concrete with high durability.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multiscale Characterization of the Caney Shale — An Emerging Play in Oklahoma

From a hydrocarbon perspective, the Caney Shale has historically been evaluated as a sealing unit, which resulted in limited studies characterizing the rock properties of the Caney Shale and its suitability for hydraulic fracturing. The objective of our research is to help bridge the current knowledge gap through the integration of multiscale laboratory techniques and to characterize the macro- and microscale rock properties of the Caney Shale. We employed an integrated approach for the characterization of the Caney using 200 ft (61 m) of Caney core from a target well in southern Oklahoma. Core observation and petrographic analysis of thin sections were combined to characterize the general rock types and associated fabrics and textures. Mineralogical composition, pore system architecture, and rock fabric were analyzed using x-ray diffraction (XRD), scanning electron microscopy/energy dispersive x-ray spectroscopy (SEM/EDS), and focused ion beam (FIB)-SEM. In addition, rebound hardness and indentation testing were carried out to determine rock hardness (brittleness) and elasticity, respectively. With the integrated multiscale characterization, three mixed carbonate-siliciclastic rock types were identified — mudstone, calcareous siltstone, and silty carbonate — likely representing a spectrum of deposition from low to relatively high energy environments in the distal portions of a ramp system. Silty carbonate contains mostly interparticle pores. The calcareous siltstones and silty mudstones contain a combination of organic matter pores and interparticle pores. Each of the rock types shows unique mineralogical compositions based on XRD. The mudstone lithofacies has the highest clay content and the least carbonate content. Calcareous siltstones show moderate carbonate and clay content. Silty carbonate indicates the highest carbonate content with the least clay content. In an order of mudstone, calcareous siltstone, and silty carbonate, rebound hardness and Young’s modulus show an increasing trend. As a result of rock-fluid interactions, there are potential scaling reactions during completion and production that could ultimately affect permeability and production rates. Overall, the proposed multiscale integration approach is critical for the geologic characterization of most rocks. However, in shale reservoirs dominated by microporosity and microstructure where engineered fractures are expected to provide permeability at a reservoir scale, successful integration is essential. An optimized, integrated geological characterization of the Caney Shale that is well aligned with the engineering designs in drilling, completing, and producing wellbores will ultimately lead to optimal production while providing safe and environmentally responsible operations.

Wang, Yulun↗
DOE OSTI2022

Physics-Informed and Data-Driven Prediction of Residual Stress in Three-Dimensional Machining

Efficient and reliable prediction of machining-induced residual stress (RS) is a key requirement for truly integrated computational materials engineering (ICME). Currently available process modeling approaches, including empirical, analytical, and numerical methodologies lack predictive power and require substantial calibration and validation data. Moreover, most model-based approaches consider only two-dimensional (2D) (i.e., orthogonal), cutting processes. Meanwhile, industrial processes such as milling, turning, and drilling are inherently three-dimensional (3D). The present work attempts to bridge the gap between 2D and 3D through careful consideration of the process physics, including geometric, kinematic, and size-effect constraints to realize robust prediction of how RS develops in 3D machining. Using a novel in-situ experimental technique and digital image correlation (DIC) to determine equivalent Hertzian contact widths, contact pressures, and friction coefficients, the proposed methodology leverages a discretized conversion algorithm that includes multi-pass shakedown effects. This paper presents a semi-analytical model to predict machining-induced RS in 3D turning operations, which are used representatively for 3D processes more generally. Rather than follow a ‘brute force’ 3D FEM approach or conduct countless experiments to train a purely data-driven machine learning algorithm, the proposed approach builds on previous 2D modeling work. Through careful consideration of the process physics, including complex geometry/kinematic considerations of 3D turning, the authors demonstrated an experimentally calibrated approach, as well as validation based on published RS data. Model predictions and previously published measurement data of RS depth profiles for turning of Inconel 718 were compared for a range of process parameters. Correlation between the proposed 3D model and validation data was found to be within the margin of experimental error for most conditions. The proposed model appears to capture the overall behavior of 3D RS depth profiles with acceptable accuracy, particularly the key metrics of near-surface stress, peak stress magnitude and location, as well as overall stress profile depth. This report presents a physics-informed, data-driven approach for efficient calibration of a 2D model for machining-induced RS through DIC analysis of in-situ characterized subsurface displacement fields.

42 ENGINEERING↗
DOE OSTI2023

Methane Partial Oxidation over Multifunctional 2-D Materials

The objective of this research is to design, synthesize, and evaluate highly selective, active, and stable multifunctional catalysts for the low temperature (< 500 Kelvin (K)) partial oxidation of methane to methanol (MTM) with molecular oxygen: CH 4 (g) + $\frac{1}{2}$O 2 (g) → CH 3 OH(g). Methane, the primary component of natural gas, is a source of energy and economic growth as well as an environmental concern. Recent developments in horizontal drilling and enhanced extraction methods have resulted in production of an estimated 62.4 trillion m 3 of ‘stranded’, or uneconomic, natural gas. Uneconomical natural gas is often flared or vented at remote oil production sites. Leaked, flared, and/or vented gas represents a "lost opportunity”, and this research project aims to maximize the value of the resource. Conventional catalysts for MTM suffer from low methanol selectivity since they exhibit ~0.55 eV higher barrier for C-H bond activation of methane compared to methanol. Without breaking these scaling relations, methanol oxidation is orders of magnitude faster than methane oxidation and it is very challenging to envision a process with economically viable single-pass yield. Here, we chose to investigate single-atom catalysts embedded and stabilized in two-dimensional materials such as graphene (GR) and "supported" on Group VIII and IB transition metals such as nickel. The electronic atomic monolayer-metal support interaction (EAMSI) present in these systems could promote methanol selectivity by breaking the scaling relations of the C-H bond activation of methane and methanol. A density functional theory (DFT) based computational study focused on predicting families of GR-based catalysts that could be active and selective for MTM. The catalyst systems predicted by the computational study were synthesized and evaluated for the gas phase MTM under relevant conditions. Unfortunately, the experimental activity and selectivity was lower than computationally predicted. The origin for the discrepancy is likely related to difficulties in synthesizing single atom catalysts in a threecomponent catalyst system at high density and with high selectivity. Future work in our groups is thus focused on reducing the system complexity to a two-component catalyst system. Finally, a techno-economic analysis (TEA) was also conducted to identify critical bottlenecks that inhibit future commercialization.

03 NATURAL GAS↗
DOE OSTI2021

Evaluation of “Tight Oil” Well Performance and Completion Practices in the Powder River Basin - “Time Slice” Analysis

The objective of this paper is to assess how well completion practices and well performance have evolved over time (“time slices”) in the tight sands and shales of the Powder River Basin (PRB). This information can provide a foundation for helping operators define more effective well completion practices and thus optimize well performance in an emerging “tight oil” basin. To start, the authors assembled a “database” containing well completion practices, production data, and geologic information for more than 800 horizontal (Hz) wells targeting three “tight oil” formations -- Turner, Frontier, and Mowry--placed on production in the past dozen years. To control for the impact of geologic and reservoir properties on well performance, the authors defined 12 geologically distinct areas (“partitions”) for the Turner, Frontier, and Mowry Shale in the Powder River Basin (4 partitions in each formation). This paper discusses the methodology for establishing “time slices” for each partition within each of the three “tight oil” formations, including (1) taking out, to the extent practical, the effects of geology by partitioning the three “tight oil” formations based on their geologic parameters; (2) using type curves and estimated ultimate recoveries (EURs) to establish a reliable measure of well performance; and (3) rigorously inspecting the production and well completion data to assure a quality dataset. Partition #2, a high thermal maturity area of the Frontier Sandstone in the western Powder River Basin, helps illustrate the results of the study. Using a dataset of 55 Hz wells, the study found that well performance has steadily improved, with oil EURs increasing from 180 MBbl in 2012-13 to 290 MBbl in 2018-19. Much of this improvement was due to the use of longer Hz laterals, increasing notably from 3,765 ft in 2012-13 to 10,260 ft in 2018-19. More intensive completion practices contributed, as well. For example, the number of frac stages more than doubled, from 17 to 35, and proppant concentrations increased from 1,060 lbs/ft to 1,320 lbs/ft. Finally, the data show diminishing returns to longer laterals and more intensive completion practices. For example, the key performance measure of oil EUR per 1,000 ft of lateral is decreasing. For 2012-13, Hz wells recovered 50 MBbl of oil per 1,000 ft of lateral while the more recent 2018-19 Hz wells only provide 30 MBbl of oil per 1,000 ft of lateral. Considerable insight can be gained by using “time slices” and geologic partitioning to better understand the relationship between changes in well drilling and completion practices and changes in well performance in emerging “tight oil” plays. The results from this study can also serve as a foundation for subsequent, more intense efforts involving data analytics for defining more optimum well completion practices targeting specified geologic settings and formations.

02 PETROLEUM↗
NASA NTRS1989

Non-LTE analysis of extremely helium-rich stars. The hot sdO stars LSE 153, 259 and 263

Results of a non-LTE fine analysis based mainly on high-resolution CASPEC spectra for three extremely helium-rich sdO stars are discussed in order to explain hydrogen deficiency in single stars. High temperature (Teff = 70,000 to 75,000 K) and a position in the log Teff - log g diagram were found close to the Eddington limit. Various abundance estimates are derived for hydrogen (upper limits only), carbon, nitrogen, and magnesium. Hydrogen is reduced to less than 10 percent by number in LSE 153 and LSE 263, and to less than 5 percent in LSE 259. The hydrogen deficiency is accompanied by nitrogen- and carbon-enrichment in LSE 153 and LSE 259 only. In LSE 263, carbon is depleted by about 1 dex. Stellar masses obtained by assuming that a core mass-luminosity relation holds for these stars, were found to be in the range 0.6-0.9 solar mass, yielding luminosities log L/L:solar = 3.7-4.5. Two of the program stars (LSE 153 and 259) appear to be possible successors of the R CrB and helium B stars, whereas the third star (LSE 263) displays a much lower carbon content in its photosphere making it an exceptional case among the known hydrogen deficient stars.

Husfeld, D.↗
DOE OSTI2024

Globally Gridded Groundwater Extraction Volumes and Costs under Six Depletion and Ponded Depth Targets

This repository contains simulated outputs from superwell – a hydro-economic tool for long-term assessment of groundwater cost and supply – providing globally gridded groundwater extractable volumes and associated unit costs ($/km³) for accessible groundwater production, based on a variety of user-defined depletion and ponded depth scenarios. Key model documentation: Niazi, H., Ferencz, S. B., Graham, N. T., Yoon, J., Wild, T. B., Hejazi, M., Watson, D. J., & Vernon, C. R. (2025). Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1. Geoscientific Model Development, 18(5), 1737-1767. https://doi.org/10.5194/gmd-18-1737-2025 Find the source code of the superwell model on GitHub: https://github.com/JGCRI/superwell Repository Overview Main output: superwell_outputs.7z contains 6 files (4.5 GB) named as superwell_py_deep_all_0.*PD_0.*DL.csv. These files present superwell outputs of global groundwater extraction volumes and cost estimates on a 0.5° scale for six scenarios with different Ponded Depth (PD; 0.3 and 0.6 m) and Depletion Limit (DL; 5%, 25%, and 40% of available volume) targets over the entire pumping lifetime of a grid cell superwell_py_deep_all_0.3PD_0.25DL_sample_100.csv contains superwell outputs for 100 data points sampled to match the global inputs' distribution superwell_py_deep_all_0.3PD_0.25DL_Grid_72548.csv contains superwell output for a single grid cell concept_v5.png provides an overview of the superwell workflow Outputs Description year_number: year of pumping depletion_limit: set depletion limit (DL) as a volume fraction of total available groundwater Mappings: continent, country, gcam_basin_id, Basin_long_name, grid_id: geographic identifiers and basin information Inputs: grid_area (km²): area of the grid cell whyclass: hydrogeological classification of the aquifer permeability (m/day), porosity (%), total_thickness (m), depth_to_water (m): aquifer properties. The geo-processed input data has been published separately: https://doi.org/10.57931/2307831 Model outputs: orig_aqfr_sat_thickness (m), aqfr_sat_thickness (m): original and remaining/instantaneous saturated thickness of the aquifer hydraulic_conductivity (m/day), transmissivity (m²/day): hydraulic properties of the aquifer radius_of_influence (m), areal_extent (km²): well radius and area of influence from the center of the well number_of_wells (-): number of wells in a grid cell determined by a ratio of well area and grid area max_drawdown (m), drawdown (m), drawdown_interference (m): well and aquifer drawdown during extraction total_head (m): total lift for the groundwater (depth to water plus drawdown) total_well_length (m): total depth of wells drilled well_yield (m³/day): pumping rate or well yield power (kW), energy (kWh): power and energy required for pumping groundwater Volume Outputs: volume_produced_perwell (m³), cumulative_vol_produced_perwell (m³): production volume metrics per well volume_produced_allwells (m³), cumulative_vol_produced_allwells (m³): aggregate extraction volumes for all wells in a grid cell available_volume (m³): available groundwater in storage for the grid cell as determined by aquifer properties depleted_vol_fraction: fraction of total volume pumped over available volumes in a grid cell (same as depletion limit) Cost Outputs: well_installation_cost ($): well installation cost based on the hydrogeological complexity of the aquifer annual_capital_cost, maintenance_cost, nonenergy_cost ($): nonenergy costs energy_cost_rate ($/kWh): electricity rate energy_cost ($): energy cost of pumping groundwater total_cost_perwell ($), total_cost_allwells ($): total annual energy and non-energy cost for each and all wells in a grid cell a unit_cost ($/m³), unit_cost_per_km3 ($/km³), unit_cost_per_acreft ($/acre-ft): total cost of pumping a unit of groundwater, indicated for different spatial units Key Resources Model documentation: Niazi, H., Ferencz, S., Graham, N., Yoon, J., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024; In-prep). Long-term Hydro-economic Assessment Tool for Evaluating Global Groundwater Cost and Supply: Superwell v1. Geoscientific Model Development. Input data: Niazi, H., Watson, D., Hejazi, M., Yonkofski, C., Ferencz, S., Vernon, C., Graham, N., Wild, T., & Yoon, J. (2024). Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins. MSD-LIVE Data repository. https://doi.org/10.57931/2307831 superwell source code: https://github.com/JGCRI/superwell Cite as Niazi, H., Ferencz, S., Yoon, J., Graham, N., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024). Globally Gridded Groundwater Extraction Volumes and Costs under Six Depletion and Ponded Depth Targets. MSD-LIVE Data repository. https://doi.org/10.57931/2307832 Contact Reach out to Hassan Niazi or Stephen Ferencz or open an issue in the superwell repository for questions or suggestions.

Earth Systems↗
DOE OSTI2024

UT-GOM2-2 Preliminary Report: Terrebonne Basin Northern Gulf of Mexico, 30 July-28 September 2023

In the summer and fall of 2023, the Gulf of Mexico Deepwater Hydrate Coring Expedition (UT-GOM2-2) drilled, cored, made downhole measurements, and analyzed samples from the seafloor to the base of the gas hydrate stability zone in one location (Site H, WR313) in the Terrebonne basin, deepwater Gulf of Mexico. Analyses of data and samples from the expedition will inform biological, geochemical, and geomechanical models to constrain the role of gas hydrates in the carbon cycle and the potential for gas hydrates as an energy resource. Pressure and conventional cores were collected continuously to a depth of 155.1 meters below the seafloor (mbsf). At deeper depths, cores were taken periodically from hydrate-bearing sands and their bounding muds to a total depth of 861.3 mbsf. 162.6 m of conventional core and 54.8 m of pressure core were obtained. Twelve temperature measurements were made between 27.1 and 144.5 mbsf to determine the geothermal gradient. At the seafloor, more than 4 m of sandy silt of unknown origin was encountered. Beneath this sand, to a depth of ~200 mbsf, the section was composed of interbedded mud and biogenic carbonate ooze. The biogenic ooze correlated to low density and high porosity intervals observed in the previously acquired logging while drilling (LWD) data and as measured. Calcareous nannofossil biostratigraphy constrains the entire record to the Pleistocene (< 0.91 million years) with a pronounced increase in sedimentation rate with depth. Beneath 200 mbsf, the section was predominantly composed of mud with two thicker, hydrate-bearing coarse-grained intervals, which are commonly known as the Blue and Orange sands. The dissolved gas concentration was quantified from pressure cores. In the shallow section, dissolved methane concentration increased below the sulfate-methane transition zone (SMTZ) and reaches saturation (the limit of solubility for methane) at 147 mbsf. Gas expansion was very common in conventional and depressurized pressure (conventionalized) cores below the SMTZ. At deeper depths, the methane concentration within muds bounding the Blue and Orange reservoirs was generally found to be less than saturation. The dissolved and hydrate gas composition is consistent with a microbial source, containing greater than 99.99% methane and only trace concentrations of ethane, propane, and butane. The methane to ethane ratio (C 1 /C 2 ) and the methane to ethane plus propane (C 1 /(C 2 +C 3 )) decrease with depth down to at least 678 mbsf, mainly driven by the increase in ethane with depth. It is unclear if this trend continues through the Orange sand interval. The δ 13 C isotopic signature of methane ranges between -69.9 and -78.5 ‰ Vienna Pee Dee Belemnite (VPDB). Pressure core recovery of all sandy intervals was poor. However, pressure core logs of the Orange sand show intervals of low density and high velocity, which are indicative of high hydrate saturation. One core from within the Orange sand was composed of interbedded graded sandy silt and mud. The sandy silts from this core are composed of mainly quartz and feldspar with some lithics. Most of the recovered pressure core samples are maintained at near in-situ pressure and temperature (within the hydrate stability field) at the University of Texas Pressure Core Center awaiting analysis. In the shallow section, samples will be used to determine the flux of organic carbon through the basin system, find the rate at which that carbon was consumed, and understand the microbial population responsible for these processes. In the deeper section, samples from in and around the hydrate reservoirs will be used to determine the petrophysical properties of the reservoir and bounding seals in these systems.

03 NATURAL GAS↗
DOE OSTI2026

Qualitative Risk Assessment of Legacy Wells within the Estimated Prairie State Generating Company Area of Review

This report details the digitization of a legacy wellbore database, including data processing assumptions, parameter estimation, and risk assessment methodology. The database, comprising 6,454 documents, was provided by ISGS. It includes valuable data from the Prairie State Generating Company (PSGC) and One Earth Energy (OEE) sites of the CarbonSAFE Phase III – Illinois Storage Corridor project. The report focuses on wells within a 15-mile radius from the Lively Grove #1 (LG#1) well at PSGC site, evaluating subsurface conditions and potential risks. A total of 4,386 wellbores within 15 miles of the LG#1 well were filtered based on depth and formation codes. LG#1 is the stratigraphic well at the PSGC site drilled in 2021. Ninety-four (94) wells penetrating the Maquoketa Shale Group (the primary confining unit) within the estimated area-of-review (AoR) for the PSGC site were evaluated using a qualitative risk assessment (QRA) methodology. The QRA developed by Arbad et al. 2022 focuses on legacy wells within the AoR and categorizes them based on well construction details. The QRA identifies wells that need immediate attention by categorizing them based on penetration depth and protection. Wells within the AoR were categorized into nine groups based on penetrations and protections. These categories range from Type 1 wells, with no documentation, to Type 9 wells, which do not penetrate the primary confining unit or storage reservoir (unit). Well accessibility within the AoR varies based on well status, including Dry & Abandoned (DA), Plugged & Abandoned (PA), Injection (INJ), Oil/Gas Producing (PROD), and Observation (Obs) wells. Accessibility levels were determined by well construction, with DA wells being the least accessible and Observation wells the most accessible, impacting gas leakage detection possibilities. Remedial action priority of wells decreases from Type 1 to Type 9 wells. Type 1 to Type 6 wells with status DA and PA require immediate attention, while Type 7 and Type 8 wells are low priority. A risk matrix used to prioritize corrective actions for legacy wells is proposed to categorize wells within an AoR based on penetrations, protections, and accessibility. The methodology involves data acquisition, well categorization into nine types, and determining CO 2 leakage pathways using well schematics and geospatial mapping. This approach is particularly useful for managing the integrity of legacy wells throughout the lifecycle of a Carbon Capture and Storage (CCS) project. A qualitative risk assessment of 94 wells within the AoR of the PSGC site identified 54 wells with high priority for corrective action due to penetration of the primary containment seal. The assessment utilizes color-coded maps to categorize well types and prioritize corrective actions, providing a comprehensive analysis. Schematics of wells penetrating the primary confining unit were drawn, and leakage pathways were identified. Details of all wells penetrating the confining zone are provided in the appendix, including information on well types, plugging, and casing status.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

GOOML - Finding Optimization Opportunities for Geothermal Operations: Preprint

Geothermal Operational Optimization with Machine Learning (GOOML) is a transferable and extensible component-based geothermal asset modeling framework that considers complex steamfield relationships and identifies optimization prospects using a data-driven approach. We have used this framework to develop digital twins that provide steamfield operators with an operational environment to analyze and understand historical and forecasted power production, explore new steamfield configuration possibilities, and seek optimal asset management for real world applications. The GOOML modeling software is built on a generic component-based systems framework that allows for both historical and forecast analysis. A GOOML model can perform historical data-assimilation using first-principal thermodynamics to create a meaningful data model. Historical production data can then be coupled with a forecast framework to train machine-learning models of steamfield components to predict future outputs. This modeling environment enables digital exploration of steamfield design configurations and operational scenarios. GOOML digital twins have been developed for steamfields in New Zealand and the United States representing differing power generation and field conditions. These digital twins have been validated by comparing hindcast predictions against historical production data. Reinforcement learning experiments were conducted to demonstrate the ability to programmatically explore the operations space using machine learning agents. Our initial results are compelling; two to five percent increases in annual energy production were demonstrated by the GOOML models with no additional infrastructure build required. GOOML offers a new approach to geothermal operations by applying state-of-the-art machine learning algorithms, comprehensive data analytics, and interaction with digital twins. Through application of these tools, operators will realize greater availability and higher net generation which will increase the cost effectiveness of geothermal energy projects.

access↗
DOE OSTI2025

Expedition UT-GOM2-2 Site H

Pressure and conventional cores were collected at Site H of the Walker Ridge Protracted Area Block 313 in the Terrebonne Basin, deepwater Gulf of America (Gulf of Mexico) during the University of Texas (UT) Deepwater Hydrate Coring Expedition (UT-GOM2-2). Pressure and conventional cores were collected continuously to a depth of 155.1 meters below the seafloor (mbsf). At deeper depths, cores were taken periodically from hydrate-bearing sands and their bounding muds to a total depth of 861.3 mbsf. 162.6 m of conventional core and 54.8 m of pressure core were recovered. Twelve temperature measurements were made between 27.1 and 144.5 mbsf to determine the geothermal gradient. At the seafloor, more than 4 m of sandy silt of unknown origin was encountered. Beneath this sand, to a depth of about 200 mbsf, the section was composed of interbedded mud and biogenic carbonate ooze. The ooze correlated to low density and high porosity intervals observed in the previously acquired logging while drilling (LWD) data and as measured. These ooze intervals also correspond to lighter sediment color, increased Ca content based on X-ray florescence (XRF) core scanning, and increased calcareous nannofossil abundance. Calcareous nannofossil biostratigraphy constrains the entire record to the Pleistocene (< 0.91 million years), with a pronounced increase in sedimentation rate with depth. Below 200 mbsf, the section was predominantly composed of mud with two thicker, hydrate-bearing coarse-grained intervals, which are commonly known as the Blue and Orange sands. The dissolved gas concentration was quantified from pressure cores. In the shallow section, dissolved methane concentration increased below the sulfate-methane transition zone (SMTZ) and reaches saturation (the limit of solubility for methane) at 147 mbsf. Gas expansion was very common in conventional and depressurized pressure (conventionalized) cores below the SMTZ. At deeper depths, the methane concentration within muds bounding the Blue and Orange reservoirs was generally found to be less than saturation. The dissolved and hydrate gas composition is consistent with a microbial source, containing greater than 99.99% methane and only trace concentrations of ethane, propane, and butane. The methane to ethane ratio (C1 /C2 ) and the methane to ethane plus propane (C1 /(C2 +C3 )) decrease with depth down to at least 678 mbsf, mainly driven by the increase in ethane with depth. It is unclear if this trend continues through the Orange sand interval. The δ13C isotopic signature of methane ranges between -69.9 and -78.5 ‰ relative to the Vienna Pee Dee Belemnite (VPDB) standard. Pressure core recovery in sandy intervals was poor. However, pressure core logs of the Orange sand show intervals of low density and high velocity, which are indicative of high hydrate saturation. One pressure core was degassed and the average hydrate saturation in the core was determined to be 24%. One core from within the Orange sand was composed of interbedded graded sandy silt and mud. The sandy silts from this core are composed of mainly quartz and feldspar with some lithics. Most of the recovered pressure core samples are maintained at near in-situ pressure and temperature (within the hydrate stability field) at the University of Texas Pressure Core Center awaiting analysis. In the shallow section, samples will be used to determine the flux of organic carbon through the basin system, find the rate at which that carbon was consumed, and understand the microbial population responsible for these processes. In the deeper section, samples from in and around the hydrate reservoirs will be used to determine the petrophysical properties of the reservoir and bounding seals in these systems.

03 NATURAL GAS↗
NASA NTRS2000

Geology and geochemistry of paleosols developed on the Hekpoort Basalt, Pretoria Group, South Africa

The Hekpoort paleosols comprise a regional paleoweathering horizon developed on 2.224 +/- 0.021 Ga basaltic andesite lavas at the top of the Hekpoort Formation of the Pretoria Group, Transvaal Supergroup, South Africa. In five separate profiles, from outcrops along road cuts near Waterval Onder and the Daspoort Tunnel and in three drill cores from the Bank Break Area (BB3, BB8, and BB14), the top of the paleosol is a sericite-rich zone. The sericite zone grades downward into a chlorite-rich zone. In core BB8 and in the road cut at the Daspoort Tunnel, we sampled the underlying or parent basaltic andesite into which the chlorite zone grades. We did not obtain samples of the parent material at Waterval Onder and in cores BB3 and BB14, but chemical analyses indicate that the chlorite and sericite zones in these profiles derive from underlying lavas similar to the ones we sampled in core BB8 and at the Daspoort Tunnel. The presence of apparent rip-up clasts of the paleosol in the overlying ironstones of the Strubenkop Formation in the cores from Bank Break makes it very unlikely that most of the alteration was a result of interactions with hydrothermal fluids. Desiccation cracks at the top of the paleosol that were filled with sand during the deposition of the overlying sediments at Waterval Onder point to a subaerial weathering origin. Very little, if any, Al, Ti, Zr, V, or Cr moved a discernible distance during weathering of any of the five profiles. The vertical distribution of Fe, Mg, Mn, Ni, and Co indicates that these elements were largely removed from the top of the soil during weathering. The overall abundance of these elements in each of the profiles indicates that a significant fraction of the complement lost from the top subsequently reprecipitated in the lower portion of the soil as constituents of an Fe2(+) -rich smectite. The loss of Fe from the top of the soil during weathering of the Hekpoort paleosols indicates that atmospheric PO2 was less than 8 x 10(-4) atm about 2.22 Ga. Fe2(+) -rich smectite should only precipitate during soil formation if atmospheric PCO2 is less than or equal to 2 x 10(-2) atm (Rye, Kuo, and Holland, 1995). Ca and Na were largely lost during weathering. Some Na was apparently added to the sericite zone in cores BB3, BB8, and BB14 after weathering. All five profiles are enriched in K and Rb, and most are enriched in Ba. The distribution of these elements indicates that they all were added during post-weathering hydrothermal metasomatism. Rb-Sr analysis of the paleosol at the Daspoort Tunnel indicates that metasomatism last affected that profile 1.925 +/- 0.032 Ga (Macfarlane and Holland, 1991).

Non-NASA Center↗
DOE OSTI2022

Using Natural Gas Liquids to Recover Unconventional Oil and Gas Resources (Final Report)

This document presents final technical findings for the project Using Natural Gas Liquids to Recover Unconventional Oil and Gas Resources (FE0031782). The project is part of the U.S. Department of Energy Oil and Gas Program to develop and advance technologies that can significantly improve the recovery efficiencies of unconventional oil and gas resources. The overall objective of this project is to improve the ultimate recovery from unconventional oil and gas (UOG) resources in the United States by developing a method for using unrefined natural gas liquids (NGLs) as treatment fluids to improve hydrocarbon production. Horizontal extended-lateral drilling coupled with high volume hydraulic fracturing has significantly increased production from UOG resources in the U.S. However, the recovery efficiency is low compared to the estimated oil and gas in place. Recent data indicate that less than 10% of the oil in the liquid-rich UOG reservoirs is produced. Alternative completion methods using NGLs could increase production (Battelle, 2016; Wan et al., 2013; Wan, 2013; Downey et al, 2021); however, field validation tests are needed to develop an approach that is economical, efficient, and compatible in the UOG setting to advance towards commercial deployment. This project aims to develop and field test a method to improve recovery of oil resources in UOG shale plays by using Y-Grade NGLs, or a similar combination of NGLs, as treatment fluids. Refined NGLs have been used as a hydraulic treatment fluid in UOG plays for decades and are shown to be particularly effective because their miscibility with oil allows oil to flow more freely; however, the use of Y-Grade (unrefined) NGLs has not been studied. The use of Y-Grade NGLs would be advantageous over refined NGLs because Y-Grade NGLs do not require infrastructure or investment in refining and are already being produced from many UOG reservoirs. The concept was tested and monitored in the field at a commercial well site owned by project partner Hopco, Ltd. The project team, which consists of multiple oil and gas operators, Linde Gas North America LLC (Linde) and the Ohio Division of Geological Survey (ODGS), has extensive experience with oil and gas production in the Appalachian basin and the ability to work together quickly to solve technical issues and research needs. A key part of the proposed work was the use of existing wells for field testing and monitoring. A total of four wells (three vertical and one horizontal) were available for this project. One of the vertical wells was utilized as the test well for the NGL treatment test. A nearby vertical well was used for microseismic monitoring. The remaining vertical well and the horizontal well provided a baseline for typical UOG production in the oil window of the Utica/Point Pleasant (UPP). Major technical tasks of the project include characterization of the geotechnical properties of the UPP with an emphasis on the field site; design and planning for the NGLs testing; field testing and monitoring; analysis and integration of field data; and economic and resource/reserve assessment. The Shoman monitoring well and Doughty NGL treatment well were successfully plugged back during September-October 2020 in preparation for treatment and monitoring. A nitrogen diagnostic fracture injection test (DFIT) was completed on the Doughty well on July 22, 2021, consisting of 133,000 scf (91 Bbl.) of nitrogen. A nitrogen foam frac was completed in the Utica-Point Pleasant interval on August 17, 2021 with funding from outside sources. A microseismic monitoring array was installed in the Shoman well and monitored microseismic activity during the Doughty well frac job. Y-Grade NGL injection commenced on August 26, 2021. A total of 215 Bbl. was injected but the job was shut down due to a small leak on the suction hose on the pump truck. The Y-Grade treatment resumed on August 27, 2021 and an additional 726 Bbl. of Y-Grade NGL was injected at a well head pressure of 3850 psi. Total volume of injected Y-Grade over the two days of injection was 941 Bbl. The well was shut-in for 17 days following injection to allow the Y-Grade NGLs to soak on the formation. Y-Grade treatment flow back commenced on 9/13/2021 on a weekly basis. Production data, including surface pressures, oil, nitrogen, natural gas, and flow times was measured and recorded. Periodic gas samples were collected and analyzed to determine composition of flowback gas. As of July 2022, the treatment well had produced 726 Bbl of oil and 2,888 mcf gas. In August 2022, tubing and packer in the well were removed and a pump was installed to enhance oil recovery. Currently, the operator is producing the well about 2 days a week for a few hours. The performance of the NGL treatment test was evaluated based on reservoir simulations of the treatment process, processing of well testing data, analysis of micro seismic monitoring data, and production data analysis. This analysis suggested that oil production in the small test would continue through 2025. Upscaling the treatment to a horizontal Utica Point-Pleasant well would allow more oil recovery, but the process would involve more investment, services, and operational support. An economic analysis was conducted for scenarios aimed at upscaling the NGL treatment process for more typical horizontal UPP wells in the Appalachian Basin.

02 PETROLEUM↗
DOE OSTI2023

Hawai‘i Play Fairway (Final Report)

Most of Hawai'i's geothermal resources are blind—their manifestations, such as hot springs and steam vents, do not appear on the ground surface because the heated water flows far below. With the exception of $K\bar{i}lauea East Rift Zone$, in most areas of Hawai'i, high lateral permeability in the first kilometer below ground surface prevents surface thermal features from developing. As a methodology for discovering these blind resources, Play Fairway Analysis (PFA) involves finding potential locations of blind hydrothermal systems and describing potential geothermal sources in rift-zone settings. Using the PFA to find Hawai'i's geothermal resources, the University of Hawai'i (UH) conducted the Hawai‘i Play Fairway Project, Hawai'i's first statewide geothermal resource assessment since 1985. Sponsored by the U.S. Department of Energy, the Hawai'i Play Fairway Project provided an updated resource assessment, a roadmap for additional exploration activities, and the identification of areas for further exploration. Benefitting from UH's core competency in earth sciences and experienced geothermal researchers, the project comprised three phases. During the first phase, the team identified, compiled, and ranked existing geologic, groundwater, and geophysical datasets relevant to subsurface heat, fluid and permeability. Using a Bayesian statistical approach, the team developed a statistical methodology to integrate these data into a resource probability map. The team evaluated the confidence in the probability value and considered development viability of areas with geothermal resources. With these analyses, the team identified 10 locations in the Hawaiian Islands for exploration activities. For the second phase, the team collected new groundwater data in 10 locations across the state and new geophysical data on $L\bar{a}na‘i, Maui$, and central Hawai'i Island and modeled topographically induced stress to better characterize subsurface permeability. Analyzing the subsurface stresses, the team evaluated the potential for fracture-induced permeability. The team inverted the MT and gravity data to produce 3D models of resistivity and density, respectively, on $L\bar{a}na‘i$, across $Haleakal\bar{a}'s$ SW rift (Maui), and surrounding Mauna Kea (Hawai‘i Island). The team developed and applied a new method for incorporating depth information about resistivity, density, and potential for fracture-induced permeability into the statistical method for computing resource probability in these three focus areas. The team incorporated the new groundwater results with the new geophysical results and the calculations of potential for fracture-induced permeability to produce updated maps of resource probability and confidence. Through combining data from the first and second phases, the team determined locations for further exploration during the third phase. For MT and gravity surveys, the team recommended $Kaua'i's$ $L\bar{i}hu'e$ $Basin$, the east rift of $Maui's$ $Haleakal\bar{a}$ volcano, and the southwest rift of Hawai'i Island's Mauna Loa volcano. The MT and gravity surveys aimed to enable improved confidence in the resource potential in these locations. For drilling deep groundwater well(s), the team recommended Southeast Mauna Kea and $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$. During the third phase, further exploration involved drilling a groundwater well in $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$ and performing more geophysical surveys. We deepened an existing water well proximal to our target area on $L\bar{a}na'i$ due to funding constraints that precluded us from spudding a new well that would exceed 1km depth. Drilling was preceded by a number of substantial elements including: writing an Environmental Assessment and the subsequent legal process, performance of deviation logging, lowering a camera down the well, coordinating site preparation with $P\bar{u}lama$ $L\bar{a}na'i$, shipping the UH-owned rig interisland, procuring supplies, and leading 3 community meetings on $L\bar{a}na'i$. Drilling occurred 24/7 the entire month of June 2019 over which time $L\bar{a}na'i$ $Well$ 10 was deepened from 427 m to 1057 m, with continuous core collected. We measured a roughly linear temperature gradient averaging 42°C/km and a maximum bottom hole temperature of 66°C. This gradient is more than twice the background for Hawai'i and within a range of gradients measured in this depth range for some exploration wells within KERZ. We consider these results encouraging for $L\bar{a}na'i's$ resource potential and recommend following with a slim hole within $L\bar{a}na'i's$ caldera (our target zone) to ~ 2 km. Further, the positive implications such results have for the island of O‘ahu are substantial - the shield stage of O'ahu's volcanoes ended 1-2 My earlier. However, O'ahu uses more electricity than the rest of the islands combined, and the utility recently called for 500-700MW of firm, dispatchable renewable electricity on O'ahu by 2033. In Phase 3, we also collected limited new encouraging groundwater data, and updated our thoughts on the probabilities of fluid and permeability at resource depths (Pr F = 1; Pr P = mostly unconstrained). Ultimately, we advocate for using our final probability of heat, and confidence in this probability, to drive the next phase of exploration. We contend further development of geothermal in Hawai‘i will enable the state to achieve its 100% renewable policy objective and Hawai'i to transition off of fossil fuels through geothermal discovery and development. The project not only produced a large amount of data and expanded the existing knowledge of Hawai'i's geothermal resources, but also produced publications, theses, presentations, core photos, datasets, media reports, television interviews, community events, and a blog. Students and new professionals benefitted from the project's hands-on research experiences and educational opportunities and earned awards and recognition.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Newly acquired core enhances geologic understanding of the northern Paradox Basin Cane Creek play, southeastern Utah

The Cane Creek interval is an emerging unconventional heterolithic play within the Pennsylvanian Paradox Formation of the northern Paradox Basin, southeastern Utah. Structural controls (basement uplifts, salt movement, faults, fractures) likely control production from the Cane Creek, where total oil production of ~10 MMBO is just a fraction of current estimates of the undiscovered resource (up to 1.2 BBOE). Oil production has been most successful from unstimulated horizontal wells in the Cane Creek Unit in the central play area, whereas areas to the north (Greentown and Gunnison Valley) have to date proven to be largely unproductive from vertical and horizontal wells drilled without the benefit of 3D seismic or the deployment of modern hydraulic fracturing techniques. However, the limited historical drilling has shown significant promise (e.g., significant shows). The lack of core data in the northern part of the play has impeded thorough assessment of reservoir quality, character, and lithological facies heterogeneity. The Cane Creek is informally divided into three distinct zones. The upper A and basal C zones are composed of fabric destructive and wavy bedded anhydrite, dolomitic mudstone, and organic-rich algal laminated mudstones (source rock). The middle B zone is generally wave rippled and burrowed, low-permeability sandstone-siltstone (reservoir) with subordinate displacive anhydrite. Recent potash cores collected from the Salt Valley anticline near Crescent Junction, Utah (but east of the main play area), show thicker reservoir sandstone/siltstone packages (~40 ft) in the B zone, and parts of the C zone, in the northern part of the Paradox Basin in comparison to the central (~35 ft) and southern (~30 ft) areas. Source rock analyses from northern core/cuttings also indicate that it is well positioned within the oil window (VRo ~0.8), with thin organic-rich mudstone beds containing up to 20 wt% TOC. Because the siliciclastic reservoir rocks have low permeabilities (0.009–0.202 mD) and variable porosities (6–17%) due to dolomite-anhydrite cement, naturally occurring and possibly stimulated fractures may be essential for hydrocarbon recovery. To better understand the northern part of the Cane Creek play, a new core was obtained from the research stratigraphic well State 16-2 (January 2021). Details from this new core will help characterize reservoir quality via core analyses, petrography, fracture orientations, geomechanics, and correlation of facies and depositional trends. Additionally, an evaluation of matrix porosity/permeability will be conducted to determine its contribution to productivity.

02 PETROLEUM↗
DOE OSTI2021

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs Results

Geothermal power plants typically show decreasing heat and power production rates over time. Mitigation strategies include optimizing the management of existing wells - increasing or decreasing the fluid flow rates across the wells - and drilling new wells at appropriate locations. The latter is expensive, time-consuming, and subject to many engineering constraints, but the former is a viable mechanism for periodic adjustment of the available fluid allocations. Data and supporting literature from a study describing a new approach combining reservoir modeling and machine learning to produce models that enable strategies for the mitigation of decreased heat and power production rates over time for geothermal power plants. The computational approach used enables translation of sets of potential flow rates for the active wells into reservoir-wide estimates of produced energy and discovery of optimal flow allocations among the studied sets. In our computational experiments, we utilize collections of simulations for a specific reservoir (which capture subsurface characterization and realize history matching) along with machine learning models that predict temperature and pressure timeseries for production wells. We evaluate this approach using an "open-source" reservoir we have constructed that captures many of the characteristics of Brady Hot Springs, a commercially operational geothermal field in Nevada, USA. Selected results from a reservoir model of Brady Hot Springs itself are presented to show successful application to an existing system. In both cases, energy predictions prove to be highly accurate: all observed prediction errors do not exceed 3.68% for temperatures and 4.75% for pressures. In a cumulative energy estimation, we observe prediction errors that are less than 4.04%. A typical reservoir simulation for Brady Hot Springs completes in approximately 4 hours, whereas our machine learning models yield accurate 20-year predictions for temperatures, pressures, and produced energy in 0.9 seconds. This paper aims to demonstrate how the models and techniques from our study can be applied to achieve rapid exploration of controlled parameters and optimization of other geothermal reservoirs. Includes a synthetic, yet realistic, model of a geothermal reservoir, referred to as open-source reservoir (OSR). OSR is a 10-well (4 injection wells and 6 production wells) system that resembles Brady Hot Springs (a commercially operational geothermal field in Nevada, USA) at a high level but has a number of sufficiently modified characteristics (which renders any possible similarity between specific characteristics like temperatures and pressures as purely random). We study OSR through CMG simulations with a wide range of flow allocation scenarios. Includes a dataset with 101 simulated scenarios that cover the period of time between 2020 and 2040 and a link to the published paper about this project, where we focus on the Machine Learning work for predicting OSR's energy production based on the simulation data, as well as a link to the GitHub repository where we have published the code we have developed (please refer to the repository's readme file to see instructions on how to run the code). Additional links are included to associated work led by the USGS to identify geologic factors associated with well productivity in geothermal fields. Below are the high-level steps for applying the same modeling + ML process to other geothermal reservoirs: 1. Develop a geologic model of the geothermal field. The location of faults, upflow zones, aquifers, etc. need to be accounted for as accurately as possible 2. The geologic model needs to be converted to a reservoir model that can be used in a reservoir simulator, such as, for instance, CMG STARS, TETRAD, or FALCON 3. Using native state modeling, the initial temperature and pressure distributions are evaluated, and they become the initial conditions for dynamic reservoir simulations 4....

15 GEOTHERMAL ENERGY↗