Refraction surveys are a well-established method of imaging subsurface velocities, both in terms of the deep crustal structure at global scales and in the shallow near surface. These surveys generally involve deploying an array of receivers on the surface (or water bottom) and recording arrivals from a seismic source initiated at or near the surface.

S U M M A R YInverting for the subsurface velocity distribution by refraction traveltime tomography is a wellaccepted imaging method by both the exploration and earthquake seismology communities. A significant drawback, however, is that the recorded traces become noisier with increasing offset from the source position, and so accurate picking of traveltimes in far-offset traces is often prevented. To enhance the signal-to-noise ratio (SNR) of the far-offset traces, we present the theory of supervirtual refraction interferometry where the SNR of far-offset head-wave arrivals can be theoretically increased by a factor proportional to √ N ; here, N is the number of receiver or source positions associated with the recording and generation of the head-wave arrival. There are two steps to this methodology: correlation and summation of the data to generate traces with virtual head-wave arrivals, followed by the convolution of the data with the virtual traces to create traces with supervirtual head-wave arrivals. This method is valid for any medium that generates head-wave arrivals recorded by the geophones. Results with both synthetic traces and field data demonstrate the feasibility of this method.There are at least four significant benefits of supervirtual interferometry: (1) an enhanced SNR of far-offset traces so the first-arrival traveltimes of the noisy far-offset traces can be more reliably picked to extend the useful aperture of the data, (2) the SNR of head waves in a trace that arrive later than the first arrival can be enhanced for accurate traveltime picking and subsequent inversion by later-arrival traveltime tomography, (3) common receiver-pair gathers can be analysed to detect the presence of diving waves in the first arrivals, which can be used to assess the nature of the refracting boundary, and (4) the source statics term is eliminated in the correlation operations so that the timing of the virtual traces is independent of the source excitation time. This suggests the possibility of applying this method to earthquake data recorded by receivers that are inline with the refraction paths and source locations.

Inverting for the subsurface velocity distribution by refraction traveltime tomography is a well-accepted imaging method by both the exploration and earthquake seismology communities. A significant drawback, however, is that the recorded traces become noisier with increasing offset from the source position, and so prevents accurate picking of traveltimes in faroffset traces. To enhance the signal-to-noise ratio of the faroffset traces, we present the theory of super-virtual refraction interferometry where the signal-to-noise ratio (SNR) of faroffset head-wave arrivals can be theoretically increased by a factor proportional to √ N; here, N is the number of receiver and source positions associated with the recording and generation of the head-wave arrival. There are two steps to this methodology: correlation and summation of the data to generate traces with virtual head-wave arrivals, followed by the convolution of the data with the virtual traces to create traces with super-virtual head-wave arrivals. This method is valid for any medium that generates head-wave arrivals. There are at least three significant benefits to this methodology: 1). enhanced SNR of far-offset traces so the first-arrival traveltimes of the noisy far-offset traces can be more reliably picked to extend the useful aperture of data, 2). the SNR of head waves in a trace that arrive after the first arrival can be enhanced for accurate traveltime picking and subsequent inversion by traveltime tomography, and 3). common receiver-pair gathers can be analyzed to detect the presence of diving waves in the first arrivals, which can be used to assess the nature of the refracting boundary.