Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “subsurface 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.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Reservoir Characterization and Techno-Economic Analysis of Enhanced Geothermal and Closed-Loop Systems in the Wattenberg Field of the Denver-Julesburg Basin, Colorado

This report provides details of reservoir characterization and techno-economic analysis for enhanced geothermal systems and closed-loop geothermal systems in the Denver-Julesburg Basin Wattenberg area. The analysis contributes to the Geothermal Limitless Approach to Drilling Efficiencies (GLADE) project, a multi-institutional research initiative focused on advancing geothermal energy development by reducing drilling costs and improving the rate of penetration. The GLADE project brings together national laboratories, academic institutions, and industry partners, including the U.S. Department of Energy's Geothermal Technologies Office, Occidental Petroleum, the National Laboratory of the Rockies (NLR), Los Alamos National Laboratory, Colorado School of Mines, Louisiana State University, Texas A&M University, and several drilling technology companies. As part of this effort, NLR developed and applied slender-body theory modeling tools and the GEOPHIRES simulator to support techno-economic analysis. The slender-body theory wellbore simulator evaluated thermal performance in U-loop, Eavor-type multilateral, and enhanced geothermal system well configurations. Results indicated that an optimized single U-loop design delivers an outlet temperature of 168.8 degrees Celsius,stabilizing at a final outlet temperature of 135.3 degrees Celsius. Heat production also improved significantly, with an average thermal output of 16.9 MWth and a final thermal output of 11.0 MWth. These results provided critical baseline data for selecting optimized systems for techno-economic analysis. Building on these outputs, NLR applied the GEOPHIRES techno-economic simulator to estimate the levelized cost of electricity for both the Eavor multilateral and enhanced geothermal systems. The sensitivity analysis showed that drilling costs exert a particularly strong impact on the levelized cost of electricity of multilateral systems. Therefore, advances in drilling technology and improved drilling rates could yield significant gains in the economic performance of these systems. In this context, the GLADE project directly supports the U.S. Department of Energy's mission to reduce the cost of geothermal energy and accelerate the deployment of next-generation geothermal technologies.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

3-D Geological Modeling for Numerical Flow Simulation Studies of Gas Hydrate Reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. Here, this study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-stepping grid were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a three-dimensional (3-D) trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

58 GEOSCIENCES↗
DOE OSTI2024

Soil Gas Survey Results Supporting Groundwater Correction Action Plan (GCAP) Development for the Moab Site

A soil gas survey was performed at the Moab Uranium Mill Tailings Remedial Action (UMTRA) Project Site during the week of November 6, 2023. Soil gas surveys are used to characterize residual subsurface sources of volatile contaminants, such as volatile organic compounds, as well as contaminants that generate a surrogate indicator gas or otherwise influence soil gas composition. The primary objective of the Moab soil gas survey was to confirm, identify, quantify, and refine secondary contaminant source area locations for uranium and ammonium/ammonia (NH 4 + /NH 3 ) in the vadose zone and shallow groundwater. The overarching goal was to provide information to assist in developing the technical basis for the Groundwater Compliance Action Plan (GCAP). Specifically, the soil gas data will support the deployment of source control technologies; e.g., where supplementary capping, removal actions, or amendments might be beneficial.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Soil Gas Survey Results Supporting Groundwater Correction Action Plan (GCAP) Development for the Moab Site

A soil gas survey was performed at the Moab Uranium Mill Tailings Remedial Action (UMTRA) Project Site during the week of November 6, 2023. Soil gas surveys are used to characterize residual subsurface sources of volatile contaminants, such as volatile organic compounds, as well as contaminants that generate a surrogate indicator gas or otherwise influence soil gas composition. The primary objective of the Moab soil gas survey was to confirm, identify, quantify, and refine secondary contaminant source area locations for uranium and ammonium/ammonia (NH 4 + /NH 3 ) in the vadose zone and shallow groundwater. The overarching goal was to provide information to assist in developing the technical basis for the Groundwater Compliance Action Plan (GCAP). Specifically, the soil gas data will support the deployment of source control technologies; e.g., where supplementary capping, removal actions, or amendments might be beneficial.

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

Hydro-EKG: In-Stream Streaming-potential (SP) Electrical Monitoring of the Lifeblood of Watersheds during Disturbance

