Improvement to seismic processing based on predictive deconvolution

摘要

A seismic processing method 200 for determining a subsurface model comprises receiving 210 seismic data recorded in a space-time domain. A gapped deconvolution prediction operator- is derived 220 from the received seismic data in a model domain which is different from said space-time domain. The deconvolution prediction operator is converted 222 to the space-time domain using a linear conversion operator. The subsurface model is determined 223 by convolving said gapped deconvolution prediction operator with the received seismic data. The linear conversion operator may be a reverse linear Radon operator, reverse parabolic Radon operator, reverse hyperbolic Radon operator, reverse migration operator, reverse curvelet transform, reverse AVP operator, reverse tilted hyperbolic Radon operator. The subsurface model may comprise primaries, multiples or both.