Search NASA⌕ Search

SEARCH · Search NASA

Results for “Aquifer 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.

At least 109 records · Page 6

Proxy-based Bayesian inversion of strain tensor data measured during well tests

Recent instrument developments have made it possible to measure the strain tensor caused by injecting or pumping fluid from aquifers or reservoirs, but the full value of these data is limited because the long runtimes of poroelastic forward models makes it impractical to use many inversion schemes. This limits the interpretation of strain data for managing the recovery of resources or storage of wastes in the subsurface. This paper describes a method of inverting deformation data using a poroelastic numerical simulator so the results can be used to manage reservoirs or aquifers. We developed a workflow designed to reduce the number of simulations sufficiently to make it feasible to use DREAMzs, an advanced Bayesian inversion method that translates the uncertainties from different sources into unbiased posterior parameter distributions and uncertainty envelopes around the field data. Using a KNN proxy model for the poroelastic simulator is key to reducing the overall computations, and the workflow includes a strategy for ensuring the proxy model results converge on the results from the simulator. The workflow is tested using an idealized example that verifies the ability to correctly identify parameters and characterize noise used to perturb the data. Field data from an injection test at an oil reservoir near Tulsa, Oklahoma, are also used to evaluate the efficacy of the workflow with a real dataset. The workflow identified 265 history matching solutions out of 1240 total simulation runs (21% acceptance ratio), where the results were used to characterize posterior parameter distribution and evaluate the prediction uncertainty. Furthermore, this workflow is significant because it enables strain tensor, or other geomechanical measurements to be interpreted to guide decision-making during energy and environmental processes in the subsurface.

42 ENGINEERING↗

Geologic Characterization of the South Georgia Rift Basin for Source Proximal CO2 Storage

The project Geologic Characterization of the South Georgia Rift Basin for Source Proximal CO2 Storage is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The South Carolina Research Foundation and partners evaluated the feasibility of CCS in the Jurassic/ Triassic (J / TR) saline formations of the buried Mesozoic South Georgia Rift (SGR) Basin that extends from South Carolina into Georgia. The J / TR sequence, based on preliminary assessment of limited geologic and geophysical data, appears to have both the appropriate areal extent and multiple horizons to permanently and safely store CO2 The presence of several igneous rock layers within the sequence may potentially provide adequate seals to prevent upward CO2 migration into the Coastal Plain aquifer systems. Approximately 81 kilometers of 2-D seismic reflection data were collected by Bay Geophysical, Inc. to explore a portion of the SGR located in southern Georgia. The 81 kilometers were divided into two lines approximately 40.5 kilometers each, with Line 1 intersecting Georgia well GGS 3457. Line 2 intersects Line 1 at the southern portion of Line 1 to maximize the extent of coverage away from GGS-3457 (a deep well drilled in the 1980s for oil and gas exploration). This well had a set of usable logs, including gamma and neutron logs that provided promising results related to CO2 storage. Results showed sandstone with porosity values greater than 10 percent and a thickness of 120 meters. The design of the seismic shot was to extrapolate information away from the well and to better define the extent of the SGR and the necessary reservoir and caprock for a successful CO2 injection. A numerical simulation model of CO2 Injection and migration was developed based on the geology log for the GGS-3457 well. The simulation model was used to investigate the feasibility of injecting 30 million metric tons of CO2 into SGR J / TA sediments and integrity of the diabase layers as seals to prevent CO2 migration.

2-D seismic↗

Subsurface microbial communities as a tool for characterizing regional-scale groundwater flow

