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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Hawai‘i Play Fairway (Final Report)

Most of Hawai'i's geothermal resources are blind—their manifestations, such as hot springs and steam vents, do not appear on the ground surface because the heated water flows far below. With the exception of $K\bar{i}lauea East Rift Zone$, in most areas of Hawai'i, high lateral permeability in the first kilometer below ground surface prevents surface thermal features from developing. As a methodology for discovering these blind resources, Play Fairway Analysis (PFA) involves finding potential locations of blind hydrothermal systems and describing potential geothermal sources in rift-zone settings. Using the PFA to find Hawai'i's geothermal resources, the University of Hawai'i (UH) conducted the Hawai‘i Play Fairway Project, Hawai'i's first statewide geothermal resource assessment since 1985. Sponsored by the U.S. Department of Energy, the Hawai'i Play Fairway Project provided an updated resource assessment, a roadmap for additional exploration activities, and the identification of areas for further exploration. Benefitting from UH's core competency in earth sciences and experienced geothermal researchers, the project comprised three phases. During the first phase, the team identified, compiled, and ranked existing geologic, groundwater, and geophysical datasets relevant to subsurface heat, fluid and permeability. Using a Bayesian statistical approach, the team developed a statistical methodology to integrate these data into a resource probability map. The team evaluated the confidence in the probability value and considered development viability of areas with geothermal resources. With these analyses, the team identified 10 locations in the Hawaiian Islands for exploration activities. For the second phase, the team collected new groundwater data in 10 locations across the state and new geophysical data on $L\bar{a}na‘i, Maui$, and central Hawai'i Island and modeled topographically induced stress to better characterize subsurface permeability. Analyzing the subsurface stresses, the team evaluated the potential for fracture-induced permeability. The team inverted the MT and gravity data to produce 3D models of resistivity and density, respectively, on $L\bar{a}na‘i$, across $Haleakal\bar{a}'s$ SW rift (Maui), and surrounding Mauna Kea (Hawai‘i Island). The team developed and applied a new method for incorporating depth information about resistivity, density, and potential for fracture-induced permeability into the statistical method for computing resource probability in these three focus areas. The team incorporated the new groundwater results with the new geophysical results and the calculations of potential for fracture-induced permeability to produce updated maps of resource probability and confidence. Through combining data from the first and second phases, the team determined locations for further exploration during the third phase. For MT and gravity surveys, the team recommended $Kaua'i's$ $L\bar{i}hu'e$ $Basin$, the east rift of $Maui's$ $Haleakal\bar{a}$ volcano, and the southwest rift of Hawai'i Island's Mauna Loa volcano. The MT and gravity surveys aimed to enable improved confidence in the resource potential in these locations. For drilling deep groundwater well(s), the team recommended Southeast Mauna Kea and $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$. During the third phase, further exploration involved drilling a groundwater well in $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$ and performing more geophysical surveys. We deepened an existing water well proximal to our target area on $L\bar{a}na'i$ due to funding constraints that precluded us from spudding a new well that would exceed 1km depth. Drilling was preceded by a number of substantial elements including: writing an Environmental Assessment and the subsequent legal process, performance of deviation logging, lowering a camera down the well, coordinating site preparation with $P\bar{u}lama$ $L\bar{a}na'i$, shipping the UH-owned rig interisland, procuring supplies, and leading 3 community meetings on $L\bar{a}na'i$. Drilling occurred 24/7 the entire month of June 2019 over which time $L\bar{a}na'i$ $Well$ 10 was deepened from 427 m to 1057 m, with continuous core collected. We measured a roughly linear temperature gradient averaging 42°C/km and a maximum bottom hole temperature of 66°C. This gradient is more than twice the background for Hawai'i and within a range of gradients measured in this depth range for some exploration wells within KERZ. We consider these results encouraging for $L\bar{a}na'i's$ resource potential and recommend following with a slim hole within $L\bar{a}na'i's$ caldera (our target zone) to ~ 2 km. Further, the positive implications such results have for the island of O‘ahu are substantial - the shield stage of O'ahu's volcanoes ended 1-2 My earlier. However, O'ahu uses more electricity than the rest of the islands combined, and the utility recently called for 500-700MW of firm, dispatchable renewable electricity on O'ahu by 2033. In Phase 3, we also collected limited new encouraging groundwater data, and updated our thoughts on the probabilities of fluid and permeability at resource depths (Pr F = 1; Pr P = mostly unconstrained). Ultimately, we advocate for using our final probability of heat, and confidence in this probability, to drive the next phase of exploration. We contend further development of geothermal in Hawai‘i will enable the state to achieve its 100% renewable policy objective and Hawai'i to transition off of fossil fuels through geothermal discovery and development. The project not only produced a large amount of data and expanded the existing knowledge of Hawai'i's geothermal resources, but also produced publications, theses, presentations, core photos, datasets, media reports, television interviews, community events, and a blog. Students and new professionals benefitted from the project's hands-on research experiences and educational opportunities and earned awards and recognition.

