SEARCH · Search NASA

Results for “subsurface flow”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Deep Learning for Simultaneous Inference of Hydraulic and Transport Properties

Abstract Identification of a heterogeneous conductivity field and reconstruction of a contaminant release history are key aspects of subsurface remediation. These two goals are achieved by combining model predictions with sparse and noisy hydraulic head and concentration measurements. Solution of this inverse problem is notoriously difficult due to, in part, high dimensionality of the parameter space and high computational cost of repeated forward solves. We use a convolutional adversarial autoencoder (CAAE) to parameterize a heterogeneous non‐Gaussian conductivity field via a low‐dimensional latent representation. A three‐dimensional dense convolutional encoder‐decoder (DenseED) network serves as a forward surrogate of the flow and transport model. The CAAE‐DenseED surrogate is fed into the ensemble smoother with multiple data assimilation (ESMDA) algorithm to sample from the Bayesian posterior distribution of the unknown parameters, forming a CAAE‐DenseED‐ESMDA inversion framework. The resulting CAAE‐DenseED‐ESMDA inversion strategy is used to identify a three‐dimensional contaminant source and conductivity field. A comparison of the inversion results from CAAE‐ESMDA with physical flow and transport simulator and from CAAE‐DenseED‐ESMDA shows that the latter yields accurate reconstruction results at the fraction of the computational cost of the former.

Zhou, Zitong↗

DOE OSTI2019

Removal of Historic Low-Level Radioactive Waste from the Highland Drive South and Alexander Street Ravines - 19155

In Port Hope, the Highland Drive South Ravine (HDSR) is the discharge area for groundwater originating from the closed Highland Drive Landfill site. The presence of low-level radioactive waste (LLRW) in the HDSR is the result of a contaminated groundwater plume flowing from the landfill site. The contaminant plume contains elevated concentrations of uranium and arsenic and extends from the landfill south to discharge groundwater to shallow soils in a wet discharge area within the ravine. The groundwater discharges to Hunt's Ponds and Highland Drive South Creek, which are immediately downstream of the wet discharge area and will be included in the remediation design. Remediation and environmental management plans for HDSR have been developed within the framework of the Port Hope Project. Recent subsurface investigations have identified LLRW at greater depths than previously anticipated. Other investigations of surrounding properties have identified LLRW that is primarily from filling activities in backyards. The steep topography of the site and the presence of the creek and seepage face present remediation challenges. Also in Port Hope, the LLRW in the Alexander Street Ravine Area is the result of end dumping, filling of private backyards and stack deposition from the former refinery. The original site extent was focused on the deeper contamination due to end dumping and included a portion of Alexander Creek. However, recent investigations have confirmed other fill areas and the extent of stack deposition that are generally contiguous with the end dumping area. Restricted access to the site will result in all remediation in this area taking place via a single access point along the narrow existing Hayward Street. While both of these sites are ravines in Port Hope containing LLRW, the origins of contamination and remediation approaches are very different. This paper will compare these two challenging remediation sites. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

The value of in-reservoir energy storage for flexible dispatch of geothermal power

Geothermal systems making use of advanced drilling and well stimulation techniques have the potential to provide tens to hundreds of gigawatts of clean electricity generation in the United States by 2050. With near-zero variable costs, geothermal plants have traditionally been envisioned as providing “baseload” power, generating at their maximum rated output at all times. However, as variable renewable energy sources (VREs) see greater deployment in energy markets, baseload power is becoming increasingly less competitive relative to flexible, dispatchable generation and energy storage. Herein we conduct an analysis of the potential for future geothermal plants to provide both of these services, taking advantage of the natural properties of confined, engineered geothermal reservoirs to store energy in the form of accumulated, pressurized geofluid and provide flexible load-following generation. We develop a linear optimization model based on multi-physics reservoir simulations that captures the transient pressure and flow behaviors within a confined, engineered geothermal reservoir. We then optimize the investment decisions and hourly operations of a power plant exploiting such a reservoir against a set of historical and modeled future electricity price series. Further, we find that operational flexibility and in-reservoir energy storage can significantly enhance the value of geothermal plants in markets with high VRE penetration, with energy value improvements of up to 60% relative to conventional baseload plants operating under identical conditions. Across a range of realistic subsurface and operational conditions, our modeling demonstrates that confined, engineered geothermal reservoirs can provide large and effectively free energy storage capacity, with round-trip storage efficiencies comparable to those of leading grid-scale energy storage technologies. Optimized operational strategies indicate that flexible geothermal plants can provide both short- and long-duration energy storage, prioritizing output during periods of high electricity prices. Sensitivity analysis assesses the variation in outcomes across a range of subsurface conditions and cost scenarios.

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

