3D finite-difference frequency-domain modeling of visco-acoustic wave propagation using a massively parallel direct solver: A feasibility study
Top Cited Papers
- 1 September 2007
- journal article
- Published by Society of Exploration Geophysicists in Geophysics
- Vol. 72 (5), SM195-SM211
- https://doi.org/10.1190/1.2759835
Abstract
We present a finite-difference frequency-domain method for 3D visco-acoustic wave propagation modeling. In the frequency domain, the underlying numerical problem is the resolution of a large sparse system of linear equations whose right-hand side term is the source. This system is solved with a massively parallel direct solver. We first present an optimal 3D finite-difference stencil for frequency-domain modeling. The method is based on a parsimonious staggered-grid method. Differential operators are discretized with second-order accurate staggered-grid stencils on different rotated coordinate systems to mitigate numerical anisotropy. An antilumped mass strategy is implemented to minimize numerical dispersion. The stencil incorporates 27 grid points and spans two grid intervals. Dispersion analysis showsthat four grid points per wavelength provide accurate simulations in the 3D domain. To assess the feasibility of the method for frequency-domain full-waveform inversion, we computed simulations in the 3D SEG/EAGE overthrust model for frequencies 5, 7, and [Formula: see text]. Results confirm the huge memory requirement of the factorization (several hundred Figabytes) but also the CPU efficiency of the resolution phase (few seconds per shot). Heuristic scalability analysis suggests that the memory complexity of the factorization is [Formula: see text] for a [Formula: see text] grid. Our method may provide a suitable tool to perform frequency-domain full-waveform inversion using a large distributed-memory platform. Further investigation is still necessary to assess more quantitatively the respective merits and drawbacks of time- and frequency-domain modeling of wave propagation to perform 3D full-waveform inversion.Keywords
This publication has 33 references indexed in Scilit:
- Efficient waveform tomography for lithospheric imaging: implications for realistic, two-dimensional acquisition geometries and low-frequency dataGeophysical Journal International, 2007
- Hybrid scheduling for the parallel solution of linear systemsParallel Computing, 2006
- Mixed-grid and staggered-grid finite-difference methods for frequency-domain acoustic wave modellingGeophysical Journal International, 2004
- Impact of reordering on the memory of a multifrontal solverParallel Computing, 2003
- Robust Ordering of Sparse Matrices using MultisectionSIAM Journal on Matrix Analysis and Applications, 1998
- Applied Numerical Linear AlgebraPublished by Society for Industrial & Applied Mathematics (SIAM) ,1997
- 3-D Modeling Project: 3rd reportGeophysics, 1995
- A perfectly matched layer for the absorption of electromagnetic wavesJournal of Computational Physics, 1994
- SEG/EAEG 3-D modeling project: 2nd updateGeophysics, 1994
- The Multifrontal Solution of Indefinite Sparse Symmetric LinearACM Transactions on Mathematical Software, 1983