Search NASA⌕ Search

SEARCH · Search NASA

Results for “Geochemical Modelling”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A Review of Lattice-Boltzmann Models Coupled with Geochemical Modeling Applied for Simulation of Advanced Waterflooding and Enhanced Oil Recovery Processes

To maintain economic profit and improve the oil production efficiency after the primary and secondary production phase, advanced waterflooding techniques such as low salinity waterflooding in carbonate reservoirs have been investigated in numerical simulations, laboratory experiments, and field pilot tests. Multiple underlying mechanisms have been proposed based on these studies, and they are still under debate. Various numerical modeling approaches are introduced, but there exists a lack of a pore-scale comprehensive modeling scheme to fully understand the processes. Lattice-Boltzmann method (LBM) is a type of numerical fluid flow modeling technique that shows capabilities and flexibilities in modeling pore-scale fluid flow to integrate physical–chemical processes within complex structures. The intrinsic feature of LBM makes it a promising framework for simulating advanced waterflooding due to its flexibility, accuracy, and parallel efficiency. LBM works either by itself for solving reactive transport problems or by coupling with a third-party reaction solver. This review mainly introduces the LBM fluid flow and reactive transport capabilities and the concept and modeling approaches to simulate advanced waterflooding techniques. Meanwhile, an evaluation of the coupled LBM models for enhanced oil recovery (EOR) simulations is discussed with future research challenges and directions concluded.

02 PETROLEUM↗
DOE OSTI2022

Geochemical Modeling of Celestite (SrSO 4 ) Precipitation and Reactive Transport in Shales

Celestite (SrSO 4 ) precipitation is a prevalent example of secondary sulfate mineral scaling issues in hydraulic fracturing systems, particularly in basins where large concentrations of naturally occurring strontium are present. Herein, we present a validated and flexible geochemical model capable of predicting celestite formation under such unconventional environments. Simulations were built using CrunchFlow and guided by experimental data derived from batch reactors. These data allowed the constraint of key kinetic and thermodynamic parameters for celestite precipitation under relevant synthetic hydraulic fracturing fluid conditions. Effects of ionic strength, saturation index, and the presence of additives were considered in the combined experimental and modeling construction. This geochemical model was then expanded into a more complex system where interactions between hydraulic fracturing fluids and shale rocks were allowed to occur subject to diffusive transport. We find that the carbonate content of a given shale and the presence of persulfate breaker in the system strongly impact the location and extent of celestite formation. The results of this study provide a novel multicomponent reactive transport model that may be used to guide future experimental design in the pursuit of celestite and other sulfate mineral scale mitigation under extreme conditions typical of hydraulic fracturing in shale formations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Unraveling Hydrogen Induced Geochemical Reaction Mechanisms through Coupled Geochemical Modeling and Machine Learning

Underground hydrogen storage (UHS) provides a promising large-scale, long-term energy storage solution. A reasonable recovery of stored hydrogen is critical for a successful storage scheme. However, in subsurface reservoirs hydrogen is subject to active geochemical reactions that might result in hydrogen loss. In this study, we implemented a geochemical modeling approach coupled with an unsupervised machine learning technique called non-negative matrix factorization (NMF) to unravel the complex brine-rock-H 2 geochemical processes responsible for hydrogen losses, with particular focus on sulfate reduction reactions. NMF is applied to modeled mineral evolution and fluid component profiles to retrieve profiles that can be interpreted to more easily assess competing processes. NMF decouples simulated competing equilibrium reactions. This facilitates separation of overlapping reaction profiles from redox processes, dissolution fronts, and secondary precipitation while considering the effects of simulation parameters such as salinity, temperature, and total H 2 pressure. NMF successfully discriminates these competing effects in nonlinear ways, allowing robust interpretation. In addition, NMF reveals subtle coupled mineral associations and reaction fronts that are invisible to conventional model analysis. This integrated approach strengthens the conceptual understanding of complex nonlinear hydrogen-brine-rock interactions and advances geochemical research on UHS systems to resolve complexities in modeled geochemical systems without the need for direct experiments or prior knowledge. Furthermore, this study highlights the efficacy of combining geochemical modeling with machine learning techniques to enhance the interpretability of the intricate geochemical simulation output through deciphering the overlapping reaction path that cannot be achieved only using conventional analysis of geochemical models alone.

08 HYDROGEN↗
DOE OSTI2023

Geochemical Modeling of the Confining Unit at a CO 2 Storage Site in the Williston Basin, North Dakota

