1887

Abstract

Summary

Frequency-domain full-waveform inversion (FWI) is suitable for long-offset stationary-recording acquisition, since reliable subsurface models can be reconstructed with a few frequencies and attenuation is easily implemented without computational overhead. In the frequency domain, wave modelling is a Helmholtz-type boundary-value problem which requires to solve a large and sparse linear system per frequency with multiple right-hand sides. This system can be solved with direct or iterative methods. While the former are suitable for FWI application on 3D dense OBC acquisitions covering spatial domains of moderate size, the later should be the approach of choice for sparse node acquisitions covering large domains (> 50 millions of unknowns). Fast convergence of iterative solvers for Helmholtz problems remains however challenging due to the non definiteness of the Helmholtz operator, hence requiring efficient preconditioners. Here, we use the Krylov subspace GMRES iterative solver combined with a multi-level domain-decomposition preconditioner. Discretization relies on finite elements on tetrahedral meshes to comply with complex geometries and adapt the size of the elements to the local wavelength (h-adaptivity). We assess the convergence and the scalability of our method with the 3D SEG/EAGE Overthrust model up to a frequency of 20Hz and discuss its efficiency for multi right-hand side processing.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.202011328
2020-12-08
2024-04-19
Loading full text...

Full text loading...

References

  1. Amestoy, P., Brossier, R., Buttari, A., L’Excellent, J.Y., Mary, T., Métivier, L., Miniussi, A. and Operto, S.
    [2016] Fast 3D frequency-domain FWI with a parallel Block Low-Rank multifrontal direct solver: application to OBC data from the North Sea.Geophysics, 81(6), R363 – R383.
    [Google Scholar]
  2. Beaudoin, G. and Ross, A.
    [2007] Field design and operation of a novel deepwater, wide-azimuth node seismic survey.The Leading Edge, 26, 494–503.
    [Google Scholar]
  3. Bonazzoli, M., Dolean, V., Graham, I.G., Spence, E.A. and Tournier, P.H.
    [2019] A 2-level domain decomposition preconditioner for the time-harmonic Maxwell’s equations.Math. Comp., 88, 2559–2604.
    [Google Scholar]
  4. Dolean, V., Jolivet, P. and Nataf, F.
    [2015] An introduction to domain decomposition methods - Algorithms, theory, and parallel implementation. SIAM.
    [Google Scholar]
  5. Duff, I.S., Erisman, A.M. and Reid, J.K.
    [1986] Direct methods for sparse matrices, second edition. Oxford Science Publications, Oxford, U. K.
    [Google Scholar]
  6. Erlangga, Y.A. and Nabben, R.
    [2008] On a multilevel Krylov method for the Helmholtz equation preconditioned by shifted Laplacian.Electronic Transactions on Numerical Analysis, 31, 403–424.
    [Google Scholar]
  7. Graham, G., Spence, E.A. and Vainikko, E.
    [2017] Recent results on domain decomposition preconditioning for the high-frequency Helmholtz equation using absorption. Springer, 3–26.
    [Google Scholar]
  8. Jolivet, P., Hecht, F., Nataf, F. and Prud’homme, C.
    [2013] Scalable Domain Decomposition Precon-ditioners for Heterogeneous Elliptic Problems. In: Proceedings of the International Conference on High Performance Computing, Networking, Storage and Analysis, SC’13. ACM, NY, USA.
    [Google Scholar]
  9. Jolivet, P. and Tournier, P.H.
    [2016] Block iterative methods and recycling for improved scalability of linear solvers. In: SC16: International Conference for High Performance Computing, Networking, Storage and Analysis. 190–203.
    [Google Scholar]
  10. Mary, T.
    [2017] Block Low-Rank multifrontal solvers: complexity, performance and scalability. Ph.D. thesis, PhD thesis, Université de Toulouse.
    [Google Scholar]
  11. Parks, M., de Sturler, E., Mackey, G., Johnson, D. and Maiti, S.
    [2006] Recycling krylov subspaces for sequences of linear systems.SIAM Journal of Scientific Computing, 28(5), 1651–1674.
    [Google Scholar]
  12. Plessix, R.E.
    [2017] Some computational aspects of the time and frequency domain formulations of seismic waveform inversion. Springer, 159–187.
    [Google Scholar]
  13. Pratt, R.G.
    [1999] Seismic waveform inversion in the frequency domain, part I: theory and verification in a physical scale model.Geophysics, 64, 888–901.
    [Google Scholar]
  14. Saad, Y.
    [2003] Iterative Methods for Sparse Linear Systems. SIAM, Philadelphia.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.202011328
Loading
/content/papers/10.3997/2214-4609.202011328
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