Nonstretch normal moveout through iterative partial correction and deconvolution

Source to receiver distances employed in seismic data acquisition have been steadily increasing and it is now common to work with data acquired with more than 10 km of offset. Subbasalt exploration and seismic undershooting are just two applications where long-offset reflections are sought. However, such reflections are often subjected to muting to suppress normal moveout (NMO) stretch artifacts, thus causing a loss of valuable information. To retrieve these portions of the recorded wavefield we develop a non-stretch NMO correction based on wavelet estimation and on an iterative procedure of partial NMO correction and deconvolution. We present this methodology using fourth-order traveltime curve approximations to increase the offset of useable reflections, but it can be adapted to traveltime curves of any order. Time- and space-variant wavelets, estimated by means of singular value decomposition (SVD) along the sought traveltimes, are used to build the desired output for the deconvolution that aims at retrieving the original shape of the partially stretched wavelets. We test our method on a synthetic gather presenting time and offset varying wavelets, on a real marine line simulating an undershooting pattern, and on a true undershooting land-marine data. These examples demonstrate that our new algorithm effectively limits the stretching associated with the NMO correction and enables the recovery of those portions of the stacked sections which are typically lost from muting in the standard NMO correction.