Search NASA⌕ Search

SEARCH · Search NASA

Results for “borehole geophysics”

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

Integrated hydrogeophysical modelling and data assimilation for geoelectrical leak detection

Time-lapse electrical resistivity tomography (ERT) measurements provide indirect observations of hydrological processes in the Earth's shallow subsurface at high spatial and temporal resolution. ERT has been used in the past decades to detect leaks and monitor the evolution of associated contaminant plumes. Specifically, inverted resistivity images allow visualization of the dynamic changes in the structure of the plume. However, existing methods do not allow the direct estimation of leak parameters (e.g. leak rate, location, etc.) and their uncertainties. We propose an ensemble-based data assimilation framework that evaluates proposed hydrological models against observed time-lapse ERT measurements without directly inverting for the resistivities. Each proposed hydrological model is run through the parallel coupled hydro-geophysical simulation code PFLOTRAN-E4D to obtain simulated ERT measurements. The ensemble of model proposals is then updated using an iterative ensemble smoother. In this paper, we demonstrate the proposed framework on synthetic and field ERT data from controlled tracer injection experiments. Our results show that the approach allows joint identification of contaminant source location, initial release time, and solute loading from the cross-borehole time-lapse ERT data, alongside with an assessment of uncertainties in these estimates. We demonstrate a reduction in site-wide uncertainty by comparing the prior and posterior plume mass discharges at a selected image plane. This framework is particularly attractive to sites that have previously undergone extensive geological investigation (e.g., nuclear sites). It is well suited to complement ERT imaging and we discuss practical issues in its application to field problems.

58 GEOSCIENCES↗
DOE OSTI2023

Monitoring Preferential Flow of Water in Sand Using Thermoacoustics Wave Imaging

Abstract Accurate predictions of fluid flow, mass transport, and reaction rates critically impact the efficiency and reliability of subsurface exploration and sustainable use of subsurface resources. Quantitative dynamical sensing and imaging can play a pivotal role in the ability to make such predictions. Geophysical thermoacoustic technology has the potential to provide the aforementioned capabilities, since it builds upon the principle that electromagnetic and mechanical wave fields can be coupled through a thermodynamic process. In this letter, we present laboratory experiments featuring the efficacy of thermoacoustic imaging in the monitoring of preferential flow of water in porous media. Our laboratory experimental equipment can be readily packaged in a form factor that fits in a borehole, and the use of multiple acoustic transducers—which can be combined with volumetric coding techniques—has the potential to provide quasi‐real‐time imaging (0.5 Hz video rate) of regions in close proximity (a few meters) of an open field well.

Geology↗
DOE OSTI2020

Hydraulic fracturing experiments at 1500 m depth in a deep mine: Highlights from the kISMET project

In support of the U.S. DOE SubTER Crosscut initiative, we established a field test facility in a deep mine and designed and carried out in situ hydraulic fracturing experiments relevant to enhanced geothermal systems (EGS) in crystalline rock to characterize the stress field, understand the effects of rock fabric on fracturing, and gain experience in monitoring using geophysical methods. The project also included pre- and post-fracturing simulation and analysis, and laboratory measurements and experiments. The kISMET (permeability (k) and Induced Seismicity Management for Energy Technologies) site was established in the West Access Drift of the Sanford Underground Research Facility (SURF) 4757 ft (1450 m) below ground (on the 4850 ft level (4850L)) in phyllite of the Precambrian Poorman Formation. We drilled and continuously cored five near-vertical boreholes in a line on 3 m (10 ft) spacing, deviating the two outermost boreholes slightly to create a five-spot pattern around the test borehole centered in the test volume 40 m below the drift invert (floor) at a total depth of ~1490 m (4890 ft). Laboratory measurements of core from the center test borehole showed P-wave velocity heterogeneity along each core indicating strong, fine-scale (~1 cm or smaller) changes in the mechanical properties of the rock. Field measurements of the stress field by hydraulic fracturing showed that the minimum horizontal stress at the kISMET site averages 21.7 MPa (3146 psi) trending approximately N-S (356 degrees azimuth) and plunging slightly NNW at 12°. The vertical and horizontal maximum stresses are similar in magnitude at 42-44 MPa (6090-6380 psi) for the depths of testing, which averaged approximately 1530 m (5030 ft). Hydraulic fractures were remarkably uniform suggesting core-scale and larger rock fabric did not play a role in controlling fracture orientation. Analytical solutions suggest that the fracture radius of the large fracture (stimulation test) was more than 6 m (20 ft), depending on the unknown amount of leak-off.

