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 667 records · Page 37

Detecting the permafrost carbon feedback: talik formation and increased cold-season respiration as precursors to sink-to-source transitions

Abstract. Thaw and release of permafrost carbon (C) due to climate change is likely to offset increased vegetation C uptake in northern high-latitude (NHL) terrestrial ecosystems. Models project that this permafrost C feedback may act as a slow leak, in which case detection and attribution of the feedback may be difficult. The formation of talik, a subsurface layer of perennially thawed soil, can accelerate permafrost degradation and soil respiration, ultimately shifting the C balance of permafrost-affected ecosystems from long-term C sinks to long-term C sources. It is imperative to understand and characterize mechanistic links between talik, permafrost thaw, and respiration of deep soil C to detect and quantify the permafrost C feedback. Here, we use the Community Land Model (CLM) version 4.5, a permafrost and biogeochemistry model, in comparison to long-term deep borehole data along North American and Siberian transects, to investigate thaw-driven C sources in NHL (> 55∘ N) from 2000 to 2300. Widespread talik at depth is projected across most of the NHL permafrost region (14 million km2) by 2300, 6.2 million km2 of which is projected to become a long-term C source, emitting 10 Pg C by 2100, 50 Pg C by 2200, and 120 Pg C by 2300, with few signs of slowing. Roughly half of the projected C source region is in predominantly warm sub-Arctic permafrost following talik onset. This region emits only 20 Pg C by 2300, but the CLM4.5 estimate may be biased low by not accounting for deep C in yedoma. Accelerated decomposition of deep soil C following talik onset shifts the ecosystem C balance away from surface dominant processes (photosynthesis and litter respiration), but sink-to-source transition dates are delayed by 20–200 years by high ecosystem productivity, such that talik peaks early (∼ 2050s, although borehole data suggest sooner) and C source transition peaks late (∼ 2150–2200). The remaining C source region in cold northern Arctic permafrost, which shifts to a net source early (late 21st century), emits 5 times more C (95 Pg C) by 2300, and prior to talik formation due to the high decomposition rates of shallow, young C in organic-rich soils coupled with low productivity. Our results provide important clues signaling imminent talik onset and C source transition, including (1) late cold-season (January–February) soil warming at depth (∼ 2 m), (2) increasing cold-season emissions (November–April), and (3) enhanced respiration of deep, old C in warm permafrost and young, shallow C in organic-rich cold permafrost soils. Our results suggest a mosaic of processes that govern carbon source-to-sink transitions at high latitudes and emphasize the urgency of monitoring soil thermal profiles, organic C age and content, cold-season CO2 emissions, and atmospheric 14CO2 as key indicators of the permafrost C feedback.

54 ENVIRONMENTAL SCIENCES↗

3D Characterization and Time-Lapse Imaging of the Desiccation Treatability Test at the Hanford BC-Cribs and Trenches Site using High Performance Electrical Resistivity Imaging

Contaminated vadose zone materials are a potential source of long-term groundwater contamination at many sites across the Department of Energy (DOE) complex. Deep vadose zone contamination presents a particularly challenging remedial problem due to the difficultly of locating contaminants and the expense of access and ex situ treatment. In situ remediation techniques, whereby remedial amendments must be delivered to contaminated soils, have been identified as a potential alternative. However, amendment delivery is typically uncertain and post delivery remedial performance is often not well understood due to limited information available from sparsely spaced boreholes. Recent advancements in electrical resistivity tomography (ERT) are being used to address these challenges at the Hanford site by providing remote, three-dimensional images of contaminant distribution and four-dimensional (three spatial dimensions plus the time dimension) images of in-situ vadose zone remediation processes. These capabilities were recently demonstrated with a large scale surface characterization effort and during a smaller scale desiccation treatability test at the Hanford BC-Cribs area. The results of these efforts demonstrate the utility of leveraging high-performance computing resources to process ERT data for 3D reconnaissance mapping of subsurface contaminants and for detailed 3D monitoring of vadose zone remediation efforts.

contamination↗

