Earth & Environmental Sciences

Previous studies have shown that underground natural gas storage (UGS) in California has served a critical role in meeting energy demands in California, and there is no immediate alternative. Therefore, it is important to ensure the safety of UGS infrastructure, especially considering that many of the UGS sites are using a combination of new and old wells, some of which were installed decades ago and re-purposed for UGS. The purpose of this project is to develop an integrated risk management and decision support system (IRMDSS) to manage risks associated with this heterogeneous subsurface infrastructure. The approach of the IRMDSS is to take advantage of the predictive capability of mechanistic models, with support from data acquired from advanced monitoring technologies, for evaluation and analysis of various incident scenarios or potential threats. In this project, we have demonstrated data collection by four advanced monitoring technologies. These include two downhole monitoring technologies, distributed temperature sensing (DTS), distributed acoustic sensing (DAS), and two surface monitoring technologies, Interferometric Synthetic Aperture Radar (InSAR), and unmanned aerial vehicle (UAV). DTS and DAS data are collected continuously, providing information related to individual wells. InSAR data are collected frequently (~every 24 days), and UAV data can be collected as frequently as is practical depending on need. Together, these subsurface and surface monitoring technologies provide near real-time information useful for risk management of UGS facilities.

Radionuclide transport in smectite clay barrier systems used for nuclear waste disposal is controlled by diffusion, with adsorption significantly retarding transport rates. While a relatively minor component of spent nuclear fuel, 79Se is a major driver of the safety case for spent fuel disposal due to its long half-life (3.3 × 105 yr) and its low adsorption to clay (KD < 10 L/kg), thus a thorough understanding of Se diffusion through clay is critical for understanding the long-term safety of spent fuel disposal systems. Through-diffusion experiments with tritiated water (HTO, conservative tracer) and Se(VI) were conducted with a well-characterized, purified montmorillonite source clay (SWy-2) under a constant ionic strength (0.1 M) and three different electrolyte compositions: Na+, Ca2+, and a Na + -Ca2+ mixture at pH 6.5 in order to probe the effects of electrolyte composition and interlayer cation composition on clay microstructure, Se(VI) aqueous speciation, and ultimately diffusion. The results were modeled using a reactive transport modeling approach to determine values of porosity (ε), De (effective diffusion coefficient), and KD (distribution coefficient for adsorption). HTO diffusive flux was higher in Ca-montmorillonite (De = 1.68 × 10−10 m2 s−1) compared to Na-montmorillonite (De = 7.83 × 10−11 m2 s−1). This increase in flux is likely due to a greater degree of clay layer stacking in the presence of Ca2+ compared to Na+, which leads to larger inter-particle pores. Overall, the Se(VI) flux was much lower than the HTO flux due to anion exclusion, with Se(VI) flux following the order Ca (De = 1.03 × 10−11 m2 s−1) > Na–Ca (De = 2.12 × 10−12 m2 s−1) > Na (De = 1.28 × 10−12 m2 s−1). These differences in Se(VI) flux are due to a combination of factors, including (1) larger accessible porosity in Ca-montmorillonite due to clay layer stacking and smaller electrostatic effects compared to Na-montmorillonite, (2) larger accessible porosity for neutral-charge CaSeO4 species which makes up 32% of aqueous Se(VI) in the pure Ca system, and (3) possibly higher Se(VI) adsorption for Ca-montmorillonite. Through a combination of experimental and modeling work, this study highlights the compounding effects that electrolyte and counterion compositions can have on radionuclide transport through clay. Diffusion models that neglect these effects are not transferable from laboratory experimental conditions to in situ repository conditions.

We evaluate if electromagnetic (EM) geophysical methods for monitoring geologic carbon storage (GCS) efforts at the Wyoming CarbonSAFE project adjacent to the Dry Fork Station power plant near Gillette, Wyoming. This first involved acquiring both electric and magnetic fields at eleven different locations ranging in distance from immediately adjacent to 4 km from the plant. Passive EM measurements were made to provide spectral EM noise measurements generated by electricity production at the plant and to determine if useful magnetotelluric (MT) data can be successfully collected in the region. The processed data indicate that useful MT data can be collected as long as the site is located more than 2km away from the power plant as well as active roads and rail lines. Controlled source EM data were collected using three different source configurations, two of which connected to steel casings used to complete the injection wells. Comparing the EM noise measurements to the CSEM data show measurable electric and magnetic field signals at all sites. Next a series of three-dimensional (3D) numerical models were built that simulate resistivity changes caused by the proposed CO2 injection at depths ranging from 2.4 to 3.0km. These models were used to simulate various EM measurement configurations. The modeling shows that casing-source CSEM monitoring can provide sensitivity to the injected CO2 if source electrodes are connected to the bottom of one or both of the injection wells.

Mollusk shells protect the animals that form and inhabit them. They are composites of minerals and organics, with diverse mesostructures, including nacre, prismatic calcite, crossed-lamellar aragonite, and foliated calcite. Twins, that is, crystals mirror symmetric with respect to their coherent interface, occurring as formation or deformation twins, are observed in all mollusk shell mesostructures but never within calcite prisms. Here, nanotwins and microwins within single calcite prisms are observed in different shells. Using Polarization-dependent Imaging Contrast (PIC) mapping with 20–60 nm resolution, twins are observed to be 0.2–3 µm thick layers of differently oriented and colored crystals with respect to the main prism crystal. Multiple twins are interspersed with the prism crystal, parallel to one another, and similarly oriented. When comparing images of calcite prisms and twins obtained by PIC mapping and by Electron Back-Scattered Diffraction (EBSD), the images correspond precisely. All twins are e-twin types, with 127° angular distance between c-axes. E-twins are the most common deformation twins in geologic calcite, as also observed here in Carrara marble. Location of all twins near the outer surface of all shells and e-twin type both suggest that twins within calcite prisms in mollusk shells result from deformation twinning.