In order to develop improved methods for characterization and monitoring of hyporheic exchange along stream and river corridors, this project evaluated and compared concentration (C) and discharge (Q) correlations between the dam-managed, intermittent Rio Grande River and the perennial East Fork Poplar Creek in Oak Ridge, Tennessee. Additionally, a multimethod geoelectric survey was implemented between January and March 2022 along a 220-meter long reach of the bedrock-lined East Fork Poplar Creek in Oak Ridge, Tennessee to identify locations of surface-water and groundwater exchange in the bedrock-lined streambed and characterize the subsurface flow paths that convey water between the stream and flood plain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

FracML: A Machine Learning Based Tool to Quantify Reservoir Scale Fracture Network for CO2 Storage

Poster on “FRACML: A Machine Learning Based Tool to Quantify Reservoir Scale Fracture Network for CO2 Storage” for the CCUS 2025 conference held in Houston, Texas March 3-5, 2025. The accurate characterization of subsurface fracture networks is essential for the secure operation of carbon capture, utilization, and storage (CCUS) projects. A thorough understanding of the spatial distribution of subsurface faults and fractures is crucial for predicting CO2 plume evolution and minimizing risks such as potential leakage into overlying formations or induced seismicity. In this context, robust fracture network quantification plays a pivotal role in reservoir management, providing the data necessary to fine-tune operational parameters, and ensure the environmental and economic viability of CCUS projects. As part of the U.S. Department of Energy’s SMART (Science-informed Machine Learning for Accelerating Real-time Decisions in Subsurface Applications) initiative, we focused on the development and application of a machine learning-based tool (FRACML) designed to quantify and map fracture networks using real-world (non-synthetic) data from an active CO2 injection site. Our objective is to demonstrate the utility of this tool in improving operational efficiency and safety across CCUS sites.

artifical intelligence / machine learning (AI/ML)↗
DOE OSTI2026

One Earth Energy Seismic Interpretation

The objectives of the Illinois Storage Corridor (ISC) project are to accelerate commercial deployment of carbon capture utilization and storage at two individual sites and receive approvals for Underground Injection Control (UIC) Class VI permits for construction at each site (ISC Project Narrative, 2020). As part of this project, and as part of the subsurface geologic characterization, 2D seismic data was acquired at both sites. This report summarizes the findings from the 2D and 3D seismic interpretation at the One Earth Energy site near Gibson City, Illinois. The seismic data confirms the stratigraphic continuity of the Mt. Simon Arkose Zone storage interval and the Eau Claire confining unit across the project area. The seismic data also indicates that there are faults that transect the Mt. Simon Arkose Zone Sandstone storage reservoir within the modeled CO 2 plume (for more detailed information, see Faults and Fractures section of One Earth Energy Class VI Permit applications). However, the seismic data also shows that there are no faults within the modeled CO 2 plume that transect the confining unit Eau Claire Formation. The faults that transect the Mt. Simon Arkose Zone Sandstone storage reservoir all tip out in the Lower Mt. Simon Formation and do not reach the overlying Eau Claire confining unit. A small 3D survey acquired around the One Earth Energy #1 characterization well confirms these findings.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2026

Prairie State Generating Company 2D Seismic Interpretation

The objectives of the Illinois Storage Corridor (ISC) project are to accelerate commercial deployment of carbon capture utilization and storage at two individual sites and receive approvals for Underground Injection Control (UIC) Class VI permits for construction at each site (ISC Project Narrative, 2020). As part of this project, and as part of the subsurface geologic characterization, two-dimensional (2D) seismic data was acquired at both sites. This report summarizes the findings from two phases of 2D acquisition and seismic interpretation at the Prairie State Generating Company site near Marissa, Illinois. The seismic data indicates the presence of three faults that completely transect the storage and confining units. The initial 2021 2D seismic acquisition revealed the presence of a feature two miles east of the Lively Grove #1 characterization well. The feature is a polyphase fault with a component of strike-slip motion, forming a small positive flower structure. The second phase of 2D seismic acquisition in 2022 constrained the maximum extent of the fault, indicating it has a relatively limited length. Other 2022 seismic lines revealed the presence of two other faults that transect the storage and confining units in the southwest and northeast portions of the 2D seismic acquisition area. The 2021 and 2022 2D seismic surveys have identified three specific locations within the project area that may have an elevated risk of out of zone CO 2 migration due to faulting. Further technical work will be needed to quantify this risk and determine how and if this risk will impact the placement of CO 2 injection wells across the project area.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2026

