-
- Mark Pilkington
- 1Geological Survey of Canada, Ottawa, Ontario, Canada. E-mail: mpilking@nrcan.gc.ca.
書誌事項
- 公開日
- 2008-12-12
- DOI
-
- 10.1190/1.3026538
- 公開者
- Society of Exploration Geophysicists
この論文をさがす
説明
<jats:title>Abstract</jats:title> <jats:p>I have developed an inversion approach that determines a 3D susceptibility distribution that produces a given magnetic anomaly. The subsurface model consists of a 3D, equally spaced array of dipoles. The inversion incorporates a model norm that enforces sparseness and depth weighting of the solution. Sparseness is imposed by using the Cauchy norm on model parameters. The inverse problem is posed in the data space, leading to a linear system of equations with dimensions based on the number of data, N. This contrasts with the standard least-squares solution, derived through operations within the M-dimensional model space (M being the number of model parameters). Hence, the data-space method combined with a conjugate gradient algorithm leads to computational efficiency by dealing with an N×N system versus an M×M one, where N≪M. Tests on synthetic data show that sparse inversion produces a much more focused solution compared with a standard model-space, least-squares inversion. The inversion of aeromagnetic data collected over a Precambrian Shield area again shows that including the sparseness constraint leads to a simpler and better resolved solution. The degree of improvement in model resolution for the sparse case is quantified using the resolution matrix.</jats:p>
収録刊行物
-
- Geophysics
-
Geophysics 74 (1), L7-L15, 2008-12-12
Society of Exploration Geophysicists
