Christopher E. Strickland

During fiscal years 2024 and 2025, the U.S. Department of Energy’s Hanford Field Office commissioned Pacific Northwest National Laboratory to conduct applied research aimed at reducing the cost, time, and uncertainty associated with in situ treatment of vadose zone contaminants at the Hanford Site. This report outlines the integration of three key research efforts into a meso-scale demonstration designed to advance field-scale solutions that aim to (1) optimize the delivery of chemical amendments to contaminated soils, (2) reduce uncertainty in amendment delivery performance assessment using advanced monitoring techniques, and (3) provide real-time insights into when and where amendment-induced precipitation reactions occur in the subsurface. To achieve these objectives, the tank-scale (~ 1 cubic meter) Geophysical Imaging of Flow and Transport (GIFT) system was developed. GIFT enables experimental testing of amendment delivery while incorporating automated multi-modal monitoring approaches, including pressure measurements, direct fluid sampling, and remote time-lapse geophysical imaging. The data generated from these monitoring techniques will serve as inputs for a generative artificial-intelligence-driven digital twin – a numerical simulation model designed to honor observed data while quantifying uncertainty in simulation accuracy. Using this simulator, researchers will refine an amendment injection strategy to maximize delivery efficiency within a low-permeability soil zone. Monitoring data will be interpreted through simulated outputs to enhance understanding of the injection process. The efficacy of this integrated approach will be evaluated through direct sampling at the conclusion of the experiment.

The GIFT (Geophysical Imaging of Flow and Transport) system will improve predictive capabilities as a platform for collecting multi-physical (MP) time-lapse datasets from well-controlled laboratory experiments. This will enable more accurate prediction of multi-phase fluid flow during delivery and transport of insitu remedial amendments as important components of the overall Hanford Site cleanup mission. Design and construction of the GIFT system, which consists of a meter-scale test cell and a network of custom gas circulation and analysis systems, was accomplished in fiscal year (FY) 2024, with plans to add liquid circulation to the system and use this new capability for multi-phase deep vadose experiments starting in FY25.

Fluid flow through fractured media is typically governed by the distribution of fracture apertures, which are in turn governed by stress. Consequently, understanding subsurface stress is critical for understanding and predicting subsurface fluid flow. Although laboratory-scale studies have established a sensitive relationship between effective stress and bulk electrical conductivity in crystalline rock, that relationship has not been extensively leveraged to monitor stress evolution at the field scale using electrical or electromagnetic geophysical monitoring approaches. In this paper we demonstrate the use time-lapse 3-dimensional (4D) electrical resistivity tomography to image perturbations in the stress field generated by pressurized borehole packers deployed during shear-stimulation attempts in a 1.25 km deep metamorphic crystalline rock formation.

A subsurface air injection at the Hanford Site’s Deep Vadose Zone Monitoring Test Bed was completed to realize a change in subsurface hydrologic conditions in accordance with a soil desiccation remedy. The injection mimicked a previous injection at the site that was performed in accordance with a vadose zone treatability study. Unlike the previous test which relied on electrical methods only, the change in hydrologic conditions during the recent test was also monitored using cross-hole seismic sensing methods to assess the ability of the seismic methods to evaluate changes in moisture conditions of desiccated sediment at the field scale. Data from in situ neutron probes indicates a reduction of soil moisture in the vicinity of the injection well due to the air injection. Similar changes were observed in the time-lapse electrical resistivity and seismic data, which indicates a loss of soil moisture over time. Tomographic inversions of the time-lapse geophysical data illustrate the 2D and 3D features of the soil moisture distribution over time. Time-lapse electrical resistivity tomography (ERT) results show a reduction in the electrical conductivity of the subsurface in the vicinity of the injection well, with most changes occurring within the screened interval. Similar patterns are observed inmore » the seismic tomography results, with both methods illustrating two lobe-shaped features of reduced soil moisture. Use of seismic and ERT technologies in tandem takes advantage of two complementary geophysical monitoring technologies, providing increased sensitivity to specific hydrologic conditions. The multiphysics approach, therefore, has the potential to improve the ability to estimate subsurface moisture conditions from sensor-based and remotely sensed geophysical data that will ultimately improve the ability of remediation contractors to evaluate remedy performance. « less

The goal of this Laboratory Directed Research and Development (LDRD) project was to develop a borehole seismic source and sensor array to enable real-time seismic imaging at scales and conditions relevant to the energy industry including both fossil-energy and geothermal. In FY21 and FY22, we designed, built, and tested both a prototype impulse source module for generating seismic energy and a sensing module for recording ground motions generated by the source module array. A pneumatically driven vibratory source was also designed. The source modules were fabricated with all high temperature components and the team has worked to incorporate the current RT-SEISMIC electronics design into a commercially available, high temperature silicon-on-insulator chip integrated circuit. Several issues were identified during fabrication and lab testing that led to redesign of several system components and subsequent retesting. The final round of testing showed that while metal/graphite-based seals worked quite well for static seals, they were unable to provide an adequate gas seal for dynamic, reciprocating part movements which necessitated a final redesign using Kalrez. This change will result in a continuous temperature rating of approximately 275 degrees C for the system. While a field test of the RT-SEISMIC system was targeted in FY22, due to the extended lab testing and redesign efforts, field testing was not achieved. As a result of this LDRD investment, several sponsors have expressed interest in RT-SEISMIC and we expect to continue towards a field demonstration of the full system in the future.

