Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Aquifer”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Machine-learning-assisted high-temperature reservoir thermal energy storage optimization

High-temperature reservoir thermal energy storage (HT-RTES) has the potential to become an indispensable component in achieving the goal of the net-zero carbon economy, given its capability to balance the intermittent nature of renewable energy generation. In this study, a machine-learning-assisted computational framework is presented to co-optimize the performance metrics of HT-RTES by combining physics-based simulation with stochastic hydrogeologic formation and thermal energy storage operation parameters, artificial neural network regression of the simulation data, and genetic algorithm-enabled multi-objective optimization. A doublet well configuration with a layered (aquitard-aquifer-aquitard) generic reservoir is simulated for cases of continuous operation and seasonal-cycle operation scenarios. Further, neural network-based surrogate models are developed for the two scenarios and applied to generate the Pareto fronts of the HT-RTES performance for four potential HT-RTES sites. The developed Pareto optimal solutions indicate the performance of HT-RTES is operation-scenario (i.e., fluid cycle) and reservoir-site dependent, and the performance metrics have competing effects for a given site and a given fluid cycle. The developed neural network models can be applied to identify suitable sites for HT-RTES, and the proposed framework sheds light on the design of resilient HT-RTES systems.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Characterizing the diverse hydrogeology underlying rivers and estuaries using new floating transient electromagnetic methodology

We report the hydrogeology below large surface water features such as rivers and estuaries is universally under-informed at the long reach to basin scales (tens of km+). This challenge inhibits the accurate modeling of fresh/saline groundwater interfaces and groundwater/surface water exchange patterns at management-relevant spatial extents. Here we introduce a towed, floating transient electromagnetic (TEM) system (i.e. FloaTEM) for rapid (up to 15 km/h) high resolution electrical mapping of the subsurface below large water bodies to depths often a factor of 10 greater than other towed instruments. The novel FloaTEM system is demonstrated at a range of diverse 4th through 6th-order riverine settings across the United States including 1) the Farmington River, near Hartford, Connecticut; 2) the Upper Delaware River near Barryville, New York; 3) the Tallahatchie River near Shellmound, Mississippi; and, 4) the Eel River estuary, on Cape Cod, near Falmouth, Massachusetts. Airborne frequency-domain electromagnetic and land-based towed TEM data are also compared at the Tallahatchie River site, and streambed geologic scenarios are explored with forward modeling. A range of geologic structures and pore water salinity interfaces were identified. Process-based interpretation of the case study data indicated FloaTEM can resolve varied sediment-water interface materials, such as the accumulation of fines at the bottom of a reservoir and permeable sand/gravel riverbed sediments that focus groundwater discharge. Bedrock layers were mapped at several sites, and aquifer confining units were defined at comparable resolution to airborne methods. Terrestrial fresh groundwater discharge with flowpaths extending hundreds of meters from shore was also imaged below the Eel River estuary, improving on previous hydrogeological characterizations of that nutrient-rich coastal exchange zone. In summary, the novel FloaTEM system fills a critical gap in our ability to characterize the hydrogeology below surface water features and will support more accurate prediction of groundwater/surface water exchange dynamics and fresh-saline groundwater interfaces.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Controls on the 36 Cl/Cl input ratio of paleo-groundwater in arid environments: New evidence from 81 Kr/Kr data

In this study, measurements of the long-lived 81 Kr and 36 Cl radioisotopes in groundwater from the Negev Desert (Israel) were used to assess the 36 Cl/Cl input ratios and Cl – contents for paleorecharge into the Nubian Sandstone Aquifer (NSA). The reconstructed Cl – content of the recharge flux was on the order of 300–400 mg/L. An initial 36 Cl/Cl ratio of 50 × 10 –15 was assessed for the groundwater replenishment in the Negev Desert since the late Pleistocene, in agreement with the 36 Cl/Cl ratios in recent local rainwater. This is despite possible changes in the climatic conditions and the 36 Cl production rates in the atmosphere over this timeframe. This similarity in values is explained by the major role played by the erosion and weathering of near-surface materials in the desert environment that dominate the hydrochemistry of rains, floods, and the consequent groundwater recharge. Spatial variation in the reconstructed initial 36 Cl/Cl ratio is accounted for by the differences in the mineral aerosol sources for specific recharge areas of the NSA. Accordingly, regional variations in the 36 Cl/Cl input in groundwater reservoirs surrounding the Mediterranean Sea indicate various processes that govern the 36 Cl/Cl system. Finally, the results of this study highlight the great advantage of integrating 81 Kr age information in evaluating the initial 36 Cl/Cl and Cl – input, which is essential for the calibration of 36 Cl radioisotope as an available long-term dating tool for a given basin.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Large seasonal fluctuations of groundwater radioiodine speciation and concentrations in a riparian wetland in South Carolina

