A computer-implemented method of seismic imaging for identification of and production from hydrocarbon reservoirs, comprising: a. receiving, at a computer processor, a seismic dataset representative of a subsurface volume of interest; b. receiving, at the computer processor, a representation of geological constraints for the subsurface volume of interest, the geological constraints comprising a geologically plausible reflector sourced independently of the seismic dataset; c. transforming, via the computer processor, the seismic dataset into a data beam set; d. receiving, at the computer processor, an earth model, the earth model being a model of the physical properties of the subsurface volume of interest; e. computing, via the computer processor, traveltime variations from the data beam set based on the geological constraints and the earth model, the traveltime variations being computed by computing traveltimes from an undesired surface of a geologically implausible reflector of the seismic dataset to a desired surface of the geologically plausible reflector sourced independently of the seismic dataset along ray-paths in directions normal to said undesired surface based on the earth model, and, based on the directions of rays, obtaining desired time shifts required to refocus from the undesired surface to the desired surface from the traveltimes; f. deriving, via the computer processor, an updated earth model based on the traveltime variations; g. performing, via the computer processor, seismic imaging of the seismic dataset using the updated earth model to generate a seismic image with correctly focused and positioned reflectors; and h. identifying geologic features based on the seismic image. a. receiving, at a computer processor, a seismic dataset representative of a subsurface volume of interest; b. receiving, at the computer processor, a representation of geological constraints for the subsurface volume of interest, the geological constraints comprising a geologically plausible reflector sourced independently of the seismic dataset; c. transforming, via the computer processor, the seismic dataset into a data beam set; d. receiving, at the computer processor, an earth model, the earth model being a model of the physical properties of the subsurface volume of interest; e. computing, via the computer processor, traveltime variations from the data beam set based on the geological constraints and the earth model, the traveltime variations being computed by computing traveltimes from an undesired surface of a geologically implausible reflector of the seismic dataset to a desired surface of the geologically plausible reflector sourced independently of the seismic dataset along ray-paths in directions normal to said undesired surface based on the earth model, and, based on the directions of rays, obtaining desired time shifts required to refocus from the undesired surface to the desired surface from the traveltimes; f. deriving, via the computer processor, an updated earth model based on the traveltime variations; g. performing, via the computer processor, seismic imaging of the seismic dataset using the updated earth model to generate a seismic image with correctly focused and positioned reflectors; and h. identifying geologic features based on the seismic image.