Ion-specific interactions at calcite-brine interfaces: a nano-scale study of the surface charge development and preferential binding of polar hydrocarbons

Atomistic Insights into Ion-Driven Interactions of Calcite/Carbonated Brine/Polar Model Oil: Implications for Carbonated Smart Waterflooding

This study investigates the fundamental ion-specific (Na+, Cl-, Mg2+, and SO42-) interactions governing a polar model oil (decane + benzoic acid) at the calcite/carbonated brine interface by adopting a fully atomistic molecular dynamics (MD) simulation. By bridging molecular-scale interactions with macroscopic mechanisms, such as interfacial tension (IFT) reduction, oil viscosity, and wettability changes, this work provides the first direct mechanistic validation of phenomena that have previously been inferred only from experimental observations in carbonated smart water flooding systems. The results demonstrate that enhanced interactions between carboxylic acids and anions at the oil/brine interface significantly influence CO2 diffusion and distribution within the oleic phase, which affects the apparent oil viscosity. While variations in brine ionic composition cause only modest changes in IFT, a pronounced reduction is observed with increased concentrations of polar molecules in the oil phase. Structural analysis reveals that divalent ions (Mg2+, SO42-) are excluded from the hydration layers near the calcite surface but alter the arrangement of Na+ and Cl- ions in the hydration layer covering the calcite surface, thereby influencing wettability. Notably, SO42- neutralizes the calcite surface positive charge and facilitates Mg2+ access to the interface, promoting desorption of benzoic acid (BA) from the surface through the Mg-BA association. This highlights the cooperative role of SO42- and Mg2+ in releasing polar species from the calcite surface. The findings underscore the dominant influence of IFT over contact angle in capillary-driven recovery and show that apparent viscosity is more sensitive to CO2 content and overall salinity than specific ions. Therefore, from an industrial perspective, maintaining seawater-like salinity enriched with divalent ions offers a practical strategy to enhance the mobilization of polar acidic components during carbonated water flooding in carbonate reservoirs, supporting the design of more efficient Enhanced Oil Recovery (EOR) formulations.

Effect of Ion-Specific Hydration Forces on the Stability of Water Films on Calcite Surfaces

The hydration force is indispensable for understanding short-range interfacial forces in aqueous systems. Perturbation of the hydration structure by ions generates an ion-specific hydration force. Surface-force measurements on calcite surfaces have suggested that Na+ decreases the repulsive hydration force by directly adsorbing the surface and disrupting the hydration layers. However, the influence of structural changes on the surface force remains unclear. We conducted molecular dynamics simulations for water films between calcite (104) surfaces and oil/water interfaces. Ion-specific hydration forces estimated by the simulations were consistent with the experimental results. Notably, the ion-specific hydration forces cannot be explained solely by the structure of water molecules because ions do not significantly change the structure of the hydration layers, such as density distributions and orientations. We propose a novel mechanism whereby ion-specific electrostatic potentials in the water films control the adhesive and repulsive nature of the interfaces. The directly adsorbed Na+ on the calcite causes the monotonically decreasing electrostatic potential from the calcite surface, thereby enhancing adhesion. Ca2+ results in a convex shape of the electrostatic potential curve, which enhances repulsion. Importantly, the shape of the electrostatic potential curve depends on the Stern layer structure and the perturbation between the surface and interfaces. This study offers important insight for interpreting surface-force measurements in aqueous systems.

Microscopic insights into the effects of interfacial dynamics and nanoconfinement on characteristics of calcium carbonate clusters within two-dimensional nanochannels

Herein, the interfacial effects on calcium carbonate clustering within two-dimensional (2D) graphene nanochannels were systematically investigated using molecular dynamics simulations. The distribution characteristics of the ions at the interface can be attributed to the ordered water layers within the 2D nanochannels. The orientation of CO32- is approximately perpendicular to the interface, which can be attributed to hydrogen bonding and its association with Ca2+ at the interface region. The results show that characteristics of CaCO3 clusters can be affected by ion dynamics at the interface and nanoconfinement, although they prefer to locate in the bulk-like region. Due to nanoconfinement, ion dynamics are slowed down, especially in the direction perpendicular to the graphene surface. Due to the distribution and orientation characteristics of CO32- in the interface region, particularly considering the hydration dynamics of Ca2+ and CO32-, the association between Ca2+ and CO32- ions in CaCO3 clusters at the interface can be promoted as Ca2+ moves from the interface region to the bulk-like region. The ion dynamics and coordination characteristics of CaCO3 near the interface region within 2D nanochannels facilitate the formation of CaCO3 clusters with highly coordinated Ca2+-CO32- structures, which might favor the nucleation of aragonite. The results provide insight into the effects of nanoconfinement and interfacial water layers on biomineral nucleation and offer theoretical insights into the new preparation methods of novel inorganic functional materials.

