Acoustic crosswell imaging using asymptotic waveforms

Seismic waveforms are inverted using an asymptotic method. The asymptotic m ethod models amplitudes correctly at the caustics and takes nonstationary r aypaths into account when computing the waveforms, and thus is an extension of geometrical ray theory. Using numerical differencing, partial derivativ es of the data with respect to the model are computed. As expected, these p artial derivatives (or sensitivity functions) are concentrated along, but n ot confined to, raypaths. The sensitivity functions enable the formulation of a waveform inversion algorithm, which is applied to a synthetic crosswel l experiment and a laboratory crosswell experiment. The synthetic experimen t shows the advantages of the waveform inversion method over conventional t raveltime inversion methods. Boundaries of anomalies are better defined, an d smearing is reduced. The waveform inversion produces a much lower misfit than the traveltime inversion. The goal of the laboratory experiment was th e detection of a nonaqueous phase liquid (NAPL) in water saturated sand. Th e sand was imaged before and after injection of the NAPL. Using the wavefor m inversion method, low-velocity anomalies were imaged that correlate well with post-experiment determination of NAPL concentrations. The low-velocity anomaly defocuses the seismic energy. However, the amplitude reduction due to the low-velocity anomaly is not enough to explain the observed low ampl itudes. We suggest that other mechanisms (such as multiple scattering, 3-D effects, or intrinsic attenuation) not included in the asymptotic waveform modeling play an important role in decreasing the amplitude.