Note, the crosscorrelation image of the point scatterer is smeared over a larger depth range than that of the Kirchhoff image. This is because the crosscorrelation of one trace with another smears the source wavelet into a longer wavelet, and also because the crosscorrelation operator has poor resolution in the depth direction (see Figure 5). Nevertheless, the crosscorrelation point-scatterer image is acceptable.
In practice, the master trace and its two nearest neighbors were muted because the direct wave migration operator in equation (4) has zero or nearly zero phase when . This is undesirable because any energy from these traces will be smeared uniformly throughout the model, not just at the exploder points. Also, a second derivative in time was applied to the crosscorrelogram traces to partly compensate for the smoothing effects of crosscorrelation and migration.
Figure 5 is the same as Figure 4 except the source wavelet is a long random time series. The crosscorrelation of traces collapses the ringy time series to an impulse-like wavelet so that the associated migration image in the lower panel of Figure 5 has good spatial resolution compared to the Kirchhoff image in the middle panel.
In the previous examples, the scatterer exploded at time zero. Figure 6 shows a similar synthetic experiment with ten scatterers, all exploding at random times with a random time series as a source wavelet. The resulting data for 1 second is shown in the top panel of Figure 6. The middle panel shows these data after crosscorrelation migration of 1 second of data, and roughly locates the location of the 10 point sources. Repeating this crosscorelation migration for fifteen data sets, each with 1 second of data generated from ten point scatterers with distinct random time histories, yields the stacked images shown in the lower panel of Figure 6. As expected, averaging the migration images tends to cancel migration noise and reinforce the energy at the location of the point exploders.
Finally, the fault-like structure denoted by stars in Figure 7 is assumed to emanate seismic energy randomly in time with random strength. This might approximate the situation where fluid is injected along a reservoir bed and seismic instruments are passively monitoring the injection front. Figure 7 shows the results after crosscorrelation migration of 1 second of data (middle panel), and 40 stacks of 1 second records (lower panel). The fault boundaries are much better delineated in the 40-stack migration image, although the resolution is much worse than that of an ordinary seismic survey.
Poor resolution of the crosscorrelation images is consistent with theoretical predictions of resolution for crosscorrelation migration operators. A possibility for improving resolution is to measure the incidence angle of energy in the correlograms and use this angle as a constraint in smearing data into the model. This strategy is similar to that of ray-map migration, but it remains to be seen if this is a practical strategy with correlograms.