Subsurface microbial community distribution patterns are influenced by biogeochemical and groundwater fluxes and may inform hydraulic connections along groundwater-flow paths. This study examined the regional-scale microbial community of the Death Valley Regional Flow System and evaluated whether subsurface communities can be used to identify groundwater-flow paths between recharge and discharge areas. Samples were collected from 36 sites in three groundwater basins: Pahute Mesa–Oasis Valley (PMOV), Ash Meadows (AM), and Alkali Flat–Furnace Creek Ranch (AFFCR). Microbial diversity within and between communities varied by location, and communities were separated into two overall groups that affiliated with the AM and PMOV/AFFCR basins. Network analysis revealed patterns between clusters of common microbes that represented groundwaters with similar geochemical conditions and largely corroborated hydraulic connections between recharge and discharge areas. Null model analyses identified deterministic and stochastic ecological processes contributing to microbial community assemblages. Most communities were more different than expected and governed by dispersal limitation, geochemical differences, or undominating processes. However, certain communities from sites located within or near the Nevada National Security Site were more similar than expected and dominated by homogeneous dispersal or selection. Overall, the (dis)similarities between the microbial communities of DVRFS recharge and discharge areas supported previously documented hydraulic connections between: (1) Spring Mountains and Ash Meadows; (2) Frenchman and Yucca Flat and Amargosa Desert; and (3) Amargosa Desert and Death Valley. However, only a portion of the flow path between Pahute Mesa and Oasis Valley could be supported by microbial community analyses, likely due to well-associated artifacts in samples from the two Oasis Valley sites. This study demonstrates the utility of combining microbial data with hydrologic, geologic, and water-chemistry information to comprehensively characterize groundwater systems, highlighting both strengths and limitations of this approach.

54 ENVIRONMENTAL SCIENCES↗

Incorporating Heterogeneity into 3D Geologic Models for the Alluvial Basin Reference Case

The unsaturated zone alluvium reference case is one of several geologic systems under consideration by the U.S. Department of Energy Office of Nuclear Energy for hosting repositories for spent nuclear fuel and associated waste (Sevougian et al., 2019). As noted by Mariner et al. (2018), the generic alluvial basin offers positive attributes that merit its consideration as a reference case by the Spent Fuel and Waste Science and Technology (SFWST) campaign. Hundreds of alluvial basins and sub-basins are scattered across the arid western United States (Figure 1-1). Precipitation and infiltration rates are relatively low with high evapotranspiration, resulting in deep regional water tables. Deep regional water tables allow for significant vertical separation between proposed repository depths and groundwater and thus longer transport paths from a hypothetical repository to an aquifer. Accumulations of alluvial sediments within these basins are typically on the order of hundreds of meters, and locally may exceed 1,000 m in thickness, as is the case for the Deming sub-basin in southern New Mexico. In alluvial basins, alluvium creates a thick geologic host medium that serves as the natural barrier system (NBS) in the conceptual model framework of a geologic disposal system and isolates the waste packages from receptors in the biosphere. Further, alluvial basin fill is typically comprised of stacked playa and lacustrine deposits along the basin axis (Perry et al., 2018). Characterized by low permeability, these layers protect the biosphere above the repository and the groundwater resources below the repository by limiting flow and transport of radionuclides.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs Results

