Water and methane in shale rocks: Flow pattern effects on fluid transport and pore structure

Water Bridges in Clay Nanopores: Mechanisms of Formation and Impact on Hydrocarbon Transport

Clays are prevalent in the earth's crust and usually deposited in the presence of water. An unusual finding in clays is that under certain conditions, water molecules can collectively form a bridge across a clay-hosted pore. However, there are relatively few studies focused on the formation mechanism of the water bridge in clay nanopores. In this work, we use molecular dynamics simulations to investigate the formation of the water bridge and its influence on fluid transport in slit-shaped illite nanopores. Two different basal illite surface chemistries are constructed: potassium-hydroxyl (P-H) and hydroxyl-hydroxyl (H-H) structures. Because pore size and water concentration are expected to control the formation of the water bridge, our simulations span a wide range of pore sizes and water concentrations. Generally, positive potassium layers and negative hydroxyl groups in P-H nanopore can induce partial charges which in return produce instant and local electric fields, favoring the formation of the water bridge. In P-H nanopores, the water bridge happens at a relatively low water concentration. However, in H-H nanopores, the water bridge only forms at high water concentrations. Additionally, smaller pore sizes favor the formation of water bridges. However, the presence of an electric field promotes the formation of a water bridge even in larger pore sizes in P-H pores. The results also indicate that in both P-H and H-H nanopores, water adsorption films initially create a smooth surface to promote the hydrocarbon flow. In P-H nanopores, further increases in the water concentration causes a sharp decline in the self-diffusion coefficients of the hydrocarbon and water due to the formation of the water bridge. The presence of electric fields in P-H pores can however weaken the confinement effect of illite and promote the hydrocarbon flow. In contrast, in H-H nanopores, the self-diffusion coefficients decline slowly with the increase of water concentration. This is because no water bridge is formed at low water concentrations in H-H nanopores.

Quantifying Pore Width Effects on Diffusivity via a Novel 3D Stochastic Approach with Input from Atomistic Molecular Dynamics Simulations

The increased production of unconventional hydrocarbons emphasizes the need to understand the transport of fluids through narrow pores. Although it is well-known that confinement affects fluids structure and transport, it is not yet possible to quantitatively predict properties such as diffusivity as a function of pore width in the range of 1-50 nm. Such pores are commonly found not only in shale rocks but also in a wide range of engineering materials, including catalysts. We propose here a novel and computationally efficient methodology to obtain accurate diffusion coefficient predictions as a function of pore width for pores carved out of common materials, such as silica, alumina, magnesium oxide, calcite, and muscovite. We implement atomistic molecular dynamics (MD) simulations to quantify fluid structure and transport within 5 nm-wide pores, with particular focus on the diffusion coefficient within different pore regions. We then use these data as input to a bespoke stochastic kinetic Monte Carlo (KMC) model, developed to predict fluid transport in mesopores. The KMC model is used to extrapolate the fluid diffusivity for pores of increasing width. We validate the approach against atomistic MD simulation results obtained for wider pores. When applied to supercritical methane in slit-shaped pores, our methodology yields data within 10% of the atomistic simulation results, with significant savings in computational time. The proposed methodology, which combines the advantages of MD and KMC simulations, is used to generate a digital library for the diffusivity of gases as a function of pore chemistry and pore width and could be relevant for a number of applications, from the prediction of hydrocarbon transport in shale rocks to the optimization of catalysts, when surface-fluid interactions impact transport.

Enhancement of oil flow in shale nanopores by manipulating friction and viscosity

Understanding the viscosity and friction of a fluid under nanoconfinement is the key to nanofluidics research. Existing work on nanochannel flow enhancement has been focused on simple systems with only one to two fluids considered such as water flow in carbon nanotubes, and large slip lengths have been found to be the main factor for the massive flow enhancement. In this study, we use molecular dynamics simulations to study the fluid flow of a ternary mixture of octane-carbon dioxide-water confined within two muscovite and kerogen surfaces. The results indicate that, in a muscovite slit, supercritical CO₂ (scCO₂) and H₂O both enhance the flow of octane due to (i) a decrease in the friction of octane with the muscovite wall because of the formation of thin layers of H₂O and scCO₂ near the surfaces; and (ii) a reduction in the viscosity of octane in nanoconfinement. Water reduces octane viscosity by weakening the interaction of octane with the muscovite surface, while scCO₂ reduces octane viscosity by weakening both octane-octane and octane-surface interactions. In a kerogen slit, water does not play any significant role in changing the friction or viscosity of octane. In contrast, scCO₂ reduces both the friction and the viscosity of octane, and the enhancement of octane flow is mainly caused by the reduction of viscosity. Our results highlight the importance of multicomponent interactions in nanoscale fluid transport. The results presented here also have a direct implication in enhanced oil recovery in unconventional reservoirs.