Many river corridor systems frequently experience rapid variations in river stage height, hydraulic head gradients, and residence times. The integrated hydrology and biogeochemistry of such systems is challenging to study, particularly in their associated hyporheic zones. Here we present an automated system to facilitate 4-dimensional study of dynamic hyporheic zones. It is based on combining real-timein-situandex-situmeasurements from sensor/sampling locations distributed in 3-dimensions. A novel dissolved oxygen (DO) sensor was integrated into the system during a small scale study. We measured several biogeochemical and hydrologic parameters at three subsurface depths in the riverbed of the Columbia River in Washington State, USA, a dynamic hydropeaked river corridor system. During the study, episodes of significant DO variations (~+/- 4 mg/l) were observed, with minor variation in other parameters (e.g., <~+/-0.15 mg/l NO₃). DO concentrations were related to hydraulic head gradients, showing both hysteretic and non-hysteretic relationships with abrupt (hours) transitions between the two types of relationships. The observed relationships provide a number of hypotheses related to the integrated hydrology and biogeochemistry of dynamic hyporheic zones. We suggest that preliminary high-frequency monitoring is advantageous in guiding the design of long term monitoring campaigns. The study also demonstrated the importance of measuring multiple parameters in parallel, where the DO sensor provided the key signal for identifying/detecting transient phenomena.

Compressional seismic (P‐wave) velocities were measured experimentally at different matric suctions and confining pressure states using a variety of sediments collected from the USDOE Hanford Site in southeastern Washington State. The P‐wave velocity was measured in variably saturated sediments using ultrasonic piezoelectric transducers, and matric suction was measured with a heat dissipation sensor. Desaturation of the samples occurred by means of evaporative drying. We propose a model to relate changes in P‐wave velocity to soil stress conditions using an effective stress formulation that includes capillary stress as well as adsorptive stress. We show that compressional wave velocities can be related to matric suction and confining stress using a model with five fitted parameters. A single set of parameters describing the relationship between effective stress and P‐wave velocity provide an excellent fit to data acquired on Hanford Formation samples. The results presented here suggest the utility of seismic methods for establishing low matric suction and flow conditions within variably saturated Hanford formation sediments, which could be applied at the field scale for evaluating the performance of environmental remediation efforts, such as soil desiccation and surface infiltration barriers.

Enhanced Geothermal Systems could provide a substantial contribution to the global energy demand if their implementation could overcome inherent challenges. Examples are insufficient created permeability, early thermal breakthrough, and unacceptable induced seismicity. Here we report on the seismic response of a meso-scale hydraulic fracturing experiment performed at 1.5 km depth at the Sanford Underground Research Facility. We have measured the seismic activity by utilizing a novel 100 kHz, continuous seismic monitoring system deployed in six 60 m-length monitoring boreholes surrounding the experimental domain in 3-D. The achieved location uncertainty was on the order of 1 m, and limited by the signal-to-noise ratio of detected events. These uncertainties were corroborated by detections of fracture intersections at the monitoring boreholes. Three intervals of the dedicated injection borehole were hydraulically stimulated by water injection at pressures up to 33 MPa and flow rates up to 5 L/min. We located 1933 seismic events during several injection periods. The recorded seismicity delineates a complex fracture network comprised of multi-strand hydraulic fractures and shear-reactivated, pre-existing planes of weakness that grew unilaterally from the point of initiation. We find that heterogeneity of stress dictates the outcome of hydraulic stimulations, even when relying on theoretically well-behaved hydraulic fractures. Once hydraulic fractures intersected boreholes, the boreholes acted as a pressure relief and fracture propagation ceased. In order to create an efficient sub-surface heat exchanger, production boreholes should not be drilled before the end of hydraulic stimulations.

River corridors, the spatial domains around rivers in which river water interacts with surrounding sediment and rock, are important components of watersheds. They comprise extremely complex ecosystems: heterogeneous at all spatial scales with strong temporal dynamics, coupled biological, geochemical, and hydrologic processes, and ubiquitous human impacts. We present several ways that our project, focused around the 75 km Hanford Reach of the Columbia River but with multiple connections to other systems, is addressing this challenge. These include 1) deployment of intensive, automated sensor networks supplemented by data from the Hanford Environmental Information System (HEIS) for hyporheic zone monitoring 2) data assimilation of these and other data into models using joint hydrologic and geophysical inversion, 3) integrating MASS2 model outputs and bathymetry data using machine learning to classify hydromorphologic features, 4) a community-based effort to develop broad understanding of organic carbon biogeochemistry and microbiomes in diverse river systems, and 5) use of multi-‘omics data to develop new biogeochemical reaction networks. These underpin the incorporation of process understanding and diverse data into high-resolution mechanistic models, and employment of those models to develop reduced-order models that can be applied at large scales while retaining the effects of local features and processes. In so doing we are contributing to reduction of uncertainties associated with major Earth system biogeochemical fluxes, thus improving predictions of environmental and human impacts on water quality and riverine ecosystems and supporting environmentally responsible management of linked energy-water systems.