Geothermal power plants typically show decreasing heat and power production rates over time. Mitigation strategies include optimizing the management of existing wells - increasing or decreasing the fluid flow rates across the wells - and drilling new wells at appropriate locations. The latter is expensive, time-consuming, and subject to many engineering constraints, but the former is a viable mechanism for periodic adjustment of the available fluid allocations. Data and supporting literature from a study describing a new approach combining reservoir modeling and machine learning to produce models that enable strategies for the mitigation of decreased heat and power production rates over time for geothermal power plants. The computational approach used enables translation of sets of potential flow rates for the active wells into reservoir-wide estimates of produced energy and discovery of optimal flow allocations among the studied sets. In our computational experiments, we utilize collections of simulations for a specific reservoir (which capture subsurface characterization and realize history matching) along with machine learning models that predict temperature and pressure timeseries for production wells. We evaluate this approach using an "open-source" reservoir we have constructed that captures many of the characteristics of Brady Hot Springs, a commercially operational geothermal field in Nevada, USA. Selected results from a reservoir model of Brady Hot Springs itself are presented to show successful application to an existing system. In both cases, energy predictions prove to be highly accurate: all observed prediction errors do not exceed 3.68% for temperatures and 4.75% for pressures. In a cumulative energy estimation, we observe prediction errors that are less than 4.04%. A typical reservoir simulation for Brady Hot Springs completes in approximately 4 hours, whereas our machine learning models yield accurate 20-year predictions for temperatures, pressures, and produced energy in 0.9 seconds. This paper aims to demonstrate how the models and techniques from our study can be applied to achieve rapid exploration of controlled parameters and optimization of other geothermal reservoirs. Includes a synthetic, yet realistic, model of a geothermal reservoir, referred to as open-source reservoir (OSR). OSR is a 10-well (4 injection wells and 6 production wells) system that resembles Brady Hot Springs (a commercially operational geothermal field in Nevada, USA) at a high level but has a number of sufficiently modified characteristics (which renders any possible similarity between specific characteristics like temperatures and pressures as purely random). We study OSR through CMG simulations with a wide range of flow allocation scenarios. Includes a dataset with 101 simulated scenarios that cover the period of time between 2020 and 2040 and a link to the published paper about this project, where we focus on the Machine Learning work for predicting OSR's energy production based on the simulation data, as well as a link to the GitHub repository where we have published the code we have developed (please refer to the repository's readme file to see instructions on how to run the code). Additional links are included to associated work led by the USGS to identify geologic factors associated with well productivity in geothermal fields. Below are the high-level steps for applying the same modeling + ML process to other geothermal reservoirs: 1. Develop a geologic model of the geothermal field. The location of faults, upflow zones, aquifers, etc. need to be accounted for as accurately as possible 2. The geologic model needs to be converted to a reservoir model that can be used in a reservoir simulator, such as, for instance, CMG STARS, TETRAD, or FALCON 3. Using native state modeling, the initial temperature and pressure distributions are evaluated, and they become the initial conditions for dynamic reservoir simulations 4....

15 GEOTHERMAL ENERGY↗

Self-Sealing Mafic Sills for Carbon and Hydrogen Storage

Tabular igneous intrusions (sills) are common features in sedimentary basins and have the potential to be useful seals for fluids (e.g., CO 2 or H 2 ) in geological storage scenarios, and may be important as the need for geologic carbon sequestration and alternative fuel storage increases. This is advantageous in regions without ready access to large-volume reservoirs in depleted hydrocarbon plays and saline aquifers, like the northeastern United States. Moreover, geological H 2 storage requires more demanding conditions than typically associated with oil and gas. An enhanced seal will have the properties of a traditional seal - competent, low permeability, and laterally extensive - with the addition of being able to self-seal pre-existing and induced fractures. Self-sealing will occur along fluid pathways like fractures where CO 2 , water, and minerals like plagioclase, olivine, and pyroxene react together. Dolerite sills from the Gettysburg Basin, Pennsylvania, have remarkably low permeability and homogeneous compositions that include minerals that will readily react with CO 2 dissolved in water. Here, we characterize the physical properties and chemical gradients within several mafic sills cored in five boreholes. In addition, preliminary CO 2 -reaction experiments on dolerite samples demonstrated rapid carbonate mineralization.

58 GEOSCIENCES↗

Geologic Mapping along the Arabia Terra Dichotomy Boundary: Mawrth Vallis and Nili Fossae, Mars: Introductory Report

Geologic mapping studies at the 1:1M-scale will be used to characterize geologic processes that have shaped the highlands along the Arabia Terra dichotomy boundary. In particular, this mapping will evaluate the distribution, stratigraphic position, and lateral continuity of compositionally distinct outcrops in Mawrth Vallis and Nili Fossae as identified by spectral instruments currently in orbit. Placing these landscapes, their material units, structural features, and unique compositional outcrops into spatial and temporal context with the remainder of the Arabia Terra dichotomy boundary will provide the ability to: 1) further test original dichotomy formation hypotheses, 2) constrain ancient paleoenvironments and climate conditions, and 3) evaluate various fluvial-nival modification processes related to past and present volatile distribution and their putative reservoirs (aquifers, lakes and oceans, surface and ground ice) and the influences of nearby volcanic and tectonic features on hydrologic processes in these regions. The result will be two 1:1M scale geologic maps of twelve MTM quadrangles (Mawrth Vallis - 20022, 20017, 20012, 25022, 25017, and 25012; and Nili Fossae - 20287, 20282, 25287, 25282, 30287, 30282).

