Method of seismic time migration using a massively parallel computer

摘要

A method for 3 dimensional, one pass migration of post stack seismic data is implemented on a massively parallel computer. The quality of the migration and speed of execution makes possible iterative 3D post stack depth migrations. The depth migration method computes in the frequency domain and downward continues one or more frequency planes through all depth levels of interest. For a single frequency, at each depth level, the method extrapolates a shifted position for each x,y position by applying a filter for 2D convolution. A processing component is assigned to each subgrid of x,y spatial positions and the processing components operate concurrently in parallel to determine filter coefficients independently for each x,y spatial position and to extrapolate the x,y shifted positions. The filter coefficients are derived independently at each x,y position by storing a limited table of filter coefficients in local memory of each processing component.