Three-dimensional finite-difference & finite-element frequency-domain wave simulation with multi-level optimized additive Schwarz domain-decomposition preconditioner: A tool for FWI of sparse node datasets

Efficient frequency-domain Full Waveform Inversion (FWI) of long-offset node data can be designed with a few discrete frequencies hence allowing for compact volume of data to be managed. Moreover, attenuation effects can be straightforwardly implemented in the forward problem without computational o...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Tournier, P. -H, Jolivet, P, Dolean, V, Aghamiry, H. S, Operto, S, Riffo, S
Format: Artikel
Sprache:eng
Schlagworte:
Online-Zugang:Volltext bestellen
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
container_end_page
container_issue
container_start_page
container_title
container_volume
creator Tournier, P. -H
Jolivet, P
Dolean, V
Aghamiry, H. S
Operto, S
Riffo, S
description Efficient frequency-domain Full Waveform Inversion (FWI) of long-offset node data can be designed with a few discrete frequencies hence allowing for compact volume of data to be managed. Moreover, attenuation effects can be straightforwardly implemented in the forward problem without computational overhead. However, 3D frequency-domain seismic modeling is challenging since it requires solving a large and sparse linear indefinite system per frequency with multiple right-hand sides. This linear system can be solved by direct or iterative methods. The former are very efficient to process multiple right-hand sides but may suffer from limited scalability for very large problems. Iterative methods equipped with a domain decomposition preconditioner provide a suitable alternative to process large computational domains for sparse node acquisition. The domain decomposition preconditioner relies on an optimized restricted additive Schwarz (ORAS) method, where a Robin or Perfectly-Matched Layer (PML) condition is implemented at the boundaries between the subdomains. The preconditioned system is solved by a Krylov subspace method, while a block low-rank Lower-Upper (LU) decomposition of the local matrices is performed at a preprocessing stage. Multiple sources are processed in group with a pseudo-block method. The accuracy, the computational cost and the scalability of the ORAS solver are assessed against several realistic benchmarks. For the considered benchmarks, a compact wavelength-adaptive 27-point finite-difference stencil on regular Cartesian grid shows better accuracy and computational efficiency than a P3 finite-element method on h-adaptive tetrahedral mesh, which remains however beneficial to comply with known boundaries such as bathymetry.
doi_str_mv 10.48550/arxiv.2110.15113
format Article
fullrecord <record><control><sourceid>arxiv_GOX</sourceid><recordid>TN_cdi_arxiv_primary_2110_15113</recordid><sourceformat>XML</sourceformat><sourcesystem>PC</sourcesystem><sourcerecordid>2110_15113</sourcerecordid><originalsourceid>FETCH-LOGICAL-a673-f6c0917d8a093efc8555d8a74afb78c6b05c0be655a5bc322683fecf8c50f6e33</originalsourceid><addsrcrecordid>eNo1kLFOwzAQhrMwoMIDMHETW0pS10nKVlUUKlVioBJjdbHP6klJHGy3pX1H3gk3wHT6_jv9On1Jcpdn42klZfaI7osP40keg1zmubhOvjc7R5RqbqnzbDtswHDH4RIZQ446RfDwn1FD8S6AcfS5j6tTqm2L3MERDwSe232DIbbAkcMOIgVOGzpQA7YP3PKZNKDWHDiev6vdEd0ZfitSTcq2vfU8FPQuYqcHIPcEcwjWxt-sg-XHCqwB36PzBJ3VBBoDegr-Jrky2Hi6_ZujZLN83ixe0_Xby2oxX6dYlCI1hcpmeakrzGaCjIpqZIRyiqYuK1XUmVRZTYWUKGslJpOiEoaUqZTMTEFCjJL739rB57Z33KI7bS9et4NX8QMrs3oI</addsrcrecordid><sourcetype>Open Access Repository</sourcetype><iscdi>true</iscdi><recordtype>article</recordtype></control><display><type>article</type><title>Three-dimensional finite-difference &amp; finite-element frequency-domain wave simulation with multi-level optimized additive Schwarz domain-decomposition preconditioner: A tool for FWI of sparse node datasets</title><source>arXiv.org</source><creator>Tournier, P. -H ; Jolivet, P ; Dolean, V ; Aghamiry, H. S ; Operto, S ; Riffo, S</creator><creatorcontrib>Tournier, P. -H ; Jolivet, P ; Dolean, V ; Aghamiry, H. S ; Operto, S ; Riffo, S</creatorcontrib><description>Efficient frequency-domain Full Waveform Inversion (FWI) of long-offset node data can be designed with a few discrete frequencies hence allowing for compact volume of data to be managed. Moreover, attenuation effects can be straightforwardly implemented in the forward problem without computational overhead. However, 3D frequency-domain seismic modeling is challenging since it requires solving a large and sparse linear indefinite system per frequency with multiple right-hand sides. This linear system can be solved by direct or iterative methods. The former are very efficient to process multiple right-hand sides but may suffer from limited scalability for very large problems. Iterative methods equipped with a domain decomposition preconditioner provide a suitable alternative to process large computational domains for sparse node acquisition. The domain decomposition preconditioner relies on an optimized restricted additive Schwarz (ORAS) method, where a Robin or Perfectly-Matched Layer (PML) condition is implemented at the boundaries between the subdomains. The preconditioned system is solved by a Krylov subspace method, while a block low-rank Lower-Upper (LU) decomposition of the local matrices is performed at a preprocessing stage. Multiple sources are processed in group with a pseudo-block method. The accuracy, the computational cost and the scalability of the ORAS solver are assessed against several realistic benchmarks. For the considered benchmarks, a compact wavelength-adaptive 27-point finite-difference stencil on regular Cartesian grid shows better accuracy and computational efficiency than a P3 finite-element method on h-adaptive tetrahedral mesh, which remains however beneficial to comply with known boundaries such as bathymetry.</description><identifier>DOI: 10.48550/arxiv.2110.15113</identifier><language>eng</language><subject>Mathematics - Analysis of PDEs</subject><creationdate>2021-10</creationdate><rights>http://creativecommons.org/licenses/by-nc-nd/4.0</rights><oa>free_for_read</oa><woscitedreferencessubscribed>false</woscitedreferencessubscribed></display><links><openurl>$$Topenurl_article</openurl><openurlfulltext>$$Topenurlfull_article</openurlfulltext><thumbnail>$$Tsyndetics_thumb_exl</thumbnail><link.rule.ids>228,230,781,886</link.rule.ids><linktorsrc>$$Uhttps://arxiv.org/abs/2110.15113$$EView_record_in_Cornell_University$$FView_record_in_$$GCornell_University$$Hfree_for_read</linktorsrc><backlink>$$Uhttps://doi.org/10.48550/arXiv.2110.15113$$DView paper in arXiv$$Hfree_for_read</backlink></links><search><creatorcontrib>Tournier, P. -H</creatorcontrib><creatorcontrib>Jolivet, P</creatorcontrib><creatorcontrib>Dolean, V</creatorcontrib><creatorcontrib>Aghamiry, H. S</creatorcontrib><creatorcontrib>Operto, S</creatorcontrib><creatorcontrib>Riffo, S</creatorcontrib><title>Three-dimensional finite-difference &amp; finite-element frequency-domain wave simulation with multi-level optimized additive Schwarz domain-decomposition preconditioner: A tool for FWI of sparse node datasets</title><description>Efficient frequency-domain Full Waveform Inversion (FWI) of long-offset node data can be designed with a few discrete frequencies hence allowing for compact volume of data to be managed. Moreover, attenuation effects can be straightforwardly implemented in the forward problem without computational overhead. However, 3D frequency-domain seismic modeling is challenging since it requires solving a large and sparse linear indefinite system per frequency with multiple right-hand sides. This linear system can be solved by direct or iterative methods. The former are very efficient to process multiple right-hand sides but may suffer from limited scalability for very large problems. Iterative methods equipped with a domain decomposition preconditioner provide a suitable alternative to process large computational domains for sparse node acquisition. The domain decomposition preconditioner relies on an optimized restricted additive Schwarz (ORAS) method, where a Robin or Perfectly-Matched Layer (PML) condition is implemented at the boundaries between the subdomains. The preconditioned system is solved by a Krylov subspace method, while a block low-rank Lower-Upper (LU) decomposition of the local matrices is performed at a preprocessing stage. Multiple sources are processed in group with a pseudo-block method. The accuracy, the computational cost and the scalability of the ORAS solver are assessed against several realistic benchmarks. For the considered benchmarks, a compact wavelength-adaptive 27-point finite-difference stencil on regular Cartesian grid shows better accuracy and computational efficiency than a P3 finite-element method on h-adaptive tetrahedral mesh, which remains however beneficial to comply with known boundaries such as bathymetry.</description><subject>Mathematics - Analysis of PDEs</subject><fulltext>true</fulltext><rsrctype>article</rsrctype><creationdate>2021</creationdate><recordtype>article</recordtype><sourceid>GOX</sourceid><recordid>eNo1kLFOwzAQhrMwoMIDMHETW0pS10nKVlUUKlVioBJjdbHP6klJHGy3pX1H3gk3wHT6_jv9On1Jcpdn42klZfaI7osP40keg1zmubhOvjc7R5RqbqnzbDtswHDH4RIZQ446RfDwn1FD8S6AcfS5j6tTqm2L3MERDwSe232DIbbAkcMOIgVOGzpQA7YP3PKZNKDWHDiev6vdEd0ZfitSTcq2vfU8FPQuYqcHIPcEcwjWxt-sg-XHCqwB36PzBJ3VBBoDegr-Jrky2Hi6_ZujZLN83ixe0_Xby2oxX6dYlCI1hcpmeakrzGaCjIpqZIRyiqYuK1XUmVRZTYWUKGslJpOiEoaUqZTMTEFCjJL739rB57Z33KI7bS9et4NX8QMrs3oI</recordid><startdate>20211028</startdate><enddate>20211028</enddate><creator>Tournier, P. -H</creator><creator>Jolivet, P</creator><creator>Dolean, V</creator><creator>Aghamiry, H. S</creator><creator>Operto, S</creator><creator>Riffo, S</creator><scope>AKZ</scope><scope>GOX</scope></search><sort><creationdate>20211028</creationdate><title>Three-dimensional finite-difference &amp; finite-element frequency-domain wave simulation with multi-level optimized additive Schwarz domain-decomposition preconditioner: A tool for FWI of sparse node datasets</title><author>Tournier, P. -H ; Jolivet, P ; Dolean, V ; Aghamiry, H. S ; Operto, S ; Riffo, S</author></sort><facets><frbrtype>5</frbrtype><frbrgroupid>cdi_FETCH-LOGICAL-a673-f6c0917d8a093efc8555d8a74afb78c6b05c0be655a5bc322683fecf8c50f6e33</frbrgroupid><rsrctype>articles</rsrctype><prefilter>articles</prefilter><language>eng</language><creationdate>2021</creationdate><topic>Mathematics - Analysis of PDEs</topic><toplevel>online_resources</toplevel><creatorcontrib>Tournier, P. -H</creatorcontrib><creatorcontrib>Jolivet, P</creatorcontrib><creatorcontrib>Dolean, V</creatorcontrib><creatorcontrib>Aghamiry, H. S</creatorcontrib><creatorcontrib>Operto, S</creatorcontrib><creatorcontrib>Riffo, S</creatorcontrib><collection>arXiv Mathematics</collection><collection>arXiv.org</collection></facets><delivery><delcategory>Remote Search Resource</delcategory><fulltext>fulltext_linktorsrc</fulltext></delivery><addata><au>Tournier, P. -H</au><au>Jolivet, P</au><au>Dolean, V</au><au>Aghamiry, H. S</au><au>Operto, S</au><au>Riffo, S</au><format>journal</format><genre>article</genre><ristype>JOUR</ristype><atitle>Three-dimensional finite-difference &amp; finite-element frequency-domain wave simulation with multi-level optimized additive Schwarz domain-decomposition preconditioner: A tool for FWI of sparse node datasets</atitle><date>2021-10-28</date><risdate>2021</risdate><abstract>Efficient frequency-domain Full Waveform Inversion (FWI) of long-offset node data can be designed with a few discrete frequencies hence allowing for compact volume of data to be managed. Moreover, attenuation effects can be straightforwardly implemented in the forward problem without computational overhead. However, 3D frequency-domain seismic modeling is challenging since it requires solving a large and sparse linear indefinite system per frequency with multiple right-hand sides. This linear system can be solved by direct or iterative methods. The former are very efficient to process multiple right-hand sides but may suffer from limited scalability for very large problems. Iterative methods equipped with a domain decomposition preconditioner provide a suitable alternative to process large computational domains for sparse node acquisition. The domain decomposition preconditioner relies on an optimized restricted additive Schwarz (ORAS) method, where a Robin or Perfectly-Matched Layer (PML) condition is implemented at the boundaries between the subdomains. The preconditioned system is solved by a Krylov subspace method, while a block low-rank Lower-Upper (LU) decomposition of the local matrices is performed at a preprocessing stage. Multiple sources are processed in group with a pseudo-block method. The accuracy, the computational cost and the scalability of the ORAS solver are assessed against several realistic benchmarks. For the considered benchmarks, a compact wavelength-adaptive 27-point finite-difference stencil on regular Cartesian grid shows better accuracy and computational efficiency than a P3 finite-element method on h-adaptive tetrahedral mesh, which remains however beneficial to comply with known boundaries such as bathymetry.</abstract><doi>10.48550/arxiv.2110.15113</doi><oa>free_for_read</oa></addata></record>
fulltext fulltext_linktorsrc
identifier DOI: 10.48550/arxiv.2110.15113
ispartof
issn
language eng
recordid cdi_arxiv_primary_2110_15113
source arXiv.org
subjects Mathematics - Analysis of PDEs
title Three-dimensional finite-difference & finite-element frequency-domain wave simulation with multi-level optimized additive Schwarz domain-decomposition preconditioner: A tool for FWI of sparse node datasets
url https://sfx.bib-bvb.de/sfx_tum?ctx_ver=Z39.88-2004&ctx_enc=info:ofi/enc:UTF-8&ctx_tim=2024-12-18T13%3A03%3A43IST&url_ver=Z39.88-2004&url_ctx_fmt=infofi/fmt:kev:mtx:ctx&rfr_id=info:sid/primo.exlibrisgroup.com:primo3-Article-arxiv_GOX&rft_val_fmt=info:ofi/fmt:kev:mtx:journal&rft.genre=article&rft.atitle=Three-dimensional%20finite-difference%20&%20finite-element%20frequency-domain%20wave%20simulation%20with%20multi-level%20optimized%20additive%20Schwarz%20domain-decomposition%20preconditioner:%20A%20tool%20for%20FWI%20of%20sparse%20node%20datasets&rft.au=Tournier,%20P.%20-H&rft.date=2021-10-28&rft_id=info:doi/10.48550/arxiv.2110.15113&rft_dat=%3Carxiv_GOX%3E2110_15113%3C/arxiv_GOX%3E%3Curl%3E%3C/url%3E&disable_directlink=true&sfx.directlink=off&sfx.report_link=0&rft_id=info:oai/&rft_id=info:pmid/&rfr_iscdi=true