Bleamaster, Leslie F., III↗

Groundwater Variability Across India, Under Contrasting Human and Natural Conditions

Characterizing local to regional scale water cycles and water resources will be crucial for achieving the United Nations' water-related Sustainable Developmental Goals. However, quantification and understanding of groundwater extraction across scales have been hampered by inadequate water usage reporting and limited information on irrigation practices. Here we analyze observations from ∼15,000 groundwater monitoring wells and the Gravity Recovery and Climate Experiment satellites together with irrigation, agricultural, and meteorological datasets to show how drought-induced coupling between natural and anthropogenic groundwater storage variations has caused sustainability challenges in India, the world's biggest consumer of groundwater for irrigation. Notably, the mechanisms and consequences of such coupling differ significantly depending on aquifer types. In Andhra Pradesh's hard rock aquifer, groundwater declines have been limited, despite the nearly constant water scarcity that its farmers face. Moreover, its free farm power policy involves an annual irrigation energy consumption of 26 billion kWh that costs US$ 2.5 billion, possibly unparalleled compared to any other part of the world of similar size (0.27 million km2). In West Bengal's highly permeable alluvial aquifer, the water table is declining rapidly (15 cm/yr) due to a policy that encourages irrigation. Situated between these two states, Odisha's aquifer shows substantial resilience to drought, owing to the state's relatively natural landscape and forest restoration policy. The findings of this study provide new insights to understand the divergent aspects of groundwater irrigation in north versus south India, which can enable development of adaptation and mitigation strategies to avert the looming water crisis.

Dileep K Panda↗

The Scientific Significance of Potential Samples From the Jezero Crater Rim

The Mars 2020 Perseverance rover has been exploring Jezero crater, Mars to characterize the geology,assess the potential for rocks to represent ancient habitable environments and/or preserve biosignatures, and collect a suite of scientifically compelling samples for return to Earth (Farleyet al.,2020).Thus far, Perseverance has collected samples that represent the diversity of igneous materials encountered on the crater floor and a suite of aqueously deposited sedimentary rocks from the fanand margin units(Simon et al., 2023; Bosak et al., 2024). When the rover ascendsthe crater rimit will encounter rocktypesthat are not included in the current sample cache. These include rocks fromthe Noachian whose mineralogical, geochemical, and geochronological recordswould greatly expand our understanding of early Mars. Some of the exposures are understood to come from subsurface aquifers that may represent the most ancient potentially habitablesubsurfaceenvironmentson Mars(Ehlmann et al., 2024). Samples ofthese rocks will add uniquescientific value tothe sample cache and to the Mars sample returnprogramby addressing outstanding questions in Mars science including potential habitability and how the climate and the interior of the planet evolved through time.

Mars 2020↗

Maximizing oil production from water alternating gas (CO 2 ) injection into residual oil zones: The impact of oil saturation and heterogeneity