Recent studies evaluating multiple years of groundwater radioiodine ( 129 I) concentration in a riparian wetland located in South Carolina, USA identified strong seasonal concentration fluctuations, such that summer concentrations were much greater than winter concentrations. These fluctuations were observed only in the wetlands but not in the upland portion of the plume and only with 129 I, and not with other contaminants of anthropogenic origin: nitrate/nitrite, strontium-90, technecium-99, tritium, or uranium. This unexplained observation was hypothesized to be the result of strongly coupled processes involving hydrology, water temperature, microbiology, and chemistry. To test this hypothesis, an extensive historical groundwater database was evaluated, and additional measurements of total iodine and iodine speciation were made from recently collected samples. During the summer, the water table decreased by as much as 0.7 m, surface water temperature increased by as much as 15 °C, and total iodine concentrations were consistently greater (up to 680%) than the following winter months. Most of the additional iodine observed in the summer could be attributed to proportional gains in organo-iodine, and not iodide or iodate. Furthermore, 129 I concentrations were observed to be two-orders-of-magnitude greater at the bottom of the upland aquifer than at the top. A coupled hydrological and biogeochemical conceptual model is proposed to tie these observations together. First, as the surface water temperature increased during the summer, microbial activity was enhanced, which in turn stimulated the formation of mobile organo-I. Hydrological processes were also likely involved in the observed iodine seasonal changes: (1) as the water table decreased in summer, the remaining upland water entering the wetland was comprised of a greater proportion of water containing elevated iodine concentrations from the low depths, and (2) water flow paths in summer changed such that the wells intercepted more of the contaminant plume and less of the diluting rainwater (due to evapotranspiration) and streamwater (as the lower levels promote a predominantly recharging system). Furthermore, these results underscore the importance of coupled processes influencing contaminant concentrations, and the need to assess seasonal contaminant variations to optimize long-term monitoring programs of wetlands.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A country-scale green energy-water-hydrogen nexus: Jordan as a case study

Many developing countries suffer from a shortage of clean energy-water production, causing significant dependency on imported fossil fuel to meet the local energy and freshwater demands. Furthermore, when planning a 100 % renewable energy system (RES) to match the energy demand on an hourly basis, the RES is usually oversized. Hence, unavoidable large amounts of excess energy would be generated during hours of high production, low demand, and fully charged energy storage systems (ESSs). Hence, using Jordan as a case study, this work proposes a novel integrated system of wind, solar photovoltaic (PV), and lithium-ion ESS to match 100 % of the country's energy demand while using the excess generated power to drive reverse osmosis water desalination plants to match the demand of freshwater as well. Furthermore, the remaining excess energy is used for producing green hydrogen. A techno-economic optimizer was developed to show that this concept is indeed feasible. Here, the model was used to scan Jordan for the fitting sites of the highest demand-supply matching and the lowest levelized cost of electricity (LCOE), simultaneously. Three different integration scenarios are presented, namely PV-ESS, wind-ESS, and hybrid PV-wind-ESS RESs. The hybrid system showed the best performance, especially compared to the wind-ESS RES, in terms of total installed capacity (33.37 GW), LCOE (0.0492 USD/kWh), specific water cost (0.3629 USD/m 3 ), and energy and water demand-supply fractions (99.47 and 96.16 %, respectively). Storing desalinated water in natural aquifers was found to be a favorable option to enhance the performance of the three systems, where stored water by the end of the year was sufficient to cover the freshwater requirements for several following years. Finally, the remaining excess energy is used to produce green hydrogen with an annual production of up to 1.37 million tons and a specific cost down to 1.08 USD/kg.

14 SOLAR ENERGY↗
DOE OSTI2020

Carbon-Negative Biofuel Production

