{"title":"Simultaneously recovering both domain and varying density in inverse gravimetry by efficient level-set methods","authors":"Wenbin Li, J. Qian","doi":"10.3934/ipi.2020073","DOIUrl":null,"url":null,"abstract":"We develop new efficient algorithms for a class of inverse problems of gravimetry to recover an anomalous volume mass distribution (measure) in the sense that we design fast local level-set methods to simultaneously reconstruct both unknown domain and varying density of the anomalous measure from modulus of gravity force rather than from gravity force itself. The equivalent-source principle of gravitational potential forces us to consider only measures of the form \\begin{document}$ \\mu = f\\,\\chi_{D} $\\end{document} , where \\begin{document}$ f $\\end{document} is a density function and \\begin{document}$ D $\\end{document} is a domain inside a closed set in \\begin{document}$ \\bf{R}^n $\\end{document} . Accordingly, various constraints are imposed upon both the density function and the domain so that well-posedness theories can be developed for the corresponding inverse problems, such as the domain inverse problem, the density inverse problem, and the domain-density inverse problem. Starting from uniqueness theorems for the domain-density inverse problem, we derive a new gradient from the misfit functional to enforce the directional-independence constraint of the density function and we further introduce a new labeling function into the level-set method to enforce the geometrical constraint of the corresponding domain; consequently, we are able to recover simultaneously both unknown domain and varying density from given modulus of gravity force. Our fast level-set method is built upon localizing the level-set evolution around a narrow band near the zero level-set and upon accelerating numerical modeling by novel low-rank matrix multiplication. Numerical results demonstrate that uniqueness theorems are crucial for solving the inverse problem of gravimetry and will be impactful on gravity prospecting. To the best of our knowledge, our inversion algorithm is the first of such for the domain-density inverse problem since it is based upon the conditional well-posedness theory of the inverse problem.","PeriodicalId":50274,"journal":{"name":"Inverse Problems and Imaging","volume":"37 1","pages":""},"PeriodicalIF":1.2000,"publicationDate":"2021-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"4","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Inverse Problems and Imaging","FirstCategoryId":"100","ListUrlMain":"https://doi.org/10.3934/ipi.2020073","RegionNum":4,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
引用次数: 4
Abstract
We develop new efficient algorithms for a class of inverse problems of gravimetry to recover an anomalous volume mass distribution (measure) in the sense that we design fast local level-set methods to simultaneously reconstruct both unknown domain and varying density of the anomalous measure from modulus of gravity force rather than from gravity force itself. The equivalent-source principle of gravitational potential forces us to consider only measures of the form \begin{document}$ \mu = f\,\chi_{D} $\end{document} , where \begin{document}$ f $\end{document} is a density function and \begin{document}$ D $\end{document} is a domain inside a closed set in \begin{document}$ \bf{R}^n $\end{document} . Accordingly, various constraints are imposed upon both the density function and the domain so that well-posedness theories can be developed for the corresponding inverse problems, such as the domain inverse problem, the density inverse problem, and the domain-density inverse problem. Starting from uniqueness theorems for the domain-density inverse problem, we derive a new gradient from the misfit functional to enforce the directional-independence constraint of the density function and we further introduce a new labeling function into the level-set method to enforce the geometrical constraint of the corresponding domain; consequently, we are able to recover simultaneously both unknown domain and varying density from given modulus of gravity force. Our fast level-set method is built upon localizing the level-set evolution around a narrow band near the zero level-set and upon accelerating numerical modeling by novel low-rank matrix multiplication. Numerical results demonstrate that uniqueness theorems are crucial for solving the inverse problem of gravimetry and will be impactful on gravity prospecting. To the best of our knowledge, our inversion algorithm is the first of such for the domain-density inverse problem since it is based upon the conditional well-posedness theory of the inverse problem.
We develop new efficient algorithms for a class of inverse problems of gravimetry to recover an anomalous volume mass distribution (measure) in the sense that we design fast local level-set methods to simultaneously reconstruct both unknown domain and varying density of the anomalous measure from modulus of gravity force rather than from gravity force itself. The equivalent-source principle of gravitational potential forces us to consider only measures of the form \begin{document}$ \mu = f\,\chi_{D} $\end{document} , where \begin{document}$ f $\end{document} is a density function and \begin{document}$ D $\end{document} is a domain inside a closed set in \begin{document}$ \bf{R}^n $\end{document} . Accordingly, various constraints are imposed upon both the density function and the domain so that well-posedness theories can be developed for the corresponding inverse problems, such as the domain inverse problem, the density inverse problem, and the domain-density inverse problem. Starting from uniqueness theorems for the domain-density inverse problem, we derive a new gradient from the misfit functional to enforce the directional-independence constraint of the density function and we further introduce a new labeling function into the level-set method to enforce the geometrical constraint of the corresponding domain; consequently, we are able to recover simultaneously both unknown domain and varying density from given modulus of gravity force. Our fast level-set method is built upon localizing the level-set evolution around a narrow band near the zero level-set and upon accelerating numerical modeling by novel low-rank matrix multiplication. Numerical results demonstrate that uniqueness theorems are crucial for solving the inverse problem of gravimetry and will be impactful on gravity prospecting. To the best of our knowledge, our inversion algorithm is the first of such for the domain-density inverse problem since it is based upon the conditional well-posedness theory of the inverse problem.
期刊介绍:
Inverse Problems and Imaging publishes research articles of the highest quality that employ innovative mathematical and modeling techniques to study inverse and imaging problems arising in engineering and other sciences. Every published paper has a strong mathematical orientation employing methods from such areas as control theory, discrete mathematics, differential geometry, harmonic analysis, functional analysis, integral geometry, mathematical physics, numerical analysis, optimization, partial differential equations, and stochastic and statistical methods. The field of applications includes medical and other imaging, nondestructive testing, geophysical prospection and remote sensing as well as image analysis and image processing.
This journal is committed to recording important new results in its field and will maintain the highest standards of innovation and quality. To be published in this journal, a paper must be correct, novel, nontrivial and of interest to a substantial number of researchers and readers.