15 GEOTHERMAL ENERGY↗

DOE OSTI2012

Three-Dimensional Bayesian Geostatistical Aquifer Characterization at the Hanford 300 Area using Tracer Test Data

Tracer testing under natural or forced gradient flow holds the potential to provide useful information for characterizing subsurface properties, through monitoring, modeling and interpretation of the tracer plume migration in an aquifer. Non-reactive tracer experiments were conducted at the Hanford 300 Area, along with constant-rate injection tests and electromagnetic borehole flowmeter (EBF) profiling. A Bayesian data assimilation technique, the method of anchored distributions (MAD) [Rubin et al., 2010], was applied to assimilate the experimental tracer test data with the other types of data and to infer the three-dimensional heterogeneous structure of the hydraulic conductivity in the saturated zone of the Hanford formation. In this study, the Bayesian prior information on the underlying random hydraulic conductivity field was obtained from previous field characterization efforts using the constant-rate injection tests and the EBF data. The posterior distribution of the conductivity field was obtained by further conditioning the field on the temporal moments of tracer breakthrough curves at various observation wells. MAD was implemented with the massively-parallel three-dimensional flow and transport code PFLOTRAN to cope with the highly transient flow boundary conditions at the site and to meet the computational demands of MAD. A synthetic study proved that the proposed method could effectively invert tracer test data to capture the essential spatial heterogeneity of the three-dimensional hydraulic conductivity field. Application of MAD to actual field data shows that the hydrogeological model, when conditioned on the tracer test data, can reproduce the tracer transport behavior better than the field characterized without the tracer test data. This study successfully demonstrates that MAD can sequentially assimilate multi-scale multi-type field data through a consistent Bayesian framework.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

A New Approach to Predict Hydrogeological Parameters Using Shear Waves from the Multichannel Analysis of Surface Waves Method

For near-surface contaminant characterization, the accurate prediction of hydrogeological parameters in anisotropic and heterogeneous environments has been a challenge since the last decades. However, recent advances in near-surface geophysics have facilitated the use of geophysical data for hydrogeological characterization in the last few years. A pseudo 3-D high resolution P-wave shallow seismic reflection survey was performed at the P Reactor Area, Savannah River Site, South Carolina in order to delineate and predict migration pathways of a large contaminant plume including trichloroethylene. This contaminant plume originates from the northwest section of the reactor facility that is located within the Upper Atlantic Coastal Plain. The data were collected with 40 Hz geophones, an accelerated weight-drop as seismic source and 1 m receiver spacing with near- and far-offsets of 0.5 and 119.5 m, respectively. In such areas with near-surface contaminants, a detailed subsurface characterization of the vadose zone hydraulic parameters is very important. Indeed, an inexpensive method of deriving such parameters by the use of seismic reflection surveys is beneficial, and our approach uses the relationship between seismic velocity and hydrogeological parameters together with empirical observations relating porosity to permeability and hydraulic conductivity. Shear wave velocity ( V s ) profiles were estimated from surface wave dispersion analysis of the seismic reflection data and were subsequently used to derive hydraulic parameters such as porosity, permeability, and hydraulic conductivity. Additional geophysical data including core samples, vertical seismic profiling, surface electrical resistivity tomography, natural gamma and electrical resistivity logs allowed for a robust assessment of the validity and geological significance of the estimated V s and hydrogeological models. The results demonstrate the usefulness of this approach for the upper 15 m of shallow unconsolidated sediments even though the survey design parameters were not optimal for surface wave analysis due to the higher than desired frequency geophones.

