In this work, we propose a novel methodology to adjust parameters of compartmental epidemiological models. It is based on solving a multi-objective optimization problem that consists in fitting some of the model outputs to real observations. First, according to the available data of the considered epidemic, we define a multi-objective optimization problem where the model parameters are the optimization variables. Then, this problem is solved by considering a particular optimization algorithm called ParWASF-GA (ParallelWeighting Achievement Scalarizing Function Genetic Algorithm). Finally, the decision maker chooses, within the set of possible solutions, the values of parameters that better suit his/her preferences. In order to illustrate the benefit of using our approach, it is applied to estimate the parameters of a deterministic epidemiological model, called Be-CoDiS (Between-Countries Disease Spread), used to forecast the possible spread of human diseases within and between countries. We consider data from different Ebola outbreaks from 2014 up to 2019. In all cases, the proposed methodology helps to obtain reasonable predictions of the epidemic magnitudes with the considered model.