Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Subsurface data models”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Three-dimensional seismic characterization and imaging of the Soda Lake geothermal field

Accurate characterization of subsurface geophysical properties and detection of the fault system are essential for geothermal energy exploration and production. The Soda Lake geothermal field is in western Nevada with a complex fault system. Previous seismic characterization only produced a low-resolution, smooth velocity model along with a simple, conceptual fault model. Using optimized correlation-based full-waveform inversion, wavefield-separation-based reverse-time migration, and automatic fault detection techniques, we present 3D seismic characterization for the Soda Lake geothermal field using 3D surface seismic data acquired with Vibroseis sources. Here, we obtain 3D high-resolution velocity, density, and acoustic impedance models, 3D seismic images with different grid spacings, and a high-resolution fault system. Consistency check between the constructed faults and currently active injection and production geothermal wells verifies that our seismic inversion and imaging results and detected faults are reliable. These results can provide valuable information for optimizing well placement and geothermal energy production at the Soda Lake geothermal field.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Stratigraphy‐Induced Localization of Microseismicity During CO 2 Injection in Illinois Basin

Abstract Subsurface fluid injection stimulates complex hydromechanical interaction, necessitating the integration of geomechanical data across spatial and temporal scales to consider the sophisticated behavior. Induced seismic response is usually associated with the complex reservoir architecture and pre‐existing features that are three‐dimensional, such as local stratigraphy, fractures, faults, and other discontinuities. This study encompasses laboratory characterization of the coupled hydromechanical response of cores extracted from rock formations in Illinois Basin: reservoir ‐ Mt. Simon sandstone, basal seal ‐ Argenta sandstone, and crystalline basement ‐ Precambrian rhyolite. High‐resolution numerical modeling allows considering the three‐dimensional complexity of the Illinois Basin Decatur Project with spatial resolution comparable to one of the active seismic surveys. A detailed reconstruction of the evolving state of stress in formations lacking direct stress measurements is achieved by numerical modeling that integrated laboratory‐derived hydromechanical properties, a porosity‐permeability relationship, active seismic data, and an inverted three‐dimensional porosity distribution. It appears that the microseismic clusters, mainly observed in the crystalline basement during the injection, are linked to zones experiencing more critically stressed conditions prior to injection. These zones have a potential for reactivation during the injection and are attributed to the specific local stratigraphy of the injection site, as well as transfer of triggering perturbations during the injection.

