Abstract: New methods for calculating acquisition illumination are computationally less expensive in comparison with conventional methods. In one such new method, source wavefield propagations are calculated and assigned to corresponding zero-offset receivers. Further, the number of non-zero-offset receivers within the coverage of the shot at the source location is decimated. Such a method is most advantages in reverse time migration, in which all source wavefield propagations are already calculated. The receiver-side illumination for each shot can be obtained by summing up all the source-side illumination with the source located within receiver coverage. All the source-side illumination and receiver-side illumination can be summed up to get the acquisition illumination for the survey. The acquisition illumination can be used to value the acquisition system and to compensate the migration images.

Abstract: New methods for calculating acquisition illumination are computationally less expensive in comparison with conventional methods. In one such new method, source wavefield propagations are calculated and assigned to corresponding zero-offset receivers. Further, the number of non-zero-offset receivers within the coverage of the shot at the source location is decimated. Such a method is most advantages in reverse time migration, in which all source wavefield propagations are already calculated. The receiver-side illumination for each shot can be obtained by summing up all the source-side illumination with the source located within receiver coverage. All the source-side illumination and receiver-side illumination can be summed up to get the acquisition illumination for the survey. The acquisition illumination can be used to value the acquisition system and to compensate the migration images.

Abstract: A method for seismic modeling on a GPU includes creating a discretized model containing a volume enclosed by boundaries with arbitrary boundary conditions, dividing the discretized model into subsets, transferring a data to the GPU, initializing the data stored in a selected number of the subsets on the GPU, performing computation related with a number of time steps for the selected number of the subsets on the GPU, and substantially simultaneously computing, off-loading a finished subsets data from the GPU, and up-loading a new subsets data to the GPU. The provided method and corresponding system can reduce the amount of data stored in GPU memory by an order of magnitude, minimize the GPU device-to-host transfer bandwidth requirements, and reduce the GPU setup time.

Abstract: In various disclosed seismic imaging systems and methods, a fast target-oriented illumination calculation technique is employed to obtain a data volume or “matrix” of approximate illumination values. These illumination values enable an image matrix of “true reflectivity” values to be found. The illumination values are derived from Green's functions which, rather than being calculated and re-calculated on a shot-by-shot basis, are calculated in multi-shot groups and combined with a rolling-sum to greatly reduce the computational overhead. As a consequence, the disclosed systems and methods can provide target region illuminations more quickly and/or with higher quality than those systems relying on conventional 3D wave-equation illumination.

Abstract: A method for seismic modeling on a GPU includes creating a discretized model containing a volume enclosed by boundaries with arbitrary boundary conditions, dividing the discretized model into subsets, transferring a data to the GPU, initializing the data stored in a selected number of the subsets on the GPU, performing computation related with a number of time steps for the selected number of the subsets on the GPU, and substantially simultaneously computing, off-loading a finished subsets data from the GPU, and up-loading a new subsets data to the GPU. The provided method and corresponding system can reduce the amount of data stored in GPU memory by an order of magnitude, minimize the GPU device-to-host transfer bandwidth requirements, and reduce the GPU setup time.

Abstract: A migration method using hybrid one-way and full-way (HOF) wave equation propagation. The HOF method extrapolates seismic wavefields in less complex media with a one-way wave equation propagator and extrapolates seismic wavefields in extremely complex media with a full-way wave equation propagator. For prestack depth migration, the HOF extrapolates source-side and receiver-side wavefields independently. Frequency-space domain and time-space domain imaging conditions are applied to the one-way and the full-way extrapolated wavefields, respectively. A suitable amplitude matching factor is introduced to combine the one-way and full-way images. The HOF method is a cost-effective migration that produces superior image quality with less noises and less computational resources.

Abstract: In various disclosed seismic imaging systems and methods, a fast target-oriented illumination calculation technique is employed to obtain a data volume or “matrix” of approximate illumination values. These illumination values enable an image matrix of “true reflectivity” values to be found. The illumination values are derived from Green's functions which, rather than being calculated and re-calculated on a shot-by-shot basis, are calculated in multi-shot groups and combined with a rolling-sum to greatly reduce the computational overhead. As a consequence, the disclosed systems and methods can provide target region illuminations more quickly and/or with higher quality than those systems relying on conventional 3D wave-equation illumination.