Abstract
The phase‐shift method of wavefield extrapolation applies a phase shift in the Fourier domain to deduce a scalar wavefield at one depth level given its value at another. The phase‐shift operator varies with frequency and wavenumber, and assumes constant velocity across the extrapolation step. We use nonstationary filter theory to generalize this method to nonstationary phase shift (NSPS), which allows the phase shift to vary laterally depending upon the local propagation velocity. For comparison, we derive an analytic form for the popular phase shift plus interpolation (PSPI) method in the limit of an exhaustive set of reference velocities. NSPS and this limiting form of PSPI can be written as generalized Fourier integrals which reduce to ordinary phase shift in the constant velocity limit. In the (x, ω) domain, these processes are the transpose of each other; however, only NSPS has the physical interpretation of forming the scaled, linear superposition of laterally‐variable impulse responses (i.e.,Huygen’s wavelets). The difference between NSPS and PSPI is clear when they are compared in the case of a piecewise constant velocity variation. Define a set of windows such that the jth window is unity when the propagation velocity is the jth distinct velocity and is zero otherwise. NSPS can be computed by applying the window set to the input data to create a set of windowed wavefields, which are individually phase‐shift extrapolated with the corresponding constant velocity, and the extrapolated set is superimposed. PSPI proceeds by phase‐shift extrapolating the input data for each distinct velocity, applying the jth window to the jth extrapolation, and superimposing. Though neither process is fully correct, PSPI has the unphysical limit that discontinuities in the lateral velocity variation cause discontinuities in the wavefield, whereas NSPS shows the expected wavefront “healing.” We then formulate a finite aperture compensation for NSPS which has the practical result of absorbing lateral boundaries for all incidence angles. Wavefield extrapolation can be regarded as the crosscorrelation of the wavefield with the expected response of a point diffractor at the new depth level. Aperture compensation simply applies a laterally varying window to the infinite, theoretical diffraction response. The crosscorrelation becomes spatially variant, even for constant velocity, and hence is a nonstationary filter. The nonstationary effects of aperture compensation can be simultaneously applied with the NSPS extrapolation through a laterally variable velocity field.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.