Experimental study of the effect of oil polarity on smart waterflooding in carbonate reservoirs

This study investigates the influence of oil polarity on interfacial tension (IFT), contact angle, oil recovery, and effluent pH in smart water and Low-salinity water injection. The results indicate that the interaction between the hydration shell of ions and the polar components (PCs) of oil is crucial. Increasing oil polarity enhances the potential for interaction with the hydration shell of ions, leading to reduced IFT, altered wettability, and improved oil recovery; which could be boosted by the contribution of a higher number of anions in the smart water bulk through the enhancement of their interaction with the PCs (especially acidic components) of oil. The study demonstrates that increasing the SO42- concentration in seawater increased oil recovery for oils with higher acid component content, as indicated by total acid number values of 0.87, 0.99, and 1.32 mgKOH/g, the tertiary oil recovery factors for these oils were 61.10%, 69.82%, and 87.09%, respectively. The effluent pH results align with the findings of contact angle and oil recovery, confirming the dominant influence of anions on oil recovery. The interaction between the PCs of oil and the hydration shell of ions is thus highlighted as a critical factor in the observed outcomes.

Salinity-Driven Structural and Viscosity Modulation of Confined Polar Oil Phases by Carbonated Brine Films: Novel Insights from Molecular Dynamics

The structural and dynamic properties of fluids under confinement in a porous medium differ from their bulk properties. This study delves into the surface structuring and hydrodynamic characteristics of oil/thin film carbonated brine two-phase within a calcite channel upon salinity variation. To this end, both equilibrium and non-equilibrium molecular dynamics simulations are utilized to unveil the effect of the carboxylic acid component (benzoic acid) in a simple model oil (decane) confined between two thin films of carbonated brine on the oil-brine-calcite characteristics. The salinity effect was scrutinized under four saline carbonated waters, deionized carbonated water (DCW), carbonated low-salinity brine (CLSB, 30,000 ppm), carbonated seawater (CSW, 60,000 ppm), and carbonated high-salinity brine (CHSB, 180,000 ppm). An electrical double layer (EDL) is observed at varying salinities, comprising a Stern-like positive layer (formed by Na+ ions) followed by a negative one (formed by Cl- ions primarily residing on top of the adsorbed sodium cations). By lowering the salinity, the Na+ ions cover the interface regions (brine-calcite and brine-oil), depleting within the brine bulk region. The lowest positive surface charge on the rock surface was found in salinity corresponding to seawater. Two distinct Na+ peaks at the oleic phase interface have been observed in the carbonated high-salinity brine system, enhancing the adsorption of polar molecules at the thin brine film interfaces. There is a pronounced EDL formation at the oleic phase interface in the case of CSW, resulting in a strong interface region containing ions and functional fractions. Likewise, the oil region confined by CSW exhibited the lowest apparent viscosity, attributed to the optimized salinity distribution and inclination of benzoic acid fractions uniformly at the brine-oil interface, acting as a slippery surface. Moreover, the results reveal that the presence of polar fractions could increase the oil phase's apparent viscosity, and introducing ions to this system reduces the polar molecules' destructive effect on the apparent viscosity of the oil region. Therefore, the fluidity of confined systems is modulated by both composition of the brine and oil phases.

Ordered and Disordered Carboxylic Acid Monolayers on Calcite (104) and Muscovite (001) Surfaces

The adsorption of carboxylic acid molecules at the calcite (104) and the muscovite (001) surface was investigated using surface X-ray diffraction. All four investigated carboxylic acid molecules, hexanoic acid, octanoic acid, lauric acid, and stearic acid, were found to adsorb at the calcite surface. Whereas the shortest two carboxylic acid molecules, hexanoic acid and octanoic acid, showed limited ordering and a flexible, disordered chain, the two longest carboxylic acid molecules form fully ordered monolayers, i.e., these form highly structured self-assembled monolayers. The latter molecules are oriented almost fully upright, with a tilt of up to 10°. The oxygen atoms of the organic molecules are found at similar positions to those of water molecules at the calcite-water interface. This suggests that in both cases, the oxygen atoms compensate for the broken bonds at the calcite surface. Under the same experimental conditions, stearic acid does not adsorb to K⁺ and Ca²⁺-functionalized muscovite mica because the neutral molecules do not engage in the ionic bonds typical for the mica interface. These differences in adsorption behavior are characteristic for the differences of the oil-solid interactions in carbonate and sandstone reservoirs.

Study of the Adsorption Behavior of Surfactants on Carbonate Surface by Experiment and Molecular Dynamics Simulation