Achievement of the 1.5 °C limit for global temperature increase relies on the large-scale deployment of carbon dioxide removal (CDR) technologies. In this article, we explore two CDR technologies: soil carbon sequestration (SCS), and carbon capture and storage (CCS) integrated with cellulosic biofuel production. These CDR technologies are applied as part of decentralized biorefinery systems processing corn stover and unfertilized switchgrass grown in riparian zones in the Midwestern United States. Cover crops grown on corn-producing lands are chosen from the SCS approach, and biogenic CO 2 in biorefineries is captured, transported by pipeline, and injected into saline aquifers. The decentralized biorefinery system using SCS, CCS, or both can produce carbon-negative cellulosic biofuels (≤-22.2 gCO 2 MJ –1 ). Meanwhile, biofuel selling prices increase by 15–45% due to CDR costs. Economic incentives (e.g., cover crop incentives and/or a CO 2 tax credit) can mitigate price increases caused by CDR technologies. Lastly, a combination of different CDR technologies in decentralized biorefinery systems is the most efficient method for greenhouse gas (GHG) mitigation, and its total GHG mitigation potential in the Midwest is 0.16 GtCO 2 year –1 .

09 BIOMASS FUELS↗
DOE OSTI2021

Respiratory Vinyl Chloride Reductive Dechlorination to Ethene in TceA-Expressing Dehalococcoides mccartyi

Bioremediation of chlorinated ethenes in anoxic aquifers hinges on organohalide-respiring Dehalococcoidia expressing vinyl chloride (VC) reductive dehalogenase (RDase). The tceA gene encoding the trichloroethene-dechlorinating RDase TceA is frequently detected in contaminated groundwater but not recognized as a biomarker for VC detoxification. We demonstrate that tceA-carrying Dehalococcoides mccartyi (Dhc) strains FL2 and 195 grow with VC as an electron acceptor when sufficient vitamin B 12 (B 12 ) is provided. Strain FL2 cultures that received 50 μg L –1 B 12 completely dechlorinated VC to ethene at rates of 14.80 ± 1.30 μM day –1 and attained 1.64 ± 0.11 × 10 8 cells per μmol of VC consumed. Strain 195 attained similar growth yields of 1.80 ± 1.00 × 10 8 cells per μmol of VC consumed, and both strains could be consecutively transferred with VC as the electron acceptor. Proteomic analysis demonstrated TceA expression in VC-grown strain FL2 cultures. Resequencing of the strain FL2 and strain 195 tceA genes identified non-synonymous substitutions, although their consequences for TceA function are currently unknown. The finding that Dhc strains expressing TceA respire VC can explain ethene formation at chlorinated solvent sites, where quantitative polymerase chain reaction analysis indicates that tceA dominates the RDase gene pool.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Contribution of Microorganisms with the Clade II Nitrous Oxide Reductase to Suppression of Surface Emissions of Nitrous Oxide

The sources and sinks of nitrous oxide, as control emissions to the atmosphere, are generally poorly constrained for most environmental systems. Initial depth-resolved analysis of nitrous oxide flux from observation wells and the proximal surface within a nitrate contaminated aquifer system revealed high subsurface production but little escape from the surface. Further, to better understand the environmental controls of production and emission at this site, we used a combination of isotopic, geochemical, and molecular analyses to show that chemodenitrification and bacterial denitrification are major sources of nitrous oxide in this subsurface, where low DO, low pH, and high nitrate are correlated with significant nitrous oxide production. Depth-resolved metagenomes showed that consumption of nitrous oxide near the surface was correlated with an enrichment of Clade II nitrous oxide reducers, consistent with a growing appreciation of their importance in controlling release of nitrous oxide to the atmosphere. Our work also provides evidence for the reduction of nitrous oxide at a pH of 4, well below the generally accepted limit of pH 5.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Molecular Structure and Thermodynamics of CO 2 and Water Adsorption on Mica