Real time explosive hazard information sensing, processing, and communication for autonomous operation

Methods, computer readable media, and apparatuses provide robotic explosive hazard detection. A robot intelligence kernel (RIK) includes a dynamic autonomy structure with two or more autonomy levels between operator intervention and robot initiative A mine sensor and processing module (ESPM) operating separately from the RIK perceives environmental variables indicative of a mine using subsurface perceptors. The ESPM processes mine information to determine a likelihood of a presence of a mine. A robot can autonomously modify behavior responsive to an indication of a detected mine. The behavior is modified between detection of mines, detailed scanning and characterization of the mine, developing mine indication parameters, and resuming detection. Real time messages are passed between the RIK and the ESPM. A combination of ESPM bound messages and RIK bound messages cause the robot platform to switch between modes including a calibration mode, the mine detection mode, and the mine characterization mode.

97 MATHEMATICS AND COMPUTING↗

Real time explosive hazard information sensing, processing, and communication for autonomous operation

Methods, computer readable media, and apparatuses provide robotic explosive hazard detection. A robot intelligence kernel (RIK) includes a dynamic autonomy structure with two or more autonomy levels between operator intervention and robot initiative A mine sensor and processing module (ESPM) operating separately from the RIK perceives environmental variables indicative of a mine using subsurface perceptors. The ESPM processes mine information to determine a likelihood of a presence of a mine. A robot can autonomously modify behavior responsive to an indication of a detected mine. The behavior is modified between detection of mines, detailed scanning and characterization of the mine, developing mine indication parameters, and resuming detection. Real time messages are passed between the RIK and the ESPM. A combination of ESPM bound messages and RIK bound messages cause the robot platform to switch between modes including a calibration mode, the mine detection mode, and the mine characterization mode.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

A hydrogeophysical framework to assess infiltration during a simulated ecosystem-scale flooding experiment

This study presents a framework to quantify changes in soil saturation in response to flooding caused by extreme hydrologic perturbation on coastal ecosystems at the interfaces and transition between terrestrial and aquatic systems. Subsurface heterogeneity limits the use of in situ measurements to quantify subsurface flow during flooding due to the spatial discontinuity in the measured data. While geophysical methods, including time-lapse electrical resistivity imaging (ERI), are increasingly used to monitor soil hydrological processes, their abilities to parameterize flow models have been underutilized. This study combines background ERI, ground penetrating radar (GPR), time-lapse ERI, soil characterization, and a numerical flow model developed using an Advanced Terrestrial Simulator (ATS) code to quantify the infiltration pathway and describe the hydrological dynamics during a simulated flooding experiment. We assessed the use of two conceptual models developed using [1] ERI and GPR data that described the stratigraphic distribution, and time-lapse ERI that mapped permeability contrast, and [2] information from a national soil database for capturing changes in saturation. Combining the ERI and GPR results with soil core data revealed the stratigraphic heterogeneity at the site with a silty clay layer from 1 to 2 m between an overlying loamy topsoil and an underlying saturated silty sand. This silty clay layer could restrict deep infiltration. The time-lapse ERI showed up to a 35% decrease in resistivity, which correlated with soil moisture data (R 2 value > 0.53) and revealed preferential infiltration zones used to inform the flow model. Numerical simulation results from both the geophysics- and soil database-informed models quantified changes in soil saturation with calculated soil moistures that agreed with field data. The geophysics-informed model captured more of the system’s variability, reflective of shallow subsurface heterogeneities. The framework presented will serve as a precursor for a robust ecohydrological model that can describe the impacts of extreme events induced by climate change on coastal ecosystems.

54 ENVIRONMENTAL SCIENCES↗

Tc(VII) and Cr(VI) Interaction with Naturally Reduced Ferruginous Smectite from a Redox Transition Zone

