ILL-POSED AND WELL-POSED FORMULATIONS OF THE REFLECTION TRAVEL TIME TOMOGRAPHY PROBLEM

被引:43
作者
DELPRATJANNAUD, F
LAILLY, P
机构
关键词
D O I
10.1029/92JB02441
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
Many studies have shown the strong dependence of the solution of reflection tomography (or travel time inversion) in the model discretization interval: such solutions are purely numerical artifacts. The cause lies in the formulation of the reflection tomography which turns out to be an ill-posed problem unless properly regularized. Adequate regularizations are those which make use of derivatives of the functions describing the model (velocity distributions and reflector geometries). The order of the derivatives to use depends on the space dimension (two or three dimension). Of practical interest is curvature regularization: it works both in two and three dimensions with minimum a priori knowledge on the substratum under survey. We give theoretical arguments to prove that, with such a regularization, the computed solution is unique, stable, and almost independent of model discretization, provided that it is fine enough. These theoretical results are confirmed by numerical tests on a two-dimensional model with a velocity field varying both vertically and horizontally: they show the practical interest of the proposed regularization. However, mathematical stability is not sufficient, and a priori information bar, to be used in order to decrease the size of the set of possible solutions. As a consequence of our theoretical study of regularization techniques, only a priori information involving derivatives of the functions describing the model yields stability. In other words, except for some pathological examples, stability is obtained only if we try to determine a smooth model of the subsurface. These conclusions remain valid for three-dimensional models and for transmission tomography. As a by-product of our study, we suggest a preconditioning which makes much easier the solution of the linear systems involved in the Gauss-Newton optimization method.
引用
收藏
页码:6589 / 6605
页数:17
相关论文
共 26 条
  • [1] [Anonymous], 1987, INVERSE PROBLEM THEO
  • [2] UNIQUENESS IN INVERSION OF INACCURATE GROSS EARTH DATA
    BACKUS, G
    GILBERT, F
    [J]. PHILOSOPHICAL TRANSACTIONS OF THE ROYAL SOCIETY OF LONDON SERIES A-MATHEMATICAL AND PHYSICAL SCIENCES, 1970, 266 (1173) : 123 - &
  • [3] RESOLVING POWER OF GROSS EARTH DATA
    BACKUS, G
    GILBERT, F
    [J]. GEOPHYSICAL JOURNAL OF THE ROYAL ASTRONOMICAL SOCIETY, 1968, 16 (02): : 169 - &
  • [4] Bartels R. H., 1986, INTRO SPLINES USE CO
  • [5] VELOCITY-DEPTH AMBIGUITY OF REFLECTION TRAVEL-TIMES
    BICKEL, SH
    [J]. GEOPHYSICS, 1990, 55 (03) : 266 - 276
  • [6] TOMOGRAPHIC DETERMINATION OF VELOCITY AND DEPTH IN LATERALLY VARYING MEDIA
    BISHOP, TN
    BUBE, KP
    CUTLER, RT
    LANGAN, RT
    LOVE, PL
    RESNICK, JR
    SHUEY, RT
    SPINDLER, DA
    WYLD, HW
    [J]. GEOPHYSICS, 1985, 50 (06) : 903 - 923
  • [7] Bois P., 1971, Geophysical Prospecting, V19, P42, DOI 10.1111/j.1365-2478.1971.tb00585.x
  • [8] Boor CD., 1978, PRACTICAL GUIDE SPLI
  • [9] Ciarlet PG, 1982, INTRO ANAL NUMERIQUE
  • [10] WHAT INFORMATION ON THE EARTH MODEL DO REFLECTION TRAVEL-TIMES PROVIDE
    DELPRATJANNAUD, F
    LAILLY, P
    [J]. JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH, 1992, 97 (B13) : 19827 - 19844