The adsorption of CO 2 and water on clay surfaces plays a key role in applications, such as gas storage in saline aquifers and depleted hydrocarbon reservoirs, but is not yet fully understood. Here, we study the adsorption of CO 2 and water vapor using Grand Canonical Monte Carlo and molecular dynamics simulations. At a bulk pressure of 100 bar, pure CO 2 adsorbs strongly on mica and forms extensive layers next to it. CO 2 adsorption is lowered substantially if introducing water vapor above mica and is largely eliminated when the relative humidity (RH) approaches about 60%. When pure water vapor is introduced above a mica surface, a subnanometer thick liquid water film develops on it to form apparent liquid–solid and liquid–vapor interfaces simultaneously. Using the identification of truly interfacial molecules (ITIM) analysis, we delineate how individual water layers develop in this film as RH increases. We highlight that the water film is spatially heterogeneous and the true liquid–vapor interface emerges only at an RH of 60–80%. Introducing 100 bar of CO 2 into the water vapor above the mica surface modulates water adsorption nonlinearly: at RH = 0.01%, the water adsorption is reduced by ∼30%; as RH increases, the reduction is weakened, and eventually, enhancement of water adsorption by about 7% occurs at RH = 90%. These variations are attributed to the interplay of film thinning by high-pressure CO 2 , competition of mica surface sites by CO 2 molecules, and energetic and entropic stabilization of interfacial water by CO 2 molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

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↗
DOE OSTI2020

Adaptive Multirate Mass Transfer (aMMT) Model: A New Approach to Upscale Regional‐Scale Transport Under Transient Flow Conditions

Abstract The long‐term evaluation of regional‐scale groundwater quality needs efficient upscaling methods for transient flow. Upscaling techniques, such as the Multirate Mass Transfer (MRMT) method with constant upscaling parameters, have been used for transport with steady‐state flow, yet the upscaling parameters (i.e., rate coefficients) may be time dependent. This study proposed and validated an adaptive MRMT (aMMT) method by allowing the mass transfer coefficients in MRMT to change with the flow field. Advective‐dispersive contaminant transport simulated in a 3‐D heterogeneous medium was used as a reference solution. Equivalent transport under homogeneous flow conditions was evaluated by applying the MRMT and aMMT models for upscaling. The relationship between mass transfer coefficients and flow rates was fitted under steady‐state flow driven by various hydraulic gradients. A power law relationship was obtained, which was then used to update the mass transfer coefficients in each stress period under transient flow conditions in the aMMT method. Results indicated that for advection‐dominated transport, both the MRMT and aMMT methods can upscale the anomalous transport dynamics affected by subgrid heterogeneity under transient flow conditions. Whereas for diffusion‐dominated systems, the MRMT model failed to capture the tails of tracer breakthrough curves after the boundary condition changed, but the results from the aMMT model were significantly improved. However, if the overall flow direction changed, both MRMT and aMMT failed to represent the breakthrough curve tail generated by the heterogeneous system. The results point toward a promising path for upscaling transport in complex aquifers with transient flow.

Guo, Zhilin↗
DOE OSTI2021

Utilizing Environmental Tracers to Reduce Groundwater Flow and Transport Model Parameter Uncertainties

Abstract Non‐uniqueness in groundwater model calibration is a primary source of uncertainty in groundwater flow and transport predictions. In this study, we investigate the ability of environmental tracer information to constrain groundwater model parameters. We utilize a pilot point calibration procedure conditioned to subsets of observed data including: liquid pressures, tritium ( 3 H), chlorofluorocarbon‐12 (CFC‐12), and sulfur hexafluoride (SF 6 ) concentrations; and groundwater apparent ages inferred from these environmental tracers, to quantify uncertainties in the heterogeneous permeability fields and infiltration rates of a steady‐state 2‐D synthetic aquifer and a transient 3‐D model of a field site located near Riverton, Wyoming (USA). To identify the relative data worth of each observation data type, the post‐calibration uncertainties of the optimal parameters for a given observation subset are compared to that from the full observation data set. Our results suggest that the calibration‐constrained permeability field uncertainties are largest when liquid pressures are used as the sole calibration data set. We find significant reduction in permeability uncertainty and increased predictive accuracy when the environmental tracer concentrations, rather than apparent groundwater ages, are used as calibration targets in the synthetic model. Calibration of the Riverton field site model using environmental tracer concentrations directly produces infiltration rate estimates with the lowest uncertainties, however; permeability field uncertainties remain similar between the environmental tracer concentration and apparent groundwater age calibration scenarios. This work provides insight on the data worth of environmental tracer information to calibrate groundwater models and highlights potential benefits of directly assimilating environmental tracer concentrations into model parameter estimation procedures.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Direct Observation of the Depth of Active Groundwater Circulation in an Alpine Watershed

