The present study aims to implement an electrolyte EoS developed in our previous work [A. Raeispour Shirazi, M. N. Lotfollahi, Fluid Phase Equilib., 502 (2019) 112289] to model the vapor-liquid equilibria (VLE) of CO2 solubility in aqueous N-methyldiethanolamine (MDEA) solution (CO2 +MDEA-H2O ternary mixture) as well as H2O+MDEA+CO2-CH4 quaternary mixture. The proposed electrolyte PC-SAFT+MSA+Born equation of state ((ePC_SAFT-MB) consists of PC-SAFT equation of state, the mean spherical approximation (MSA) term and the Born term. First, an appropriate association scheme for water is specified via performing VLE calculations for H2O-CO2 binary mixture and without any adjustment in binary interaction parameter (BIP). By determining the proper association scheme for water (4C), the BIPs of H2O-CO2 , MDEA-CH4 and CH4-H2O binary mixtures are adjusted with binary VLE experimental data. In what follows, the nonelectrolyte PC-SAFT EoS is applied considering the obtained BIPs to predict 162 selected data points of bubble pressure for H2O-MDEA-CO2 ternary system in a wide range of temperature (298-413 K), MDEA concentration (18.7-48.8 wt.%) and CO2 loading (0.0001-1.5). The average absolute relative deviation (AAD%) for the calculation results is obtained 25.9%. In order to consider the ions interactions in the solution and to improve the precision of calculations, Born and MSA contributions are added to PC-SAFT EoS. By adjusting two BIPs (k(MDEAH+ --CO2) and k(MDEAH+) (-HCO3-)) and without any readjustment of other BIPs, the AAD% is obtained 15.3%. Finally, the obtained BIPs between the molecular and ionic species for proposed ePC_SAFT-MB EoS are employed to predict the VLE of H2O-MDEA-CO2-CH4 quaternary mixture without any readjustment in BIPs. For this case, the AAD% of 28.05% indicates the reliability of proposed ePC_SAFT-MB EoS for the prediction of H2O-MDEA-CO2-CH4 quaternary VLE results. (C) 2020 Elsevier B.V. All rights reserved.