Exploration geophysics
Geophysical methods for investigation of terrestrial crust and mineral prospecting
GEOPHYSICS↗
SEARCH · Engineering Papers
Results for “geophysical methods”
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.
Integration of Electromagnetic Geophysics Forward Simulation in Coupled Flow and Geomechanics for Monitoring a Gas Hydrate Deposit Located in the Ulleung Basin, East Sea, Korea
We investigate the feasibility of electromagnetic (EM) geophysics methods to detect the dissociation of gas hydrate specifically from a gas hydrate deposit located in the Ulleung Basin, East Sea, Korea via an integrated flow-geomechanics-EM geophysics simulation. To this end, coupled flow and geomechanics simulation is first performed with the multiple porosity model employed, where a mixed formulation with the finite volume (FV) and finite element (FE) methods are taken for the flow and geomechanics, respectively. From the saturation and porosity fields obtained from the coupled flow and geomechanics, the electrical conductivity model is established for the EM simulation. Solving the partial differential equation of electrical diffusion which is linearized using the 3D finite element method (FEM), the EM fields are then computed. For numerical experiments, particularly two approaches in the configuration for the EM methods are compared in this contribution: the surface-to-surface and the surface-to-borehole methods. When the surface-to-surface EM method is employed, the EM is found to be less sensitive, implying low detectability. Especially for the short term of production, the low detectability is attributed to the similarity of electrical resistivity between the dissociated gas (CH4) and hydrate as well as the specific dissociation pattern within the intercalated composites of the field. On the other hand, when the surface-to-borehole EM method is employed, its sensitivity to capture the produced gas flow is improved, confirming its detectability in monitoring gas flow. Hence, the EM geophysics simulation integrated with coupled flow and geomechanics can be a potential tool for monitoring gas hydrate deposits.
depressurization↗
Using Controlled-Source Electromagnetic (CSEM) for CO 2 Storage Monitoring in the North Dakota CarbonSAFE Project
Conference presentation for International Meeting for Applied Geophysics & Energy (IMAGE), September 26 – October 1, 2021, Denver, CO. Geophysical methods are key for characterizing the geologic formations to store CO 2 and monitor the injected CO 2 over time to ensure containment. In the integrated multimeasurement geophysical approach considered for this project, it is expected that the controlled-source electromagnetic (CSEM) method is a strong contributor to mapping the CO 2 movement. A feasibility study of the CSEM method, including 1D and 3D modeling and a field noise test, was conducted to determine its effectiveness in monitoring CO 2 in the Broom Creek and Deadwood Formations. The study results demonstrate that the CSEM method can be used for CO 2 storage monitoring in the study area.
20 FOSSIL-FUELED POWER PLANTS↗
Sensitivity of geophysical techniques for monitoring secondary CO 2 storage plumes
For geologic carbon storage, the ability to detect secondary CO 2 plumes—defined as those CO 2 plumes accumulating outside the intended storage reservoir—is fundamental to preventing unexpected CO 2 migration into groundwater resources and for risk and liability management. Understanding the sensitivity of various geophysical methods to secondary plumes is crucial for designing cost-effective monitoring schemes. We use several modeling scenarios to demonstrate the process of assessing sensitivities and detection thresholds of three primary geophysical techniques—surface seismic, borehole-to-surface electromagnetic (EM), and surface and borehole gravity—for early detection of secondary CO 2 plumes in the post-injection phase. While seismic reflection methods are often considered in monitoring strategies to track the evolution of CO 2 plumes, they are also the most expensive. Due to cost considerations, especially for long-term post-injection monitoring, other techniques complement seismic monitoring when designing an adaptive monitoring network. Borehole-to-surface EM or surface gravity surveys are feasible for time-lapse monitoring of deep secondary CO 2 plumes. Furthermore, these surveys could be carried at intervals defined by site-specific conditions. If time-lapse EM and/or gravity surveys detect any signal responses beyond the expected change, it would trigger a need for the higher resolution seismic survey.
58 GEOSCIENCES↗
Capturing soil moisture and salinity changes in flooded coastal forests using electrical resistivity and induced polarization
Geophysical methods provide high-resolution spatial measurements of physical quantities sensitive to changes in soil moisture and salinity and can be used to monitor soil hydrological responses to flooding. However, extracting quantitative hydrological information from geophysical data remains challenging. In this study, we extended existing petrophysical models to estimate soil moisture and salinity from electrical measurements to address this challenge. We manipulated two hydrologically isolated 2000 m 2 experimental plots by simultaneously inundating them with 265 m 3 of either freshwater or estuarine water. Repeated electrical resistivity and induced polarization measurements were used to image the water and solute infiltration along two transects that are 100 and 42 m in length. Petrophysical models derived from laboratory multi-salinity electrical measurements were used to estimate changes in soil moisture and fluid salinity from field measurements of real and imaginary conductivity during the ecosystem-scale flooding experiment. The real conductivity increased by ∼100% in the freshwater plot and ∼570% in the saltwater plot. The change in imaginary conductivity in the freshwater plot was <1 mS/m, whereas that of the estuarine water plot was ∼5 mS/m. Real conductivity shows dependence on soil moisture content with a coefficient of determination (R 2 ) > 0.7, while the imaginary conductivity shows a dependence on soil salinity with R 2 > 0.6. The results validate the use of electrical resistivity for estimating changes in soil moisture content in response to flooding. Combining electrical resistivity imaging with induced polarization measurements provides the possibility to account for changes in pore fluid conductivity.
Adebayo, Moses B. [Univ. of Toledo, OH (United Sta↗
Material studies related to lunar surface exploration. Volume 2 - Application of geophysical and geotechnical methods to lunar sites exploration Final report, 6 Mar. 1967 - 30 Jun. 1968
Geophysical and geotechnical methods for lunar landing site soil property studies
Carmichael, I. S. E.↗
Advanced High-Performance Computational Modeling of the Seismic Response of High-Hazard and/or Nuclear Facilities and Critical Infrastructure at the NNSS
New methods for predicting the amplitude and variability of ground shaking from earthquakes (and explosions) are needed for seismic hazard analysis for buildings, nuclear power plants, and critical infrastructure at the NNSS. We are comparing existing 1-D and new 3-D geophysical methods for estimating the shear-wave velocity structure in the upper 30 meters of the ground surface (Vs30), which plays a major role in ground motion amplification and seismic response of buildings. We evaluate the performance of these methodologies at the U1a Complex at the NNSS and develop simple 1-D and high-resolution 3-D Vs30 models. We then emplace these high-resolution models into a background seismic velocity model. We will collaborate with Lawrence Livermore National Laboratory (LLNL) to conduct numerical modeling of the ground shaking at the NNSS using their high-performance computing technology and state-of-the-art ground motion simulation methodology. The primary work that was completed in FY 2019 was to acquire the seismic systems and familiarize staff at the NNSS with their use. We also worked on developing a collection plan with the Device Assembly Facility (DAF) at the NNSS, but due to time constraints and other ongoing projects at the DAF, we had to use U1a as a backup. We were able to coordinate the seismic survey, and we will complete the collection of seismic data in FY 2020. Additionally, during FY 2019, we completed the geologic framework model (GFM) for the U1a Complex and modeled the Yucca fault. LLNL worked with us through FY 2019 to prepare the data files for their modeling software and tested the software for reliability. In FY 2020 we will develop the end-to-end capability so that any facility could easily be modeled and the expected shaking from a local earthquake understood. The work in FY 2020 will include building fault models from the GFM and finalizing the velocity model analysis. The final simulations will be run for multiple rupture models, and final assessments will demonstrate the seismic hazard at the U1a Complex.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Geophysical Signatures of Crack Network Coalescence in Rocks at Multiple Length Scales
The main goal of the research project was to identify the geophysical signatures of fracture growth in natural rocks by utilizing novel geophysical techniques. The research objectives were to (a) investigate the potential for geophysical methods to determine when cracks initiate, the types and locations of propagated cracks, and the coalescence of networks of cracks in natural rocks at multiple scales, (b) determine how damage at the microscale evolved into damage at the macroscale and then link the microscopic and macroscopic observations, (c) quantify crack coalescence in rocks under realistic stress conditions using coupled mechanical-geophysical-optical visualization, and (d) identify the precursors in geophysical signals to crack coalescence. The following research thrusts were explored to achieve the research objectives: (1) uniaxial compression testing of rock specimens with and without a set of pre‐existing flaws and (2) triaxial compression testing of natural rock specimens. These thrusts allowed for exploring fracturing in rocks under realistic in situ environments and at multiple scales. This project provided educational opportunities for nine graduate and undergraduate students and resulted in 27 peer-reviewed publications. This first research thrust focused on investigating the micromechanics of fractures in rocks through uniaxial compression testing combined with advanced geophysical and imaging techniques, specifically acoustic emission (AE) monitoring, ultrasonic imaging, and 2-dimensional Digital Image Correlation (2D-DIC). By examining damage processes under time-independent and time-dependent loading conditions, insights into damage localization, crack initiation, and fracturing mechanisms were gained. It was observed that the AE signals and the strain-based measurements directly reflect the state of damage in the rock specimen and could be used to identify the cracking levels, such as the crack initiation (CI) and crack damage (CD), and the mode of deformation. A novel calibration apparatus was developed to enhance the accuracy of AE sensors, allowing for the estimation of key parameters such as magnitude, source dimension, stress drop, and radiated seismic energy associated with the fractures. The findings highlighted significant variations in the temporal evolution of AE source parameters during the primary, secondary, and tertiary stages of creep, identifying tensile cracking as the primary deformation mode. The second research thrust focused on enhancing the understanding of fracturing processes in natural rocks through triaxial compression testing, real-time AE monitoring, and ultrasonic monitoring. We investigated the impact of various factors such as fracture propagation regimes, injection parameters, rock types, and pre-existing conditions on the hydraulic fracture (HF) behavior using scaled true-triaxially loaded specimens of Barre granite and Lyons sandstone. Custom sensor housing facilitated concurrent active and passive monitoring to analyze hydro-mechanical responses and microseismicity associated with different HF scenarios. A coupled investigation of passive microseismicity and active signal attributes permitted a detailed comprehension of the various HF processes (aseismic deformation, fracture initiation and propagation, fluid permeation, and leak-off) and their dependence on the specific rock type. The findings of this research demonstrated the effectiveness of AE monitoring techniques in providing valuable insights into the impact of various factors on the behavior and dynamics of HF processes. The advancements in monitoring techniques, offering a more thorough and precise approach, represent a significant step towards optimizing HF practices and ensuring sustainable resource extraction.
58 GEOSCIENCES↗
Field Test Bed for Vadose-Zone Monitoring Approaches - 20404
In the Central Plateau at the U.S. Department of Energy Hanford Site, a large inventory of contaminants resides in unsaturated sediments within the approximately 100-meter-thick vadose zone, posing a potential continuing risk to groundwater. Vadose zone remedies used to address these contaminants will require performance monitoring to provide feedback during implementation and for long-term verification that remedial action objectives have been met. Passive approaches may also need long-term monitoring to demonstrate that the flux of contaminants from the vadose zone to the groundwater are below thresholds established to meet groundwater protection goals. Collection of physical (e.g., groundwater or sediment) samples is a common method for identifying contaminant concentration distributions, but this approach is limited by the number of locations and the frequency with which data can be collected. In situ vadose zone measurements have evolved over the past few years to include key measurements of water content, soil water pressure, temperature, and chemical concentration. However, the current generation of sensors is designed for relatively short-term use in near-surface soils or sediments. Geophysical methods have been evolving but are also limited in that they have not been designed for the specific long-term vadose zone monitoring needs at the Hanford Site. Overall, monitoring under unsaturated conditions can be difficult due to the need to install and maintain instrumentation over a large area and depth and the need to identify preferential flow pathways due to geologic and chemical heterogeneities over long periods. Thus, a vadose zone monitoring test bed was initiated to address these challenges and identify cost-effective approaches for implementation and postclosure monitoring of the deep vadose zone. The monitoring test bed is expected to provide valuable field-scale information for the design of vadose zone monitoring systems. (authors)
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NETL Plastic Pipes Project (Final Report)
Plastic or composite pipelines have been the bane of the utility locating industry because they are neither conductive nor magnetic which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipelines where resistive cover allows for adequate penetration of radar energy. However, GPR has limited applicability in areas where the soil cover is conductive due to significant clay and/or salt content. This study examines complementary near-surface geophysical methods that are potentially useful for locating buried plastic/composite pipelines, either singly or in combination. Specifically, this modeling study used computational numerical methods to forward model the response of GPR, resistivity, seismic, gravity gradiometry, and photoacoustic/thermoacoustic imaging methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft. to 4 ft.; (3) various degrees in contrast in physical properties (i.e., electrical permittivity, elasticity, resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar and seismic source frequencies, electrode spacing). Numerical modeling performed herein reconfirmed that GPR is the preferred method for detecting/locating plastic pipelines. A caveat for GPR detection is that the material covering the plastic pipe (trench fill material and adjacent soil) must be sufficiently resistive to allow the two-way propagation to the required depth of investigation and back to the surface. GPR was the only method modeled in this study that can be used to directly detect plastic pipelines of 2-in.-diameter and larger when buried 3-ft-deep. GPR data processing and imaging also can determine pipe depth, pipe diameter, trench dimensions, and moisture conditions. Seismic modeling results suggest that direct detection of a 12-in.-diameter plastic pipe at 3-ft.-depth may be possible under favorable conditions; however, the associated signature would be weak (e.g., surface- to S-wave, backscattered surface-waves, and/or forward scattered surface-waves to S-wave). Direct pipe detection under field conditions with noise and strong lateral geologic heterogeneity is doubtful. Numerical modeling also suggests that plastic pipelines can be indirectly located by detecting the trench in which they are buried. GPR, direct current (DC) resistivity, and seismic methods have the potential to locate the pipeline trench if there is sufficient contrast between the trench-wall and trench-fill materials for the physical property being measured by each method (i.e., electrical permittivity for GPR; resistivity for DC resistivity; or density, compressional velocity, or shear velocity for seismic). Modeling also indicated that currently available (commercial) gravity gradiometers would be unable to directly detect/locate plastic pipelines ≤ 8-in.-diameter when buried 3-ft.-deep given the typical instrument noise floor for field surveying as well as the expected density variations due to geologic heterogeneity. The numerical modeling performed in this project did not identify a universal geophysical technology that can locate buried plastic pipelines in all parts of the United States (although GPR is suggested for all areas with resistive cover). However, the project results suggest that a towed land streamer simultaneously acquiring multiple geophysical data types including multi-offset GPR, multi-channel DC resistivity, seismic geophone- and/or distributed acoustic sensing (DAS), and potentially photoacoustic/thermoacoustic data would be an appropriate platform for locating buried plastic pipeline. Moreover, the complementary multiphysics data acquired by a towed land streamer would permit the use of joint and/or cooperative inversion frameworks for a more rigorous and consistent data interpretation.
42 ENGINEERING↗
A Flux Detection Probe to Quantify Dynamic Groundwater-Surface Water Exchange in the Hyporheic Zone
A new probe was designed to quantify groundwater-surface water exchange in the hyporheic zone under dynamic stage condition. Current methods focus on either vertical pore water velocity or Darcy flux measurements. Both parameters must be understood to evaluate residence time and mass flux of constituents. Furthermore, most instruments are not well suited for monitoring instantaneous velocity or flux under dynamic exchange conditions. For this reason, the flux detection probe (FDP) was designed that employs electrogeophysical measurements to estimate in situ sediment porosity, which can be used to convert pore water velocity to Darcy flux. Dynamic pore water velocity is obtained by monitoring fluid conductivity and temperature along the FDP probe. Pressure sensors deployed at the top and bottom of the probe provide the additional information necessary to estimate vertical permeability. This study focuses on the use of a geophysical method to estimate pore water velocity, porosity, and permeability within a controlled soil column where simulated river water displaces simulated groundwater. The difference between probe derived and theoretical pore water velocity using natural tracers such as electrical conductivity and temperature was -4.9 and 3.9% for downward flow and 1.1 and 12.8% for upward flow, respectively. The difference in porosity calculated from mass and volume packed in the soil column and probe measure porosity ranged between -3.2% and 1.5%. Also, the calculated hydraulic conductivity differed from probe derived values by -8.9%.
54 ENVIRONMENTAL SCIENCES↗
Development of an ERT‐Based Framework for Bentonite Buffers Monitoring From Laboratory Tests: 1. Characterizing Thermal–Hydrological–Mechanical Processes
Abstract Bentonite clay is widely used in engineered barrier systems for the permanent disposal of high‐level radioactive waste due to its low permeability, high swelling capacity, and thermal stability. However, the complex thermal‐hydrological‐mechanical (THM) processes induced by heating from decaying radioactive waste and hydration from surrounding rock can lead to heterogeneous changes that are difficult to measure and predict. This study develops an Electrical Resistivity Tomography (ERT)‐based framework for monitoring THM processes, progressing from sample‐scale to bench‐scale tests, to inform field‐scale applications. Sample‐scale tests analyzed small bentonite samples under controlled variations in water content, temperature, and porosity to establish fundamental resistivity relationships. Bench‐scale tests involved larger bentonite columns subjected to heating (up to 200°C) and hydration under controlled pressure, simulating repository conditions. ERT measurements, complemented by X‐ray CT imaging, temperature monitoring, and tracing sensors, revealed coupled THM processes, such as hydration‐induced compression, swelling, and thermal gradients, leading to complex resistivity patterns. The results demonstrate the potential of ERT for capturing THM‐induced resistivity changes, though challenges remain in upscaling and quantitative analysis. This study evaluates laboratory test capabilities and proposes future improvements for understanding THM‐induced resistivity responses. A conceptual framework for ERT implementation in field‐scale monitoring is presented, synthesizing findings from both scales and exploring how ERT data can inform long‐term modeling and reduce prediction uncertainties. Overall, this ERT‐based framework offers a robust method for monitoring bentonite buffers, aiding in early issue detection and supporting the safe long‐term disposal of radioactive waste in geological repositories, while highlighting the need for future development. Plain Language Summary Bentonite clay is crucial in engineered barrier systems (EBS) for containing high‐level radioactive waste due to its ability to absorb water, swell, seal and remain stable under high temperatures. When bentonite absorbs water and heats up from radioactive decay, it experiences complex changes in its physical and mechanical properties. Understanding these changes is important for ensuring the long‐term safety and effectiveness of EBS. This study used Electrical Resistivity Tomography (ERT), a non‐invasive method that measures electrical conductivity to monitor these changes during laboratory experiments. The ERT data revealed significant variations in resistivity corresponding to changes in water content, temperature, and density, providing detailed spatial and temporal insights into the behavior of bentonite. These findings enhance our ability to predict the long‐term performance of bentonite barriers, ensuring the safe containment of radioactive waste. By improving our understanding of bentonite's behavior, this research supports the development of more reliable and effective barrier systems for radioactive waste disposal, protecting the environment and public health. Key Points ERT monitoring was employed to capture resistivity changes in bentonite during controlled heating and hydration experiments, providing insights into THM processes ERT data reveal significant resistivity changes correlated with water content, temperature, and mechanical effects, enhancing the understanding of THM dynamics in bentonite This study explores the potential of the framework for application in field‐scale EBS monitoring, emphasizing the need for integrating additional geophysical methods for comprehensive subsurface imaging
Chen, Hang↗
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↗
Exploring the Use of Non‐Invasive Drone‐Based Ground‐Penetrating Radar (GPR) to Characterize Biogenic Gas Dynamics in Subtropical Peat Soils
Peat soils are a critical component of the global carbon cycle as natural producers of biogenic greenhouse gases (e.g., methane and carbon dioxide) that accumulate within the soil and are released to the atmosphere. Previous studies have showed the ability of ground-based minimally-invasive geophysical methods such as ground-penetrating radar (GPR) to characterize carbon dynamics in peat soils. However, ground-based GPR is limited by scale of measurement and soil disturbance potentially altering gas releases during deployment. Here, we explore the potential of drone-based GPR for identification of hot spots and hot moments of gas accumulation and release in subtropical soils. Here, we collected drone-based GPR data sets across two grids (∼17,500 m 2 ) in the Everglades during January (dry season), September, and November (wet season) of 2023 to characterize peat thickness and seasonal variability of gas content. Results show that drone-based GPR is effective and efficient for: (a) capturing the temporal variation of in situ biogenic gas content in peat soils with changes between 1% and 25 % volumetric gas content over repeatable grids; (b) inferring a total peat thickness between 0.8 and 1.2 m; and (c) estimating flux releases of 63 and 135 mg CH 4 m −2 day −1 for specific locations and periods that are strikingly consistent with our coincident gas trap measurements. This work also indicates that (a) spatial distribution of gas content in the Everglades is strongly controlled by landscape morphology such as ridges and sloughs and (b) the temporal variation of gas content is seasonal with increased gas production during the wet season.
54 ENVIRONMENTAL SCIENCES↗
Geomechanical and Hydrogeological Evaluation of a Shallow Hydraulic Fracture at the Devine Fracture Pilot Site, Medina County, Texas
UT-Austin’s Devine Fracture Pilot Site, 50 miles southwest of San Antonio, Texas, has been targeted for a comprehensive, multidisciplinary development of fracture diagnostic techniques that are cross-validated by ground-truth data acquisition near a recently created, 175-ft-deep, horizontal hydraulic fracture (Ahmadian et al. 2018 Demonstration of proof of concept of electromagnetic geophysical methods for high resolution illumination of induced fracture networks. In Proceedings of the SPE Hydraulic Fracturing Technology Conference and Exhibition, The Woodlands, Texas, USA, 23–25 January 2018. SPE-189858-MS.). To evaluate the fracture diagnostic methods at this site, we conducted injection tests with a predefined volumetric flow-rate profile, resembling a diagnostic fracture injection test on September 2020. Subsequently, we developed hydrogeological and geomechanical models based on flow-rate and bottomhole-pressure measurements. History-matching efforts using a simplified layer-cake hydrogeological model resulted in the field-scale formation permeability of 9.87 × 10 –15 m 2 (10 mD) and Darcy-scale fracture permeability. The analysis of the bottomhole pressure and injection-rate history showed that (1) the newly created horizontal fracture was closed adjacent to the injection well pre-injection and (2) the initial pump-pressure increase at a nominal volumetric injection rate led to near-well fracture reopening, fluid conductivity increase, and abrupt injection-rate increase. To overcome hydrogeological-model limitations of predicting fracture reopening throughout injection, we extended the modeling to a finite-element, poroelastic analysis of horizontal-fracture growth using a cohesive-zone model. Using this fracture-reopening model, we improved the history match of the transient-pressure response during the experiment by adjusting the hydromechanical properties. Furthermore, post-injection pressure transient analyses helped reduce uncertainty in the overburden-stress gradient, and the initial hydraulic-fracturing simulation verified the plausibility of achieving the surveyed propped fracture area.
02 PETROLEUM↗
Three-dimensional fracture continuum characterization aided by surface time-domain electromagnetics and hydrogeophysical joint inversion—proof-of-concept
Efficient and safe production of hydraulically fractured reservoirs benefits from the prediction of their geometrical attributes. Geophysical methods have the potential to provide data that are sensitive to fracture geometries, alleviating the typically sparse nature of in situ reservoir observations. Moreover, surface-based methods can be logistically and economically attractive since they avoid operational interference with the injection well infrastructure. Here, we investigate the potential of the surface-based time-domain electromagnetic (EM) method. EM methods can play an important role owing to their sensitivity to injection-induced fluid property changes. Two other advantageous factors are the EM signal-enhancing effect of vertical steel-cased wells and the fact that injected proppants can be enhanced to produce a stronger electrical conductivity contrast with the reservoir’s connate fluid. Nevertheless, an optimal fracture characterization will no doubt require the integration of EM and reservoir injection and production data. We hence carry out our investigations within a hydrogeophysical parameter estimation framework where EM data and injection flow rates are combined in a fully coupled way. Given the interdisciplinary nature of coupled hydrogeophysical inverse modeling, we dedicate one section to laying out key aspects in a didactic manner.
58 GEOSCIENCES↗
Integration of Large-Scale Electrical Imaging into Geological Framework Development and Refinement
Geologic framework models (GFMs) are critical to the construction of reliable simulation models of groundwater flow and contaminant transport. To support GFM development, direct information (e.g., core samples, fluid samples, hydraulic testing) tends to be sparse and separated by large distances relative to the spatial scales of aquifer heterogeneity. There are additional challenges associated with highly contaminated legacy waste sites, where drilling is particularly costly, and invasive sampling requires specialized handling and disposal of hazardous materials. At these sites in particular, non-invasive geophysical imaging can play an important role in filling spatial gaps between boreholes and reducing characterization costs by optimizing and minimizing the number of necessary boreholes. Here this paper presents a case study demonstrating the use of large-scale (> 30 km 2 ) electrical mapping to identify hydrostratigraphy and potential paleochannels at the Hanford Site, located in Washington State, USA. In two field campaigns, over 36 line-kilometers of electrical resistivity tomography (ERT) data were collected along 14 transects. ERT surveys were sited and performed to image critical aspects (e.g., paleochannels, stratigraphic contacts) of the subsurface, demonstrating a general workflow for integrating ERT with GFM development. Inconsistencies between the GFM and ERT were catalogued to provide a basis for future site characterization using complementary geophysical methods and (or) direct sampling.
58 GEOSCIENCES↗
Deep learning to estimate permeability using geophysical data
Time-lapse electrical resistivity tomography (ERT) is a popular geophysical method to estimate three-dimensional (3D) permeability fields from electrical potential difference measurements. Traditional inversion and data assimilation methods are used to ingest this ERT data into hydrogeophysical models to estimate permeability. Due to ill-posedness and the curse of dimensionality, existing inversion strategies provide poor estimates and low resolution of the 3D permeability field. Recent advances in deep learning provide us with powerful algorithms to overcome this challenge. This paper presents a deep learning (DL) framework to estimate the 3D subsurface permeability from time-lapse ERT data. To test the feasibility of the proposed framework, we train DL-enabled inverse models on simulation data. Each measurement in both synthetic and field data is standardized by removing the mean and scaling the time-series to unit variance. This pre-processing step is necessary to bring simulation data closer to field observations. Subsurface process models based on hydrogeophysics are used to generate this synthetic data. Training performed on limited simulation data resulted in the DL model over-fitting. An advanced data augmentation based on mixup is implemented to generate additional training samples to overcome this issue. This mixup technique creates weakly labeled (low-fidelity) samples from strongly labeled (high-fidelity) data. The weakly labeled training data is then used to develop DL-enabled inverse models and reduce over-fitting. As both time-lapse ERT (1133048 features/realization) and 3D permeability (585453 features/realization) data samples are from a high-dimensional space, principal component analysis (PCA) is employed to reduce dimensionality. Encoded ERT and encoded permeability are generated using the trained PCA estimators. A deep neural network is then trained to map the encoded ERT to encoded permeability. This mixup training and unsupervised learning allowed us to build a fast and reasonably accurate DL-based inverse model under limited simulation data. Results show that proposed weak supervised learning can capture salient spatial features in the 3D permeability field. Quantitatively, the average mean squared error (in terms of the natural log) on the strongly labeled training, validation, and test datasets is less than 0.5. The R 2 -score (global metric) is greater than 0.75, and the percent error in each cell (local metric) is less than 10%. Finally, an added benefit in terms of computational cost is that the proposed DL-based inverse model is at least O(10 4 ) times faster than running a forward model once it is trained. Data generation, DL model training, and hyperparameter tuning to identify optimal neural network architectures utilized high-performance computing resources while the DL inference is performed on a standard laptop. Approximately, O(10 5 ) processor hours are used for generating data and DL tuning and training. We acknowledge that the data generation and DL model development are expensive. But once a DL model is trained, it can be re-used for inversion rapidly for the given system, with set physics and domain. Note that traditional inversion may require multiple forward model simulations (e.g., in the order of 10 to 1000), which are very expensive. This computational savings ≈ O(10 5 ) – O(10 7 )) makes the proposed DL-based inverse model attractive for subsurface imaging and real-time ERT monitoring applications due to fast and yet reasonably accurate estimations of permeability field.
58 GEOSCIENCES↗