Path-integral seismic imaging

A new type of seismic imaging, based on Feynman path integrals for waveform modelling, is capable of producing accurate subsurface images without any need for a reference velocity model. Instead of the usual optimization for traveltime curves with maximal signal semblance, a weighted summation over all representative curves avoids the need for velocity analysis, with its common difficulties of subjective and time-consuming manual picking. The summation over all curves includes the stationary one that plays a preferential role in classical imaging schemes, but also multiple stationary curves when they exist. Moreover, the weighted summation over all curves also accounts for non-uniqueness and uncertainty in the stacking/migration velocities. The path-integral imaging can be applied to stacking to zero-offset and to time and depth migration. In all these cases, a properly defined weighting function plays a vital role: To emphasize contributions from traveltime curves close to the optimal one and to suppress contributions from unrealistic curves. The path-integral method is an authentic macromodel-independent technique in the sense that there is strictly no parameter optimization or estimation involved. Development is still in its initial stage, and several conceptual and implementation issues are yet to be solved. However, application to synthetic and real data examples shows that it has the potential for becoming a fully automatic imaging technique.