Hybridized method of pseudopotential lattice Boltzmann and cubic-plus-association equation of state assesses thermodynamic characteristics of associating fluids

Impact of nanopore confinement on phase behavior and enriched gas minimum miscibility pressure in asphaltenic tight oil reservoirs

Miscible gas injection in tight/shale oil reservoirs presents a complex problem due to various factors, including the presence of a large number of nanopores in the rock structure and asphaltene and heavy components in crude oil. This method performs best when the gas injection pressure exceeds the minimum miscibility pressure (MMP). Accordingly, accurate calculation of the MMP is of special importance. A critical issue that needs to be considered is that the phase behavior of the fluid in confined nanopores is substantially different from that of conventional reservoirs. The confinement effect may significantly affect fluid properties, flow, and transport phenomena characteristics in pore space, e.g., considerably changing the critical properties and enhancing fluid adsorption on the pore wall. In this study, we have investigated the MMP between an asphaltenic crude oil and enriched natural gas using Peng-Robinson (PR) and cubic-plus-association (CPA) equations of state (EoSs) by considering the effect of confinement, adsorption, the shift of critical properties, and the presence of asphaltene. According to the best of our knowledge, this is the first time a model has been developed considering all these factors for use in porous media. We used the vanishing interfacial tension (VIT) method and slim tube test data to calculate the MMP and examined the effects of pore radius, type/composition of injected gas, and asphaltene type on the computed MMP. The results showed that the MMP increased with an increasing radius of up to 100 nm and then remained almost constant. This is while the gas enrichment reduced the MMP. Asphaltene presence changed the trend of IFT reduction and delayed the miscibility achievement so that it was about 61% different from the model without the asphaltene precipitation effect. However, the type of asphaltene had little impact on the MMP, and the controlling factor was the amount of asphaltene in the oil. Moreover, although cubic EoSs are particularly popular for their simplicity and accuracy in predicting the behavior of hydrocarbon fluids, the CPA EoS is more accurate for asphaltenic oils, especially when the operating pressure is within the asphaltene precipitation range.

Improving the estimation accuracy of confined vapor-liquid equilibria by fine-tuning the pure component parameter in the PC-SAFT equation of state

We propose a thermodynamic model that combines the Young-Laplace equation and perturbed chain-statistical associating fluid theory (PC-SAFT) equation of state to estimate capillary condensation pressure in microporous and mesoporous sorbents. We adjust the PC-SAFT dispersion-energy parameter when the pore size becomes comparable to the molecular dimension. This modelling framework is applied to diverse systems containing associating and non-associating gases, various sorbents, and a wide range of temperatures. Our simulation results show that under extreme confinement, a higher value of the dispersion-energy parameter (ε) is required. Furthermore, using the experimental saturation pressure data for 18 different associating and non-associating confined fluids, we find that the shift in the PC-SAFT dispersion energy correlates with the ratio of the sorbent mean pore size to the PC-SAFT segment size (rp/σ). By fitting to the capillary condensation data, the relative deviation between the confined and bulk PC-SAFT dispersion energy parameter is only 0.1% at rp/σ = 15; however, this deviation starts to increase exponentially as rp/σ decreases. For a sorbent with large pores, when rp/σ > 15, the capillary condensation pressure results from our model are similar to the predictions from the Kelvin equation. Using a dataset containing 235 saturation pressure data points composed of 18 pure gases and 4 binary mixtures, the overall AARD% from our model is 12.26%, which verifies the good accuracy of our model. Because the mean sorbent pore radius (rp), the PC-SAFT energy parameter (ε), and segment size (σ) are known a priori, our model estimates the corrected energy parameter for small pores and, thus, extends its applicability.

Central-Moments-Based Lattice Boltzmann for Associating Fluids: A New Integrated Approach

Dynamic and thermodynamic behaviors of associating fluids play a crucial role in various science and engineering disciplines. Cubic plus association equation of state (CPA EOS) is implemented in a central-moments-based lattice Boltzmann method (LBM) in order to mimic the thermodynamic behavior of associating fluids. The pseudopotential approach is selected to model the multiphase thermodynamic characteristics such as reduced density of associating fluids. The priority of central-moments-based approach over multiple-relaxation-time collision operator is highlighted by performing double shear layers. The integration of central-moments-based LBM and CPA EOS is useful to simulate the dynamic and thermodynamic characteristics of associating fluids at high flow rate conditions, which is extended to high-density ratio scenarios by increasing the anisotropy order of gradient operator. In order to increase the stability of the model, a higher anisotropy order of the gradient operator is implemented; about 34 present reduction in spurious velocities is noticed in some cases. The type of gradient operator considerably affects the model thermodynamic consistency. Finally, the model is validated by observing a straight line in the Laplace law test. Prediction of thermodynamic behaviors of associating fluids is of significance in various applications including biological processes as well as fluid flow in porous media.