geoPFA: A Python-Based Open-Source Software for 3D Geothermal PFA

This work presents a novel Python-based framework, geoPFA, for conducting 3D play fairway analysis (PFA) tailored to superhot geothermal systems. The workflow has been applied to the Nesjavellir field in Iceland, a candidate site for the third Iceland Deep Drilling Project's superhot production scenarios. This application demonstrates the value of modular, transparent, and extensible workflows for integrating geological, geophysical, and simulation-derived datasets in high-enthalpy environments. Preliminary results indicate favorable zones consistent with known hydrothermal activity. The geoPFA library will soon be publicly available, offering a scalable and reproducible approach to geothermal exploration across varied geological contexts.

15 GEOTHERMAL ENERGY↗
NASA NTRS1988

Airborne in-situ spectral characterization and concentration estimates of fluorescent organics as a function of depth

The primary purpose of many in-situ airborne light scattering experiments in natural waters is to spectrally characterize the subsurface fluorescent organics and estimate their relative concentrations. This is often done by shining a laser beam into the water and monitoring its subsurface return signal. To do this with the proper interpretation, depth must be taken into account. If one disregards depth dependence when taking such estimates, both their spectral characteristics and their concentrations estimates can be rather ambiguous. A simple airborne lidar configuration is used to detect the subsurface return signal from a particular depth and wavelength. Underwater scatterometer were employed to show that in-situ subsurface organics are very sensitive to depth, but they also require the use of slow moving boats to cover large sample areas. Also, their very entry into the water disturbs the sample it is measuring. The method described is superior and simplest to any employed thus far.

Tittle, R. A.↗
NASA NTRS1990

Glass sample characterization

The development of in-house integrated optical performance modelling capability at MSFC is described. This performance model will take into account the effects of structural and thermal distortions, as well as metrology errors in optical surfaces to predict the performance of large an complex optical systems, such as Advanced X-Ray Astrophysics Facility. The necessary hardware and software were identified to implement an integrated optical performance model. A number of design, development, and testing tasks were supported to identify the debonded mirror pad, and rebuilding of the Technology Mirror Assembly. Over 300 samples of Zerodur were prepared in different sizes and shapes for acid etching, coating, and polishing experiments to characterize the subsurface damage and stresses produced by the grinding and polishing operations.

Ahmad, Anees↗
NASA NTRS1992

Glass sample preparation and performance investigations

This final report details the work performed under this delivery order from April 1991 through April 1992. The currently available capabilities for integrated optical performance modeling at MSFC for large and complex systems such as AXAF were investigated. The Integrated Structural Modeling (ISM) program developed by Boeing for the U.S. Air Force was obtained and installed on two DECstations 5000 at MSFC. The structural, thermal and optical analysis programs available in ISM were evaluated. As part of the optomechanical engineering activities, technical support was provided in the design of support structure, mirror assembly, filter wheel assembly and material selection for the Solar X-ray Imager (SXI) program. As part of the fabrication activities, a large number of zerodur glass samples were prepared in different sizes and shapes for acid etching, coating and polishing experiments to characterize the subsurface damage and stresses produced by the grinding and polishing operations. Various optical components for AXAF video microscope and the x-ray test facility were also fabricated. A number of glass fabrication and test instruments such as a scatter plate interferometer, a gravity feed saw and some phenolic cutting blades were fabricated, integrated and tested.

Johnson, R. Barry↗
NASA NTRS1992

Characterization of Martian near-subsurface materials by determination of cohesion and angle of internal friction

Back-analysis (reconstruction) of the stability of thirty avalanche chutes was performed in the very limited areas where high resolution imaging overlapped with available 1:500 K topographic map coverage. A new technique was developed to incorporate the third dimension (width) of an avalanche chute in stability back-analysis in order to yield unambiguous values of cohesion and angle of internal friction. The procedure is based upon extending the ordinary method of slices to three dimensions, in order to construct avalanche chute cross-sections whose widths and depths vary as a function of gradient, gravity, density of material, and phi and c. Applying the technique to the well documented slide at Lodalen, Norway as a test produces excellent correspondence with reality. Generally, the technique reveals that the width:depth ratio of any avalanche chute decreases with increasing contrast between the average slope angle and the angle of internal friction. Applying this technique to the martian avalanche chute yields results consistent with indications from earlier work, but with greater certainty. Values of cohesion and angle of internal friction identify the materials at the time of failure as moderately cohesive debris. If Sharp's identification of these features as avalanche chutes is correct, then the results here imply that weathering processes have had a significant effect to depths of tens of meters (where failure has occured) below the martian surface. It is also implied that on relatively steep slopes within Valles Marineris, sizable, unaltered, unmantled bedrock exposures for high resolution spectral and spatial scanning by Mars Observer may be scarce.