Techno-economic feasibility of geothermal energy production using inactive oil and gas wells for district heating and cooling systems in Tuttle, Oklahoma

Oil and gas wells have been drilled in the United States and often encounter hot co-produced water possibly suitable for geothermal direct-use applications. Here this study evaluated geothermal resources, heating and cooling demand, and techno-economic potential of four oil and gas wells repurposed for geothermal energy production to serve two public schools and 250 nearby houses in Tuttle, Oklahoma, USA. Subsurface geology in the study area primarily consists of Permian to Mississippian sandstone and limestone formations approximately from 2 km to 3.5 km depth. With a geothermal gradient of 25.3 °C/km, the temperature of geothermal energy production was estimated as 68.2 °C at 2.1 km. The potential of calcite scaling within boreholes and pipes was minimal. Using the characterized reservoir properties, techno-economic analyses were conducted for evaluating levelized costs of geothermal energy production for direct-use heating and cooling and district heating systems with five different production scenarios. Heating and cooling demand in two schools and 250 houses and additional 0.3 MW and 0.6 MW base heating loads for potential geothermal direct-use applications in the study area were also incorporated in the modeling. Results indicated that the levelized cost of heating of geothermal district heating systems utilizing existing wells were significantly lower than those of conventional geothermal energy systems that include well drillings. The geothermal energy production with doublet and quartet configurations was generally sufficient to supply total annual heating demand in the district, while a peaking boiler was used to supply peak loads. Doublet configurations showed higher production temperature with a lower thermal drawdown than the quartet configurations. The doublet system at 2.1 km production depth with 9.3 kg/s flow rate represented the optimal levelized cost of $71/MWh and 91.2% utilization factor. These results imply that the Tuttle geothermal district heating system is techno-economically feasible and comparable to geothermal district heating systems in the United States and the regional natural gas rate for residentials by recycling inactive oil and gas wells with no drillings.

15 GEOTHERMAL ENERGY↗

DOE OSTI2019

Spatially tracking carbon through the root-rhizosphere-soil system using laser ablation-IRMS

The intimate relationships between plant roots, rhizosphere, and soil are fostered by the release of organic compounds from the plant into soil through various forms of rhizodeposition and the simultaneous harvestingof inorganic nutrients from the soil to the plant. Here we present a method to spatially track and map the migration of plant-derived carbon (C) through roots into the rhizosphere and surrounding soil using laser ablation-isotope ratio mass spectrometry (LA-IRMS). Weused switchgrass microcosms containing soil from field plots at the Kellogg Biological Station (Hickory Corners, Michigan, USA), which have been cropped with switchgrass since 2008. We used a 13 CO 2 tracerto isotopically label switchgrass plants for two diel cycles and tracked subsequent movement of labeled C using the spatially specific (<100 µm resolution) δ 13 C analysis enabled by LA-IRMS.This approach permitted assessment of variable C flow through different roots and enabled mapping of spatial variability of C allocation to the rhizosphere. Highly 13 C-enriched C (consistent with production during the 13 CO 2 application period) extended ~0.5 – 1 mm from the root into the soil, suggesting that the majority of recent plant-derived C was within this distance of the root after 48 hours. Tracking the physical extent of root exudation into the rhizosphere can help evaluate the localization of plant-microbe interactions in highly variable subsurface environments,and the use of the isotopic label can differentiate freshly fixed C (presumably from root exudates) from other types of subsurface C (e.g., plant necromass and microbial turnover). The LA-IRMS technique mayalso serve as a valuable screening technique to identify areas of high activity for additional microbial or geochemicalassays.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Effect of CO 2 -brine-rock reactions on pore architecture and permeability in dolostone: Implications for CO 2 storage and EOR

