Publisher's copyright statement:This document is the Accepted Manuscript version of a Published Work that appeared in nal form in Journal of physical chemistry C, copyright c 2015 American Chemical Society after peer review and technical editing by the publisher. To access the nal edited and published work see http://dx.doi.org/10.1021/acs.jpcc.5b00555Additional information: Use policyThe full-text may be used and/or reproduced, and given to third parties in any format or medium, without prior permission or charge, for personal research or study, educational, or not-for-pro t purposes provided that:• a full bibliographic reference is made to the original source • a link is made to the metadata record in DRO • the full-text is not changed in any way The full-text must not be sold in any format or medium without the formal permission of the copyright holders.Please consult the full DRO policy for further details. AbstractEnhanced oil recovery is becoming commonplace in order to maximise recovery from oilfields. One of these methods, low-salinity enhanced oil recovery (EOR) has shown promise, however the fundamental underlying chemistry requires elucidating. Here, three mechanisms proposed to account for low-salinity enhanced oil recovery in sandstone reservoirs are investigated using molecular dynamic simulations. The mechanisms probed are electric double layer expansion, multicomponent ionic exchange and pH effects arising at clay mineral surfaces. Simulations of smectite basal planes interacting with uncharged non-polar decane, uncharged polar decanoic acid and charged Nadecanoate model compounds are used to this end. Various salt concentrations of NaCl are modelled: 0% , 1% , 5% and 35% to determine the role of salinity upon the three separate mechanisms. Furthermore, the initial oil/water-wetness of the clay surface is modeled. Results show that electric double layer expansion is not able to fully explain the effects of low-salinity enhanced oil recovery. The pH surrounding a clays basal plane, and hence the protonation and charge of acid molecules is determined to be one of the dominant effects driving low-salinity EOR. Further, results present that the presence of calcium cations can drastically alter the oil wettability of a clay ⇤To whom correspondence should be addressed mineral surface. Replacing all divalent cations with monovalent cations through multicomponent cation exchange dramatically increases the water wettability of a clay surface, and will increase EOR.

Fine-grained sediments and sedimentary rocks play important roles in a variety of modern energy technologies from petroleum geology to geological carbon sequestration and radioactive waste management. However, despite their utility and ubiquity, many of their properties remain poorly understood. In particular, the ability to predict the permeability and mechanics of these media remains a persistent fundamental challenge in the geosciences. In the present work, we show how large-scale classical molecular dynamics (MD) simulations can help interpret the properties of fine-grained sedimentary material. All-atom simulations containing 30 discrete clay particles are utilized to understand the evolution of a clay nanoparticle suspension during its progressive dehydration. Microstructural (pore size distribution, tortuosity, anisotropy), thermodynamic (enthalpy and free energy of hydration, anion exclusion), mechanical (total suction), and transport properties (diffusion coefficient tensors of water and sodium) are calculated and compared to the experiment. Overall, our results provide new insight into the coupled chemistry, mechanics, and transport properties of disordered nanoparticle assemblages and shed light upon the important role of water films in controlling these properties.

No abstract

In this study, classical molecular dynamic simulations have been used to examine the molecular properties of the water-alkane interface at various NaCl salt concentrations (up to 3.0 mol/kg). A variety of different force field combinations have been compared against experimental surface/interfacial tension values for the water-vapour, decane-vapour and water-decane interfaces. Six different force fields for water (SPC, SPC/E, TIP3P, TIP3Pcharmm, TIP4P & TIP4P2005), and three further force fields for alkane (TraPPE-UA, CGenFF & OPLS) have been compared to experimental data. CGenFF, OPLS-AA and TraPPE-UA all accurately reproduce the interfacial properties of decane. The TIP4P2005 (four-point) water model is shown to be the most accurate water model for predicting the interfacial properties of water. The SPC/E water model is the best three-point parameterisation of water for this purpose. The CGenFF and TraPPE parameterisations of oil accurately reproduce the interfacial tension with water using either the TIP4P2005 or SPC/E water model. The salinity dependence on surface/interfacial tension is accurately captured using the Smith & Dang parameterisation of NaCl. We observe that the Smith & Dang model slightly overestimates the surface/interfacial tensions at higher salinities (>1.5 mol/kg). This is ascribed to an overestimation of the ion exclusion at the interface.

Many important properties of clay minerals are defined by the species of chargebalancing cation. Phenomena such as clay swelling and cation exchange depend upon the cation species present and it is therefore important to understand how the cations bind with the mineral surface at a fundamental level. In the present study, the binding affinities of several different charge-balancing cations with the basal surface of the smectite clay montmorillonite have been calculated using molecular dynamics in conjunction with the well-tempered metadynamics algorithm. The results follow a Hofmeister series of preferred ion adsorption to the smectite basal surfaces of the form: