Earth & Environmental Sciences

This study presents the results from in-well hydraulic fracture monitoring within a horizontal well in an unconventional reservoir utilizing Distributed Fiber Optic Sensing (DFOS). An in-house-developed Brillouin-based Distributed Strain Sensing (DSS) interrogator was deployed to obtain strain measurements, complemented by a commercial Raman-based Distributed Temperature Sensing (DTS) interrogator for temperature measurements and a commercial Rayleigh-based Low-Frequency Distributed Acoustic Sensing (LF-DAS) interrogator for strain-rate measurements. Examined over a ten-day period, the spatio-temporal distribution of temperature-compensated strain obtained from DSS and DTS revealed distinct signatures of the multi-stage hydraulic fracturing process. These signatures were analyzed with respect to fracture width growth and closure, residual strain effects, and fracture conductivity near the wellbore. Fracture widths within the fracture zone were estimated for individual stages. The findings were assessed with LF-DAS measurements for further evaluation. This work integrates DFOS-measured strain, temperature, and strain-rate data for monitoring in-well hydraulic fracturing, with the goal of supporting future studies in interpreting DFOS measurements for improved understanding of hydraulic fracturing in unconventional reservoirs.

Executive Summary Scientists at Lawrence Berkeley National Laboratory (Berkeley Lab) have teamed up with the University of Hawai‘i at Manoa (UH Manoa) through the U.S. Department of Energy’s Energy Technology Innovation Partnership Project to evaluate the technological and market feasibility of shallow geothermal heat exchanger (GHE) technology. UH requested this analysis to evaluate opportunities in building cooling, energy efficiency, and emissions reduction applications in Hawai‘i. UH has an abundance of geologic and geothermal data and is looking to the national labs’ expertise to execute this analysis. UH is also interested in investigating policy, regulatory, and business conditions advantageous for implementation of a pilot project and more broad deployment of this technology in Hawai‘i. In many locations around the world, the demands for heating and cooling are roughly balanced over the course of the year, so GHEs do not cause significant long-term changes in subsurface temperature. This is not the case in Hawai’i, where the demand for heating is very small, meaning that, over time, GHEs will add heat to the subsurface. If temperatures increase significantly, GHE systems will not work as designed. Regional groundwater flow has the potential to sweep heated water away from boreholes, thereby maintaining the functionality of the GHE system. Significant regional groundwater flow requires two things: a sufficiently large driving hydraulic head gradient (usually closely related to surface topography), and sufficient porosity and permeability to enable groundwater to flow in large enough quantities to enable near-borehole temperatures to be maintained at ambient values. Hawai‘i’s volcanic terrain offers ample surface topographic variation. The lava itself shows an extremely large range of porosity and permeability, so sites with large enough values of these properties must be selected. Numerical modeling of coupled groundwater and heat flow can be used to determine how large is large enough. Primarily, closed-loop systems have been investigated. Other options considered are open-loop systems and using cool seawater as the chilling source. Project work investigated the feasibility of GHE technology at two scales. At the island scale, GIS layers of various attributes relevant for GHE were combined to develop an overall favorability map for employing GHE in Hawai‘i. At the local scale, a hydrogeologic model for the subsurface component of a closed-loop system was developed for the Stan Sheriff Center at the UH Manoa campus. This site is considered promising because the rock below and immediately downgradient of the borefield is highly permeable, consisting of a subsurface karst system (limestone containing high-permeability open channels), which is underlain by a thick, high-permeability fractured basalt. Moreover, the site is near the base of the Ko‘olau Range, providing a large hydraulic head gradient. Thus, groundwater flow through the site is expected to be large, enabling efficient removal of heated groundwater. A full-GHE-system model of the site was also developed, with a simplified representation of the subsurface, in which groundwater flow is not considered and heat transfer is purely by conduction. Using the building cooling load data provided by UH, simulation results show that with groundwater flow present, a GHE can operate successfully for at least 10 years, but with no groundwater flow, the subsurface begins to heat up after only one year of operation, making the GHE unviable within 2-6 years. The team also developed a techno-economic model for this site to compare the cost of cooling using a GHE system with the costs of operating the current air-conditioning system. The GHE system is advantageous economically if favorable tax incentives and interest rates can be obtained.