Engineering↗

DOE OSTI2026

Wyoming CarbonSAFE Phase III: Site Characterization and Permitting Commercial-Scale Carbon Storage Complex Feasibility Study at Dry Fork Station, Wyoming

This report presents the findings of the technical and non-technical site characterization and permitting activities (“Phase III”) conducted under the Wyoming CarbonSAFE: Accelerating CCUS Commercialization and Deployment at Dry Fork Power Station (DFS) and the Wyoming Integrated Test Center project (“Wyoming CarbonSAFE”). Wyoming CarbonSAFE is part of the U.S. Department of Energy’s (DOE) National Energy Technology Laboratory (NETL) Carbon Storage Assurance Facility Enterprise (“CarbonSAFE”) Initiative. The results of Phase III demonstrate that the Wyoming CarbonSAFE project - referred to as the Northern Powder River Basin Carbon Sequestration Hub (NPRB-CSH) - meets the technical, regulatory, and commercial feasibility requirements necessary to advance toward commercial development and construction. The activities completed under this project Phase make the NPRB-CSH one of the region’s most commercially ready carbon storage sites. Completion of this Phase included the finalization of all site surface and subsurface characterization activities, completion and testing of two Class VI standard wells, 10 draft Class VI permits-to construct to address the future needs of a storage complex, finalized NEPA assessments, transportation and capture FEED studies, and a full commercialization strategy with economic modeling, operation and site closure strategies. The NPRB-CSH meets all requirements to progress to a CarbonSAFE Phase IV program or advance to full commercial operations under the development of a business partner.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2018

Geophysical Imaging of Geothermal Systems Spanning Various Geologic Settings

Despite the wide range of geologic settings in which geothermal systems reside, many systems share characteristics that can be targeted by geophysical techniques. Over the last few years, the authors have used potential field and magnetotelluric (MT) methods to jointly image and characterize numerous geothermal systems in various geologic and tectonic settings. This includes systems in the Basin and Range, Washington Cascade Range, Northeastern Oregon, and the Snake River Plain. The principal geothermal targets are areas of elevated temperatures near fault interactions that promote permeability and fluid flow. Gravity, magnetic, and MT data are all sensitive to different physical parameters that together form a complementary set of tools for characterizing subsurface structure and identifying favorable geothermal targets. In general, potential fields delineate important geologic structures and MT locates fluids and zones of hydrothermal alteration. Joint interpretation of these complementary methods, along with other geologic information, provides a robust framework for geothermal prospecting and understanding a hydrothermal system. Specifically, common geothermal targets lie along potential field gradients and are associated with an electrically conductive anomaly that resides in the less-dense, less-magnetic material.

15 GEOTHERMAL ENERGY↗

DOE OSTI2019

Extending accurate spectral induced polarization measurements into the kHz range: modelling and removal of errors from interactions between the parasitic capacitive coupling and the sample holder

SUMMARY Spectral induced polarization (SIP), describing the measurement of the frequency domain electrical impedance magnitude and phase of porous materials, has been widely used to characterize subsurface hydrological/biogeochemical properties and processes. SIP data collected at frequencies higher than 100 Hz are expected to describe the polarization of small particles providing insights into the physicochemical properties of clays, nanoparticles and microorganisms. However, the phase measurements at these high frequencies are often contaminated by errors due to the parasitic capacitive coupling of the SIP instrument, especially for lower conductivity samples. We developed a model showing the measured phase is the sum of the true sample phase and an error term ωCinZx, where ω is the angular frequency; Cin is the instrument input capacitance and Zx is a measurable impedance function related to the sample holder properties and the reference resistor. Based on this model, a new phase correction method is proposed that results in highly accurate SIP data up to 20 kHz as well as the determination of Cin. We tested the correction method using electric circuits, NaCl fluids and three unconsolidated samples (sand, sand-clay and sand-pyrite mixtures). The corrected phase for the circuit and NaCl fluid experiments showed excellent agreement with the theoretical phase response across the studied frequency range (errors <1 mrad). For unconsolidated samples, removal of errors results in phase spectra more consistent with expected polarization mechanisms, as based on phase peaks recorded for small pyrite and clay particles at high frequencies. These phase peaks could not be identified in the uncorrected data. Our approach can substantially enhance the value of the SIP method for the characterization of fine-grained sediments and rocks.

