Search NASA⌕ Search

SEARCH · Search NASA

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 127 records · Page 7

Survey Plan For Characterization of the Subsurface Underlying the National Aeronautics and Space Administration's Marshall Space Flight Center in Huntsville, Alabama

Topic considered include: survey objectives; technologies for non-Invasive imaging of subsurface; cost; data requirements and sources; climatic condition; hydrology and geology; chemicals; magnetometry; electrical(resistivity, potential); optical-style imaging; reflection/refraction seismics; gravitometry; photo-acoustic activation;well drilling and borehole analysis; comparative assessment matrix; ground sensors; choice of the neutron sources; logistic of operations; system requirements; health and safety plans.

Source record↗

Characterization of a Subsurface Biosphere in a Massive Sulfide Deposits at Rio Tinto, Spain: Implications for Extant Life on Mars

The recent discovery of abundant sulfate minerals, particularly Jarosite by the Opportunity Rover at Sinus Merdiani on Mars has been interpreted as evidence for an acidic lake or sea on ancient Mars [1,2], since the mineral Jarosite is soluble in liquid water at pH above 4. The most likely mechanism to produce sufficient protons to acidify a large body of liquid water is near surface oxidation of pyrite rich deposits [3]. The acidic waters of the Rio Tinto, and the associated deposits of Hematite, Goethite, and Jarosite have been recognized as an important chemical analog to the Sinus Merdiani site on Mars [4]. The Rio Tinto is a river in southern Spain that flows 100 km from its source in the Iberian pyrite belt, one of the Earth s largest Volcanically Hosted Massive Sulfide (VHMS) provinces, into the Atlantic ocean. The river originates in artesian springs emanating from ground water that is acidified by the interaction with subsurface pyrite ore deposits. The Mars Analog Rio Tinto Experiment (MARTE) has been investigating the hypothesis that a subsurface biosphere exists at Rio Tinto living within the VHMS deposit living on chemical energy derived from sulfur and iron minerals. Reduced iron and sulfur might provide electron donors for microbial metabolism while in situ oxidized iron or oxidants entrained in recharge water might provide electron acceptors.

Stoker, C. R.↗

Characterization of a Subsurface Biosphere in a Massive Sulfide Deposit At Rio Tinto, Spain: Implications For Extant Life On Mars

The recent discovery of abundant sulfate minerals, particularly Jarosite by the Opportunity Rover at Sinus Merdiani on Mars has been interpreted as evidence for an acidic lake or sea on ancient Mars [1,2], since the mineral Jarosite is soluble in liquid water at pH above 4. The most likely mechanism to produce sufficient protons to acidify a large body of liquid water is near surface oxidation of pyrite rich deposits [3]. The acidic waters of the Rio Tinto, and the associated deposits of Hematite, Goethite, and Jarosite have been recognized as an important chemical analog to the Sinus Merdiani site on Mars [4]. The Rio Tinto is a river in southern Spain that flows 100 km from its source in the Iberian pyrite belt, one of the Earth's largest Volcanically Hosted Massive Sulfide (VHMS) provinces, into the Atlantic ocean. The river originates in artesian springs emanating from ground water that is acidified by the interaction with subsurface pyrite ore deposits. The Mars Analog Rio Tinto Experiment (MARTE) has been investigating the hypothesis that a subsurface biosphere exists at Rio Tinto living within the VHMS deposit living on chemical energy derived from sulfur and iron minerals. Reduced iron and sulfur might provide electron donors for microbial metabolism while in situ oxidized iron or oxidants entrained in recharge water might provide electron acceptors.

Stoker, C. R.↗

Forward and inverse modeling of fault transmissibility in subsurface flows

Characterizing physical properties of faults, such as their transmissibility, is crucial for performing predictive numerical simulation of subsurface flows, such as those encountered in petroleum engineering and remediation of subsurface contamination. Here, this paper provides a complete investigation of the inverse problem for fault transmissibility in subsurface flow models, under appropriate assumptions on fault structure. In particular, the following aspects are considered: 1) fault modeling and well-posedness of the forward problem; 2) finite element (FEM) discretizations of the forward problem and their rigorous a priori convergence analysis; 3) Well-posedness of the Bayesian inverse problem, FEM discretization of the infinite dimensional Bayesian inverse formulation, and its rigorous a priori analysis. Moreover, computation of the maximum a posteriori (MAP) point via fast inexact Newton-conjugate gradient optimization and a Laplace approximation of the Bayesian posterior are also presented. Numerical results illustrate the use of the proposed fault model in forward and inverse problems for subsurface flows in two dimensional domains with multiple faults.

97 MATHEMATICS AND COMPUTING↗

Geology of the One Earth Energy Site