Geologic carbon sequestration (GCS) is considered a feasible technology for storing substantive volumes of greenhouse gases in subsurface geological formations. In the reservoir, far from carbon dioxide (CO 2 ) injection wells or in post-injection scenarios, diffusion dominates over advection. This condition conjoins with spatially distributed geochemical reactions to induce heterogeneous changes in pore architecture, i.e. pore body and throat sizes or surface roughness. These changes can affect CO 2 transport properties and storage capacity. In this work, we investigated mineral dissolution and precipitation in dolomite samples saturated with a CO 2 -saturated brine at 93 °C and 34.5 MPa, aged without flow. Two rock types samples, i.e. intergranular- and vuggy-dominant, were selected to investigate changes in pore size, porosity and permeability under reactive conditions. Mineral dissolution and precipitation were characterized using scanning electron microscopy. Changes in pore size were quantified via time-domain nuclear magnetic resonance (TD-NMR) transverse relaxation time (T 2 ) and diffusion coefficient (D) distributions. We show that mineral dissolution likely occurs in highly permeable pathways. These observations are confirmed through analysis of (T 2 ) and diffusion coefficient (D) distributions. In contrast to results during CO 2 -enriched brine continuous injection, mineral precipitation was observed in micropores. The leftward shift of the T 2 peaks, corresponding to micropores, also evidenced mineral precipitation in lowpermeability zones. However, microscale alterations resulted only in a subtle increase in porosity and permeability. Results in this study shed light on effects of geochemical reactions on alteration of rock properties in diffusion-dominated regions during CO 2 storage.

58 GEOSCIENCES↗

DOE OSTI2023

Developing Methods to Assess Changes in Mechanical Properties of Shale Modified by Engineered Mineral Precipitation

