Search

The phase behavior of hydrocarbon fluids confined in porous media has been reported to deviate significantly from that in the bulk environment due to the existence of sub-10nm pores. Though experiments and simulations have measured the bubble/dew points and sorption isotherms of hydrocarbons confined in both natural and synthetic nanopores, the confinement effects in terms of the strength of fluid-pore interactions tuned by surface wettability and chemistry have received comparably less discussion. More importantly, the underlying physics of confinement-induced phenomena remain obfuscated. In this work, we studied the phase behavior and capillary condensation of n-hexane to understand the effects of confinement at the molecular level. To systematically investigate the pore effects, we constructed two types of wall confinements; one is a structureless virtual wall described by the Steele potential and the other one is an all-atom amorphous silica structure with surface modified by hydroxyl groups. Our numerical results demonstrated the importance of fluid-pore interaction, pore size, and pore morphology effects in mediating the pressure-volume-temperature (PVT) properties of hydrocarbons. The most remarkable finding of this work was that the saturation pressure predicted from the van der Waals-type adsorption isothermal loop could be elevated or suppressed relative to the bulk phase, as illustrated in the graphical abstract. As the surface energy (i.e., fluid-pore interaction) decreased, the isothermal vapor pressure increased, indicating a greater preference for the fluid to exist in the vapor state. Sufficient reduction of the fluid-pore interactions could even elevate the vapor pressure above that of the bulk fluid.

Accurate and efficient flow models for hydrocarbons are important in the development of enhanced geotechnical engineering for energy source recovery and carbon capture & storage in low-porosity, low-permeability rock formations. This work reports an atomistically-validated, mesoscopic model for heptane based on a many-body dissipative particle dynamics (mDPD) method. In this model, each heptane molecule is coarse-grained in one mDPD bead and the mDPD model parameters are calibrated with a rigorous approach using reference data, including experimental measurements and/or molecular dynamics (MD) simulations. Results show that this mDPD model accurately predicts the bulk pressure-density relation of heptane and surface tension. Notice that our approach can be used to calibrate the mDPD model for other hydrocarbons as well, though heptane is chosen as a representative source fluid for its abundance in source rocks. Further, our timing test indicates that the mDPD model is three orders of magnitude faster than its MD counterpart for simulations of bulk heptane in equivalent volumes. Overall, this work serves as a key prerequisite for the development of accurate and efficient mesoscale models for the flow of hydrocarbons confined in mesoporous rock formations., Comment: 16 pages, 7 figures

Mesoscopic simulations of hydrocarbon flow in source shales are challenging, in part due to the heterogeneous shale pores with sizes ranging from a few nanometers to a few micrometers. Additionally, the sub-continuum fluid-fluid and fluid-solid interactions in nano- to micro-scale shale pores, which are physically and chemically sophisticated, must be captured. To address those challenges, we present a GPU-accelerated package for simulation of flow in nano- to micro-pore networks with a many-body dissipative particle dynamics (mDPD) mesoscale model. Based on a fully distributed parallel paradigm, the code offloads all intensive workloads on GPUs. Other advancements, such as smart particle packing and no-slip boundary condition in complex pore geometries, are also implemented for the construction and the simulation of the realistic shale pores from 3D nanometer-resolution stack images. Our code is validated for accuracy and compared against the CPU counterpart for speedup. In our benchmark tests, the code delivers nearly perfect strong scaling and weak scaling (with up to 512 million particles) on up to 512 K20X GPUs on Oak Ridge National Laboratory's (ORNL) Titan supercomputer. Moreover, a single-GPU benchmark on ORNL's SummitDev and IBM's AC922 suggests that the host-to-device NVLink can boost performance over PCIe by a remarkable 40\%. Lastly, we demonstrate, through a flow simulation in realistic shale pores, that the CPU counterpart requires 840 Power9 cores to rival the performance delivered by our package with four V100 GPUs on ORNL's Summit architecture. This simulation package enables quick-turnaround and high-throughput mesoscopic numerical simulations for investigating complex flow phenomena in nano- to micro-porous rocks with realistic pore geometries.