SEARCH · Search NASA
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.
This report analyzes experimental data from Test Plans TP 08-02, TP 12-02, and TP 20-01 to add new log K values and Pitzer interaction parameters for Fe, Pb, Mg, Nd and B reactions to the WIPP geochemical thermodynamic database, data0.fm 1.
Within the Haynesville-Bossier Shale complex, the Bossier Shale has not been extensively studied by either industry and academia, despite it being an unconventional gas reservoir and a potential caprock for carbon storage in the underlaying Haynesville Shale. The lack of knowledge of the complex pore structures and fluid-rock interactions hinders the effective extraction of gas and the characterization of fluid reservoirs and sealing capacity. Integrated experimental studies of pore structure and fluid-rock interactions were conducted in seven Bossier Shale core samples collected in eastern Texas. Petrographic, geochemical, and petrophysical properties such as mineral composition, organic richness, thermal maturity, porosity, pore/pore throat diameter distribution, water-accessible pores, liquid water imbibition, and water vapor adsorption were characterized using complementary approaches of thin-section petrography, scanning electron microscopy, X-ray diffraction, total organic matter, pyrolysis, mercury intrusion porosimetry, nuclear magnetic resonance, (Ultra-) small angle X-rays scattering as well as small angle neutron scattering with deuterated liquids and contrast variation. Further, the results show that the thermally mature Bossier Shales are composed of mixed argillaceous mudstone, mixed mudstone, and mixed carbonate mudstone. The shale contains both organic and inorganic pores, with porosities of 3.24-9.37 %, pore-to-throat ratios of 1.65 to 19.4, and water-accessible pores accounting for 28.7-72.6 % of total pores. Approaches of liquid water imbibition and water vapor adsorption, with and without direct contact of water with shale samples, indicate that liquid water first enters the nano-sized pores under high capillary pressures, and water vapor adsorption is mainly controlled by both clay minerals and pores with diameters less than 10 nm. These findings contribute to a better understanding of pore structures and water-shale interactions and their controlling factors in the Bossier Shale.
Reactions at mineral interfaces with aqueous solutions control many geochemical and biogeochemical processes in the Earth’s critical zone. At the molecular level, insights into important properties such as the structure of the electrical double layer (EDL) at specific mineral interfaces continue to improve due to the increasing fidelity of laboratory instrumentation and computational approaches. However, molecular simulation approaches suffer from limited reach into relevant scales of time, length, and system complexity. To span this gap, a novel hybrid approach that couples first principle plane-wave density functional theory (DFT) with classical DFT (cDFT) is demonstrated and applied to calcite (104) interfaces with various electrolytes. In this approach, a region of interest described using DFT interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. Benchmarking against experimental observations and entirely first principle DFT simulations demonstrates that this hybrid model efficiently encompasses the key short-range and collective interactions in the EDL. Simulations of calcite (104)/solution interfaces reveal the key static and dynamic polarization interactions that give rise to structuring of ions and water. Ion hydration interactions have the strongest effect on the depth of the first minimum in the density distribution of counterions at the surface, and the position and width of the first density peak is largely determined by the strength of ion-correlation forces. Finer details of ion distributions are controlled by mutual polarization of the calcite surface and interfacial electrolyte. Finally, this new ability to efficiently and rigorously predict EDL structure at mineral surfaces in contact with complex solutions paves the way to accurately modeling sorption, nucleation, dissolution, and growth in realistic systems.
Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.
The aquatic interfaces exposing terrestrial soils to oxic-anoxic regime shifts represent biogeochemical “hotspots” that are extremely sensitive to climate and environmental change. However, processes and interaction across theses aquatic interfaces are poorly understood and underrepresented in current Earth system models. In this project, we aim to develop predictive understanding of the feedbacks between microbial systems and geochemical environments that determine emergent ecosystem behaviors and resilience in response to disturbances. We use experimental, mechanistic modeling and meta-analysis tools to elucidate interactions among soil, water, geomorphology and microbiology that regulate the molecular transformations and fluxes of carbon, nutrients, and redox-sensitive compounds across aquatic interfaces.
This closeout presentation summarizes a Department of Energy–funded study that evaluated whether large amounts of hydrogen can be safely stored underground in the North Dakota portion of the Williston Basin. The project combined lab testing, computer simulations, and basin‑wide analysis to assess saline formations, depleted oil and gas reservoirs, and salt formations for hydrogen storage capacity, recovery, and risks. Results show that underground hydrogen storage is technically feasible across multiple formation types, with depleted oil and gas reservoirs offering higher recovery and saline formations providing large long‑term storage potential. The study also identifies key challenges—such as wellbore material durability and gas purity management—and recommends pilot projects and further site‑specific studies to support future commercialization.
The Mount Simon formation in the Midwest region of the U.S. is considered as a potential candidate host reservoir for carbon storage. Obtaining knowledge of possible geochemically induced changes to the permeability and porosity of host CO2 storage sandstone will enable us to gain a deeper insight of the long-term reservoir behavior under the CO2 storage conditions. An experimental study of the interaction of CO2/brine/rock on saline formations in a static system under CO2 storage conditions was conducted.
It has been reported that among the various geologic storage options, deep saline aquifers have the largest estimated capacity for CO 2 storage. Obtaining knowledge of possible artificially geochemically induced changes to the permeability and porosity of host CO 2 storage deposits will enable us to gain insight on long-term reservoir behavior under CO 2 storage conditions. An experimental study of the interaction of CO 2 /brine/rock on saline formations was conducted in a static system under CO 2 storage conditions. Chemical interactions in the Cedar Keys-Lawson Formation carbonate during exposure to CO 2 and brine under sequestration conditions were studied. Samples were exposed to the simulated in-situ reaction conditions for one and six months. The samples were exposed to a model brine at 55 °C and CO 2 pressure of 23.8 MPa (3,500 psig). Computed tomography (CT), x-ray diffraction (XRD), scanning electron microscopy (SEM)-energy dispersive x-ray spectroscopy (EDS), brine composition, core porosity, and core permeability analyses were conducted prior to and after the exposure experiments. Preliminary permeability measurements obtained from the core samples showed changes after they were exposed to CO 2 -saturated brine for one and six months. This observation suggests that mineral dissolution and mineral precipitation could occur in the host deposit altering its characteristics for CO 2 storage over time. The 3D images of the pore space clearly illustrate the degree of dissolution that occurred during exposure. It is noted that the dissolution that occurred during the six-month exposure could have enhanced the connectivity between voids. This may contribute the increase of permeability after the CO 2 /brine exposure. In addition, the primary minerals comprising the core are dolomite and gypsum. Both minerals could dissolve in the CO 2 /brine environment resulting in the increase of permeability after the six-month exposure.
ABSTRACT: Studies were conducted on cores and rock cuttings recovered from two wells drilled in Caney shale in the Southern Oklahoma. Mineralogical compositions of these samples were obtained by X-Ray Diffraction (XRD) measurements using two methods of sample preparation: conventional powder (by crushing and grinding of rock fragments) and spray-drying of micronized rock powder method. Microstructure and microanalysis were acquired by Scanning Electron Microscopy/Energy Dispersive Spectroscopy (SEM/EDS). Batch experiments were then conducted using selected rock samples and fracturing fluid whilst deionized water was used as control sample for fluid/rock interaction. These were conducted at 95°C and ambient pressure for various time periods to assess the geochemical reactivity of the Caney shale. The rock samples were subsequently recovered from fluid, dried and analysis repeated using identical XRD and SEM/EDS measurements, following the same conditions. Preliminary experimental results show rock mineralogical compositions are predominantly quartz with relatively high clay content. The clay minerals are mainly illite. Carbonate content is less than 10% by weight for most samples, with other non-clay minerals pyrite, albite, dolomite, and some apatite. In general, sprayed samples showed more intense peaks and less full width at half maxima (FWHM) relative to conventional rock-powder samples. Post-experimental mineralogical changes were observed in samples with clay portions breaking down to amorphous entities. As expected, dissolution of carbonates (calcite and dolomite) was observed and newly detected dissolution of biogenic micro quartz resulted in formation of amorphous silica. 1.0 INTRODUCTION Formation of sedimentary rocks involves processes where rock sediments are emplaced with fluids starting interactions that persist over long periods. During the lithification process, sediments and fluids react and attain equilibrium state. This equilibrium is disturbed during drilling, completion, hydraulic fracturing and other operations that involve the introduction of engineered fluids into the subsurface. To return to equilibrium, dissolution and precipitation of minerals may occur. These reactions have implications on the petrophysical properties of the reservoir. Whiles dissolution and precipitation of minerals can enhance porosity and permeability, the reverse can also be true. For example, when acidic fluids are introduced into formations with high carbonate contents, carbonate minerals dissolve into solution thus increasing the porosity and permeability of the formation. In contrast, when alkaline fluids with pH above 11 are introduced into an illite-rich shale formation, the rock-fluid interaction leads to de-flocculation of illite platelets causing substantial increase in clay fines migration, consequently significant reduction in permeability (Salles et al., 2008).
Here, this study investigates the impacts of geomechanical and geochemical changes on carbon storage in a partially depleted oil reservoir, using results from four different coupled simulation models. Models were used to examine the relative importance of storage mechanisms, and how changing reservoir parameters might affect these mechanisms through time. The study uses data from a Morrowan sandstone reservoir in the Farnsworth Unit (FWU), Ochiltree County, Texas which is currently undergoing CO 2 enhanced oil recovery (EOR). Partially depleted oil reservoirs such as the FWU offer attractive carbon utilization and/or storage targets because of existing infrastructure and economic benefits from incremental oil recovery as well as tax credits. However, prediction of storage capacity or long-term fluid migration in these fields can be difficult because of the wide variation in formation fluids and operational histories that may have undergone. CO 2 injection can cause complex geomechanical and geochemical responses in a reservoir as a result of interplay between dynamic changes in pore pressure, reservoir temperature, fluid composition, and interactions between formation fluids, CO 2 , and reservoir rock. Thus, multiple coupled numerical simulation models must be developed and used to more precisely understand what CO 2 storage mechanisms are most significant, as well as the long-term fate of the stored CO 2 . Our study used results from hydrodynamic, coupled hydro-geomechanical, coupled hydro-geochemical, and coupled hydro-geomechanical-geochemical models to examine how changes in geomechanical and geochemical properties can impact the injectivity or storage capacity of CO 2 . Models simulated historical field operations and then forward-modeled a water-alternate-gas (WAG) operation for 20 years, followed by a 1000-year post-injection monitoring. The work demonstrates that in this specific reservoir, geomechanical impacts appear to be more significant than any geochemical processes, which showed minimal impact on reservoir properties. Furthermore, these simulation scenarios corroborated other studies showing more than 90% of CO 2 injected are being stored and will remain stored in the reservoir. The study was able to draw on a rich set of data from this currently operational project for the work. Because the Morrow is a widespread oil-bearing reservoir in the midcontinent of the United States and is fairly typical of many sandstone reservoirs elsewhere, the study has significance beyond the boundaries of our field area.
ABSTRACT: In this comparative study, we analyzed the changes in elemental concentrations of hydraulic fracturing fluids after interaction with the Marcellus and Caney Shale formations. The focus was on assessing the inherent risks and environmental implications associated with flowback waters, including their impact on soil, ground, and drinking water quality, and human health safety. The chemical compositions of effluents were determined through Inductively Coupled Plasma Optical Emission Spectroscopy (ICP-OES) and Mass spectrometry (ICP-MS). The Marcellus Shale, showed concentrations of Cd, averaging 0.380 ppm far exceeding safe water thresholds. Significant levels of As, Se, B, and Pb were detected in both shales, raising concerns about soil and water contamination. Analyzing the effluents from representative samples of sections of the Marcellus (S2 and S7) and Caney (R1 and R2) indicates different geochemical responses over 4 weeks of experiments. This comparison underscores the chemical changes and environmental considerations linked to hydraulic fracturing across shale formations, suggesting the value of tailored monitoring and regulatory measures for each type. 1. INTRODUCTION The Caney and Marcellus shales, differing in geology and geochemistry, represent distinct unconventional reservoirs. The Caney shale, is more ductile, with higher produced water volumes, (Smith et al., 2022) contrasts with the brittle Marcellus shale known for lower brine production but significant data availability. This study aims to elucidate the possible environmental health and safety issues that may arise from the hydraulic fracturing processes. The interaction between fracturing fluids and clays presents a significant challenge. The primary base of these fluids is water, which, when introduced to clay, can induce swelling and constrict flow pathways. This phenomenon is attributed to water molecules infiltrating the layers of clay, particularly in 2:1-type clays, leading to an increased distance between layers. To counteract this, clay stabilizers are employed (Awejori et al., 2021). Despite their effectiveness, these stabilizers are considered temporary solutions. Upon completion of the fracturing process, a concomitant amount of contaminated waters (flow back), with varied content is collected at the surface. From these, we can infer the geochemical reactions and the environmental challenges associated with these waters. Flowback waters can be reinjected or used for other purposes such as irrigation. This requires adequate screening and treatment for safe use.
Pahute Mesa (Nevada, USA) was the site of 85 underground nuclear tests between 1965 and 1992 whose residual radio-chemical inventory poses a contaminant threat to local groundwater resources. To assess the fate and transport of these radionuclides in groundwater, several geochemical studies, including this one, have been conducted to understand long-term groundwater flow patterns at, and downgradient of, Pahute Mesa. Groundwater samples from the Pahute Mesa area at the Nevada National Security Site were analyzed for noble gases ( 3 He/ 4 He isotopic ratio and He, Ne, Kr, Ar and Xe concentrations) to infer groundwater recharge temperatures and provide additional evidence in support of groundwater flow patterns and the presence of late-pluvial Pleistocene groundwater. These data confirm recent findings (based on groundwater geochemistry and isotope compositions) that most recharge in the area occurred under temperature conditions much cooler than are found today, except near major drainages where ephemeral run of continues to recharge shallow groundwater. These noble gas data indicate that noble gas recharge temperatures and terrigenic He components have been reset where groundwater has flowed upward over barriers formed by caldera margins, or near the regional discharge area (Oasis Valley) where shallow groundwater has interacted with shallow soil gas. This evidence from groundwater noble gas measurements, in combination with other isotopic and geochemical studies, supports the current understanding of the unusual hydrogeology and recharge history of Pahute Mesa.
Rhizosphere arguably embodies the most diverse microbial ecosystem on the planet, yet it is largely a functional ‘black box’ of belowground plant-microbiome interactions. The rhizosphere is the primary site of entry for subsurface injection of fixed carbon (C) into soil with impacts on local to global scale C biogeochemistry and ultimately Earth’s climate. While spatial organization of rhizosphere is central to its function, small scale and steep microbial and geochemical gradients within this dynamic region make it easily disrupted by sampling. The significant challenge presented by sampling blocks elucidation of discreet functions, drivers, and interactions within rhizosphere ecosystems. Here, we describe a non-destructive sampling method linked to metaproteomic analysis in order to measure temporal shifts in the microbial composition and function of rhizosphere. A robust, non-destructive method of sampling microbial hotspots within rhizosphere provides an unperturbed window into the elusive functional interactome of this system over time and space.
Abstract. Soil pore water (SPW) chemistry can vary substantially across multiple scales in Arctic permafrost landscapes. The magnitude of these variations and their relationship to scale are critical considerations for understanding current controls on geochemical cycling and for predicting future changes. These aspects are especially important for Arctic change modeling where accurate representation of sub-grid variability may be necessary to predict watershed-scale behaviors. Our research goal is to characterize intra- and inter-watershed soil water geochemical variations at two contrasting locations in the Seward Peninsula of Alaska, USA. We then attempt to identify the key factors controlling concentrations of important pore water solutes in these systems. The SPW geochemistry of 18 locations spanning two small Arctic catchments was examined for spatial variability and its dominant environmental controls. The primary environmental controls considered were vegetation, soil moisture and/or redox condition, water–soil interactions and hydrologic transport, and mineral solubility. The sampling locations varied in terms of vegetation type and canopy height, presence or absence of near-surface permafrost, soil moisture, and hillslope position. Vegetation was found to have a significant impact on SPW NO3- concentrations, associated with the localized presence of nitrogen-fixing alders and mineralization and nitrification of leaf litter from tall willow shrubs. The elevated NO3- concentrations were, however, frequently equipoised by increased microbial denitrification in regions with sufficient moisture to support it. Vegetation also had an observable impact on soil-moisture-sensitive constituents, but the effect was less significant. The redox conditions in both catchments were generally limited by Fe reduction, seemingly well-buffered by a cache of amorphous Fe hydroxides, with the most reducing conditions found at sampling locations with the highest soil moisture content. Non-redox-sensitive cations were affected by a wide variety of water–soil interactions that affect mineral solubility and transport. Identification of the dominant controls on current SPW hydrogeochemistry allows for qualitative prediction of future geochemical trends in small Arctic catchments that are likely to experience warming and permafrost thaw. As source areas for geochemical fluxes to the broader Arctic hydrologic system, geochemical processes occurring in these environments are particularly important to understand and predict with regards to such environmental changes.
Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance-based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain-scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. Furthermore, these results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems.
Highlights: • Concept to relate microstructural changes, mass transport, and chemical reactions in cement-based materials is presented. • Multi-physics modeling approach utilizing an extended version of the PNP equation and the geochemical code IPHREEQC. • Case studies indicate significant changes in mass transport properties due to alterations in the microstructure. • Microstructural changes depend on exposure and interactions between chemical reactions, pore structure, and mass transport. A generic modeling framework is presented to relate microstructural changes, i.e. changes in the pore volume distribution and phase assemblage, and mass transport in cement-based materials. The modeling framework accounts for mass transport and chemical equilibrium between ions in the pore solution and solid hydrates by means of an extended version of the Poisson-Nernst-Planck equation, in which the chemical equilibrium is solved by the external geochemical code IPHREEQC. Results of numerical studies concerning carbonation and leaching of a cement-based material indicate significant changes in mass transport properties due to chemical reactions leading to alterations in the microstructure. The resulting highly non-linear (both spatial and temporal) microstructural changes are thereby depending on the exposure and indicate that interactions between chemical reactions, pore structure changes, and mass transport have to be accounted for as a whole.
Geogenic arsenic has become a globally-distributed groundwater contaminant, liberated from the weathering of arsenic-bearing sulfide minerals and often transported to aquifer sediments adsorbed to iron oxides. Among the iron oxides, goethite (α-FeOOH) is uniquely important for the fate of arsenic because of its widespread abundance, stability, and high affinity for binding arsenic. Goethite is ubiquitous in soils and sediments and often contains substituted elements, including manganese. Structural manganese may affect the surface reactivity and redox capacity of goethite and alter the mechanisms of recrystallization catalyzed by dissolved Fe(II), potentially affecting arsenic adsorption. Here, this study examined the fate of As(V) during the interactions between dissolved Fe(II) and Mn-substituted goethites at pH 4 and 7 as well as associated changes in arsenic speciation. At pH 7, the addition of dissolved Fe(II) initially increases the adsorption of As(V) onto Mn-bearing and Mn-free goethites. For the Mn-substituted goethites, the adsorbed As(V) slowly releases to solution at longer aging times. Fe(II) addition at pH 4 slightly increases As(V) uptake by Mn-substituted goethites, with differences in total sorption correlating with the Mn content in goethite. The addition of Fe(II) releases substantial dissolved manganese but the amount solubilized is higher at pH 4 compared to 7, suggesting that the presence of adsorbed As(V) may substantially promote the Mn release at pH 4. X-ray absorption fine structure spectroscopy shows that arsenic is stabilized as As(V) in all the samples and adsorbed on goethite via a bidentate binuclear mechanism. Fitting results show that the binding distance and coordination numbers are stable in Mn-free goethite and Mn-substituted goethite samples; the effect of substituted Mn on the surface complex structure is minor. High resolution transmission electron microscopy and X-ray diffraction confirm that no secondary ferrous arsenate minerals precipitate under both pH conditions. This study improves our understanding of the Fe(II)-As(V) interactions on iron oxides, and demonstrates that the substituted cations such as manganese may quantitatively alter the geochemical fate of arsenic during the reaction of dissolved Fe(II) with Fe(III) oxides.
Hydrogen (H 2 ) has the potential to be a transformative technology as an enabler to a low-carbon future and promoter of renewable energy. When H 2 is injected and stored in the subsurface, it has the potential to interact with the caprock (usually shale) which overlies and seals the storage reservoir. Here, this study examines geochemical reactions or changes in surface morphology to Eagle Ford Shale, a proxy for caprock, upon exposure to H 2 at 50°C and 10.3 MPa. Reactions were also performed with N 2 to provide an experimental control. Scanning Electron Microscopy with Energy Dispersive Spectroscopy (SEM-EDS), Atomic Force Microscopy (AFM), and Optical Photothermal-Infrared (O-PTIR) Spectroscopy were applied to quantify changes on the microscale and nanoscale level. Fluid chemistry changes were monitored with Ion Chromatography (IC) and Inductively Coupled Plasma Mass Spectrometry (ICP-MS). Exposure of Eagle Ford Shale to H 2 gas alone did not result in any alterations to the shale chemically or any changes in the surface morphology. Exposure of Eagle Ford Shale to both H 2 and water as well as N 2 and water resulted in changes to the surface morphology because of gypsum dissolution and reprecipitation, thus indicating that H 2 is not necessary to promote changes. Pure pyrite was the most reactive with H 2 possibly resulting in a reduction to pyrrhotite. These initial studies suggest that the extent of reactions activated by hydrogen with caprock are minor under the temperature and pressure conditions that would represent underground hydrogen storage.