Abstract The depth of active groundwater circulation is a fundamental control on stream flows and chemistry in mountain watersheds, yet it remains challenging to characterize and is rarely well constrained. We collected hydraulic conductivity, hydraulic head, temperature, chemical, noble gas, and 3 H/ 3 He groundwater age data from discrete levels in two boreholes 46 and 81 m deep in an alpine watershed, in combination with chemical and age data from shallow groundwater discharge, to discern groundwater flow rates at different depths and directly observe active and inactive groundwater. Vertical head gradients are steep (average of 0.4) and thermal profiles are consistent with typical linear conductive continental geotherms. Groundwater deeper than ∼20 m is distinct from shallow groundwater and creek water in its chemistry, noble gas signature, and age (dominantly >65 years compared to <9 years). Together these results suggest low vertical groundwater flow velocities and a relatively shallow active circulation depth of ∼20 m. This hypothesis is tested with a simple 2‐D numerical fluid flow and heat transport model representing a hillslope transect through the two boreholes. The modeling indicates that the subhorizontally bedded sedimentary rocks underlying the basin are highly anisotropic with low vertical hydraulic conductivity, and at most ∼10% of bedrock recharge (equivalent to <2% of stream baseflow) flows below a depth of 20 m. The study demonstrates the considerable value of discrete‐depth hydrogeologic, chemical, and age data for determining active circulation depth, and illustrates an approach for maximizing the utility of individual boreholes drilled for mountain bedrock aquifer characterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Numerical and Analytical Modeling of Flow Partitioning in Partially Saturated Fracture Networks

Infiltration processes in fractured-porous media remain a crucial, yet not very well understood component of recharge and vulnerability assessment. Under partially-saturated conditions flows in fractures, percolating fracture networks, and fault zones contribute to the fastest spectrum of infiltration velocities via preferential pathways. Specifically, the partitioning dynamics at fracture intersections determine the magnitude of flow fragmentation into vertical and horizontal components and hence the bulk flow velocity and dispersion of fracture networks. In this work, we derive an approximate analytical solution for the partitioning process and validate it using smoothed particle hydrodynamics simulations. The developed transfer function is conceptually based on simulation results and laboratory experiments carried out in previous works. It allows to efficiently simulate flow through fracture networks with simple cubic structure and arbitrary number of fractures and aperture sizes via linear response theory and convolution of a given input signal. We derive a non-dimensional bulk flow velocity ($\widetilde{v}$) and dispersion coefficient ($\widetilde{D}$) to characterize fracture networks in terms of dimensionless horizontal and vertical time scales $\tau_m$ and $\tau_0$. The dispersion coefficient is shown to strongly depend on the horizontal time scale and converges towards a constant value of $0.08$ within reasonable fluid and geometrical parameter ranges, while the non-dimensional velocity exhibits a characteristic $\widetilde{v} \sim \tau_m^{-1/2}$ scaling. Given that hydraulic information is often only available at limited places within (fractured-porous) aquifer system, such as boreholes or springs, our study intends to provide a rudimentary analytical concept to potentially reconstruct internal fracture network geometries from external boundary information, e.g., the dispersive properties of discharge (groundwater level fluctuations).

42 ENGINEERING↗
DOE OSTI2021

Precambrian Crystalline Basement Properties From Pressure History Matching and Implications for Induced Seismicity in the US Midcontinent

Injection-induced seismicity across the US midcontinent has almost exclusively occurred in the crystalline basement that underlies the Arbuckle Group aquifer and its equivalents, the primary wastewater disposal zone in this region. However, the properties of the basement are not well known. Newly compiled data, from Class I wells in Kansas, provide a unique record of pressures in the Arbuckle and an opportunity to constrain the reservoir-scale properties of the basement such as permeability, diffusivity, and specific storage. Constraints on these parameters are critical for modeling fluid flow and pressures across the entire Arbuckle-basement system, and are necessary for accurate evaluation and prediction of injection-induced earthquakes. Here, we present a detailed, three-dimensional geological and pressure history-matched numerical model for the Arbuckle and basement, based on data from >400 wells covering a large region in south-central Kansas, where injection-induced seismicity has been concentrated since 2014. Simulations of dynamic data from 319 wells indicate that Arbuckle pressures have increased by 1.1 MPa in high injection rate areas and an overpressure of <0.1 MPa may be the cause of seismicity in the basement. Pressure-history matching also yields the likely range in porosity (0.3%–7%), permeability (0.1–0.7 mD), and diffusivity (0.004–0.07 m2 /s) for the basement. The resulting estimates suggest reservoir-scale properties of the basement are enhanced by faults and fractures. Importantly, the diffusivities determined in this study are lower than estimates derived from Kansas earthquake triggering fronts, and suggest that such seismicity-based techniques may have limitations, particularly where spacetime patterns between injection and seismicity are complex.

