1887

Abstract

Summary

Complex topography, free surface boundary condition and inelastic properties of media should be well considered for onshore geophysical prospecting. Thus an appropriate and accurate forward modeling engine is very important. Unlike the time-domain implementation of many seismic imaging techniques, the counterpart in the frequency domain is rarely studied, despite of having many advantages, for example, only limited number of frequencies is needed for the inversion process, and solving the multiple-source problem is quite cheap if a direct solver is used. In this study, the spectral element method is applied to discretize the 3D frequency-domain anisotropic elastic wave modeling and the parallel direct solver MUMPS is used to solve the linear system. The structure and building process of the impedance matrix is thoroughly explained. We validate the numerical results by comparing with analytical solutions. A hybrid implementation of MPI and OpenMP for MUMPS is shown to be more efficient in flops and memory cost during the factorization. The influence of complex topography on MUMPS performance is negligible. With the available resources, the largest scale modeling, 30 wavelength in each dimension, is achieved. Other direct solvers and different low-rank techniques will also be investigated to reduce the flops and memory cost.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.201901534
2019-06-03
2024-03-29
Loading full text...

Full text loading...

References

  1. Deville, M., Fischer, P. and Mund, E.
    [2002] High Order Methods for Incompressible Fluid Flow. Cambridge University Press.
    [Google Scholar]
  2. Karypis, G.
    [2013] METIS - A Software Package for Partitioning Unstructured Graphs, Partitioning Meshes, and Computing Fill-Reducing Orderings of Sparse Matrices - Version 5.1.0.University of Minnesota.
    [Google Scholar]
  3. Komatitsch, D., Barnes, C. and Tromp, J.
    [2000] Simulation of anisotropic wave propagation based upon a spectral element method. Geophysics, 65(4), 1251–1260.
    [Google Scholar]
  4. Li, Y., Métivier, L., Brossier, R., Han, B. and Virieux, J.
    [2015] 2D and 3D frequency-domain elastic wave modeling in complex media with a parallel iterative solver. Geophysics, 80(3), T101–T118.
    [Google Scholar]
  5. MUMPS-team
    [2017] MUMPS - MUltifrontal Massively Parallel Solver users' guide - version 5.1.2 (October, 2017). ENSEEIHT-ENS Lyon, http://www.enseeiht.fr/apo/MUMPS/ or http://graal.ens-lyon.fr/MUMPS.
    [Google Scholar]
  6. Operto, S. and Miniussi, A.
    [2018] On the role of density and attenuation in 3D multi-parameter visco-acoustic VTI frequency-domain FWI: an OBC case study from the North Sea. Geophysical Journal International, 213, 2037–2059.
    [Google Scholar]
  7. Shi, L., Zhou, Y., Wang, J.M, Zhuang, M., Liu, N. and Liu, Q.H
    [2016] Spectral element method for elastic and acoustic waves in frequency domain. Journal of Computational Physics, 327(Supplement C), 19 – 38.
    [Google Scholar]
  8. Trinh, P.T, Brossier, R., Métivier, L., Tavard, L. and Virieux, J.
    [2019] Efficient 3D time-domain elastic and viscoelastic Full Waveform Inversion using a spectral-element method on flexible Cartesian-based mesh. Geophysics, 84(1), R75–R97.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.201901534
Loading
/content/papers/10.3997/2214-4609.201901534
Loading

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error