1. A method for modeling shallow water multiples, the method comprising: reconstructing a seafloor reflectivity image using seismic data; evaluating a seismic wave equation using, as an input, said reconstructed seafloor reflectivity image to determine surface multiples; and subtracting said determined surface multiples from said seismic data.