The energy hub concept is gradually getting popular due to its capability to integrate renewable energy technologies into the energy grid with a minimum impact. Designing distributed energy hubs is a challenging task due to the coupling of optimum dispatch and energy system sizing problems. Operation of the system needs to be considered for 8760 time steps and energy system sizing optimization can take several days to complete. This time must be shortened in order to make it easier to optimize the multi-energy hubs connected to multi energy grids where there is a strong coupling among energy hubs and network. This study introduces a novel optimization algorithm, coupling an existing energy hub model with a hybrid surrogate model in order to reduce computational time in the optimization process.