Residual oil zones (ROZs) are widespread reservoirs, characterized by oil at residual saturation, either underlying oil fields (brownfield) or lateral (greenfield) to such fields. These reservoirs have the potential to produce volumes of oil sufficiently significant to make appreciable impacts on the US’s oil reserves and associated incidental CO 2 sequestration. The objective of this study is to improve our understanding the impact of heterogeneous and low oil saturations, in brownfield ROZs, on the effectiveness of water alternating gas (WAG) injection strategies. ROZs occur in the Permian Basin and elsewhere, and operators are using CO 2 injection for enhanced oil recovery (EOR) in these zones. The consensus model for the formation of ROZs is that they were formed by the effect of faster regional aquifer flow, acting over millions of years. Both the magnitude of oil saturation and the spatial distribution of oil differ from water-flooded main pay zones (MPZs). To explore the most effective injection strategies, we conducted simulations of CO 2 injection into synthetic geologic reservoirs. These simulations focused on injection into reservoirs subject to either man-made waterflooding or long-term natural waterflooding. By exploring the impact of varying: oil saturation; well patterns; reservoir heterogeneity; and permeability anisotropy, we attempt to quantify the factors that most influence the effectiveness of WAG injection. WAG ratios (the ratio of injected water and CO 2 , in reservoir volumes) of interest are those that either minimize the net CO 2 utilization ratios or maximize oil production rates. In general, the most effective WAG ratios for ROZs, are consistently less than those observed undergoing CO 2 injection in the same geologic reservoir models after traditional (man-made) waterflooding. This work demonstrates that most favorable WAG ratios for oil production in ROZs are different from those in traditional MPZs because of oil saturation differences. Therefore, CO 2 injection into both zones or directly copying WAG injection designs from MPZs to ROZs might not maximize oil production.

04 OIL SHALES AND TAR SANDS↗

The stable isotopes of natural waters at the Marcell Experimental Forest

Here, we present a new data set from the Marcell Experimental Forest (MEF) that compiles water isotope measurements from multiple research catchments, some of which have been studied since the 1960s. The MEF is located in northern Minnesota, USA, and is home to heavily studied and monitored forests, streams, bogs, and fens. Peat-forming systems (bogs and fens) are an important component of the MEF landscape and have a profound impact on the water cycle in these catchments. Within the last decade, analysis of stable isotopes of water (expressed as δD and δ 18 O) has been implemented to characterize the different components of the water budget, and to allow researchers to look at catchment and peatland-specific hydrologic effects in the watershed. This δD and δ 18 O data set of natural waters from MEF catchments is primarily composed of measurements from three peatlands (S1, S2, S6) during an 11-year period. More recently collection and analysis were expanded to also include samples from the Spruce and Peatland Responses Under Changing Environments (SPRUCE) project in the S1 bog, peatlands S3, S4, S5, as well as nearby lakes. We establish a local meteoric water line by analyzing the isotopic composition of precipitation, which fills a void in regional meteoric water lines for Minnesota. Furthermore, we establish baseline isotopic composition for bog outlet streams, bog porewater, aquifer groundwater, overland flow, subsurface stormflow, and snowpack, as well as runoff from the SPRUCE experimental chambers. These data are publicly available and will be expanded upon in the future.

54 ENVIRONMENTAL SCIENCES↗

Application of Modified Meshgraphnets for Subsurface Prediction during CO2 Sequestration

In the face of the increasingly dire consequences of anthropogenic climate change, capturing and storing carbon dioxide is paramount. However, several impediments exist to the safe and effective subsurface storage of CO2, such as cost of transport, identification of suitable sites for subsurface storage, and assessment of long-term risk from storage in subsurface aquifers. Accurate subsurface modeling is necessary to ensure that CO2 storage is both safe and effective. Still, such modeling has traditionally required either substantial time and computational power (numerical simulation) or a substantial amount of pre-existing data for training (machine learning models). Additionally, these models lack flexibility in dealing with both changes in discretization of the input data and generalizability beyond the data on which they are trained. In order to address these issues, this research applies graph neural networks (GNNs) to predict subsurface saturation and pressure during CO₂ injection in a model of the Illinois Basin-Decatur Project (IBDP). GNNs provide a flexible, intuitive method for representing and manipulating complex unstructured data, which is often found in many practical domain problems such as fluid flow and subsurface characterization. These unstructured grids are easily represented in GNNs by representing spatially-localized features such as permeability, porosity, saturation, and pressure as nodes in a graph and relationships between these properties as edges connecting these nodes. This research applies a specific GNN model called MeshGraphNets (MGN) to model the change in CO2 saturation and pressure over a 50-month time period (36 months of injection, 14 months post-injection). The MGN model leverages a message passing process that allows the network to learn both the spatial and temporal dynamics of this system simultaneously. Additionally, training on a limited dataset (64 realizations, 20 time points each) resulted in a high degree of accuracy in saturation prediction both within the same timeframe as the training (20 months, 0.039 average RMSE) and when projecting out to the end of injection (36 months, 0.053 average RMSE). Temporal predictions such as those generated by MGNs and other similar models are prone to accumulated error over time; in order to address this, a multi-step rollout (MSR) training process was applied to calculate training loss. This method mimics the forward prediction during inference by “rolling out” multiple time points in a single training step using the previous prediction as input to the MGN model. By calculating the loss several time steps forward from the current prediction, the model is forced to find a more stable state over time. Application of MSR to the MGN model resulted in an average 15% reduction in inference error over time during forward prediction. This study showcases the immense potential of GNNs as a game-changing methodology for predicting pressure and saturation evolution in CCS projects, ultimately paving the way for more sustainable and effective carbon storage solutions. Presentation prepared for the 2024 AiChE Annual Meeting, October 27 to November 1 2024, San Diego, CA.

