SEARCH · Search NASA
Results for “Subsurface characterization”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Estimating biofuel contaminant concentration from 4D ERT with mixing models
Here, we present the results of a lab-scaled feasibility study to assess the performance of electrical resistivity tomography for detection, characterization, and monitoring of fuel grade ethanol releases to the subsurface. Further, we attempt to determine the concentration distribution of the ethanol from the electrical resistivity tomography data using mixing-models. Ethanol is a renewable fuel source as well as an oxygenate fuel additive currently used to replace the known carcinogen methyl tert-butyl ether; however, ethanol is preferentially biodegraded and a cosolvent. When introduced to areas previously impacted by nonethanol-based fuels, it will facilitate the persistence of carcinogenic fuel compounds like benzene and ethylbenzene, as well as remobilize them to the ground water. These compounds would otherwise be retained in the soil column undergoing active or passive remediation processes such as soil vapor extraction or natural attenuation. Here, we introduce ethanol to a saturated Ottawa sand in a tank instrumented for four-dimensional geoelectrical measurements. Forward model results suggest pure phase ethanol released into a water saturated silica sand should present a detectable target for electrical resistivity tomography relative to a saturated silica sand only. We observe the introduction of ethanol to the closed hydraulic system and subsequent migration over the duration of the experiment. One-dimensional and three–dimensional temporal data are assessed for the detection, characterization, and monitoring of the ethanol release. Results suggest one-dimensional geoelectrical measurements may be useful for monitoring a release, while three-dimensional geoelectrical field imaging would be useful to characterize, monitor, and design effective remediation approaches for an ethanol release, assuming field conditions do not preclude the application of geoelectrical methods. We then attempt to use predictive mixing models to calculate the distribution of ethanol concentration within the measurement domain. For this study we examine four different models: a nested parallel mixing model, a nested cubic mixing model, the complex refractive index model (CRIM), and the Lichtenecker-Rother (L-R) model. The L-R model, modified to include an electrical formation factor geometry term, provided the best agreement with expected EtOH concentrations.
Spectral Induced Polarization-Biogeochemical Relationships for Remediation Amendment Monitoring
Geophysical tools such as electrical resistivity (ER) can indirectly monitor subsurface changes in response to remedial injections. These methods exhibit relatively low spatial resolution compared to sediment core characterization but are advantageous due to the ability to collect measurements non-intrusively over time across large volumes of the subsurface. Moreover, along with confirmatory groundwater or core sampling, geophysical methods can be used during active biogeochemical remedies to monitor short-term contaminant transformations and mobility, as well as part of an overall strategy for long-term monitoring of subsurface contaminated sites. The use of alternating current spectral induced polarization (SIP) provides significantly more information than conventional geophysical methods like direct current ER. SIP allows for monitoring changes in both solution and surface conductivity by separation of real and imaginary conductivity, respectively, as well as surface capacitance. In principle, SIP can measure indicators of remedy progression such as precipitation reactions that sequester contaminants, potentially providing a better indication of amendment delivery and reactivity as compared to conventional ER methods. However, multiple processes and material properties have overlapping (interacting) electrical responses across a range of frequencies. Hence, the purpose of this scoping study was to evaluate if SIP can be used to monitor (a) amendment delivery and (b) precipitation and reactivity of amendments under consideration at Hanford. SIP measurements were collected in fully saturated columns packed with sand and Hanford formation sediments containing (a) amendments that were highly conductive with significant capacitance (zero valent iron – ZVI, sulfur modified iron – SMI) and (b) amendments that exhibited low electrical conductivity with a small capacitance (calcite, apatite, bismuth). The sand was a quartz material with homogenous particle size that exhibited a relatively low surface conductivity. It was used as a control for comparison with the sediments from the Hanford Site, which have a greater surface conductivity because of their complex mineralogy and heterogeneous particle size distribution and may have complex interactions with amendments. The amendment mass fraction was varied to represent the different stages and subsurface locations associated with the delivery of a remedy. The primary objective was to identify the solution and solid surface changes associated with the delivery amendments and their secondary reactions within the subsurface. The figure below summarizes results for the amendments tested through this project. The SIP phase shift or imaginary conductivity change for the high conductivity amendments was more than 10 times that of the low conductivity amendments, highlighting the relative ease of detection of ZVI and SMI independent of the background signal from sand or Hanford sediments. The ZVI phase shift and imaginary conductivity changes occur primarily at high frequency (> 100 Hz) whereas SMI changes were at low frequency (0.01 to 10 Hz). SIP signals of SMI also increased over time and the maximum shifted to lower frequencies. The low conductivity amendments (calcite, abiotic and biotic apatite, bismuth subnitrate) exhibited relatively small phase and imaginary conductivity changes when added to sediments. The changes were above minimum detection limits (0.5 mrad for phase shift, 0.03 µS/cm for imaginary conductivity) for the highest concentration except for the commercial bismuth material. The lowest amendment concentration that can be detected is likely sediment specific, as minerals in sediments (clays, magnetite, Fe-oxides) have some capacitance and, therefore, exhibit a variable phase shift.
Characterization of natural organic matter in low-carbon sediments: Extraction and analytical approaches
Sediment-associated natural organic matter (NOM) is a complex assemblage of organic molecules with a wide range of sizes, functional groups, and structures, which is intricately associated with mineral particles. Organic carbon (OC) concentrations in subsurface sediments are typically 10 to 200 times lower than in surface soils, posing a distinct challenge for characterization. A range of chemical extractions were evaluated for extraction of NOM, and a NOM extraction scheme was developed using a combination of sequential extraction with water (MQ) and sodium pyrophosphate at pH 10 (PP), and purification by dialysis and solid phase extraction in order to isolate different fractions of sediment-associated NOM. Analysis of these different NOM fractions was then carried out by Fourier transform infrared spectroscopy (FTIR), ultraviolet-visible (UV-Vis) spectroscopy, and electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (ESI-FTICR-MS). The OC extraction efficiency of the tested extractions followed the order PP>NaOH>HCl=hydroxylamine hydrochloride>dithionite. Greater concentrations of OC in PP vs NaOH extracts suggest that metal complexation and/or ligand exchange plays an important role in OC stabilization. Characterization of different pools of extraction NOM by FITR shows that the water soluble fraction has a higher fraction of aliphatic and carboxylic groups, while the PP fractions have higher C=C groups. This trend from aliphatic to more aromatic is also supported by the UV-Vis and ESI-FTICR-MS data.
Defining Modeling and Analysis Complexity Needs for Remediation - 18425
Many subsurface remedy decisions rely on model predictions or other types of data analyses to support site characterization and remedy selection, design, implementation, performance assessment, and closure. Thus, determining the appropriate level of complexity for modeling and data analysis is important to facilitate subsurface remedy decisions. Multiple factors are important in selecting the appropriate level of complexity in modeling and data analysis. For instance, the type and number of inputs required typically increases with increased complexity of analysis, and the ability to obtain or verify those inputs needs to be considered. The approach needs to establish the relationship between parameters or processes and contaminant transport and the detail of the description needed to meet the analysis objectives. Assessing the modeling complexity associated with the appropriate level of detail can be facilitated by considering complexity of the 1) modeling approach, 2) description of process, and 3) description of heterogeneity. As demonstrated by case studies, using a framework to guide determination of modeling and analysis complexity facilitates selection of the appropriate approach for the site situation. (authors)
Machine-Learning of Nonlocal Kernels for Anomalous Subsurface Transport from Breakthrough Curves
Anomalous behavior is ubiquitous in subsurface solute transport due to the presence of high degrees of heterogeneity at different scales in the media. Although fractional models have been extensively used to describe the anomalous transport in various subsurface applications, their application is hindered by computational challenges. Simpler nonlocal models characterized by integrable kernels and finite interaction length represent a computationally feasible alternative to fractional models; yet, the informed choice of their kernel functions still remains an open problem. We propose a general data-driven framework for the discovery of optimal kernels on the basis of very small and sparse data sets in the context of anomalous subsurface transport. Using spatially sparse breakthrough curves recovered from fine-scale particle-density simulations, we learn the best coarse-scale nonlocal model using a nonlocal operator regression technique. Predictions of the breakthrough curves obtained using the optimal nonlocal model show good agreement with fine-scale simulation results even at locations and time intervals different from the ones used to train the kernel, confirming the excellent generalization properties of the proposed algorithm. A comparison with trained classical models and with black-box deep neural networks confirms the superiority of the predictive capability of the proposed model.
Groundwater characterization and monitoring at a complex industrial waste site using electrical resistivity imaging
A contaminated industrial waste site in Washington State (USA) containing buried, metallic-waste storage tanks, pipes, and wells, was evaluated to determine the feasibility of monitoring groundwater remediation activities associated with an underlying perched aquifer system using electrical resistivity tomography. The perched aquifer, located ~65 m below ground surface and ~10 m above the regional water table, contains high concentrations of nitrate, uranium, and other contaminants of concern from past tank leaks and intentional releases of wastes to surface disposal sites. The extent of the perched water aquifer is not well known, and the effectiveness of groundwater extraction for contaminant removal is uncertain, so supplemental characterization and monitoring technologies are being evaluated. Numerical simulations of subsurface flow and contaminant transport were performed with a highly resolved model of the hydrogeologic system and waste site infrastructure, and these simulations were used as the physical basis for electrical resistivity tomography modeling. The modeling explicitly accounted for metallic infrastructure at the site. The effectiveness of using surface electrodes versus surface and horizontal subsurface electrodes, for imaging groundwater extraction from the perched water aquifer, was investigated. Although directional drilling is a mature technology, its use for electrode emplacement in the deep subsurface under a complex industrial waste site via horizontal wells has not yet been demonstrated. Results from this study indicate that using horizontal subsurface electrode arrays could significantly improve the ability of electrical resistivity tomography to image deep subsurface features and monitor remediation activities under complex industrial waste sites.
Reactive Transport Modeling with Physics-Informed Machine Learning for Critical Minerals Applications
This study presents a physics-informed neural network (PINN) framework for reactive transport modeling for simulating fast bimolecular reactions in porous media. Accurate characterization of cAhemical interactions and product formation in surface and subsurface environments is essential for advancing critical mineral extraction and related geoscience applications. The proposed methodology sequentially addresses the flow and diffusion–reaction subproblems. The flow field is computed using a mixed formulation, while the diffusion–reaction system is modeled via two uncoupled tensorial diffusion equations reformulated in terms of chemical invariants. PINNs are employed to solve the governing equations, enabling data-efficient, mesh-free prediction of chemical concentration fields. The framework is validated through a series of benchmark problems involving flow in heterogeneous porous media. Initial verification is conducted using patch tests for the flow field, followed by validation of the transport problem with emphasis on preserving non-negativity of concentrations. The complete fast bimolecular reaction scenario is then solved, yielding spatial distributions of reactants and product species. Results demonstrate that the PINNs-based approach effectively captures sharp, mixing-limited reaction fronts and dispersive mixing behavior, offering reliable predictions of reactive plume evolution. These capabilities are crucial for evaluating long-term subsurface behavior in applications such as fluid storage, energy extraction, and efficient extraction of critical minerals.
Understanding Subsurface Fracture Evolution Dynamics Using Time‐Lapse Full Waveform Inversion of Continuous Active‐Source Seismic Monitoring Data
Abstract Predicting the behavior, geometry, and flow properties of subsurface fractures remains a challenging problem. Seismic models that can characterize fractures usually suffer from low spatiotemporal resolution. Here, we develop a correlative double‐difference time‐lapse full waveform inversion of continuous active source seismic monitoring data for determining high‐spatiotemporal‐resolution time‐lapse V p models of in‐situ fracture evolution at a shallow contamination site in Wyoming, USA. Assisted by rock physics modeling, we find that (a) rapidly increasing pore pressure initializes and grows the fracture, increasing the porosity slightly (from ∼13.7% to ∼14.6%) in the tight clay formation, thus decreasing V p (∼50 m/s); (b) the fluid injection continues decreasing V p , likely through the introduction of gas bubbles in the injectate; and (c) final V p reductions reach over ∼150 m/s due to a posited ∼4.5% gas saturation. Our results demonstrate that high‐resolution V p changes are indicative of mechanical and fluid changes within the fracture zone during hydrofracturing.
Real-Time Monitoring of Fracture Dynamics with a Contrast Agent-Assisted Electromagnetic Method
In collaboration with the Advanced Energy Consortium, our team has previously demonstrated that the placement of electrically active proppants (EAPs) in a hydraulic fracture surveyed by electromagnetic (EM) methods can enhance the imaging of the stimulated reservoir volumes during hydraulic fracturing. That work culminated in constructing a well-characterized EAP-filled fracture anomaly at the Devine field pilot site (DFPS). In subsequent laboratory studies, we observed that the electrical conductivity of our EAP correlates with changes in pressure, salinity, and flow. Thus, we postulated that the EAP could be used as an in-situ sensor for the remote monitoring of these changes in previously EAP-filled fractures. This paper presents our latest field data from the DFPS to demonstrate such correlations at an intermediate pilot scale. We conducted surface-based EM surveys during freshwater (200 ppm) and saltwater (2,500 ppm) slug injections while running surfaced-based EM surveys. Simultaneously, we measured the following: 1) bottomhole pressure and salinity in five monitoring wells; 2) injection rate using high-precision data loggers; 3) distributed acoustic sensors in four monitoring wells; and 4) tiltmeter data on the survey area. We demonstrated that injections into an EAP-filled fracture could be successfully coupled with real-time electric field measurements on the surface, leading to remote monitoring of dynamic changes within the EAP-filled fracture. Furthermore, by comparing the electrical field traces with the bottomhole pressure, flow rate, and salinity, we concluded that the observed electric field in our study is influenced by fracture dilation and flow rate. Salinity effect was observed when saltwater was injected. EM simulations solely based on assumptions of fracture conductivity changes during injection did not reproduce all of the measured electric field magnitudes. Preliminary estimates showed that including streaming potential in our geophysical model may be needed to reduce the simulation mismatch. The methods developed and demonstrated during this study will lead to a better understanding of the extent of fracture networks, formation stress states, fluid leakoff and invasion, characterizations of engineered fracture systems, and other applications where monitoring subsurface flow tracking is deemed important.
Experimental Investigation of Buoyant Flow in Realistic Bedforms With Heterogeneous Wettability
Submeter-scale geologic heterogeneity greatly affects CO 2 plume migration and retention. In this work, we present meter-scale laboratory experiments that can capture the impact of realistic submeter-scale geologic heterogeneity on multiphase flow and trapping. We produce realistic sedimentary formations consisting of ripple deposits with varying grain size contrast and wettability in a meter-scale slab chamber. Then, we conduct multiphase flow experiments with analog fluids through these structures and measure the saturation patterns, capillary heterogeneity trapping (CHT), and overall trapping performance. When we alter the ripple bedform architecture, variations in trapped saturation and CHT (10–20%) increment are exhibited. Similar growth in trapping performance is also observed when grain size contrast increases. Finally, wettability changes (water- to oil-wet) can increase nonwetting saturation and CHT up to 5% and 10–20%, respectively. These results emphasize the importance of correctly characterizing the impact of small-scale heterogeneities and wettability changes. We believe this is the first time that multiphase flow experiments were conducted in meter-scale domains with realistic ripple bedforms and heterogeneous wettability to investigate plume migration and trapping.
Phase I Investigation at Bladen, Nebraska (Final Work Plan)
The village of Bladen is a town of population approximately 237 in the northwest part of Webster County, Nebraska, 30 mi southwest of Hastings and 140 mi southwest of Lincoln, Nebraska. In 2000, the fumigant-related compound carbon tetrachloride was detected in public water supply well PWS 68-1, at a trace level. Low-level contamination, below the maximum contamination level (MCL) of 5.0 μg/L, has been detected intermittently in well PWS 68-1 since 2000, including in the last sample taken in July 2013. In 2006, the village installed a new well, PWS 2006-1, that remains free of contamination. Because the carbon tetrachloride found in well PWS 68-1 might be linked to historical use of fumigants containing carbon tetrachloride at grain storage facilities, including its former facility in Bladen, the CCC/USDA is proposing an investigation to (1) delineate the source and extent of the carbon tetrachloride contamination potentially associated with its former facility, (2) characterize pathways and controlling factors for contaminant migration in the subsurface, and (3) establish a basis for estimating potential health and environmental risks. The work will be performed in accordance with the Intergovernmental Agreement established between the NDEQ and the Farm Service Agency of the USDA. The site investigation at Bladen will be implemented in phases, so that data collected and interpretations developed during each phase can be evaluated to determine if a subsequent phase of investigation is warranted and, if warranted, to provide effective guidance for the subsequent investigation activities. This Work Plan identifies the specific technical objectives and defines the scope of work proposed for the Phase I investigation by compiling and evaluating historical data. The proposed investigation activities will be performed on behalf of the CCC/USDA by the Environmental Science Division of Argonne National Laboratory. Argonne is a nonprofit, multidisciplinary research institute operated by UChicago Argonne, LLC, for the U.S. Department of Energy.
Distributed Acoustic Sensing of Strain at Earth Tide Frequencies
The solid Earth strains in response to the gravitational pull from the Moon, Sun, and other planetary bodies. Measuring the flexure of geologic material in response to these Earth tides provides information about the geomechanical properties of rock and sediment. Such measurements are particularly useful for understanding dilation of faults and fractures in competent rock. A new approach to measuring earth tides using fiber optic distributed acoustic sensing (DAS) is presented here. DAS was originally designed to record acoustic vibration through the measurement of dynamic strain on a fiber optic cable. Here, laboratory experiments demonstrate that oscillating strain can be measured with DAS in the microHertz frequency range, corresponding to half-day (M2) lunar tidal cycles. Although the magnitude of strain measured in the laboratory is larger than what would be expected due to earth tides, a clear signal at half-day period was extracted from the data. With the increased signal-to-noise expected from quiet field applications and improvements to DAS using engineered fiber, earth tides could potentially be measured in deep boreholes with DAS. Because of the distributed nature of the sensor (0.25 m measurement interval over kilometres), fractures could be simultaneously located and evaluated. Such measurements would provide valuable information regarding the placement and stiffness of open fractures in bedrock. Characterization of bedrock fractures is an important goal for multiple subsurface operations such as petroleum extraction, geothermal energy recovery, and geologic carbon sequestration.
Advancements in 3D Structural Analysis of Geothermal Systems
Robust geothermal activity in the Great Basin, USA is a product of both anomalously high regional heat flow and active fault-controlled extension. Elevated permeability associated with some fault systems provides pathways for circulation of geothermal fluids. Constraining the local-scale 3D geometry of these structures and their roles as fluid flow conduits is crucial in order to mitigate both the costs and risks of geothermal exploration and to identify blind (no surface expression) geothermal resources. Ongoing studies have indicated that much of the robust geothermal activity in the Great Basin is associated with high density faulting at structurally complex fault intersection/interaction areas, such as accommodation/transfer zones between discrete fault systems, step-overs or relay ramps in fault systems, intersection zones between faults with different strikes or different senses of slip, and horse-tailing fault terminations. These conceptualized models are crucial for locating and characterizing geothermal systems in a regional context. At the local scale, however, pinpointing drilling targets and characterizing resource potential within known or probable geothermal areas requires precise 3D characterization of the system. Employing a variety of surface and subsurface data sets, we have conducted detailed 3D geologic analyses of two Great Basin geothermal systems. Using EarthVision (Dynamic Graphics Inc., Alameda, CA) we constructed 3D geologic models of both the actively producing Brady’s geothermal system and a ‘greenfield’ geothermal prospect at Astor Pass, NV. These 3D models allow spatial comparison of disparate data sets in 3D and are the basis for quantitative structural analyses that can aid geothermal resource assessment and be used to pinpoint discrete drilling targets. The relatively abundant data set at Brady’s, ~80 km NE of Reno, NV, includes 24 wells with lithologies interpreted from careful analysis of cuttings and core, a 1:24,000 scale detailed geologic map and cross-sections, 2D seismic reflection profiles and other geophysical data, and downhole temperature data. The 3D geologic model based on these data consists of 61 fault planes, 25 distinct stratigraphic units, and 2 intrusive bodies. Geothermal fluids are produced from a left step-over/relay ramp within the Brady’s Fault Zone (BFZ). Under local stress conditions, fault segments that strike NNE-to-NE are most likely to slip and/or dilate, and therefore transmit geothermal fluids. The 3D model defines the locations of discrete fault intersections within the BFZ and indicates that the densest zones of structurally controlled fracture permeability are ~10-to-10s of meters in diameter and plunge ~55° NW-NNW beneath the heart of the BFZ step over. The locations of high intersection density, high fault slip and dilation tendency, high subsurface temperature, and lithologies known to support high fracture permeability are combined to produce 3D ‘fairway’ maps useful in both assessments of geothermal resource potential and for defining drilling targets. Astor Pass is located on the Pyramid Lake Paiute Reservation, ~80 km north of Reno, NV. It is a prospective ‘greenfield’ geothermal area, and thus subsurface data are relatively sparse. Available data include: two relatively deep wells (~1400 m) and one shallower well (~500 m) with lithologies interpreted from drill cuttings, several 2D seismic reflection profiles, a 1:24,000 scale geologic map and cross-section, a shallow temperature survey, and downhole temperature data. 3D modeling based on these data has defined 19 distinct fault planes and 16 stratigraphic units. Based on the stress field calculated from borehole breakouts, drilling induced tensile cracks and petal-centerline cracks in the two relatively deep wells, 3D slip and dilation tendency analysis indicates that northerly striking fault segments are most likely to slip and/or dilate, and therefore transmit geothermal fluids. Analysis of fault intersection density indicates that the highest density of structurally controlled permeability within the field lies in a narrow (10-to-10s of m) zone plunging moderately (~35°) to the NNW beneath Pleistocene tufa deposits. This zone of increased fracture density, which we interpret as the primary upflow zone, is controlled by the intersection of N-to-NNW striking normal faults and a WNW striking dextral fault zone and represents the most promising target for future drilling. Construction of a 3D geologic model involves integration of a variety of data into an internally consistent framework. A robust model allows for spatial comparison between the various types of data (structural, stratigraphic, geophysical, temperature, etc.) that are commonly used independently to site geothermal wells. Furthermore, highly detailed 3D geologic models provide the basis for additional quantitative analysis, including 3D fault slip and dilation tendency analysis and the precise location of structurally controlled permeability pathways. These analyses provide detailed information relating to the internal dynamics of geothermal systems and can mitigate the costs and risks of geothermal exploration and development by contributing to better well targeting and more accurate evaluations of resource potential.
Moisture Monitoring at Area G, Technical Area 54, Los Alamos National Laboratory
Hydrological characterization and moisture monitoring activities provide data required for evaluating the transport of subsurface contaminants in the unsaturated and saturated zones beneath Area G, and for the Area G Performance Assessment and Composite Analysis. These activities have been ongoing at Area G, Technical Area 54 of the Los Alamos National Laboratory since waste disposal operations began in 1957. This report summarizes the hydrological characterization and moisture monitoring activities conducted at Area G. It includes moisture monitoring data collected from 1986 through 2014 from numerous boreholes and access tubes with neutron moisture meters, as well as data collected by automated dataloggers for water content measurement sensors installed in a waste disposal pit cover, and buried beneath the floor of a waste disposal pit. This report is an update of the Jennings and French (2009) moisture monitoring report and includes all data from Jennings and French (2009) and the Draft 2010 Addendum moisture monitoring report (Jennings and French, 2010). For this updated report, all neutron logging data, including neutron probe calibrations, were investigated for quality and pedigree. Some data have been recalculated using more defensible calibration data. Therefore, some water content profiles are different from those in the Jennings and French (2009) report. Monitoring and characterization data generally indicate that some areas of the Area G vadose zone are consistent with undisturbed conditions, with water contents of less than five percent by volume in the top two layers of the Bandelier tuff at Area G. These data also indicate that other areas of the vadose zone are affected by waste disposal activities that have been ongoing at Area G since 1957, a period of nearly 60 years. In some areas, water content profiles indicate increases in water content to depths of tens of meters, especially in areas covered by asphalt and structures.
Moisture Monitoring at Area G, Technical Area 54, Los Alamos National Laboratory, 2016 Status Report
Hydrological characterization and moisture monitoring activities provide data required for evaluating the transport of subsurface contaminants in the unsaturated and saturated zones beneath Area G, and for the Area G Performance Assessment and Composite Analysis. These activities have been ongoing at Area G, Technical Area 54 of the Los Alamos National Laboratory since waste disposal operations began in 1957. This report summarizes the hydrological characterization and moisture monitoring activities conducted at Area G. It includes moisture monitoring data collected from 1986 through 2016 from numerous boreholes and access tubes with neutron moisture meters, as well as data collected by automated dataloggers for water content measurement sensors installed in a waste disposal pit cover, and buried beneath the floor of a waste disposal pit. This report is an update of a nearly identical report by Levitt et al., (2015) that summarized data collected through early 2015; this report includes additional moisture monitoring data collected at Pit 31 and the Pit 38 extension through December, 2016. It also includes information from the Jennings and French (2009) moisture monitoring report and includes all data from Jennings and French (2009) and the Draft 2010 Addendum moisture monitoring report (Jennings and French, 2010). For the 2015 version of this report, all neutron logging data, including neutron probe calibrations, were investigated for quality and pedigree. Some data were recalculated using more defensible calibration data. Therefore, some water content profiles are different from those in the Jennings and French (2009) report. All of that information is repeated in this report for completeness. Monitoring and characterization data generally indicate that some areas of the Area G vadose zone are consistent with undisturbed conditions, with water contents of less than five percent by volume in the top two layers of the Bandelier tuff at Area G. These data also indicate that other areas of the vadose zone are affected by waste disposal activities that have been ongoing at Area G since 1957, a period of nearly 60 years. In some areas, water content profiles indicate increases in water content to depths of tens of meters, especially in areas covered by asphalt and structures.
Beyond the Surface: Non-Invasive Low-Field NMR Analysis of Microbially-Induced Calcium Carbonate Precipitation in Shale Fractures
Microbially-induced calcium carbonate precipitation (MICP) is a biological process in which microbially-produced urease enzymes convert urea and calcium into solid calcium carbonate (CaCO 3 ) deposits. MICP has been demonstrated to reduce permeability in shale fractures under elevated pressures, raising the possibility of applying this technology to enhance shale reservoir storage safety. For this and other applications to become a reality, non-invasive tools are needed to determine how effectively MICP seals shale fractures at subsurface temperatures. In this study, two different MICP strategies were tested on 2.54 cm diameter and 5.08 cm long shale cores with a single fracture at 60 °C. Flow-through, pulsed-flow MICP-treatment was repeatedly applied to Marcellus shale fractures with and without sand (“proppant”) until reaching approximately four orders of magnitude reduction in apparent permeability, while a single application of polymer-based “immersion” MICP-treatment was applied to an Eagle Ford shale fracture with proppant. Low-field nuclear magnetic resonance (LF-NMR) and X-Ray computed microtomography (micro-CT) techniques were used to assess the degree of biomineralization. With the flow-through approach, these tools revealed that while CaCO 3 precipitation occurred throughout the fracture, there was preferential precipitation around proppant. Without proppant, the same approach led to premature sealing at the inlet side of the core. In contrast, immersion MICP-treatment sealed off the fracture edges and showed less mineral precipitation overall. This study highlights the use of LF-NMR relaxometry in characterizing fracture sealing and can help guide NMR logging tools in subsurface remediation efforts.