Oldenburg, C↗
DOE OSTI2024

Interpretation of large‐scale, long‐term electrical geophysical monitoring guided by a process simulation

Abstract Surface electrical resistivity tomography (ERT) was used at a waste site to monitor vadose zone changes in electrical properties as a proxy for contaminant flux over a span of 17 years. The BC Cribs and Trenches (BCCT) site at the Hanford site contains 20 disposal trenches and six disposal cribs. Wastes include a large inventory of technetium‐99 and large masses of nitrate and uranium‐238. ERT data were collected along 41 profiles in 2005 to characterize regions of elevated bulk electrical conductivity (BEC) associated with past liquid waste discharges. Previous analyses performed on samples from four boreholes showed a high correlation between nitrate concentration and BEC. In 2022, ERT data were re‐collected along the same profiles and six additional profiles in an area not previously surveyed. Compared to background uncontaminated areas, BEC was higher in contaminated areas at the waste sites. Given the correlation between nitrate concentration and BEC previously found at this site, ERT images show the spatial distribution and relative ionic concentration of vadose zone contaminants at BCCT. Between 2005 and 2022, ERT difference images showed a decrease in BEC surrounding most waste sites, with exceptions where there were known anthropogenic surface changes. An evaluation of recharge‐driven nitrate migration using synthetic flow and transport simulations showed that downward migration causes a decrease in BEC from the decrease in ionic strength at the trailing end of the plume where contaminants migrated downward. From this, we interpret ERT difference images as showing the predominant regions of downward ion flux.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

High-resolution Bayesian spatial autocorrelation (SPAC) quasi-3-D Vs model of Utah FORGE site with a dense geophone array

SUMMARY We expand the application of spatial autocorrelation (SPAC) from typical 1-D Vs profiles to quasi-3-D imaging via Bayesian Monte Carlo inversion (BMCI) using a dense nodal array (49 nodes) located at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site. Combinations of 4 and 9 geophones in subarrays provide for 36 and 25 1-D Vs profiles, respectively. Profiles with error bars are determined by calculating coherency functions that fit observations in a frequency range of 0.2–5 Hz. Thus, a high-resolution quasi-3-D Vs model from the surface to 2.0 km depth is derived and shows that surface-parallel sedimentary strata deepen to the west, consistent with a 3-D seismic reflection survey. Moreover, the resulting Vs profile is consistent with a Vs profile derived from distributed acoustic sensing (DAS) data located in a borehole at the FORGE site. The quasi-3-D velocity model shows that the base of the basin dips ∼22° to the west and topography on the basement interface coincident with the Mag Lee Wash suggests that the bedrock interface is an unconformity.

Geochemistry & Geophysics↗
DOE OSTI2023

Autonomous Inversion of In Situ Deformation Measurement Data for Injection-Induced Stress Change