Holcomb, Paul↗

Subsurface H 2 Storage: A Williston Basin Commercial-Scale Resource Study

Poster for the 2024 NETL Resource Sustainability Project Review Meeting, Pittsburgh, Pennsylvania, April 2-4, 2024. This poster presents a commercial‑scale assessment of subsurface hydrogen storage potential in the North Dakota portion of the Williston Basin, evaluating saline formations, depleted oil and gas reservoirs, and salt formations. The study integrates laboratory characterization, reservoir simulation, and basinwide analysis to assess storage capacity, injectivity, recovery, and risks related to geochemical, microbial, and wellbore interactions. Results support the feasibility of large‑volume, secure hydrogen storage and provide a framework to guide future hydrogen commercialization and infrastructure development.

08 HYDROGEN↗

Effect of Formation Heterogeneity on CO 2 Dissolution in Subsurface Porous Media

Dissolution trapping is one of the most dominant mechanisms for the secure storage of CO 2 injected in porous subsurface formations saturated with brine. This trapping mechanism is enhanced by convective mixing, which occurs due to the gravitational instability between the different fluid layers in the aquifer. The reservoir permeability also plays a crucial role in the dissolution rate and overall fluid flow dynamics during the density-driven convection in porous media. Herein this study investigates the role of complex heterogeneity, i.e., irregular permeability distribution in CO 2 dissolution, using a novel experimental approach to create medium permeability heterogeneity in Hele–Shaw cells. Complex subsurface transport phenomena such as a preferential dissolution path, CO 2 sweep efficiency, changes in finger morphology, and CO 2 concentration distribution are visualized by creating heterogeneous media. Experimental results showed that reservoir permeability heterogeneity causes significant channeling effects and poor sweep efficiency. A scaling relationship between average finger growth rate (Gr) and permeability (k) was obtained as Gr [m s –1 ] = 266.8k [m 2 ] + 1.20 × 10 –6 . Furthermore, the mass of CO 2 dissolved is calculated using the spectrophotometric method to characterize the convective instability. The convective flux was analyzed by comparing the experimental dissolution flux with the theoretical diffusion flux, calculating a maximum Sherwood number of 6.8. The study’s findings improve the current understanding of the CO 2 convection morphology in heterogeneous media, allowing better assessment of long-term CO 2 storage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Kilometer-Scale Hydrologic Exchange Flows in a Gravel Bed River Corridor and Their Implications to Solute Migration