58 GEOSCIENCES↗

DOE OSTI2016

High-stress triaxial direct-shear fracturing of Utica shale and in situ X-ray microtomography with permeability measurement

The challenge of characterizing subsurface fluid flow has motivated extensive laboratory studies, yet fluid flow through rock specimens in which fractures are created and maintained at high-stress conditions remains underinvestigated at this time. The studies of this type that do exist do not include in situ fracture geometry measurements acquired at stressed conditions, which would be beneficial for interpreting the flow behavior. Therefore, this study investigates the apparent permeability induced by direct-shear fracture stimulation through Utica shale (a shale gas resource and potential caprock material) at high triaxial stress confinement and for the first time relates these values to simultaneously acquired in situ X-ray radiography and microtomography images. Change in fracture geometry and apparent permeability was also investigated at additional reduced stress states. Finite element and combined finite-discrete element modeling were used to evaluate the in situ observed fracturing process. Results from this study indicate that the increase in apparent permeability through fractures created at high-stress (22.2 MPa) was minimal relative to the intact rock (less than 1 order of magnitude increase), while fractures created at low stress (3.4 MPa) were significantly more permeable (2 to 4 orders of magnitude increase). Here, this study demonstrates the benefit of in situ X-ray observation coupled with apparent permeability measurement to analyze fracture creation in the subsurface. Our results show that the permeability induced by fractures through shale at high stress can be minor and therefore favorable in application to CO 2 sequestration caprock integrity but unfavorable for hydrocarbon recovery from unconventional reservoirs.

58 GEOSCIENCES↗

DOE OSTI2016

A computationally efficient parallel Levenberg-Marquardt algorithm for highly parameterized inverse model analyses