Geologic carbon storage (GCS) is likely to play a key part of the global effort to dramatically reduce CO2 emissions and perhaps even reduce atmospheric CO2 concentrations through carbon negative operations. A critical part of effort to commercialize and widely deploy this technology is developing the capability to rapidly assimilate real-time monitoring data into a form that will enable site operators to make decisions to manage the safe and efficient operations. Two of the risks associate with GCS are the risk of inducing fractures in the sealing formations that can create leakage pathways and the risk of inducing earthquakes of sufficient magnitude to cause public concern, property damage, or safety risks. To properly manage these risks the site operator needs to know the initial state of stress, the change in stress induced by injection, and the relationship between operational parameters such as injection rate and pressure and the change in stress. Current methods of estimating the change in stress require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the characterization phase, with little feedback from operational observations to validate or refine these choices. These characterization methods interrogate the geologic formations using length scales, loading rates or magnitudes that are quite different from those encountered by the actual storage system. It is shown that errors in the assumed constitutive response, even when informed by laboratory tests on core samples, are likely to be common, large, and underestimate the magnitude of stress change caused by injection. Recent advances in borehole-based strain instruments and borehole and surface-based tilt and displacement instruments have now enabled monitoring of the deformation of the storage system throughout its operational lifespan. This data can enable validation and refinement of the knowledge of the geomechanical properties and state of the system, but brings with it a challenge to transform the raw data into actionable knowledge. We demonstrate a method that uses automatic differentiation and a finite-element based geomechanical model perform a gradient-based deterministic inversion of geomechanical monitoring data. This approach allows autonomous integration of the instrument data without the need for time consuming manual interpretation and selection of updated model parameters. Furthermore, only isotropic linear elasticity is considered in this paper, the approach presented is very flexible as to what type of geomechanical constitutive response can be used. The approach is easily adaptable to nonlinear physics-based constitutive models to account for common rock behaviors such as creep and plasticity. The approach also enables training of machine learning-based constitutive models by allowing back propagation of errors through the finite element calculations. This enables strongly enforcing known physics, such as conservation of momentum and continuity, while allowing data-driven models to learn the truly unknown physics such as the constitutive or petrophysical responses.

Burghardt, Jeffrey A.↗
DOE OSTI2021

EGS Collab Experiment 2: Testbed Well Orientation Data

This dataset contains the north seeking gyro data for each of the 11 boreholes drilled at the Experiment 2 testbed on the 4100 foot level of the SURF (Sanford Underground Research Facility). Each gyro file is in individual folders for each well. A single folder called 'E2 All Well Trajectories' contains the trajectories for all the wells but in the mine coordinate system. An image of the well layout is provided in a PowerPoint.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Peak Ground Velocity Spatial Variability Revealed by Dense Seismic Array in Southern California

There are an infinite number of locations where seismic waves interact with Earth's surface, yet seismologists can only record a finite set of ground motions. This requires assumptions about how known measurements or predictions can be extrapolated to nearby locations. We test to what extent extrapolation is viable using an array of 1,108 geophones installed 10–20 m apart within a 0.36 km 2 study area for 1 month atop the San Jacinto Fault in Southern California. Measuring ground motions from 38 small local earthquakes, we observe peak ground velocity (PGV) variations up to 22% across the array. We find the basin structure consistently traps energy, amplifying the seismic velocities, and that the damaged region within the fault reduces the PGVs. Comparing surface and borehole measurements (148 m deep), we find PGVs recorded at the surface are 3–10 times greater than the measurements at depth. We attribute these differences to amplification from soft sediments within the basin. We show that by leveraging information from >1,000 measurements per earthquake the ground motion velocities can be tightly constrained. These results suggest the spatial extrapolation of PGVs in regions of known faults and basins, even across a small area, should be done with caution.

58 GEOSCIENCES↗
DOE OSTI2020

A New Capability of E4D For 3D Parallel Joint Inversion of DC Resistivity And Traveltime Data on Unstructured Mesh

A major challenge in interpreting geophysical data is how to derive consistent three-dimensional (3D) earth models of different physical properties from spatially and temporally limited measurements. Joint inversion with cross-gradient constraints is an approach to find such models by imposing structural similarities between different physical parameters. We have developed a parallel distributed-memory joint inversion code for direct-current (DC) resistivity and traveltime data using the cross-gradient constraint on unstructured mesh. The code utilizes existing E4D framework for parallel forward simulation, distributed storage and computation of the Jacobian matrix of forward operator, and parallel execution of matrix-vector multiplication during inversion. Besides, the joint inversion is solved by nonlinear conjugate gradient algorithm parallelized for DC resistivity and traveltime data. The joint inversion capability of E4D was tested using synthetic data from cross-borehole DC resistivity and traveltime data. The results indicate that the shape and size of the anomalies from the joint inversion are more reliable than those from separate inversions.

