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

Study on the Effects of Wettability and Pressure in Shale Matrix Nanopore Imbibition during Shut-in Process by Molecular Dynamics Simulations

Shut-in after fracturing is generally adopted for wells in shale oil reservoirs, and imbibition occurring in matrix nanopores has been proven as an effective way to improve recovery. In this research, a molecular dynamics (MD) simulation was used to investigate the effects of wettability and pressure on nanopore imbibition during shut-in for a typical shale reservoir, Jimsar. The results indicate that the microscopic advancement mechanism of the imbibition front is the competitive adsorption between "interfacial water molecules" at the imbibition front and "adsorbed oil molecules" on the pore wall. The essence of spontaneous imbibition involves the adsorption and aggregation of water molecules onto the hydroxyl groups on the pore wall. The flow characteristics of shale oil suggest that the overall push of the injected water to the oil phase is the main reason for the displacement of adsorbed oil molecules. Thus, shale oil, especially the heavy hydrocarbon component in the adsorbed layer, tends to slip on the walls. However, the weak slip ability of heavy components on the wall surface is an important reason that restricts the displacement efficiency of shale oil during spontaneous imbibition. The effectiveness of spontaneous imbibition is strongly dependent on the hydrophilicity of the matrix pore's wall. The better hydrophilicity of the matrix pore wall facilitates higher levels of adsorption and accumulation of water molecules on the pore wall and requires less time for "interfacial water molecules" to compete with adsorbed oil molecules. During the forced imbibition process, the pressure difference acts on both the bulk oil and the boundary adsorption oil, but mainly on the bulk oil, which leads to the occurrence of wetting hysteresis. Meanwhile, shale oil still existing in the pore always maintains a good, stratified adsorption structure. Because of the wetting hysteresis phenomenon, as the pressure difference increases, the imbibition effect gradually increases, but the actual capillary pressure gradually decreases and there is a loss in the imbibition velocity relative to the theoretical value. Simultaneously, the decline in hydrophilicity further weakens the synergistic effect on the imbibition of the pressure difference because of the more pronounced wetting hysteresis. Thus, selecting an appropriate well pressure enables cost savings and maximizes the utilization of the formation's natural power for enhanced oil recovery (EOR).

Surface Interactions and Nanoconfinement of Methane and Methane plus CO2 Revealed by High-Pressure Magic Angle Spinning NMR Spectroscopy and Molecular Dynamics

This study explores the fundamental, molecular- to microscopic-level behavior of methane gas confined into nanoporous silica proxies with different pore diameters and surface-to-volume (S/V) ratios. Surfaces and pore walls of nanoporous silica matrices are decorated with hydroxyl (-OH) groups, resembling natural heterogeneity. High-pressure MAS NMR was utilized to characterize the interactions between methane and the engineered nanoporous silica proxies under various temperature and pressure regimes. There was a change in the chemical shift position of confined methane slightly in the mixtures with nanoporous silica up to 393 K, as shown by high-pressure ¹³C-NMR. The ¹³C-NMR chemical shift of methane was changed by pressure, explained by the densification of methane inside the nanoporous silica materials. The influence of pore diameter and S/V of the nanoporous silica materials on the behaviors and dynamics of methane were studied. The presence of CO₂ in mixtures of silica and methane needs analysis with caution because CO₂ in a supercritical state and gaseous CO₂ change the original structure of nanoporous silica and change surface area and pore volume. According to simulation, the picosecond scale dynamics of methane confined in larger pores of amorphous silica is faster. In the 4 nm pore, the diffusivity obtained from MD simulations in the pore with a higher S/V ratio is slower due to the trapping of methane molecules in adsorbed layers close to the corrugated pore surface. In contrast, relaxation measured with NMR for smaller pores (higher S/V) exhibits larger T₁, indicating slower relaxation.

Breakthrough pressure of oil displacement by water through the ultra-narrow kerogen pore throat from the Young-Laplace equation and molecular dynamic simulations

As one common unconventional reservoir, shale plays a pivotal role to compensate the depletion of conventional oil resources. There are numerous nanoscale pores and ultra-narrow pore throats (sub 2 nm) in shale media. To displace oil through ultra-narrow pore throats by water, one needs to overcome excessively-high capillary pressure. Understanding the water-oil two-phase displacement process through pore throats is critical to numerical simulation on tight/shale oil exploitation and ultimate oil recovery estimation. In this work, we use molecular dynamics simulations to investigate oil (represented by n-octane) displacement by water through a ~2 nm kerogen (represented by Type II-C kerogen) pore throat. Besides, the applicability of the Young-Laplace equation to the ultra-narrow kerogen pore throat has been assessed. We find that although the Type II-C kerogen is generally oil-wet, water has an excellent displacement efficiency without the oil film on the substrate, thanks to the hydrogen bonding formed between water and heteroatoms (such as O, N, and S) on the kerogen surface. Unlike previous studies, the capillary pressure obtained from the widely used Young-Laplace equation shows a good agreement with the breakthrough pressure obtained from MD simulations for the ∼2 nm kerogen pore throat. Our work indicates that explicitly considering intermolecular interactions as well as atomistic and molecular level characteristics is imperative to study the two-phase displacement process through ultra-narrow pore throats.

Effects of pore connectivity and tortuosity on the dynamics of fluids confined in sub-nanometer pores

Dynamical behavior of fluids under nano-pore confinement is studied extensively as it has important implications for several industrial as well as geological processes. Pore network in many porous materials exhibits a varied degree of inter connections. The extent of this pore connectivity may affect the structural and dynamical behavior of the confined fluid. However, studies of fluid confinement addressing these effects systematically are lacking. Here, we report molecular dynamics simulation studies addressing the effects of pore connectivity on the dynamics of two representative fluids - CO₂ and ethane in silicalite by systematically varying the degree of pore connectivity through selectively blocking some pore space with immobile methane molecules. By selectively turning off the pore spaces in the shape of straight, or tortuous zigzag channels, we also probe the effects of pore tortuosity. In general, pore connectivity is found to facilitate both the translational as well as rotational dynamics of both fluids, while the intermolecular modes of vibration in both fluids remain largely unaffected. The effects of providing connections between a set of straight or zigzag channel-like pores are however more nuanced. Pore tortuosity facilitates the rotational motion, but suppresses the translational motion of CO₂, while its effects on the rotational and translational motion of ethane are less pronounced. The intermolecular vibrational modes of both fluids shift to higher energies with an increase in the number of tortuous pores. The results reported here provide a detailed molecular level understanding of the effects of pore connectivity on the dynamics of fluids and thus have implications for applications like fluid separation.

Using Computationally-Determined Properties for Machine Learning Prediction of Self-Diffusion Coefficients in Pure Liquids

The ability to predict transport properties of liquids quickly and accurately will greatly improve our understanding of fluid properties both in bulk and complex mixtures, as well as in confined environments. Such information could then be used in the design of materials and processes for applications ranging from energy production and storage to manufacturing processes. As a first step, we consider the use of machine learning (ML) methods to predict the diffusion properties of pure liquids. Recent results have shown that Artificial Neural Networks (ANNs) can effectively predict the diffusion of pure compounds based on the use of experimental properties as the model inputs. In the current study, a similar ANN approach is applied to modeling diffusion of pure liquids using fluid properties obtained exclusively from molecular simulations. A diverse set of 102 pure liquids is considered, ranging from small polar molecules (e.g., water) to large nonpolar molecules (e.g., octane). Self-diffusion coefficients were obtained from classical molecular dynamics (MD) simulations. Since nearly all the molecules are organic compounds, a general set of force field parameters for organic molecules was used. The MD methods are validated by comparing physical and thermodynamic properties with experiment. Computational input features for the ANN include physical properties obtained from the MD simulations as well as molecular properties from quantum calculations of individual molecules. Fluid properties describing the local liquid structure were obtained from center of mass radial distribution functions (COM-RDFs). Feature sensitivity analysis revealed that isothermal compressibility, heat of vaporization, and the thermal expansion coefficient were the most impactful properties used as input for the ANN model to predict the MD simulated self-diffusion coefficients. The MD-based ANN successfully predicts the MD self-diffusion coefficients with only a subset (2 to 3) of the available computationally determined input features required. A separate ANN model was developed using literature experimental self-diffusion coefficients as model targets. Although this second ML model was not as successful due to a limited number of data points, a good correlation is still observed between experimental and ML predicted self-diffusion coefficients.

Consequences of Cylindrical Pore Geometry for Interfacial Phenomena in Reversed-Phase Liquid Chromatography

The interfacial phenomena behind analyte separation in a reversed-phase liquid chromatography column take place nearly exclusively inside the silica mesopores. Their cylindrical geometry can be expected to shape the properties of the chromatographic interface with consequences for the analyte density distribution and diffusivity. To investigate this topic through molecular dynamics simulations, we introduce a cylindrical pore inside a slit pore configuration, where the inner curved and outer planar silica surface bear the same bonded phase. The present model replicates an average-sized (9 nm) mesopore in an endcapped C₁₈ column equilibrated with a mobile phase of 70/30 (v/v) water/acetonitrile. Simulations performed for ethylbenzene and acetophenone show that the surface curvature shifts the bonded phase and analyte density toward the pore center, decreases the solvent density in the bonded-phase region, increases the acetonitrile excess in the interfacial region, and considerably enhances the surface diffusivity of both analytes. Overall, the cylindrical pore provides a more hydrophobic environment than the slit pore. Ethylbenzene density is decidedly increased in the cylindrical pore, whereas acetophenone density is nearly equally distributed between the cylindrical and slit pore. The cylindrical pore geometry thus sharpens the discrimination between the apolar and moderately polar analytes while enhancing the mass transport of both.

