{"title":"通过有限差分频域全波形反演 GPR 数据进行多参数成像:沉积构造建模指南","authors":"Mrinal Kanti Layek, Probal Sengupta","doi":"10.1007/s00024-024-03520-1","DOIUrl":null,"url":null,"abstract":"<div><p>The need for reconstruction of the distribution of physical properties like dielectric permittivity and electrical conductivity of shallow subsurface sedimentary architecture leads to the development of an optimum strategy of GPR data inversion. In this paper, we present finite difference frequency domain (FDFD) full waveform inversion (FWI) method to get high-resolution subsurface model using GPR data. FWI is an optimization technique which involves in search of the minima between recorded and predicted data. The inversion process includes the quasi-Newton method and simultaneous frequency sampling strategy of irregular sampling. The Hessian term in quasi-Newton algorithm is approximated using preconditioned-LBFGS consideration and the search directions are also optimized after following the Wolfe conditions. At the end of each iteration during inversion, permittivity and conductivity models were updated and became ready to be the initial model for the next iteration. The goals of this research were to develop a robust framework for sedimentary-GPR data inversion and to evaluate the efficacy of the novel grid strategy introduced by Layek and Sengupta (2021) proposed for FWI. This paper presents a comparative analysis between conventional and newly proposed technique from Layek and Sengupta (2021), supported by numerical experiments conducted through our own MATLAB programming. Numerical tests conducted on a benchmark from previously published article, established the fact that new grid formulation produces a faster converging rate and required less computation time. This approach demonstrates remarkable efficacy when applied to a comprehensive sedimentary model comprising a lossy medium.</p></div>","PeriodicalId":21078,"journal":{"name":"pure and applied geophysics","volume":"181 7","pages":"2107 - 2130"},"PeriodicalIF":1.9000,"publicationDate":"2024-06-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Multi-parameter Imaging by Finite Difference Frequency Domain Full Waveform Inversion of GPR Data: A Guide for Sedimentary Architecture Modeling\",\"authors\":\"Mrinal Kanti Layek, Probal Sengupta\",\"doi\":\"10.1007/s00024-024-03520-1\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<div><p>The need for reconstruction of the distribution of physical properties like dielectric permittivity and electrical conductivity of shallow subsurface sedimentary architecture leads to the development of an optimum strategy of GPR data inversion. In this paper, we present finite difference frequency domain (FDFD) full waveform inversion (FWI) method to get high-resolution subsurface model using GPR data. FWI is an optimization technique which involves in search of the minima between recorded and predicted data. The inversion process includes the quasi-Newton method and simultaneous frequency sampling strategy of irregular sampling. The Hessian term in quasi-Newton algorithm is approximated using preconditioned-LBFGS consideration and the search directions are also optimized after following the Wolfe conditions. At the end of each iteration during inversion, permittivity and conductivity models were updated and became ready to be the initial model for the next iteration. The goals of this research were to develop a robust framework for sedimentary-GPR data inversion and to evaluate the efficacy of the novel grid strategy introduced by Layek and Sengupta (2021) proposed for FWI. This paper presents a comparative analysis between conventional and newly proposed technique from Layek and Sengupta (2021), supported by numerical experiments conducted through our own MATLAB programming. Numerical tests conducted on a benchmark from previously published article, established the fact that new grid formulation produces a faster converging rate and required less computation time. This approach demonstrates remarkable efficacy when applied to a comprehensive sedimentary model comprising a lossy medium.</p></div>\",\"PeriodicalId\":21078,\"journal\":{\"name\":\"pure and applied geophysics\",\"volume\":\"181 7\",\"pages\":\"2107 - 2130\"},\"PeriodicalIF\":1.9000,\"publicationDate\":\"2024-06-14\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"pure and applied geophysics\",\"FirstCategoryId\":\"89\",\"ListUrlMain\":\"https://link.springer.com/article/10.1007/s00024-024-03520-1\",\"RegionNum\":4,\"RegionCategory\":\"地球科学\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q2\",\"JCRName\":\"GEOCHEMISTRY & GEOPHYSICS\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"pure and applied geophysics","FirstCategoryId":"89","ListUrlMain":"https://link.springer.com/article/10.1007/s00024-024-03520-1","RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"GEOCHEMISTRY & GEOPHYSICS","Score":null,"Total":0}
Multi-parameter Imaging by Finite Difference Frequency Domain Full Waveform Inversion of GPR Data: A Guide for Sedimentary Architecture Modeling
The need for reconstruction of the distribution of physical properties like dielectric permittivity and electrical conductivity of shallow subsurface sedimentary architecture leads to the development of an optimum strategy of GPR data inversion. In this paper, we present finite difference frequency domain (FDFD) full waveform inversion (FWI) method to get high-resolution subsurface model using GPR data. FWI is an optimization technique which involves in search of the minima between recorded and predicted data. The inversion process includes the quasi-Newton method and simultaneous frequency sampling strategy of irregular sampling. The Hessian term in quasi-Newton algorithm is approximated using preconditioned-LBFGS consideration and the search directions are also optimized after following the Wolfe conditions. At the end of each iteration during inversion, permittivity and conductivity models were updated and became ready to be the initial model for the next iteration. The goals of this research were to develop a robust framework for sedimentary-GPR data inversion and to evaluate the efficacy of the novel grid strategy introduced by Layek and Sengupta (2021) proposed for FWI. This paper presents a comparative analysis between conventional and newly proposed technique from Layek and Sengupta (2021), supported by numerical experiments conducted through our own MATLAB programming. Numerical tests conducted on a benchmark from previously published article, established the fact that new grid formulation produces a faster converging rate and required less computation time. This approach demonstrates remarkable efficacy when applied to a comprehensive sedimentary model comprising a lossy medium.
期刊介绍:
pure and applied geophysics (pageoph), a continuation of the journal "Geofisica pura e applicata", publishes original scientific contributions in the fields of solid Earth, atmospheric and oceanic sciences. Regular and special issues feature thought-provoking reports on active areas of current research and state-of-the-art surveys.
Long running journal, founded in 1939 as Geofisica pura e applicata
Publishes peer-reviewed original scientific contributions and state-of-the-art surveys in solid earth and atmospheric sciences
Features thought-provoking reports on active areas of current research and is a major source for publications on tsunami research
Coverage extends to research topics in oceanic sciences
See Instructions for Authors on the right hand side.