Representation of bi-parameter singular integrals by dyadic operators