This multi-lab, DOE-funded project addresses the significant energy and water consumption and cost to cool information technology (IT) equipment in data centers by utilizing subsurface thermal energy storage systems, more specifically, reservoir thermal energy storage (RTES). The project was augmented by an industrial advisory group (IAG), including experts from both the data center and subsurface energy storage sectors, to provide feedback. A scenario-based method was applied to perform techno-economic feasibility analysis based on three types of data centers covering a range of sizes and in three geographical locations. The techno-economic analysis (TEA) was performed to compare RTES scenarios with commonly used or most competitive non-RTES cooling scenarios. The main conclusions from the investigation are that all RTES systems studied are technically feasible and sustainable for at least a period of 20 years without major modifications of the RTES and IT cooling systems. Within the context of the assumptions made by this study, the key factor to make RTES for data center cooling economically feasible and attractive in the right location includes: 1) a shallow non-potable water-bearing geological formation with large transmissivity (thick and high permeability formation) to maximize storability and minimize the number and depth of wells needed; and 2) potential to use free (compressorless) or inexpensive cooling. Compressorless cooling can be provided by dry coolers in mild climates (although that is not the only option,) and inexpensive cooling can utilize compressor cooling when power costs are very low or negative (e.g., excessive renewable energy production). Future studies should further consider using chillers for RTES cooling (in addition to dry coolers) when there is a significant grid value to do so (large difference between peak and off-peak power cost). Additionally, system optimization should be performed for a specific site to maximize the benefit of using RTES for cooling when deployed. Additional benefits, such as resiliency during high heat events, are often not captured in traditional TEA studies, and should be considered.

Molecular hydrogen (H2) is found in a variety of settings on and in the Earth from low-temperature sediments to hydrothermal vents, and is actively being considered as an energy resource for the transition to a green energy future. The hydrogen isotopic composition of H2, given as D/H ratios or δD, varies in nature by hundreds of per mil from ∼−800 ‰ in hydrothermal and sedimentary systems to ∼+450 ‰ in the stratosphere. This range reflects a variety of processes, including kinetic isotope effects associated with formation and destruction and equilibration with water, the latter proceeding at fast (order year) timescales at low temperatures (<100 °C). At isotopic equilibrium, the D/H fractionation factor between liquid water and hydrogen (DαH2O(l)-H2(g)) is a function of temperature and can thus be used as a geothermometer for H2 formation or re-equilibration temperatures. Multiple studies have produced theoretical calculations for hydrogen isotopic equilibrium between H2 and water vapor. However, only three published experimental calibrations used in geochemistry exist for the H2O-H2 system: two between 51 and 742 °C for H2O(g)-H2(g) (Suess, 1949; Cerrai et al., 1954), and one in the H2O(l)-H2(g) system for temperatures <100 °C (Rolston et al., 1976). Despite these calibrations existing, there is uncertainty on their accuracy at low temperatures (<100 °C; e.g., Horibe and Craig, 1995). Here we present a new experimental calibration of the equilibrium hydrogen isotopic fractionation factor for liquid water and molecular hydrogen from 3 to 90 °C. Equilibration was achieved using platinum catalysts and verified via experimental bracketing by approaching final values of DαH2O(l)-H2(g) at a given temperature from both higher (top-bracket) and lower (bottom-bracket) initial Dα values. Our calibration yields the following equation: [Formula presented] Where T is in Kelvin. We find that our calibrations differ from prior experimental calibrations by, on average, up to 20 ‰ and prior theoretical results by up to, on average, 25 ‰. Good agreement with theoretical results (<11 ‰ differences) is found for calculations that consider both anharmonic effects and the Diagonal Born-Oppenheimer correction.

Direct lithium extraction (DLE) of brines after geothermal power production offers opportunities to produce environmentally benign “green” lithium; however, some environmental impact is inevitable. We examined solid waste production at geothermal power plants in southern California that are also locations for planned DLE facilities. Currently, the geothermal plants in this region produce approximately 79,800 metric tons (wet weight) per year of solid waste, which represents about 28 metric tons per GWh of net electricity production or approximately 500 mg solids per kg geothermal brine. Approximately 15% of this waste requires management as hazardous waste. Solids produced during power production represent about 0.2% of the total dissolved solids in the brine. Lithium production will require the removal of silica, iron, and other metals as part of the DLE process. Using a mass balance approach, we calculate that precipitation of silica and metals could produce up to an additional 6800 mg solids per kg brine. Calcium occurs at very high concentrations, and the amount of solids disposed in landfills will be dependent on the amount of calcium removed during lithium recovery. Our analysis shows that evaluation of brine chemistry in the context of the DLE process is useful for evaluating the potential solid waste impacts of producing lithium from brines.