Surfactants adsorption onto carbonate reservoirs would cause surfactants concentrations decrease in surfactant flooding, which would decrease surfactant efficiency in practical applications of enhanced oil recovery (EOR) processes. Different surfactants could be classified as cationic surfactants, anionic surfactants, non-ionic surfactants according to the main charge, or be classified as chemical surfactant and bio-surfactant according to the surfactant origin. However, the research on different type surfactants adsorption on carbonate reservoirs surface differences was few. Therefore, five representative surfactants (CTAB, SDS, TX-100, sophorolipid, rhamonilipid) adsorption effect onto carbonate reservoirs surface was studied. Owing to the fact that the salinity and temperature in underground carbonate reservoirs were high during the EOR process, it is vital to study the salinity effect and temperature effect on surfactant adsorption. In this study, different surfactants species, temperature and salinity adsorption onto carbonate reservoirs were studied. The adsorption isotherms were fitted by Langmuir, Freundlich, Temkin and Linear models, and the first three models fitting effect were good. The results showed that cationic surfactants adsorption quantity was higher than anionic surfactants, and the non-ionic surfactants adsorption quantity was the lowest. When the temperature increased, the surfactants adsorption would decrease, because the adsorption process was exothermic process, and increasing temperature would inhibit the adsorption. The higher salinity would increase surfactants adsorption because higher salinity could compress electric double layer. In order to decrease surfactants adsorption, SiO2 nanoparticles and TiO2 nanoparticles were added to surfactants solutions, and then surfactants could adsorb onto nanoparticles surface, then the steric hindrance between surfactant molecules would increase, which could decrease surfactants adsorption. Contact angle results indicated that surfactants adsorption made the carbonate reservoir wettability alteration. In the end, surfactants (with or without SiO2 nanoparticles) adsorption onto carbonate reservoirs mechanism were studied by molecular dynamics simulation. The simulation results indicated that the surfactants molecules could adsorb onto SiO2 nanoparticles surface, and then the surfactants adsorption quantity onto carbonate rocks would decrease, which was in accordance with the experiments results.

A systematic and critical review of application of molecular dynamics simulation in low salinity water injection

Low Salinity Water Injection (LSWI) has been a well-researched EOR method, with several experimental and theoretical scientific papers reported in the literature over the past few decades. Despite this, there is still an ongoing debate on dominant mechanisms behind this complex EOR process, and some issues remain elusive. Part of the complexity arises from the scale of investigation, which spans from sub-pore scale (atomic and electronic scale) to pore scale, core scale, and reservoir scale. Molecular Dynamics (MD) simulation has been used as a research tool in the past decade to investigate the nano-scale interactions among reservoir rock (e.g., calcite, silica), crude oil, and brine systems in presence of some impurities (e.g., clay minerals) and additives (e.g., nanoparticles). In this paper, fundamental concepts of MD simulation and common analyses driven by MD are briefly reviewed. Then, an overview of molecular models of the most common minerals encountered in petroleum reservoirs: quartz, calcite, and clay, with their most common types of potential function, is provided. Next, a critical review and in depth analysis of application of MD simulations in LSWI process in both sandstone and carbonate reservoirs in terms of sub-pore scale mechanisms, namely electrical double layer (EDL) expansion, multi-ion exchange (MIE), and cation hydration, is presented to scrutinize role of salinity, ionic composition, and rock surface chemistry from an atomic level. Some inconsistencies observed in the literature are also highlighted and the reasons behind them are explained. Finally, a future research guide is provided after critically discussing the challenges and potential of the MD in LSWI to shed more light on governing mechanisms behind LSWI by enhancing the reliability of MD outcomes in future researches. Such insights can be used for design of new MD researches with complementary experimental studies at core scale to capture the main mechanisms behind LSWI.

Atomistic insight into salinity dependent preferential binding of polar aromatics to calcite/brine interface: implications to low salinity waterflooding

This paper resolve the salinity-dependent interactions of polar components of crude oil at calcite-brine interface in atomic resolution. Molecular dynamics simulations carried out on the present study showed that ordered water monolayers develop immediate to a calcite substrate in contact with a saline solution. Carboxylic compounds, herein represented by benzoic acid (BA), penetrate into those hydration layers and directly linking to the calcite surface. Through a mechanism termed screening effect, development of hydrogen bonding between –COOH functional groups of BA and carbonate groups is inhibited by formation of a positively-charged Na⁺ layer over CaCO₃ surface. Contrary to the common perception, a sodium-depleted solution potentially intensifies surface adsorption of polar hydrocarbons onto carbonate substrates; thus, shifting wetting characteristic to hydrophobic condition. In the context of enhanced oil recovery, an ion-engineered waterflooding would be more effective than injecting a solely diluted saltwater.