Fe(II)-rich clay minerals found in subsurface redox transition zones (RTZs) can serve as important source of electron equivalents limiting the transport of redox active contaminants. While most laboratory reactivity studies are based on reduced model clays, the reactivity of naturally reduced clays in field samples remains poorly explored. Characterization of the clay size fraction of a fine-grained unit from RTZ interface at the Hanford site, Washington, including mineralogy, crystal chemistry, and Fe(II)/(III) content, indicates that ferruginous montmorillonite is the dominant mineralogical component. Oxic and anoxic fractions differ significantly in Fe(II) concentration, but FeTOTAL remains constant demonstrating no Fe loss during reduction-oxidation cycling. At its native pH of 8.6, the anoxic fraction despite its significant Fe(II) (~23% of FeTOTAL), exhibits minimal reactivity with TcO4- and CrO42- and much slower reaction kinetics than that measured in studies with biologically/chemically reduced model clays. Reduction capacity is enhanced by added Fe(II) (if Fe(II)SORBED >8% clay Fe(II)LABILE), however the kinetics of this conceptually surface-mediated reaction remain sluggish. Surface-sensitive Fe L-edge X-ray absorption spectroscopy shows that Fe(II)SORBED and the resulting reducing equivalents are not available in the outermost few nanometers of clay surfaces. Slow kinetics thus appear related to diffusion-limited access to electron equivalents retained within clay mineral.

Environmental Molecular Sciences Laboratory↗

Role of Surface Termination in the Metal–Insulator Transition of V 2 O 3 (0001) Ultrathin Films

Surface termination is known to play an important role in determining the physical properties of materials. It is crucial to know how surface termination affects the metal–insulator transition (MIT) of V 2 O 3 films for both fundamental understanding and its applications. By changing growth parameters, we achieved a variety of surface terminations in V 2 O 3 films that are characterized by low-energy electron diffraction (LEED) and photoemission spectroscopy techniques. Depending upon the terminations, our results show that MIT can be partially or fully suppressed near the surface region due to the different fillings of the electrons at the surface and subsurface layers and the change of screening length compared to the bulk. Across MIT, a strong redistribution of spectral weight and its transfer from a high-to-low-binding energy regime is observed in a wide energy scale. Our results show that the total spectral weight in the low-energy regime is not conserved across MIT, indicating a breakdown of the “sum rules of spectral weight”, signature of a strongly correlated system. Such a change in spectral weight is possibly linked to the change in hybridization, lattice volume (i.e., effective carrier density), and the spin degree of freedom in the system that occurs across MIT. We find that MIT in this system is strongly correlation-driven, where the electron–electron interactions play a pivotal role. Furthermore, our results provide better insight into the understanding of the electronic structure of strongly correlated systems and highlight the importance of accounting for surface effects during interpretation of the physical property data mainly using surface-sensitive probes, such as surface resistivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Fluorescent Functionalized Mesoporous Silica for Radioactive Material Extraction

Mesoporous silica with covalently bound salicylic acid molecules incorporated in the structure was synthesized with a one-pot, co-condensation reaction at room temperature. The as-synthesized material has a large surface area, uniform particle size, and an ordered pore structure as determined by characterization with transmission electron microscopy, thermal gravimetric analysis, and infrared spectra, etc. Using the strong fluorescence and metal coordination capability of salicylic acid, functionalized mesoporous silica (FMS) was developed to track and extract radionuclide contaminants, such as uranyl [U(VI)] ions encountered in subsurface environments. Adsorption measurements showed a strong affinity of the FMS toward U(VI) with a Kd value of 105 mL/g, which is four orders of magnitude higher than the adsorption of U(VI) onto most of the sediments in natural environments. The new materials have a potential for synergistic environmental monitoring and remediation of the radionuclide U(VI) from contaminated subsurface environments.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Ambient‐Pressure Relithiation of Degraded Li x Ni 0.5 Co 0.2 Mn 0.3 O 2 (0 < x < 1) via Eutectic Solutions for Direct Regeneration of Lithium‐Ion Battery Cathodes

