The Nanoscale Basis of CO2 Trapping for Geologic Storage

Nanoscale Understanding on CO2 Diffusion and Adsorption in Clay Matrix Nanopores: Implications for Carbon Geosequestration

Carbon capture and storage (CCS) in subsurface reservoirs represents a highly promising and viable strategy for mitigating global carbon emissions. In the context of CCS implementation, it is particularly crucial to understand the complex molecular diffusive and adsorptive behaviors of anthropogenic carbon dioxide (CO2) in the subsurface at the nanoscale. Yet, conventional molecular models typically represent only single-slit pores and overlook the complexity of interconnected nanopores. In this work, finite kaolinite lamellar assemblages with abundant nanopores (r < 2 nm) were used. Molecular dynamics simulations were performed to quantify the spatial distribution correlations, adsorption preference, diffusivity, and residence time of the CO2 molecules in kaolinite nanopores. The movement of the CO2 molecules primarily occurs in the central and proximity regions of the siloxane surfaces, progressing from larger to smaller nanopores. CO2 prefers smaller nanopores over larger ones. The diffusion coefficients increase, while residence times decrease, with the pore size increasing, differing from typical slit-pore models due to the pore shape and interconnectivity. The perspectives in this study, which would be challenging in conventional slit-pore models, will facilitate our comprehension of the CO2 molecular behaviors in the complex subsurface clay sediments for developing quantitative estimation techniques throughout the CCS project durations.

Artificial subsurface lithoautotrophic microbial ecosystems and gas storage in deep subsurface

Over the next few years, it is planned to convert all or part of the underground gas storage (UGS) facilities used for natural gas (salt caverns, depleted hydrocarbon reservoirs, and deep aquifers) into underground dihydrogen (H2) storage reservoirs. These deep environments host microbial communities, some of which are hydrogenotrophic (sulfate reducers, acetogens, and methanogens). The current state of microbiological knowledge is thus presented for the three types of UGS facilities. In the mid-1990s, the concept of anaerobic subsurface lithoautotrophic microbial ecosystems, or SLiMEs, emerged. It is expected that the large-scale injection of H2 into subsurface environments will generate new microbial ecosystems called artificial SLiMEs, which could persist over time. These artificial SLiMEs could lead to H2 loss, an intense methanogenic activity, a degradation of gas quality and a risk to installations through sulfide production. However, recent studies on salt caverns and deep aquifers suggest that hydrogenotrophic microbial activity also leads to alkalinization (up to pH 10), which can constrain hydrogenotrophy. Therefore, studying and understanding these artificial SLiMEs is both a necessity for the development of the H2 industry and presents an opportunity for ecologists to monitor the evolution of deep environments in real time.

Phase quantification using deep neural network processing of XRD patterns

Mineral identification and quantification are key to the understanding and, hence, the capacity to predict material properties. The method of choice for mineral quantification is powder X-ray diffraction (XRD), generally using a Rietveld refinement approach. However, a successful Rietveld refinement requires preliminary identification of the phases that make up the sample. This is generally carried out manually, and this task becomes extremely long or virtually impossible in the case of very large datasets such as those from synchrotron X-ray diffraction computed tomography. To circumvent this issue, this article proposes a novel neural network (NN) method for automating phase identification and quantification. An XRD pattern calculation code was used to generate large datasets of synthetic data that are used to train the NN. This approach offers significant advantages, including the ability to construct databases with a substantial number of XRD patterns and the introduction of extensive variability into these patterns. To enhance the performance of the NN, a specifically designed loss function for proportion inference was employed during the training process, offering improved efficiency and stability compared with traditional functions. The NN, trained exclusively with synthetic data, proved its ability to identify and quantify mineral phases on synthetic and real XRD patterns. Trained NN errors were equal to 0.5% for phase quantification on the synthetic test set, and 6% on the experimental data, in a system containing four phases of contrasting crystal structures (calcite, gibbsite, dolomite and hematite). The proposed method is freely available on GitHub and allows for major advances since it can be applied to any dataset, regardless of the mineral phases present.

Water Film Structure and Wettability of Different Quartz Surfaces: Hydrogen Bonding Across Various Cutting Planes

Quartz is ubiquitous in subsurface formations. The crystal faces have different atomic arrangements. Knowledge of the molecular structures on the surface of quartz is key in many processes. Molecular dynamics simulations are conducted to investigate the atomic arrangement effect on the water film structure, ion adsorption, and wettability at three different α-quartz surfaces. The interfacial structures depend on the quartz surface. Intrasurface hydrogen bonding between surface silanols differs in the α-quartz surface. At the (0001) surface, the OH density is 9.58 nm-2, consisting of Q2 units with two hydroxyl groups per silicone atom. At the (101̅0)-β surface, the OH density is 7.54 nm-2, consisting of Q3 units with one hydroxyl group per silicone atom; there is significant intrasurface hydrogen bonding. At the (101̅0)-α surface, the OH density is 7.54 nm-2, consisting of Q2 units; however, there is little intrasurface hydrogen bonding. The intrasurface hydrogen bonding results in the exposure of hydrogen-bond acceptors to the aqueous phase, causing water molecules to have an H-up (hydrogen toward surface) orientation. This orientation can be found at the (0001) and (101̅0)-β surfaces; it is related to the degree of ordering at the surface. The ordering at the (0001) and (101̅0)-β surfaces is higher than that at the (101̅0)-α surface. In aqueous systems with ions, cation adsorption is the most dominant at the (0001) surface due to the largest surface density of the intrasurface hydrogen bonding, providing interaction sites for cations to be adsorbed. We observe a pronounced decrease in water film thickness from the ions at the (0001) surface only, likely due to significant cation adsorption. In this work, we demonstrate that the hydrogen-bond network, which varies from the plane cut, affects the water film structure and ion adsorption. The contact is nearly zero despite the changes in the film thickness and molecular structure at the temperature of 318 K.

Effect of the Temperature on Interfacial Properties of CO2/H2 Mixtures Contacting with Brine and Hydrophilic Silica by Molecular Dynamics Simulations

Underground H2 storage (UHS) is a promising technology to achieve large-scale, long-term H2 storage. Using CO2 as a cushion gas to maintain the pressure of the reservoir and withdraw stored H2 in the saline aquifer simultaneously enables the implementation of UHS and underground CO2 storage (UCS). The difference in the molecular properties of CO2 and H2 leads to distinct interfacial behavior when in contact with the brine and rock, thereby affecting the flow patterns and trapping mechanisms of gases in geological formations. Accurate prediction of the interfacial properties of CO2, H2, and the mixtures when interacting with brine and rock is crucial to minimizing the uncertainties in UHS and UCS projects. In this study, molecular dynamics (MD) simulations are performed to predict the interfacial tension, surface excess, bubble evolution, and contact angle of CO2, H2, and the mixtures at 10 MPa and 300-400 K. The MD results show that the interaction of CO2 with H2O and hydrophilic silica is considerably stronger than that of H2. The interfacial tension reduces linearly with the temperature in H2-dominated mixture systems, and the surface adsorption of H2 can diminish in a CO2-dominated system or at high-temperature conditions. The hydrophilic silica is more CO2-wet than H2-wet, and the attached CO2 bubble is more easily disconnected. Ions and the temperature play different roles in the contact angle.

