Abstract:Prestack reverse time migration is a popular method for imaging underground structures in seismic exploration. Its imaging condition requires to gather data from the source wavefield and receiver wavefield samultaneously, which means that we have to store all source wavefield data at all times of the forward simulation and huge storage demand. The source wavefield reconstruction technologies can be used to solve this problem at the cost of introducing additional computational complexity. In this paper we introduce the POD finite element method and apply it to the reconstruction of viscous source wavefield, here the proper orthogonal decomposition (POD) method is used to reduce the storage demand by decreasing the dimensionality of data while keeping high enough accuracy. Numerical examples show that, compared with the traditional finite element methods, the POD finite element method can save more computer memory and speed up the reconstruction.