58 GEOSCIENCES↗
DOE OSTI2021

Impacts of Pore Network-Scale Wettability Heterogeneity on Immiscible Fluid Displacement: A Micromodel Study

The mixed-wet nature of reservoir formations imposes a wide range of rock wettability from strong resident-fluid wetting to strong invading-fluid wetting. The characteristics of two-phase flow in porous media composed of mixed-wetting surfaces remain poorly understood. In this study, we investigated the displacement of resident ethylene glycol (EG) by hexane in two mixed-wet micromodels of identical 2.5-D geometry heterogeneity, with uniformly or heterogeneously distributed patches strongly wetting to hexane. These patches are mixed among pores with unaltered EG-wetting surfaces. Along with control tests in the originally EG-wet micromodel, here we show the classic fingering and transitions in flow regimes at $logCa$ (capillary number) from -7.2 to -3.9. Moreover, pore-scale distributions of wettability and their spatial correlation influence displacement efficiency. In the two mixed-wet micromodels, we found (a) an increase of steady-state hexane saturation at the end of experiments by up to 0.12 in the capillary fingering regime and a decrease of at most by 0.06 in the viscous fingering regime, compared to the EG-wet micromodel, and (b) dispersed and fragmented hexane distribution after displacement. Brine drainage during supercritical CO 2 (scCO 2 ) injections in these micromodels occurs with lower wettability contrasts, and under similar viscosity ratios and interfacial tensions resulted in higher displacement efficiency relative to displacement of EG by hexane. While mixed-wettability can enhance displacement efficiency compared to uniform wettability, the dynamics of immiscible fluids in strong mixed-wet reservoirs are expected to be less pronounced in contributing to the efficiency of geological CO 2 sequestration, oil recovery, and remediation of hydrocarbon-contaminated aquifers.

2.5-D micromodel↗
DOE OSTI2022

Geologic Heterogeneity Controls on Trapping and Migration of CO 2

When CO 2 is injected into sedimentary basins, the fate of the buoyant CO 2 will be determined by how it interacts with the underlying geologic heterogeneities. Here we geologically print a cross-bedded pattern using two different bead sizes, and then conduct flow experiments to observe buoyant flow pathways, migration speed, and immobilized volumes. The amount of buoyant phase trapped under the heterogeneities is observed to vary by two orders of magnitude, with the controlling parameter being the size contrast between the coarse and fine grains that make up the structure. Furthermore, the results allow an estimate of heterogeneity trapping, and show how CO 2 trapping in aquifers can be experimentally evaluated before injection.

58 GEOSCIENCES↗
DOE OSTI2022

Transition Zone Morphology Dynamics of Dissolved Oxygen (DO) in a Salinity-Impacted Hyporheic Zone

Freshwater salinization is a common yet underappreciated environmental problem in rivers, yet how this process impacts the transport of dissolved oxygen (DO) in the hyporheic zone remains unknown. Here, using flume experiments and numerical simulations, we have demonstrated that the morphological dynamics of the DO transition zone are controlled by the dimensionless number Da·Ra –2 , a measure of the rate of aerobic respiration to free convection. At low Da·Ra –2 , the oxic front breaks up into fingers which subsequently grow, and leads to accelerating of DO transport and shrinking of the anaerobic zone in the underlying sediment. As Da·Ra –2 increases, the growth rate of the instability decreases, and the DO plume fingers are suppressed and delayed compared with the saltwater-freshwater interface, the mixing area between saltwater and freshwater. These results indicate that the freshwater salinization syndrome can have significant impacts on the functioning of aquifer ecosystems unless regulated and managed effectively.

54 ENVIRONMENTAL SCIENCES↗