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
Clemens Elster et al 2007 Metrologia 44 L31
C J Lundy and R E Olson 1996 J. Phys. B: At. Mol. Opt. Phys. 29 1723
Igor V Meglinski and Stephen J Matcher 2002 Physiol. Meas. 23 741
Archana Pai et al 2002 Class. Quantum Grav. 19 1477
J Schurr et al 2009 Metrologia 46 619
S Carlip 1998 Class. Quantum Grav. 15 2629
B Hohler and H Schreyer 1982 J. Phys. F: Met. Phys. 12 857
R Stadler et al 2003 Nanotechnology 14 138
D T A Blair et al 1989 J. Phys. D: Appl. Phys. 22 755