Conference presentation for American Association of Petroleum Geologists (AAPG) Carbon Capture, Utilization, and Storage (CCUS) Conference, Houston, TX, April 25–27, 2023. This geochemical modeling study investigates the evolution of geochemical reactions among CO 2 , formation water, and rocks of the two confining units above and below the CO 2 storage unit at one of the permitted CO 2 storage sites in the Williston Basin. The model is built on PHREEQC-coded chemical equilibrium thermodynamics.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

Geochemical Modeling Using OLI Systems Mixed Solvent Electrolyte Model in Support of the Waste Isolation Pilot Plant

Aqueous electrolyte thermodynamic models were developed for synthetic Salado (GWB) and Castile (ERDA-6) brines to inform the experimental design and to facilitate laboratory data interpretation in geochemical studies supporting disposal in the Waste Isolation Pilot Plant (WIPP). The thermodynamic models, which include the relevant vapor-liquid-solid equilibria for the oxidation-reduction conditions of interest, enable computation of system parameters such as activity-based pH, E h [vs. standard hydrogen electrode (SHE)], ionic strength, total dissolved solids, and solution density; the chemical speciation of all phases; saturation tendencies for thermodynamically feasible solids; and changes in brine chemistry because of the addition of H 2 (g), argon, MgO, Fe(0), and other waste form components. The impacts of oxidation-reduction potential (ORP) and pH on aqueous speciation and dominant solid phases were also assessed with the aid of real-solution Pourbaix diagrams for key elements of interest. In addition to the OLI modeling results, limited measurement data for pH and ORP standards are provided, analyzed, and discussed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Integrated Experimental and Modeling Study of Geochemical Reactions of Simple Fracturing Fluids with Caney Shale

