A time-domain SGFD-FK hybrid method for 2D teleseismic elastic wave modeling and inversion

IF 2.3 4区 地球科学
Mauricio del Valle-Rosales, Francisco José Chávez-García
{"title":"A time-domain SGFD-FK hybrid method for 2D teleseismic elastic wave modeling and inversion","authors":"Mauricio del Valle-Rosales,&nbsp;Francisco José Chávez-García","doi":"10.1007/s11600-024-01335-1","DOIUrl":null,"url":null,"abstract":"<div><p>Full waveform inversion (FWI) has proved to be a reliable tool for high-resolution imaging of lithospheric structures at various depths down to the upper mantle. However, when the size of the model is large, the computational burden is significative and applications are restricted to low frequencies. To tackle this issue, we developed a new 2D time-domain hybrid method to simulate high-frequency teleseismic body waves propagating through a local heterogeneous elastic Earth model: The frequency-wavenumber (FK) integration method is coupled with the staggered grid finite difference method (SGFD). The FK method is used to compute the wavefield due to obliquely incident plane P and SV waves in a 1D multilayered half space that excites a local heterogeneous region. Inside this region, the velocity-stress staggered grid FD method (SGFD) is used to accurately deal with wave propagation in heterogeneous media. Spurious waves that might be generated at the boundaries of the local region are avoided using convolutional perfectly matched layers (CPML). This new hybrid method inherits the low-memory requirements of the FK method and the accuracy, efficiency and easy implementation of the SGFD. The new hybrid method is benchmarked against the analytical FK method for some canonical models and shows good agreement with analytical solutions. Subsequently, our modeling tool is incorporated into a full waveform inversion algorithm adapted for teleseismic configurations to invert the incident P wave and its coda. The inversion is carried out using a gradient approach that is efficiently implemented via the adjoint-state method. The results suggest that our hybrid method and FWI algorithm represent a valuable tool for 2D forward and inverse regional applications using teleseismic data sets.</p></div>","PeriodicalId":6988,"journal":{"name":"Acta Geophysica","volume":"72 6","pages":"3903 - 3915"},"PeriodicalIF":2.3000,"publicationDate":"2024-04-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s11600-024-01335-1.pdf","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Acta Geophysica","FirstCategoryId":"89","ListUrlMain":"https://link.springer.com/article/10.1007/s11600-024-01335-1","RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0

Abstract

Full waveform inversion (FWI) has proved to be a reliable tool for high-resolution imaging of lithospheric structures at various depths down to the upper mantle. However, when the size of the model is large, the computational burden is significative and applications are restricted to low frequencies. To tackle this issue, we developed a new 2D time-domain hybrid method to simulate high-frequency teleseismic body waves propagating through a local heterogeneous elastic Earth model: The frequency-wavenumber (FK) integration method is coupled with the staggered grid finite difference method (SGFD). The FK method is used to compute the wavefield due to obliquely incident plane P and SV waves in a 1D multilayered half space that excites a local heterogeneous region. Inside this region, the velocity-stress staggered grid FD method (SGFD) is used to accurately deal with wave propagation in heterogeneous media. Spurious waves that might be generated at the boundaries of the local region are avoided using convolutional perfectly matched layers (CPML). This new hybrid method inherits the low-memory requirements of the FK method and the accuracy, efficiency and easy implementation of the SGFD. The new hybrid method is benchmarked against the analytical FK method for some canonical models and shows good agreement with analytical solutions. Subsequently, our modeling tool is incorporated into a full waveform inversion algorithm adapted for teleseismic configurations to invert the incident P wave and its coda. The inversion is carried out using a gradient approach that is efficiently implemented via the adjoint-state method. The results suggest that our hybrid method and FWI algorithm represent a valuable tool for 2D forward and inverse regional applications using teleseismic data sets.

Abstract Image

用于二维远震弹性波建模和反演的时域 SGFD-FK 混合方法
事实证明,全波形反演(FWI)是对上地幔不同深度岩石圈结构进行高分辨率成像的可靠工具。然而,当模型规模较大时,计算负担会很大,应用也仅限于低频。为了解决这个问题,我们开发了一种新的二维时域混合方法来模拟通过局部异质弹性地球模型传播的高频远震体波:频率-波数(FK)积分法与交错网格有限差分法(SGFD)相结合。FK 方法用于计算斜入射平面 P 波和 SV 波在激发局部异质区域的一维多层半空间中产生的波场。在该区域内,使用速度-应力交错网格 FD 方法 (SGFD) 精确处理波在异质介质中的传播。利用卷积完全匹配层(CPML)避免了可能在局部区域边界产生的杂波。这种新的混合方法继承了 FK 方法的低内存要求和 SGFD 的精确、高效和易于实施的特点。针对一些典型模型,新的混合方法与分析 FK 方法进行了基准测试,结果表明新方法与分析解法具有良好的一致性。随后,我们的建模工具被纳入了一个适用于远震配置的全波形反演算法,以反演入射 P 波及其尾波。反演采用梯度方法,通过邻接态方法有效实现。结果表明,我们的混合方法和全波反演算法是利用远震数据集进行二维正演和反演区域应用的重要工具。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
求助全文
约1分钟内获得全文 求助全文
来源期刊
Acta Geophysica
Acta Geophysica GEOCHEMISTRY & GEOPHYSICS-
CiteScore
3.80
自引率
13.00%
发文量
251
期刊介绍: Acta Geophysica is open to all kinds of manuscripts including research and review articles, short communications, comments to published papers, letters to the Editor as well as book reviews. Some of the issues are fully devoted to particular topics; we do encourage proposals for such topical issues. We accept submissions from scientists world-wide, offering high scientific and editorial standard and comprehensive treatment of the discussed topics.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
确定
请完成安全验证×
copy
已复制链接
快去分享给好友吧!
我知道了
右上角分享
点击右上角分享
0
联系我们:info@booksci.cn Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。 Copyright © 2023 布克学术 All rights reserved.
京ICP备2023020795号-1
ghs 京公网安备 11010802042870号
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术官方微信