Modelling CO2-Brine Interfacial Tension using Density Gradient Theory

Abstract

Knowledge regarding carbon dioxide (CO2)-brine interfacial tension (IFT) is important for petroleum industry and Carbon Capture and Storage (CCS) strategies. In petroleum industry, CO2-brine IFT is especially importance for CO2 – based enhanced oil recovery strategy as it affects phase behavior and fluid transport in porous media. CCS which involves storing CO2 in geological storage sites also requires understanding regarding CO2-brine IFT as this parameter affects CO2 quantity that could be securely stored in the storage site. Several methods have been used to compute CO2-brine interfacial tension. One of the methods employed is by using Density Gradient Theory (DGT) approach. In DGT model, IFT is computed based on the component density distribution across the interface. However, current model is only applicable for modelling low to medium ionic strength solution. This limitation is due to the model only considers the increase of IFT due to the changes of bulk phases properties and does not account for ion distribution at interface.

In this study, a new modelling strategy to compute CO2-brine IFT based on DGT was proposed. In the proposed model, ion distribution across interface was accounted for by separating the interface to two sections. The saddle point of tangent plane distance where ( ) was defined as the boundary separating the two sections of the interface. Electrolyte is assumed to be present only in the second section which is connected to the bulk liquid phase side. Numerical simulations were performed using the proposed approach for single and mixed salt solutions for three salts (NaCl, KCl, and CaCl2), for temperature (298 K to 443 K), pressure (2 MPa to 70 MPa), and ionic strength (0.085 mol·kg-1 to 15 mol·kg-1). The simulation result shows that the tuned model was able to predict with good accuracy CO2-brine IFT for all studied cases. Comparison with current DGT model showed that the proposed approach yields better match with the experiment data. In this study, the thermodynamic properties were computed using Cubic Plus Association (CPA) equation of state, and the electrolyte contribution was accounted for by adding Debye-Huckel activity coefficient in the thermodynamic properties computation.