Interferometry examines the general interference phenomena between pairs of signals in order to gain useful information about the subsurface. Seismic interferometry (SI) utilizes the crosscorrelation of signal pairs to reconstruct the impulse response of a given media. Jon Claerbout's initial (1968) conjecture provided the framework upon which modern theory is based.
A signal at a location A can be crosscorrelated with a signal at a location B to reproduce a virtual source-receiver pair using seismic interferometry. Crosscorrelation is often considered the key mathematical operation in this approach, but it is also possible to use convolution to come up with a similar result. The crosscorrelation of passive noise measured at a free surface reproduces the subsurface impulse response. As such, it is possible to obtain information about the subsurface with no need for an active seismic source. This method, however, is not limited to passive sources, and can be extended for use with active sources and computer–generated waveforms.
As of 2006 the field of seismic interferometry was beginning to change the way geophysicists view noise. Seismic interferometry uses this previously–ignored noise in models of the shallow subsurface. Potential applications include both research and industry.
History and Development
Claerbout (1968) first considered seismic interferometry for investigating the shallow subsurface. He later predicted that seismic interferometry could be applied to real world media, which was later proven. The long term average of random ultrasound waves can reconstruct the impulse response between two points on an aluminum block. However, they had assumed random diffuse noise, limiting interferometry in real world conditions. In a similar case, it was shown that the expressions for uncorrelated noise sources reduce to a single crosscorrelation of observations at two receivers. The interferometric impulse response of the subsurface can be reconstructed using only an extended record of background noise, initially only for the surface and direct wave arrivals.
Crosscorrelations of seismic signals from both active and passive sources at the surface or in the subsurface can be used to reconstruct a valid model of the subsurface. Seismic interferometry can produce a result similar to traditional methods without limitations on the diffusivity of the wavefield or ambient sources. In a drilling application, it is possible to utilize a virtual source to image the subsurface adjacent to a downhole location. This application is increasingly utilized particularly for exploration in subsalt settings.
Mathematical and Physical Explanation
Seismic interferometry provides for the possibility of reconstructing the subsurface reflection response using the crosscorrelations of two seismic traces. Following Claerbout's conjecture, recent work has mathematically demonstrated applications of crosscorrelation for reconstructing Green's function using wave field reciprocity theorem in a lossless, 3D heterogeneous medium. Traces are most often extended records of passive background noise, but it is also possible to utilize active sources depending on the objective. Seismic interferometry essentially exploits the phase difference between adjacent receiver locations to image the subsurface.
Seismic interferometry consists of simple crosscorrelation and stacking of actual receiver responses to approximate the impulse response as if a virtual source was placed at the location of the applicable receiver. Crosscorrelation of continuous functions in the time domain is presented as Equation 1.
(f 1∗f 2)(t)= ∫ f 1(λ)f 2(λ-t)dt
Where the functions are integrated as a function of time at different lag values. In fact, crosscorrelation can be understood conceptually as the traveltime lag associated with waveforms in two discrete receiver locations. Crosscorrelation is similar to convolution where the second function is folded relative to the first.
Seismic interferometry is fundamentally similar to the optical interferogram produced by the interference of a direct and reflected wave passing through a glass lens where intensity is primarily dependent upon the phase component.
I = 1+2R2 cos[ω(λAr+λrB)]+R^4
Where: Intensity is related to the magnitude of the reflection coefficient (R) and the phase component ω(λAr+λrB). An estimate of the reflectivity distributions can be obtained through the crosscorrelation of the direct wave at a location A with the reflection recorded at a location B where A represents the reference trace. The multiplication of the conjugate of the trace spectrum at A and the trace spectrum at B gives:
ФAB =Re^iω(λAr+λrB) + o.t.
Where: ФAB = product spectrum o.t. = additional terms, e.g. correlations of directdirect, etc. As in the previous case, the product spectrum is a function of phase.
Key: Changes in reflector geometry lead to changes in the correlation result and the reflector geometry can be recovered through the application of a migration kernel. Interpretation of raw interferograms is not normally attempted; crosscorrelated results are generally processed using some form of migration.
In the simplest case, consider a rotating drill bit at depth radiating energy that is recorded by geophones on the surface. It is possible to assume that the phase of the source wavelet at a given position is random and utilize the crosscorrelation of the direct wave at a location A with a ghost reflection at a location B to image a subsurface reflector without any knowledge regarding the source location. The crosscorrelation of traces A and B in the frequency domain simplifies as:
Ф(A, B) = −(Wiω)^2 Re^iω(λArλrB)+o.t.
Where: Wi(ω) = frequency domain source wavelet (ith wavelet)
The crosscorrelation of the direct wave at a location A with a ghost reflection at a location B removes the unknown source term where:
This form is equivalent to a virtual source configuration at a location A imaging hypothetical reflections at a location B. Migration of these correlation positions removes the phase term and yields a final migration image at position x where:
m(x) = Σø(A,B,λAx+λxB) Where: ø(A,B,t) = temporal correlation between locations A and B with lag time t
This model has been applied to simulate subsurface geometry in West Texas using simulated models including a traditional buried source and a synthetic (virtual) rotating drill bit source to produce similar results. A similar model demonstrated the reconstruction of a simulated subsurface geometry. In this case, the reconstructed subsurface response correctly modeled the relative positions of primaries and multiples. Additional equations can be derived to reconstruct signal geometries in a wide variety of cases.
Seismic interferometry is currently utilized primarily in research and academic settings. In one example, passive listening and the crosscorrelation of long noise traces was used to approximate the impulse response for shallow subsurface velocity analysis in Southern California. Seismic interferometry provided a result comparable to that indicated using elaborate inversion techniques. Seismic interferometry is most often used for the examination of the near surface and is often utilized to reconstruct surface and direct waves only. As such, seismic interferometry is commonly used to estimate ground roll to aid in its removal. Seismic interferometry simplifies estimates of shear wave velocity and attenuation in a standing building. Seismic interferometry has been applied to image the seismic scattering  and velocity structure  of volcanoes.
Exploration and production
Increasingly, seismic interferometry is finding a place in exploration and production. SI can image dipping sediments adjacent to salt domes. Complex salt geometries are poorly resolved using traditional seismic reflection techniques. An alternative method calls for the use of downhole sources and receivers adjacent to subsurface salt features. It is often difficult to generate an ideal seismic signal in a downhole location. Seismic interferometry can virtually move a source into a downhole location to better illuminate and capture steeply dipping sediments on the flank of a salt dome. In this case, the SI result was very similar to that obtained using an actual downhole source. Seismic interferometry can locate the position of an unknown source and is often utilized in hydrofrac applications to map the extent of induced fractures. It is possible that interferometric techniques can be applied to timelapse seismic monitoring of subtle changes in reservoir properties in the subsurface.
Seismic interferometry applications are currently limited by a number of factors. Real world media and noise represent limitations for current theoretical development. For example, for interferometry to work noise sources must be uncorrelated and completely surround the region of interest. In addition, attenuation and geometrical spreading are largely neglected and need to be incorporated into more robust models. Other challenges are inherent to seismic interferometry. For example, the source term only drops out in the case of the crosscorrelation of a direct wave at a location A with a ghost reflection at a location B. The correlation of other waveforms can introduce multiples to the resulting interferogram. Velocity analysis and filtering can reduce but not eliminate the occurrence of multiples in a given dataset.
Although there have been many advancements in seismic interferometry challenges still remain. One of the biggest remaining challenges is extending the theory to account for real world media and noise distributions in the subsurface. Natural sources typically do not comply with mathematical generalizations and may in fact display some degree of correlation. Additional problems must be addressed before applications of seismic interferometry can become more widespread.
- Curtis et al. 2006
- Draganov, Wapenaar & Thorbecke 2006
- Weaver & Lobkis 2001
- Wapenaar 2004
- Campillo & Paul 2003
- Schuster et al. 2004
- Bakulin & Calvert 2004
- Wapenaar & Fokkema 2006
- Yu, Followill & Schuster 2003
- Snieder & Safak 2006
- Chaput et al., 2012
- Brenguier et al. 2007
- Hornby & Yu 2007
- Lu et al. 2006
- Bakulin, A.; Calvert, R. (2004). "Virtual source: new method for imaging and 4D below complex overburden". SEG Expanded Abstracts: 24772480.
- Brenguer, F.; Shapiro, N. (2007). "3D surface wave tomography of the Piton de la Fournaise volcano using seismic noise correlations". Geophysical Research Letters. 34. Bibcode:2007GeoRL..34.2305B. doi:10.1029/2006gl028586.
- Campillo, H.; Paul, A. (2003). "Longrange correlations in the diffuse seismic coda". Science. 299: 547549. doi:10.1126/science.1078551.
- Chaput, J.; Zandomeneghi, D.; Aster, R.; Knox, H.A.; Kyle, P.R. (2012). "Imaging of Erebus volcano using body wave seismic interferometry of Strombolian eruption coda". Geophysical Research Letters. 39. Bibcode:2012GeoRL..39.7304C. doi:10.1029/2012gl050956.
- Curtis, A.; Gerstoft, P.; Sato, H.; Snieder, R.; Wapenaar, K. (2006). "Seismic interferometry turning noise into signal". The Leading Edge. 25: 10821092. doi:10.1190/1.2349814.
- Draganov, D.; Wapenaar, K.; Thorbecke, J. (2006). "Seismic interferometry: Reconstructing the earth's reflection response". Geophysics. 71: S161S170. doi:10.1190/1.2209947.
- Hornby, B.; Yu, J. (2007). "Interferometric imaging of a salt flank using walkaway VSP data". The Leading Edge. 26: 760763. doi:10.1190/1.2748493.
- Lu, R.; Willis, M.; Campman, X.; Franklin, J.; Toksoz, M. (2006). "Imaging dipping sediments at a salt dome flankVSP seismic interferometry and reversetime migration". SEG Expanded Abstracts: 21912195.
- Schuster, G.; Yu, J.; Sheng, J.; Rickett, J. (2004). "Interferometric/daylight seismic imaging". Geophysical Journal International. 157: 838852. doi:10.1111/j.1365-246x.2004.02251.x.
- Yu, J.; Followill, F.; Schuster, G. (2003). "Autocorrelogram migration of IVSPWD data: Field data test". Geophysics. 68: 297307. doi:10.1190/1.1543215.
- Weaver, R.; Lobkis, O. (2001). "Ultrasonics without a source: Thermal fluctuation correlations at MHz frequencies". Physical Review Letters. 87. Bibcode:2001PhRvL..87m4301W. doi:10.1103/PhysRevLett.87.134301.