Subsurface Geophysical Imaging
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
SUMMARY Active or man-made seismic sources are often used to delineate subsurface geological structures via seismic imaging techniques. Although conventional controlled seismic sources have provided high-resolution subsurface images, the high cost of data acquisition necessitates improved use of cheaper alternative seismic sources for subsurface imaging, such as traffic noise, quarry blasts and earthquakes. Here, we explore the potential use of quarry blasts as an active seismic source for subsurface imaging. The goal of our study is to demonstrate the utility of quarry blasts as seismic sources for investigating subsurface structure with a field experiment in Mill Creek, southern Oklahoma. We deployed a 7455-m-long receiver line of 72 receivers (105 m spacing) extending away from a quarry, with a recording time of 28 d. We analyse similarities between blasts, wave types, origin time and excitation location of quarry blasts, which is important information for using such signal as an active source. Given a relative homogeneous geological condition in the study area, we could not observe clear reflected waves. However, the use of P- and S-wave velocities led to the delineation of a major basement damage zone in the study area, which further demonstrates the viability of quarry blasts as a useful seismic source for subsurface geological imaging.
The 3013-container-package consists of a convenience, inner and outer container and is used for long-term storage of plutonium-bearing materials. A destructive examination (DE) protocol has been developed to examine the container package visually and with microscopic aid to find any corrosion conditions that could result in the loss of the integrity of the container package over its lifetime. The DE protocol contains three main steps: initial container examination, helium (He) leak testing, and detailed imaging and analysis of the inner container closure weld region (ICCWR). The ICCWR has been determined to be bounding, defined as exhibiting worst case conditions for stress corrosion cracking (SCC) of the inner container. To assess SCC in the ICCWR, the inner container lid is cut into quarters and the weld and He-leak testing gasket are removed. Then a citric acid wash is performed to remove adherent chlorides from the ICCWR. The wash is then sent for analysis to determine the concentration of chlorides in the ICCWR. While the analysis for chloride concentration is being performed, the quarter sections are further sectioned into 1/8th subsections by cutting each section in half. These subsections are washed using nitric acid to remove corrosion products. Then each subsection is imaged using a Wide Angle 3-D Measurement System (WAMS). After analysis of microscope images for potential SCC, additional imaging can be performed, including subsurface imaging. After review, a determination is made of whether the container integrity may potentially affect the safe storage of the material.
Muon tomography has been used to seek hidden chambers in Egyptian pyramids and image subsurface features in volcanoes. It seemed likely that it could be used to image injected, supercritical carbon dioxide as it is emplaced in porous geological structures being used for carbon sequestration, and also to check on subsequent leakage. It should work equally well in any other application where there are two fluids of different densities, such as water and oil, or carbon dioxide and heavy oil in oil reservoirs. Continuous monitoring of movement of oil and/or flood fluid during enhanced oil recovery activities for managing injection is important for economic reasons. Checking on leakage for geological carbon storage is essential both for safety and for economic purposes. Current technology (for example, repeat 3D seismic surveys) is expensive and episodic. Muons are generated by high- energy cosmic rays resulting from supernova explosions, and interact with gas molecules in the atmosphere. This innovation has produced a theoretical model of muon attenuation in the thickness of rock above and within a typical sandstone reservoir at a depth of between 1.00 and 1.25 km. Because this first simulation was focused on carbon sequestration, the innovators chose depths sufficient for the pressure there to ensure that the carbon dioxide would be supercritical. This innovation demonstrates for the first time the feasibility of using the natural cosmic-ray muon flux to generate continuous tomographic images of carbon dioxide in a storage site. The muon flux is attenuated to an extent dependent on, amongst other things, the density of the materials through which it passes. The density of supercritical carbon dioxide is only three quarters that of the brine in the reservoir that it displaces. The first realistic simulations indicate that changes as small as 0.4% in the storage site bulk density could be detected (equivalent to 7% of the porosity, in this specific case). The initial muon flux is effectively constant at the surface of the Earth. Sensitivity of the method would be decreased with increasing depth. However, sensitivity can be improved by emplacing a greater array of particle detectors at the base of the reservoir.
This research program aims to address critical challenges in subsurface exploration and monitoring of anthropogenic CO 2 storage by advancing quantitative dynamic sensing and imaging technologies. Two primary applications—petrophysical assessment of hydrocarbon reservoirs (C1) and CO 2 storage monitoring (C2)—require significant innovation to overcome three barriers (B1-B3). These barriers include the need for enhanced understanding of interactions between physical wave fields (thermoacoustic, electromagnetic, acoustic/seismic, and X-ray) with fluid-filled porous media, the development of multi-sensor data fusion for real-time imaging, and upscaling microscopic quantum effects for macroscopic observations. The project’s objective was to establish a unified 4D coded sensing and imaging approach, integrating EM, AC/S, and TA fields for multi-scale and multi-physics material characterization and subsurface imaging.
VelocityGAN is a software to solve subsurface imaging problems. It leverages a generative adversarial network to model the correspondence from seismic data to subsurface velocity structures.
InversionNet is a software to solve subsurface imaging problems. It leverages a convolutional neural network with an encoder-decoder structure to model the correspondence from seismic data to subsurface velocity structures.
This IR is a follow-up conference Paper to the IR approved abstract; PNNL-SA-194096. --- Verification that equipment is authentic and not changed; even at the circuit card assembly (CCA) level will likely be an important component of future arms control treaties. This effort was an initial evaluation of the potential for laser-based ultrasound (LBU) as an inspection tool for Unique Identification (UID) of Circuit Card Assembly (CCA) boards and CCA components. The LBU system used a laser pulse for ultrasound generation and an Optical Microphone for ultrasound detection to image subsurface structures without physical contact and in a dry state. This paper described the selected CCA surrogate, CCA components that were examined, the LBU system, LBU images of the CCA components, and the initial development of UID algorithms. Receiver operating characteristics (ROC) curves indicated good performance for two algorithms as a potential means for UID of CCA boards and CCA components.
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.
The advent and subsequent proliferation of autonomous airborne, waterborne, and groundbased vehicles (i.e., “drones”) promises to broadly transform the geosciences and associated industries, including fossil energy exploration and development, mineral resource exploration and development, water-resource management, and environmental remediation. For geophysical characterization and monitoring, the prospect of programming highly repeatable and low-cost drone missions for subsurface imaging will allow for deployments in hazardous and previously inaccessible areas. Coupled with autonomous workflows for data processing, management, and visualization, drone-based geophysical characterization and monitoring will enable unprecedented, real-time insight into diverse subsurface properties and processes of scientific and engineering importance. Toward this end, the objectives of this Lab Directed Research and Development (LDRD) project were to develop new (1) instrumentation for dronebased electromagnetic induction (EMI) geophysical imaging, including separated transmitter and receivers and associated electronics, (2) software for real-time data telemetry, processing, management, and visualization. Although EMI has been previously deployed using unoccupied aerial systems (UASs), these applications failed to capitalize on the game-changing capabilities of drone platforms. Whereas drone-based data acquisition allows for collection of rich, three-dimensional (3D) multi-offset/multi-angle configurations between transmitters and receivers, past efforts have relied on conventional instrumentation that was designed for ground-based data collection with the transmitter and a single receiver housed in the same unit; nor did these previous applications demonstrate real-time delivery of results to support rapid management decisions in the field. In this 1-year project, we (1) designed and constructed new lightweight independent transmitter and receiver antenna platforms that communicate with a laptop computer; (2) developed software to control data acquisition, manage/transfer data, and visualize data as its collected; and (3) demonstrated the operation of the new hardware and software systems in a ground-based field test. Our work entails major technological advances for EMI and established a foundation on which to build a new drone-based, real-time geophysical EMI imaging capability to support diverse challenges facing the nation.
The SESAR-LITE (Space Exploration Synthetic Aperture Radar - Lunar Investigations Targeted Experiment) instrument is a compact P-band (70 cm wavelength) polarimetric synthetic aperture radar under development at the NASA Goddard Space Flight Center to measure the surface and subsurface of the Moon at full polarimetry and at meter-scale resolution. The radar will use a compact deployable antenna, distributed RF electronics, and multi-channel digital processing system to enable a set of focused mission goals for small payload opportunities. The instrument development is leveraging proven technology advancements recently developed and demonstrated at NASA Goddard Space Flight center for SESAR (Space Exploration Synthetic Aperture Radar), a flagship version of the instrument that was tailored for larger orbital missions. The development of SESAR-LITE addresses accommodation flexibility on multiple launch vehicle families that require small packages while providing unprecedented surface and subsurface imaging of the Moon called for NASA’s Artemis program.
The SESAR-LITE (Space Exploration Synthetic Aperture Radar - Lunar Investigations Targeted Experiment) instrument is a compact P-band (70 cm wavelength) polarimetric synthetic aperture radar under development at the NASA Goddard Space Flight Center (GSFC) to measure the surface and upper subsurface of the Moon at full polarimetry and at meter-scale resolution. The radar will use a compact deployable antenna, distributed RF electronics, and multi-channel digital processing system to enable a set of focused mission goals for small payload opportunities. The instrument development leverages proven technology advancements recently developed and demonstrated at NASA GSFC for SESAR (Space Exploration Synthetic Aperture Radar), a flagship version of the instrument that was tailored for larger orbital missions. The development of SESAR-LITE addresses accommodation flexibility on multiple launch vehicle families that require small packages while providing unprecedented surface and subsurface imaging of the Moon as required by NASA’s Artemis program.
Abstract Subsurface heterogeneity influences watershed hydrology strongly but remains difficult to characterize at catchment scales with sparse and costly field data. Geophysical surveys such as electromagnetic induction (EMI) provide local spatial subsurface images yet scaling them to watershed scales and converting EMI‐derived resistivity into hydraulic properties remains a challenge. We present a Model–Experiment (ModEx) framework that integrates limited EMI data with machine learning (ML) and hydrologic modeling to improve process representation and guide field investigations. Sparse EMI surveys were scaled to the catchment scale using a Random Forest model, and the resulting resistivity fields were combined with nearby borehole constraints to parameterize a hydrologic model. The EMI‐informed hydrological simulations improved predictions of streamflow sustained by subsurface flow and shallow saturation patterns. By combining EMI data and ML with hydrologic modeling, the ModEx framework guides future subsurface surveys, providing a transferable and efficient strategy for data–model integration across diverse watersheds. Plain Language Summary Mapping the underground network of soil and rock that controls water is essential for predicting floods and droughts, but seeing underground is difficult and expensive. We cannot drill everywhere, so scientists use geophysical tools to scan broad areas. There are two key challenges: these geophysical scans are often sparse across the whole watershed, and the geophysical data is hard to translate into water‐related properties. We used artificial intelligence to solve these problems. We taught a computer to find patterns linking the limited geophysical data to the land surface properties. This allowed it to fill in the gaps and create a complete, useful subsurface map for the entire watershed. This new map improves hydrologic simulations, leading to more accurate predictions of water movement in the watershed. It also helps scientists build better models with less data and generates a priority map showing where to measure next, making future investigations more efficient. Key Points Limited EMI scaled with ML improves catchment‐scale subsurface parameterization for hydrologic models The framework integrates hydrologic modeling with limited geophysical data to support subsurface investigation design ModEx framework offers a transferable data–model integration strategy that quantifies and reduces uncertainty guiding watershed studies
Materials whose luminescence can be switched by optical stimulation drive technologies ranging from superresolution imaging, nanophotonics, and optical data storage, to targeted pharmacology, optogenetics, and chemical reactivity. Furthermore, these photoswitchable probes, including organic fluorophores and proteins, can be prone to photodegradation and often operate in the ultraviolet or visible spectral regions. Colloidal inorganic nanoparticles can offer improved stability, but the ability to switch emission bidirectionally, particularly with near-infrared (NIR) light, has not, to our knowledge, been reported in such systems. Here, we present two-way, NIR photoswitching of avalanching nanoparticles (ANPs), showing full optical control of upconverted emission using phototriggers in the NIR-I and NIR-II spectral regions useful for subsurface imaging. Employing single-step photodarkening and photobrightening, we demonstrate indefinite photoswitching of individual nanoparticles (more than 1,000 cycles over 7 h) in ambient or aqueous conditions without measurable photodegradation. Critical steps of the photoswitching mechanism are elucidated by modelling and by measuring the photon avalanche properties of single ANPs in both bright and dark states. Unlimited, reversible photoswitching of ANPs enables indefinitely rewritable two-dimensional and three-dimensional multilevel optical patterning of ANPs, as well as optical nanoscopy with sub-Å localization superresolution that allows us to distinguish individual ANPs within tightly packed clusters.
Part of the DEEPEN (DE-risking Exploration of geothermal Plays in magmatic ENvironments) project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology, detailed in the journal article below, is based on the method proposed by Poux & O'brien (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., (2022). For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty.
DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. Part of the DEEPEN project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology is based on the method proposed by Poux et al. (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., 2022. For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. More information on this process is described in the following sections. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty. There are two different versions of the Leapfrog model and associated favorability models: - v1.0: The first release in June 2023 - v2.1: The second release, with improvements made to the earthquake catalog (included additional identified events, removed duplicate events), to the temperature model (fixed a deep BHT), and to the index models (updated the seismicity-heat source index models for supercritical and EGS, and the resistivity-insulation index models for all three play types). Also uses the jet color map rather than the magma color map for improved interpretability. - v2.1.1: Updated to include v2.0 uncertainty results (see below for uncertainty model versions) There are two different versions of the associated uncertainty models: - v1.0: The first release in June 2023 - v2.0: The se...
This is the final technical report for the Virginia Tech Center for Photonics Technology (VT-CPT) research project entitled “Fully Distributed Acoustic and Magnetic Field Monitoring Via a Single Fiber Line for Optimized Production of Unconventional Resource Plays”. In the the four-year effort, VT-CPT collaborated with the Departments of Mathmatics at VT and Sentek Instrument, LLC (Sentek) to develop a fiber-optic sensing system capable of real-time simultaneous distributed measurement of multiple subsurface, drilling, and production parameters. The ultra-sensitive fiber optic distributed acoustic sensing technology, picoDAS, developed by Sentek was integrated with a novel multi-material optical sensing fiber fabricated by CPT-VT to obtain distributed acoustic and magnetic field measurements with ultrahigh sensitivity and high spatial resolution. The sensing technology successfully developed and demonstrated under this research program is truly unique has application beyond subsurface imaging to include carbon storage and electrical grid monitoring, as well as healthcare and nuclear fusion. The project, sponsored by the Advanced Technology Program, began in 2019 and had the original goal of developing the next generation of harsh environment sensing systems. Exhaustive theoretical modeling was performed to optimize the multi-material optical fiber design. VT developed the processing techniques employed to successfully fabricate single mode fibers with metal cladding wires. Sensing fibers with nickel and Metglas cladding wires were fully characterized, inscribed with sensors, and fully integrated with the picoDAS system. Laboratory scale testing was performed to demonstrate performance upon exposure to lateral and transverse magnetic fields. Preliminary field trials were performed for prototype magnetic and acoustic sensing systems upon near surface deployment at VT. The distributed acoustic and magnetic field sensing system was advanced from a (Technology Readiness Level) TRL=2 to a TRL=5 via laboratory scale system validation in relevant environments. The truly one-of-its kind distributed magnetic and acoustic sensing system is expected to find immediate applications that require subsurface monitoring, to include carbon storage site monitoring, and has the potential to be a distributive technology in healthcare and electric grid markets. The technologies developed by Virginia Tech’s Center for Photonics Technology will support the mission of the National Energy Technology Laboratory (NETL) to drive innovation and deliver solutions for an environmentally sustainable and prosperous energy future by ensuring affordable, abundant, and reliable energy that drives a robust economy and national security. Technical accomplishments during the program are summarized briefly here and described in detail in the remainder of the report.
Here, we present deep learning (DL) networks for three-dimensional (3D) joint inversion of active seismic full waveform and passive seismic traveltime data to image reservoirs and their properties and quantify imaging uncertainties. Active seismic full-waveform data can provide high-resolution monitoring images but are collected only intermittently because of their high acquisition cost. In contrast, passive seismic data can be gathered at relatively low cost between regular active surveys, although their imaging quality can be compromised by factors such as low signal-to-noise ratios and limited ray coverage of the target. Although these datasets are routinely acquired together at CO 2 storage sites, their combined inversion within a 3D DL framework has not been previously demonstrated. To our knowledge, this is the first study to address this gap, combining the strength of both data types. For efficient data storage and DL training with large 3D seismic datasets, we use a 3D data matrix in which a random number of passive seismic traveltime data are stored as parabolic envelopes using one-hot encoding and a 3D full-waveform data matrix in which multiple shot gathers are summed. Two network architectures are evaluated: a single-encoder U-Net for single-data type inversion and a dual-encoder U-Net for joint inversion of active and passive seismic data. We also evaluate the single-encoder U-Net for joint inversion by concatenating full-waveform data and traveltime data. We propose a systematic approach for selecting an optimal dropout rate that balances regularization during training and Monte Carlo dropout-based uncertainty quantification during prediction by examining the correlation coefficient between standard deviation and prediction error, along with the training misfit, across a range of dropout rates. 3D DL inversion experiments include five different network configurations, with evaluations under ideal, noisy and dropout-enabled conditions. Both model and data uncertainties are assessed, as well as their combined effects. Across all conditions, the networks consistently predict accurate CO 2 saturation models with low prediction errors, such as a structural similarity index of 0.993 and CO 2 difference of 1.1%. Uncertainty estimates show strong spatial correlation with prediction errors, confirming the effectiveness of the proposed dropout selection approach. The results demonstrate that our DL approach, utilizing compact data representations and appropriate uncertainty quantification, yields accurate subsurface images under various inversion conditions and provides valuable insights into the reliability of predictions.