Oxide- and Silicate-Water Interfaces and Their Roles in Technology and the Environment

Interfacial reactions drive all elemental cycling on Earth and play pivotal roles in human activities such as agriculture, water purification, energy production and storage, environmental contaminant remediation, and nuclear waste repository management. The onset of the 21st century marked the beginning of a more detailed understanding of mineral aqueous interfaces enabled by advances in techniques that use tunable high-flux focused ultrafast laser and X-ray sources to provide near-atomic measurement resolution, as well as by nanofabrication approaches that enable transmission electron microscopy in a liquid cell. This leap into atomic- and nanometer-scale measurements has uncovered scale-dependent phenomena whose reaction thermodynamics, kinetics, and pathways deviate from previous observations made on larger systems. A second key advance is new experimental evidence for what scientists hypothesized but could not test previously, namely, interfacial chemical reactions are frequently driven by "anomalies" or "non-idealities" such as defects, nanoconfinement, and other nontypical chemical structures. Third, progress in computational chemistry has yielded new insights that allow a move beyond simple schematics, leading to a molecular model of these complex interfaces. In combination with surface-sensitive measurements, we have gained knowledge of the interfacial structure and dynamics, including the underlying solid surface and the immediately adjacent water and aqueous ions, enabling a better definition of what constitutes the oxide- and silicate-water interfaces. This critical review discusses how science progresses from understanding ideal solid-water interfaces to more realistic systems, focusing on accomplishments in the last 20 years and identifying challenges and future opportunities for the community to address. We anticipate that the next 20 years will focus on understanding and predicting dynamic transient and reactive structures over greater spatial and temporal ranges as well as systems of greater structural and chemical complexity. Closer collaborations of theoretical and experimental experts across disciplines will continue to be critical to achieving this great aspiration.

Influence of stearic acid and alumina nanofluid on CO2 wettability of calcite substrates: Implications for CO2 geological storage in carbonate reservoirs

HYPOTHESIS

Atmospheric CO₂ emissions trigger global warming and climate change challenges. Thus, geological CO₂ storage appears to be the most viable choice to mitigate CO₂ emissions in the atmosphere. However, the adsorption capacity of reservoir rock in the presence of diverse geological conditions, including organic acids, temperature, and pressure, can cause reduced certainty for CO₂ storage and injection problems. Wettability is critical in measuring the adsorption behavior of rock in various reservoir fluids and conditions.

EXPERIMENT

We systematically evaluated the CO₂-wettability of calcite substrates at geological conditions (323 K and 0.1, 10, and 25 MPa) in the presence of stearic acid (a replicate realistic reservoir organic material contamination). Similarly, to reverse the effects of organics on wettability, we treated calcite substrates with various alumina nanofluid concentrations (0.05, 0.1, 0.25, and 0.75 wt%) and evaluated the CO₂-wettability of calcite substrates at similar geological conditions.

FINDINGS

Stearic acid profoundly affects the contact angle of calcite substrates where wettability shifts from intermediate to CO₂-wet conditions, reducing the CO₂ geological storage potential. The treatment of organic acid-aged calcite substrates with alumina nanofluid reversed the wettability to a more hydrophilic state, increasing CO₂ storage certainty. Further, the optimum concentration displaying the optimum potential for changing the wettability in organic acid-aged calcite substrates was 0.25 wt%. The effect of organics and nanofluids should be augmented to improve the feasibility of CO₂ geological projects at the industrial scale for reduced containment security.

Assessment of the oxygen reactivity in a gas storage facility by multiphase reactive transport modeling of field data for air injection into a sandstone reservoir in the Paris Basin, France

The first objective of this study is to present unique field data on a three-year pilot test during which air containing 8 mol% O₂(g) was injected as a cushion gas into a natural gas reservoir, a carbonate-cemented sandstone aquifer located in the Paris Basin (France). 10-year system survey showed that: the oxygen was fully depleted several months after injection completion, meanwhile CO₂(g) was detected around 2-6 mol%; the pH decreased from 8 to 6, while reducing conditions shifted to mildly oxidizing ones with increasing concentration of sulfates in equilibrium with gypsum. 3 years after injection completion, the pH gradually returned to its near initial state and sulfates were reduced by 2 to 3 times. The second objective is to develop a multiphase reactive transport model based on the field data. Simulations were constructed using the HYTEC reactive transport code, progressing from 0D-batch to 2D-reservoir configurations. The model reproduced the gas-water-rock reactive sequence: 1/ full depletion of the injected O₂(g) due to pyrite oxidation, 2/ leading to acidity production and dissolved sulfates, 3/ acidity buffering by calcite dissolution, 4/ followed by gypsum precipitation and CO₂(g) exsolution. The model demonstrated that pyrite kinetics was the most significant factor governing not only the amount of O₂(g), CO₂(g) and dissolved minerals, but also the spatial extent of these chemical reactions and, hence, the gas spread inside the reservoir. The formulated advective Damköhler number for oxygen consumption indicated advection- and reaction-dominant regimes explaining the gas composition and extension. The developed field-based model could be used as a workflow for other gas storage facilities, e.g. biomethane, compressed air, and CO₂. For underground biomethane storage, the O₂(g) contents recommended in Europe, i.e. the EASEE-gas specification 2005-001-02, should have a low impact on gas composition and reservoir geochemistry when the reservoir contains efficient pH-buffers such as calcite.

Quantifying microstructures of earth materials using higher-order spatial correlations and deep generative adversarial networks

The key to most subsurface processes is to determine how structural and topological features at small length scales, i.e., the microstructure, control the effective and macroscopic properties of earth materials. Recent progress in imaging technology has enabled us to visualise and characterise microstructures at different length scales and dimensions. However, one limitation of these technologies is the trade-off between resolution and sample size (or representativeness). A promising approach to this problem is image reconstruction which aims to generate statistically equivalent microstructures but at a larger scale and/or additional dimension. In this work, a stochastic method and three generative adversarial networks (GANs), namely deep convolutional GAN (DCGAN), Wasserstein GAN with gradient penalty (WGAN-GP), and StyleGAN2 with adaptive discriminator augmentation (ADA), are used to reconstruct two-dimensional images of two hydrothermally rocks with varying degrees of complexity. For the first time, we evaluate and compare the performance of these methods using multi-point spatial correlation functions-known as statistical microstructural descriptors (SMDs)-ultimately used as external tools to the loss functions. Our findings suggest that a well-trained GAN can reconstruct higher-order, spatially-correlated patterns of complex earth materials, capturing underlying structural and morphological properties. Comparing our results with a stochastic reconstruction method based on a two-point correlation function, we show the importance of coupling training/assessment of GANs with higher-order SMDs, especially in the case of complex microstructures. More importantly, by quantifying original and reconstructed microstructures via different GANs, we highlight the interpretability of these SMDs and show how they can provide valuable insights into the spatial patterns in the synthetic images, allowing us to detect common artefacts and failure cases in training GANs.

