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...
Gespeichert in:
Hauptverfasser: | , , , , , |
---|---|
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 & 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 & 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 & 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 & 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 |