Surrogate model for gravitational waveforms of spin-aligned binary black holes with eccentricities