Modelling the asymmetry of the halo cross-correlation function with relativistic effects at quasi-linear scales