3D Gravity Fast Inversion Based on Krylov Subspace Methods

IF 1.6 3区 地球科学 Q3 GEOCHEMISTRY & GEOPHYSICS
Min Yang, Xinqiang Xu, Wanyin Wang, Dongming Zhao, Wei Zhou
{"title":"3D Gravity Fast Inversion Based on Krylov Subspace Methods","authors":"Min Yang, Xinqiang Xu, Wanyin Wang, Dongming Zhao, Wei Zhou","doi":"10.1093/jge/gxad091","DOIUrl":null,"url":null,"abstract":"Abstract Mapping the density contrast through the 3D gravity inversion can help detect the goals under the subsurface. However, it is a challenge to accurately and efficiently solve the 3D gravity inversion. Krylov subspace method is commonly used for large linear problems due to its high computational efficiency and low storage requirement. In this study, two classical algorithms of Krylov subspace method, namely the Generalized Minimum Residual method and the Conjugate Gradient method, are applied to 3D gravity inversion. Based on the recovered models of the deep mineral and the shallow L-shaped tunnel models, it was found that the Generalized Minimum Residual method provided similar density contrast results as the Conjugate Gradient method. The obtained inversion results of density contrast corresponded well to the position of the deep mineral resources model and the L-shaped tunnel model. The 3D distribution of Fe content underground was obtained by inverting the measured gravity data from Olympic Dam in Australia. The recovered results correspond well with the distribution of Fe content in the geological profile collected. The accuracy of inversion using the Generalized Minimum Residual method was similar to that of the Conjugate Gradient method under the same conditions. However, the Generalized Minimum Residual method had a faster convergence speed and increased inversion efficiency by about 90%, greatly reducing the inversion time and improves the inversion efficiency.","PeriodicalId":54820,"journal":{"name":"Journal of Geophysics and Engineering","volume":"132 45","pages":"0"},"PeriodicalIF":1.6000,"publicationDate":"2023-11-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Geophysics and Engineering","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/10.1093/jge/gxad091","RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q3","JCRName":"GEOCHEMISTRY & GEOPHYSICS","Score":null,"Total":0}
引用次数: 0

Abstract

Abstract Mapping the density contrast through the 3D gravity inversion can help detect the goals under the subsurface. However, it is a challenge to accurately and efficiently solve the 3D gravity inversion. Krylov subspace method is commonly used for large linear problems due to its high computational efficiency and low storage requirement. In this study, two classical algorithms of Krylov subspace method, namely the Generalized Minimum Residual method and the Conjugate Gradient method, are applied to 3D gravity inversion. Based on the recovered models of the deep mineral and the shallow L-shaped tunnel models, it was found that the Generalized Minimum Residual method provided similar density contrast results as the Conjugate Gradient method. The obtained inversion results of density contrast corresponded well to the position of the deep mineral resources model and the L-shaped tunnel model. The 3D distribution of Fe content underground was obtained by inverting the measured gravity data from Olympic Dam in Australia. The recovered results correspond well with the distribution of Fe content in the geological profile collected. The accuracy of inversion using the Generalized Minimum Residual method was similar to that of the Conjugate Gradient method under the same conditions. However, the Generalized Minimum Residual method had a faster convergence speed and increased inversion efficiency by about 90%, greatly reducing the inversion time and improves the inversion efficiency.
基于Krylov子空间方法的三维重力快速反演
摘要通过三维重力反演绘制密度对比图,有助于探测地下目标。然而,如何准确、高效地求解三维重力反演是一个难题。Krylov子空间方法具有计算效率高、存储容量小等优点,是求解大型线性问题的常用方法。本研究将Krylov子空间方法中的两种经典算法,即广义最小残差法和共轭梯度法应用于三维重力反演。基于深部矿物恢复模型和浅层l形隧道模型,发现广义最小残差法与共轭梯度法的密度对比结果相似。所得密度对比反演结果与深部矿产资源模型和l型隧道模型的位置吻合较好。通过对澳大利亚奥林匹克大坝实测重力数据的反演,得到了地下铁含量的三维分布。恢复结果与所收集的地质剖面中铁含量的分布吻合较好。在相同条件下,广义最小残差法反演的精度与共轭梯度法相当。而广义最小残差法收敛速度较快,反演效率提高约90%,大大缩短了反演时间,提高了反演效率。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
求助全文
约1分钟内获得全文 求助全文
来源期刊
Journal of Geophysics and Engineering
Journal of Geophysics and Engineering 工程技术-地球化学与地球物理
CiteScore
2.50
自引率
21.40%
发文量
87
审稿时长
4 months
期刊介绍: Journal of Geophysics and Engineering aims to promote research and developments in geophysics and related areas of engineering. It has a predominantly applied science and engineering focus, but solicits and accepts high-quality contributions in all earth-physics disciplines, including geodynamics, natural and controlled-source seismology, oil, gas and mineral exploration, petrophysics and reservoir geophysics. The journal covers those aspects of engineering that are closely related to geophysics, or on the targets and problems that geophysics addresses. Typically, this is engineering focused on the subsurface, particularly petroleum engineering, rock mechanics, geophysical software engineering, drilling technology, remote sensing, instrumentation and sensor design.
×
引用
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学术官方微信