Molecular-scale mechanisms of CO2 mineralization in nanoscale interfacial water films

The calamitous impacts of unabated carbon emission from fossil-fuel-burning energy infrastructure call for accelerated development of large-scale CO₂ capture, utilization and storage technologies that are underpinned by a fundamental understanding of the chemical processes at a molecular level. In the subsurface, rocks rich in divalent metals can react with CO₂, permanently sequestering it in the form of stable metal carbonate minerals, with the CO₂-H₂O composition of the post-injection pore fluid acting as a primary control variable. In this Review, we discuss mechanistic reaction pathways for aqueous-mediated carbonation with carbon mineralization occurring in nanoscale adsorbed water films. In the extreme of pores filled with a CO₂-dominant fluid, carbonation reactions are confined to angstrom to nanometre-thick water films coating mineral surfaces, which enable metal cation release, transport, nucleation and crystallization of metal carbonate minerals. Although seemingly counterintuitive, laboratory studies have demonstrated facile carbonation rates in these low-water environments, for which a better mechanistic understanding has come to light in recent years. The overarching objective of this Review is to delineate the unique underlying molecular-scale reaction mechanisms that govern CO₂ mineralization in these reactive and dynamic quasi-2D interfaces. We highlight the importance of understanding unique properties in thin water films, such as how water dielectric properties, and consequently ion solvation and hydration behaviour, can change under nanoconfinement. We conclude by identifying important frontiers for future work and opportunities to exploit these fundamental chemical insights for decarbonization technologies in the twenty-first century.

Study on the Alteration of Pore Parameters of Shale with Different Natural Fractures under Supercritical Carbon Dioxide Seepage

Supercritical CO2 can reduce formation fracture pressure, form more complex fractures in the near-well zone, and replace methane to complete carbon sequestration, which is an important direction for the efficient development of deep shale gas with carbon sequestration. In this paper, based on the scCO2 fracturing field test parameters and the characteristics of common shale calcite filled natural fractures, we simulated the porosity change in shale with three kinds of fractures (no fracture, named NF; axial natural fracture, named AF; and transversal natural fracture, named TF) under scCO2 seepage, and carried out the experimental verification of shale under supercritical CO2 seepage. It was found that: (1) At the same pressure, when the temperature is greater than the critical temperature, the shale porosity of three kinds of fractures gradually increases with the injection of CO2, and the higher the temperature, the more obvious the increase in porosity. (2) At the same temperature and different pressures, the effect of pressure change on the porosity of shale specimens was more obvious than that of temperature. (3) Multi-field coupling experiments of shale under supercritical CO2 seepage revealed that the porosity of all three shale specimens at the same temperature and pressure increased after CO2 injection, and the relative increase in shale porosity measured experimentally was basically consistent with the numerical simulation results. This paper reveals the mechanism of the effect of different temperatures and pressures of scCO2 and different natural fractures on the change in shale porosity, which can be used to optimize the CO2 injection in supercritical CO2 fracturing and carbon sequestration.

Physical characterization of fault rocks within the Opalinus Clay formation

Near-surface disposal of radioactive waste in shales is a promising option to safeguard the population and environment. However, natural faults intersecting these geological formations can potentially affect the long-term isolation of the repositories. This paper characterizes the physical properties and mineralogy of the internal fault core structure intersecting the Opalinus Clay formation, a host rock under investigation for nuclear waste storage at the Mont Terri Laboratory (Switzerland). We have performed porosity, density, microstructural and mineralogical measurements in different sections of the fault, including intact clays, scaly clays and fault gouge. Mercury intrusion porosimetry analysis reveal a gouge that has a pore network dominated by nanopores of less than 10 nm, yet a high-porosity (21%) and low grain density (2.62 g/cm³) when compared to the intact rock (14.2%, and 2.69 g/cm³). Thus, a more permeable internal fault core structure with respect to the surrounding rock is deduced. Further, we describe the OPA fault gouge as a discrete fault structure having the potential to act as a preferential, yet narrow, and localized channel for fluid-flow if compared to the surrounding rock. Since the fault gouge is limited to a millimetres-thick structure, we expect the barrier property of the geological formation is almost not affected.

Thin Water Films Enable Low-Temperature Magnesite Growth Under Conditions Relevant to Geologic Carbon Sequestration

Injecting supercritical CO₂ (scCO₂) into basalt formations for long-term storage is a promising strategy for mitigating CO₂ emissions. Mineral carbonation can result in permanent entrapment of CO₂; however, carbonation kinetics in thin H₂O films in humidified scCO₂ is not well understood. We investigated forsterite (Mg₂SiO₄) carbonation to magnesite (MgCO₃) via amorphous magnesium carbonate (AMC; MgCO₃·xH₂O, 0.5 < x < 1), with the goal to establish the fundamental controls on magnesite growth rates at low H₂O activity and temperature. Experiments were conducted at 25, 40, and 50 °C in 90 bar CO₂ with a H₂O film thickness on forsterite that averaged 1.78 ± 0.05 monolayers. In situ infrared spectroscopy was used to monitor forsterite dissolution and the growth of AMC, magnesite, and amorphous SiO₂ as a function of time. Geochemical kinetic modeling showed that magnesite was supersaturated by 2 to 3 orders of magnitude and grew according to a zero-order rate law. The results indicate that the main drivers for magnesite growth are sustained high supersaturation coupled with low H₂O activity, a combination of thermodynamic conditions not attainable in bulk aqueous solution. This improved understanding of reaction kinetics can inform subsurface reactive transport models for better predictions of CO₂ fate and transport.

Swelling layered minerals applications: A solid state NMR overview

Swelling layered clay minerals form an important sub-group of the phyllosilicate family. They are characterized by their ability to expand or contract in the presence or absence of water. This property makes them useful for a variety of applications, ranging from environmental technologies to heterogeneous catalysis, and including pharmaceutical and industrial applications. Solid State Nuclear Magnetic Resonance (SS-NMR) has been extensively applied in the characterization of these materials, providing useful information on their dynamics and structure that is inaccessible using other characterization methods such as X-ray diffraction. In this review, we present the key contributions of SS-NMR to the understanding of the mechanisms that govern some of the main applications associated to swelling clay minerals. The article is divided in two parts. The first part presents SS-NMR conventional applications to layered clay minerals, while the second part comprises an in-depth review of the information that SS-NMR can provide about the different properties of swelling layered clay minerals.