58 GEOSCIENCES↗
DOE OSTI2020

Integration of Soft Data Into Geostatistical Simulation of Categorical Variables

Uncertain or indirect “soft” data, such as geologic interpretation, driller’s logs, geophysical logs or imaging, offer potential constraints or “soft conditioning” to stochastic models of discrete categorical subsurface variables in hydrogeology such as hydrofacies. Previous bivariate geostatistical simulation algorithms have not fully addressed the impact of data uncertainty in formulation of the (co) kriging equations and the objective function in simulated annealing (or quenching). This paper introduces the geostatistical simulation code tsim-s, which accounts for categorical data uncertainty through a data “hardness” parameter. In generating geostatistical realizations with tsim-s, the uncertainty inherent to soft conditioning is factored into both 1) the data declustering and spatial correlation functions in cokriging and 2) the acceptance probability for change of category in simulated quenching. The degree or sensitivity to which soft data conditions a realization as a function of hardness can be quantified by mapping category probabilities derived from multiple realizations. In addition to point or borehole data, arrays of data (e.g., as derived from a depth-dependency function, probability map, or “prior realization”) can be used as soft conditioning. The tsim-s algorithm provides a theoretically sound and general framework for integrating datasets of variable location, resolution, and uncertainty into geostatistical simulation of categorical variables. A practical example shows how tsim-s is capable of generating a large-scale three-dimensional simulation including curvilinear features.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Utah FORGE: Seismic Velocity Models, February 2021

This dataset contains a map, showing the Utah FORGE seismic stations, and seismic velocity model data. There are 61 1-D velocity models which are in a compressed TAR file. A paper is referenced at the end of this description which discusses the use of these data in 3D modelling. The paper summary follows: We expand the application of spatial autocorrelation (SPAC) from typical 1-D Vs profiles to quasi-3-D imaging via Bayesian Monte Carlo inversion (BMCI) using a dense nodal array (49 nodes) located at the Utah Frontier Observatory for Research in Geothermal Energy (FORGE) site. Combinations of 4 and 9 geophones in subarrays provide for 36 and 25 1-D Vs profiles, respectively. Profiles with error bars are determined by calculating coherency functions that fit observations in a frequency range of 0.2-5 Hz. Thus, a high-resolution quasi-3-D Vs model from the surface to 2.0 km depth is derived and shows that surface-parallel sedimentary strata deepen to the west, consistent with a 3-D seismic reflection survey. Moreover, the resulting Vs profile is consistent with a Vs profile derived from distributed acoustic sensing (DAS) data located in a borehole at the FORGE site. The quasi-3-D velocity model shows that the base of the basin dips ~22 degrees to the west and topography on the basement interface coincident with the Mag Lee Wash suggests that the bedrock interface is an unconformity. Reference: Zhang, H. and K. L. Pankow (2021). High-resolution Bayesian spatial auto-correlation (SPAC) pseudo-3D Vs model of Utah FORGE site with a dense geophone array, Geophys. Res. Int, https://doi.org/10.1093/gji/ggab049

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Shallow EGS Regional Resource Potential and Map - Snake River Plain

SMU Geothermal Lab developed a methodology to estimate shallow (1 km to 4 km) Enhanced Geothermal Systems (EGS) resource potential using an approach that utilizes recent geology and geophysical research along with new well data to improve the thermal conductivity model, mitigate impacts from groundwater flow in the thermal model, and examine radioactivity data variations. By incorporating the results of the most recent projects with the SMU shallow methodology, we developed a more accurate, updated resource estimate for the Snake River Plain (SRP). The resulting maps and resource estimates can be used by the National Renewable Energy Lab (NREL), Bureau of Land Management (BLM), and the public to determine how best to move forward with future project development.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Cubic-meter scale laboratory fault re-activation experiments to improve the understanding of induced seismicity risks

