Multiphase equilibrium calculations are important in primary and enhanced oil recovery production mechanisms such as natural depletion, gas injection, and steam injection. Predicting the physical properties of phases under different conditions is necessary for designing optimum recovery techniques. In this work, we report a cubic-plus-association equation of state (CPA EoS) parameterization of vapor−liquid−liquid equilibrium (VLLE) of water/bitumen/solvent (C 1 , C 2 , C 3 , and n-C 4 ) systems. In this approach, first, the unknown parameters of the EoS are tuned versus the solvent mole fraction in the bitumen-rich phase and saturated bitumen density data. Next, the binary interaction parameters (BIPs) between bitumen and hydrocarbons (C 1 , C 2 , C 3 , and n-C 4 ) and nonhydrocarbon (CO 2 and N 2 ) solvents are adjusted versus experimental data. The absolute average relative deviation (AARD) values for solubility of