×

A microlocal analysis of migration. (English) Zbl 1074.74582

Summary: The mathematics of the propagation of seismic energy relevant in seismic reflection experiments is assumed to be governed by the linear acoustic wave equation, in which the coefficient is the speed of sound in the subsurface. If the speed of sound suddenly changes, the seismic signal is (partly) reflected. To find the position of these changes, one considers first the so-called forward map, which sends the coefficient of the wave equation to its solution at the surface of the earth, where the recording equipment is positioned. This map is highly nonlinear. For inversion one therefore usually takes its formal derivative, which leads to a linearized inverse problem.
It is well known that this linearized forward map corresponds in a high-frequency approximation to a Fourier integral operator. The construction of a parametrix for this Fourier integral operator requires the computation of the so-called normal operator, i.e. the composition of the linearized forward map with its adjoint. An important result by G. Beylkin [Commun. Pure Appl. Math. 37, 579–599 (1984; Zbl 0577.44003)] states that, if there are no caustics in the medium, this normal operator is an elliptic pseudo-differential operator.
In many practical situations, however, the no-caustics assumption is violated. In this paper a microlocal analysis of this more general case will be presented. We will show that for spatial dimensions less than or equal to 3, the normal operator remains a Fourier integral operator, be it not a pseudo-differential operator anymore. Instead, it is the sum of an elliptic pseudo-differential operator and a more general Fourier integral operator of lower order than the pseudo-differential part.
We will also formulate a mild injectivity condition on the traveltime function under which Beylkin’s result remains true, i.e. under which the normal operator is purely pseudo-differential. This injectivity condition includes various types of multi-valued traveltimes, which occur frequently in practice. Its geometrical interpretation will be discussed.
Finally, we derive an approximate explicit formula for the inverse, which is suitable for numerical evaluation.

MSC:

74J25 Inverse problems for waves in solid mechanics
35A27 Microlocal methods and methods of sheaf theory and homological algebra applied to PDEs
35R30 Inverse problems for PDEs
86A15 Seismology (including tsunami modeling), earthquakes

Citations:

Zbl 0577.44003
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Rakesh, A linearized inverse problem for the wave equation, Comm. PDE, 13, 573 (1988) · Zbl 0671.35078
[2] Beylkin, G., The inversion problem and applications of the generalized Radon transform, Comm. Pure and Appl. Math., 37, 579 (1984) · Zbl 0577.44003
[3] Duistermaat, J. J.; Guillemin, V., The spectrum of positive elliptic operators and periodic bicharacteristics, Inv. Math., 29, 39 (1975) · Zbl 0307.35071
[4] Hansen, S., Solution of a hyperbolic inverse problem by linearization, Comm. PDE, 16, 291 (1991) · Zbl 0737.35149
[5] Maslov, V. P., Théorie des Perturbations et des Méthodes Asymptotiques (1972), Dunod: Dunod Paris · Zbl 0247.47010
[6] Maslov, V. P.; Fedoriuk, M. V., Semi-Classical Approximation in Quantum Mechanics (1981), Reidel: Reidel Dordrecht
[7] Duistermaat, J. J., Fourier Integral Operators (1995), Birkhäuser: Birkhäuser Boston · Zbl 0272.47028
[8] Treves, F., (Introduction to Fourier Integral Operators, vol. 2 (1982), Plenum Press: Plenum Press New York)
[9] Guillemin, G., On some results of Gelfand integral geometry, (Proc. Symp. Pure Math., 43 (1985)), 149
[10] Burridge, R.; Beylkin, G., On double integrals over spheres, Inverse Problems, 4, 1 (1988) · Zbl 0652.44004
[11] ten Kroode, A. P.E.; Smit, D.-J.; Verdel, A. R., Linearized scattering in the presence of caustics, Shell Report P1295 (1995), Rijswijk
[12] Nolan, C. J.; Symes, W. W., Global solution of a linearized inverse problem for the wave equation, Comm. PDE, 22, 919 (1997) · Zbl 0889.35122
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.