Fast solution of geophysical inversion using adaptive mesh, space-filling curves and wavelet compression

Abstract

S U M M A R YMany geophysical inverse problems involve large and dense coefficient matrices that often exceed the limitations of physical memory in commonly available computers. The repeated multiplications of such matrices to vectors during processing or inversion require an immense amount of computing power. These two factors pose a significant challenge to solving largescale inverse problems in practice and can render many realistic problems intractable. To overcome these limitations, we develop a new computational approach for this class of problems by combining an adaptive quadtree or octree model discretization and wavelet transforms on reordered parameter sets. The adaptive mesh discretizes the model region according to the required resolutions based on localized anomalies. Hilbert space-filling curves and similar ordering of the reduced parameter set then enable a higher compression of the coefficient matrix by forming its sparse representation in the 1-D wavelet domain. This combination can reduce the storage requirement by 100 to 1000 times and, therefore, also speeds up the computation during the processing stage by the same factor. As a result, problems can now be solved that were computationally prohibitive. We present the algorithm and illustrate its effectiveness with an example from equivalent source construction in potential-field processing.

Toward improved prediction of the bedrock depth underneath hillslopes: Bayesian inference of the bottom‐up control hypothesis using high‐resolution topographic data

“…Alternatively, one can use so‐called model‐reduction techniques and let the actual regolith depth observations determine simultaneously the pattern and values of the DTB model parameters. Examples of such model‐reduction approaches include the discrete cosine transform [ Linde and Vrugt , ; Lochbühler et al ., ], wavelet transform [ Davis and Li , ; Jafarpour , ], and singular value decomposition [ Laloy et al ., ; Oware et al ., ]. We have tested this alternative approach in the present study but found little improvements in the quality of fit of the DTB model (not shown).…”

A new stable downward continuation of airborne magnetic data based on Wavelet deconvolution

“…Low‐pass filtering suppresses the noise but blurs the signal and defeats the enhancement of multi‐source anomalies (Daubechies et al . ; Davis and Li ; Moreau et al . ; Trompat et al .…”

Markov chain Monte Carlo simulation using the DREAM software package: Theory, concepts, and MATLAB implementation

“…This allows for the use of standard closed-form prior distributions for the reduced set of parameters. Examples of such approaches include the discrete cosine transform (Jafarpour et al, 2009(Jafarpour et al, , 2010Linde and Vrugt, 2013;Lochbühler et al, 2015), wavelet transform (Davis and Li, 2011;Jafarpour, 2011), and singular value decomposition Oware et al, 2013).…”