Source location using time‐reverse imaging

We present the chain of time‐reverse modeling, image space wavefield decomposition and several imaging conditions as a migration‐like algorithm called time‐reverse imaging. The algorithm locates subsurface sources in passive seismic data and diffractors in active data. We use elastic propagators to capitalize on the full waveforms available in multicomponent data, although an acoustic example is presented as well. For the elastic case, we perform wavefield decomposition in the image domain with spatial derivatives to calculate P and S potentials. To locate sources, the time axis is collapsed by extracting the zero‐lag of auto and cross‐correlations to return images in physical space. The impulse response of the algorithm is very dependent on acquisition geometry and needs to be evaluated with point sources before processing field data. Band‐limited data processed with these techniques image the radiation pattern of the source rather than just the location. We present several imaging conditions but we imagine others could be designed to investigate specific hypotheses concerning the nature of the source mechanism. We illustrate the flexible technique with synthetic 2D passive data examples and surface acquisition geometry specifically designed to investigate tremor type signals that are not easily identified or interpreted in the time domain.