Interactions between rock minerals and hydraulic fluids directly impact the geochemical and geomechanical properties of shale formations. However, the mechanisms of geochemical reactions in shale unconventional reservoirs remain poorly understood. Here, to investigate the geochemical reactions between shale and hydraulic fracturing fluids, a series of batch reactor experiments were undertaken. Three rock samples with different mineralogical compositions and three fluid samples of different compositions (deionized water, deionized water + 2% potassium chloride (KCl), and deionized water + 0.5% choline chloride (C 5 H 14 ClNO) were used. Experiments were undertaken at reservoir temperature and atmospheric pressure. Elemental compositions of effluents after 1, 3, 7, 14, 28 days were analyzed using Inductively Couple Plasma Mass Spectrometry (ICP-MS). Medical Computed Tomography (CT) scan and X-ray Fluorescence (XRF) were conducted on the entire core run to help upscale results obtained from rock-fluid interaction experiments. Geochemical modeling using a reactive simulator, TOUGHREACT, was undertaken to corroborate experimental results. Results show that lower pH triggered high dissolution rates in the rock samples, especially the carbonate components. As pH increased, the rate of dissolution declined significantly, though for most cases dissolution still continued. Observed dissolved silica concentrations were much higher than the quartz solubility, suggesting that much of the silica originates from more soluble silica polymorphs and possibly desorption from clay mineral exchange sites. Concentration of most elemental species in solution increased but aluminium and magnesium concentrations declined rapidly following initial entry into solution. Geochemical modeling corroborated the conclusions regarding mineral dissolution and precipitation observed from experiments, notably; the dissolution of calcite and oxidation of pyrite in reacted shale samples, the likely presence of silica polymorphs such as opal, chalcedony or amorphous silica in these samples, and the reduction of Al and Mg concentrations in solution by precipitation of secondary aluminosilicate phases. The de-flocculation of clay minerals during reaction implies fines migration after hydraulic fracturing. This is detrimental to reservoir productivity as clay fines are displaces and lodged within the micro and nano-fractures created during fracturing. The immediate consumption of aluminium and magnesium also has implications on blockage of hydrocarbon pathways due to precipitation of new minerals in these locations.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2022

Stable carbon isotope values of syndepositional carbonate spherules and micrite record spatial and temporal changes in photosynthesis intensity

Marine and lacustrine carbonate minerals preserve carbon cycle information, and their stable carbon isotope values (δ 13 C) are frequently used to infer and reconstruct paleoenvironmental changes. However, multiple processes can influence the δ 13 C values of bulk carbonates, confounding the interpretation of these values in terms of conditions at the time of mineral precipitation. Co-existing carbonate forms may represent different environmental conditions, yet few studies have analyzed δ 13 C values of syndepositional carbonate grains of varying morphologies to investigate their origins. Here, we combine stable isotope analyses, metagenomics, and geochemical modeling to interpret δ 13 C values of syndepositional carbonate spherules (>500 μm) and fine-grained micrite (<63 μm) from a ~1600-year-long sediment record of a hypersaline lake located on the coral atoll of Kiritimati, Republic of Kiribati (1.9°N, 157.4°W). Petrographic, mineralogic, and stable isotope results suggest that both carbonate fractions precipitate in situ with minor diagenetic alterations. The δ 13 C values of spherules are high compared to the syndepositional micrite and cannot be explained by mineral differences or external perturbations, suggesting a role for local biological processes. We use geochemical modeling to test the hypothesis that the spherules form in the surface microbial mat during peak diurnal photosynthesis when the δ 13 C value of dissolved inorganic carbon is elevated. In contrast, we hypothesize that the micrite may precipitate more continuously in the water as well as in sub-surface, heterotrophic layers of the microbial mat. Both metagenome and geochemical model results support a critical role for photosynthesis in influencing carbonate δ 13 C values. The down-core spherule–micrite offset in δ 13 C values also aligns with total organic carbon values, suggesting that the difference in the δ13C values of spherules and micrite may be a more robust, inorganic indicator of variability in productivity and local biological processes through time than the δ 13 C values of individual carbonate forms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

ATS (Advanced Terrestrial Simulator) integrated hydrology and reactive transport model output in Copper Creek, Colorado

This dataset is generated using the ATS (Advanced Terrestrial Simulator) model at Copper Creek, Colorado, the largest catchment in the East River watershed. ATS is an integrated hydrology and reactive transport model to simulate the Concentration-Discharge (C-Q) relationship, and is used to quantify the geochemical export and understand how watersheds respond to climate disturbances. Key hydrogeochemical variables, including discharge, precipitation, evapotranspiration, and geochemical concentrations, were simulated and output at a daily frequency at the Copper Creek outlet. A group of hydrology and geochemical model variables over the spatial domain are output at a bi-monthly frequency. The model used 1-km resolution daily frequency Daymet meteorological forcing, and mineral composition from previous studies. The simulations are performed on NERSC. The datasets include the input files for running the ATS (Advanced Terrestrial Simulator) simulation, including the mesh file, input script and geochemical model files. The output files from the ATS (Advanced Terrestrial Simulator) are the daily time-series discharge and geochemical concentration simulated at the outlet of Copper Creek, and the bi-monthly spatial output over the simulation domain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

UV-Vis spectrophotometric determination of rare earth elements (REE) speciation at near-neutral to alkaline pH. Part I: m-cresol purple properties from 25-75 °C and Er hydrolysis

The speciation and mobility of rare earth elements (REE) strongly depends on pH which controls the formation of charged aqueous hydroxyl species. The latter potentially play an important role in controlling heavy REE adsorption on clay minerals in near-neutral to alkaline waters such as in regolith-hosted REE mineral deposits. However, accurate REE hydrolysis constants are needed for developing geochemical models that can predict the role of these charged species in natural systems. Here, we develop a robust experimental UV-Vis spectrophotometric method using m-cresol purple to determine in situ pH from 25 to 75 °C. This method is used to derive the average ligand number and hydrolysis constants of erbium (Er) at 25 °C in aqueous solutions with low ionic strength (≤ 0.001 mol/L) at pH from ~7 to 9.5 and in the presence of Er concentrations from 0 to 0.057 mM. The average ligand number ranges between 1 and 3 indicating that Er(OH) 2+ , Er(OH) 2 + and Er(OH) 3 0 control speciation in the experiments. The logarithm of the Er hydrolysis constants (log*β n ° , n= 1 to 3) derived at infinite dilution for the reaction Er 3+ + nH 2 O = Er(OH) n 3-n + nH + are: *β 1 ° = –7.22 ± 0.10, *β 2 ° = –14.52 ± 0.08, *β 3 ° = –23.24 ± 0.04. Implementation of these experimental data into a geochemical model indicates that the Er(OH) 2 + and Er(OH) 3 0 species are both stable in a much wider pH range than previously predicted. Consequently, the positively charged REE hydroxyl complexes can potentially control the fractionation of light vs. heavy REE via adsorption as observed in the formation of certain regolith-hosted REE deposits.

58 GEOSCIENCES↗
DOE OSTI2024

Hydrothermal solubility of Dy hydroxide as a function of pH and stability of Dy hydroxyl aqueous complexes from 25 to 250 °C

The rare earth elements (REE) have important applications in green energy technologies. The formation of mineral deposits in geologic systems commonly involves hydrothermal fluids which can mobilize the REE. However, the REE speciation is not well known as a function of pH. The thermodynamic properties of REE hydroxyl complexes used in geochemical models are based on the Helgeson-Kirkham-Flowers (HKF) equation of state parameters which were derived by extrapolation of low temperature experimental and estimated data. In this study, Dy hydroxide solubility experiments are combined with available literature data to improve these models from 25 to 250 °C and optimize the thermodynamic properties of Dy 3+ and Dy hydroxyl complexes using GEMSFITS. Batch-type solubility experiments were conducted from 150 to 250 °C and at saturated water vapor pressure in perchloric acid solutions with initial pH values of 2 to 5 in 0.5 pH unit increments. The measured solubility of Dy hydroxide is retrograde with temperature and decreases with pH. The logarithm of total dissolved Dy molality ranges from –2.3 to –5.3 at 150 °C (pH 4.7–5.5), from –2.4 to –5.6 at 200 °C (pH 3.9–5.1), and from –3.7 to –6.9 at 250 °C (pH of 3.4 and 5.0). The optimized standard partial molal Gibbs energies of formation (Δ f G° T ) derived for Dy 3+ and DyOH 2+ display a close to linear relationship with temperature, fitting with previous optimizations based on DyPO 4 solubility data in the literature. A comparison of the optimized ΔfG°T values for aqueous Dy species with predictions from available HKF parameters indicates significant differences ranging from +11 to –26 kJ/mol between 25 and 250 °C. The experimental fits are used to derive the Dy hydroxide solubility products (K s0 ) and formation constants for the hydrolysis of Dy (β n with n = 1 to 3; Dy 3+ + nOH – = DyOH n 3-n ) as a function of temperature. The optimization method presented yields accurate thermodynamic properties for the Dy 3+ aqua ions and the DyOH 2+ species at the acidic to mildly acidic pH studied whereas more experimental work is needed at near-neutral and alkaline conditions to better constrain the other hydroxyl complexes. Furthermore, the optimized thermodynamic data have a significant impact on geochemical modeling of the mobility and solubility of REE minerals in acidic hydrothermal fluids.

58 GEOSCIENCES↗
DOE OSTI2023

Coupled Lattice Boltzmann Modeling Framework for Pore-Scale Fluid Flow and Reactive Transport

In this paper, we propose a modeling framework for pore-scale fluid flow and reactive transport based on a coupled lattice Boltzmann model (LBM). We develop a modeling interface to integrate the LBM modeling code parallel lattice Boltzmann solver and the PHREEQC reaction solver using multiple flow and reaction cell mapping schemes. The major advantage of the proposed workflow is the high modeling flexibility obtained by coupling the geochemical model with the LBM fluid flow model. Consequently, the model is capable of executing one or more complex reactions within desired cells while preserving the high data communication efficiency between the two codes. Meanwhile, the developed mapping mechanism enables the flow, diffusion, and reactions in complex pore-scale geometries. We validate the coupled code in a series of benchmark numerical experiments, including 2D single-phase Poiseuille flow and diffusion, 2D reactive transport with calcite dissolution, as well as surface complexation reactions. The simulation results show good agreement with analytical solutions, experimental data, and multiple other simulation codes. In addition, we design an AI-based optimization workflow and implement it on the surface complexation model to enable increased capacity of the coupled modeling framework. Compared to the manual tuning results proposed in the literature, our workflow demonstrates fast and reliable model optimization results without incorporating pre-existing domain knowledge.

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

AquaMEND: Reconciling multiple impacts of salinization on soil carbon biogeochemistry

Soil salinization, exacerbated by climate change, poses a global threat to coastal ecosystem function and soil quality. Salinity influences carbon cycling through direct effects on microbial activity and indirect alterations to soil physicochemical properties including cation exchange, pH, and soil organic carbon availability. Current models inadequately represent these complexities, relying on linear reduction functions that overlook specific physicochemical changes induced by salinity. To address this gap, we propose an integrated model framework, AquaMEND, that combines microbial-explicit carbon decomposition and geochemical models. This model allows cation exchange and surface complexation processes to capture solute chemistry and nutrient availability in soils upon saltwater intrusion. Using response functions that capture salinity impacts on both salt-sensitive and salt-resistant microbial processes, AquaMEND simulates how the abiotic and biotic mechanisms work individually and collectively to regulate organic and inorganic pools and fluxes. Here, the parallel structure of aqueous and non-aqueous phases, together with microbial functions, result in a versatile model for solving dynamic coupling of organics, minerals and microbes under various environmental settings.

54 ENVIRONMENTAL SCIENCES↗