Yu Zhang et al 2003 Inverse Problems 19 1113 doi:10.1088/0266-5611/19/5/307
Yu Zhang1, Guanquan Zhang2 and Norman Bleistein3
Show affiliationsOne-way wave operators are powerful tools for use in forward modelling and inversion. Their implementation, however, involves introduction of the square root of an operator as a pseudo-differential operator. Furthermore, a simple factoring of the wave operator produces one-way wave equations that yield the same travel times as the full wave equation, but do not yield accurate amplitudes except for homogeneous media and for almost all points in heterogeneous media. Here, we present augmented one-way wave equations. We show that these equations yield solutions for which the leading order asymptotic amplitude as well as the travel time satisfy the same differential equations as the corresponding functions for the full wave equation. Exact representations of the square-root operator appearing in these differential equations are elusive, except in cases in which the heterogeneity of the medium is independent of the transverse spatial variables. Here, we address the fully heterogeneous case. Singling out depth as the preferred direction of propagation, we introduce a representation of the square-root operator as an integral in which a rational function of the transverse Laplacian appears in the integrand. This allows us to carry out explicit asymptotic analysis of the resulting one-way wave equations. To do this, we introduce an auxiliary function that satisfies a lower dimensional wave equation in transverse spatial variables only. We prove that ray theory for these one-way wave equations leads to one-way eikonal equations and the correct leading order transport equation for the full wave equation. We then introduce appropriate boundary conditions at z = 0 to generate waves at depth whose quotient leads to a reflector map and an estimate of the ray theoretical reflection coefficient on the reflector. Thus, these true amplitude one-way wave equations lead to a 'true amplitude wave equation migration' (WEM) method. In fact, we prove that applying the WEM imaging condition to these newly defined wavefields in heterogeneous media leads to the Kirchhoff inversion formula for common-shot data when the one-way wavefields are replaced by their ray theoretic approximations. This extension enhances the original WEM method. The objective of that technique was a reflector map, only. The underlying theory did not address amplitude issues. Computer output obtained using numerically generated data confirms the accuracy of this inversion method. However, there are practical limitations. The observed data must be a solution of the wave equation. Therefore, the data over the entire survey area must be collected from a single common-shot experiment. Multi-experiment data, such as common-offset data, cannot be used with this method as currently formulated. Research on extending the method is ongoing at this time.
03.65.Ge Solutions of wave equations: bound states
02.60.Lj Ordinary and partial differential equations; boundary value problems
35S30 Fourier integral operators
35R30 Inverse problems (undetermined coefficients, etc.) for PDE
81Q20 Semiclassical techniques including WKB and Maslov methods
Issue 5 (October 2003)
Received 24 April 2003
Published 5 September 2003
Yu Zhang et al 2003 Inverse Problems 19 1113
M A Rodríguez-Valverde et al 2003 Eur. J. Phys. 24 159
K B Sundaram and G K Bhagavat 1981 J. Phys. D: Appl. Phys. 14 333
Sheng-Jui Chen et al 2006 J. Phys.: Conf. Ser. 32 244
John Robertson 2006 Rep. Prog. Phys. 69 327
Arnd Bäcker and Roman Schubert 2002 J. Phys. A: Math. Gen. 35 539
Peng-Fei Hao et al 2006 J. Micromech. Microeng. 16 1397
J M Aguirregabiria et al 2004 Eur. J. Phys. 25 555
K F Freed 1979 J. Phys. C: Solid State Phys. 12 L17
S Baillet et al 2001 Phys. Med. Biol. 46 77