Interfacial Tension–Temperature–Pressure–Salinity Relationship for the Hydrogen–Brine System under Reservoir Conditions: Integration of Molecular Dynamics and Machine Learning

Sina Omrani,Mehdi Ghasemi,Mrityunjay Singh,Saeed Mahmoodpour,Tianhang Zhou,Masoud Babaei,Vahid Niasar
DOI: https://doi.org/10.1021/acs.langmuir.3c01424
IF: 3.9
2023-08-31
Langmuir
Abstract:Hydrogen (H<sub>2</sub>) underground storage has attracted considerable attention as a potentially efficient strategy for the large-scale storage of H<sub>2</sub>. Nevertheless, successful execution and long-term storage and withdrawal of H<sub>2</sub> necessitate a thorough understanding of the physical and chemical properties of H<sub>2</sub> in contact with the resident fluids. As capillary forces control H<sub>2</sub> migration and trapping in a subsurface environment, quantifying the interfacial tension (IFT) between H<sub>2</sub> and the resident fluids in the subsurface is important. In this study, molecular dynamics (MD) simulation was employed to develop a data set for the IFT of H<sub>2</sub>-brine systems under a wide range of thermodynamic conditions (298-373 K temperatures and 1-30 MPa pressures) and NaCl salinities (0-5.02 mol·kg<sup>-1</sup>). For the first time to our knowledge, a comprehensive assessment was carried out to introduce the most accurate force field combination for H<sub>2</sub>-brine systems in predicting interfacial properties with an absolute relative deviation (ARD) of less than 3% compared with the experimental data. In addition, the effect of the cation type was investigated for brines containing NaCl, KCl, CaCl<sub>2</sub>, and MgCl<sub>2</sub>. Our results show that H<sub>2</sub>-brine IFT decreases with increasing temperature under any pressure condition, while higher NaCl salinity increases the IFT. A slight decrease in IFT occurs when the pressure increases. Under the impact of cation type, Ca<sup>2+</sup> can increase IFT values more than others, i.e., up to 12% with respect to KCl. In the last step, the predicted IFT data set was used to provide a reliable correlation using machine learning (ML). Three white-box ML approaches of the group method of data handling (GMDH), gene expression programming (GEP), and genetic programming (GP) were applied. GP demonstrates the most accurate correlation with a coefficient of determination (<i>R</i><sup>2</sup>) and absolute average relative deviation (AARD) of 0.9783 and 0.9767%, respectively.
chemistry, multidisciplinary, physical,materials science
What problem does this paper attempt to address?