In the supercritical CO2-water-mineral systems relevant to subsurface CO2 sequestration, interfacial processes at the supercritical fluid-mineral interface will strongly affect core- and reservoir-scale hydrologic properties. Experimental and theoretical studies have shown that water films will form on mineral surfaces in supercritical CO2, but will be thinner than those that form in vadose zone environments at any given matric potential. The theoretical model presented here allows assessment of water saturation as a function of matric potential, a critical step for evaluating relative permeabilities the CO2 sequestration environment. The experimental water adsorption studies, using Quartz Crystal Microbalance and Fourier Transform Infrared Spectroscopy methods, confirm the major conclusions of the adsorption/condensation model. Additional data provided by the FTIR study is that CO2 intercalation into clays, if it occurs, does not involve carbonate or bicarbonate formation, or significant restriction of CO2 mobility. We have shown that the water film that forms in supercritical CO2 is reactive with common rock-forming minerals, including albite, orthoclase, labradorite, and muscovite. The experimental data indicate that reactivity is a function of water film thickness; at an activity of water of 0.9, the greatest extent of reaction in scCO2 occurred in areas (step edges, surface pits) where capillary condensation thickened the water films. This suggests that dissolution/precipitation reactions may occur preferentially in small pores and pore throats, where it may have a disproportionately large effect on rock hydrologic properties. Finally, a theoretical model is presented here that describes the formation and movement of CO2 ganglia in porous media, allowing assessment of the effect of pore size and structural heterogeneity on capillary trapping efficiency. The model results also suggest possible engineering approaches for optimizing trapping capacity and for monitoring ganglion formation in the subsurface.
This report summarizes the work completed under the Laboratory Directed Research and Development (LDRD) project 09-1351, 'Computational Investigation of Thermal Gas Separation for CO{sub 2} Capture'. Thermal gas separation for a binary mixture of carbon dioxide and nitrogen is investigated using the Direct Simulation Monte Carlo (DSMC) method of molecular gas dynamics. Molecular models for nitrogen and carbon dioxide are developed, implemented, compared to theoretical results, and compared to several experimental thermophysical properties. The molecular models include three translational modes, two fully excited rotational modes, and vibrational modes, whose degree of excitation depends on the temperature. Nitrogen has one vibrational mode, and carbon dioxide has four vibrational modes (two of which are degenerate). These models are used to perform a parameter study for mixtures of carbon dioxide and nitrogen confined between parallel walls over realistic ranges of gas temperatures and nominal concentrations of carbon dioxide. The degree of thermal separation predicted by DSMC is slightly higher than experimental values and is sensitive to the details of the molecular models.