Propane simulated in silica pores: Adsorption isotherms, molecular structure, and mobility

Adsorption Behaviors of Different Components of Shale Oil in Quartz Slits Studied by Molecular Simulation

Understanding the adsorption state and molecular behavior of the diverse components of shale oil in shale slits is of critical importance for exploring novel enhanced shale oil recovery techniques, but it is hard to be achieved by experimental measurements. In this paper, molecular dynamics (MD) simulations are performed to quantitatively describe the microbehavior of shale oil mixtures containing different kinds of hydrocarbon components, including asphaltene, in quartz slits. The spatial distributions of all the presenting components are given, the interaction energy between the components and quartz is analyzed, and the diffusion coefficients of all the components are calculated. It was found that asphaltene molecules play a vitally important role in restricting the detachment and diffusion movement of all hydrocarbon components, which is actually a key problem limiting the recovery efficiency of shale oil. The effects of temperature, slit aperture, and the appearance of CO₂ on the adsorption behavior of the different shale oil components are examined; the results suggest that the light and medium components are the fractions with the most potential in thermal exploitation, while injection of CO₂ is beneficial for the extraction of all the components, especially the medium components. This work gives insights into the effect of asphaltene on shale oil recovery in quartz slits and might provide guidance on the utilization of thermal and CO₂-enhanced enhanced oil recovery (EOR) techniques in shale oil production.

Effects of salinity and shear stress on clay deformation: A molecular dynamics study

The deformation of clay minerals is an important phenomenon that is relevant to many problems, particularly those that occur in subsurface geological formations. The salinity of the formations and external shear stress applied to them are two important factors that contribute to the deformation of such porous media. To gain a deeper understanding of such phenomena, we have carried out extensive molecular dynamics simulations using the Na-montmorillonite (Na-MMT) structure as the model of clay minerals and have studied the effect of salt concentration on its swelling. As the NaCl concentration increases, so also does the basal spacing. We demonstrate the effect of the coupling between the applied shear stress and NaCl salinity on the swelling behavior of Na-MMT, namely, deformation of the interlayer space that results in swelling. According to the results, the extent of Na-MMT deformation depends on both the brine salinity and the shear rate.

The role of the potential field on occurrence and flow of octane in quartz nanopores

Occurrence and flow of hydrocarbons in nanopores are two important issues in the effective exploitation of shale oil reservoirs. In this study, molecular dynamics simulations are employed to investigate the mechanisms about occurrence and flow of octane in slit-shaped quartz nanopores. We show that the occurrence state of octane and, therefore, its flow behavior are profoundly affected by the potential field from quartz walls and adsorption layers if the nanopore width w becomes less than 50 Å. Two main adsorption layers are always formed, adjacent to the walls and independent of w, due to two potential wells generated by the attractive potentials of the walls. Each pair of symmetrical adsorption layers, each of which can be considered as a solid-like surface, forms a confined environment similar to a nano-slit. The attractive potentials from them are found to be the cause for the formation of the adsorption layers between them. The obvious bulk phase of octane is formed in the pore of w = 50 Å due to the wide zero potential barrier induced by the innermost two adsorption layers. The nonlinear dependence of flow rate on pressure gradient shows that Darcy's law fails to describe the flow in the nanopore. The non-Darcy behavior mainly arises from adsorption effects from the walls and the adsorption layers, slippage between octane and walls and between adjacent two adsorption layers, and the molecular exchange between adsorption layers. A modified microscopic model is established to predict the dependence of flow rate on potential field, pressure gradient and w, which is in a good agreement with our simulation results and verified by the dodecane flow through the nanopore. Our work can be of great importance for revealing the mechanisms of occurrence and transport and guiding the estimation and exploitation of shale oil resources.

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.

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.

Effects of Confinement and Pressure on the Vibrational Behavior of Nano-Confined Propane

Fluids confined in nanopores exhibit significant deviations in their structure and dynamics from the bulk behavior. Although phase, structural, and diffusive behaviors of confined fluids have been investigated and reported extensively, confinement effects on the vibrational properties are less understood. We study the vibrational behavior of propane confined in 1.5 nm nanopores of MCM-41-S using inelastic neutron scattering (INS) and molecular dynamics (MD) simulations. Vibrational spectra have been obtained from INS data as functions of temperature and pressure. At ambient pressure, a strong quasielastic signal observed in the INS spectrum at 80 K suggests that confined propane remains liquid below the bulk phase melting point of 85 K. The quasielastic signal is heavily suppressed when either the pressure is increased to 1 kbar or the temperature is lowered to 30 K, indicating solidification of pore-confined propane. Confinement in MCM-41-S pores results in a glass-like state of propane that exhibits a relatively featureless low-energy vibrational spectrum compared to that of the bulk crystalline propane. Increasing the pressure to 3 kbar results in hardening of the intermolecular and methyl torsional modes. The INS data are used for estimating the isochoric specific heat of confined propane, which is compared with the specific heat of bulk propane reported in literature. Data from MD simulations are used to calculate the vibrational power spectra that agree qualitatively with the experimental data. Simulation data also suggest a reduction of the structural ordering (positional, orientational, and intramolecular) of propane under confinement.

Quasi-Two-Dimensional Phase Transition of Methane Adsorbed in Cylindrical Silica Mesopores

Using Monte Carlo and molecular dynamics simulations, we examine the adsorption of methane in cylindrical silica mesopores in an effort to understand a possible phase transition of adsorbed methane in MCM-41 and SBA-15 silica that was previously identified by an unexpected increase in the adsorbed fluid density following capillary condensation, as measured by small-angle neutron scattering (SANS) [Chiang, W-S., et al., Langmuir 2016, 32, 8849]. Our initial simulation results identify a roughly 10 % increase in the density of the liquidlike adsorbed phase for either an isotherm with increasing pressure or an isobar with decreasing temperature and that this densification is associated with a local maximum in the isosteric enthalpy of adsorption. Subsequent analysis of the simulated fluid, via computation of bond-orientational order parameters of specific annular layers of the adsorbed fluid, showed that the layers undergo an ordering transition from a disordered, amorphous state to one with two-dimensional hexagonal structure. Furthermore, this two-dimensional restructuring of the fluid occurs at the same thermodynamic state points as the aforementioned densification and local maximum in the isosteric enthalpy of adsorption. We thus conclude that the densification of the fluid is the result of structural reorganization, which is signaled by the maximum in the isosteric enthalpy. Owing to the qualitative similarity of the structural transitions in the simulated and experimental methane fluids, we propose this hexagonal reorganization as a plausible explanation of the densification observed in SANS measurements. Lastly, we speculate how this structural transition may impact the transport properties of the adsorbed fluid.

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.

Molecular dynamics simulations of propane in slit shaped silica nano-pores: direct comparison with quasielastic neutron scattering experiments

MD simulations reveal the origin of anomalous pressure dependence of propane diffusion in silica mesopores.