Effects of water on the stochastic motions of propane confined in MCM-41-S pores

Quasi-elastic neutron scattering (QENS) and molecular dynamics simulations (MDS) reveal the effects of water on the structure and dynamics of propane confined in 1.5 nm wide cylindrical pores of MCM-41-S.

Invasion of gas into mica nanopores: a molecular dynamics study

The invasion of gas into liquid-filled nanopores is encountered in many engineering problems but is not yet well understood. We report molecular dynamics simulations of the invasion of methane gas into water-filled mica pores with widths of 2-6 nm. Gas invades into a pore only when the pressure exceeds a breakthrough pressure and a thin residual water film is left on the mica wall as the gas phase moves deeper into the pore. The gas breakthrough pressure of pores as narrow as 2 nm can be modeled reasonably well by the capillary pressure if the finite thickness of residual liquid water film and the liquid-gas interface are taken into account. The movement of the front of the liquid meniscus during gas invasion can be quantitatively described using the classical hydrodynamics when the negative slip length on the strongly hydrophilic mica walls is taken into account. Understanding the molecular mechanisms underlying the gas invasion in the system studied here will form the foundation for quantitative prediction of gas invasion in practical porous media.

Enhanced Ion Adsorption on Mineral Nanoparticles

Classical molecular dynamics simulation was used to study the adsorption of Na⁺, Ca²⁺, Ba²⁺, and Cl^- ions on gibbsite edge (1 0 0), basal (0 0 1), and nanoparticle (NP) surfaces. The gibbsite NP consists of both basal and edge surfaces. Simulation results indicate that Na⁺ and Cl^- ions adsorb on both (1 0 0) and (0 0 1) surfaces as inner-sphere species (i.e., no water molecules between an ion and the surface). Outer-sphere Cl^- ions (i.e., one water molecule between an ion and the surface) were also found on these surfaces. On the (1 0 0) edge, Ca²⁺ ions adsorb as inner-sphere and outer-sphere complexes, whereas on the (0 0 1) surface, outer-sphere Ca²⁺ ions are the dominant species. Ba²⁺ ions were found as inner-sphere and outer-sphere complexes on both surfaces. Calculated ion surface coverages indicate that, for all ions, surface coverages are always higher on the basal surface compared to those on the edge surface. More importantly, surface coverages for cations on the gibbsite NP are always higher than those calculated for the (1 0 0) and (0 0 1) surfaces. This enhanced ion adsorption behavior for the NP is due to the significant number of inner-sphere cations found at NP corners. Outer-sphere cations do not contribute to the enhanced surface coverage. In addition, there is no ion adsorption enhancement observed for the Cl^- ion. Our work provides a molecular-scale understanding of the relative significance of ion adsorption onto gibbsite basal versus edge surfaces and demonstrates the corner effect on ion adsorption on NPs.

Atomistic Structure of Mineral Nano-aggregates from Simulated Compaction and Dewatering

The porosity of clay aggregates is an important property governing chemical reactions and fluid flow in low-permeability geologic formations and clay-based engineered barrier systems. Pore spaces in clays include interlayer and interparticle pores. Under compaction and dewatering, the size and geometry of such pore spaces may vary significantly (sub-nanometer to microns) depending on ambient physical and chemical conditions. Here we report a molecular dynamics simulation method to construct a complex and realistic clay-like nanoparticle aggregate with interparticle pores and grain boundaries. The model structure is then used to investigate the effect of dewatering and water content on micro-porosity of the aggregates. The results suggest that slow dewatering would create more compact aggregates compared to fast dewatering. Furthermore, the amount of water present in the aggregates strongly affects the particle-particle interactions and hence the aggregate structure. Detailed analyses of particle-particle and water-particle interactions provide a molecular-scale view of porosity and texture development of the aggregates. The simulation method developed here may also aid in modeling the synthesis of nanostructured materials through self-assembly of nanoparticles.

Quantifying the anisotropy and tortuosity of permeable pathways in clay-rich mudstones using models based on X-ray tomography

