Algorithms and Software for Large-Scale Geophysical Reconstructions

Christian Eske Bruun, Trine Brandt Nielsen

AbstractThe main focus of the thesis is algorithms and software for large-scale geophysical reconstructions. Throughout the project we deal with a special Toeplitz structure of the coeffecient matrix that enables a significant loss-less compression.

The geophysical surveying problems dealt with in the thesis are by nature ill-posed. For this reason we use regularization methods to perform the reconstuctions. The methods used in the paper are TSVD, Tikhonov, and CGLS. We will describe the constraints in the surveying problems that need to be met in order to achieve the Toeplitz structures. Aside from enabling compression the Toeplitz structures makes it possible to use a FFT matrix-vector multiplication to achieve fast multiplications in the regularization methods. Multi-level observation sets are used to meet the constraints and a Fourier based upward continuation method can be used to achieve the multi-levels, however, for forward problems we are able to compute levels at different altitudes. The performance of the upward continuation method is investigated.

The result of this thesis is a Matlab object-oriented package implemented within the frameworks of MOOReTools and expands the existing GravMagTools package. The package is tested and reconstructions are performed.
KeywordsToeplitz structures, Large-scale inversions, Ill-posed problems, Compression, FFT matrix-vector multiplication, Regularization methods, Upward continuation, Object-oriented implementation, MOOReTools
TypeMaster's thesis [Academic thesis]
Year2007
PublisherInformatics and Mathematical Modelling, Technical University of Denmark, DTU
AddressRichard Petersens Plads, Building 321, DK-2800 Kgs. Lyngby
SeriesIMM-Thesis-2007-29
NoteSupervised by Professor Per Christian Hansen, IMM, DTU.
Electronic version(s)[pdf]
BibTeX data [bibtex]
IMM Group(s)Scientific Computing