Concentrated aqueous sodium chloride solution in clays at thermodynamic conditions of hydraulic fracturing: Insight from molecular dynamics simulations

Abstract

A relationship between the effect of uni-univalent electrolytes on the structure of water and on its volatility The Journal of Chemical Physics 148, 222807 (2018) To address a high salinity of flow-back water during hydraulic fracturing, we use molecular dynamics (MD) simulations and study the thermodynamics, structure, and diffusion of concentrated aqueous salt solution in clay nanopores. The concentrated solution results from the dissolution of a cubic NaCl nanocrystal, immersed in an aqueous NaCl solution of varying salt concentration and confined in clay pores of a width comparable to the crystal size. The size of the nanocrystal equals to about 18 Å which is above a critical nucleus size. We consider a typical shale gas reservoir condition of 365 K and 275 bar, and we represent the clay pores as pyrophyllite and Na-montmorillonite (Na-MMT) slits. We employ the Extended Simple Point Charge (SPC/E) model for water, JoungCheatham model for ions, and CLAYFF for the slit walls. We impose the pressure in the normal direction and the resulting slit width varies from about 20 to 25 Å when the salt concentration in the surrounding solution increased from zero to an oversaturated value. By varying the salt concentration, we observe two scenarios. First, the crystal dissolves and its dissolution time increases with increasing salt concentration. We describe the dissolution process in terms of the number of ions in the crystal, and the crystal size and shape. Second, when the salt concentration reaches a system solubility limit, the crystal grows and attains a new equilibrium size; the crystal comes into equilibrium with the surrounding saturated solution. After crystal dissolution, we carry out canonical MD simulations for the concentrated solution. We evaluate the hydration energy, density profiles, orientation distributions, hydrogen-bond network, radial distribution functions, and in-plane diffusion of water and ions to provide insight into the microscopic behaviour of the concentrated aqueous sodium chloride solution in interlayer galleries of the slightly hydrophobic pyrophyllite and hydrophilic Na-MMT pores. Published by AIP Publishing. https://doi

Methane and carbon dioxide in dual‐porosity organic matter: Molecular simulations of adsorption and diffusion

“…51 force fields, respectively. The OPLS, SPC/E, and EPM2 force fields have been shown to satisfactorily reproduce the thermodynamic and vapor-liquid equilibrium properties, [52][53][54] additionally, the force fields are compatible with the CVFF. 47 The LJ crossinteractions were given by the Lorentz-Berthelot combining rules: 11…”

“…For the hydrocarbons, water, and carbon dioxide, we employed the Optimized Potential for Liquid Simulations (OPLS), 48 the flexible Extended Simple Point Charge (SPC/E), 49,50 and the flexible Elementary Physical Model (EPM2) 51 force fields, respectively. The OPLS, SPC/E, and EPM2 force fields have been shown to satisfactorily reproduce the thermodynamic and vapor–liquid equilibrium properties, 52‐54 additionally, the force fields are compatible with the CVFF 47 . The LJ cross‐interactions were given by the Lorentz‐Berthelot combining rules: 11 and σ ij = ( σ ii + σ jj )/2.…”

The Influence of Nanoporosity on the Behavior of Carbon-Bearing Fluids

“…For molecular simulations to reliably predict the properties of confined electrolytes, it is important that adequate force fields describe the interactions between water, electrolytes, and pore constituents, as well as all the cross-interactions between these components. 139,140 The reliability of the force fields used to describe the electrolytes depends on the model chosen to simulate water, and it is essential to tune the ion-ion interactions if one seeks to replicate the experimental salt solubility. 141,142 Svoboda and Lísal 140 recently employed the grand canonical Monte Carlo algorithm to predict the solubility of NaCl in montmorillonite pores.…”

Molecular Dynamics Simulation of Calcium Sulfate Nucleation in Homogeneous and Heterogeneous Crystallization Conditions: An Application in Water Flooding

“…The same amount of diffusivity for calcium and sulfate ions is due to strong ion association in a cluster structure. Svoboda and Lísal showed that ions’ self-diffusion coefficients in a solution change for different salinities. , In low-salinity cases, the small ions would experience higher hydration than the larger ions; therefore, they are more affected by surrounding water; hence, their mobility decreases. While in higher salinity solutions, ions’ self-diffusivities become similar because of the probability of ion association enhancement.…”