· · · Institutional Repository

It has been shown by many authors that the cross correlation of two recordings of a diffuse wave field at different receivers yields the Green’s function between these receivers. Recently the theory has been extended for situations where time-reversal invariance does not hold (e.g., in attenuating media) and where source-receiver reciprocity breaks down (in moving fluids). Here we present a unified theory for Green’s function retrieval that captures all these situations and, because of the unified form, readily extends to more complex situations, such as electrokinetic Green’s function retrieval in poroelastic or piezoelectric media. The unified theory has a wide range of applications in ‘‘remote sensing without a source.’’

We show that the two-point cross-correlation of self-potential field recordings is equal to the electric resistivity between the two points. This holds under the condition that spatially and temporally uncorrelated noise sources exist throughout the volume. These sources should have a known amplitude spectrum and their correlated strengths should be proportional to the dissipative medium property function. Natural fluctuations, such as thermal noise, may occur that satisfy the necessary conditions. When these fluctuations are random deviations from a state of thermal equilibrium, the fluctuation-dissipation theorem can be used to describe these sources. Other types of sources may exist, such as the ones creating the self-potential field through coupling with fluctuations in pressure, temperature and chemical potential gradients.

We compare two approaches for deriving the fact that the Green’s function in an arbitrary inhomogeneous open system can be obtained by cross correlating recordings of the wave field at two positions. One approach is based on physical arguments, exploiting the principle of time-reversal invariance of the acoustic wave equation. The other approach is based on Rayleigh’s reciprocity theorem. Using a unified notation, we show that the result of the time-reversal approach can be obtained as an approximation of the result of the reciprocity approach.

With seismic interferometry one can retrieve the response to a virtual source inside an unknown medium, if there is a receiver at the position of the virtual source. Using inverse scattering theory, we demonstrate that, for a 1D medium, the requirement of having an actual receiver inside the medium can be circumvented, going beyond seismic interferometry. In this case, the wavefield can be focused inside an unknown medium with independent variations in velocity and density using reflection data only.

Turning noise into useful data—every geophysicist's dream? And now it seems possible. The field of seismic interferometry has at its foundation a shift in the way we think about the parts of the signal that are currently filtered out of most analyses—complicated seismic codas (the multiply scattered parts of seismic waveforms) and background noise (whatever is recorded when no identifiable active source is emitting, and which is superimposed on all recorded data). Those parts of seismograms consist of waves that reflect and refract around exactly the same subsurface heterogeneities as waves excited by active sources. The key to the rapid emergence of this field of research is our new understanding of how to unravel that subsurface information from these relatively complex-looking waveforms. And the answer turned out to be rather simple. This article explains the operation of seismic interferometry and provides a few examples of its application.

Seismic interferometry is a technique for estimating the Green's function that accounts for wave propagation between receivers by correlating the waves recorded at these receivers. We present a derivation of this principle based on the method of stationary phase. Although this derivation is intended to be educational, applicable to simple media only, it provides insight into the physical principle of seismic interferometry. In a homogeneous medium with one horizontal reflector and without a free surface, the correlation of the waves recorded at two receivers correctly gives both the direct wave and the singly reflected waves. When more reflectors are present, a product of the singly reflected waves occurs in the crosscorrelation that leads to spurious multiples when the waves are excited at the surface only. We give a heuristic argument that these spurious multiples disappear when sources below the reflectors are included. We also extend the derivation to a smoothly varying heterogeneous background medium.

Interferometry by multidimensional deconvolution applied to Controlled-Source Electromagnetic data replaces the medium above the receivers by a homogeneous halfspace, suppresses the direct field and redatums the source positions to the receiver locations. In that sense, the airwave and any other interactions of the signal with the air-water interface and the water layer are suppressed and the source uncertainty is reduced. Interferometry requires grid data and cannot be applied to line data unless the source is infinitely long in the crossline direction. To create such a source, a set of source lines is required. We use an iterative algorithm to determine the optimal locations of these source lines and show that more source lines are required if the source is towed closer to the sea bottom and closer to the receivers.