In this work, a well-characterized field site along a major river corridor was used to investigate the dynamic pathways and impacts of subsurface hydrogeologic structures on kilometer-scale hydrologic exchange flows between river water and groundwater. A aqueous uranium (U aq ) plume exists at the site that discharges to the Columbia River. We performed temporally intensive monitoring of specific conductance (SpC) and U aq concentrations within the plume at increasing distances from the river shoreline and at locations both within and outside a presumed subsurface pathway of lateral hydrologic exchange. Specific conductance and U aq were utilized as in-situ tracers of hydrologic exchange and associated groundwater-surface water mixing. Monitoring data was collected for an approximate two year period during which seasonal river elevations varied by more than 2 m causing distinct events of river water intrusion and retreat from the near-shore aquifer, and large temporal variations in well-water analyte concentrations. Statistical clustering analysis was applied to identify well groupings exhibiting common time series behaviors. PFLOTRAN was implemented within a modeling domain of 3.3 x 1.9 km with dynamic flow boundary conditions and accurate subsurface topography for the basal aquitard underlying the aquifer. In-silico calculations of river water tracer behavior were performed to assess predominant flow directions and velocities influenced by river stage variations and variable aquitard topography, and their collective influences on SpC and U aq concentrations. We demonstrate complicated multidirectional flow behaviors at the kilometer scale resulting from hydrologic exchange and subsurface structure that influence plume dynamics and determine the overall system behavior.

54 ENVIRONMENTAL SCIENCES↗

Risk-based area of review estimation in overpressured reservoirs to support injection well storage facility permit requirements for CO 2 storage projects

This paper by the Energy & Environmental Research Center presents a workflow and modeling approach for delineating a risk-based area of review (AOR) to support a U.S. Environmental Protection Agency (EPA) Class VI permit for a carbon dioxide (CO 2 ) storage project. The approach combines semianalytical solutions for estimating formation fluid leakage through a hypothetical leaky wellbore with the results of numerical reservoir simulations to define the AOR. The modeling utilizes 1) semianalytical solutions from the peer-reviewed literature for formation fluid leakage through abandoned wellbores by Raven (1990) and Avci (1994), 2) a FORTRAN model compiled and described in Cihan et al. (2011, 2012) called ASLMA (Analytical Solution for Leakage in Multilayered Aquifers), and 3) a computational framework for estimating a risk-based AOR first proposed by Oldenburg et al. (2014, 2016). Therefore, the approach builds upon well-established research and underlying hydrogeological principles that have been upheld for nearly three decades. Moreover, the ASLMA model has been broadly applied to an array of storage projects. The work presented herein extends these earlier works using a custom wrapper written in the software environment, R (R Core Team, 2020), which was developed to perform multiple runs of the ASLMA model using given ranges for one or more input parameters. In addition, the current work simulates the pressure buildup within the storage reservoir in response to CO 2 injection using a compositional simulator to better accommodate the temporospatial evolution of pressure buildup within the storage reservoir that is more accurately modeled using a heterogeneous geologic model and a compositional simulator that accounts for the multiphase interactions. The workflow is demonstrated using a case study for a 180,000-metric-ton-per-year storage project located in the PCOR (Plains CO 2 Reduction) Partnership region. For the storage project evaluated here, under the scenario where the leaky wellbore is open to a saline aquifer (thief zone) between the overlying seal (cap rock) and the underground sources of drinking water (USDW), the risk-based AOR essentially collapses to the areal extent of the CO 2 plume in the storage reservoir because the pressure buildup in the storage reservoir beyond the CO 2 plume is insufficient to drive formation fluids up a hypothetical leaky wellbore into the USDW. However, even under the conservative assumption that the leaky wellbore is not open to a thief zone, beyond the areal extent of the CO 2 plume, the incremental leakage is less than 400 m 3 over 20 years, which represents ~0.0001% or less of the total volume of water contained within the USDW rock volume. As discussed in the text, the threshold criterion for defining the risk-based AOR is site-specific and should be informed by the results of the sensitivity analysis and available site characterization data. The approach outlined in this paper is designed to be protective of USDWs and, therefore, comply with the Safe Drinking Water Act requirements and provisions for the U.S. EPA Class VI Underground Injection Control (UIC) Program (Class VI Rule) and North Dakota Administrative Code Chapter 43-05-01.

54 ENVIRONMENTAL SCIENCES↗

Sedimentary Geothermal Resources in Nevada, Utah, Colorado, and Texas

