Time-reversal imaging is a powerful method for the localization of microseismic events. Conventional time-reversal imaging methods such as autocorrelation imaging or grouped crosscorrelation imaging may suffer from imaging artifacts, e.g., caused by less refined velocity models, noise-contaminated data, and data acquired from sparse receiver networks. These artifacts typically reduce imaging quality and may cause subsequent misinterpretations leading to false positives. To address these issues, we develop a new imaging condition that comprises three steps for each time step. First, we divide the back-propagated wavefield into parts according to their maximum absolute amplitudes. Second, the amplitudes of the back-propagated wavefield are weighted by a Gaussian function with the spatial extent of the prevailing wavelength of the event, centered at the absolute maximum of that part of the wavefield. Finally, we zero-lag crosscorrelate these weighted wavefields at each space point to obtain an image for this time step. The final image is gained by summing the images for each time step. This process collects all energy concentrations along the back projection process, and the energy on the wavefront overlaps and collapses at the hypocenter leading to high-resolution images displaying little to no imaging artifacts. Numerical examples using the Marmousi-II and the 3D SEG overthrust models and a 3D field data example indicate the performance of our method. High-resolution low-noise source images allow unique identification of sources even for source clusters, noisy data, and sparse acquisitions. The source localization errors are smaller than the dominant wavelength of the signal, where a smooth model with a mean velocity error of approximately 5% was considered in the synthetic examples, and a homogeneous model was used in the field data example.
Read full abstract