Regular Article
A Gibbs free energy minimization based model for liquid–liquid equilibrium calculation of a system containing oil, brine, and surfactant
^{1}
Department of Petroleum Engineering, Faculty of Chemical Engineering, Tarbiat Modares University, P.O. Box: 141154838, Tehran, Iran
^{2}
Discipline of Chemical Engineering, School of Engineering, University of KwaZuluNatal, Howard College Campus, King George V Avenue, Durban
4041, South Africa
^{*} Corresponding author: amir_h_mohammadi@yahoo.com
Received:
23
December
2019
Accepted:
13
February
2020
Accurate and reliable phase equilibrium calculations of microemulsion systems are of great importance. This study deals with the thermodynamic modeling of Liquid–Liquid Equilibrium (LLE) of a system including oil (ndecane), brine (containing CaCl_{2} salt), and ionic surfactant (sodium dodecyl sulfonate). Two models of UNIQUAC and UNIQUAC + Debye–Hückel were used for thermodynamic calculations. The LLE experimental data were utilized to estimate the binary interaction parameters of UNIQUAC model and the adjustable parameter, b, of the Debye–Hückel model. The thermodynamic model calculates the microemulsion phase’s compositions by minimizing the Gibbs free energy of the LLE system using a combination of genetic algorithm and fmincon function in order to prevent local minima. The thermodynamic modeling results show an appropriate agreement with the experimental data. Accordingly, the presented model of this study can be used as a suitable method to investigate the liquid–liquid equilibrium of systems containing oil, water, and surfactant.
© M. Hosseini & A.H. Mohammadi, published by IFP Energies nouvelles, 2020
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Petroleum products (oil, gas, and fuels), colloidal systems, and those which result from biomass conversion are among the fluids that have received much attention in various industries [1, 2]. Because of the high demand for oil around the world, more attention is being paid to Enhanced Oil Recovery (EOR) methods [3]. Given the pressure drop in conventional reservoirs, many methods are used to recover the residual oil in reservoirs [4]. For instance, water flooding has been extensively considered as an effective method in the secondary oil recovery [5]; however, given the high capillary forces and heterogeneity of reservoir, nearly 70% original oil in place remains in the reservoir [6]. Therefore, the importance of tertiary oil recovery methods is specified. In these methods, different chemicals including surfactants, polymers, and so on along with seawater are injected to oil reservoirs [7–9]. To this end, different surfactants, including cationic, anionic or even a mixture of these two, are widely used [10–13]. As a result of the surfactants injection, the surfactant is initially located due to its structure between the oil and water phases, which leads to the formation of a threephase system including aqueous brine, oil as well as, surfactant, and also reduces the interfacial tension [14–17]. Increasing the surfactants’ concentration results in the formation of a continuous phase, or the microemulsion. Finally, a twophase system consisting of oil and microemulsion phases, which are in equilibrium, are formed [17, 18].
According to the previous researches, one can conclude that surfactant plays a very important role in the LLE of oil, brine, and surfactant systems. Therefore, accurate and reliable phase equilibrium calculations of microemulsion systems to facilitate the EOR conditions are of great importance so that the exact number of phases and the values of the components are determined in equilibrium [19–22]. As indicated in the other researches, there are two methods for calculating liquid–liquid equilibrium: in the first method, the isoactivity equations are solved under mass balance constraints [23, 24], and in the second method, the Gibbs free energy minimization is addressed [25]. In the second method, the results are dependent on initial estimates, and poor estimates may lead to wrong solutions [26]. For instance, in 1999, Yushan and Zhihong proposed an algorithm for liquid–liquid phase equilibrium calculation. They utilized tangent plane distance function criterion for stability analysis so that this approach was considered as sufficient and necessary condition. Moreover, activity coefficient models, namely UNIQUAC and NRTL were used to formulate Gibbs free energy functions [27]. In 2018, Dadmohammadi et al. employed the NRTL, mNRTL1, and mNRTL2 activity coefficient models for the characterization of LLE systems [28]. Their study showed the efficiency of the modified NRTL models in characterization of LLE systems. Also in 2016, Li et al. unlike the use of minimizing the Gibbs free energy (coupled with stability test), solved isoactivity equations for their LLE calculations using NRTL model, and conducted their study on twophase LLE systems. They concluded that this approach could be applied to multiphase systems [26]. It should be noted that both Equation of State (EOS) and activity coefficient models have been implemented in different researches for characterization of the phase behavior of LLE [29–35] and vapor–liquid equilibrium [36] systems. Riazi and Moshfeghian presented a thermodynamic approach using NRTL activity coefficient model. They used Debye–Hückel theory to calculate the longrange interactions. They investigated the LLE of a system containing water, oil, and surfactant. For the adjustable parameter of the NRTL model, the value of 0.1 was chosen. In addition, only two interaction parameters were optimized and others considered equal to zero [37].
In this study, the LLE of oil, brine, and surfactant systems was studied through Gibbs free energy minimization method on the data provided in Riazi and Moshfeghian’s research [37]. To do so, two models of UNIQUAC and UNIQUAC + Debye–Hückel were used, and consequently the results were compared with the previous research.
2 Experimental data
As mentioned earlier, the experimental (phase equilibrium) data of a system including oil (ndecane), brine (containing CaCl_{2} salt), and an ionic surfactant, named sodium dodecyl sulfonate, were obtained from Riazi and Moshfeghian’s research [37]. The studied system includes two phases of oil and microemulsion, so that in each of the two phases, all three components of oil, brine, and surfactant are available. It should be noted that the salt presents in water is CaCl_{2} with a concentration of 500 mg/L (Ca^{++}), which is equal to 0.138 wt%. In addition, the experiments were carried out at a temperature of 20 °C and atmospheric pressure. Three concentrations of surfactants with values of 0.5 wt%, 1 wt%, and 2 wt% were used in the experiments, and in each set of experiments ten different oil volumes ranging from 5% to 90% of the total volume of system were utilized. The initial mixture compositions for three sets of experiments as well as the physical properties and mole fractions of the components in microemulsion phase are presented in Tables 1–3, respectively. It is noteworthy that the ϕ parameter represents the molar ratio of the oil phase to the total number of mole in the system, the amount of which alters from 0 to 1. It should be noted that numbers 1–3 are considered for oil, water, and surfactant molecules, respectively.
The structural parameters for the functional groups are reported in Table 4 [38, 39].
3 Thermodynamic model
3.1 The UNIQUAC activity coefficient model
The molar excess Gibbs and activity coefficient functions of UNIQUAC model are as follows [40]:(1) (2) (3) (4) (5) (6) (7) (8) (9)where r, q, and q′ indicate molecular structure parameters of a pure component, expressing the volume and surface areas of molecule, respectively, where with the exception of water and some small alcohols, q = q′. Also, the subscript k and superscript i in equations (7) and (8) refer to the group k and the occurrence of the group k in molecule i, respectively. The parameter Z, which implies the coordination number, is assumed as a fixed value of 10. The variable Φ* represents the segment fraction. In addition, θ and θ′ denote the area fractions of a simple and a hydrogen bonded molecule, respectively.
3.2 Debye–Hückel theory
The ionic activity coefficient for electrolyte is as follows [41]:(10)where Z _{+} and Z _{−} are the positive and negative charges of the electrolyte. The parameter b is an adjustable parameter, and I is the ionic strength which is defined as follows:(11)where m _{ i } and Z _{ i } stand for the molality of ion i (mol/kg of water) and the valances of ions, respectively.
In Debye–Hückel theory, the A _{ γ } coefficient is calculated as follows:(12)
In the calculations, it is assumed that the surfactant is fully dissociated into positive and negative ions just in the microemulsion phase, and the dissociation does not occur in the oil phase. The experimental data of this study were measured at 20 °C [37]. The constant values utilized in equation (12) are presented in Table 5. The calculated value for A _{ γ } is equal to 1.1667.
The structural parameters for oil, water, and surfactant molecules in UNIQUAC activity coefficient model.
The regressed binary parameter for UNIQUAC model using binary LLE data.
3.3 Gibbs free energy minimization
Calculations were performed based on the Gibbs free energy minimization [19]. Total Gibbs free energy of a mixture of n _{ c } components and n _{ p } phases is obtained from the following equations:(13)where,(14) (15)where n ^{ k }, , F, and z _{ i } are the total number of moles in phase k, the number of moles of component i phase k, the total number of moles in input feed, and the molar fraction of component i in the input feed, respectively. In the case of using the Gibbs energy functions for the LLE calculations, the following equation is utilized to calculate the Gibbs energy of the mixture:(16)where the superscript, L(k) denotes to liquid phase k, and γ _{ i } represents the activity coefficient. Through the minimization of the above equation, equilibrium concentrations are calculated in each phase. In the above equation, the mole number of each component in each phase is unknown that can be calculated by minimizing. To achieve this goal, change of variables, (for i = 1, 2, …, n _{ c }; k = 1, 2, …, n _{ p }−1), is used which is as follows:(17) (18) (19)
The unknown variables are , which are calculated by minimizing. The following limitation must be applied to avoid negative or imaginary answers:(20)
3.4 Liquid–liquid equilibrium calculation
In this study, to find the accurate values of in the minimization of Gibbs energy and the binary interaction parameters of UNIQUAC and UNIQUAC + Debye–Hückel models, a combination of genetic algorithm and fmincon function (fmincon function within MATLAB environment) was utilized which leads to calculation of logical answers in all points.
One of the methods to solve constrained and unconstrained optimization problems is to use genetic algorithm based on natural selection, a process that leads to biological evolution. The way it works is that the genetic algorithm frequently modifies a population of individual solutions. At each stage, the genetic algorithm randomly chooses individuals from the current population as parents and also it utilizes them to produce children for the next generation. Consequently, throughout consecutive generations, the population “evolves” toward an optimal solution. Genetic algorithms can be applied to solve different optimization problems that cannot be implemented using standard optimization algorithms. These include problems where the objective function is stochastic, discontinuous, nonlinear, or nondifferentiable [42, 43]. In addition, the fmincon function is used to find a constrained minimum of a scalar function whose variables start from an initial estimation.
As we know, in the calculations of the liquid–liquid equilibrium, the equilibrium condition is as follows:(21)where O and ME denote to the oil and microemulsion phases, respectively. In equation (21), the activity coefficients are calculated based on the mole fraction. The problem with this is that the Debye–Hückel model is based on molality. Therefore, the conversion of the unit of activity coefficient from molality to mole fraction is required. To do so, the equation provided by Panah et al. can be used [44], which is as follows:(22)where x, γ ^{(x)}, m, γ ^{(m)} are the mole fraction, mole fraction based activity coefficient, molality, and molality based activity coefficient, respectively. In the above equation, the term of “1000/18.01528” is the mole number of water contained in 1 kg.
When using the UNIQUAC + Debye–Hückel model for the surfactant, the UNIQUAC model is used to calculate the activity in the oil phase and the Debye–Hückel model is applied for the microemulsion phase. Therefore, the normalization of the Debye–Hückel model is necessary and it is proposed as follows:(23)where and are the calculated activity coefficient of the surfactant by the Debye–Hückel model and the surfactant activity coefficient at infinite dilution, respectively. In the calculations, the value of was obtained by fitting which is reported in Table 8.
The regressed binary parameter for UNIQUAC + Debye–Hückel model using binary LLE data.
To find the binary interaction parameters of UNIQUAC activity coefficient model, the adjustable parameter of Debye–Hückel model (b), as well as the surfactant activity coefficient at infinite dilution (), equation (21) should be solved. Afterwards the Gibbs free energy minimization will be performed.
In UNIQUAC + Debye–Hückel model, the activity coefficient of ndecane and water in both phases are calculated using UNIQUAC equations. The activity coefficient of the surfactant in the oil phase is also calculated using UNIQUAC equations, but the activity coefficient in the microemulsion phase is calculated using Debye–Hückel equations.
In the last step of the present study, Absolute Average Deviation (AAD) is used to find the quality of the predicted data compared to the experimental data. The AAD% is defined as below:(24)where N is the number of tie lines and X is the variable whose error is calculated. In addition, the superscripts “exp” and “cal” refer to experimental and calculated data, respectively.
4 Results and discussion
The structural parameters for oil, water, and surfactant molecules in UNIQUAC activity coefficient model are calculated using equations (7) and (8). These parameters are presented in Table 6.
As Teh and Rangaiah mentioned in their study, in order to have precise phase equilibrium calculations, an accurate global minimum is required [19]. As we know, although the genetic algorithm optimizes the parameters globally, this optimization is not accurate, since different initial guesses for parameters give different answers. Moreover, functions such as fmincon optimize the parameters locally. To solve these problems, in the present study, the combination of these two algorithms was used so that both binary interaction and parameters are first optimized by genetic algorithm and then by fmincon function, this procedure leads to accurate results.
The binary interaction parameters for UNIQUAC and UNIQUAC + Debye–Hückel models are presented in Tables 7 and 8, respectively.
The UNIQUAC and UNIQUAC + Debye–Hückel models have 6 and 7 binary interaction parameters, respectively. In UNIQUAC + Debye–Hückel model, full dissociation of the surfactant into two positive and negative ions in the microemulsion phase is assumed, but in UNIQUAC model, no dissociation is considered.
The calculated values of the compositions (%) in the microemulsion phase as well as ϕ values (%) for UNIQUAC and UNIQUAC + Debye–Hückel models are presented in Table 9. In addition, for a better comparison, for example the predicted values of the present study for surfactant composition as well as ϕ values are shown in Figures 1 and 2, respectively, as compared to the experimental data and the model presented by Riazi and Moshfeghian [37].
Fig. 1 Prediction of surfactant composition in the microemulsion phase using UNIQUAC and UNIQUAC + Debye–Hückel in comparison with the experimental data and the presented NRTL model by Riazi and Moshfeghian [37] for oil/brine/surfactant system. 
Fig. 2 Prediction of Φ values using UNIQUAC and UNIQUAC + Debye–Hückel in comparison with the experimental data and the presented NRTL model by Riazi and Moshfeghian [37] for oil/brine/surfactant system. 
Calculated molar composition (%) of microemulsion phase using UNIQUAC and UNIQUAC + Debye–Hückel models.
It should be mentioned that the mixture of water and CaCl_{2} is not considered as a pseudo component. The concentration of salt in the initial mixture is very low and its presence is ignored in the oil phase.
The AAD% of the present modeling results are presented in Table 10 in comparison with the results acquired by Riazi and Moshfeghian [37].
The AAD% values.
Due to the low concentrations of surfactants, the addition of electrolyte term or, in other words, the consideration of longrange interactions, has little effect on the improvement of results. Previously, this conclusion was reported by Riazi and Moshfeghian [37]. Moreover, the comparison between the results obtained by Riazi and Moshfeghian [37] with the results of this study shows a much better performance of the UNIQUAC and UNIQUAC + Debye–Hückel models, especially in predicting the oil composition in the microemulsion phase.
It is necessary to mention that all three models have poor performances in predicting the molar fraction of the oil component. However, UNIQUAC and UNIQUAC + Debye–Hückel models predict the molar fraction of the oil component with a higher accuracy than the NRTL model used by Riazi and Moshfeghian [37]. Nevertheless, in the analysis of the obtained results, this point should be considered that the NRTL model requires four fitting parameters to describe the proposed LLE system; while the obtained lower errors by the UNIQUAC and UNIQUAC + Debye–Hückel models were originated from more fitting values and it is logical to present better performances than the NRTL model. Regarding the molar fraction of the water component, all three models give the same results. The NRTL model employed by Riazi and Moshfeghian [37] offers a more successful performance in predicting the molar fraction of the surfactant component, and there is an appropriate agreement with the experimental values. The UNIQUAC and UNIQUAC + Debye–Hückel models represent an excellent performance in predicting Φ. However, in the NRTL model used by Riazi and Moshfeghian [37], a very high error is observed.
One of the goals pursued in modeling studies could be to reduce the fitting parameters to the lowest possible number as the high number of fitting parameters can reduce the validity of the model. In 2015, Dadfar and Biria studied the LLE of the system with the same experimental data used in the present study [21]. They used 35 fitting parameters. However, the error obtained by them were not considerably better.
It should be mentioned that, as the oil composition is not reported in the literature [37], the oil is considered as ndecane. Moreover, in modeling study carried out on the same system, it is common to consider oil as ndecane [21].
5 Conclusion
Gibbs free energy minimization based modeling of the LLE conditions for a system containing oil, brine, and surfactant was investigated in this work. This modeling approach employs UNIQUAC and UNIQUAC + Debye–Hückel activity coefficients models. A combination of genetic algorithm and fmincon function was used which leads to calculation of reasonably accurate results. According to Table 10, one can note that the investigated model in this study displays an acceptable overall performance in predicting the molar fractions of components as compared to a previous study reported in literature so that the experimental and computational values were slightly different. It is concluded that the addition of electrolyte term has little effect on the improvement of results given the low concentrations of surfactants. In fact, considering the electrolyte term, the AAD% value for the mole fractions of the components and ϕ changes less than 1%. The model can be used as a robust method to investigate the liquid–liquid equilibrium of systems containing oil, water, and surfactant.
References
 Mahi M.R., Ouaar F., Negadi A., Bahadur I., Negadi L. (2018) Excess/deviation properties of binary mixtures of 2,5dimethylfuran with furfuryl alcohol, methyl isobutyl ketone, 1butanol and 2butanol at temperature range of (293.15–323.15) K, Oil Gas Sci. Technol.  Rev. IFP Energies nouvelles 73, 64. [CrossRef] [Google Scholar]
 Jafari A., Hasani M., Hosseini M., Gharibshahi R. (2019) Application of CFD technique to simulate enhanced oil recovery processes: Current status and future opportunities, Pet. Sci. 1–23. https://doi.org/10.1007/s12182019003637 [Google Scholar]
 Hirasaki G.J., Miller C.A., Puerto M. (2008) Recent advances in surfactant EOR, in: IPTC 2008: International Petroleum Technology Conference, 3–5 December, Kuala Lumpur, Malaysia. [Google Scholar]
 Chen X., Feng Q., Liu W., Sepehrnoori K. (2017) Modeling preformed particle gel surfactant combined flooding for enhanced oil recovery after polymer flooding, Fuel 194, 42–49. [CrossRef] [Google Scholar]
 Craig F.F. (1971) The reservoir engineering aspects of waterflooding, Vol. 3, HL Doherty Memorial Fund of AIME, New York. [Google Scholar]
 Kumar N., Mandal A. (2018) Surfactant stabilized oilinwater nanoemulsion: Stability, interfacial tension, and rheology study for enhanced oil recovery application, Energy Fuels 32, 6452–6466. [Google Scholar]
 Shi S., Wang Y., Wang L., Jin Y., Wang T., Wang J. (2015) Potential of spontaneous emulsification flooding for enhancing oil recovery in hightemperature and highsalinity oil reservoir, J. Dispers. Sci. Technol. 36, 660–669. [Google Scholar]
 Liu S., Zhang D., Yan W., Puerto M., Hirasaki G.J., Miller C.A. (2008) Favorable attributes of alkalinesurfactantpolymer flooding, SPE J. 13, 5–16. [CrossRef] [Google Scholar]
 Manrique E.J., Thomas C.P., Ravikiran R., Izadi Kamouei M., Lantz M., Romero J.L., Alvarado V. (2010) EOR: current status and opportunities, in:SPE Improved Oil Recovery Symposium, 24–28 April, Tulsa, OK. Society of Petroleum Engineers. [Google Scholar]
 Puerto M.C. (2001) Surfactants: Fundamentals and Applications in the Petroleum IndustryCambridge University Press, 2000, pp. 621, @ 85.00 (US 140.00)(hardback), ISBN 0521640679, Chem. Eng. J. 83, 1, 63. [Google Scholar]
 Li Y., Zhang W., Kong B., Puerto M., Bao X., Sha O., Shen Z., Yang Y., Liu Y., Gu S. (2016) Mixtures of anionic/cationic surfactants: A new approach for enhanced oil recovery in lowsalinity, hightemperature sandstone reservoir, SPE J. 21, 1–164. [CrossRef] [Google Scholar]
 Ahmadi S., Hosseini M., Tangestani E., Mousavi S.E., Niazi M. (2020) Wettability alteration and oil recovery by spontaneous imbibition of smart water and surfactants into carbonates, Pet. Sci. 1–10. https://doi.org/10.1007/s12182019004121. [Google Scholar]
 Dehaghani A.H.S., Badizad M.H. (2019) Impact of ionic composition on modulating wetting preference of calcite surface: Implication for chemically tuned water flooding, Colloids Surf. A Physicochem. Eng. Asp. 568, 470–480. [Google Scholar]
 Sheng J.J. (2010) Modern chemical enhanced oil recovery: theory and practice, Gulf Professional Publishing, Oxford, UK. [Google Scholar]
 Flaaten A., Nguyen Q.P., Pope G.A., Zhang J. (2008) A systematic laboratory approach to lowcost, highperformance chemical flooding, in: SPE Symposium on Improved Oil Recovery, 20–23 April, Tulsa, OK. Society of Petroleum Engineers. [Google Scholar]
 Hirasaki G., Zhang D.L. (2004) Surface chemistry of oil recovery from fractured, oilwet, carbonate formations, SPE J. 9, 151–162. [CrossRef] [Google Scholar]
 Rosen M.J. (2004) Emulsification by surfactants, in: Surfactants Interfacial Phenomena, 3rd edn., Wiley, Hoboken, NJ, pp. 303–331. [CrossRef] [Google Scholar]
 Winsor P.A. (1948) Hydrotropy, solubilisation and related emulsification processes, Trans. Faraday Soc. 44, 376–398. [CrossRef] [Google Scholar]
 Teh Y.S., Rangaiah G.P. (2002) A study of equationsolving and Gibbs free energy minimization methods for phase equilibrium calculations, Chem. Eng. Res. Des. 80, 745–759. [Google Scholar]
 Ashrafizadeh S.N., Motaee E., Hoshyargar V. (2012) Emulsification of heavy crude oil in water by natural surfactants, J. Pet. Sci. Eng. 86, 137–143. [Google Scholar]
 Dadfar B., Biria D. (2015) Application of group contribution–NRTL model with closure to predict LLE behavior of an oil/brine/surfactant system, J. Chem. Eng. Data 60, 2575–2584. [Google Scholar]
 Jin L., Budhathoki M., Jamili A., Li Z., Luo H., Ben Shiau B.J., Delshad M., Harwell J.H. (2017) Predicting microemulsion phase behavior using physics based HLDNAC equation of state for surfactant flooding, J. Pet. Sci. Eng. 151, 213–223. [Google Scholar]
 Privat R., Jaubert J.N., Privat Y. (2013) A simple and unified algorithm to solve fluid phase equilibria using either the gamma–phi or the phi–phi approach for binary and ternary mixtures, Comput. Chem. Eng. 50, 139–151. [Google Scholar]
 Panah H.S. (2018) Modeling binary vapor–liquid equilibrium data containing perfluorocarbons using the PengRobinson and the PCSAFT equations of state, Int. J. Refrig. 85, 13–26. [Google Scholar]
 Li Z., Mumford K.A., Shang Y., Smith K.H., Chen J., Wang Y., Stevens G.W. (2014) Analysis of the nonrandom twoliquid model for prediction of liquid–liquid equilibria, J. Chem. Eng. Data 59, 2485–2489. [Google Scholar]
 Li Z., Mumford K.A., Smith K.H., Chen J., Wang Y., Stevens G.W. (2016) Solution structure of isoactivity equations for liquidliquid equilibrium calculations using the nonrandom twoliquid model, Ind. Eng. Chem. Res. 55, 2852–2859. [Google Scholar]
 Yushan Z.H.U., Zhihong X.U. (1999) A reliable method for liquidliquid phase equilibrium calculation and global stability analysis, Chem. Eng. Commun. 176, 133–160. [Google Scholar]
 Dadmohammadi Y., Gebreyohannes S., Neely B.J., Gasem K.A.M. (2018) Application of modified NRTL models for binary LLE phase characterization, Ind. Eng. Chem. Res. 57, 7282–7290. [Google Scholar]
 Boshkov L.Z., Yelash L.V. (1997) Closedloops of liquidliquid immiscibility in binary mixtures predicted from the RedlichKwong equation of state, Fluid Phase Equilib. 141, 105–112. [Google Scholar]
 Sofyan Y., Ghajar A.J., Gasem K.A.M. (2003) A systematic method to predict cloud point temperature and solid precipitation, Pet. Sci. Technol. 21, 409–424. [Google Scholar]
 Domańska U., Morawski P., Wierzbicki R. (2006) Phase diagrams of binary systems containing nalkanes, or cyclohexane, or 1alkanols and 2, 3pentanedione at atmospheric and high pressure, Fluid Phase Equilib. 242, 154–163. [Google Scholar]
 Giovanoglou A., Adjiman C.S., Jackson G., Galindo A. (2009) Fluid phase stability and equilibrium calculations in binary mixtures: Part II: Application to singlepoint calculations and the construction of phase diagrams, Fluid Phase Equilib. 275, 95–104. [Google Scholar]
 Królikowska M., Karpińska M. (2013) Phase equilibria study of the (Noctylisoquinolinium thiocyanate ionic liquid + aliphatic and aromatic hydrocarbon, or thiophene) binary systems, J. Chem. Thermodyn. 63, 128–134. [Google Scholar]
 Królikowska M., Karpińska M., Zawadzki M. (2013) Phase equilibria study of (ionic liquid + water) binary mixtures, Fluid Phase Equilib. 354, 66–74. [Google Scholar]
 ReynelÁvila H.E., BonillaPetriciolet A., TapiaPicazo J.C. (2019) An artificial neural networkbased NRTL model for simulating liquidliquid equilibria of systems present in biofuels production, Fluid Phase Equilib. 483, 153–164. [Google Scholar]
 Nasrifar K., Rahmanian N. (2018) Equations of state with group contribution binary interaction parameters for calculation of twophase envelopes for synthetic and real natural gas mixtures with heavy fractions, Oil Gas Sci. Technol.  Rev. IFP Energies nouvelles 73, 7. [CrossRef] [Google Scholar]
 Riazi M.R., Moshfeghian M. (2009) A thermodynamic model for LLE behavior of oil/brine/ionicsurfactant/alcohol cosurfactant systems for EOR processes, J. Pet. Sci. Eng. 67, 75–83. [Google Scholar]
 Cheng H., Kontogeorgis G.M., Stenby E.H. (2005) Correlation and prediction of environmental properties of alcohol ethoxylate surfactants using the UNIFAC method, Ind. Eng. Chem. Res. 44, 7255–7261. [Google Scholar]
 Qiu T., Li S., Li S., Wu Y. (2009) Liquid–liquid phase equilibria of the ternary system of water/1, 4dioxane/dihydromyrcene, Fluid Phase Equilib. 280, 84–87. [Google Scholar]
 Prausnitz J.M., Lichtenthaler R.N., De Azevedo E.G. (1998) Molecular thermodynamics of fluidphase equilibria, Pearson Education, London, UK. [Google Scholar]
 Kontogeorgis G.M., Folas G.K. (2009) Thermodynamic models for industrial applications: From classical and advanced mixing rules to association theories, John Wiley & Sons, Hoboken, NJ. [Google Scholar]
 Sivanandam S.N., Deepa S.N. (2008) Genetic algorithms, in:Introduction to genetic algorithms, Springer, Berlin, Germany, pp. 15–37. [CrossRef] [Google Scholar]
 Gen M., Cheng R., Lin L. (2008) Network models and optimization: Multiobjective genetic algorithm approach, Springer Science & Business Media, Berlin, Germany. [Google Scholar]
 Panah H.S., Mohammadi A.H., Ramjugernath D. (2016) Development of a novel approach for modeling acid gas solubility in alkanolamine aqueous solution, J. Nat. Gas Sci. Eng. 34, 112–123. [Google Scholar]
All Tables
The structural parameters for oil, water, and surfactant molecules in UNIQUAC activity coefficient model.
The regressed binary parameter for UNIQUAC + Debye–Hückel model using binary LLE data.
Calculated molar composition (%) of microemulsion phase using UNIQUAC and UNIQUAC + Debye–Hückel models.
All Figures
Fig. 1 Prediction of surfactant composition in the microemulsion phase using UNIQUAC and UNIQUAC + Debye–Hückel in comparison with the experimental data and the presented NRTL model by Riazi and Moshfeghian [37] for oil/brine/surfactant system. 

In the text 
Fig. 2 Prediction of Φ values using UNIQUAC and UNIQUAC + Debye–Hückel in comparison with the experimental data and the presented NRTL model by Riazi and Moshfeghian [37] for oil/brine/surfactant system. 

In the text 