The objectives of this project were to (1) perform a literature review of sedimentary geothermal resources, (2) identify data sources and develop data-collection methodologies that characterize selected resources, (3) screen sedimentary basins and formations for sedimentary geothermal potential, and (4) evaluate the technical feasibility of one or more selected locations. Numerous publications have characterized geothermal resources within sedimentary basins. A literature search reviewed publications describing resources located in Colorado, Louisiana, Nevada, Texas, Utah, and Wyoming. The most attractive resources have high temperature gradients, low drilling costs, and reservoir permeabilities greater than 10 millidarcies (mD). Prospects in Colorado, Nevada, Texas, and Utah exhibit attractive characteristics and were chosen for further analysis. Sedimentary resources in Nevada and Utah are most attractive, followed by tested resources in Texas and untested resources in Colorado. The identified resources in Wyoming and Louisiana had lower geothermal gradients and were not evaluated. Reservoir modeling and techno-economic analysis were performed at Marys River Basin–North in Nevada. Geothermal energy production at this location is expected to have a levelized cost of energy (LCOE) ranging between 10 and 20 cents/kWh. Additional work may result in lower LCOE estimates at this location and at other attractive prospects in these three regions. The Great Basin carbonate and alluvial aquifer system of eastern Nevada and western Utah includes a lower carbonate aquifer unit, which has the potential for hosting both conduction- and convection-dominated geothermal systems. Mapping has identified lateral thickness variability expressed as a roughly 120–150-km-wide central corridor, which hosts the thickest and most continuous formations and extends from near Las Vegas north to the Idaho border. Purchase and analysis of privately held legacy seismic data could potentially compensate for the lack of sufficient data documenting measured depth to the lower carbonate aquifer unit. Multiple orogenies, extension episodes, and intrusive events have deformed and displaced the target formations of the lower carbonate aquifer unit. This structural complexity and potential for dipping reservoirs emphasizes the need for detailed geological, geophysical, and reservoir modeling. Heat flow within three Colorado sedimentary basins reviewed as part of this study was calculated in targeted studies by the Colorado Geologic Survey and Colorado School of Mines. These calculations are based on bottom-hole temperature data sets with significant limitations and some variability but produce values consistently higher than the global continental average of 65 mW/m2 for all three basins. Heat flow in the Raton Basin is the highest; however, permeability measurements from specific sedimentary formations with high heat flow have not been obtained. Promising formations for sedimentary geothermal systems were found in all three regions studied—Nevada-Utah, Colorado, and Texas. The next steps for developing sedimentary geothermal resources vary due to differences in available data and resource uncertainties. Additional scopes of work are recommended for identified basins in these three regions.

15 GEOTHERMAL ENERGY↗

Temporal flow variations interact with spatial physical heterogeneity to impact solute transport in managed river corridors

The interactions between surface water and groundwater in river corridors lead to temporal fluctuations in subsurface water fluxes which have a critical role on solute transport dynamics. Here, we develop a framework to analyze the relative impacts of different temporal frequencies of the flow field in a spatially heterogeneous aquifer on solute transport. Our analysis indicates that the advection-dispersion equation behaves as a low-pass filter by wiping out the effect of high-frequency velocity fluctuations on the first two spatial moments of the solute plume, namely its center of mass and spreading. The concepts discussed in the theoretical analysis are then applied to understand solute transport dynamics at the 300 Area of the Hanford site (USA) adjacent to the Columbia River. We examine the temporal behavior of the solute plume's spatial moments for different temporal frequencies utilizing geostatistical parameters estimated in the 300 Area. Due to the proximity to the Columbia river, groundwater fluxes at the Hanford site are highly dynamic resulting in a large range of characteristic temporal frequencies. Nonetheless, similar to the theoretical analysis, our results show that the effect of high-frequency fluctuations is filtered, with most of the solute transport dynamics being controlled by fluctuations characterized by a large characteristic period.

54 ENVIRONMENTAL SCIENCES↗