Structure of ice confined in silica nanopores

Observed anomalous thermodynamic properties of confined water such as deviations in the melting point and freezing point motivate the determination of the structure of confined water as a function of pore size and temperature. In this study, we investigate the dynamic evolution of the structure of confined ice in SBA-15 porous materials with pore diameters of 4 nm, 6 nm, and 8 nm at temperatures ranging from 183 K to 300 K using in operando Wide-Angle X-Ray Scattering (WAXS) measurements, X-Ray Partial Distribution Function (PDF) measurements, and classical Molecular Dynamics (MD) simulations. Formation of hexagonal ice structures is noted in all the three pore sizes. In silica nanopores with diameters of 4 nm, cubic ice formation is noted in addition to hexagonal ice. Longer lasting hydrogen bonds and longer residence times of the water molecules in the first coordination shell contribute to observed crystalline organization of ice in confinement. Self-diffusion coefficients of confined liquid water, predicted from classical MD simulations, are four orders of magnitude higher compared to ice formed in confinement. These experimental and simulation results provide comprehensive insights underlying the organization of confined water and ice in silica nanopores and the underlying physico-chemical interactions that contribute to the observed structures.

Diffusion, Interactions, and Disparate Kinetic Trapping of Water-Hydrocarbon Mixtures in Nanoporous Solids

Mixed fluids confined in porous solid hosts present challenges for the accurate characterization of individual-component behavior. NMR diffusometry with chemical resolution is used to identify unexpected loading- and composition-dependent anomalous diffusion in water/cyclohexane mixtures confined to solid nanoporous glass (NPG) hosts. Diffusion NMR results indicate that data obtained on pure-component liquids in confinement cannot be extrapolated to their nonideal liquid mixtures confined in the same solid host. Loading-dependent data must be obtained on each component in the confined mixture in order to determine which of the liquid components exhibits chemical affinity for the host and, conversely, which of the components exhibits anomalous diffusivity. Most notably, NMR diffusometry revealed that cyclohexane diffusivity varied by 2 orders of magnitude in a water-rich mixture depending on the total fluid loading in the NPG host, ranging from anomalously high diffusivities that significantly exceeded that for pure cyclohexane in NPG at low fluid loadings to kinetically trapped sequestration at high fluid loadings. NMR diffusometry indicates that nonideal solution behavior in fluids confined within nanoporous hosts may have practical implications for enhanced oil recovery methods. Specifically, kinetic trapping of hydrocarbons in water-flooding regimes can result from complex liquid-vapor equilibrium that is significantly perturbed from that which exists in bulk or microporous confinement.

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.

A Novel Modeling Approach to Stochastically Evaluate the Impact of Pore Network Geometry, Chemistry and Topology on Fluid Transport

Fine-grained sandstones, siltstones, and shales have become increasingly important to satisfy the ever-growing global energy demands. Of particular current interest are shale rocks, which are mudstones made up of organic and inorganic constituents of varying pore sizes. These materials exhibit high heterogeneity, low porosity, varying chemical composition and low pore connectivity. Due to the complexity and the importance of such materials, many experimental, theoretical and computational efforts have attempted to quantify the impact of rock features on fluids diffusivity and ultimately on permeability. In this study, we introduce a stochastic kinetic Monte Carlo approach developed to simulate fluid transport. The features of this approach allow us to discuss the applicability of 2D vs 3D models for the calculation of transport properties. It is found that a successful model should consider realistic 3D pore networks consisting of pore bodies that communicate via pore throats, which however requires a prohibitive amount of computational resources. To overcome current limitations, we present a rigorous protocol to stochastically generate synthetic 3D pore networks in which pore features can be isolated and varied systematically and individually. These synthetic networks do not correspond to real sample scenarios but are crucial to achieve a systematic evaluation of the pore features on the transport properties. Using this protocol, we quantify the contribution of the pore network’s connectivity, porosity, mineralogy, and pore throat width distribution on the diffusivity of supercritical methane. A sensitivity analysis is conducted to rank the significance of the various network features on methane diffusivity. Connectivity is found to be the most important descriptor, followed by pore throat width distribution and porosity. Based on such insights, recommendations are provided on possible technological approaches to enhance fluid transport through shale rocks and equally complex pore networks. The purpose of this work is to identify the significance of various pore network characteristics using a stochastic KMC algorithm to simulate the transport of fluids. Our findings could be relevant for applications that make use of porous media, ranging from catalysis to radioactive waste management, and from environmental remediation to shale gas production.

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.