Abstract To understand fluid induced seismicity, we have designed a large-scale laboratory experiment consisting of a one-cubic-meter sandstone with an artificial fault cut and fluid-injection boreholes. The sandstone block is assembled in a true triaxial loading frame and equipped with 38 piezoelectric sensors to locate and characterise acoustic emission events. The differential stress on the artificial fault is increased in stages to bring it towards a critically stressed state. After each stage of differential stress increase, fluids are injected at low pressures through boreholes to test the potential of fault re-activation. In addition, a high-pressure injection was conducted that created a hydraulic fracture from the injection borehole towards the artificial fault. The newly generated fluid pathway resulted in an activation of the complete block through a stick–slip movement. We compare acoustic emission measurements from the laboratory experiment with seismicity observations from the field-scale CO 2 injection at Decatur, Illinois, U.S., and conclude that the existence of fluid pathways plays a decisive role for the potential of induced seismicity.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Estimation of the orientation of stress in the Earth’s crust without earthquake or borehole data

Abstract Mechanical stress acting in the Earth’s crust is a fundamental property that is important for a wide range of scientific and engineering applications. The orientation of maximum horizontal compressive stress can be estimated by inverting earthquake source mechanisms and measured directly from borehole-based measurements, but large regions of the continents have few or no observations. Here we present an approach to determine the orientation of maximum horizontal compressive stress by measuring stress-induced anisotropy of nonlinear susceptibility, which is the derivative of elastic modulus with respect to strain. Laboratory and Earth experiments show that nonlinear susceptibility is azimuthally dependent in an anisotropic stress field and is maximum in the orientation of maximum horizontal compressive stress. We observe this behavior in the Earth—in Oklahoma and New Mexico, U.S.A, where maximum nonlinear susceptibility coincides with the orientation of maximum horizontal compressive stress measured using traditional methods. Our measurements use empirical Green’s functions and solid-earth tides and can be applied at different temporal and spatial scales.

58 GEOSCIENCES↗
DOE OSTI2024

Stress Field Dynamics and Fault Slip Potential in the Paradox Basin

Abstract The Paradox Basin, straddling Utah, Colorado, Arizona, and New Mexico is characterized by an intricate amalgamation of evaporites and clastic layers and is dominated by prominent salt walls and related subsurface structures. Our research offers a new examination of the stress distribution across the basin, deriving from continuous and discrete stress measurements conducted in boreholes in the region and focal mechanism analysis, emphasizing variations over salt structures. Integrating Coulomb failure criteria with probabilistic methods, we assess potential fault movements resulting from fluid pressure alterations. Our approach provides a comprehensive understanding of the Paradox Basin's state of stress, showing a continuous change of the maximum horizontal stress orientation from N‐S at the Wasatch Fault Zone to WNW‐ESE in the northern part of the Paradox Basin and to WSW‐ENE in the southern part of the basin. Further East, into the Colorado Plateau and the Uncompahgre Uplift, the S H max orientation becomes E‐W. Decoding stress orientation dynamics has enabled critical insights into fault slip potential, especially in the basin's northern region. The salt wall faults are less likely to slip, and the Paradox Formation's evaporite and clastic rock sequence can serve as a potential low seismic risk target for carbon storage and hydrocarbon extraction.

Geochemistry & Geophysics↗
DOE OSTI2021

Identifying geological structures through microseismic cluster and burst analyses complementing active seismic interpretation

