Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We develop a 2-D pore scale model of coupled fluid flow, reactive transport, and calcium carbonate (CaCO 3) precipitation and dissolution. The model is used to simulate transient experimental results of CaCO 3 precipitation and dissolution under supersaturated conditions in a microfluidic pore network (i.e., micromodel) in order to improve understanding of coupled reactive transport systems perturbed by geological CO 2 injection. In the micromodel, precipitation is induced by transverse mixing along the centerline in pore bodies. The reactive transport model includes the impact of pH upon carbonate speciation and a CaCO 3 reaction rate constant, the effect of changing reactive surface area upon the reaction, and the impact of pore blockage from CaCO 3 precipitation on diffusion and flow. Overall, the pore scale model qualitatively captured the precipitate morphology, precipitation rate, and maximum precipitation area using parameter values from the literature. In particular, we found that proper estimation of the effective diffusion coefficient (D eff) and the reactive surface area is necessary to adequately simulate precipitation and dissolution rates. In order to match the initial phase of fast precipitation, it was necessary to consider the top and bottom of the micromodel as additional reactive surfaces. In order to match a later phase when dissolution occurred, it was necessary to increase the dissolution rate compared to the precipitation rate, but the simulated precipitate area was still higher than the experimental results after ∼30 min, highlighting the need for future study. The model presented here allows us to simulate and mechanistically evaluate precipitation and dissolution of CaCO 3 observed in a micromodel pore network. This study leads to improved understanding of the fundamental physicochemical processes of CaCO 3 precipitation and dissolution under far-from-equilibrium conditions. Copyright 2012 by the American Geophysical Union.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This paper describes the modeling efforts undertaken during a recently completed feasibility study of a generic shale repository for disposal of high-level radioactive waste within the United States. A coupled thermal-hydrological-mechanical-chemical analysis of the shale repository was performed using the SIERRA Mechanics code developed at Sandia National Laboratories. Because U.S. efforts have focused on the volcanic tuff site at Yucca Mountain, radioactive waste disposal in U.S. shale formations has not been considered for many years. However, advances in multi-physics computational modeling and research into clay mineralogy continue to improve the scientific basis for assessing nuclear waste repository performance in such formations. Disposal of high-level radioactive waste in suitable shale formations is attractive because the material is essentially impermeable and self-sealing, conditions are chemically reducing, and sorption tends to prevent radionuclide transport. Vertically and laterally extensive shale and clay formations exist in multiple locations in the contiguous 48 states. © 2011 ARMA, American Rock Mechanics Association.

Abstract not provided.

Abstract not provided.

Mudstone pore networks are strong modifiers of sedimentary basin fluid dynamics and have a critical role in the distribution of hydrocarbons and containment of injected fluids. Using core samples from continental and marine mudstones, we investigate properties of pore types and networks from a variety of geologic environments, together with estimates of capillary beam- scanning electron microscopy, suggest seven dominant mudstone pore types distinguished by geometry and connectivity. A dominant planar pore type occurs in all investigated mudstones and generally has high coordination numbers (i.e., number of neighboring connected pores). Connected networks of pores of this type contribute to high mercury capillary pressures due to small pore throats at the junctions of connected pores and likely control most matrix transport in these mudstones. Other pore types are related to authigenic (e.g., replacement or pore-lining precipitation) clay minerals and pyrite nodules; pores in clay packets adjacent to larger, more competent clastic grains; pores in organic phases; and stylolitic and microfracture-related pores. Pores within regions of authigenic clay minerals often form small isolated networks (<3 μm). Pores in stringers of organic phases occur as tubular pores or slit- and/or sheet-like pores. These form short, connected lengths in 3D reconstructions, but appear to form networks no larger than a few microns in size. Sealing efficiency of the studied mudstones increases with greater distal depositional environments and greater maximum depth of burial. © 2011 Geological Society of America.

Abstract not provided.

Abstract not provided.

This document summarizes research performed under the SNL LDRD entitled - Computational Mechanics for Geosystems Management to Support the Energy and Natural Resources Mission. The main accomplishment was development of a foundational SNL capability for computational thermal, chemical, fluid, and solid mechanics analysis of geosystems. The code was developed within the SNL Sierra software system. This report summarizes the capabilities of the simulation code and the supporting research and development conducted under this LDRD. The main goal of this project was the development of a foundational capability for coupled thermal, hydrological, mechanical, chemical (THMC) simulation of heterogeneous geosystems utilizing massively parallel processing. To solve these complex issues, this project integrated research in numerical mathematics and algorithms for chemically reactive multiphase systems with computer science research in adaptive coupled solution control and framework architecture. This report summarizes and demonstrates the capabilities that were developed together with the supporting research underlying the models. Key accomplishments are: (1) General capability for modeling nonisothermal, multiphase, multicomponent flow in heterogeneous porous geologic materials; (2) General capability to model multiphase reactive transport of species in heterogeneous porous media; (3) Constitutive models for describing real, general geomaterials under multiphase conditions utilizing laboratory data; (4) General capability to couple nonisothermal reactive flow with geomechanics (THMC); (5) Phase behavior thermodynamics for the CO2-H2O-NaCl system. General implementation enables modeling of other fluid mixtures. Adaptive look-up tables enable thermodynamic capability to other simulators; (6) Capability for statistical modeling of heterogeneity in geologic materials; and (7) Simulator utilizes unstructured grids on parallel processing computers.