Enhanced Sorption of Supercritical CO2 and CH4 in the Hydrated Interlayer Pores of Smectite

Understanding the long-term confinement of supercritical fluids in the clay pores of subsurface rocks is important for many geo-energy technologies, including geological CO₂ storage. However, the adsorption properties of hydrated clay minerals remain largely uncertain because competitive adsorption experiments of supercritical fluids in the presence of water are difficult. Here, we report on the sorption properties of four source clay minerals-Ca-rich montmorillonite (STx-1b), Na-rich montmorillonite (SWy-2), illite-smectite mixed layer (ISCz-1), and illite (IMt-2)-for water at 20 °C up to relative humidity of 0.9. The measurements unveil the unsuitability of physisorption analysis by N₂ (at 77 K) and Ar (at 87 K) gases to quantify the textural properties of clays because of their inability to probe the interlayers. We further measure the sorption of CO₂ and CH₄ on swelling STx-1b and nonswelling IMt-2, both in the absence (dehydrated at 200 °C) and the presence of sub-1W preadsorbed water (following dehydration) up to 170 bar at 50 °C. We observe enhanced sorption of CO₂ and CH₄ in STx-1b (50 and 65% increase at 30 bar relative to dry STx-1b, respectively), while their adsorption on IMt-2 remains unchanged, indicating the absence of competition with water. By describing the supercritical adsorption isotherms on hydrated STx-1b with the lattice density functional theory model, we estimate that the pore volume has expanded by approximately 6% through the formation of sub-nanometer pore space. By presenting a systematic approach of quantifying the smectite clay mineral's hydrated state, this study provides an explanation for the conflicting literature observations of gas uptake capacities in the presence of water.

Assessment of Geochemical Limitations to Utilizing CO2 as a Cushion Gas in Compressed Energy Storage Systems

Compressed energy storage (CES) of air, CO₂, or H₂ in porous formations is a promising means of energy storage to abate the intermittency of renewable energy production. During operation, gas is injected during times of excess energy production and extracted during excess demands to drive turbines. Storage in saline aquifers using CO₂ as a cushion or working gas has numerous advantages over typical air storage in caverns. However, interactions between CO₂ and saline aquifers may result in potential operational limitations and have not been considered. This work utilizes reactive transport simulations to evaluate the geochemical reactions that occur during injection and extraction operational cycles for CES in a porous formation using CO₂ as a cushion gas. Simulation results are compared with similar simulations considering an injection-only flow regime of geologic CO₂ storage. Once injected, CO₂ creates conditions favorable for dissolution of carbonate and aluminosilicate minerals. However, the dissolution extent is limited in the cyclic flow regime where significantly smaller dissolution occurs after the first cycle such that CO₂ is a viable choice of cushion gas. In the injection-only flow regime, larger extents of dissolution occur as the fluid continues to be undersaturated with respect to formation minerals throughout the study period and porosity increased uniformly from 24.84% to 33.6% throughout the simulation domain. For the cyclic flow conditions, porosity increases nonuniformly to 31.1% and 25.8% closest and furthest from the injection well, respectively.

Crystallization in Confinement

Many crystallization processes of great importance, including frost heave, biomineralization, the synthesis of nanomaterials, and scale formation, occur in small volumes rather than bulk solution. Here, the influence of confinement on crystallization processes is described, drawing together information from fields as diverse as bioinspired mineralization, templating, pharmaceuticals, colloidal crystallization, and geochemistry. Experiments are principally conducted within confining systems that offer well-defined environments, varying from droplets in microfluidic devices, to cylindrical pores in filtration membranes, to nanoporous glasses and carbon nanotubes. Dramatic effects are observed, including a stabilization of metastable polymorphs, a depression of freezing points, and the formation of crystals with preferred orientations, modified morphologies, and even structures not seen in bulk. Confinement is also shown to influence crystallization processes over length scales ranging from the atomic to hundreds of micrometers, and to originate from a wide range of mechanisms. The development of an enhanced understanding of the influence of confinement on crystal nucleation and growth will not only provide superior insight into crystallization processes in many real-world environments, but will also enable this phenomenon to be used to control crystallization in applications including nanomaterial synthesis, heavy metal remediation, and the prevention of weathering.

Functionalized multiscale visual models to unravel flow and transport physics in porous structures

The fluid flow, species transport, and chemical reactions in geological formations are the chief mechanisms in engineering the exploitation of fossil fuels and geothermal energy, the geological storage of carbon dioxide (CO₂), and the disposal of hazardous materials. Porous rock is characterized by a wide surface area, where the physicochemical fluid-solid interactions dominate the multiphase flow behavior. A variety of visual models with differences in dimensions, patterns, surface properties, and fabrication techniques have been widely utilized to simulate and directly visualize such interactions in porous media. This review discusses the six categories of visual models used in geological flow applications, including packed beds, Hele-Shaw cells, synthesized microchips (also known as microfluidic chips or micromodels), geomaterial-dominated microchips, three-dimensional (3D) microchips, and nanofluidics. For each category, critical technical points (such as surface chemistry and geometry) and practical applications are summarized. Finally, we discuss opportunities and provide a framework for the development of custom-built visual models.

Evidence of Facilitated Transport in Crowded Nanopores

Fluid transport in nature often occurs through crowded nanopores, where a number of phenomena can affect it, because of fluid-fluid and fluid-solid interactions, as well as the presence of organic compounds filling the pores and their structural fluctuations. Employing molecular dynamics, we probe here the transport of fluid mixtures (CO₂-CH₄ and H₂S-CH₄) through silica nanopores filled with benzene. Both CO₂ and H₂S are strongly adsorbed within the organic-filled pore, partially displacing benzene. Unexpectedly, CO₂/H₂S adsorption facilitates CH₄ transport. Analysis of the trajectories suggests that both CO₂ and H₂S act as vehicle-like carriers and might swell benzene, generating preferential transport pathways within the crowded pore. The results are useful for identifying unexpected transport mechanisms and for developing engineering approaches that could lead to storage of CO₂ in caprocks.

Structure and dynamics of ethane confined in silica nanopores in the presence of CO2