Abstract With the rapid growth of the lithium‐ion battery (LIBs) market, recycling and re‐use of end‐of‐life LIBs to reclaim lithium (Li) and transition metal (TM) resources (e.g., Co, Ni), as well as eliminating pollution from disposal of waste batteries, has become an urgent task. Here, for the first time the ambient‐pressure relithiation of degraded LiNi 0.5 Co 0.2 Mn 0.3 O 2 (NCM523) cathodes via eutectic Li + molten‐salt solutions is successfully demonstrated. Combining such a low‐temperature relithiation process with a well‐designed thermal annealing step, NCM523 cathode particles with significant Li loss (≈40%) and capacity degradation (≈50%) can be successfully regenerated to achieve their original composition and crystal structures, leading to effective recovery of their capacity, cycling stability, and rate capability to the levels of the pristine materials. Advanced characterization tools including atomic resolution electron microscopy imaging and electron energy loss spectroscopy are combined to demonstrate that NCM523's original layered crystal structure is recovered. For the first time, it is shown that layer‐to‐rock salt phase change on the surfaces and subsurfaces of the cathode materials can be reversed if lithium can be incorporated back to the material. The result suggests the great promise of using eutectic Li + molten–salt solutions for ambient‐pressure relithiation to recycle and remanufacture degraded LIB cathode materials.

Shi, Yang↗

Characterization of Most Promising Sequestration Formations in the Rocky Mountain Region

The project Characterization of Most Promising Sequestration Formations in the Rocky Mountain Region is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The Rocky Mountain Carbon Capture and Storage (RMCCS) project investigated multiple geologic formations and characterized a local site on the Colorado Plateau for future CCS opportunities. The RMCCS project focused on the Cretaceous Dakota, Jurassic Entrada, and Pennsylvanian Weber Sandstones, the three largest regional formations. All formations in this project are potential CO2 storage resources for future power plants, natural gas processing plants, cement plants, and oil shale development projects. The area adjacent to Craig, Colorado, (Sand Wash Basin) was the area selected for detailed geologic characterization on the RMCCS project. The basin was selected in part because the geology can be extrapolated to other sites on the Colorado Plateau. Field mapping and seismic surveys were conducted to identify and evaluate the basin's structural configuration. A 9,745-foot deep characterization well was drilled to collect 131 feet of core and a suite of geophysical well log data. Petrophysical tests on samples of core were used to calibrate geophysical log data, which can be used to obtain storage resource estimates and evaluate associated uncertainty as well as simulate the hydrologic behavior of injected CO2. A detailed analysis of the primary formations (Dakota, Entrada and Weber sandstones) yielded a more accurate CO2 storage resource assessment for these formations within the Colorado Plateau; RMCCS estimates indicate a total CO2 storage resource of more than 38,000 million metric tons. The characterization of the Sand Wash Basin (2-D seismic surveys, multiple well logs and lithological, petrophysical and geochemical analyses) allowed for a detailed 3-D model to be constructed. The model served as the framework for analyses ranging from CO2 storage resource, injectivity, and subsurface flow to uncertainty estimates to evaluation of risk.

2-D seismic↗

Interpretation of Hydraulic Fractures Based on Microseismic Response in the Marcellus Shale, Monongalia County, West Virginia, USA: Implications for Shale Gas Production

Summary Hydraulic fracturing is critical for extracting shale gas in the subsurface. The treatment technique of multistage hydraulic fracturing is widely used to maximize production. In multistage hydraulic fracturing, not all pumping stages make the same contribution for production, although the designed stimulation process is almost same in every pumping stage. In this study, we characterize the microseismic responses of hydraulic fracturing. Significant variations of microseismic characteristics are observed among the different pumping stages. Pre-existing natural fractures are examined along the horizontal well in selected stages. Combined with the Mohr’s circle analysis, results of the fracture study show that induced hydraulic fractures can be captured by pre-existing natural fractures. Induced hydraulic fractures are simulated by the unconventional fracture model (UFM), and the result reveals that the stress-shadow effect diverts the direction of hydraulic fractures. The diverted hydraulic fracture affects the development of the hydraulic fracture network, which has influence on production. The treatment of two-step pumping is investigated by application of hydraulic diffusivity. Hydraulic fracturing performance and production could be optimized by the treatment of two-step pumping in a single stage. The second pumping creates new fractures and fills the fractures with additional proppants to maintain production for a long duration.