The one-dimensional Marchenko equation forms the basis for inverse scattering problems in which the scattering object is accessible from one side only. Here we derive a three-dimensional (3D) Marchenko equation which relates the single-sided reflection response of a 3D inhomogeneous medium to a field inside the medium. We show that this equation is solved by a 3D iterative data-driven focusing method, which yields the 3D Green’s function with its virtual source inside the medium. The 3D single-sided Marchenko equation and its iterative solution method form the basis for imaging of 3D strongly scattering inhomogeneous media that are accessible from one side only.

It has recently been shown that the response to a virtual source in the subsurface can be derived from reflection data at the surface and an estimate of the direct arrivals between the virtual source and the surface. Hence, unlike for seismic interferometry, no receivers are needed inside the medium. This new method recovers the complete wavefield of a virtual source, including all internal multiple scattering. Because no actual receivers are needed in the medium, the virtual source can be placed anywhere in the subsurface. With some additional processing steps (decomposition and multidimensional deconvolution) it is possible to obtain a redatumed reflection response at any depth level in the subsurface, from which all the overburden effects are eliminated. By applying standard migration between these depth levels, a true amplitude image of the subsurface can be obtained, free from ghosts due to internal multiples. The method is non-recursive and therefore does not suffer from error propagation. Moreover, the internal multiples are eliminated by deconvolution, hence no adaptive prediction and subtraction is required.

Controlled‐source electromagnetics (CSEM) has been used as a de‐risking tool in the hydrocarbon exploration industry. Although there have been successful applications of CSEM, this technique is still not widely used in the industry because the limited types of hydrocarbon reservoirs CSEM can detect. In this paper, we apply the concept of synthetic aperture to CSEM data. Synthetic aperture allows us to design sources with specific radiation patterns for different purposes. The ability to detect reservoirs is dramatically increased after forming an appropriate synthetic aperture antenna. Consequently, the types of hydrocarbon reservoirs that CSEM can detect are significantly extended. Because synthetic apertures are constructed as a data processing step, there is no additional cost for the CSEM acquisition. Synthetic aperture has potential for simplifying and reducing the cost of CSEM acquisition. We show a data example that illustrates the increased sensitivity obtained by applying synthetic aperture CSEM source.

We have analyzed the far-field approximation of the Green's function representation for seismic interferometry. By writing each of the Green's functions involved in the correlation process as a superposition of a direct wave and a scattered wave, the Green's function representation is rewritten as a superposition of four terms. When the scattered waves are modeled with the Born approximation, it appears that a three-term approximation of the Green's function representation (omitting the term containing the crosscorrelation of the scattered waves) yields a nearly exact retrieval, whereas the full four-term expression leads to a significant nonphysical event. This is because the Born approximation does not conserve energy and therefore is an insufficient model to explain all aspects of seismic interferometry. We use the full four-term expression of the Green's function representation to derive the generalized optical theorem. Unlike other recent derivations, which use stationary phase analysis, our derivation uses reciprocity theory. From the generalized optical theorem, we derive the nonlinear scattering matrix of a point scatterer. This nonlinear model accounts for primary and multiple scattering at the point scatterer and conforms with well-established scattering theory of classical waves. The model is essential to explain fully the results of seismic interferometry, even when it is applied to the response of a single point scatterer. The nonlinear scattering matrix also has implications for modeling, inversion, and migration.

In the 1990s, the method of time-reversed acoustics was developed. This method exploits the fact that the acoustic wave equation for a lossless medium is invariant for time reversal. When ultrasonic responses recorded by piezoelectric transducers are reversed in time and fed simultaneously as source signals to the transducers, they focus at the position of the original source, even when the medium is very complex. In seismic interferometry the time-reversed responses are not physically sent into the earth, but they are convolved with other measured responses. The effect is essentially the same: The time-reversed signals focus and create a virtual source which radiates waves into the medium that are subsequently recorded by receivers. A mathematical derivation, based on reciprocity theory, formalizes this principle: The crosscorrelation of responses at two receivers, integrated over differ-ent sources, gives the Green's function emitted by a virtual source at the position of one of the receivers and observed by the other receiver. This Green's function representation for seismic interferometry is based on the assumption that the medium is lossless and nonmoving. Recent developments, circumventing these assumptions, include interferometric representations for attenuating and/or moving media, as well as unified representations for waves and diffusion phenomena, bending waves, quantum mechanical scattering, potential fields, elastodynamic, electromagnetic, poroelastic, and electroseismic waves. Significant improvements in the quality of the retrieved Green's functions have been obtained with interferometry by deconvolution. A trace-by-trace deconvolution process compensates for complex source functions and the attenuation of the medium. Interferometry by multidimensional deconvolution also compensates for the effects of one-sided and/or irregular illumination.