Fundamental understanding of the subcritical/supercritical behavior of key hydrocarbon species inside nano-porous matrices at elevated pressure and temperature is less developed compared to bulk fluids, but this knowledge is of great importance for chemical and energy engineering industries. This study explores in detail the structure and dynamics of ethane (C2H6) fluid confined in silica nanopores, with a focus on the effects of pressure and different ratios of C2H6 and CO2 at non-ambient temperature. Quasi-elastic neutron scattering (QENS) experiments were carried out for the pure C2H6, C2H6:CO2 = 3:1, and 1:3 mixed fluids confined in 4-nm cylindrical silica pores at three different pressures (30 bars, 65 bars, and 100 bars) at 323 K. Two Lorentzian functions were required to fit the spectra, corresponding to fast and slow translational motions. No localized motions (rotations and vibrations) were detected. Higher pressures resulted in hindrances of the diffusivity of C2H6 molecules in all systems investigated. Pore size was found to be an important factor, i.e., the dynamics of confined C2H6 is more restricted in smaller pores compared to the larger pores used in previous studies. Molecular dynamics simulations were performed to complement the QENS experiment at 65 bars, providing supportive structure information and comparable dynamic information. The simulations indicate that CO2 molecules are more strongly attracted to the pore surface compared to C2H6. The C2H6 molecules interacting with or near the pore surface form a dense first layer (L1) close to the pore surface and a second less dense layer (L2) extending into the pore center. Both the experiments and simulations revealed the role that CO2 molecules play in enhancing C2H6 diffusion ("molecular lubrication") at high CO2:C2H6 ratios. The energy scales of the two dynamic components, fast and slow, quantified by both techniques, are in very good agreement. Herein, the simulations identified the fast component as the main contributor to the dynamics. Molecule motions in the L2 region are mostly responsible for the dynamics (fast and slow) that can be detected by the instrument.

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

Assessment of the risks and environmental impacts of carbon geosequestration requires knowledge about the wetting behavior of mineral surfaces in the presence of CO₂ and the pore fluids. In this context, the interfacial tension (IFT) between CO₂ and the aqueous fluid and the contact angle, θ, with the pore mineral surfaces are the two key parameters that control the capillary pressure in the pores of the candidate host rock. Knowledge of these two parameters and their dependence on the local conditions of pressure, temperature, and salinity is essential for the correct prediction of structural and residual trapping. We have performed classical molecular dynamics simulations to predict the CO₂-water IFT and the CO₂-water-calcite contact angle. The IFT results are consistent with previous simulations, where simple point charge water models have been shown to underestimate the water surface tension, thus affecting the simulated IFT values. When combined with the EPM2 CO₂ model, the SPC/Fw water model indeed underestimates the IFT in the low-pressure region at all temperatures studied. On the other hand, at high pressure and low temperature, the IFT is overestimated by ∼5 mN/m. Literature data regarding the CO₂/water/calcite contact angle on calcite are contradictory. Using our new set of force field parameters, we performed NVT simulations at 323 K and 20 MPa to calculate the contact angle of a water droplet on the calcite {10.4} surface in a CO₂ atmosphere. We performed simulations for both spherical and cylindrical droplet configurations for different initial radii to study the size dependence of the water contact angle on calcite in the presence of CO₂. Our results suggest that the contact angle of a cylindrical droplet, is independent of droplet size, for droplets with a radius of 50 Å or more. On the contrary, spherical droplets make a contact angle that is strongly influenced by their size. At the largest size explored in this study, both spherical and cylindrical droplets converge to the same contact angle, 38°, indicating that calcite is strongly wetted by water.

Supercritical CO2 and CH4 Uptake by Illite-Smectite Clay Minerals

Clay minerals abound in sedimentary formations and the interaction of reservoir gases with their submicron features have direct relevance to many geoenergy applications. The quantification of gas uptake over a broad range of pressures is key toward assessing the significance of these physical interactions on enhancing storage capacity and gas recovery. We report a systematic investigation of the sorption properties of three source clay minerals-Na-rich montmorillonite (SWy-2), illite-smectite mixed layer (ISCz-1), and illite (IMt-2)-using CO₂ and CH₄ up to 30 MPa at 25-115 °C. The textural characterization of the clays by gas physisorption indicates that micropores are only partly accessible to N₂ (77 K) and Ar (87 K), while larger uptakes are measured with CO₂ (273 K) in the presence of illite. The supercritical excess sorption experiments confirm these findings while revealing differences in uptake capacities that originate from the clay-specific pore size distribution. The lattice density functional theory model describes accurately the measured sorption isotherms by using a distribution of properly weighted slit pores and clay-specific solid-fluid interaction energies, which agree with isosteric heats of adsorption obtained experimentally. The model indicates that the maximum degree of pore occupancy is universal to the three clays and the two gases, and it depends solely on temperature, reaching values near unity at the critical temperature. These observations greatly support the model's predictive capability for estimating gas adsorption on clay-bearing rocks and sediments.

Anomalously low activation energy of nanoconfined MgCO3 precipitation

Experimental study of nanoconfined MgCO₃ nucleation and growth processes reveals elevated kinetics due to less strongly hydrated Mg²⁺.

Effects of Na+ and K+ Exchange in Interlayers on Biotite Dissolution under High-Temperature and High-CO2-Pressure Conditions

Cations in formation brine can affect CO₂-induced dissolution of minerals during geologic CO₂ sequestration (GCS), affecting the GCS performance. This study investigated the dissolution of biotite with 0-4 M Na⁺ and 0-10 mM K⁺ under high temperature and high CO₂ pressure (i.e., 95 °C and 100 bar CO₂). At <0.5 M Na⁺ concentration, Na⁺ replaced K⁺ in the biotite interlayer and enhanced the biotite dissolution. In >0.5 M Na⁺, however, the enhancing effect of Na⁺ was mitigated by an inhibition caused by competing sorption between Na⁺ and protons. With 0.5 M Na⁺ concentration, coexisting K⁺ significantly inhibited the biotite dissolution with high sensitivity at even lower K⁺ concentrations, such as 0.1-0.5 mM. In this study, we also reported the dissolution of Na-treated biotite, mimicking biotite naturally equilibrated with Na⁺-abundant brine. Na-treated biotite dissolved faster than natural K-containing biotite, and during the dissolution, it transformed to vermiculite. Aqueous Na⁺ inhibited the dissolution of Na-treated biotite by suppressing the release of interlayer Na⁺, and aqueous K⁺ inhibited the dissolution of Na-treated biotite by replacing the interlayer Na⁺. These findings contribute to better understanding of biotite dissolution in the presence of potassium-containing clay-swelling inhibitors and different salinities at GCS sites.

Consumer Attitudes towards Industrial CO2 Capture and Storage Products and Technologies