The permeability of shales is important, because it controls where oil and gas resources can migrate to and where in the Earth hydrocarbons are ultimately stored. Shales have a well-known anisotropic directional permeability that is inherited from the depositional layering of sedimentary laminations, where the highest permeability is measured parallel to laminations and the lowest permeability is perpendicular to laminations. We combine state of the art laboratory permeability experiments with high-resolution X-ray computed tomography and for the first time can quantify the three-dimensional interconnected pathways through a rock that define the anisotropic behaviour of shales. Experiments record a physical anisotropy in permeability of one to two orders of magnitude. Two- and three-dimensional analyses of micro- and nano-scale X-ray computed tomography illuminate the interconnected pathways through the porous/permeable phases in shales. The tortuosity factor quantifies the apparent decrease in diffusive transport resulting from convolutions of the flow paths through porous media and predicts that the directional anisotropy is fundamentally controlled by the bulk rock mineral geometry. Understanding the mineral-scale control on permeability will allow for better estimations of the extent of recoverable reserves in shale gas plays globally.

Propane-Water Mixtures Confined within Cylindrical Silica Nanopores: Structural and Dynamical Properties Probed by Molecular Dynamics

Despite the multiple length and time scales over which fluid-mineral interactions occur, interfacial phenomena control the exchange of matter and impact the nature of multiphase flow, as well as the reactivity of C-O-H fluids in geologic systems. In general, the properties of confined fluids, and their influence on porous geologic phenomena are much less well understood compared to those of bulk fluids. We used equilibrium molecular dynamics simulations to study fluid systems composed of propane and water, at different compositions, confined within cylindrical pores of diameter ∼16 Å carved out of amorphous silica. The simulations are conducted within a single cylindrical pore. In the simulated system all the dangling silicon and oxygen atoms were saturated with hydroxyl groups and hydrogen atoms, respectively, yielding a total surface density of 3.8 -OH/nm2. Simulations were performed at 300 K, at different bulk propane pressures, and varying the composition of the system. The structure of the confined fluids was quantified in terms of the molecular distribution of the various molecules within the pore as well as their orientation. This allowed us to quantify the hydrogen bond network and to observe the segregation of propane near the pore center. Transport properties were quantified in terms of the mean square displacement in the direction parallel to the pore axis, which allows us to extract self-diffusion coefficients. The diffusivity of propane in the cylindrical pore was found to depend on pressure, as well as on the amount of water present. It was found that the propane self-diffusion coefficient decreases with increasing water loading because of the formation of water bridges across the silica pores, at sufficiently high water content, which hinder propane transport. The rotational diffusion, the lifespan of hydrogen bonds, and the residence time of water molecules at contact with the silica substrate were quantified from the simulated trajectories using the appropriate autocorrelation functions. The simulations contribute to a better understanding of the molecular phenomena relevant to the behavior of fluids in the subsurface.

Confinement Effects on Carbon Dioxide Methanation: A Novel Mechanism for Abiotic Methane Formation

An important scientific debate focuses on the possibility of abiotic synthesis of hydrocarbons during oceanic crust-seawater interactions. While on-site measurements near hydrothermal vents support this possibility, laboratory studies have provided data that are in some cases contradictory. At conditions relevant for sub-surface environments it has been shown that classic thermodynamics favour the production of CO2 from CH4, while abiotic methane synthesis would require the opposite. However, confinement effects are known to alter reaction equilibria. This report shows that indeed thermodynamic equilibrium can be shifted towards methane production, suggesting that thermal hydrocarbon synthesis near hydrothermal vents and deeper in the magma-hydrothermal system is possible. We report reactive ensemble Monte Carlo simulations for the CO2 methanation reaction. We compare the predicted equilibrium composition in the bulk gaseous phase to that expected in the presence of confinement. In the bulk phase we obtain excellent agreement with classic thermodynamic expectations. When the reactants can exchange between bulk and a confined phase our results show strong dependency of the reaction equilibrium conversions, [Formula: see text], on nanopore size, nanopore chemistry, and nanopore morphology. Some physical conditions that could shift significantly the equilibrium composition of the reactive system with respect to bulk observations are discussed.

Surface hydration drives rapid water imbibition into strongly hydrophilic nanopores

Surface hydration-driven imbibition of water into strongly hydrophilic pores follows a diffusive scaling law and exhibits effective diffusion coefficients much higher than water molecules.

Modelling aqueous solubility of sodium chloride in clays at thermodynamic conditions of hydraulic fracturing by molecular simulations

To address the high salinity of flow-back water during hydraulic fracturing, we have studied the equilibrium partitioning of NaCl and water between the bulk phase and clay pores.