Controlled-source electromagnetics (CSEM) has been used as a derisking tool in the hydrocarbon exploration industry. We apply the concept of synthetic aperture to the lowfrequency electromagnetic field in CSEM. Synthetic aperture sources have been used in radar imaging for many years. Using the synthetic aperture concept, big synthetic sources can be constructed by adding the response to small sources (building blocks) in different ways, and consequently, big sources with different radiation patterns can be created. We show that the detectability of hydrocarbons is significantly enhanced by applying synthetic aperture to CSEM data. More challenging targets such as deep reservoirs (4km belowsea floor) can be detected. The synthetic aperture technique also increases the sensitivity of the field to subsurface targets in the towing streamer acquisition.We also show that a pseudovertical source (orthogonally distributed dipole pairs) can be constructed synthetically, and that the detection capability of this pseudovertical source is increased by applying field steering. The synthetic aperture concept opens a new line of research in CSEM, with the freedom to design suitable synthetic aperture sources for a given purpose.

With interferometry applied to controlled-source electromagnetic data, the direct field and the airwave and all other effects related to the air-water interface can be suppressed in a data-driven way. Interferometry allows for retreival of the scattered field Green’s function of the subsurface or, in other words, the subsurface reflection response. This reflection response can then be further used to invert for the subsurface conductivity distribution. To perform interferometry in 3D, measurements on an areal grid are necessary. We discuss 3D interferometry by multidimensional deconvolution in the frequency-wavenumber and in the frequency-space domains and provide examples for a layered earth model. We use the synthetic aperture source concept to damp the signal at high wavenumbers to allow large receiver sampling distances. Interferometry indeed increases the detectability of a subsurface reservoir. Finally, we discuss the dependency of the accuracy of the retrieved reflection response on the two crucial parameters: the conductivity of the seabed at the receiver location and the stabilization parameter of the least-squares inversion.

Seismic interferometry involves the crosscorrelation of responses at different receivers to obtain the Green's function between these receivers. For the simple situation of an impulsive plane wave propagating along the x-axis, the crosscorrelation of the responses at two receivers along the x-axis gives the Green's function of the direct wave between these receivers. When the source function of the plane wave is a transient (as in exploration seismology) or a noise signal (as in passive seismology), then the crosscorrelation gives the Green's function, convolved with the autocorrelation of the source function. Direct-wave interferometry also holds for 2D and 3D situations, assuming the receivers are surrounded by a uniform distribution of sources. In this case, the main contributions to the retrieved direct wave between the receivers come from sources in Fresnel zones around stationary points. The main application of direct-wave interferometry is theretrieval of seismic surface-wave responses from ambient noise and the subsequent tomographic determination of the surface-wave velocity distribution of the subsurface. Seismic interferometry is not restricted to retrieving direct waves between receivers. In a classic paper, Claerbout shows that the autocorrelation of the transmission response of a layered medium gives the plane-wave reflection response of that medium. This is essentially 1D reflected-wave interferometry. Similarly, the crosscorrelation of the transmission responses, observed at two receivers, of an arbitrary inhomogeneous medium gives the 3D reflection response of that medium. One of the main applications of reflected-wave interferometry is retrieving the seismic reflection response from ambient noise and imaging of the reflectors in the subsurface. A common aspect of direct- and reflected-wave interferometry is that virtual sources are created at positions where there are only receivers without requiring knowledge of the subsurface medium parameters or of the positions of the actual sources.