Energy & Fuels↗

Pore-scale visualization of natural hydrate-bearing sediments

Accurate modeling of gas hydrate reservoir productivity and geomechanical risks associated with subsurface dissociation of natural gas hydrates (NGH) requires the determination of model parameters through physical testing on natural hydrate-bearing sediments (HBS). This involves investigating the hydro-mechanical behavior of undisturbed hydrate samples from nature under in situ conditions using pressure core characterization and analysis, which provides a unique opportunity for research. By employing state-of-the-art micro computed tomography imagery on cryogenically preserved, hydrate-bearing sediment samples, we can determine hydrate saturation as well as permeability with and without the presence of hydrates in the sediment. Furthermore, utilizing a machine learning based image segmentation technique, it is possible to extract pore space and grain information. Subsections of the entire image volume were used to determine anisotropic permeabilities using a finite-difference method Stokes solver (FDMSS). Additionally, permeability measurements on whole pressure and temperature preserved hydrate-bearing core were analyzed by utilizing the National Energy Technology Laboratory’s (NETL) Pressure Core Characterization and X-ray CT Visualization Tool (PCXT) to manipulate, cut, and analyze pressure preserved sediment. Permeabilities were measured under a broad range of vertical stress states to simulate expected pressure changes during production scenarios, and the results show that permeabilities derived from images are in agreement with those from traditional core derived experiments. The collected stress-dependent permeability, permeability anisotropy, and corresponding gas hydrate saturations provide valuable input into numerical simulations of reservoir productivity. These properties have been proven to be key parameters determining a long-term reservoir response under depressurization.