This paper discusses and elicits consumer attitudes towards industrial carbon dioxide (CO2) capture and storage (ICCS) products and technologies. It presents a comprehensive review of the relevant research literature on consumer attitudes towards ICCS represented by the willingness-to-pay (WTP) and willingness-to-accept (WTA) negative externalities and outcomes of the carbon capture and storage (the so-called “not-in-my-backyard” (NIMBY) approach). In addition, it employs a concise empirical model that uses the data from the online questionnaire survey conducted in 7 European Union (EU) countries with and without ICSS sites. Our results demonstrate that having at least one ICCS site significantly reduces the WTA for the ICCS products and technologies. It is shown that further increase of ICCS sites, including those in the neighboring regions and countries, leads to the increase of negative consumer attitudes to the ICCS technologies and renewable energy policies. It becomes apparent that the majority of consumers are willing to support industrial CO2 capture and storage only if it happens far away from their dwellings. The outcomes of this paper might be informative for the EU local industries and policy-makers who are planning the location of ICCS sites and optimizing the public support for their endeavors. Moreover, they might be relevant for the stakeholders dealing with the threat of climate change and the necessity for the decarbonization of the economy.

Scrutinising the Gap between the Expected and Actual Deployment of Carbon Capture and Storage—A Bibliometric Analysis

For many years, carbon capture and storage (CCS) has been discussed as a technology that may make a significant contribution to achieving major reductions in greenhouse gas emissions. At present, however, only two large-scale power plants capture a total of 2.4 Mt CO2/a. Several reasons are identified for this mismatch between expectations and realised deployment. Applying bibliographic coupling, the research front of CCS, understood to be published peer-reviewed papers, is explored to scrutinise whether the current research is sufficient to meet these problems. The analysis reveals that research is dominated by technical research (69%). Only 31% of papers address non-technical issues, particularly exploring public perception, policy, and regulation, providing a broader view on CCS implementation on the regional or national level, or using assessment frameworks. This shows that the research is advancing and attempting to meet the outlined problems, which are mainly non-technology related. In addition to strengthening this research, the proportion of papers that adopt a holistic approach may be increased in a bid to meet the challenges involved in transforming a complex energy system. It may also be useful to include a broad variety of stakeholders in research so as to provide a more resilient development of CCS deployment strategies.

Effect of CO2 Phase States and Flow Rate on Salt Precipitation in Shale Caprocks-A Microfluidic Study

Fracture networks inside the caprock for CO₂ storage reservoirs may serve as leakage pathways. Fluid flow through fractured caprocks and bypass conduits, however, can be restrained or diminished by mineral precipitations. This study investigates precipitation of salt crystals in an artificial fracture network as a function of pressure-temperature conditions and CO₂ phase states. The impact of CO₂ flow rate on salt precipitation was also studied. The primary research objective was to examine whether salt precipitation can block potential CO₂ leakage pathways. In this study, we developed a novel microfluidic high-pressure high-temperature vessel to house geomaterial micromodels. A fracture network was laser-scribed on the organic-rich shales of the Draupne Formation, the primary caprock for the Smeaheia CO₂ storage in Norway. Experimental observations demonstrated that CO₂ phase states influence the magnitude, distribution, and precipitation patterns of salt accumulations. The CO₂ phase states also affect the relationship between injection rate and extent of precipitated salts due to differences in solubility of water in CO₂ and density of different CO₂ phases. Injection of gaseous CO₂ resulted in higher salt precipitation compared to liquid and supercritical CO₂. It is shown that micrometer-sized halite crystals have the potential to partially or entirely clog fracture apertures.

Two-phase flow visualization under reservoir conditions for highly heterogeneous conglomerate rock: A core-scale study for geologic carbon storage

Geologic storage of carbon dioxide (CO₂) is considered a viable strategy for significantly reducing anthropogenic CO₂ emissions into the atmosphere; however, understanding the flow mechanisms in various geological formations is essential for safe storage using this technique. This study presents, for the first time, a two-phase (CO₂ and brine) flow visualization under reservoir conditions (10 MPa, 50 °C) for a highly heterogeneous conglomerate core obtained from a real CO₂ storage site. Rock heterogeneity and the porosity variation characteristics were evaluated using X-ray computed tomography (CT). Multiphase flow tests with an in-situ imaging technology revealed three distinct CO₂ saturation distributions (from homogeneous to non-uniform) dependent on compositional complexity. Dense discontinuity networks within clasts provided well-connected pathways for CO₂ flow, potentially helping to reduce overpressure. Two flow tests, one under capillary-dominated conditions and the other in a transition regime between the capillary and viscous limits, indicated that greater injection rates (potential causes of reservoir overpressure) could be significantly reduced without substantially altering the total stored CO₂ mass. Finally, the capillary storage capacity of the reservoir was calculated. Capacity ranged between 0.5 and 4.5%, depending on the initial CO₂ saturation.

Geo-material surface modification of microchips using layer-by-layer (LbL) assembly for subsurface energy and environmental applications

A key constraint in the application of microfluidic technology to subsurface flow and transport processes is the surface discrepancy between microchips and the actual rocks/soils. This research employs a novel layer-by-layer (LbL) assembly technology to produce rock-forming mineral coatings on microchip surfaces. The outcome of the work is a series of 'surface-mimetic micro-reservoirs (SMMR)' that represent multi-scales and multi-types of natural rocks/soils. For demonstration, the clay pores of sandstones and mudrocks are reconstructed by representatively coating montmorillonite and kaolinite in polydimethylsiloxane (PDMS) microchips in a wide range of channel sizes (width of 10-250 μm, depth of 40-100 μm) and on glass substrates. The morphological and structural properties of mineral coatings are characterized using a scanning electron microscope (SEM), optical microscope and profilometer. The coating stability is tested by dynamic flooding experiments. The surface wettability is characterized by measuring mineral oil-water contact angles. The results demonstrate the formation of nano- to micro-scale, fully-covered and stable mineral surfaces with varying wetting properties. There is an opportunity to use this work in the development of microfluidic technology-based applications for subsurface energy and environmental research.

Soft x-ray spectroscopy of high pressure liquid

We describe a new experimental technique that allows for soft x-ray spectroscopy studies (∼100-1000 eV) of high pressure liquid (∼100 bars). We achieve this through a liquid cell with a 100 nm-thick Si₃N₄ membrane window, which is sandwiched by two identical O-rings for vacuum sealing. The thin Si₃N₄ membrane allows soft x-rays to penetrate, while separating the high-pressure liquid under investigation from the vacuum required for soft x-ray transmission and detection. The burst pressure of the Si₃N₄ membrane increases with decreasing size and more specifically is inversely proportional to the side length of the square window. It also increases proportionally with the membrane thickness. Pressures > 60 bars could be achieved for 100 nm-thick square Si₃N₄ windows that are smaller than 65 μm. However, above a certain pressure, the failure of the Si wafer becomes the limiting factor. The failure pressure of the Si wafer is sensitive to the wafer thickness. Moreover, the deformation of the Si₃N₄ membrane is quantified using vertical scanning interferometry. As an example of the performance of the high-pressure liquid cell optimized for total-fluorescence detected soft x-ray absorption spectroscopy (sXAS), the sXAS spectra at the Ca L edge (∼350 eV) of a CaCl₂ aqueous solution are collected under different pressures up to 41 bars.

