skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A 2.5D Reactive Transport Model for Fracture Alteration Simulation

 [1];  [1];  [1];  [2];  [1];  [1];  [1]
  1. Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
  2. Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States; Earth and Planetary Science, University of California, Berkeley, Berkeley, California 94720, United States
Publication Date:
Research Org.:
Energy Frontier Research Centers (EFRC) (United States). Center for Nanoscale Control of Geologic CO2 (NCGC)
Sponsoring Org.:
USDOE Office of Science (SC), Basic Energy Sciences (BES) (SC-22)
OSTI Identifier:
DOE Contract Number:
Resource Type:
Journal Article
Resource Relation:
Journal Name: Environmental Science and Technology; Journal Volume: 50; Journal Issue: 14; Related Information: NCGC partners with Lawrence Berkeley National Laboratory (lead); University of California, Davis; Lawrence Livermore National Laboratory; Massachusetts Institute of Technology; Ohio State University; Oak Ridge National Laboratory; Washington University, St. Louis
Country of Publication:
United States
bio-inspired, mechanical behavior, carbon sequestration

Citation Formats

Deng, Hang, Molins, Sergi, Steefel, Carl, DePaolo, Donald, Voltolini, Marco, Yang, Li, and Ajo-Franklin, Jonathan. A 2.5D Reactive Transport Model for Fracture Alteration Simulation. United States: N. p., 2016. Web. doi:10.1021/acs.est.6b02184.
Deng, Hang, Molins, Sergi, Steefel, Carl, DePaolo, Donald, Voltolini, Marco, Yang, Li, & Ajo-Franklin, Jonathan. A 2.5D Reactive Transport Model for Fracture Alteration Simulation. United States. doi:10.1021/acs.est.6b02184.
Deng, Hang, Molins, Sergi, Steefel, Carl, DePaolo, Donald, Voltolini, Marco, Yang, Li, and Ajo-Franklin, Jonathan. 2016. "A 2.5D Reactive Transport Model for Fracture Alteration Simulation". United States. doi:10.1021/acs.est.6b02184.
title = {A 2.5D Reactive Transport Model for Fracture Alteration Simulation},
author = {Deng, Hang and Molins, Sergi and Steefel, Carl and DePaolo, Donald and Voltolini, Marco and Yang, Li and Ajo-Franklin, Jonathan},
abstractNote = {},
doi = {10.1021/acs.est.6b02184},
journal = {Environmental Science and Technology},
number = 14,
volume = 50,
place = {United States},
year = 2016,
month = 7
  • Recent advances in microbiology have enabled the quantitative simulation of microbial metabolism and growth based on genome-scale characterization of metabolic pathways and fluxes. We have incorporated a genome-scale metabolic model of the iron-reducing bacteria Geobacter sulfurreducens into a pore-scale simulation of microbial growth based on coupling of iron reduction to oxidation of a soluble electron donor (acetate). In our model, fluid flow and solute transport is governed by a combination of the Navier-Stokes and advection-diffusion-reaction equations. Microbial growth occurs only on the surface of soil grains where solid-phase mineral iron oxides are available. Mass fluxes of chemical species associated withmore » microbial growth are described by the genome-scale microbial model, implemented using a constraint-based metabolic model, and provide the Robin-type boundary condition for the advection-diffusion equation at soil grain surfaces. Conventional models of microbially-mediated subsurface reactions use a lumped reaction model that does not consider individual microbial reaction pathways, and describe reactions rates using empirically-derived rate formulations such as the Monod-type kinetics. We have used our pore-scale model to explore the relationship between genome-scale metabolic models and Monod-type formulations, and to assess the manifestation of pore-scale variability (microenvironments) in terms of apparent Darcy-scale microbial reaction rates. The genome-scale model predicted lower biomass yield, and different stoichiometry for iron consumption, in comparisonto prior Monod formulations based on energetics considerations. We were able to fit an equivalent Monod model, by modifying the reaction stoichiometry and biomass yield coefficient, that could effectively match results of the genome-scale simulation of microbial behaviors under excess nutrient conditions, but predictions of the fitted Monod model deviated from those of the genome-scale model under conditions in which one or more nutrients were limiting. The fitted Monod kinetic model was also applied at the Darcy scale; that is, to simulate average reaction processes at the scale of the entire pore-scale model domain. As we expected, even under excess nutrient conditions for which the Monod and genome-scale models predicted equal reaction rates at the pore scale, the Monod model over-predicted the rates of biomass growth and iron and acetate utilization when applied at the Darcy scale. This discrepancy is caused by an inherent assumption of perfect mixing over the Darcy-scale domain, which is clearly violated in the pore-scale models. These results help to explain the need to modify the flux constraint parameters in order to match observations in previous applications of the genome-scale model at larger scales. These results also motivate further investigation of quantitative multi-scale relationships between fundamental behavior at the pore scale (where genome-scale models are appropriately applied) and observed behavior at larger scales (where predictions of reactive transport phenomena are needed).« less
  • One-dimensional transport of solutes in a saturated porous medium has been simulated using a numerical model CHEMTRN. This model includes dispersion/diffusion, advection, ion exchange, formation of complexes in the aqueous phase, and the dissociation of water. The mass action, transport, and site constraint equations are expressed in a differential/algebraic form and solved simultaneously. In the mass action equations describing complexation and sorption processes, the activities of the species in the aqueous phase are approximated using the Davies equation. A solubility constraint can be set at the inner boundary. The model was applied to a study of the migration of strontiummore » in groundwater. Results show that surface site competition is important when the strontium concentration is comparable to the background electrolyte concentration, in which case a constant k/sub D/ model does not simulate the transport of chemical species adequately. In addition, the pH level was buffered to a level of approximately 10 near the wasteform due to the incorporation of the solubility constraint. For the particular case studied, the incorporation of aqueous complexation did not significantly enhance strontium mobility. 24 references, 8 figures, 1 table.« less
  • Abstract not provided.
  • We investigated the ability of Pseudomonas stutzeri strain DCP-Ps1 to drive CaCO 3 biomineralization in a microfluidic flowcell (i.e., micromodel) that simulates subsurface porous media. Results indicate that CaCO 3 precipitation occurs during NO 3 reduction with a maximum saturation index (SI calcite) of ~1.56, but not when NO 3 was removed, inactive biomass remained, and pH and alkalinity were adjusted to SI calcite ~ 1.56. CaCO 3 precipitation was promoted by metabolically active cultures of strain DCP-Ps1, which at similar values of SIcalcite, have a more negative surface charge than inactive strain DCP-Ps1. A two-stage NO 3more » reduction (NO 3 → NO 2 → N 2) pore-scale reactive transport model was used to evaluate denitrification kinetics, which was observed in the micromodel as upper (NO 3 reduction) and lower (NO 2 reduction) horizontal zones of biomass growth with CaCO 3 precipitation exclusively in the lower zone. Our model results are consistent with two biomass growth regions and indicate that precipitation occurred in the lower zone because the largest increase in pH and alkalinity is associated with NO 2 reduction. CaCO 3 precipitates typically occupied the entire vertical depth of pores and impacted porosity, permeability, and flow. This study provides a framework for incorporating microbial activity in biogeochemistry models, which often base biomineralization only on SI (caused by biotic or abiotic reactions) and, thereby, underpredict the extent of this complex process. Furthermore, these results have wide-ranging implications for understanding reactive transport in relevance to groundwater remediation, CO 2 sequestration, and enhanced oil recovery.« less