Information about evapotranspiration (ET) and its components, that is, evaporation and transpiration, is crucial for a wide range of water and ecosystem management applications. However, partitioning ET into its two components is often challenging because of their spatiotemporal variabilities and lack of process understanding. This study developed a machine learning (ML) framework to shed light on ET processes and assess the relative importance of different drivers by incorporating hydrometeorology and biomass productivity variables. The Shapley Additive Explanations (SHAP) approach was applied to enhance explainability and rank the importance of ET drivers and their components. A total of 62 variables covering hydrometeorological and biomass productivity dimensions were considered from the Reynolds Creek Critical Zone Observatory (CZO) station in Idaho. The variable importance assessment identified the leading drivers individually for evaporation, transpiration and ET (soil water content for evaporation, vapour pressure deficit for transpiration and soil water content for ET). The results further highlighted the value of combining hydrometeorological and biomass productivity variables to achieve better predictability of ET processes.

Land surface models consider the exchange of water, energy, and carbon along the soil-canopy-atmosphere continuum, which is challenging to model due to their complex interdependency and associated challenges in representing and parameterizing them. Differentiable modeling provides a new opportunity to capture these complex interactions by seamlessly hybridizing process-based models with deep neural networks (DNNs), benefiting both worlds, that is, the physical interpretation of process-based models and the learning power of DNNs. Here, we developed a differentiable land model, JAX-CanVeg. The new model builds on the legacy CanVeg by incorporating advanced functionalities through JAX in the graphic processing unit support, automatic differentiation, and integration with DNNs. We demonstrated JAX-CanVeg's hybrid modeling capability by applying the model at four flux tower sites with varying aridity. To this end, we developed a hybrid version of the Ball-Berry equation that emulates the water stress impact on stomatal closure to explore the capability of the hybrid model in (a) improving the simulations of latent heat fluxes (Formula presented.) and net ecosystem exchange (Formula presented.), (b) improving the optimization trade-off when learning observations of both (Formula presented.) and (Formula presented.), and (c) benefiting a multi-layer canopy model setup. Our results show that the proposed hybrid model improved the simulations of (Formula presented.) and (Formula presented.) at all sites, with an improved optimization trade-off over the process-based model. Additionally, the multi-layer canopy set benefited hybrid modeling at some sites. Anchored in differentiable modeling, our study provides a new avenue for modeling land-atmosphere interactions by leveraging the benefits of both data-driven learning and process-based modeling.

A model for gas migration through clay-based buffer material is developed for modeling gas generation and migration associated with deep geologic nuclear waste disposal. The model is based on a multiphase fluid flow and geomechanics simulator that is adapted to consider enhanced gas flow when gas pressure is high enough to approach the confining stress magnitude. A key feature in the model is a direct coupling between gas permeability and stress, through a non-linear stress-dependent permeability function. The model was first tested and calibrated by modelling two different laboratory gas migration tests on Wyoming (MX-80) bentonite samples. The calibrated model was then applied to model gas migration through a bentonite buffer of a large-scale gas injection test (Lasgit) conducted at the Äspö Hard Rock Laboratory in Sweden. Observed preferential gas migration along interfaces (between compacted blocks and along the canister surface) required explicit representation of such interfaces in the model. The model with the stress-dependent gas permeability accurately captured observed experimental responses in terms of gas breakthrough time, peak gas pressure, and cumulative gas flow rates. The calibrated model was finally applied to simulate migration of hydrogen gas generated within a breached nuclear waste canister over 10,000 years, involving migration of much larger gas volumes. For the considered gas generation rate and host rock properties, the generated gas could migrate through the bentonite buffer and released into the surrounding host rock at a maximum gas pressure somewhat higher than the initial total stress, though a significant amount of hydrogen remained within the buffer. This modelling sets the stage for further detailed analysis of the impact of hydrogen gas generation on the long-term performance of nuclear waste repositories.

It is challenging to quantitatively predict shearing of intersecting fractures/faults because of dynamic frictional contacts accompanied by possible nonlinear rock deformation. To address such challenges, a new conceptual model—the simplified DFN model—was proposed and validated by Hu et al.46 to use major paths (MPs) to represent complicated DFNs for calculation of shearing. In this work, we conducted a benchmark study for three examples that involve different levels of complexity of intersecting fractures, and correspondingly different numbers of MPs. The codes and software that were used in the benchmark cover a range of continuum, discontinuum and hybrid numerical methods: NMM (LBNL), FLAC3D (LBNL), GBDEM (KIGAM), FRACOD (DynaFrax), and CASRock (CAS). The general consistency between DFN and MP cases as predicted by all the codes/software demonstrates that major paths can be used to simplify the geometry of DFNs in a wide range of software. Disagreement in results made by some software and potential future improvements are discussed. We show that (1) shearing of one or multiple major fractures can be reduced if there are multiple smaller intersecting fractures in that area, which is a useful basis for understanding and controlling induced seismicity and merits further analysis, and (2) the agreement achieved in the benchmark examples provide confidence that the simplified DFN model is a promising conceptual model that can be used for different types of numerical approaches and software for simplifying the analysis of the shearing of intersecting fractures and faults.