Knowledge Resource Center for Ecological Environment in Arid Area
DOI | 10.1111/1365-2478.12548 |
Linear geophysical inversion via the discrete cosine pseudo-inverse: application to potential fields | |
Fernandez-Martinez, J. L.1; Fernandez-Muniz, Z.1; Pallero, J. L. G.1,2; Bonvalot, S.3 | |
通讯作者 | Fernandez-Martinez, J. L. |
来源期刊 | GEOPHYSICAL PROSPECTING
![]() |
ISSN | 0016-8025 |
EISSN | 1365-2478 |
出版年 | 2017 |
卷号 | 65页码:94-111 |
英文摘要 | In this paper, we present a methodology to perform geophysical inversion of large-scale linear systems via a covariance-free orthogonal transformation: the discrete cosine transform. The methodology consists of compressing the matrix of the linear system as a digital image and using the interesting properties of orthogonal transformations to define an approximation of the Moore-Penrose pseudo-inverse. This methodology is also highly scalable since the model reduction achieved by these techniques increases with the number of parameters of the linear system involved due to the high correlation needed for these parameters to accomplish very detailed forward predictions and allows for a very fast computation of the inverse problem solution. We show the application of this methodology to a simple synthetic two-dimensional gravimetric problem for different dimensionalities and different levels of white Gaussian noise and to a synthetic linear system whose system matrix has been generated via geostatistical simulation to produce a random field with a given spatial correlation. The numerical results show that the discrete cosine transform pseudo-inverse outperforms the classical least-squares techniques, mainly in the presence of noise, since the solutions that are obtained are more stable and fit the observed data with the lowest root-mean-square error. Besides, we show that model reduction is a very effective way of parameter regularisation when the conditioning of the reduced discrete cosine transform matrix is taken into account. We finally show its application to the inversion of a real gravity profile in the Atacama Desert (north Chile) obtaining very successful results in this non-linear inverse problem. The methodology presented here has a general character and can be applied to solve any linear and non-linear inverse problems (through linearisation) arising in technology and, particularly, in geophysics, independently of the geophysical model discretisation and dimensionality. Nevertheless, the results shown in this paper are better in the case of ill-conditioned inverse problems for which the matrix compression is more efficient. In that sense, a natural extension of this methodology would be its application to the set of normal equations. |
英文关键词 | Linear inversion Model reduction Discrete cosine transform (DCT) Pseudo-inverse |
类型 | Article |
语种 | 英语 |
国家 | Spain ; France |
收录类别 | SCI-E |
WOS记录号 | WOS:000418757400007 |
WOS关键词 | GRAVITY-DATA |
WOS类目 | Geochemistry & Geophysics |
WOS研究方向 | Geochemistry & Geophysics |
来源机构 | French National Research Institute for Sustainable Development |
资源类型 | 期刊论文 |
条目标识符 | http://119.78.100.177/qdio/handle/2XILL650/199292 |
作者单位 | 1.Univ Oviedo, Grp Problemas Inversos Optimizac & Aprendizaje Au, Dept Matemat, Oviedo 33007, Spain; 2.Univ Politecn Madrid, ETSI Topog Geodesia & Cartog, E-28040 Madrid, Spain; 3.Univ Toulouse, Lab GET, CNRS, IRD,CNES,BGI, F-31401 Toulouse, France |
推荐引用方式 GB/T 7714 | Fernandez-Martinez, J. L.,Fernandez-Muniz, Z.,Pallero, J. L. G.,et al. Linear geophysical inversion via the discrete cosine pseudo-inverse: application to potential fields[J]. French National Research Institute for Sustainable Development,2017,65:94-111. |
APA | Fernandez-Martinez, J. L.,Fernandez-Muniz, Z.,Pallero, J. L. G.,&Bonvalot, S..(2017).Linear geophysical inversion via the discrete cosine pseudo-inverse: application to potential fields.GEOPHYSICAL PROSPECTING,65,94-111. |
MLA | Fernandez-Martinez, J. L.,et al."Linear geophysical inversion via the discrete cosine pseudo-inverse: application to potential fields".GEOPHYSICAL PROSPECTING 65(2017):94-111. |
条目包含的文件 | 条目无相关文件。 |
除非特别说明,本系统中所有内容都受版权保护,并保留所有权利。