Sullivan, R. J.↗
NASA NTRS2005

Constraining Oceanic Sources and Sinks of CO2 by Inverse Modeling: A First Step Toward Assimilation of In Situ and Remote Data

The primary goal of this collaborative project, lead by Nicholas Gruber at UCLA, is to better quantify the air-sea fluxes of pre-industrial and anthropogenic carbon dioxide CO2 for the global ocean with a newly developed inverse method. The method involves using ocean general circulation models (OGCMs) to characterize the subsurface distribution of tracers injected for different regions of the ocean surface. By linearly combining a series of these basis (or Greens) functions, one should be able to recreate any arbitrary ocean tracer field. In our case, we are using the global database of ocean dissolved inorganic carbon (DIC) from the WOCE/JGOFS global CO2 survey from the 1990s. The observed DIC field has been partitioned into pre-industrial and anthropogenic components using the empirical C* technique. We then solve for an optimal combination of the regional basis functions that best match the two sets observations in a least-squares sense with singular value decomposition. This in turn provides solutions for the regional air-sea CO2 fluxes (pre-industrial and anthropogenic) that are independent of more traditional estimates derived from surface water disequilibria in the partial pressure of CO2 and wind speed gas exchange parameterizations. We are apply the method to a suite of nine different OGCMs to test the robustness of the derived surface fluxes fields to differences in ocean circulation. Together with a collegue Keith Lindsay at NCAR, S. Doney's role in this project was to create the appropriate steady-state and temporally evolving basis functions from the NCAR CSM-1 ocean circulation model. Because of a slow start across the project in defining a common set of experimental protocols, we requested and were given a one year, no cost extension. Both the steady-state and temporally evolving basis functions were delivered to the UCLA group for further analysis. Additional information is included in the original extended abstract.

Doney, Scott↗
NASA NTRSDate not supplied

Interpretable Machine Learning Models for Autonomous Characterization of Analogue Ocean World Seawater Chemistry and Biosignature Potential Using Isotope Ratio Data

Background: Future missions to ocean worlds, such as Enceladus and Europa, will attempt to characterize the subsurface seawater chemistry and assess the potential for life. Such missions will be equipped with capabilities to precisely measure volatile isotopes in plumes, atmospheres, and exospheres. Motivation: While large isotopic fractionations can indicate a biological source, there are signatures resulting from abiotic geochemical processes that mimic isotopic biosignatures. While machine learning (ML) has the potential to disentangle competing effects and biotic mimicry, high-dimensional isotope ratio mass spectrometry (IRMS) data is likely to contain noise/irrelevant features and involve complex statistical interactions that make human inference and interpretation difficult. Further, ML predictions with as far-reaching implications as an extraterrestrial biosignature on an ocean world requires the use of interpretable models (i.e., not “black box” models) with physically and mathematically meaningful feature spaces along with false positive diagnostics. Methods: We use volatile CO2 IRMS data of analogue ocean world seawaters to validate an ML approach to provide biogeochemical context for biosignature detection. We employ a feature selection method called nearest-neighbor projected distance regression (NPDR) that detects statistical interactions and helps elucidate the mechanisms of the Random Forest classification models. Results: We train and validate predictive ML models on volatile CO2 IRMS data of analogue ocean world seawaters to predict major salt components (e.g., MgSO4, NaHCO3), pH, ionic strength, and the presence of biosignatures. Features derived from IRMS measurements are augmented with extracted time-series features. Our results show high test accuracy and interpretability, which is increased by interaction network visualization, sample-wise variable importance scores, and single-sample class probability estimates. We demonstrate an ML mission software solution that triggers autonomous data transmission and biogeochemical sample prediction.

geochemistry↗
DOE OSTI2024

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University: Preprint

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ("CUBO") was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