Parallel algorithm for 3D modeling of monochromatic acoustic field by the method of integral equations


  • Mikhail S. Malovichko Moscow Institute of Physics and Technology, Moscow
  • Nikolay E. Khokhlov Moscow Institute of Physics and Technology, Moscow
  • Nikolay B. Yavich Moscow Institute of Physics and Technology, Moscow
  • Michael S. Zhdanov The University of Utah, Salt Lake City



We present a parallel algorithm for solution of the three-dimensional Helmholtz equation in the frequency domain by the method of volume integral equations.
The algorithm is applied to seismic forward modeling.
The method of integral equations reduces the size of the problem by dividing the geologic model into the anomalous and background parts, but leads to a dense system matrix.
A tolerable memory consumption and numerical complexity were achieved by applying an iterative solver, accompanied by an effective matrix-vector multiplication operation, based on the fast Fourier transform.
We used OpenMP to speed up the matrix-vector multiplication, while MPI was used to speed up the equation system solver, and also for parallelizing across multiple sources.
Practical examples and efficiency tests are presented.


Abubakar, A., T.M. Habashy. Three-dimensional visco-acoustic modeling using a renormalized integral equation iterative solver ,Journal of computational physics, vol. 249, pp.1-12, 2013.

Aki, K. and P. G. Richards. Quantitative seismology. 1980. W. R. Freeman and Co.

Aminzadeh, F., Brac, J., Kunz, T. 3-D Salt and Overthrust Models. SEG/EAGE 3-D Modeling Series No.1., Soc. Explor. Geophysicists, Tulsa, 1997.

Freter, H. An Integral Equation Method for Seismic Modelling// in Inversion Theory and Practice of Geophysical Data Inversion, vol. 5, 1992. Vieweg+Teubner Verlag.

Fu, Li-Yun. Numerical study of generalized Lippmann-Schwinger integral equation including surface topography. //Geophysics, vol. 68, no. 2, pp.665-671, 2003.

Fu, Li-Yun, Yong-Guang Mu, Huey-Ju Yang. Forward problem of nonlinear Fredholm integral equation in reference medium via velocity-weighted wavefield function. //Geophysics, vol. 62, no. 2, pp.650-656, 1997.

Grama, A., A. Gupta, G. Karypis, V. Kumar Introduction to parallel computing, 2nd ed. Addison Wesley, 2003.

Johnson, S.A., Y. Zhou, M.J. Berggren, M.L. Tracy. Acoustic Inverse Scattering Solutions by Moment Methods and Back Propagation // Conference on inverse scattering: theory and applications. 1983. SIAM.

Morse, P. M. and Feshbach, H. Methods in theoretical physics. 1953, McGraw-Hill Book Company, inc.

Saad, Y. Iterative methods for sparse linear systems, 2nd ed., SIAM, Philadelphia, PA, USA. 2003.

Zhang, Rongfeng , Tadeusz J. Ulrych. Seismic forward modeling by integral equation and some practical considerations. //SEG Technical Program Expanded Abstracts 2000, chp. 593, pp. 2329-2332, SEG, 2000.




How to Cite

Malovichko, M. S., Khokhlov, N. E., Yavich, N. B., & Zhdanov, M. S. (2016). Parallel algorithm for 3D modeling of monochromatic acoustic field by the method of integral equations. Supercomputing Frontiers and Innovations, 3(4), 74–78.