Non-linear inversion using general measures of data misfit and model structure
Tóm tắt
We investigate the use of general, non-l 2 measures of data misfit and model structure in the solution of the non-linear inverse problem. Of particular interest are robust measures of data misfit, and measures of model structure which enable piecewise-constant models to be constructed. General measures can be incorporated into traditional linearized, iterative solutions to the non-linear problem through the use of an iteratively reweighted least-squares (IRLS) algorithm. We show how such an algorithm can be used to solve the linear inverse problem when general measures of misfit and structure are considered. The magnetic stripe example of Parker (1994) is used as an illustration. This example also emphasizes the benefits of using a robust measure of misfit when outliers are present in the data. We then show how the IRLS algorithm can be used within a linearized, iterative solution to the non-linear problem. The relevant procedure contains two iterative loops which can be combined in a number of ways. We present two possibilities. The first involves a line search to determine the most appropriate value of the trade-off parameter and the complete solution, via the IRLS algorithm, of the linearized inverse problem for each value of the trade-off parameter. In the second approach, a schedule of prescribed values for the trade-off parameter is used and the iterations required by the IRLS algorithm are combined with those for the linearized, iterative inversion procedure. These two variations are then applied to the 1-D inversion of both synthetic and field time-domain electromagnetic data.
Từ khóa
Tài liệu tham khảo
5E. Andersonet al.1995 LAPACK Users' Guide2nd edn, SIAM, Philadelphia, PA.
Chave A.D., 1989, J. geophys. Res., 134, 215
13S.E. Dosso ,1990 .Inversion and appraisal for the one-dimensional magnetotellurics problem . PhD thesisUniversity of British Columbia, Vancouver, Canada.
16H. Ekblom ,1987 .TheL1-estimate as limiting case of anLp- or Huber-estimate, in . Statistical Data Analysis Based on the L1-Norm and Related Methods pp. 109 -116 , ed. Dodge, Y., Elsevier, Amsterdam.
20P.W. Holland, and R.E. Welsch ,1977 .Robust regression using iteratively reweighted least-squares . Communications in Statistics-Theory and MethodsA6,813 -827 .
LaD. Brecque, G. Morelli, B. Daily, A. Ramirez, and P. Lundegard ,1997 .Occam's inversion of 3D ERT data, in . Three-Dimensional Electromagnetics eds Oristaglio, M.L. & Spies, B.R., Soc. Expl. Geophysicists, Tulsa, OK, in press.
Nabighian M.N., 1991, Electromagnetic Methods in Applied Geophysics, 134, 427, 10.1190/1.9781560802686.ch6
27D.W. Oldenburg ,1984 .An introduction to linear inverse theory . IEEE Trans. Geosci. Remote SensingGE-22,665 -674 .
30R.L. Parker ,1994 .Geophysical Inverse Theory Princeton University Press, Princeton, NJ.
Sacchi M.D., 1995, J. seis. Expl., 134, 315
37G.A. Watson ,1980 .Approximation Theory and Numerical Methods John Wiley, Chichester.