Rate equations for positronium formation at metal surfaces