Quasi-ray Gaussian beam algorithm for time-harmonic two-dimensional scattering by moderately rough interfaces