1
|
Chen X, Tang L, Jia C, Yue P, Zhang Z, Liu W. Phase Behavior of Hydrocarbon Fluids in Shale Reservoirs, Considering Pore Geometries, Adsorption, and Water Film. ACS OMEGA 2024; 9:2104-2112. [PMID: 38250389 PMCID: PMC10795131 DOI: 10.1021/acsomega.3c03601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/23/2023] [Revised: 10/01/2023] [Accepted: 10/04/2023] [Indexed: 01/23/2024]
Abstract
Phase behavior of hydrocarbon fluids in nanopores is different from that observed in a PVT cell due to the confinement effect. While scholars have established various models for studying the phase behavior in nanopores, the authors often ignore the effect of pore geometries, which can significantly affect the critical fluid properties in shale nanopores. In this study, we extend the Soave-Redlich-Kwong equation of state (SRK EOS) using potential theory and establish models of critical property shift, considering pore geometries, adsorption, and water film. Our research shows that the critical property shifts, considering fluid adsorption, begin at rp ≤ 10 nm and are seriously strengthened with nanopore radius reduction. The extended SRK EOS is applied to compute phase diagrams of the 50% C1-50% nC10 mixture at different pore sizes and find that the thickness of adsorption and water film causes a depression in the P-T diagram and that the bubble point pressure is lower in cylindrical pores. At pressures above 6 MPa, the irreducible water saturation and pore geometries greatly impact the vapor-liquid ratio. This study is significant for evaluating residual oil distribution and studying fluid flow laws in shale reservoirs.
Collapse
Affiliation(s)
- Xiaofan Chen
- State
Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest Petroleum University, Chengdu 610500, China
| | - Liandong Tang
- State
Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest Petroleum University, Chengdu 610500, China
| | - Chunsheng Jia
- State
Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest Petroleum University, Chengdu 610500, China
| | - Ping Yue
- State
Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest Petroleum University, Chengdu 610500, China
| | - Zhenzu Zhang
- Huanqing
Oil Production Plant of PetroChina Yumen Oilfield Branch, Jiuquan 735202, Gansu, China
| | - Wei Liu
- Huanqing
Oil Production Plant of PetroChina Yumen Oilfield Branch, Jiuquan 735202, Gansu, China
| |
Collapse
|
2
|
Li Z, Turner J, Snurr RQ. Computational investigation of hysteresis and phase equilibria of n-alkanes in a metal-organic framework with both micropores and mesopores. Commun Chem 2023; 6:90. [PMID: 37156883 PMCID: PMC10167368 DOI: 10.1038/s42004-023-00889-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2022] [Accepted: 04/19/2023] [Indexed: 05/10/2023] Open
Abstract
Adsorption hysteresis is a phenomenon related to phase transitions that can impact applications such as gas storage and separations in porous materials. Computational approaches can greatly facilitate the understanding of phase transitions and phase equilibria in porous materials. In this work, adsorption isotherms for methane, ethane, propane, and n-hexane were calculated from atomistic grand canonical Monte Carlo (GCMC) simulations in a metal-organic framework having both micropores and mesopores to better understand hysteresis and phase equilibria between connected pores of different size and the external bulk fluid. At low temperatures, the calculated isotherms exhibit sharp steps accompanied by hysteresis. As a complementary simulation method, canonical (NVT) ensemble simulations with Widom test particle insertions are demonstrated to provide additional information about these systems. The NVT+Widom simulations provide the full van der Waals loop associated with the sharp steps and hysteresis, including the locations of the spinodal points and points within the metastable and unstable regions that are inaccessible to GCMC simulations. The simulations provide molecular-level insight into pore filling and equilibria between high- and low-density states within individual pores. The effect of framework flexibility on adsorption hysteresis is also investigated for methane in IRMOF-1.
Collapse
Affiliation(s)
- Zhao Li
- Department of Chemical and Biological Engineering, Northwestern University, 2145 Sheridan Road, Evanston, IL, 60208, USA
| | - Jake Turner
- Department of Chemical and Biological Engineering, Northwestern University, 2145 Sheridan Road, Evanston, IL, 60208, USA
| | - Randall Q Snurr
- Department of Chemical and Biological Engineering, Northwestern University, 2145 Sheridan Road, Evanston, IL, 60208, USA.
| |
Collapse
|
3
|
Machine-Learned Free Energy Surfaces for Capillary Condensation and Evaporation in Mesopores. ENTROPY 2022; 24:e24010097. [PMID: 35052123 PMCID: PMC8774451 DOI: 10.3390/e24010097] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/15/2021] [Revised: 12/29/2021] [Accepted: 01/05/2022] [Indexed: 12/04/2022]
Abstract
Using molecular simulations, we study the processes of capillary condensation and capillary evaporation in model mesopores. To determine the phase transition pathway, as well as the corresponding free energy profile, we carry out enhanced sampling molecular simulations using entropy as a reaction coordinate to map the onset of order during the condensation process and of disorder during the evaporation process. The structural analysis shows the role played by intermediate states, characterized by the onset of capillary liquid bridges and bubbles. We also analyze the dependence of the free energy barrier on the pore width. Furthermore, we propose a method to build a machine learning model for the prediction of the free energy surfaces underlying capillary phase transition processes in mesopores.
Collapse
|
4
|
Desgranges C, Delhommelle J. Free energy calculations along entropic pathways. III. Nucleation of capillary bridges and bubbles. J Chem Phys 2017. [DOI: 10.1063/1.4982943] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Caroline Desgranges
- Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| | - Jerome Delhommelle
- Department of Chemistry, University of North Dakota, Grand Forks, North Dakota 58202, USA
| |
Collapse
|
5
|
Tan S(J, Do DD, Nicholson D. An efficientmethod to determine chemical potential of mixtures in the isothermal and isobaric bulk phase with kineticMonte Carlo simulation. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1090634] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
| | - D. D. Do
- School of Chemical Engineering, University of Queensland, St. Lucia, Australia
| | - D. Nicholson
- School of Chemical Engineering, University of Queensland, St. Lucia, Australia
| |
Collapse
|
6
|
Liu Z, Nguyen VT, Do D, Nicholson D. A Monte Carlo study of equilibrium transition in finite cylindrical pores. MOLECULAR SIMULATION 2013. [DOI: 10.1080/08927022.2013.829230] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
7
|
Miyahara MT, Tanaka H. Determination of phase equilibria in confined systems by open pore cell Monte Carlo method. J Chem Phys 2013; 138:084709. [PMID: 23464174 DOI: 10.1063/1.4792715] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a modification of the molecular dynamics simulation method with a unit pore cell with imaginary gas phase [M. Miyahara, T. Yoshioka, and M. Okazaki, J. Chem. Phys. 106, 8124 (1997)] designed for determination of phase equilibria in nanopores. This new method is based on a Monte Carlo technique and it combines the pore cell, opened to the imaginary gas phase (open pore cell), with a gas cell to measure the equilibrium chemical potential of the confined system. The most striking feature of our new method is that the confined system is steadily led to a thermodynamically stable state by forming concave menisci in the open pore cell. This feature of the open pore cell makes it possible to obtain the equilibrium chemical potential with only a single simulation run, unlike existing simulation methods, which need a number of additional runs. We apply the method to evaluate the equilibrium chemical potentials of confined nitrogen in carbon slit pores and silica cylindrical pores at 77 K, and show that the results are in good agreement with those obtained by two conventional thermodynamic integration methods. Moreover, we also show that the proposed method can be particularly useful for determining vapor-liquid and vapor-solid coexistence curves and the triple point of the confined system.
Collapse
Affiliation(s)
- Minoru T Miyahara
- Department of Chemical Engineering, Kyoto University, Katsura, Nishikyo, Kyoto 615-8510, Japan.
| | | |
Collapse
|
8
|
Liu Z, Do DD, Nicholson D. A thermodynamic study of the mid-density scheme to determine the equilibrium phase transition in cylindrical pores. MOLECULAR SIMULATION 2012. [DOI: 10.1080/08927022.2011.613383] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
9
|
Rasmussen CJ, Vishnyakov A, Neimark AV. Calculation of chemical potentials of chain molecules by the incremental gauge cell method. J Chem Phys 2011; 135:214109. [DOI: 10.1063/1.3657438] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
10
|
Do DD, Nicholson D, Fan C. Development of equations for differential and integral enthalpy change of adsorption for simulation studies. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2011; 27:14290-9. [PMID: 22013949 DOI: 10.1021/la203531j] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/25/2023]
Abstract
We present equations to calculate the differential and integral enthalpy changes of adsorption for their use in Monte Carlo simulation. Adsorption of a system of N molecules, subject to an external potential energy, is viewed as one of transferring these molecules from a reference gas phase (state 1) to the adsorption system (state 2) at the same temperature and equilibrium pressure (same chemical potential). The excess amount adsorbed is the difference between N and the hypothetical amount of gas occupying the accessible volume of the system at the same density as the reference gas. The enthalpy change is a state function, which is defined as the difference between the enthalpies of state 2 and state 1, and the isosteric heat is defined as the negative of the derivative of this enthalpy change with respect to the excess amount of adsorption. It is suitable to determine how the system behaves for a differential increment in the excess phase adsorbed under subcritical conditions. For supercritical conditions, use of the integral enthalpy of adsorption per particle is recommended since the isosteric heat becomes infinite at the maximum excess concentration. With these unambiguous definitions we derive equations which are applicable for a general case of adsorption and demonstrate how they can be used in a Monte Carlo simulation. We apply the new equations to argon adsorption at various temperatures on a graphite surface to illustrate the need to use the correct equation to describe isosteric heat of adsorption.
Collapse
Affiliation(s)
- D D Do
- School of Chemical Engineering, University of Queensland, St. Lucia, Queensland 4072, Australia.
| | | | | |
Collapse
|
11
|
Liu Z, Herrera L, Nguyen VT, Do DD, Nicholson D. A Monte Carlo scheme based on mid-density in a hysteresis loop to determine equilibrium phase transition. MOLECULAR SIMULATION 2011. [DOI: 10.1080/08927022.2011.578135] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
12
|
Pharmacological and toxicological target organelles and safe use of single-walled carbon nanotubes as drug carriers in treating Alzheimer disease. NANOMEDICINE-NANOTECHNOLOGY BIOLOGY AND MEDICINE 2010; 6:427-41. [DOI: 10.1016/j.nano.2009.11.007] [Citation(s) in RCA: 199] [Impact Index Per Article: 14.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/20/2009] [Revised: 11/03/2009] [Accepted: 11/18/2009] [Indexed: 01/01/2023]
|
13
|
Adsorption of light alkanes and alkenes onto single-walled carbon nanotube bundles: Langmuirian analysis and molecular simulations. Colloids Surf A Physicochem Eng Asp 2010. [DOI: 10.1016/j.colsurfa.2009.09.002] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
14
|
|