Automatic phase correction of common-midpoint stacked data

The residual wavelet on a processed seismic section is often not zero phase despite all efforts to make it so. In this paper we adopt the convolutional model for the processed seismogram, assume that the residual phase shift can be approximated by a frequency-independent constant, and use the varimax norm to generate an algorithm to estimate the residual phase directly. Application of our algorithm to reflectivities from well logs suggests that it should work in the majority of cases so long as the reflectivity is non-Gaussian. An application of our algorithm to stacked data enhances the interpretability of the seismic section and leads to an improved match between the recovered relative acoustic impedance and a measured velocity log.