A computer-implemented method comprising: receiving seismic data for a geologic region of the Earth, wherein the seismic data comprise primaries and multiples acquired using a seismic survey; building a velocity model that characterizes the geologic region of the Earth pre-processing the received seismic data by applying one or more multiple attenuation techniques to preserve primaries in the seismic data; generating a first image representative of coverage and quality of the primaries based on the pre-processed data; identifying a portion of the geologic region of the earth in the first image that comprises a deficiency in coverage, wherein the deficiency in coverage corresponds to one or more gaps between adjacent sail lines of the seismic survey in a shallow subsurface region; pre-processing the received seismic data a second time by applying one or more primary attenuation techniques to preserve one or more modes of multiples in the seismic data; selecting at least one mode of multiple and corresponding travel time data from a data storage , wherein the travel time data correspond to at least one complex ray signature in the geologic region of the earth and are based at least in part on the velocity model; performing Kirchhoff depth migration on the pre-processed seismic data with preserved modes of multiples, using at least the selected travel time data to generate migrated seismic data; and rendering a second image of the geologic region of the Earth to a display for use in determining a location of hydrocarbons, wherein the second image comprises at least a multiple image based at least in part on the migrated seismic data that improves the coverage at the shallow subsurface region receiving seismic data for a geologic region of the Earth, wherein the seismic data comprise primaries and multiples acquired using a seismic survey; building a velocity model that characterizes the geologic region of the Earth pre-processing the received seismic data by applying one or more multiple attenuation techniques to preserve primaries in the seismic data; generating a first image representative of coverage and quality of the primaries based on the pre-processed data; identifying a portion of the geologic region of the earth in the first image that comprises a deficiency in coverage, wherein the deficiency in coverage corresponds to one or more gaps between adjacent sail lines of the seismic survey in a shallow subsurface region; pre-processing the received seismic data a second time by applying one or more primary attenuation techniques to preserve one or more modes of multiples in the seismic data; selecting at least one mode of multiple and corresponding travel time data from a data storage , wherein the travel time data correspond to at least one complex ray signature in the geologic region of the earth and are based at least in part on the velocity model; performing Kirchhoff depth migration on the pre-processed seismic data with preserved modes of multiples, using at least the selected travel time data to generate migrated seismic data; and rendering a second image of the geologic region of the Earth to a display for use in determining a location of hydrocarbons, wherein the second image comprises at least a multiple image based at least in part on the migrated seismic data that improves the coverage at the shallow subsurface region <EOS>