We present calculations on energy- and time-resolved two-photon photoemission spectra of images states in Cu(100) and Cu(111) surfaces. The surface is modeled by a 1D effective potential and the states are propagated within a real-space, real-time method. To obtain the energy resolved spectra we employ a geometrical approach based on a subdivision of space into two regions. We treat electronic inelastic effects by taking into account the scattering rates calculated within a GW scheme. To get further insight into the decaying mechanism we have also studied the effect of the variation of the classical Hartree potential during the excitation. This effect turns out to be small.