Molecular dynamics of wetting layer formation and forced water invasion in angular nanopores with mixed wettability

Molecular Dynamics Insight into the CO2 Flooding Mechanism in Wedge-Shaped Pores

Because of the growing demand for energy, oil extraction under complicated geological conditions is increasing. Herein, oil displacement by CO2 in wedge-shaped pores was investigated by molecular dynamics simulation. The results showed that, for both single and double wedge-shaped models, pore Ⅱ (pore size from 3 to 8 nm) exhibited a better CO2 flooding ability than pore Ⅰ (pore size from 8 to 3 nm). Compared with slit-shaped pores (3 and 8 nm), the overall oil displacement efficiency followed the sequence of 8 nm > double pore Ⅱ > single pore Ⅱ > 3 nm > double pore Ⅰ > single pore Ⅰ, which confirmed that the exits of the wedge-shaped pores had determinant effects on CO2 enhanced oil recovery over their entrances. “Oil/CO2 inter-pore migration” and “siphoning” phenomena occurred in wedge-shaped double pores by comparing the volumes of oil/CO2 and the center of mass. The results of the interaction and radial distribution function analyses indicate that the wide inlet and outlet had a larger CO2−oil contact surface, better phase miscibility, higher interaction, and faster displacement. These findings clarify the CO2 flooding mechanisms in wedge-shaped pores and provide a scientific basis for the practical applications of CO2 flooding.

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.

Surfactant-induced wettability reversal on oil-wet calcite surfaces: Experimentation and molecular dynamics simulations with scaled-charges

HYPOTHESIS

Surfactant flooding is the leading approach for reversing the wettability of oil-wet carbonate reservoirs, which is critical for the recovery of the remaining oil. Combination of molecular dynamics (MD) simulations with experiments on simplified model systems can uncover the molecular mechanisms of wettability reversal and identify key molecular properties for systematic design of new, effective chemical formulations for the enhanced oil recovery.

EXPERIMENTS/SIMULATIONS

Wettability reversal by a series of surfactant solutions was studied experimentally using contact angle measurements on aged calcite chips, and a novel MD simulation methodology with scaled-charges that provides superior description of the ionic interactions in aqueous solutions.

FINDINGS

The MD simulation results were in excellent agreement with the experiments. Cationic surfactants were the most effective in reversing the calcite wettability, resulting in complete detachment of the oil from the surface. Some nonionic surfactants also altered the wettability, but to a lesser degree, while the amphoteric and anionic surfactants had no effect. From the tested cationic surfactants, the double-tailed one was the least effective, but the experiments were inconclusive due to its poor solubility. Contributions of specific interactions to the wettability reversal process and implications for the design and optimization of surfactants for the enhanced oil recovery are discussed.

Influence of organic molecules on wetting characteristics of mica/H2/brine systems: Implications for hydrogen structural trapping capacities

HYPOTHESIS

Actualization of the hydrogen (H₂) economy and decarbonization goals can be achieved with feasible large-scale H₂ geo-storage. Geological formations are heterogeneous, and their wetting characteristics play a crucial role in the presence of H₂, which controls the pore-scale distribution of the fluids and sealing capacities of caprocks. Organic acids are readily available in geo-storage formations in minute quantities, but they highly tend to increase the hydrophobicity of storage formations. However, there is a paucity of data on the effects of organic acid concentrations and types on the H₂-wettability of caprock-representative minerals and their attendant structural trapping capacities.

EXPERIMENT

Geological formations contain organic acids in minute concentrations, with the alkyl chain length ranging from C₄ to C₂₆. To fully understand the wetting characteristics of H₂ in a natural geological picture, we aged mica mineral surfaces as a representative of the caprock in varying concentrations of organic molecules (with varying numbers of carbon atoms, lignoceric acid C₂₄, lauric acid C₁₂, and hexanoic acid C₆) for 7 days. To comprehend the wettability of the mica/H₂/brine system, we employed a contact-angle procedure similar to that in natural geo-storage environments (25, 15, and 0.1 MPa and 323 K).

FINDINGS

At the highest investigated pressure (25 MPa) and the highest concentration of lignoceric acid (10^-2 mol/L), the mica surface became completely H₂ wet with advancing (θ_a= 106.2°) and receding (θ_r=97.3°) contact angles. The order of increasing θ_a and θ_r with increasing organic acid contaminations is as follows: lignoceric acid > lauric acid > hexanoic acid. The results suggest that H₂ gas leakage through the caprock is possible in the presence of organic acids at higher physio-thermal conditions. The influence of organic contamination inherent at realistic geo-storage conditions should be considered to avoid the overprediction of structural trapping capacities and H₂ containment security.

Chemical computational approaches for optimization of effective surfactants in enhanced oil recovery

The surfactant flooding becomes an attractive method among several Enhanced Oil Recovery (EOR) processes to improve the recovery of residual oil left behind in the reservoir after secondary oil recovery process. The designing of a new effective surfactant is a comparatively complex and often time consuming process as well as cost-effective due to its dependency on the crude oil and reservoir properties. An alternative chemical computational approach is focused in this article to optimize the performance of effective surfactant system for EOR. The molecular dynamics (MD), dissipative particle dynamics (DPD) and density functional theory (DFT) simulations are mostly used chemical computational approaches to study the behaviour in multiple phase systems like surfactant/oil/brine. This article highlighted a review on the impact of surfactant head group structure on oil/water interfacial property like interfacial tensions, interface formation energy, interfacial thickness by MD simulation. The effect of entropy in micelle formation has also discussed through MD simulation. The polarity, dipole moment, charge distribution and molecular structure optimization have been illustrated by DFT. A relatively new coarse-grained method, DPD is also emphasized the phase behaviour of surfactant/oil/brine as well as polymer-surfactant complex system.