Clay, Water, and Salt: Controls on the Permeability of Fine-Grained Sedimentary Rocks

The ability to predict the permeability of fine-grained soils, sediments, and sedimentary rocks is a fundamental challenge in the geosciences with potentially transformative implications in subsurface hydrology. In particular, fine-grained sedimentary rocks (shale, mudstone) constitute about two-thirds of the sedimentary rock mass and play important roles in three energy technologies: petroleum geology, geologic carbon sequestration, and radioactive waste management. The problem is a challenging one that requires understanding the properties of complex natural porous media on several length scales. One inherent length scale, referred to hereafter as the mesoscale, is associated with the assemblages of large grains of quartz, feldspar, and carbonates over distances of tens of micrometers. Its importance is highlighted by the existence of a threshold in the core scale mechanical properties and regional scale energy uses of shale formations at a clay content X_clay ≈ 1/3, as predicted by an ideal packing model where a fine-grained clay matrix fills the gaps between the larger grains. A second important length scale, referred to hereafter as the nanoscale, is associated with the aggregation and swelling of clay particles (in particular, smectite clay minerals) over distances of tens of nanometers. Mesoscale phenomena that influence permeability are primarily mechanical and include, for example, the ability of contacts between large grains to prevent the compaction of the clay matrix. Nanoscale phenomena that influence permeability tend to be chemomechanical in nature, because they involve strong impacts of aqueous chemistry on clay swelling. The second length scale remains much less well characterized than the first, because of the inherent challenges associated with the study of strongly coupled nanoscale phenomena. Advanced models of the nanoscale properties of fine-grained media rely predominantly on the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory, a mean field theory of colloidal interactions that accurately predicts clay swelling in a narrow range of conditions (low salinity, low compaction, Na⁺ counterion). An important feature of clay swelling that is not predicted by these models is the coexistence, in most conditions of aqueous chemistry and dry bulk density, of two types of pores between parallel smectite particles: mesopores with a pore width of >3 nm that are controlled by long-range interactions (the osmotic swelling regime) and nanopores with a pore width <1 nm that are controlled by short-range interactions (the crystalline swelling regime). Nanogeochemical characterization and simulation techniques, including coarse-grained and all-atom molecular dynamics simulations, hold significant promise for the development of advanced constitutive relations that predict this coexistence and its dependence on aqueous chemistry.

Effect of Mineral Dissolution/Precipitation and CO2 Exsolution on CO2 transport in Geological Carbon Storage

Geological carbon sequestration (GCS) in deep saline aquifers is an effective means for storing carbon dioxide to address global climate change. As the time after injection increases, the safety of storage increases as the CO₂ transforms from a separate phase to CO₂(aq) and HCO₃^- by dissolution and then to carbonates by mineral dissolution. However, subsequent depressurization could lead to dissolved CO₂(aq) escaping from the formation water and creating a new separate phase which may reduce the GCS system safety. The mineral dissolution and the CO₂ exsolution and mineral precipitation during depressurization change the morphology, porosity, and permeability of the porous rock medium, which then affects the two-phase flow of the CO₂ and formation water. A better understanding of these effects on the CO₂-water two-phase flow will improve predictions of the long-term CO₂ storage reliability, especially the impact of depressurization on the long-term stability. In this Account, we summarize our recent work on the effect of CO₂ exsolution and mineral dissolution/precipitation on CO₂ transport in GCS reservoirs. We place emphasis on understanding the behavior and transformation of the carbon components in the reservoir, including CO₂(sc/g), CO₂(aq), HCO₃^-, and carbonate minerals (calcite and dolomite), highlight their transport and mobility by coupled geochemical and two-phase flow processes, and consider the implications of these transport mechanisms on estimates of the long-term safety of GCS. We describe experimental and numerical pore- and core-scale methods used in our lab in conjunction with industrial and international partners to investigate these effects. Experimental results show how mineral dissolution affects permeability, capillary pressure, and relative permeability, which are important phenomena affecting the input parameters for reservoir flow modeling. The porosity and the absolute permeability increase when CO₂ dissolved water is continuously injected through the core. The MRI results indicate dissolution of the carbonates during the experiments since the porosity has been increased after the core-flooding experiments. The mineral dissolution changes the pore structure by enlarging the throat diameters and decreasing the pore specific surface areas, resulting in lower CO₂/water capillary pressures and changes in the relative permeability. When the reservoir pressure decreases, the CO₂ exsolution occurs due to the reduction of solubility. The CO₂ bubbles preferentially grow toward the larger pores instead of toward the throats or the finer pores during the depressurization. After exsolution, the exsolved CO₂ phase shows low mobility due to the highly dispersed pore-scale morphology, and the well dispersed small bubbles tend to merge without interface contact driven by the Ostwald ripening mechanism. During depressurization, the dissolved carbonate could also precipitate as a result of increasing pH. There is increasing formation water flow resistance and low mobility of the CO₂ in the presence of CO₂ exsolution and carbonate precipitation. These effects produce a self-sealing mechanism that may reduce unfavorable CO₂ migration even in the presence of sudden reservoir depressurization.

Experimental Investigation on the Behavior of Supercritical CO2 during Reservoir Depressurization

CO₂ sequestration in saline aquifers is a promising way to address climate change. However, the pressure of the sequestration reservoir may decrease in practice, which induces CO₂ exsolution and expansion in the reservoir. In this study, we conducted a core-scale experimental investigation on the depressurization of CO₂-containing sandstone using NMR equipment. Three different series of experiments were designed to investigate the influence of the depressurization rate and the initial CO2 states on the dynamics of different trapping mechanisms. The pressure range of the depressurization was from 10.5 to 4.0 MPa, which covered the supercritical and gaseous states of the CO₂ (named as CO₂(sc) and CO₂(g), respectively). It was found that when the aqueous phase saturated initially, the exsolution behavior strongly depended on the depressurization rate. When the CO₂ and aqueous phase coexisting initially, the expansion of the CO₂(sc/g) contributed to the incremental CO₂ saturation in the core only when the CO₂ occurred as residually trapped. It indicates that the reservoir depressurization has the possibility to convert the solubility trapping to the residual trapping phase, and/or convert the residual trapping to mobile CO₂.

Modeling CO2-Water-Mineral Wettability and Mineralization for Carbon Geosequestration

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₂-H₂O 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.

Dissolution Trapping of Carbon Dioxide in Heterogeneous Aquifers

