Modeling CO₂–Water–Mineral Wettability and Mineralization for Carbon Geosequestration

Abstract

Carbon dioxide (CO) capture and storage (CCS) is an important climate change mitigation option along with improved energy efficiency, renewable energy, and nuclear energy. CO geosequestration, that is, to store CO under the subsurface of Earth, is feasible because the world's sedimentary basins have high capacity and are often located in the same region of the world as emission sources. How CO interacts with the connate water and minerals is the focus of this Account. There are four trapping mechanisms that keep CO in the pores of subsurface rocks: (1) structural trapping, (2) residual trapping, (3) dissolution trapping, and (4) mineral trapping. The first two are dominated by capillary action, where wettability controls CO and water two-phase flow in porous media. We review state-of-the-art studies on CO/water/mineral wettability, which was found to depend on pressure and temperature conditions, salt concentration in aqueous solutions, mineral surface chemistry, and geometry. We then review some recent advances in mineral trapping. First, we show that it is possible to reproduce the CO/water/mineral wettability at a wide range of pressures using molecular dynamics (MD) simulations. As the pressure increases, CO gas transforms into a supercritical fluid or liquid at ∼7.4 MPa depending on the environmental temperature. This transition leads to a substantial decrease of the interfacial tension between CO and reservoir brine (or pure water). However, the wettability of CO/water/rock systems depends on the type of rock surface. Recently, we investigated the contact angle of CO/water/silica systems with two different silica surfaces using MD simulations. We found that contact angle increased with pressure for the hydrophobic (siloxane) surface while it was almost constant for the hydrophilic (silanol) surface, in excellent agreement with experimental observations. Furthermore, we found that the CO thin films at the CO-hydrophilic silica and CO-HO interfaces displayed a linear correlation, which can in turn explain the constant contact angle on the hydrophilic silica surface. In view of the literature and our study results, a few recommendations seem necessary to construct a molecular system suitable to study wettability with MD simulations. Future work should be conducted to determine the influence of brine salinity on the wettability of minerals with high cation exchange capacity. Mineral trapping is believed to be an extremely slow process, likely taking thousands of years. However, a recent pilot study demonstrated that CO mineralization occurs within 2 years in highly reactive basalt reservoirs. A first-principles MD study has also shown that carbonation reactions occur rapidly at the surface oxygen sites of a reactive mineral. We observed carbonate ions on both a newly cleaved quartz surface (without hydrolysis), and a basalt andesine surface after hydrolysis in a CO-rich environment. Future work should consider the influence of water, gas impurities, and mineral cation type on carbonation.

Wetting Properties of the CO2–Water–Calcite System via Molecular Simulations: Shape and Size Effects

“…According to Li et al, 28 the deviation at high pressures does not result from density deviation but can be due to other simulation parameters such as the system size. 32 (open blue pentagons) and Liang et al 36 (closed blue pentagons), 323 K, this study (red diamonds), experimental data from Chiquet et al 10 (closed green diamonds) and Kvamme et al 26 (open green diamonds) and simulation data from Liu et al 29 (blue diamonds) and 383 K, this study (red circles), experimental data from Chiquet et al 10 (green circles) and simulations results of Nielsen et al 8 (open blue circles) and Liang et al 36 (closed blue circles). The black symbols represent the water surface tension calculated using the IAPWS recommended interpolating equation for the surface tension of ordinary substances 71 at the given temperature.…”

Crystallization in Confinement

“…The crystallization of minerals in heterogeneous porous media is a widespread phenomenon in geochemistry and building materials, and is observed in a wide range of settings and environments. [350][351][352] With relevance to processes including the degradation of rocks in natural systems, and of porous building materials such as stone, brick, cement and concrete, [11] the formation of ore deposits, [353] the storage of carbon dioxide [13,14] and sequestration of toxic metals in minerals, [15] and the extraction of oil, [354][355][356][357] this topic attracts enormous interest and is the subject of an extensive body of literature. Predicting the outcome of crystallization within these porous media is an extremely complex problem that depends on an ability to model both transport processes and crystal nucleation and growth within confined volumes.…”

Wetting Properties of the CO2–Water–Calcite System via Molecular Simulations: Shape and Size Effects

“…According to Li et al, 28 the deviation at high pressures does not result from density deviation but can be due to other simulation parameters such as the system size. 36 (closed blue circles). The black symbols represent the water surface tension calculated using the IAPWS recommended interpolating equation for the surface tension of ordinary substances 71 at the given temperature.…”

“…12,16,17 Classical molecular dynamics (MD) simulation has proven to be a powerful tool for investigating wetting phenomena in CO 2 -watermineral systems. MD has been used to predict the interfacial tension between CO 2 and water (or saline solution) 8,[25][26][27][28][29][30][31][32][33][34][35][36] and the contact angle of CO 2 -water-mineral surfaces 27,30,32,[35][36][37][38][39][40][41][42] as a function of pressure, temperature and salinity. However, while the accuracy of the predicted interfacial tension only depends on the ability of the force fields to reproduce the interfacial energies, the contact angle determination presents a few more challenges.…”