The One Earth Energy site is one of two sites in the Illinois Storage Corridor (ISC) project. The objectives of the ISC project is 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. At the One Earth Energy site, an extensive data collection program was undertaken, which included the drilling of a test well (One Earth Energy #1 [OEE #1]), four 2D seismic lines, and a small 3D seismic survey. The OEE #1 well was drilled in 2022 and acquired extensive core, log, and testing data to characterize the subsurface geology of the site. Coring was focused on the storage interval, the Mt. Simon Sandstone, and the confining interval, the Eau Claire Formation. The core and log data were used to evaluate the sedimentology and sequence stratigraphy, as well as to develop the conceptual geologic model. This report includes the geological summaries of the Mt. Simon Sandstone and the Eau Claire Formation. The extensive analysis of the log data is included in the petrophysical section, showing ranges of porosity, estimated pore size, and the mineral content of selected zones in the well. The separate petrographic technical report entitled “Petrographic and Advanced Geologic Characterization Report on One Earth Energy #1 (API# 1211325373)”, report number DOE-UIUC-0031892-04, details thin section point-counting analysis that includes mineralogical and pore space analysis, including grain size analysis, annotated thin section photomicrographs, scanning electron microscopy (SEM) with energy dispersive X-ray spectroscopy (EDS), and statistics of grain size analysis on Mt. Simon thin sections from OEE #1. The final OEE #1 well data to be included in this geology report is the routine core analysis of both whole core plugs and rotary sidewall core plugs. In addition to the OEE #1 well, four 2D seismic lines and a small 3D survey were acquired as part of the overall subsurface geological characterization. This geology report references the seismic interpretation report, entitled “One Earth Energy Site Seismic Interpretation Task 5.0”, report number DOE-UIUC-0031892-07. This report details the stratigraphic and structural interpretation of the 2D and 3D seismic data acquired at the One Earth Energy site. The 2D seismic data was acquired in 2019 and 2021, and the 3D survey was acquired in 2022. The objectives of the seismic programs were to contribute to the subsurface characterization of the Mt. Simon-Eau Claire Storage Complex by evaluating the continuity of potential storage reservoirs and containment intervals across the project area, and to determine if any geologic features are present that would increase containment risk to the proposed carbon storage project.

09 BIOMASS FUELS↗

3D Seismic Experimentation and Advanced Processing/Inversion Development for Investigations of the Shallow Subsurface

The work plan for this project was to develop and apply advanced seismic reflection and wide-angle processing and inversion techniques to high resolution seismic data for the shallow subsurface to seismically characterize the shallow subsurface at hazardous waste sites as an aid to containment and cleanup activities. We proposed to continue work on seismic data that we had already acquired under a previous DoE grant, as well as to acquire additional new datasets for analysis. The project successfully developed and/or implemented the use of 3D reflection seismology algorithms, waveform tomography and finite-frequency tomography using compressional and shear waves for high resolution characterization of the shallow subsurface at two waste sites. These two sites have markedly different near-surface structures, groundwater flow patterns, and hazardous waste problems. This is documented in the list of refereed documents, conference proceedings, and Rice graduate theses, listed below.

54 ENVIRONMENTAL SCIENCES↗

Applications of Ensemble-based Data Assimilation Techniques for Aquifer Characterization using Tracer Data at Hanford 300 Area

Subsurface aquifer characterization often involves high parameter dimensionality and requires tremendous computational resources if employing a full Bayesian approach. Ensemble-based data assimilation techniques, including filtering and smoothing, are computationally efficient alternatives. Despite the increasing number of applications of ensemble-based methods in assimilating flow and transport related data for subsurface aquifer charaterization, most are limited to either synthetic studies or two-dimensional problems. In this study, we applied ensemble-based techniques for assimilating field tracer experimental data obtained from the Integrated Field Research Challenge (IFRC) site at the Hanford 300 Area. The forward problem was simulated using the massively-parallel three-dimensional flow and transport code PFLOTRAN to effectively deal with the highly transient flow boundary conditions at the site and to meet the computational demands of ensemble-based methods. This study demonstrates the effectiveness of ensemble-based methods for characterizing a heterogeneous aquifer by sequentially assimilating multiple types of data. The necessity of employing high performance computing is shown to enable increasingly mechanistic non-linear forward simulations to be performed within the data assimilation framework for a complex system with reasonable turnaround time.

data assimilation, ensemble kalman filter, inverse↗

Consensus Equilibrium for Subsurface Delineation

Heterogeneity and insufficient site characterization limit our knowledge of the subsurface. Inversion techniques, which minimize the mismatch between observations and model predictions, have become an essential tool of subsurface characterization. Most optimization-based approaches fail to incorporate various implicit priors and capture the geological complexity. We overcome these limitations by deploying the plug-and-play and consensus equilibrium (CE) strategies, which provide a flexible framework for image reconstruction. Our CE methodology for spatial delineation of geologic formations consists of an image denoiser and a variational auto-encoder (deep learning-based emulator). The former ameliorates the reconstruction noise, yielding well-defined geological structures; its mathematical equivalence with the proximal operator allows the deployment of advanced denoisers (e.g., CNN-based denoiser) that do not correspond to a regularization objective. The latter defines a geology prior that imposes a geological constraint, for example, continuity and shape of geological features, onto the reconstructed image. Here, we conduct a series of numerical experiments dealing with transient two-dimensional flow driven by a pumping well and natural hydraulic head gradient. They demonstrate the CE framework's ability to delineate, both probabilistically and deterministically, complex subsurface environments with sufficient quality.

42 ENGINEERING↗

Sequential geophysical and flow inversion to characterize fracture networks in subsurface systems

Subsurface applications, including geothermal, geological carbon sequestration, and oil and gas, typically involve maximizing either the extraction of energy or the storage of fluids. Fractures form the main pathways for flow in these systems, and locating these fractures is critical for predicting flow. However, fracture characterization is a highly uncertain process, and data from multiple sources, such as flow and geophysical are needed to reduce this uncertainty. We present a nonintrusive, sequential inversion framework for integrating data from geophysical and flow sources to constrain fracture networks in the subsurface. In this framework, we first estimate bounds on the statistics for the fracture orientations using microseismic data. These bounds are estimated through a combination of a focal mechanism (physics-based approach) and clustering analysis (statistical approach) of seismic data. Then, the fracture lengths are constrained using flow data. In conclusion, the efficacy of this inversion is demonstrated through a representative example.

97 MATHEMATICS AND COMPUTING↗

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗

Subsurface microbial communities as a tool for characterizing regional-scale groundwater flow

Subsurface microbial community distribution patterns are influenced by biogeochemical and groundwater fluxes and may inform hydraulic connections along groundwater-flow paths. This study examined the regional-scale microbial community of the Death Valley Regional Flow System and evaluated whether subsurface communities can be used to identify groundwater-flow paths between recharge and discharge areas. Samples were collected from 36 sites in three groundwater basins: Pahute Mesa–Oasis Valley (PMOV), Ash Meadows (AM), and Alkali Flat–Furnace Creek Ranch (AFFCR). Microbial diversity within and between communities varied by location, and communities were separated into two overall groups that affiliated with the AM and PMOV/AFFCR basins. Network analysis revealed patterns between clusters of common microbes that represented groundwaters with similar geochemical conditions and largely corroborated hydraulic connections between recharge and discharge areas. Null model analyses identified deterministic and stochastic ecological processes contributing to microbial community assemblages. Most communities were more different than expected and governed by dispersal limitation, geochemical differences, or undominating processes. However, certain communities from sites located within or near the Nevada National Security Site were more similar than expected and dominated by homogeneous dispersal or selection. Overall, the (dis)similarities between the microbial communities of DVRFS recharge and discharge areas supported previously documented hydraulic connections between: (1) Spring Mountains and Ash Meadows; (2) Frenchman and Yucca Flat and Amargosa Desert; and (3) Amargosa Desert and Death Valley. However, only a portion of the flow path between Pahute Mesa and Oasis Valley could be supported by microbial community analyses, likely due to well-associated artifacts in samples from the two Oasis Valley sites. This study demonstrates the utility of combining microbial data with hydrologic, geologic, and water-chemistry information to comprehensively characterize groundwater systems, highlighting both strengths and limitations of this approach.

54 ENVIRONMENTAL SCIENCES↗

Geoelectric Characterization of Hyporheic Exchange Flow in the Bedrock‐Lined Streambed of East Fork Poplar Creek, Oak Ridge, Tennessee

Abstract A multimethod geoelectric survey was implemented between January and March 2022 along a 220‐m long reach of the bedrock‐lined streambed of East Fork Poplar Creek in Oak Ridge, Tennessee to identify locations of surface‐water and groundwater exchange and characterize the subsurface flow paths that convey water between the stream and flood plain. A waterborne self‐potential (WaSP) survey was completed in January 2022 to measure the electric streaming‐potential field in the stream. Electric resistivity tomography (ERT) was performed in March 2022 on the flood plain adjacent to the WaSP survey reach to map the electric resistivity distribution and characterize the hydrogeology and subsurface flow paths that facilitate surface‐water and groundwater exchange in the bedrock‐lined stream. The combination of WaSP and ERT data support the qualitative interpretation that surface‐water and groundwater exchange likely occurs along fractures in outcropping bedrock and along two fault lines that intersect the limestone creek bed.

54 ENVIRONMENTAL SCIENCES↗

Technical Report on Subsurface Monitoring of the Brady Hot Spring Geothermal Site, Nevada, based upon Full Waveform Inversion

Abilities to accurately characterize the subsurface in a geothermal setting is key to assess and support production. An important element of geothermal reservoir monitoring is also the ability to investigate fluid transport within fracture network. This report focuses on improving subsurface imaging and monitoring in geothermal settings using full waveform inversion based on the adjoint method and time-lapse imaging. To assess our method, we rely on a dense seismic dataset collected in 2016 at the Brady Hot Springs geothermal site in Nevada for the DOE-funded project Poroelastic Tomography by Adjoint Inverse Modeling of Data from Seismology, Geodesy, and Hydrology. This dataset captures subsurface changes across four stages of geothermal power plant operations, which involve varying rates of fluid injection and extraction. Two velocity models were previously derived from this dataset using different methods: one based on travel times and another on sweep interferometry. Our first step is to refine these models using adjoint tomography, which has been applied successfully at global and regional-scales but is less common at the reservoir-scale. Two approaches are then explored for time-lapse analysis: directly comparing refined tomographic models from different stages or backpropagating waveform differences relative to a baseline tomographic model. The main take away is that both approaches highlight similar reservoir behaviors, but the latter approach is more computationally effective in capturing small-scale changes in subsurface properties. For this work, we leverage the use of Salvus (www.mondaic.com), an end-to-end seismic imaging solution, relying on the spectral element method to compute forward and adjoint simulations, and developed by Mondaic Ltd. It includes integrated workflow management that handles waveform and metadata, launches simulations, computes waveform misfits and adjoint sources, and iterates for model updates by nonlinear optimization.

15 GEOTHERMAL ENERGY↗

Subsurface Spectroscopy of Thermal Degradation Inside an Inert Plastic Bonded Explosive (PBX) Simulant Using Feedback-Assisted Wavefront Shaping

We characterize the subsurface thermal degradation of an inert analog of high-explosive molecular crystals (Eu:Y(acac) 3 (DPEPO)) (EYAD) embedded inside of a plastic bonded explosive simulant using feedback-assisted wavefront shaping-based fluorescence and Raman spectroscopies. This technique utilizes wavefront shaping to focus pump light inside a heterogeneous material onto a target particle, which significantly improves its spectroscopic signature. We find that embedding the EYAD crystals in the heterogeneous polymer results in improved thermal stability, relative to bare crystal measurements, with the crystal remaining fluorescent to >612 K inside of the heterogeneous material, while the bare crystal’s fluorescence is fully quenched by 500 K. We hypothesize that this improvement is due to the polymer restricting the effects of EYAD melting, which occurs at 400 K and is the primary mechanism for spectroscopic changes in the temperature range explored.

Anderson, Benjamin R.↗

Can MARSIS Measure the Low-Altitude Components of the Mars Magnetic Field?

Measuring the magnetic field anomaly of Mars at low altitudes (e.g. 100-200 km) can be an interesting application of Mars Advance Radar for Subsurface and Ionospheric Sounder (MARSIS). Due to a low HF operation frequency, the radio wave propagating in the ionosphere of Mars, over the magnetic anomaly regions, will be affected and distorted by the localized magnetic field. This distortion in the sounder signal is due to the Faraday rotation and provides information about the strength of the magnetic field. MARSIS is especially sensitive to the radial magnetic field at altitudes where the electron density in the ionosphere peaks (i.e. 100-200 km). Consequently, MARSIS is potentially capable of providing measurements for the radial component of the magnetic field at altitudes between 100 to 200 km that are normally out of reach for orbital magnetometers (with the exception of the aero-braking phase). Such low-altitude measurements would be complementary to already existing measurements at 400 km by MAG-ER on Mars Global Surveyor. This paper will explain the sensitivity of MARSIS as a magnetometer and the method envisioned to measure the radial magnetic field component. MARSIS (Picardi et al.), the first major planetary radar sounder, is the result of an international collaboration between NASA, the Italian Space Agency (ASI), and European Space Agency (ESA), and will arrive at Mars in early 2004 for a two-year mission. MARSIS has a frequency range between 0.1-5.5 MHz and is designed to penetrate the subsurface to a depth of a few kilometers. MARSIS primary objective is to map and characterize the subsurface geological structure of Mars, and search for subsurface liquid water reservoirs. The secondary objective of MARSIS is to study the ionosphere of Mars providing the most extensive amount of data on Martian ionosphere to date. In addition to MARSIS, a second radar sounder named SHARAD (SHallow RADar) with operation frequency of 15-25 MHz is under development. SHARAD is an Italian instrument (Seu et. al) that will fly on NASA s Mars Reconnaissance orbiter in 2005. SHARAD can also provide magnetic measurements, however, it is not expected to be as sensitive as MARSIS to magnetic field variations.

A Safaeinili↗