A computer-implemented method of generating a seismic image and earth model related to a subsurface region of the earth, the method comprising: storing in a computer storage media, a data beam set generated from recorded seismic data that samples the subsurface region; and utilizing a computer system comprising one or more processors configured to communicate with the computer storage media and execute one or more computer programs configured for performing operations comprising: migrating the data beam set, comprising a plurality of data beams, using an initial earth model having an initial velocity model representative of the subsurface region, to generate a seismic image of the subsurface region; cross-correlating a modeled beam, derived from the seismic image and earth model, with a corresponding data beam within the data beam set; storing in the computer storage media beam pairs of cross-correlated modeled beams and data beams that for a measured time shift have a cross-correlation value that exceeds a determined threshold, and repeating the cross-correlating and storing operations until a determined number of corresponding data beams have been modeled; and inverting the measured time shifts to generate an updated earth model, having an updated velocity model, wherein the resulting updated earth model of the subsurface region is capable of generating an updated seismic image that accurately describes the propagation of the data beams through the subsurface region. storing in a computer storage media, a data beam set generated from recorded seismic data that samples the subsurface region; and utilizing a computer system comprising one or more processors configured to communicate with the computer storage media and execute one or more computer programs configured for performing operations comprising: migrating the data beam set, comprising a plurality of data beams, using an initial earth model having an initial velocity model representative of the subsurface region, to generate a seismic image of the subsurface region; cross-correlating a modeled beam, derived from the seismic image and earth model, with a corresponding data beam within the data beam set; storing in the computer storage media beam pairs of cross-correlated modeled beams and data beams that for a measured time shift have a cross-correlation value that exceeds a determined threshold, and repeating the cross-correlating and storing operations until a determined number of corresponding data beams have been modeled; and inverting the measured time shifts to generate an updated earth model, having an updated velocity model, wherein the resulting updated earth model of the subsurface region is capable of generating an updated seismic image that accurately describes the propagation of the data beams through the subsurface region. migrating the data beam set, comprising a plurality of data beams, using an initial earth model having an initial velocity model representative of the subsurface region, to generate a seismic image of the subsurface region; cross-correlating a modeled beam, derived from the seismic image and earth model, with a corresponding data beam within the data beam set; storing in the computer storage media beam pairs of cross-correlated modeled beams and data beams that for a measured time shift have a cross-correlation value that exceeds a determined threshold, and repeating the cross-correlating and storing operations until a determined number of corresponding data beams have been modeled; and inverting the measured time shifts to generate an updated earth model, having an updated velocity model, wherein the resulting updated earth model of the subsurface region is capable of generating an updated seismic image that accurately describes the propagation of the data beams through the subsurface region. <EOS>