Document Type : Original Research Paper

Authors

1 Institute of Geophysics, Unversity of Tehran, Iran

2 Institute of Geophysics, Unversity of Tehran, Iran.

Abstract

        In experimental sciences we often need to solve inverse problems. That is, we want to obtain information about the internal structure of a physical system from indirect noisy observations. Information about the errors in the observations is essential to solve any inverse problem; otherwise it is impossible to say when a feature ‘fits the data’. In practice, however, one seldom has a direct estimate of the data errors. Here, we exploit the trade-off between data prediction and model or data structure to determine model based estimates of the noise characteristics from a single realization of the data. Noise estimates are then used to characterize the set of reasonable models that fit the data. By intersecting set of prior model parameter constraints with the set of data fitting models, we obtain a set of models that fit the data and are in agreement with prior constraints. This prior information can also be used to set bounds on the bias. We illustrate our methods with a synthetic example of vertical seismic profiling (VSP).

Keywords

References
Gouveia, W. and Scales, J. A., 1998- Bayesian seismic waveform inversion: parameter estimation and uncertainty analysis: J. Geophys. Res., 103, 2759–2779.
Green, P. J. and Silverman, B. W., 1994 - Nonparametric regression and generalized linear models: A roughness penalty approach, Chapman and Hall, London.
Lawson, C. L. and Hanson, R., J., 1974 - Solving least squares problems, Prentice Hall, Englewood Clifts, NJ.
Li, Y. and Oldenburg, D. W., 1999- 3-D inversion of DC resistivity data using an L-curve criterion, in Annual Meeting Abstracts, pp. 251–254, Society of Exploration Geophysicists.
Scales, J. A., 1987- Tomographic inversion via the conjugate gradient method: Geophysics, 52, 179–185.
Scales, J. A., Smith, M., L., 2001- Introductory geophysical inverse theory, Samizdat Press, Colorado School of Mine, Denver.
Scales, J. A. and Snieder, R. K., 1998- What is noise?: Geophysics, 63, 1122–1124.
Scales, J. A. and Tenorio, L., 2001- Prior information and uncertainty in inverse problems: Geophysics, 66, 389–397.
Stark, P., 1992- Inference in infinite dimensional inverse problems: Discretization and duality: J. Geophys. Res., 97, 14055–14082.
Tikhonov, A. N., and Arsenin, V. Y., 1977- Solutions of Ill-posed problems,Winston, John Wiley and Sons, New York.
Trampert, J. and Snieder, R. K., 1996- Model estimations biased by truncated expansion: Possible artifacts in seismic tomography: Science, 271, 1257–1260.
Van Wijk, K., Scales, J. A., Navidi, W. and Roy-Chowdhury, K., 1998- Estimating Data uncertainties for least squares optimization, in Annual Project Review, Vol. CWP283, Center for Wave Phenomena, Colorado School of Mines, Golden.
Van Wijk, K., Scales, J. A., Navidi, W. and Tenorio, L., 2002- data and model uncertainty estimation for linear inversion: Geophys. J. Int., 149, 625–632
Xu, P., 1998- Svd methods for linear ill-posed problems: Geophys. J. Int., 135, 505–514.