Non-linear inversion using general measures of data misfit and model structure

Geophysical Journal International - Tập 134 Số 1 - Trang 213-227 - 1998
Colin G. Farquharson1, Douglas W. Oldenburg1
1UBC-Geophysical Inversion Facility, Department of Earth & Ocean Sciences, University of British Columbia, Vancouver, BC, Canada. E-mail: [email protected]

Tóm tắt

Abstract

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

10.1145/355887.355894

10.1145/355887.355895

10.1109/78.277854

10.1190/1.1443015

5E. Andersonet al.1995 LAPACK Users' Guide2nd edn, SIAM, Philadelphia, PA.

10.1190/1.1442960

Chave A.D., 1989, J. geophys. Res., 134, 215

Chave A.D., 1987, J. geophys. Res., 134, 633, 10.1029/JB092iB01p00633

10.1190/1.1440378

10.1190/1.1442303

10.1190/1.1442864

10.1190/1.1442813

13S.E. Dosso ,1990 .Inversion and appraisal for the one-dimensional magnetotellurics problem . PhD thesisUniversity of British Columbia, Vancouver, Canada.

Egbert G.D., 1986, Geophys. J. R. astr. Soc., 134, 173, 10.1111/j.1365-246X.1986.tb04552.x

Ekblom H., 1973, BIT, 134, 292, 10.1007/BF01951940

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.

Ellis R.G., 1993, J. Geomag. Geoelectr., 134, 1001, 10.5636/jgg.45.1001

Farquharson C.G., 1993, Geophys. J. Int., 134, 433, 10.1111/j.1365-246X.1993.tb06977.x

10.1190/1.1442095

20P.W. Holland, and R.E. Welsch ,1977 .Robust regression using iteratively reweighted least-squares . Communications in Statistics-Theory and MethodsA6,813 -827 .

Huber P.J., 1964, Ann. math. Stat., 134, 73, 10.1214/aoms/1177703732

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.

Larsen J.C., 1996, Geophys. J. Int., 134, 801, 10.1111/j.1365-246X.1996.tb05639.x

10.1190/1.1441261

Nabighian M.N., 1991, Electromagnetic Methods in Applied Geophysics, 134, 427, 10.1190/1.9781560802686.ch6

Oldenburg D.W., 1984, Geophys. J. R. astr. Soc., 134, 139, 10.1111/j.1365-246X.1984.tb06476.x

27D.W. Oldenburg ,1984 .An introduction to linear inverse theory . IEEE Trans. Geosci. Remote SensingGE-22,665 -674 .

10.1016/0031-9201(93)90130-2

O'Leary D.P., 1990, SIAM J. Matrix Anal. Appl., 134, 466, 10.1137/0611032

30R.L. Parker ,1994 .Geophysical Inverse Theory Princeton University Press, Princeton, NJ.

Parker R.L., 1980, J. geophys. Res., 134, 4429, 10.1029/JB085iB08p04429

Sacchi M.D., 1995, J. seis. Expl., 134, 315

10.1190/1.1444033

Scales J.A., 1988, J. comput. Phys., 134, 314, 10.1016/0021-9991(88)90115-5

10.1190/1.1442438

10.1190/1.1443012

37G.A. Watson ,1980 .Approximation Theory and Numerical Methods John Wiley, Chichester.

Wolke R., 1988, SIAM J. Sci. Stat. Comput., 134, 907, 10.1137/0909062