Fractures in subsurface shale formations serve multiple purposes, for example, in the recovery of resources in hydraulic fracturing or as potential harmful leakage passages through caprocks that may contribute undesired fluids to the atmosphere or functional groundwater aquifers. A proposed method to seal or influence fracture properties is Ureolysis-Induced Calcium Carbonate Precipitation (UICP), a bio-mineralization technology driven by the enzymatic hydrolysis of urea, resulting in the formation of calcium carbonate. The resulting calcium carbonate can bridge the gaps in fractured shale and reduce fluid flow through fractures. This study represents the first step toward determining the influence of UICP treatment on shale material and its subsequent mechanical strength properties. The goal of this preliminary work is twofold: first, we aim to identify a method to test tensile strength along a core axis and second, we seek to assess the effect of temperature on the tensile strength of intact, unfractured shale cores (2.54 cm (1 in) diameter, 5.08 cm (2 in) long for comparison with future fractured and UICP-treated cores. A modified Brazilian indirect tensile strength test successfully measured splitting tensile strength of shale cores from Eagle Ford and Wolfcamp formations at room temperature and 60°C.

clastic rock↗

DOE OSTI2022

Conservative solute transport processes and associated transient storage mechanisms: Comparing streams with contrasting channel morphologies, land use and land cover

Abstract Land use within a watershed impacts stream channel morphology and hydrology and, therefore, in‐stream solute transport processes and associated transient storage mechanisms. This study evaluated transport processes in two contrasting stream sites where channel morphology was influenced by the surrounding land use, land cover, climate and geologic controls: Como Creek, CO, a relatively undisturbed, high gradient, forested stream with a gravel bed and complex channel morphology, and Clear Creek, IA, an incised, low‐gradient stream with low‐permeability substrate draining an agricultural landscape. We performed conservative stream tracer injections at these sites to address the following questions: (1) How does solute transport vary between streams with differing morphologies? and (2) How does solute transport at each stream site change as a function of discharge? We analysed in‐stream tracer time series data and compared results quantifying solute attenuation in surface and subsurface transient storage zones. Significant trends were observed in these metrics with varying discharge conditions at the forested site but not at the agricultural site. There was a broad range of transport mechanisms and evidence of substantial exchange with both surface and hyporheic transient storage in the relatively undisturbed, forested stream. Changing discharge conditions activated or deactivated different solute transport mechanisms in the forested site and greatly impacted advective travel time. Conversely in the simplified agricultural stream, there was a narrow range of solute transport behaviour across flows and predominantly surface transient storage at all measured discharge conditions. These results demonstrate how channel simplification inhibits available solute transport mechanisms across varying discharge conditions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Plutonium reactive transport in fractured granite: Multi-species experiments and simulations

Plutonium (Pu) in the subsurface environment can transport in different oxidation states as an aqueous solute or as colloidal particles. The transport behavior of Pu is affected by the relative abundances of these species and can be difficult to predict when they simultaneously exist. This study investigates the concurrent transport of Pu intrinsic colloids, Pu(IV) (aq) and Pu(V-VI) (aq) through a combination of controlled experiments and semi-analytical dual-porosity transport modeling. Pu transport experiments were conducted in a fractured granite at high and low flow rates to elucidate sorption processes and their scaling behavior. In the experiments, Pu(IV) (aq) was the least mobile of the Pu species, Pu(V-VI) (aq) had intermediate mobility, and the colloidal Pu, which consisted mainly of precipitated and/or hydrolyzed Pu(IV), was the most mobile. The semi-analytical modeling revealed that the sorption of each Pu species was rate-limited, as the sorption could not be described by assuming local equilibrium in the experiments. The model was able to describe the sorption of the different Pu species that occurring either on fracture surfaces, in the pores of the rock matrix, or simultaneously in both locations. While equally good fits to the data could be achieved using any of these assumptions, a fracture-dominated process was considered to be the most plausible because it provided the most reasonable estimates of sorption rate constants. Importantly, a key result of this work is that the sorption rate constant of all Pu species tends to decrease with increasing time scales, which implies that Pu will tend to be more mobile at longer time scales than observations at shorter time scales suggest. This result has important implications for predicting the environmental impacts of Pu in the safety assessments of geologic repositories for radioactive waste disposal, and we explore potential mechanistic bases for upscaling the sorption rate constants to time and distance scales that cannot be practically evaluated in experiments.

58 GEOSCIENCES↗

DOE OSTI2024

Dynamic Surface Incorporation of Pb 2+ Ions at the Actively Dissolving Calcite (104) Surface

The reaction of dissolved Pb 2+ with calcite surfaces at near-equilibrium conditions involves adsorption of Pb 2+ and precipitation of secondary heteroepitaxial Pb-carbonate minerals. A more complex behavior is observed under far-from-equilibrium conditions, including strong inhibition of calcite dissolution, development of microtopography, and near-surface incorporation of multiple monolayers (ML) of Pb 2+ without precipitation of secondary phases [where 1 ML ≡ 1 Ca/20.2 Å 2 , the crystallographic site density of the calcite (104) lattice plane]. However, the mechanistic controls governing far-from-equilibrium reactivity are not well understood. Here, in this study, we observe the interfacial incorporation of dissolved Pb 2+ during the dissolution of calcite (104) surfaces at pH ~3.7 in a flow-through reaction cell, revealing the formation of a ~1 nm thick Pb-rich calcite layer with a total Pb coverage of ~1.4 ML. These observations of the sorbed Pb distribution used resonant anomalous X-ray reflectivity, X-ray fluorescence, and nanoinfrared atomic force microscopy. We propose that this altered surface layer represents a novel sorption mode that is stabilized by conditions of sustained disequilibrium. This behavior may significantly impact the transport of dissolved metals during disequilibrium processes occurring in acid mine drainage and subsurface CO 2 injection and, if appropriately accounted for, could improve the predictive capability of geochemical reactive-transport models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Transient Pressure Interference During CO2 Injection in Saline Aquifers

Abstract CO2 injection in subsurface geological formations (e.g. deep saline aquifers) causes pressure perturbations over a large area surrounding the injection well. Observation wells are widely considered in geologic CO2 storage (GCS) projects where the pressure perturbation induced by CO2 injection is measured. In this work, we use analytical and numerical modeling tools along with field data to examine the pressure behavior in GCS projects before and after CO2 arrival at an observation well. Prior to CO2 arrival, a baseline pressure trend is established which corresponds to single-phase brine flow across the observation well (approximated by Theis solution). Therefore, analysis of early-time pressure data is straightforward, provides the single-phase flow characteristics (mobility and storativity), and helps establishing a baseline pressure change that can be extended beyond the single-phase flow period at the observation well. Upon CO2 arrival, a departure from this baseline trend is expected. For the pressure to detect the CO2 arrival at an observation well, the departure from baseline pressure behavior must be significant and well above the background noise levels. We use existing analytical models to determine the strength of the expected pressure departure signal from the baseline trend upon CO2 arrival. The strength of the expected pressure departure is found to be directly proportional to the mobility ratio. Accordingly, we establish a criterion to determine whether the pressure at an observation well can detect the CO2 arrival. We present an analysis approach through application to synthetic and field data and show the characteristic pressure behavior before and after CO2 arrival. We show that while generally the pressure can be either above or below the expected baseline pressure trend, it would be likely above the baseline upon CO2 arrival. This is because the mobility ratio becomes less than unity after CO2 arrival. We show that depending on the reservoir characteristics, changes in the pressure trend may or may not be sufficient to detect the CO2 arrival.

Engineering↗

NASA NTRS2010

Tectonics of the Outer Planet Satellites

Tectonic features on the satellites of the outer planets range from the familiar, such as clearly recognizable graben on many satellites, to the bizarre, such as the ubiquitous double ridges on Europa, the twisting sets of ridges on Triton, or the isolated giant mountains rising from Io's surface. All of the large and middle-sized outer planet satellites except Io are dominated by water ice near their surfaces. Though ice is a brittle material at the cold temperatures found in the outer solar system, the amount of energy it takes to bring it close to its melting point is lower than for a rocky body. Therefore, some unique features of icy satellite tectonics may be influenced by a near-surface ductile layer beneath the brittle surface material, and several of the icy satellites may possess subsurface oceans. Sources of stress to drive tectonism are commonly dominated by the tides that deform these satellites as they orbit their primary giant planets. On several satellites, the observed tectonic features may be the result of changes in their tidal figures, or motions of their solid surfaces with respect to their tidal figures. Other driving mechanisms for tectonics include volume changes due to ice or water phase changes in the interior, thermoelastic stress, deformation of the surface above rising diapirs of warm ice, and motion of subsurface material toward large impact basins as they fill in and relax. Most satellites exhibit evidence for extensional deformation, and some exhibit strike-slip faulting, whereas contractional tectonism appears to be rare. Io s surface is unique, exhibiting huge isolated mountains that may be blocks of crust tilting and foundering into the rapidly emptying interior as the surface is constantly buried by deposits from hyperactive volcanoes. Of the satellites, diminutive Enceladus is spectacularly active; its south polar terrain is a site of young tectonism, copious heat flow, and tall plumes venting into space. Europa's surface is pervasively tectonized, covered with a diverse array of exotic and incompletely understood tectonic features. The paucity of impact craters on Europa suggests that its tectonic activity is ongoing. Geysers on Triton show that some degree of current activity, while tectonic features that cross sparsely cratered terrain indicate that it may also be tectonically active. Ganymede and Miranda both exhibit ancient terrains that have been pulled apart by normal faulting. On Ganymede these faults form a global network, while they are confined to regional provinces on Miranda. Ariel, Dione, Tethys, Rhea, and Titania all have systems of normal faults cutting across their surfaces, though the rifting is less pronounced than it is on Ganymede and Miranda. Iapetus exhibits a giant equatorial ridge that has defied simple explanation. The rest of the large and middle-sized satellites show very little evidence for tectonic features on their surfaces, though the exploration of Titan's surface has just begun.

McKinnon, W. B.↗

DOE OSTI2021

Origin of water masses in Floridan Aquifer System revealed by 81 Kr

Groundwater in coastal aquifers serves as an essential resource in densely populated areas throughout the world. However, hydrological connection to the ocean leaves it at risk of salinization via changes in climate and hydrological cycles. Therefore, an accurate hydrological characterization of coastal aquifers is of the utmost importance. Although localized salinization of shallow aquifers has commonly been studied, geochemical constraints on the time scale of regional freshwater and seawater in coastal aquifers are limited especially in carbonate aquifers. Such information is not only crucial for water resource management, but also for understanding solute flux over the land-ocean boundary and for interpreting the geochemical signatures that each water mass carries. Toward this goal, we report tracer-based, subsurface residence times of groundwater in the Floridan Aquifer System in southern Florida, which is one of the most productive aquifers on Earth. The first application of Kr 81-an emerging tracer of old groundwater-to this aquifer has identified freshwater recharge during the last glacial period as anticipated from previous studies using other geochemical tools including radiocarbon, stable isotopes and noble gas concentrations. Moreover, freshwater of Holocene age was also detected at 100 km-distance from the recharge area, suggesting the possibility of an active flow system in the upgradient region. A contribution of fossil seawater that predates the last glacial maximum was locally identified, suggesting slow seawater circulation and limited but significant solute flux from the ocean into the aquifer for dolomitization. This also implies that the deep saline coastal groundwater can potentially serve as a paleo-seawater archive.

58 GEOSCIENCES↗

NASA NTRSDate not supplied

Imdr Regio: The Geology of A Possible Active Hot Spot on Venus and A Target for Future Exploration By the ESA Envision Mission.

With Venus being the focus of future new missions like ESA’s EnVision, we need to look for target sites for future investigations on Venus. The different instruments will help to study the planet in a unprecedented detail. In particular, the VenSAR instrument will map the 20% of the surface of Venus at a resolution of 5 to 10 meters/pixel, so it gets crucial to select what areas of Venus deserve to be targeted for such a high-resolution mapping. We propose here Imdr Regio as a good candidate location for combined high resolution SAR mapping and for study with the rest of the mission instruments (e.g., Subsurface Radar Sounder). Imdr Regio is a large topographic rise that extends approximately from 35ºS-50ºS and 195ºE-225ºE between the volcanic plains of Helen, Nsomeka and Wawalag Planitiae. Imdr Regio was classified as a volcano-dominated topographic rise with a minimum-maximum diameter of 1200-1400 km and a swell height of 1.6 km. The southeast of Imdr Regio is dominated by Idunn Mons (46.5ºS/214.5ºE). Studies using the 1 µm (derived) surface emissivity from the Venus Express mission of the volcanic flows surrounding Idunn Mons suggest that high emissivity values in the volcano flanks are related to relatively unweathered basaltic rock and therefore indicative of a recent or even ongoing volcanic activity. We have carried out geologic mapping in the area to constrain volcanic and tectonic structures, and the geologic history of the large topographic rise. This geologic mapping reveals that different styles of volcanism are present across Imdr Regio and that this volcanic activity takes place in all the hot spot in close relation with the formation of a rift throughout all time represented by in the area. The first stages on the evolution of this large volcano are characterized by the formation of a radial fracture system and contemporaneous large sheet flows. These large sheet flows are locally difficult to distinguish from regional plains and are also deformed by regional N-trending wrinkle ridges. After this initial phase, multiple overlapping digitate flow units form the flanks and summit of the volcano and are contemporaneous with NW-SE trending fractures and graben of Olapa Chasma, a rift system that cross the topographic rise. To the northwest another large volcano also presents a system of radial fractures but lacks clear large sheet flows. Numerous pit chains are associated to fractures and graben, suggesting that transport of magma under the surface is important. The geologic history based on the mapping suggests that activity in Imdr Regio started with the emplacement of a plume/diapir in the southeast, which resulted in the formation of a radial dyke system and the emplacement of large sheet flows in Idunn Mons. The presence of other large volcanoes to the northwest suggests the presence of another magmatic source (plume or diapir). After this initial stage, volcanism continues in Idunn Mons and activity in the topographic rise is strongly related to the formation of Olapa Chasma. Rift-related volcanic flows postdate these other volcanoes and are contemporaneous with the late activity in Idunn Mons. Monogenetic volcanism is present in the rift and in all the units that are in the area.

I. López↗

DOE OSTI2021

Injection data analysis using material balance time for CO 2 storage capacity estimation in deep closed saline aquifers

Estimating the ultimate storage capacity of deep saline aquifers is important to address the formation potential to store the envisioned large volumes of CO 2 . Injection data (i.e. injection rate, bottomhole pressure, and cumulative injected volume of CO 2 ) are routinely recorded during storage operations. These data contain valuable information on the subsurface (e.g. the reservoir pore volume and the formation storage capacity) that can be extracted. In this paper, we present a two-step graphical technique to infer the pore volume and the ultimate storage capacity of closed saline aquifers by analyzing the available injection data. First, the pore volume is inferred through adapting the concept of the material balance time. Material balance time is an approximate superposition time function developed to interpret production data from oil and gas wells operating at variable pressure/rate conditions during the boundary-dominated flow period. Using material balance techniques, the ultimate storage capacity is then estimated through linear extrapolation of the average pressure trend to the maximum allowable pressure the formation can withstand. The average pressure is not available in practice, but is can be obtained from the injection data. Two approaches are presented in this study to calculate the average pressure; namely the rigorous and the approximate approaches. Unlike the rigorous approach, the approximate approach does not require a prior knowledge of some reservoir properties (e.g. relative permeability, absolute permeability, formation porosity and thickness) to calculate the average pressure. To investigate its potential and reliability in analyzing CO 2 injection data, the proposed technique is applied to four synthetic cases representing different well operating conditions. Results indicate that the approximate approach consistently overestimates the actual (simulated) storage capacity as compared to the rigorous approach. The agreement - between the inferred and the simulated reservoir pore volume, and between the analytical and numerical estimates of storage capacity - validates the potential application of the technique to CO 2 storage in closed saline aquifers. The technique is further substantiated through application to a field data set utilized from a commercial-scale geological storage (CGS) project. Finally, field data interpretation shows that the proposed technique can be utilized to identify the degree of hydraulic continuity and reservoir compartmentalization within a target formation by interpreting the corresponding pressure and rate responses.

02 PETROLEUM↗

DOE OSTI2020

Creation of a Mixed-Mode Fracture Network at Mesoscale Through Hydraulic Fracturing and Shear Stimulation

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

58 GEOSCIENCES↗

DOE OSTI2026

Numerical Analysis of Offshore Wind-Farm-Induced Drag Effects on Coastal Upwelling Dynamics

Wind farms extract momentum from the atmospheric flow, generating wind-speed deficits both within the plant, and extending downstream. When located offshore, these deficits modulate air-sea coupling, potentially impacting coastal upwelling in sensitive regions. We investigate impacts of wind farms on coastal upwelling using kilometre-scale, three-way-coupled simulations with the coupled ocean-atmosphere-wave-sediment transport system for the US West Coast. Wind-farm effects are represented by a generalised turbine drag formulation, an idealised, height-dependent body force whose magnitude is systematically varied. This approach isolates the leading-order fluid-dynamical response in a realistic coastal configuration. The atmospheric adjustment exhibits an approximately linear relation between drag force and wind-speed deficit, with wakes that expand downstream and increase in magnitude as drag increases. An empirical orthogonal function analysis of sea-surface-temperature anomalies reveals the emergence of a canonical dipole pattern under strong drag forcing. Subsurface diagnostics show consistent shoaling of the mixed layer and suppressed upward velocities in areas near wind-farm region, accompanied by compensating enhancements of shoaling closer to the coast. These results identify turbine drag as a control parameter in assessing interactions between wind-farm wake and coastal upwelling and provide scaling relationships for understanding offshore wind-farm effects on the coastal circulation dynamics.

16 TIDAL AND WAVE POWER↗

DOE OSTI2018

Simulating Surface and Subsurface Water Balance Changes Due to Burn Severity

Forest fires have a significant impact on hydrology, such as reduced infiltration rates leading to increased flooding. However, post-fire water balance changes and the competing hydrologic response of increased runoff and evapotranspiration as a function of burn severity are not well understood. Comparing pre- and post-fire water balance changes is challenging because measurements of fire-disturbed landscapes with the previously undisturbed character are impractical due to non-repetitive observational conditions. We used a physically based modeling experiment to incorporate burn severity data from the Las Conchas fire to approximate model parameterization to evaluate continuous water balance progression for pre- and post-fire simulations using the same forcing conditions. Fire disturbance decreased evapotranspiration and increased overland flow response to precipitation events. The reduction of evapotranspiration often dominated the new water balance compared with the increase in overland flow, resulting in higher soil moisture. However, this modeling experiment also identified a tipping point where increased overland flow from high burn severity sites eclipses the effect of reduced evapotranspiration on the water balance, causing comparatively drier post-fire soils. In particular, high burn severity sites approach a threshold that results in larger changes to overland flow than changes in evapotranspiration, potentially moving the site to an overland flow dominated regime. The shifts in water balance components have implications for how site conditions will change under a range of burn severity scenarios.

54 ENVIRONMENTAL SCIENCES↗