Geological CO2 sequestration (GCS) can help mitigate global warming and enhance methane recovery from coal beds. However, few studies have linked the effects of CO2 to surface chemistry changes controlling wetting behavior in deep coal beds. Contact angles (CAs) of CO2/N2-high volatile bituminous coal-water systems were measured under different temperatures and pressures. The surface chemistry and physical structure of coals were characterized to investigate changes in physicochemical properties and their relations with wettability after reactions. For N2 treatment, the time-dependence of static and dynamic CAs were insignificant, ranging within 4°. For gaseous CO2 treatment, the static CAs and the average advancing angles increased slightly. With supercritical (sc) CO2, both the static and dynamic CAs increased significantly, and θ adv changed to intermediate-wet (92°). Reactions with minerals exposed to scCO2 resulted in greater surface roughness and heterogeneity, greater contact angle hysteresis and more surface sites occupied by scCO2 rather than H2O. Increases in hydrophobic functional groups and decreases in hydrophilicity were shown by FTIR spectra, reflecting the shedding of polar oxygen-containing functional groups, reduction of hydrogen bonds, and increasing percentage of hydrocarbons. XRD patterns obtained following scCO2-treatment showed that crystallite growth and molecular polymerization were higher toward graphite-like. The calculated structural parameters of functional groups and crystallites both showed elevated coal rank. Changes in crystallite structure, notably higher carbon content and decreased negative surface charge, are unfavorable for water-wetting. This study contributes to understanding surface chemistry changes responsible for decreased wettability during CO2-enhanced coal bed methane recovery and GCS in coal reservoirs.

Abstract: Coupled thermo-hydro-mechanical (THM) processes in fractured rocks have been a topic of intense scientific research for more than 30 years. The present paper takes a look into the past and highlights some scientific advances which are of an unusual “out-of-the-box” nature, and then looks forward and discusses possible directions of future research in this interesting field of study. Concerning future research directions, we see a trend from a focus on coupled THM processes in single fractures or a few interacting fractures, to the study of coupled THM behavior in complex fracture network systems where the fractures act collectively giving rise to local stress concentration points and points of large pressure gradients. Three examples of future research directions are presented. First is an effort towards identifying characterizing parameters of a fracture network that play a direct controlling role in major coupled THM phenomena (such as induced seismicity and flow channeling), rather than parameters of stochastic distributions of fractures in the network. The second example of research direction is accounting for the heterogeneity and hierarchy of fractures in a fault or fracture zone which has been associated with major THM events in a number of geo-energy projects. The third example is at the opposite end of the first; here it is recognized that in some cases, the coupled THM processes in fractured rocks may be controlled dominantly by only a few key bridges. Identification, characterization, and evaluation of these key bridges should be one of the important research directions in the coming days.

Optimization of aquifer thermal energy storage (ATES) performance in a building system is an important topic for maximizing the seasonal offset between energy demand and supply and minimizing the building's primary energy consumption. To evaluate ATES performance with bidirectional operation, this study develops an analytical solution-based model to simulate the spatiotemporal thermal response in an aquifer. The model consists of three temperature response functions, similar to the G functions in borehole thermal energy storage (BTES), to estimate the transient temperature profile in the aquifer during seasonally varying injection and extraction of hot/cold water. Applying machine learning (ML) based data classification and regression techniques to the results of a series of finite element (FE) benchmark simulations of typical ATES configurations, model input parameters are linked to the subsurface thermal, hydrogeological, and ATES operational properties. Compared to the benchmark simulation results, the errors of the proposed model in estimating the annual energy storage and locating the thermally affected area are about 3 % and 1 %, respectively. The model was applied to a previous short-term case study, and the error in the transient production temperature estimation is about 1 %. The long-term heat recovery ratio estimated from the model also compares well to those calculated from the previous study and the validated numerical model. Because of its fast computation, the proposed model can be coupled with the individual building system simulation and used for preliminary ATES design, and this will allow for greater exploration of ATES operational space and, therefore, better choices of ATES operating conditions. The proposed model can also be coupled with the district heating and cooling network simulation for computationally efficient city-scale long-term ATES potential assessment.

Shearing of fractures and faults is important because it can result in permeability change or even induce seismicity—both are keys for efficient and safe energy recovery and storage in Earth systems. Quantitative analysis of shearing of intersecting fractures and faults is challenging because it can involve dynamic frictional contacts that are complicated by deformation of the rock matrix. To predict the shearing of intersecting fractures/faults, we attempt to answer the question of how intersections impact the shearing of a fracture network and whether we can simplify the description as compared to classical discrete fracture network (DFN) models. To answer these questions, we conducted a series of numerical simulations on scenarios for variable numbers of intersecting fractures. All these examples yield consistent results: the results of using DFNs are consistent with those of using hypothetical major paths. This leads to a new model, which we name simplified discrete fracture network model, to analyze shearing of intersecting fractures/faults using major path(s). We found that the intersections of fractures do not fundamentally change the shearing of two intersecting fractures if the intersecting angles are small. Furthermore, increasing the number of fractures/faults may relax the stress as more fractures/faults become available for shearing and distributing the stress. The simplified DFN model, which can capture efficiently the shearing behavior of each major paths from a large number of intersecting fractures/faults, will be a promising conceptual model that is complementary to existing equivalent continuum and discrete fracture models to analyze shearing of intersecting fractures/faults.