Accurate 3D frequency-domain seismic wave modeling with the wavelength-adaptive 27-point finite-difference stencil: a tool for full waveform inversion - Archive ouverte HAL Accéder directement au contenu
Article Dans Une Revue Geophysics Année : 2022

Accurate 3D frequency-domain seismic wave modeling with the wavelength-adaptive 27-point finite-difference stencil: a tool for full waveform inversion

Hossein Aghamiry
Ali Gholami
Laure Combe
  • Fonction : Auteur
  • PersonId : 752975
  • IdHAL : laure-combe
Stéphane Operto

Résumé

Efficient frequency-domain Full Waveform Inversion (FWI) of long-offset node data can be performed with a few frequencies. The seismic response of these frequencies can be computed with compact finite-difference stencils on regular Cartesian grid with direct or hybrid direct/iterative methods. Compactness, which is necessary to mitigate the fill-in induced by the Lower-Upper (LU) factorization, is implemented with second-order stencils while accuracy is achieved by building a consistent mass matrix and a compound stiffness matrix with optimal weights so that the stencil covers the eight cells surrounding the central point, leading to a 27-point stencil. Classical approaches estimate constant weights by jointly minimizing numerical dispersion in homogeneous media for several numbers of grid points per wavelength ( G). Then, the impedance matrix is built by using these constant weights at each central point covering the heterogeneous subsurface model, leading to non-uniform wavefield accuracy. Instead, we estimate G-dependent weights once and for all by minimizing dispersion separately for each value of 1/ G in a range found in FWI applications. Then, we build the impedance matrix without computational overhead by selecting at each central point the weights corresponding to the local 1/ G, hence leading to a wavelength-adaptive stencil. This separate approach with adaptive weights makes wavefield accuracy uniform.We benchmark wavefields, which are computed in several 3D large-scale subsurface models with a sparse multifrontal direct solver and the non-adaptive/adaptive stencils, against analytical solutions when available and the highly accurate discretization-free convergent Born series (CBS) method. Each benchmark reveals the higher accuracy of the adaptive stencil relative to the non-adaptive one. In the presence of sharp contrasts, the adaptive stencil is also more accurate than a classical O(Δ t 2 ,Δ h 8 ) finite-different time-domain staggered-grid stencil. The relevance of the adaptive stencil is finally illustrated with a 3D FWI case study in the 3.5–13 Hz frequency ba

Dates et versions

hal-03592166 , version 1 (01-03-2022)

Identifiants

Citer

Hossein Aghamiry, Ali Gholami, Laure Combe, Stéphane Operto. Accurate 3D frequency-domain seismic wave modeling with the wavelength-adaptive 27-point finite-difference stencil: a tool for full waveform inversion. Geophysics, 2022, pp.1-66. ⟨10.1190/geo2021-0606.1⟩. ⟨hal-03592166⟩
17 Consultations
0 Téléchargements

Altmetric

Partager

Gmail Facebook X LinkedIn More