Source-receiver wave field interferometry

Next-generation seismic model of the Australian crust from synchronous and asynchronous ambient noise imaging

The proliferation of seismic networks in Australia has laid the groundwork for high-resolution probing of the continental crust. Here we develop an updated 3D shear-velocity model using a large dataset containing nearly 30 years of seismic recordings from over 1600 stations. A recently-developed ambient noise imaging workflow enables improved data analysis by integrating asynchronous arrays across the continent. This model reveals fine-scale crustal structures at a lateral resolution of approximately 1-degree in most parts of the continent, highlighted by 1) shallow low velocities (<3.2 km/s) well correlated with the locations of known sedimentary basins, 2) consistently faster velocities beneath discovered mineral deposits, suggesting a whole-crustal control on the mineral deposition process, and 3) distinctive crustal layering and improved characterization of depth and sharpness of the crust-mantle transition. Our model sheds light on undercover mineral exploration and inspires future multi-disciplinary studies for a more comprehensive understanding of the mineral systems in Australia.

Wave-field representations with Green's functions, propagator matrices, and Marchenko-type focusing functions

Classical acoustic wave-field representations consist of volume and boundary integrals, of which the integrands contain specific combinations of Green's functions, source distributions, and wave fields. Using a unified matrix-vector wave equation for different wave phenomena, these representations can be reformulated in terms of Green's matrices, source vectors, and wave-field vectors. The matrix-vector formalism also allows the formulation of representations in which propagator matrices replace the Green's matrices. These propagator matrices, in turn, can be expressed in terms of Marchenko-type focusing functions. An advantage of the representations with propagator matrices and focusing functions is that the boundary integrals in these representations are limited to a single open boundary. This makes these representations a suitable basis for developing advanced inverse scattering, imaging and monitoring methods for wave fields acquired on a single boundary.

Seismic Interferometry from Correlated Noise Sources

It is a well-established principle that cross-correlating seismic observations at different receiver locations can yield estimates of band-limited inter-receiver Green’s functions. This principle, known as Green’s function retrieval or seismic interferometry, is a powerful technique that can transform noise into signals which enable remote interrogation and imaging of the Earth’s subsurface. In practice it is often necessary and even desirable to rely on noise already present in the environment. Theory that underpins many applications of ambient noise interferometry assumes that the sources of noise are uncorrelated in time. However, many real-world noise sources such as trains, highway traffic and ocean waves are inherently correlated in space and time, in direct contradiction to the these theoretical foundations. Applying standard interferometric techniques to recordings from correlated energy sources makes the Green’s function liable to estimation errors that so far have not been fully accounted for theoretically nor in practice. We show that these errors are significant for common noise sources, always perturbing or entirely obscuring the phase one wishes to retrieve. Our analysis explains why stacking may reduce the phase errors, but also shows that in commonly encountered circumstances stacking will not remediate the problem. This analytical insight allowed us to develop a novel workflow that significantly mitigates effects arising from the use of correlated noise sources. Our methodology can be used in conjunction with already existing approaches, and improves results from both correlated and uncorrelated ambient noise. Hence, we expect it to be widely applicable in ambient noise studies.

Automatic identification of multiply diffracted waves and their ordered scattering paths

An automated algorithm uses recordings of acoustic energy across a spatially-distributed array to derive information about multiply scattered acoustic waves in heterogeneous media. The arrival time and scattering-order of each recorded diffracted acoustic wave, and the exact sequence of diffractors encountered by that wave, are estimated without requiring an explicit model of the medium through which the wave propagated. Individual diffractors are identified on the basis of their unique single-scattering relative travel-time curves (move-outs) across the array, and secondary (twice-scattered) waves are detected using semblance analysis along temporally offset primary move-outs. This information is sufficient to estimate travel times and scattering paths of all multiply diffracted waves of any order, and for these events to be identified in recorded data. The algorithm is applied to synthetic acoustic data sets from a variety of media, including different numbers of point-diffractors and a medium with strong heterogeneity and non-hyperbolic move-outs.

Physical and non-physical energy in scattered wave source-receiver interferometry

Source-receiver interferometry allows Green's functions between sources and receivers to be estimated by means of convolution and cross-correlation of other wavefields. Source-receiver interferometry has been observed to work surprisingly well in practical applications when theoretical requirements (e.g., complete enclosing boundaries of other sources and receivers) are contravened: this paper contributes to explain why this may be true. Commonly used inter-receiver interferometry requires wavefields to be generated around specific stationary points in space which are controlled purely by medium heterogeneity and receiver locations. By contrast, application of source-receiver interferometry constructs at least kinematic information about physically scattered waves between a source and a receiver by cross-convolution of scattered waves propagating from and to any points on the boundary. This reduces the ambiguity in interpreting wavefields generated using source-receiver interferometry with only partial boundaries (as is standard in practical applications), as it allows spurious or non-physical energy in the constructed Green's function to be identified and ignored. Further, source-receiver interferometry (which includes a step of inter-receiver interferometry) turns all types of non-physical or spurious energy deriving from inter-receiver interferometry into what appears to be physical energy. This explains in part why source-receiver interferometry may perform relatively well compared to inter-receiver interferometry when constructing scattered wavefields.

Source-receiver interferometry for seismic wavefield construction and ground-roll removal

Seismic interferometry describes the construction of unmeasured wavefield responses (or Green's functions) between two or more points by applying cross-correlation, deconvolution, or convolution to seismic data recordings. The practical implications are that, applying inter-receiver interferometry, a “virtual” (imaginary) source of energy can be created at the location to a real receiver by using energy recorded from surrounding sources. Similarly, by using inter-source interferometry, a virtual receiver can be created at the location of a real source by using energy recorded at surrounding receivers (Curtis et al., 2009). These two methods can be combined to create the new technique of source-receiver interferometry (Curtis, 2009; Curtis and Halliday, 2010; Halliday and Curtis, 2010) which synthesizes real-source to real-receiver Green's function estimates using only energy recorded at a surrounding boundary of receivers and from an additional surrounding boundary of sources. The boundary sources in each case can be active (such as an explosive or vibrating source) or passive (such as noise from anthropogenic activity or ocean waves). This paper describes the first real-data application of source-receiver interferometry, and demonstrates its potential to enhance existing methods of interferometric ground-roll suppression.