Liu, Mengwei [Oak Ridge Institute for Science and ↗

Fracture Network Quantification for Enhanced Reservoir Characterization in CO2 Storage Sites

Presentation slides on “Fracture Network Quantification for Enhanced Reservoir Characterization in CO2 Storage Sites” for the CCUS 2025 conference held in Houston, Texas March 3-5, 2025. Accurate reservoir characterization is fundamental for the safe and efficient operation of carbon capture, utilization, and storage (CCUS) projects. This process is not only critical during initial site selection but also plays a key role in continuously updating reservoir geomodels throughout injection operations to effectively manage subsurface storage sites. Our primary objective is to develop a comprehensive understanding of fracture networks within CO2 injection reservoirs, which are critical for predicting the behavior of injected fluids, ensuring storage integrity post-injection, and mitigating risks such as leakage or induced seismicity.

artifical intelligence / machine learning (AI/ML)↗

Fully Distributed Acoustic and Magnetic Field Monitoring Via a Single Fiber Line for Optimized Production of Unconventional Resource Plays

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.

02 PETROLEUM↗

Nanoscale Editing of Multi and Single Layer Tungsten Disulfide via Gas‐Assisted Focused Electron Beam Induced Etching for Device Prototyping

Focused electron beam induced etching (FEBIE) with XeF 2 (xenon difluoride) precursor is conducted on multi-layer exfoliated WS 2 (tungsten disulfide) and monolayer WS 2 grown by chemical vapor deposition (CVD). The films are characterized by atomic force microscopy (AFM) and Raman and photoluminescence (PL) spectroscopy post-etching. The etch rates/efficiencies are reported as a function of electron beam energy, current, dwell time, and XeF 2 pressure. Bulk film Raman spectra are unchanged post-FEBIE, indicating minimal subsurface damage. Monolayer WS 2 shows a decrease in Raman and PL intensity post-FEBIE, with a dose-to-clear of ≈2 nC µm −2 . The study reveals regimes affected by the various mass transport contributions such as refresh time and the ratio of electrons/XeF 2 . Spontaneous etching was discovered during FEBIE of large patterned areas due to the long frame/refresh times. Density functional theory and ab initio molecular dynamics simulations compares desorption of SF x and WF x molecules from pristine WS 2 basal planes and pore edges, revealing the spontaneous etching is consistent with etching of partially etched monolayers during each frame. Single-line etching width of 21 nm, and patterning flakes into 100 nm wide channels are demonstrated. In conclusion, this work demonstrates the possibility of editing WS 2 flakes into electronic devices of arbitrary dimensions for semiconductor applications.

2D materials↗

Liquid-Rich Shale Potential of Utah’s Uinta and Paradox Basins: Reservoir Characterization and Development Optimization

The enclosed report is the culmination of a multi-year and multi-faceted research project investigating Utah’s unconventional tight oil potential. From the beginning, the project team focused efforts on two different plays: (1) the basal Green River Formation’s (GRF) Uteland Butte unconventional play in the Uinta Basin and (2) the more established but understudied Cane Creek shale play in the Paradox Basin. The 2009-2014 high price of crude oil, coupled with lower natural gas prices, generated renewed interest in exploration and development of liquid hydrocarbon reserves. Following the success of the mid-2000s shale gas boom and employing many of the same well completion techniques, petroleum companies started exploring for liquid petroleum in shale formations. In fact, many shales targeted for natural gas include areas in which the shale is more prone to liquid production. In Utah, organic-rich shales in the Uinta and Paradox Basins have been the source of significant hydrocarbon generation, with companies traditionally targeting the interbedded sands or carbonates for their conventional resource recovery. Because of the advances in horizontal drilling and hydraulic fracturing techniques, operators in these basins started to explore the petroleum production potential of the shale units themselves. The GRF in the Uinta Basin has been studied for over 50 years, since the first hydrocarbon discoveries. However, those studies focused on the many conventional sandstone reservoirs currently producing oil and gas. In contrast, less information was available about the more unconventional crude oil production potential of thinner carbonate/shale units, most notably the basal Uteland Butte member. The Cane Creek shale of the Paradox Basin has been a target for exploration periodically since the 1960s and produces oil from several small fields. The play generated much interest in the early 1990s with the successful use of horizontal drilling. Recently, the USGS assessed the undiscovered oil resource in the Cane Creek shale of the Paradox Basin at 103 million barrels at a 95 percent confidence level and 198 million barrels at a 50 percent confidence level. Nonetheless, limited research was available or published to further define the play and the reservoir characteristics. The specific objectives of the enclosed research were to (1) characterize geologic, geochemical, and geomechanical rock properties of target zones in the two designated basins by compiling data and by analyzing available cores, cuttings, and well logs; (2) describe outcrop reservoir analogs of GRF plays (Cane Creek shale is not exposed) and compare them to subsurface data; (3) map major regional trends for targeted intervals and identify “sweet spots” that have the greatest oil potential; (4) reduce exploration costs and drilling risks, especially in environmentally sensitive areas; (5) improve drilling and fracturing effectiveness by determining optimal well completion design; and (6) reduce field development costs, maximize oil recovery, and increase reserves. These objectives are all addressed in a series of nine publications that resulted from this extensive research project. Each publication is included in this report as an independent appendix.

02 PETROLEUM↗

Physics-embedded inverse analysis with algorithmic differentiation for the earth’s subsurface

Abstract Inverse analysis has been utilized to understand unknown underground geological properties by matching the observational data with simulators. To overcome the underconstrained nature of inverse problems and achieve good performance, an approach is presented with embedded physics and a technique known as algorithmic differentiation. We use a physics-embedded generative model, which takes statistically simple parameters as input and outputs subsurface properties (e.g., permeability or P-wave velocity), that embeds physical knowledge of the subsurface properties into inverse analysis and improves its performance. We tested the application of this approach on four geologic problems: two heterogeneous hydraulic conductivity fields, a hydraulic fracture network, and a seismic inversion for P-wave velocity. This physics-embedded inverse analysis approach consistently characterizes these geological problems accurately. Furthermore, the excellent performance in matching the observational data demonstrates the reliability of the proposed method. Moreover, the application of algorithmic differentiation makes this an easy and fast approach to inverse analysis when dealing with complicated geological structures.

54 ENVIRONMENTAL SCIENCES↗