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.
The National Water, Energy and Carbon Sequestration simulation model (WECSsim) is being developed to address the question, 'Where in the current and future U.S. fossil fuel based electricity generation fleet are there opportunities to couple CO{sub 2} storage and extracted water use, and what are the economic and water demand-related impacts of these systems compared to traditional power systems?' The WECSsim collaborative team initially applied this framework to a test case region in the San Juan Basin, New Mexico. Recently, the model has been expanded to incorporate the lower 48 states of the U.S. Significant effort has been spent characterizing locations throughout the U.S. where CO{sub 2} might be stored in saline formations including substantial data collection and analysis efforts to supplement the incomplete brine data offered in the NatCarb database. WECSsim calculates costs associated with CO{sub 2} capture and storage (CCS) for the power plant to saline formation combinations including parasitic energy costs of CO{sub 2} capture, CO{sub 2} pipelines, water treatment options, and the net benefit of water treatment for power plant cooling. Currently, the model can identify the least-cost deep saline formation CO{sub 2} storage option for any current or proposed coal or natural gas-fired power plant in the lower 48 states. Initial results suggest that additional, cumulative water withdrawals resulting from national scale CCS may range from 676 million gallons per day (MGD) to 30,155 MGD depending on the makeup power and cooling technologies being utilized. These demands represent 0.20% to 8.7% of the U.S. total fresh water withdrawals in the year 2000, respectively. These regional and ultimately nation-wide, bottom-up scenarios coupling power plants and saline formations throughout the U.S. can be used to support state or national energy development plans and strategies.
This report describes the progress in fiscal year 2010 in developing the Waste Integrated Performance and Safety Codes (IPSC) in support of the U.S. Department of Energy (DOE) Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Campaign. The goal of the Waste IPSC is to develop an integrated suite of computational modeling and simulation capabilities to quantitatively assess the long-term performance of waste forms in the engineered and geologic environments of a radioactive waste storage or disposal system. The Waste IPSC will provide this simulation capability (1) for a range of disposal concepts, waste form types, engineered repository designs, and geologic settings, (2) for a range of time scales and distances, (3) with appropriate consideration of the inherent uncertainties, and (4) in accordance with robust verification, validation, and software quality requirements. Waste IPSC activities in fiscal year 2010 focused on specifying a challenge problem to demonstrate proof of concept, developing a verification and validation plan, and performing an initial gap analyses to identify candidate codes and tools to support the development and integration of the Waste IPSC. The current Waste IPSC strategy is to acquire and integrate the necessary Waste IPSC capabilities wherever feasible, and develop only those capabilities that cannot be acquired or suitably integrated, verified, or validated. This year-end progress report documents the FY10 status of acquisition, development, and integration of thermal-hydrologic-chemical-mechanical (THCM) code capabilities, frameworks, and enabling tools and infrastructure.
This work uses a bifurcation approach to develop theoretical predictions for deformation band formation for a suite of true triaxial tests on Castlegate sandstone. In particular, the influence of the intermediate principal stress on strain localization is examined. Using common simplifying assumptions (localization occurs at peak stress, and the failure surface is similar to the yield surface), theoretical predictions captured the overall trends observed experimentally. However, agreement between predicted and observed band orientations for individual specimens was varied. This highlights the importance of detailed data analyses to accurately determine key material parameter values at the inception of localization.
Mudstone mechanical testing is often limited by poor core recovery and sample size, preservation and preparation issues, which can lead to sampling bias, damage, and time-dependent effects. A micropillar compression technique, originally developed by Uchic et al. 2004, here is applied to elasto-plastic deformation of small volumes of mudstone, in the range of cubic microns. This study examines behavior of the Gothic shale, the basal unit of the Ismay zone of the Pennsylvanian Paradox Formation and potential shale gas play in southeastern Utah, USA. Precision manufacture of micropillars 5 microns in diameter and 10 microns in length are prepared using an ion-milling method. Characterization of samples is carried out using: dual focused ion - scanning electron beam imaging of nano-scaled pores and distribution of matrix clay and quartz, as well as pore-filling organics; laser scanning confocal (LSCM) 3D imaging of natural fractures; and gas permeability, among other techniques. Compression testing of micropillars under load control is performed using two different nanoindenter techniques. Deformation of 0.5 cm in diameter by 1 cm in length cores is carried out and visualized by a microscope loading stage and laser scanning confocal microscopy. Axisymmetric multistage compression testing and multi-stress path testing is carried out using 2.54 cm plugs. Discussion of results addresses size of representative elementary volumes applicable to continuum-scale mudstone deformation, anisotropy, and size-scale plasticity effects. Other issues include fabrication-induced damage, alignment, and influence of substrate.
Subsurface containment of CO2 is predicated on effective caprock sealing. Many previous studies have relied on macroscopic measurements of capillary breakthrough pressure and other petrophysical properties without direct examination of solid phases that line pore networks and directly contact fluids. However, pore-lining phases strongly contribute to sealing behavior through interfacial interactions among CO2, brine, and the mineral or non-mineral phases. Our high resolution (i.e., sub-micron) examination of the composition of pore-lining phases of several continental and marine mudstones indicates that sealing efficiency (i.e., breakthrough pressure) is governed by pore shapes and pore-lining phases that are not identifiable except through direct characterization of pores. Bulk X-ray diffraction data does not indicate which phases line the pores and may be especially lacking for mudstones with organic material. Organics can line pores and may represent once-mobile phases that modify the wettability of an originally clay-lined pore network. For shallow formations (i.e., < {approx}800 m depth), interfacial tension and contact angles result in breakthrough pressures that may be as high as those needed to fracture the rock - thus, in the absence of fractures, capillary sealing efficiency is indicated. Deeper seals have poorer capillary sealing if mica-like wetting dominates the wettability. We thank the U.S. Department of Energy's National Energy Technology Laboratory and the Office of Basic Energy Sciences, and the Southeast and Southwest Carbon Sequestration Partnerships for supporting this work.
Injection of CO2 into underground rock formations can reduce atmospheric CO2 emissions. Caprocks present above potential storage formations are the main structural trap inhibiting CO2 from leaking into overlying aquifers or back to the Earth's surface. Dissolution and precipitation of caprock minerals resulting from reaction with CO2 may alter the pore network where many pores are of the micrometer to nanometer scale, thus altering the structural trapping potential of the caprock. However, the distribution, geometry and volume of pores at these scales are poorly characterized. In order to evaluate the overall risk of leakage of CO2 from storage formations, a first critical step is understanding the distribution and shape of pores in a variety of different caprocks. As the caprock is often comprised of mudstones, we analyzed samples from several mudstone formations with small angle neutron scattering (SANS) and high-resolution transmission electron microscopy (TEM) imaging to compare the pore networks. Mudstones were chosen from current or potential sites for carbon sequestration projects including the Marine Tuscaloosa Group, the Lower Tuscaloosa Group, the upper and lower shale members of the Kirtland Formation, and the Pennsylvanian Gothic shale. Expandable clay contents ranged from 10% to approximately 40% in the Gothic shale and Kirtland Formation, respectively. During SANS, neutrons effectively scatter from interfaces between materials with differing scattering length density (i.e., minerals and pores). The intensity of scattered neutrons, I(Q), where Q is the scattering vector, gives information about the volume and arrangement of pores in the sample. The slope of the scattering data when plotted as log I(Q) vs. log Q provides information about the fractality or geometry of the pore network. On such plots slopes from -2 to -3 represent mass fractals while slopes from -3 to -4 represent surface fractals. Scattering data showed surface fractal dimensions for the Kirtland formation and one sample from the Tuscaloosa formation close to 3, indicating very rough surfaces. In contrast, scattering data for the Gothic shale formation exhibited mass fractal behavior. In one sample of the Tuscaloosa formation the data are described by a surface fractal at low Q (larger pores) and a mass fractal at high Q (smaller pores), indicating two pore populations contributing to the scattering behavior. These small angle neutron scattering results, combined with high-resolution TEM imaging, provided a means for both qualitative and quantitative analysis of the differences in pore networks between these various mudstones.
Mercury intrusion porosimetry (MIP) is an often-applied technique for determining pore throat distributions and seal analysis of fine-grained rocks. Due to closure effects, potential pore collapse, and complex pore network topologies, MIP data interpretation can be ambiguous, and often biased toward smaller pores in the distribution. We apply 3D imaging techniques and lattice-Boltzmann modeling in interpreting MIP data for samples of the Cretaceous Selma Group Chalk. In the Mississippi Interior Salt Basin, the Selma Chalk is the apparent seal for oil and gas fields in the underlying Eutaw Fm., and, where unfractured, the Selma Chalk is one of the regional-scale seals identified by the Southeast Regional Carbon Sequestration Partnership for CO2 injection sites. Dual focused ion - scanning electron beam and laser scanning confocal microscopy methods are used for 3D imaging of nanometer-to-micron scale microcrack and pore distributions in the Selma Chalk. A combination of image analysis software is used to obtain geometric pore body and throat distributions and other topological properties, which are compared to MIP results. 3D data sets of pore-microfracture networks are used in Lattice Boltzmann simulations of drainage (wetting fluid displaced by non-wetting fluid via the Shan-Chen algorithm), which in turn are used to model MIP procedures. Results are used in interpreting MIP results, understanding microfracture-matrix interaction during multiphase flow, and seal analysis for underground CO2 storage.
Shales and other mudstones are the most abundant rock types in sedimentary basins, yet have received comparatively little attention. Common as hydrocarbon seals, these are increasingly being targeted as unconventional gas reservoirs, caprocks for CO2 sequestration, and storage repositories for waste. The small pore and grain size, large specific surface areas, and clay mineral structures lend themselves to rapid reaction rates, high capillary pressures, and semi-permeable membrane behavior accompanying changes in stress, pressure, temperature and chemical conditions. Under far from equilibrium conditions, mudrocks display a variety of spatio-temporal self-organized phenomena arising from nonlinear thermo-mechano-chemo-hydro coupling. Beginning with a detailed examination of nano-scale pore network structures in mudstones, we discuss the dynamics behind such self-organized phenomena as pressure solitons in unconsolidated muds, chemically-induced flow self focusing and permeability transients, localized compaction, time dependent well-bore failure, and oscillatory osmotic fluxes as they occur in clay-bearing sediments. Examples are draw from experiments, numerical simulation, and the field. These phenomena bear on the ability of these rocks to serve as containment barriers.
Dissolved CO2 from geological CO2 sequestration may react with dissolved minerals in fractured rocks or confined aquifers and cause mineral precipitation. The overall rate of reaction can be limited by diffusive or dispersive mixing, and mineral precipitation can block pores and further hinder these processes. Mixing-induced calcite precipitation experiments were performed by injecting solutions containing CaCl2 and Na2CO3 through two separate inlets of a micromodel (1-cm x 2-cm x 40-microns); transverse dispersion caused the two solutions to mix along the center of the micromodel, resulting in calcite precipitation. The amount of calcite precipitation initially increased to a maximum and then decreased to a steady state value. Fluorescent microscopy and imaging techniques were used to visualize calcite precipitation, and the corresponding effects on the flow field. Experimental micromodel results were evaluated with pore-scale simulations using a 2-D Lattice-Boltzmann code for water flow and a finite volume code for reactive transport. The reactive transport model included the impact of pH upon carbonate speciation and calcite dissolution. We found that proper estimation of the effective diffusion coefficient and the reaction surface area is necessary to adequately simulate precipitation and dissolution rates. The effective diffusion coefficient was decreased in grid cells where calcite precipitated, and keeping track of reactive surface over time played a significant role in predicting reaction patterns. Our results may improve understanding of the fundamental physicochemical processes during CO2 sequestration in geologic formations.