Effects of Surfactant Charge and Molecular Structure on Wettability Alteration of Calcite: Insights from Molecular Dynamics Simulations

Wettability alteration of oil-wet calcite by surfactants was studied by means of molecular dynamics (MD) simulations. The simulations use the recently developed model for positively charged calcite surface, whose oil-wet state originates from binding of negatively charged carboxylate molecules contained in the oil, consistently with the bulk of the available experimental data. The ability to alter the surface wettability, which can be directly quantified by the release of the surface-bound carboxylates, is tested for nine different surfactants of all charge types-cationic, anionic, nonionic, and zwitterionic-and compared to that of brine. It was found that only the cationic surfactants are able to detach the organic carboxylates more efficiently than brine, while the neutral and anionic surfactants do not seem to have any measurable effect on the wettability. The outperformance of the cationic surfactants is generally consistent with the majority of previously published experimental observations. The data also point toward a consistently better performance of single-tailed cationic surfactants over the two-tailed structure. Molecular mechanism of the wettability alteration by different types of surfactants is discussed, along with the implications of the results for the design of new surfactant formulations for the enhanced oil recovery.

Properties of a novel imbibition polymer with ultra-high wetting ability as a fracturing fluid system

A fracturing fluid with wetting ability was obtained by a crosslinking strategy. The ultra-high wettability transformed a sandstone surface from oil-wet to water-wet.

Atomistic Molecular Dynamics Simulations of Crude Oil/Brine Displacement in Calcite Mesopores

Unconventional reservoirs such as hydrocarbon-bearing shale formations and ultratight carbonates generate a large fraction of oil and gas production in North America. The characteristic feature of these reservoirs is their nanoscale porosity that provides significant surface areas between the pore walls and the occupying fluids. To better assess hydrocarbon recovery from these formations, it is crucial to develop an improved insight into the effects of wall-fluid interactions on the interfacial phenomena in these nanoscale confinements. One of the important properties that controls the displacement of fluids inside the pores is the threshold capillary pressure. In this study, we present the results of an integrated series of large-scale molecular dynamics (MD) simulations performed to investigate the effects of wall-fluid interactions on the threshold capillary pressures of oil-water/brine displacements in a calcite nanopore with a square cross section. Fully atomistic models are utilized to represent crude oil, brine, and calcite in order to accommodate electrostatic interactions and H-bonding between the polar molecules and the calcite surface. To this end, we create mixtures of various polar and nonpolar organic molecules to better represent the crude oil. The interfacial tension between oil and water/brine and their contact angle on calcite surface are simulated. We study the effects of oil composition, water salinity, and temperature and pressure conditions on these properties. The threshold capillary pressure values are also obtained from the MD simulations for the calcite nanopore. We then compare the MD results against those generated using the Mayer-Stowe-Princen (MSP) method and explain the differences.

Pore diameter effects on phase behavior of a gas condensate in graphitic one-and two-dimensional nanopores

Molecular dynamics (MD) simulations were performed on a hydrocarbon mixture representing a typical gas condensate composed mostly of methane and other small molecules with small fractions of heavier hydrocarbons, representative of mixtures found in tight shale reservoirs. The fluid was examined both in bulk and confined to graphitic nano-scale slits and pores. Numerous widths and diameters of slits and pores respectively were examined under variable pressures at 300 K in order to find conditions in which the fluid at the center of the apertures would not be affected by capillary condensation due to the oil-wet walls. For the bulk fluid, retrograde phase behavior was verified by liquid volumes obtained from Voronoi tessellations. In cases of both one and two-dimensional confinement, for the smallest apertures, heavy molecules aggregated inside the pore space and compression of the gas outside the solid structure lead to decreases in density of the confined fluid. Normal density/pressure relationships were observed for slits having gaps of above 3 nm and pores having diameters above 6 nm. At 70 bar, the minimum gap width at which the fluid could pass through the center of slits without condensation effects was predicted to be 6 nm and the corresponding diameter in pores was predicted to be 8 nm. The models suggest that in nanoscale networks involving pores smaller than these limiting dimensions, capillary condensation should significantly impede transmission of natural gases with similar composition.

Hydraulic transport across hydrophilic and hydrophobic nanopores: Flow experiments with water and n-hexane

We experimentally explore pressure-driven flow of water and n-hexane across nanoporous silica (Vycor glass monoliths with 7- or 10-nm pore diameters, respectively) as a function of temperature and surface functionalization (native and silanized glass surfaces). Hydraulic flow rates are measured by applying hydrostatic pressures via inert gases (argon and helium, pressurized up to 70 bar) on the upstream side in a capacitor-based membrane permeability setup. For the native, hydrophilic silica walls, the measured hydraulic permeabilities can be quantitatively accounted for by bulk fluidity provided we assume a sticking boundary layer, i.e., a negative velocity slip length of molecular dimensions. The thickness of this boundary layer is discussed with regard to previous capillarity-driven flow experiments (spontaneous imbibition) and with regard to velocity slippage at the pore walls resulting from dissolved gas. Water flow across the silanized, hydrophobic nanopores is blocked up to a hydrostatic pressure of at least 70 bar. The absence of a sticking boundary layer quantitatively accounts for an enhanced n-hexane permeability in the hydrophobic compared to the hydrophilic nanopores.