An extended Gauss-Newton method for full waveform inversion

GEOPHYSICS Pub Date : 2024-02-15 DOI:10.1190/geo2022-0673.1
Ali Gholami
{"title":"An extended Gauss-Newton method for full waveform inversion","authors":"Ali Gholami","doi":"10.1190/geo2022-0673.1","DOIUrl":null,"url":null,"abstract":"Full waveform inversion (FWI) is a large-scale nonlinear ill-posed problem for which computationally expensive Newton-type methods can become trapped in undesirable local minima, particularly when the initial model lacks a low-wavenumber component and the recorded data lacks low-frequency content. A modification to the Gauss-Newton (GN) method is proposed to address these issues. The standard GN system for multisource multireceiver FWI is reformulated into an equivalent matrix equation form, with the solution becoming a diagonal matrix rather than a vector as in the standard system. The search direction is transformed from a vector to a matrix by relaxing the diagonality constraint, effectively adding a degree of freedom to the subsurface offset axis. The relaxed system can be explicitly solved with only the inversion of two small matrices that deblur the data residual matrix along the source and receiver dimensions, which simplifies the inversion of the Hessian matrix. When used to solve the extended source FWI objective function, the Extended GN (EGN) method integrates the benefits of both model and source extension. The EGN method effectively combines the computational effectiveness of the reduced FWI method with the robustness characteristics of extended formulations and offers a promising solution for addressing the challenges of FWI. It bridges the gap between these extended formulations and the reduced FWI method, enhancing inversion robustness while maintaining computational efficiency. The robustness and stability of the EGN algorithm for waveform inversion are demonstrated numerically.","PeriodicalId":509604,"journal":{"name":"GEOPHYSICS","volume":null,"pages":null},"PeriodicalIF":0.0000,"publicationDate":"2024-02-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"GEOPHYSICS","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/10.1190/geo2022-0673.1","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0

Abstract

Full waveform inversion (FWI) is a large-scale nonlinear ill-posed problem for which computationally expensive Newton-type methods can become trapped in undesirable local minima, particularly when the initial model lacks a low-wavenumber component and the recorded data lacks low-frequency content. A modification to the Gauss-Newton (GN) method is proposed to address these issues. The standard GN system for multisource multireceiver FWI is reformulated into an equivalent matrix equation form, with the solution becoming a diagonal matrix rather than a vector as in the standard system. The search direction is transformed from a vector to a matrix by relaxing the diagonality constraint, effectively adding a degree of freedom to the subsurface offset axis. The relaxed system can be explicitly solved with only the inversion of two small matrices that deblur the data residual matrix along the source and receiver dimensions, which simplifies the inversion of the Hessian matrix. When used to solve the extended source FWI objective function, the Extended GN (EGN) method integrates the benefits of both model and source extension. The EGN method effectively combines the computational effectiveness of the reduced FWI method with the robustness characteristics of extended formulations and offers a promising solution for addressing the challenges of FWI. It bridges the gap between these extended formulations and the reduced FWI method, enhancing inversion robustness while maintaining computational efficiency. The robustness and stability of the EGN algorithm for waveform inversion are demonstrated numerically.
用于全波形反演的扩展高斯-牛顿法
全波形反演(FWI)是一个大规模的非线性问题,对于这个问题,计算昂贵的牛顿型方法可能会陷入不理想的局部极小值,特别是当初始模型缺乏低频分量和记录数据缺乏低频内容时。为了解决这些问题,我们提出了对高斯-牛顿(GN)方法的一种改进。用于多源多接收器 FWI 的标准 GN 系统被重新表述为等效矩阵方程形式,解成为对角矩阵,而不是标准系统中的矢量。通过放宽对角约束,搜索方向从矢量转变为矩阵,从而有效地为地下偏移轴增加了一个自由度。放宽后的系统只需对两个小矩阵进行反演,就能显式求解,这两个小矩阵可使数据残差矩阵沿声源和接收器维度发生模糊,从而简化了赫塞斯矩阵的反演。当用于求解扩展源 FWI 目标函数时,扩展 GN(EGN)方法综合了模型和源扩展的优点。EGN 方法有效地结合了简化 FWI 方法的计算有效性和扩展公式的鲁棒性特点,为解决 FWI 面临的挑战提供了一个前景广阔的解决方案。它在这些扩展公式和简化 FWI 方法之间架起了一座桥梁,在保持计算效率的同时增强了反演的鲁棒性。数值演示了 EGN 算法在波形反演中的鲁棒性和稳定性。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
求助全文
约1分钟内获得全文 求助全文
来源期刊
自引率
0.00%
发文量
0
×
引用
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学术官方微信