Method based on graphics processor for computing direct pre-stack reverse-time migration based on irregular topography

Abstract

The invention discloses a method based on a graphics processor for computing direct pre-stack reverse-time migration based on an irregular topography. Migration imaging is directly performed from the irregular topography by applying a reverse-time migration method, and solving is performed by adopting a finite difference method, and the computation speed and efficiency is obviously improved by utilizing the concurrent computation characteristic of the graphics processor.

Description

technical field [0001] The invention relates to a method for calculating and calculating direct pre-stack reverse time migration of undulating land surface, in particular to a method for calculating and calculating direct pre-stack reverse time migration of undulating land surface based on a graphic processor. Background technique [0002] With the continuous development of seismic exploration technology, the focus of oil and gas seismic exploration is shifting to areas with undulating surface conditions and complex geological conditions, such as mountains, beaches, and swamps, which pose new challenges to seismic exploration and data processing. Traditional seismic imaging methods mainly have two problems in areas with complex surface and complex structures: one is that the surface has large undulations, and the surface velocity structure is very complex; the other is that the underground structure is complex, such as strong folds, fault development, steep structures, and st...

AI Technical Summary

Problems solved by technology

Traditional seismic imaging methods mainly have two problems in areas with complex surface and complex structures: one is that the surface has large undulations, and the surface velocity structure is very complex; the other is that the underground structure is complex, such as strong folds, fault development, steep structures, and formation changes Big and so on

[0004] For the solution of the two-way wave wave equation on the undulating surface, the finite element method is more convenient to deal with the irregular surface, but the calculation efficiency is low, the calculation accuracy is difficult to improve, and it is not convenient to use

The combination of finite element and other methods, such as discrete wavenumber and domain splitting, has improved the calculation efficiency to a certain extent, but the problems of low calculation accuracy and inconvenient use have not been fundamentally solved.

The semi-analytical nature of the boundary integral or boundary element method determines that this method cannot be applied to the case of large surface velocity changes, which limits its application in actual exploration

Images