Parallel Numerical Simulation of Seismic Waves Propagation with Intel Math Kernel Library

This paper describes the implementation of parallel computing to model seismic waves in heterogeneous media based on Laguerre transform with respect to time. The main advantages of the transform are a definite sign of the spatial part of the operator and its independence of the parameter of separation. This property allows one to efficiently organize parallel computations by means of decomposition of the computational domain with successive application of the additive Schwarz method. At each step of the Schwarz alternations, a system of linear algebraic equations in each subdomain is resolved independently of all the others. A proper choice of Domain Decomposition reduces the size of matrices and ensures the use of direct solvers, in particular, the ones based on LU decomposition. Thanks to the independence of the matrix of the parameter of Laguerre transform with respect to time, LU decomposition for each subdomain is done only once, saved in the memory and used afterwards for different right-hand sides. A software is being developed for a cluster using hybrid OpenMP and MPI parallelization. At each cluster node, a system of linear algebraic equations with different right-hand sides is solved by the direct sparse solver PARDISO from Intel Math Kernel Library (Intel MKL). The solver is extensively parallelized and optimized for the high performance on many core systems with shared memory. A high performance parallel algorithm to solve the problem has been developed. The algorithm scalability and efficiency is investigated. For a two-dimensional heterogeneous medium, describing a realistic geological structure, which is typical of the North Sea, the results of numerical modeling are presented.

[1]  Satoshi Tsujimoto,et al.  Classical Orthogonal Polynomials and Bannai-Ito Polynomials , 2012 .

[2]  Chrysoula Tsogka,et al.  Application of the PML absorbing layer model to the linear elastodynamic problem in anisotropic hete , 1998 .

[3]  G. Reshetova,et al.  Numerical Modeling of Transient Seismic Fields in Viscoelastic Media Based on the Laguerre Spectral Method , 2003 .

[4]  H. Fossen,et al.  Structural geology of the Gullfaks Field, northern North Sea , 1998, Geological Society, London, Special Publications.

[5]  Olaf Schenk,et al.  Two-level dynamic scheduling in PARDISO: Improved scalability on shared memory multiprocessing systems , 2002, Parallel Comput..

[6]  T. Chan,et al.  Domain decomposition algorithms , 1994, Acta Numerica.

[7]  Wolfgang Fichtner,et al.  Efficient Sparse LU Factorization with Left-Right Looking Strategy on Shared Memory Multiprocessors , 2000 .

[8]  J. Virieux.  P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method , 1986 .

[9]  Tony F. Chan,et al.  Domain decomposition preconditioners for convection diffusion problems , 1994 .

[10]  C. Tsogka,et al.  Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media , 2001 .

[11]  René-Édouard Plessix,et al.  A Helmholtz iterative solver for 3D seismic-imaging problems , 2007 .

[12]  Martin J. Gander,et al.  Optimized Schwarz Methods , 2006, SIAM J. Numer. Anal..

[13]  Vipin Kumar,et al.  A Fast and High Quality Multilevel Scheme for Partitioning Irregular Graphs , 1998, SIAM J. Sci. Comput..

[14]  J. Zahradník,et al.  Heterogeneous formulations of elastodynamic equations and finite-difference schemes , 1995 .

[15]  Hongyuan Zha,et al.  Low-Rank Approximations with Sparse Factors I: Basic Algorithms and Error Analysis , 2001, SIAM J. Matrix Anal. Appl..

[16]  S. Nepomnyaschikh,et al.  Domain Decomposition Methods , 2007 .

[17]  Phillip Colella,et al.  Performance and scaling of locally-structured grid methods for partial differential equations - eScholarship , 2008 .