At the Decatur carbon capture and storage site (IL, USA) CO 2 has been injected from 2011–2014 and from 2017 to present near the base of the Lower Mt. Simon Sandstone saline reservoir, resulting in microseismicity. Microseismicity is mainly located in the basement and distributed in distinct spatial clusters. The lack of significant impedance contrasts within the basement makes the interpretation of active-source seismic reflection data challenging, however, recent reprocessing allowed to resolve faults above and at the top of the basement. These faults generally do not coincide with the location of microseismic events and their continuation to the general depth of the seismic events cannot be assumed. This paper shows how the interpretation of the microseismicity can complement structural interpretations of active-source seismic reflection data. In particular, we analyze clusters and bursts (abrupt increases) of microseismicity, identify unresolved, smaller-scale weaknesses and extract statistical parameters. These parameters allow comparisons with the interpreted faults, and with fracture sets intercepted by boreholes. During injection at the Decatur site, the injection pressure was kept far below fracture pressure, nevertheless, seismic events were induced and spread far beyond the expected extent of the CO 2 plume. We argue that local stress transfers related to the CO 2 injection reactivated pre-existing fractures within the critically stressed basement. Finally, we conducted a slip tendency analysis for faults interpreted from active seismic, selected cluster, bursts and nodal planes from focal mechanisms to determine if the interpreted structures are optimally oriented with respect to the stress regime. Our results suggest that the orientation of fractures close to the injection well, generally shows slight deviations from the optimal orientation for slip. This might indicate either slight local deviations of the maximum horizontal stress azimuth from the average direction used in the analysis, or the lack of optimally oriented fractures at this location.

58 GEOSCIENCES↗
DOE OSTI2020

EXPERIMENTAL EVALUATION OF ULTRASONIC VELOCITIES AND ANISOTROPY IN THE TUSCALOOSA MARINE SHALE FORMATION

Tuscaloosa Marine Shale (TMS) formation is a clay- and organic-rich emerging shale play with a considerable amount of hydrocarbon resources. Despite the substantial potential, there have been only a few wells drilled and produced in the formation over the recent years. The analyzed TMS samples contain an average of 50 wt% total clay, 27 wt% quartz and 14 wt% calcite and the mineralogy varies considerably over the small intervals. The high amount of clay leads to pronounced anisotropy and the frequent changes in mineralogy result in the heterogeneity of the formation. We studied the compressional (V P ) and shear-wave (V S ) velocities to evaluate the degree of anisotropy and heterogeneity, which impact hydraulic fracture growth, borehole instabilities, and subsurface imaging. The ultrasonic measurements of P- and S-wave velocities from five TMS wells are the best fit to the linear relationship with R 2 = 0.84 in the least-squares criteria. We observed that TMS S-wave velocities are relatively lower when compared to the established velocity relationships. Most of the velocity data in bedding-normal direction lie outside constant V P /V S lines of 1.6–1.8, a region typical of most organic-rich shale plays. For all of the studied TMS samples, the S-wave velocity anisotropy exhibits higher values than P-wave velocity anisotropy. In the samples in which the composition is dominated by either calcite or quartz minerals, mineralogy controls the velocities and V P /V S ratios to a great extent. Additionally, the organic content and maturity account for the velocity behavior in the samples in which the mineralogical composition fails to do so. The results provide further insights into TMS Formation evaluation and contribute to a better understanding of the heterogeneity and anisotropy of the play.

Geochemistry & Geophysics↗
DOE OSTI2020

Joint Focal Mechanism Inversion Using Downhole and Surface Monitoring at the Decatur, Illinois, CO2 Injection Site

ABSTRACT The three-year CO2 injection period at the Illinois Basin - Decatur Project site (Decatur, Illinois, United States) produced a number of microseismic events distributed in very distinct spatiotemporal clusters with different orientations. Further characterization of the microseismicity encompasses the determination of the event source mechanisms. Initially, the microseismic monitoring network consisted solely of borehole sensors, but has been extended with surface sensors, thereby significantly improving the data coverage over the focal sphere. This article focuses on 23 events from the northernmost microseismic cluster (about 2 km from the injection point) and takes advantage of both, surface and downhole, recordings. The resulting strike-slip east–west-oriented focal planes are all consistent with the east–west orientation of the cluster in map view. The injection-related increase of pore pressure is far below the formation fracture pressure; however, small stress-field changes associated with the pore-pressure increase may reach as far as to the investigated cluster location. Monte Carlo modeling of the slip reactivation potential within this cluster showed that the observed maximum stress-field orientation of N068° is the optimum orientation for fault reactivation of the east–west-oriented cluster. Our results suggest that the east–west orientation of the investigated cluster is the main reason for its activation, even though the cluster is about 2 km away from the low-pressure injection point.

Geochemistry & Geophysics↗