Multiresolution parameterization for geophysical inverse problems
- 1 January 2003
- journal article
- Published by Society of Exploration Geophysicists in Geophysics
- Vol. 68 (1), 199-209
- https://doi.org/10.1190/1.1543207
Abstract
A model parameterization based on a multiresolution wavelet representation is proposed for geophysical inverse problems in this study. The matrix equation for the wavelet representation of the model is constructed by dual wavelet transforms of a readily built matrix based on the pixel parameterization. We demonstrate the merits of the proposed parameterization using the simple examples of geophysical data gridding and 3D seismic tomography, and comparing the results with those obtained by the conventional pixel parameterization. It is shown that for the conventional damped least‐squares solutions of pixel‐parameterized models, regularization by norm damping deteriorates the underlying model correlation due to sparse and heterogeneous sampling, whereas roughness damping schemes impose a priori arbitrary correlation length irrespective of the heterogeneous sampling. The scale hierarchy embedded in the proposed multiresolution parameterization facilitates the compromise between the dual spatial and scale resolution. Depending on the local richness of the data constraints across different scales at a site, model variation of various scale spectra can be robustly resolved through the scale hierarchy. There is thus no need to invoke additional smoothness regularization.Keywords
This publication has 34 references indexed in Scilit:
- High‐ and low‐resolution images of the Earth's mantle: Implications of different approaches to tomographic modelingJournal of Geophysical Research, 1999
- Closing the gap between regional and global travel time tomographyJournal of Geophysical Research, 1998
- Resolution of seismic waveform inversion: Bayes versus OccamInverse Problems, 1997
- A high‐resolution P wave model for the top 1200 km of the mantleJournal of Geophysical Research, 1996
- Model Estimations Biased by Truncated Expansions: Possible Artifacts in Seismic TomographyScience, 1996
- Are our parameter estimators biased? The significance of finite-difference regularization operatorsInverse Problems, 1995
- A theory for multiresolution signal decomposition: the wavelet representationIeee Transactions On Pattern Analysis and Machine Intelligence, 1989
- Biharmonic spline interpolation of GEOS‐3 and SEASAT altimeter dataGeophysical Research Letters, 1987
- Decomposition of Hardy Functions into Square Integrable Wavelets of Constant ShapeSIAM Journal on Mathematical Analysis, 1984
- A REVIEW OF LEAST‐SQUARES INVERSION AND ITS APPLICATION TO GEOPHYSICAL PROBLEMS*Geophysical Prospecting, 1984