The geologic architecture in sedimentary reservoirs affects the behavior of density-driven flow and the dispersion of CO₂-rich brine. The spatial organization and connectivity of facies types play an important role. Low-permeability facies may suppress fingering and reduce vertical spreading, but may also increase transverse mixing. This is more pronounced when geologic structures create preferential flow pathways through connected facies types. We perform high-resolution simulations of three-dimensional (3D) heterogeneous formations whose connectivity cannot be represented in two-dimensional models consistent with percolation theory. This work focuses on the importance of 3D facies-based heterogeneity and connectivity on advection-diffusion transport of dissolved CO₂. Because the dissolution of CO₂ and the subsequent density increase of brine are the driving force for gravitational instabilities, we model the phase behavior with the accurate cubic-plus-association equation-of-state, which accounts for the self-association of polar water molecules and the cross-association between CO₂ and water. Our results elucidate how the spatial organization of facies affects the dynamics of CO₂ convective mixing. Scaling relations for the evolution of a global dispersion-width provide insights that can be universally applied. The results suggest that the long-term evolution and scaling of dispersion are surprisingly similar for homogeneous and (binary and multiscale) heterogeneous porous media.

Molecular Modeling of Thermodynamic and Transport Properties for CO2 and Aqueous Brines

Molecular simulation techniques using classical force-fields occupy the space between ab initio quantum mechanical methods and phenomenological correlations. In particular, Monte Carlo and molecular dynamics algorithms can be used to provide quantitative predictions of thermodynamic and transport properties of fluids relevant for geologic carbon sequestration at conditions for which experimental data are uncertain or not available. These methods can cover time and length scales far exceeding those of quantum chemical methods, while maintaining transferability and predictive power lacking from phenomenological correlations. The accuracy of predictions depends sensitively on the quality of the molecular models used. Many existing fixed-point-charge models for water and aqueous mixtures fail to represent accurately these fluid properties, especially when descriptions covering broad ranges of thermodynamic conditions are needed. Recent work on development of accurate models for water, CO₂, and dissolved salts, as well as their mixtures, is summarized in this Account. Polarizable models that can respond to the different dielectric environments in aqueous versus nonaqueous phases are necessary for predictions of properties over extended ranges of temperatures and pressures. Phase compositions and densities, activity coefficients of the dissolved salts, interfacial tensions, viscosities and diffusivities can be obtained in near-quantitative agreement to available experimental data, using relatively modest computational resources. In some cases, for example, for the composition of the CO₂-rich phase in coexistence with an aqueous phase, recent results from molecular simulations have helped discriminate among conflicting experimental data sets. The sensitivity of properties on the quality of the intermolecular interaction model varies significantly. Properties such as the phase compositions or electrolyte activity coefficients are much more sensitive than phase densities, viscosities, or component diffusivities. Strong confinement effects on physical properties in nanoscale media can also be directly obtained from molecular simulations. Future work on molecular modeling for CO₂ and aqueous brines is likely to be focused on more systematic generation of interaction models by utilizing quantum chemical as well as direct experimental measurements. New ion models need to be developed for use with the current generation of polarizable water models, including ion-ion interactions that will allow for accurate description of dense, mixed brines. Methods will need to be devised that go beyond the use of effective potentials for incorporation of quantum effects known to be important for water, and reactive force fields developed that can handle bond creation and breaking in systems with carbonate and silicate minerals. Another area of potential future work is the integration of molecular simulation methods in multiscale models for the chemical reactions leading to mineral dissolution and flow within the porous media in underground formations.

Effect of Layer Charge on CO2 and H2O Intercalations in Swelling Clays

The effect of layer charge on the intercalation of supercritical carbon dioxide (scCO₂)-H₂O mixture in Na-montmorillonite clay interlayers under T = 323 K and P = 90 bar geologic sequestration conditions has been further investigated. This effect includes the charge amount and its location (within either octahedral or tetrahedral layers due to isomorphic substitutions). Two clay models with different layer charges are used in this study. Simulation results show that the increase of charge amount shifts the monolayer-to-bilayer (1W-to-2W) hydration transition toward the lower relative humidity (RH), increasing water sorption at the expense of reducing the overall sorption amount of CO₂ in the clay interlayer. However, the combination of the influence of charge amount and charge location leads to insignificant changes in equilibrium basal spacings of the high- and low-charge clays. Molecular dynamics simulations show that the CO₂ dimers, which are frequently seen in low-charge clay interlayers, vanish in high-charge clay interlayers even at low RH of 30%.

Anorthite Dissolution under Conditions Relevant to Subsurface CO2 Injection: Effects of Na+, Ca2+, and Al3

Supercritical CO₂ is injected into subsurface environments during geologic CO₂ sequestration and CO₂-enhanced oil recovery. In these processes, the CO₂-induced dissolution of formation rocks, which contain plagioclase, can affect the safety and efficiency of the subsurface operation. In subsurface brines, Na⁺ and Ca²⁺ are naturally abundant, and Al³⁺ concentration increases due to acidification by injected CO₂. However, our current understanding of cation effects on plagioclase dissolution does not provide sufficiently accurate prediction of plagioclase dissolution at such high salinities. This study investigated the effects of up to 4 M Na⁺, 1 M Ca²⁺, and 200 μM Al³⁺ on anorthite (as a representative mineral of Ca-containing plagioclase) dissolution under conditions closely relevant to subsurface CO₂ injection. For the first time, we elucidated the inhibition effects of Al³⁺ on anorthite dissolution in far-from-equilibrium systems, and found that the Al³⁺ effects were enhanced at elevated temperature. Interestingly, Na⁺ inhibited anorthite dissolution as well, and the effects of Na⁺ were 50% stronger at 35 °C than at 60 °C. Ca²⁺ had similar effects to those of Na⁺, and the Ca²⁺ effects did not suppress Na⁺ effects when they coexisted. These findings can contribute to better predicting plagioclase dissolution in geologic formations and will also be helpful in improving designs for subsurface CO₂ injection.

The dynamic nature of crystal growth in pores

The kinetics of crystal growth in porous media controls a variety of natural processes such as ore genesis and crystallization induced fracturing that can trigger earthquakes and weathering, as well as, sequestration of CO₂ and toxic metals into geological formations. Progress on understanding those processes has been limited by experimental difficulties of dynamically studying the reactive surface area and permeability during pore occlusion. Here, we show that these variables cause a time-dependency of barite growth rates in microporous silica. The rate is approximately constant and similar to that observed on free surfaces if fast flow velocities predominate and if the time-dependent reactive surface area is accounted for. As the narrower flow paths clog, local flow velocities decrease, which causes the progressive slowing of growth rates. We conclude that mineral growth in a microporous media can be estimated based on free surface studies when a) the growth rate is normalized to the time-dependent surface area of the growing crystals, and b) the local flow velocities are above the limit at which growth is transport-limited. Accounting for the dynamic relation between microstructure, flow velocity and growth rate is shown to be crucial towards understanding and predicting precipitation in porous rocks.