Robust optimization of energy hubs operation based on extended affine arithmetic