The cubic-plus-association equation of state was applied to model the phase behaviour of bitumen/nalkane systems including saturation pressures, liquid-liquid boundaries, yields, and phase compositions. Yield is defined here as the mass of bitumen in the heavy phase divided by the mass of bitumen in the feed. To implement the model, the bitumen was divided into a set of pseudo-components based on a distillation assay and either the n-pentane insoluble content of the oil (CPA-C5 approach) or the propane insoluble content (CPA-C3 approach). The pseudo-components in the solvent insoluble part of the oil were defined as self-associating components, all other pseudo-components were non-associating. The critical properties and acentric factor for each pseudo-component were determined from established correlations. A set of CPA parameters was then developed to fit the available phase behavior data. The self-associating pseudo-components were assigned a distribution of self-association energies in order to capture the sequential partitioning of asphaltenes to the heavy phase upon solvent addition or in different solvents. Both approaches matched the phase behavior data for mixtures of bitumen with n-pentane and higher carbon number n-alkanes almost to within the experimental error. The CPA-C3 approach also matched the phase behavior data for mixtures of propane and bitumen. The CPA-C5 approach could not match the yield data for propane diluted bitumen but was more straightforward to implement and was less computationally intensive because it employed fewer self-associating components. To apply either approach to another oil, only the self-association energy of the self-associating pseudo-components need be adjusted. The cross-association energy between the solvent and the self-associating pseudo-components must be tuned for any new solvent. (C) 2019 Elsevier B.V. All rights reserved.