Inverse modeling seeks model parameters given a set of observations. However, for practical problems because the number of measurements is often large and the model parameters are also numerous, conventional methods for inverse modeling can be computationally expensive. We have developed a new, computationally-efficient parallel Levenberg-Marquardt method for solving inverse modeling problems with a highly parameterized model space. Levenberg-Marquardt methods require the solution of a linear system of equations which can be prohibitively expensive to compute for moderate to large-scale problems. Our novel method projects the original linear problem down to a Krylov subspace, such that the dimensionality of the problem can be significantly reduced. Furthermore, we store the Krylov subspace computed when using the first damping parameter and recycle the subspace for the subsequent damping parameters. The efficiency of our new inverse modeling algorithm is significantly improved using these computational techniques. We apply this new inverse modeling method to invert for random transmissivity fields in 2D and a random hydraulic conductivity field in 3D. Our algorithm is fast enough to solve for the distributed model parameters (transmissivity) in the model domain. The algorithm is coded in Julia and implemented in the MADS computational framework (http://mads.lanl.gov). By comparing with Levenberg-Marquardt methods using standard linear inversion techniques such as QR or SVD methods, our Levenberg-Marquardt method yields a speed-up ratio on the order of ~10 1 to ~10 2 in a multi-core computational environment. Furthermore, our new inverse modeling method is a powerful tool for characterizing subsurface heterogeneity for moderate- to large-scale problems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

An integrated approach to derive relative permeability from capillary pressure

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.

58 GEOSCIENCES↗

DOE OSTI2017

High-Dimensional Intrinsic Interpolation Using Gaussian Process Regression and Diffusion Maps

This article considers the challenging task of estimating geologic properties of interest using a suite of proxy measurements. The current work recast this task as a manifold learning problem. In this process, this article introduces a novel regression procedure for intrinsic variables constrained onto a manifold embedded in an ambient space. The procedure is meant to sharpen high-dimensional interpolation by inferring non-linear correlations from the data being interpolated. The proposed approach augments manifold learning procedures with a Gaussian process regression. It first identifies, using diffusion maps, a low-dimensional manifold embedded in an ambient high-dimensional space associated with the data. It relies on the diffusion distance associated with this construction to define a distance function with which the data model is equipped. This distance metric function is then used to compute the correlation structure of a Gaussian process that describes the statistical dependence of quantities of interest in the high-dimensional ambient space. The proposed method is applicable to arbitrarily high-dimensional data sets. Here, it is applied to subsurface characterization using a suite of well log measurements. The predictions obtained in original, principal component, and diffusion space are compared using both qualitative and quantitative metrics. Considerable improvement in the prediction of the geological structural properties is observed with the proposed method.

58 GEOSCIENCES↗

DOE OSTI2014

Krypton-81 in groundwater of the Culebra Dolomite near the Waste Isolation Pilot Plant, New Mexico

TheWaste Isolation Pilot Plant (WIPP) in New Mexico is the first geologic repository for disposal of transuranic nuclear waste fromdefense-related programs of the US Department of Energy. It is constructed within halite beds of the Permian-age Salado Formation. The Culebra Dolomite, confined within Rustler Formation evaporites overlying the Salado Formation, is a potential pathway for radionuclide transport from the repository to the accessible environment in the human-disturbed repository scenario. Although extensive subsurface characterization and numerical flow modeling of groundwater has been done in the vicinity of theWIPP, few studies have used natural isotopic tracers to validate the flow models and to better understand solute transport at this site. The advent of Atom-Trap Trace Analysis (ATTA) has enabled routine measurement of cosmogenic 81Kr (half-life 229,000 yr), a near-ideal tracer for long-term groundwater transport. We measured 81Kr in saline groundwater sampled from two Culebra Dolomite monitoring wells near the WIPP site, and compared 81Kr model ages with reverse particle-tracking results of well-calibrated flow models. The 81Kr model ages are ~130,000 and ~330,000 yr for high-transmissivity and low-transmissivity portions of the formation, respectively. Compared with flow model results which indicate a relatively young mean hydraulic age (~32,000 yr), the 81Kr model ages imply substantial physical attenuation of conservative solutes in the Culebra Dolomite and provide limits on the effective diffusivity of contaminants into the confining aquitards.

Sturchio, Neil C.↗

DOE OSTI2022

Photothermal Atomic Force Microscopy Coupled with Infrared Spectroscopy (AFM-IR) Analysis of High Extinction Coefficient Materials: A Case Study with Silica and Silicate Glasses

Photothermal atomic force microscopy coupled with infrared spectroscopy (AFM-IR) brings significant value as a spatially resolved surface analysis technique for disordered oxide materials such as glasses, but additional development and fundamental understanding of governing principles is needed to interpret AFM-IR spectra, since the existing theory described for organic materials does not work for materials with high extinction coefficients for infrared (IR) absorption. This paper describes theoretical calculation of a transient temperature profile inside the IR-absorbing material considering IR refraction at the interface as well as IR adsorption and heat transfer inside the sample. Furthermore, this calculation explains the differences in peak positions and amplitudes of AFM-IR spectra from those of specular reflectance and extinction coefficient spectra. It also addresses the information depth of the AFM-IR characterization of bulk materials. AFM-IR applied to silica and silicate glass surfaces has demonstrated novel capability of characterizing subsurface structural changes and surface heterogeneity due to mechanical stresses from physical contacts, as well as chemical alterations manifested in surface layers through aqueous corrosion.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

DeFault: DEep‐Learning‐Based FAULT Delineation Using the IBDP Passive Seismic Data at the Decatur CO2 Storage Site

Abstract The carbon capture, utilization, and storage (CCUS) framework is an essential component in reducing greenhouse gas emissions, with its success hinging on the comprehensive knowledge of subsurface geology and geomechanics. Passive seismic event relocation and fault detection offer vital insights into subsurface structures and the ability to monitor fluid migration pathways. Accurate identification and localization of seismic events, however, face significant challenges, including the necessity for high‐quality seismic data and advanced computational methods. To address these challenges, we introduce a novel deep learning method, , specifically designed for passive seismic source relocation and fault delineating for passive seismic monitoring projects. By leveraging data domain‐adaptation, allows us to train a neural network with labeled synthetic data and apply it directly to field data. Using , the passive seismic sources are automatically clustered based on their recording time and spatial locations, and subsequently, faults and fractures are delineated accordingly. We demonstrate the efficacy of on a field case study involving injection related microseismic data from Decatur, Illinois area. Our approach accurately and efficiently relocated passive seismic events, identified faults and could aid in potential damage induced by seismicity. Our results highlight the potential of as a valuable tool for passive seismic monitoring, emphasizing its role in ensuring CCUS project safety. This research bolsters the understanding of subsurface characterization in CCUS, illustrating machine learning’s capacity to refine these methods. Ultimately, our work has significant implications for CCUS technology deployment, an essential strategy in combating climate change. Plain Language Summary In our quest to tackle climate change, we use a strategy known as carbon capture, utilization, and storage (CCUS) to keep greenhouse gases out of the atmosphere. This strategy relies heavily on our ability to understand what's happening deep under the earth's surface. To make sure we store super critical safely, we need to accurately map out the geological structure, especially faults, but this is tough without high‐quality data and complex computer programs. We've developed a new tool called “DeFault,” which uses advanced machine learning to improve how we find and map these underground features. “DeFault” is smart enough to learn from numerically simulated data and then apply what it’s learned to real‐world situations. It groups together seismic activity—tiny tremors and shifts in the earth—based on when and where they happen, which helps us spot where there might be cracks or faults. We tested “DeFault” in Illinois, where CO 2 is injected underground, and it successfully pinpointed where these tremors occurred and mapped out the faults, helping to prevent accidents accurately in the future. Our study shows that “DeFault” will be a powerful ally in making CCUS safer and more effective, especially for the Illinois Basin Decatur Project. Key Points Faults and fractures introduced by carbon storage can be monitored by passive seismicity DeFault algorithm enables an automatic process for accurate and efficient passive seismic event locating and clustering

58 GEOSCIENCES↗

DOE OSTI2020

The emergence of small-scale self-affine surface roughness from deformation

Most natural and man-made surfaces appear to be rough on many length scales. There is presently no unifying theory of the origin of roughness or the self-affine nature of surface topography. One likely contributor to the formation of roughness is deformation, which underlies many processes that shape surfaces such as machining, fracture, and wear. Using molecular dynamics, we simulate the biaxial compression of single-crystal Au, the high-entropy alloy Ni 36.67 Co 30 Fe 16.67 Ti 16.67 , and amorphous Cu 50 Zr 50 and show that even surfaces of homogeneous materials develop a self-affine structure. By characterizing subsurface deformation, we connect the self-affinity of the surface to the spatial correlation of deformation events occurring within the bulk and present scaling relations for the evolution of roughness with strain. These results open routes toward interpreting and engineering roughness profiles.

42 ENGINEERING↗

DOE OSTI2016

The direct-current response of electrically conducting fractures excited by a grounded current source

Hydraulic fracture stimulation of low permeability reservoir rocks is an established and cross–cutting technology for enhancing hydrocarbon production in sedimentary formations and increasing heat exchange in crystalline geothermal systems. Whereas the primary measure of success is the ability to keep the newly generated fractures sufficiently open, long–term reservoir management requires a knowledge of the spatial extent, morphology, and distribution of the fractures — knowledge primarily informed by microseismic and ground deformation monitoring. To minimize the uncertainty associated with interpreting such data, we investigate through numerical simulation the usefulness of direct-current (DC) resistivity data for characterizing subsurface fractures with elevated electrical conductivity by considering a geophysical experiment consisting of a grounded current source deployed in a steel cased borehole. In doing so, the casing efficiently energizes the fractures with steady current. Finite element simulations of this experiment for a horizontal well intersecting a small set of vertical fractures indicate that the fractures manifest electrically in (at least) two ways: (1) a local perturbation in electric potential proximal to the fracture set, with limited farfield expression and (2) an overall reduction in the electric potential along the borehole casing due to enhanced current flow through the fractures into the surrounding formation. The change in casing potential results in a measurable effect that can be observed far from fractures themselves. Under these conditions, our results suggest that farfield, timelapse measurements of DC potentials can be interpreted by simple, linear inversion for a Coulomb charge distribution along the borehole path, including a local charge perturbation due to the fractures. As a result, this approach offers an inexpensive method for detecting and monitoring the time-evolution of electrically conducting fractures while ultimately providing an estimate of their effective conductivity — the latter providing an important measure independent of seismic methods on fracture shape, size, and hydraulic connectivity.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2022

Observing and modeling the effects of production infrastructure in electromagnetic surveys

Electromagnetic (EM) methods are among the original techniques for subsurface characterization in exploration geophysics because of their particular sensitivity to the earth electrical conductivity, a physical property of rocks distinct yet complementary to density, magnetization, and strength. However, this unique ability also makes them sensitive to metallic artifacts — infrastructure such as pipes, cables, and other forms of cultural clutter — the EM footprint of which often far exceeds their diminutive stature when compared to that of bulk rock itself. In the hunt for buried treasure or unexploded ordnance, this is an advantage; in the long-term monitoring of mature oil fields after decades of production, it is quite troublesome indeed. In this study we consider the latter through the lens of an evolving energy industry landscape in which the traditional methods of EM characterization for the exploration geophysicist are applied toward emergent problems in well-casing integrity, carbon capture and storage, and overall situational awareness in the oil field. We introduce case studies from these exemplars, showing how signals from metallic artifacts can dominate those from the target itself and impose significant burdens on the requisite simulation complexity. Furthermore, we show how recent advances in numerical methods mitigate the computational explosivity of infrastructure modeling, providing feasible and real-time analysis tools for the desktop geophysicist. Lastly, we demonstrate through comparison of field data and simulation results that incorporation of infrastructure into the analysis of such geophysical data is, in a growing number of cases, a requisite but now manageable step.

58 GEOSCIENCES↗

DOE OSTI2020

GBCGE Subsurface Database Explorer and APIs

This submission defines a DOI for the Great Basin Center for Geothermal Energy's (GBCGE) Subsurface Database Explorer web application and underlying data services, and acknowledges the INGENIOUS project as a major source of funding for data compilation and quality assurance. The GBCGE Subsurface Database Explorer is an interactive web mapping application that provides public access to the GBCGE Subsurface Database, and its collection of datasets pertinent to geothermal exploration, oil and gas exploration, critical mineral exploration, and other subsurface characterization for the Great Basin Region, western US. This is a living database, and will be continuously updated with new data and datasets as funding and motivations allow. The underlying database views that populate the web application are on an automated refresh schedule. Data sources and acknowledgements: We thank our partners with the Nevada Division of Minerals (NDOM), the Southern Methodist University (SMU), and Great Basin State Geological Surveys for their active efforts in data curation, schema design, and quality assurance. We also thank contributors among the USGS, Oregon Institute of Technology, State Divisions of Water Resources, State Divisions of Oil, Gas, and Minerals, and State Geological Surveys for open data availability and direct contributions made under the National Geothermal Data System (NGDS).

15 GEOTHERMAL ENERGY↗

DOE OSTI2015

Inverse problems in heterogeneous and fractured media using peridynamics

The following work presents an adjoint-based methodology for solving inverse problems in heterogeneous and fractured media using state-based peridynamics. We show that the inner product involving the peridynamic operators is self-adjoint. The proposed method is illustrated for several numerical examples with constant and spatially varying material parameters as well as in the context of fractures. We also present a framework for obtaining material parameters by integrating digital image correlation (DIC) with inverse analysis. This framework is demonstrated by evaluating the bulk and shear moduli for a sample of nuclear graphite using digital photographs taken during the experiment. The resulting measured values correspond well with other results reported in the literature. Lastly, we show that this framework can be used to determine the load state given observed measurements of a crack opening. Furthermore, this type of analysis has many applications in characterizing subsurface stress-state conditions given fracture patterns in cores of geologic material.

97 MATHEMATICS AND COMPUTING↗