Full text loading...
-
Parallel Integral Equation Method and Algorithm for 3D Seismic Modelling
- Publisher: European Association of Geoscientists & Engineers
- Source: Conference Proceedings, 79th EAGE Conference and Exhibition 2017, Jun 2017, Volume 2017, p.1 - 5
Abstract
We present a novel parallel algorithm for acoustic field modeling in the frequency domain using the method of integral equation (IE). The IE method reduces the size of the problem by dividing the geologic model into the anomalous and background layered models. However, the matrix of the corresponding system of linear equations is dense. We use an efficient implementation of the matrix-vector multiplication, based on the fast Fourier transform (FFT), which made possible an application of an iterative solver to this dense system matrix. We have developed a parallel algorithm and studied its scalability on shared- and distributed-memory compute systems. The method was benchmarked by a comparison with a finite-difference solution. SEG/EAGE Overthrust model was used to evaluate the applicability of the algorithm to realistic seismic models.