Bondarenko, N. [University of Illinois Urbana‐Cham↗
DOE OSTI2023

Complete Performance Comparison Between the Optimized Image Construction Algorithm (U-MBIR) and the Existing Reconstruction Algorithm for Detecting Defects and Damage in Concrete

Reinforced concrete (RC) is a composite material subjected to mechanical, thermal, and chemical loads throughout its service life. Because of these external stressors and the susceptibility of RC structural members to shrinkage and microcracking, the material degrades throughout its life cycle. This deterioration can lead to a decrease in member capacity and, ultimately, poses a threat to the structural integrity. Thus, it is crucial that the damage caused by aging and degradation be monitored and assessed at regular intervals throughout the material’s service life. Since coring of the material is typically not feasible for in-service structural systems, non-destructive evaluation (NDE) methodologies are used to assess remaining structural capacity. NDE methods enable surface and subsurface examination without damaging or degrading the medium. Moreover, RC is a critical component of nuclear power plants; thus, its safety and reliability must be thoroughly examined throughout the life cycle of the structural system. Ultrasonic measurements have been an industry standard for both surface and subsurface inspections. To this end, Oak Ridge National Laboratory (ORNL) has researched and developed advanced image reconstruction algorithms to capture internal damage. The results and discussion presented herein summarize the current state of the ultrasonic model–based iterative reconstruction (U-MBIR) algorithm developed at ORNL. More specifically, this report presents a comparison between reconstruction images produced via a widely employed ultrasonic NDE technique—the synthetic aperture focusing technique (SAFT)—and the ORNL-developed U-MBIR algorithm. These NDE methodologies are demonstrated using ultrasonic data collected from four concrete specimens. Overall, the U-MBIR algorithm eliminates artifacts and noise that are typically present within the SAFT reconstructions, and it shows defects and anomalies more clearly than the SAFT images. In conclusion, this algorithm is suitable for identifying concrete defects, although more improvements and optimization could be done to better define internal defects.

36 MATERIALS SCIENCE↗
NASA NTRS1978

Noble gases in Luna 24 core soils

Isotopic abundances of the noble gases have been determined in grain size separates of six soils from different depths in the Luna 24 core. Concentrations of cosmogenic Ne-21 in these soils show a complex pattern as a function of soil depth, and are inconsistent with a single stage, in situ production in a static regolith or with a simple model of slow accumulation of previously unirradiated soils followed by a static period. All soils have apparently had a complex pre-irradiation history. Some soils show significant differences in the concentration of cosmogenic Ne-21 between different grain sizes, which is consistent with suggestions that they formed as a mixture of soils with different maturities. The exact subsurface depth relationships of the soils are not known, and several possibilities are discussed. Because of this factor and the complex nature of the soils, the cosmogenic Ne data do not define the depositional-irradiational history of the core. Luna 24 core soils may have been deposited any time during the past 300 million years.

Bogard, D. D.↗
NASA NTRS1993

Observations of time variation in the sun's rotation

Observations of solar p-mode frequency splittings obtained at Big Bear Solar Observatory in 1986 and during 1988-90 reveal small (about 1 percent) changes in the sun's subsurface angular velocity with solar cycle. An asymptotic inversion of the splitting data yields the latitude dependence of the rotation rate and shows that the largest changes in the angular velocity, about 4 nanoHz, occurred between 1986 and the later years, at high (about 60 deg) solar latitudes. Earlier helioseismic observations suggest that solar cycle changes in the ratio of magnetic to turbulent pressure in the solar convection zone are large enough to account for the magnitude of the observed angular velocity variations, but a detailed model of the phenomenon does not exist.

Woodard, M. F.↗
DOE OSTI2021

Stepwise Dynamic Calibration of a Hydromechanical Simulation Using Time-Lapse Vertical Seismic Profile

This study aims to develop a methodology for calibrating subsurface stress changes through time-lapse Vertical Seismic Profiling (VSP) integration. The selected study site is the 13-10A injector well within the ongoing CO2-EOR operation of the Farnsworth Field Unit. The Time-lapse VSP dataset carries the combined effects of fluid substitution and mean effective stress changes, thereby providing a dataset amenable for the calibration of production and injection-induced stress changes. The concept is similar to calibrating a reservoir simulation model in that the process honor real field data to set up an inverse problem. The solution optimizes the independent and impactful geomechanical parameters that replicate the observed time-lapse seismic velocity changes. This stress calibration is enabled by 4D geomechanical modeling and the VSP Integration workflow. This calibration benefits from extensive geological, geophysical and geomechanical characterization through 3D seismic data, geophysical well logs, and core assessed as part of the 1D MEM conducted on the 13-10A subject well. These data are used to develop a site-specific rock physics model. The Biot Gassmann workflow combines rock physics and reservoir simulation outputs to determine the fluid substitution contribution to seismic velocity change. Additionally, modeled seismic velocity attributed to mean effective stress are determined from the geomechanical simulation outputs, and the stress-velocity relationship developed from the ultrasonic seismic velocity measurements on the extracted Morrow B core. A penalty function is then formed between the modeled seismic velocities and the observed time-lapse VSP dataset. Four independent and impactful geomechanical parameters have been determined. These are the bulk modulus and shear modulus for zero porosity and the shear and compressional seismic velocity to mean effective stress derivatives. The dataset of numerous coupled hydromechanical- geomechanical simulation realizations is built by combining variations of the four stated geomechanical parameters. A machine learning-assisted workflow comprised of an artificial neural network and a particle swarm optimizer are used to converge on the optimal geomechanical parameters. The successful execution of this workflow has affirmed the suitability of acoustic time-lapse measurements for 4D-VSP geomechanical stress calibration pending measurable stress sensitivities within the anticipated effective stress changes and the availability of suitable and reliable datasets for petroelastic modeling.

02 PETROLEUM↗
NASA NTRS2003

Reliability of CHAMP Anomaly Continuations

CHAMP is recording state-of-the-art magnetic and gravity field observations at altitudes ranging over roughly 300 - 550 km. However, anomaly continuation is severely limited by the non-uniqueness of the process and satellite anomaly errors. Indeed, our numerical anomaly simulations from satellite to airborne altitudes show that effective downward continuations of the CHAMP data are restricted to within approximately 50 km of the observation altitudes while upward continuations can be effective over a somewhat larger altitude range. The great unreliability of downward continuation requires that the satellite geopotential observations must be analyzed at satellite altitudes if the anomaly details are to be exploited most fully. Given current anomaly error levels, joint inversion of satellite and near- surface anomalies is the best approach for implementing satellite geopotential observations for subsurface studies. We demonstrate the power of this approach using a crustal model constrained by joint inversions of near-surface and satellite magnetic and gravity observations for Maude Rise, Antarctica, in the southwestern Indian Ocean. Our modeling suggests that the dominant satellite altitude magnetic anomalies are produced by crustal thickness variations and remanent magnetization of the normal polarity Cretaceous Quiet Zone.

vonFrese, Ralph R. B.↗
DOE OSTI2021

Multichannel Analysis of Surface Waves Accelerated (MASWAccelerated): Software for efficient surface wave inversion using MPI and GPUs

Multichannel Analysis of Surface Waves (MASW) is a technique frequently used in geotechnical engineering and engineering geophysics to infer 1D layered models of seismic shear wave velocities in the top tens to hundreds of meters of the subsurface. We aim to accelerate MASW calculations by capitalizing on modern computer hardware available in the workstations of most engineers: multiple cores and graphics processing units (GPUs). We propose new parallel and GPU accelerated algorithms for computing 1D MASW inversion, and provide software implementations in C using Message Passing Interface (MPI) and CUDA. These algorithms take advantage of sparsity that arises in the problem, and the work balance between processes considers typical data trends. We compare our methods to an existing open source Matlab MASW tool. Our serial C implementation achieves a 2x speedup over the Matlab software, and we continue to see improvements by parallelizing the problem with MPI. Here we see nearly perfect strong and weak scaling for uniform data, and improve strong scaling for realistic data by repartitioning the problem to process mapping. By utilizing GPUs available on most modern workstations, we observe an additional 1.3x speedup over the serial C implementation on the first use of the method. We typically repeatedly evaluate theoretical dispersion curves as part of an optimization procedure, and on the GPU the kernel can be cached for faster reuse on later runs. We observe a 3.2x speedup on the cached GPU runs compared to the serial C runs. This work is the first open-source parallel or GPU-accelerated software tool for MASW imaging, and should enable geotechnical engineers to fully utilize all computer hardware at their disposal.

58 GEOSCIENCES↗
DOE OSTI2025

Advancing the Representation of Human Actions in Large‐Scale Hydrological Models: Challenges and Future Research Directions

Characterizing the impact of human actions on terrestrial water fluxes and storages at multi-basin, continental, and global scales has long been on the agenda of scientists engaged in climate science, hydrology, and water resources systems analysis. This need has resulted in a variety of modeling efforts focused on the representation of water infrastructure operations. Yet, the representation of human-water interactions in large-scale hydrological models is still relatively crude, fragmented across models, and often achieved at coarse resolutions (~10–100 km) that cannot capture local water management decisions. In this commentary, we argue that the concomitance of four drivers and innovations is poised to change the status quo: “hyper-resolution” hydrological models (~0.1–1 km), multi-sector modeling, satellite missions able to monitor the outcome of human actions, and machine learning are creating a fertile environment for human-water research to flourish. We then outline four challenges that chart future research in hydrological modeling: (a) creating hyper-resolution global data sets of water management practices, (b) improving the characterization of anthropogenic interventions on water quantity, stream temperature, and sediment transport, (c) improving model calibration and diagnostic evaluation, and (d) reducing the computational requirements associated with the successful exploration of these challenges. Overcoming them will require addressing modeling, computational, and data development needs that cut across the hydrology community, thereby requiring a major communal effort.

catchment hydrology↗
DOE OSTI2024

A new perspective on density and strength loss profiles at the surface of thermally oxidized nuclear graphite

Oxidation of graphite components could influence their designed life in a high-temperature nuclear reactor. The oxidized regions could potentially lower the allowed stress capacity. The American Society of Mechanical Engineers rules for the design and construction of graphite-moderated reactors recommend that subsurface regions that might become excessively damaged by oxidation during reactor operation be identified and excluded from geometry and stress calculations. Identification of oxidation-affected regions is possible, in principle, through complex modeling exercises of reactor behavior during hypothetical accident scenarios coupled with graphite oxidation models, but this procedure may not have the precision needed for informed decisions. Here, this paper proposes an alternate method, based on interpretation of a series of well-designed oxidation experiments, which could augment the designer's tools. The procedure is illustrated by data on oxidation by air of several graphite grades (NBG-18, PCEA, IG-110, R4-650) that are corroborated with independent literature information, when available. The Wichner model for graphite oxidation used for this analysis provides conservative results that could be quickly implemented in the design process.

36 MATERIALS SCIENCE↗
NASA NTRS2020

A Fresh Look at the Stratigraphy of Northern Australe

The roughly circular collection of mare deposits centered at ~38.9°S, 93°E is often re- ferred to as Mare Australe. It is located outside of the Procellarum KREEP Terrain. The circular arrangement of Australe’s mare patches has suggested an ancient, heavily degraded or relaxed impact basin roughly 900 km in diameter. The mare deposits are generally thought to have erupted into smaller post- basin craters. The type, volume, and distribution of mare eruptions potentially resembles the early stages of basin-filling mare events, but which are preserved in Australe and some farside locations. Gravity data suggest that if there was a basin, it is much smaller than originally proposed (now ~600 km) and located in the northern part of Mare Australe, between Humboldt, Milne, and Jenner craters. As a whole, Mare Australe lacks the topography typically associated with a basin; however, northern Australe has a slight topographic depression that roughly corresponds to the basin-like Bouguer gravity signature in the same area. The compositions exposed in Humboldt crater suggest that a preexisting basin might have excavated deeper crustal material. However, the underlying cause of the circularity of Mare Australe’s deposits, particularly those extending outside of the potential impact basin setting, is not yet understood. Thus, Australe may preserve fundamental information about mare volcanism potentially uncoupled from basin formation and structure. The objectives of this study are to use new high- resolution data (images, gravity, topography, and com- position) to reassess Australe’s mare deposits, deter- mine the timing and style of volcanism, identify discrete basalt deposits, and to further characterize the evolution of magmatism and subsurface structure in this area. Here, we focus on the northern Australe deposits (between Humboldt, Jenner, and Milne). As originally noted by Whitford-Stark (1979), Humboldt crater and its ejecta make an excellent stratigraphic marker that can be traced across much of the Australe region. The ejecta serves as a stratigraphic constraint for absolute model ages (AMAs) derived from crater size-frequency distributions (CSFDs).

Stopar, J. D.↗
NASA NTRS1993

A new MAP for Mars

A Mars Aerial Platform (MAP) mission capable of generating thousands of very-high-resolution (20 cm/pixel) pictures of the Martian surface is considered. The MAP entry vehicle will map the global circulation of the planet's atmosphere and examine the surface and subsurface. Data acquisition will use instruments carried aboard balloons flying at nominal altitude of about 7 km over the Martian surface. The MAP balloons will take high- and medium-resolution photographs of Mars, sound its surface with radar, and provide tracking data to chart its winds. Mars vehicle design is based on the fourth-generation NTP, NEP, SEP vehicle set that provides a solid database for determining transportation system costs. Interference analysis and 3D image generation are performed using manual system sizing and sketching in conjunction with precise CAD modeling.

Zubrin, Robert↗
NASA NTRS2003

Polar and Non-Polar Layers on Mars: A Single Mechanism for Formation?

The recent discovery of vast quantities of near-subsurface ice in both polar regions of Mars by the Mars Odyssey Gamma Ray Spectrometer (GRS) has presented us with an interesting quandary. On one hand, these deposits, found poleward of 60 deg in both hemispheres, are consistent with thermal models suggesting ice will be best protected in these regions during periods of high obliquity. On the other hand, the current paradigm regarding the placement of these deposits, i.e., diffusive deposition of water vapor, appears to be inconsistent with the large volume mixing ratios (approx. 90%) inferred from the GRS data. This incongruity argues that diffusion alone cannot be the primary mechanism for the creation of these reservoirs, and that an alternate, large-scale process should be considered.

Mischna, M. A.↗
NASA NTRS2017

Connecting Projects to Complete the In Situ Resource Utilization Paradigm

Terrain Identify specifics such as slope, rockiness, traction parameters Identify what part of ISRU needs each Physical Geotechnical Hardness, density, cohesion, etc. Identify what part of ISRU needs each (e.g., excavation needs to know hardness, density; soil processing needs to know density, cohesion; etc.)Mineral Identify specifics Identify what part of ISRU needs each Volatile Identify specifics Identify what part of ISRU needs each Atmosphere Identify specifics Identify what part of ISRU needs each Environment Identify specifics Identify what part of ISRU needs each Resource Characterization What: Develop an instrument suite to locate and evaluate the physical, mineral, and volatile resources at the lunar poles Neutron Spectrometer Near Infrared (IR) to locate subsurface hydrogen surface water Near IR for mineral identification Auger drill for sample removal down to 1 m Oven with Gas Chromatograph Mass Spectrometer to quantify volatiles present ISRU relevance: Water volatile resource characterization and subsurface material access removal Site Evaluation Resource Mapping What: Develop and utilize new data products and tools for evaluating potential exploration sites for selection and overlay mission data to map terrain, environment, and resource information e.g., New techniques applied to generate Digital Elevation Map (DEMs) at native scale of images (1mpxl)ISRU relevance: Resource mapping and estimation with terrain and environment information is needed for extraction planning Mission Planning and Operations What: Develop and utilize tools and procedures for planning mission operations and real time changes Planning tools include detailed engineering models (e.g., power and data) of surface segment systems allows evaluation of designs ISRU relevance: Allows for iterative engineering as a function of environment and hardware performance.

Paradigm↗
DOE OSTI2023

Harnessing Uncertainty through Functional Data Analysis in Gas Breakthrough Data

Detecting subsurface explosions from radionuclide gas migration through rock fractures is an effective way to identify nuclear activity. Los Alamos National Laboratory (LANL) has developed simulation methods, based on data from the 1962 Hardhat underground nuclear test, to predict gas breakthrough times at the surface. However, these methods rely on an imperfect understanding of the relationship between rock damage and fracture permeability. Our clinic project studies methods for predicting breakthrough curves that characterize total mass produced as a function of time, as well as quantifying the uncertainty associated with these predictions. The model that is currently employed to relate damage to permeability uses an empirically motivated power-law expression, with a range of parameter values that are compatible with the experimental Hardhat data. We develop emulators, built from functional data analysis techniques and trained on simulation data, that rapidly predict the gas breakthrough curve given a damage field and given the parameter values of the power-law equation. Using Bayesian regression, we address the problem of uncertainty quantification in our emulators. Finally, in order to test the robustness of the model, we further validate it on a damage field representing different physical conditions.

58 GEOSCIENCES↗
DOE OSTI2022

Flexible and Accessible 4D Subsurface Visualization Using a Web-Based Platform

4D subsurface visualization using a web platform can provide improved communication, education and outreach to non-experts and stakeholders. It allows for the improved understanding of complex relationships and interactions that occur in inaccessible locations of which many experimental testbeds exist in. The use of a web-based visualization tool lessens the cognitive impact further by reducing the need of specialized software. Users can define their own exploration of a 3D or 4D scene, adding and removing data as needed and are able to make informed decisions based on accurate model visualizations, while providing a collaborative tool that can be accessed on any platform using a browser window. The tool is also able to point to real time streams to display up-to-the-second data as an experiment is ongoing. This level of latency can also provide operators with essential information that can direct an experiment’s progress. Additionally, this tool can leverage augmented and virtual reality (AR, VR) capabilities of certain mobile devices and head-mounted displays, providing further engaging visualization possibilities.

Pratt, Martin J.↗
NASA NTRS1983

The stability of ground ice in the equatorial region of Mars

The lifetime of an unreplenished layer of ground ice lying within 30 deg of the Martian equator was examined within the context of the existing data base on Martian regolith and climate. Data on the partial pressure of H2O in the Martian atmosphere and the range of mean annual temperatures indicated the ground ice would be restricted to latitudes poleward of 40 deg. However, the ground ice near the poles may be a relic from early Martian geologic times held in place by a thin layer of regolith. Consideration of twelve model pore size distributions, similar to silt- and clay-type earth soils, was combined with a parallel pore model of gaseous diffusion to calculate the flux of H2O molecules escaping from the subsurface ground ice layer. Martian equatorial ground ice was found to be influenced by the soil structure, the magnitude of the geothermal gradient, the climatic desorption of CO2 from the regolith. It is concluded that equatorial ground ice is present on Mars only if a process of replenishment is active.

Clifford, S. M.↗
DOE OSTI2020

Overland flow numerical model prediction, Lower Triangle Region in East River Watershed, Colorado, 3 days

This data package contains numerical simulation results of surface flow variables such as flow velocity and water depth in Lower Triangle Region in East River Watershed, Colorado. The surface flow is a consequence of a high intensity rainfall event with a total duration of 3 days, available at a resolution of 10 minutes. The results are computed on triangular multiresolution meshes with resolutions ranging from 10 meter to 80 meter. The data package also contains a simulation on a uniform triangular mesh with a resolution of 10 meter. The simulations consider surface flow only and neglect subsurface flow, infiltration, and evapotranspiration. The purpose of the data is to assess the quality of a mesh refinement strategy.

54 ENVIRONMENTAL SCIENCES↗