Menu
Home Explore People Places Arts History Plants & Animals Science Life & Culture Technology
On this page
Seismic interferometry

Seismic interferometry utilizes the crosscorrelation of signal pairs to reconstruct the impulse response of the subsurface, enabling virtual source-receiver pairs without active sources. Originating from foundational work by Keiiti Aki, Géza Kunetz, and Jon Claerbout, this technique uses ambient seismic noise—once ignored—to extract Green’s functions that reveal subsurface properties. Besides crosscorrelation, convolution can also achieve similar results. Applying this method transforms the understanding of geophysics by framing subsurface imaging as an inverse problem. Its applications span from large-scale continental studies to monitoring natural hazards, industry, and environmental concerns, all without necessarily requiring active seismic sources, leveraging the powerful information contained in seismic noise.

Related Image Collections Add Image
We don't have any YouTube videos related to Seismic interferometry yet.
We don't have any PDF documents related to Seismic interferometry yet.
We don't have any Books related to Seismic interferometry yet.
We don't have any archived web articles related to Seismic interferometry yet.

History and development

Claerbout (1968) developed a workflow to apply existing interferometry techniques to investigating the shallow subsurface, although it was not proven until later that seismic interferometry could be applied to real world media.89 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.1011

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.12 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.13

Mathematical and Physical Explanation

Seismic interferometry provides for the possibility of reconstructing the subsurface reflection response using the crosscorrelations of two seismic traces.1415 Recent work16 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.

The conditions for the method to be valid, meaning to retrieve the Green's function from correlated signals, are given as follows:1718

  • sources are uncorrelated in time,
  • sources are located all around the receivers to reconstruct surface waves,
  • the wavefield is equipartitioned, meaning it comprises both compressional and shear waves.

The last two conditions are hard to meet directly in nature. However, thanks to the wave scattering, the waves are converted, which satisfies the equipartition condition. The equal distribution of sources is met thanks to the fact, that the waves are scattered in every direction.19

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.20 Crosscorrelation of continuous functions in the time domain is presented as Equation 1.

Equation 1

( f 1 ∗ f 2 ) ( t ) = ∫ f 1 ( λ ) f 2 ( λ − t ) d λ {\displaystyle (f_{1}*f_{2})(t)=\int f_{1}(\lambda )f_{2}(\lambda -t)d\lambda }

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.21

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.

Equation 2

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).2223 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.24 The multiplication of the conjugate of the trace spectrum at A and the trace spectrum at B gives:

Equation 3

ФAB =Re^iω(λAr+λrB) + o.t.

Where: ФAB = product spectrum o.t. = additional terms, e.g. correlations of direct-direct, 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.2526 Interpretation of raw interferograms is not normally attempted; crosscorrelated results are generally processed using some form of migration.27

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.28 The crosscorrelation of traces A and B in the frequency domain simplifies as:

Equation 4

Ф(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:

Equation 5

Ф(A,B)≈Re^iω(λArλrB)

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.2930 A similar model demonstrated the reconstruction of a simulated subsurface geometry.31 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.

Applications

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.32 Seismic interferometry simplifies estimates of shear wave velocity and attenuation in a standing building.33 Seismic interferometry has been applied to image the seismic scattering 34 and velocity structure 35 of volcanoes.

Exploration and production

Increasingly, seismic interferometry is finding a place in exploration and production.36 SI can image dipping sediments adjacent to salt domes.37 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.3839 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.40 It is possible that interferometric techniques can be applied to timelapse seismic monitoring of subtle changes in reservoir properties in the subsurface.41

Limitations

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.42 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.43

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.44 Additional problems must be addressed before applications of seismic interferometry can become more widespread.

Notes

References

  1. Curtis et al. 2006 - Curtis, A.; Gerstoft, P.; Sato, H.; Snieder, R.; Wapenaar, K. (2006). "Seismic interferometry turning noise into signal". The Leading Edge. 25 (9): 1082–1092. Bibcode:2006LeaEd..25.1082C. doi:10.1190/1.2349814. http://resolver.tudelft.nl/uuid:bc099609-47bd-4ad8-91ad-2a4a9a949c0f

  2. Aki, Keiiti (1957). "Space and time spectra of stationary stochastic waves, with special reference to microtremors". Bulletin of the Earthquake Research Institute. 35: 415–457. hdl:2261/11892. /wiki/Hdl_(identifier)

  3. Claerbout, Jon F. (April 1968). "Synthesis of a Layered Medium from ITS Acoustic Transmission Response". Geophysics. 33 (2): 264–269. Bibcode:1968Geop...33..264C. doi:10.1190/1.1439927. /wiki/Bibcode_(identifier)

  4. Snieder, Roel; Wapenaar, Kees (2010-09-01). "Imaging with ambient noise". Physics Today. 63 (9): 44–49. Bibcode:2010PhT....63i..44S. doi:10.1063/1.3490500. ISSN 0031-9228. https://physicstoday.scitation.org/doi/10.1063/1.3490500

  5. Draganov, Wapenaar & Thorbecke 2006 - Draganov, D.; Wapenaar, K.; Thorbecke, J. (2006). "Seismic interferometry: Reconstructing the earth's reflection response". Geophysics. 71 (4): SI61 – SI70. Bibcode:2006Geop...71SI.61D. CiteSeerX 10.1.1.75.113. doi:10.1190/1.2209947. https://ui.adsabs.harvard.edu/abs/2006Geop...71SI.61D

  6. Curtis et al. 2006 - Curtis, A.; Gerstoft, P.; Sato, H.; Snieder, R.; Wapenaar, K. (2006). "Seismic interferometry turning noise into signal". The Leading Edge. 25 (9): 1082–1092. Bibcode:2006LeaEd..25.1082C. doi:10.1190/1.2349814. http://resolver.tudelft.nl/uuid:bc099609-47bd-4ad8-91ad-2a4a9a949c0f

  7. Curtis et al. 2006 - Curtis, A.; Gerstoft, P.; Sato, H.; Snieder, R.; Wapenaar, K. (2006). "Seismic interferometry turning noise into signal". The Leading Edge. 25 (9): 1082–1092. Bibcode:2006LeaEd..25.1082C. doi:10.1190/1.2349814. http://resolver.tudelft.nl/uuid:bc099609-47bd-4ad8-91ad-2a4a9a949c0f

  8. Curtis et al. 2006 - Curtis, A.; Gerstoft, P.; Sato, H.; Snieder, R.; Wapenaar, K. (2006). "Seismic interferometry turning noise into signal". The Leading Edge. 25 (9): 1082–1092. Bibcode:2006LeaEd..25.1082C. doi:10.1190/1.2349814. http://resolver.tudelft.nl/uuid:bc099609-47bd-4ad8-91ad-2a4a9a949c0f

  9. Weaver & Lobkis 2001 - Weaver, R.; Lobkis, O. (2001). "Ultrasonics without a source: Thermal fluctuation correlations at MHz frequencies". Physical Review Letters. 87 (13): 134301. Bibcode:2001PhRvL..87m4301W. doi:10.1103/PhysRevLett.87.134301. PMID 11580591. https://ui.adsabs.harvard.edu/abs/2001PhRvL..87m4301W

  10. Wapenaar 2004 harvnb error: no target: CITEREFWapenaar2004 (help)

  11. Campillo & Paul 2003 - Campillo, H.; Paul, A. (2003). "Longrange correlations in the diffuse seismic coda". Science. 299 (5606): 547–549. Bibcode:2003Sci...299..547C. doi:10.1126/science.1078551. PMID 12543969. S2CID 22021516. https://ui.adsabs.harvard.edu/abs/2003Sci...299..547C

  12. Schuster et al. 2004 - Schuster, G.; Yu, J.; Sheng, J.; Rickett, J. (2004). "Interferometric/daylight seismic imaging". Geophysical Journal International. 157 (2): 838–852. Bibcode:2004GeoJI.157..838S. doi:10.1111/j.1365-246x.2004.02251.x. https://doi.org/10.1111%2Fj.1365-246x.2004.02251.x

  13. Bakulin & Calvert 2004 - Bakulin, A.; Calvert, R. (2004). "Virtual source: new method for imaging and 4D below complex overburden". SEG Expanded Abstracts: 2477–2480. doi:10.1190/1.1845233. https://doi.org/10.1190/1.1845233

  14. Curtis et al. 2006 - Curtis, A.; Gerstoft, P.; Sato, H.; Snieder, R.; Wapenaar, K. (2006). "Seismic interferometry turning noise into signal". The Leading Edge. 25 (9): 1082–1092. Bibcode:2006LeaEd..25.1082C. doi:10.1190/1.2349814. http://resolver.tudelft.nl/uuid:bc099609-47bd-4ad8-91ad-2a4a9a949c0f

  15. Draganov, Wapenaar & Thorbecke 2006 - Draganov, D.; Wapenaar, K.; Thorbecke, J. (2006). "Seismic interferometry: Reconstructing the earth's reflection response". Geophysics. 71 (4): SI61 – SI70. Bibcode:2006Geop...71SI.61D. CiteSeerX 10.1.1.75.113. doi:10.1190/1.2209947. https://ui.adsabs.harvard.edu/abs/2006Geop...71SI.61D

  16. Wapenaar & Fokkema 2006 harvnb error: no target: CITEREFWapenaarFokkema2006 (help)

  17. Curtis et al. 2006 - Curtis, A.; Gerstoft, P.; Sato, H.; Snieder, R.; Wapenaar, K. (2006). "Seismic interferometry turning noise into signal". The Leading Edge. 25 (9): 1082–1092. Bibcode:2006LeaEd..25.1082C. doi:10.1190/1.2349814. http://resolver.tudelft.nl/uuid:bc099609-47bd-4ad8-91ad-2a4a9a949c0f

  18. Larose, Eric; Carrière, Simon; Voisin, Christophe; Bottelin, Pierre; Baillet, Laurent; Guéguen, Philippe; Walter, Fabian; Jongmans, Denis; Guillier, Bertrand; Garambois, Stéphane; Gimbert, Florent (2015-05-01). "Environmental seismology: What can we learn on earth surface processes with ambient noise?". Journal of Applied Geophysics. 116: 62–74. Bibcode:2015JAG...116...62L. doi:10.1016/j.jappgeo.2015.02.001. ISSN 0926-9851. https://www.sciencedirect.com/science/article/pii/S0926985115000403

  19. Larose, Eric; Carrière, Simon; Voisin, Christophe; Bottelin, Pierre; Baillet, Laurent; Guéguen, Philippe; Walter, Fabian; Jongmans, Denis; Guillier, Bertrand; Garambois, Stéphane; Gimbert, Florent (2015-05-01). "Environmental seismology: What can we learn on earth surface processes with ambient noise?". Journal of Applied Geophysics. 116: 62–74. Bibcode:2015JAG...116...62L. doi:10.1016/j.jappgeo.2015.02.001. ISSN 0926-9851. https://www.sciencedirect.com/science/article/pii/S0926985115000403

  20. Curtis et al. 2006 - Curtis, A.; Gerstoft, P.; Sato, H.; Snieder, R.; Wapenaar, K. (2006). "Seismic interferometry turning noise into signal". The Leading Edge. 25 (9): 1082–1092. Bibcode:2006LeaEd..25.1082C. doi:10.1190/1.2349814. http://resolver.tudelft.nl/uuid:bc099609-47bd-4ad8-91ad-2a4a9a949c0f

  21. Animation http://mathworld.wolfram.com/Convolution.html

  22. Draganov, Wapenaar & Thorbecke 2006 - Draganov, D.; Wapenaar, K.; Thorbecke, J. (2006). "Seismic interferometry: Reconstructing the earth's reflection response". Geophysics. 71 (4): SI61 – SI70. Bibcode:2006Geop...71SI.61D. CiteSeerX 10.1.1.75.113. doi:10.1190/1.2209947. https://ui.adsabs.harvard.edu/abs/2006Geop...71SI.61D

  23. Wapenaar & Fokkema 2006 harvnb error: no target: CITEREFWapenaarFokkema2006 (help)

  24. Schuster et al. 2004 - Schuster, G.; Yu, J.; Sheng, J.; Rickett, J. (2004). "Interferometric/daylight seismic imaging". Geophysical Journal International. 157 (2): 838–852. Bibcode:2004GeoJI.157..838S. doi:10.1111/j.1365-246x.2004.02251.x. https://doi.org/10.1111%2Fj.1365-246x.2004.02251.x

  25. Curtis et al. 2006 - Curtis, A.; Gerstoft, P.; Sato, H.; Snieder, R.; Wapenaar, K. (2006). "Seismic interferometry turning noise into signal". The Leading Edge. 25 (9): 1082–1092. Bibcode:2006LeaEd..25.1082C. doi:10.1190/1.2349814. http://resolver.tudelft.nl/uuid:bc099609-47bd-4ad8-91ad-2a4a9a949c0f

  26. Schuster et al. 2004 - Schuster, G.; Yu, J.; Sheng, J.; Rickett, J. (2004). "Interferometric/daylight seismic imaging". Geophysical Journal International. 157 (2): 838–852. Bibcode:2004GeoJI.157..838S. doi:10.1111/j.1365-246x.2004.02251.x. https://doi.org/10.1111%2Fj.1365-246x.2004.02251.x

  27. Schuster et al. 2004 - Schuster, G.; Yu, J.; Sheng, J.; Rickett, J. (2004). "Interferometric/daylight seismic imaging". Geophysical Journal International. 157 (2): 838–852. Bibcode:2004GeoJI.157..838S. doi:10.1111/j.1365-246x.2004.02251.x. https://doi.org/10.1111%2Fj.1365-246x.2004.02251.x

  28. Schuster et al. 2004 - Schuster, G.; Yu, J.; Sheng, J.; Rickett, J. (2004). "Interferometric/daylight seismic imaging". Geophysical Journal International. 157 (2): 838–852. Bibcode:2004GeoJI.157..838S. doi:10.1111/j.1365-246x.2004.02251.x. https://doi.org/10.1111%2Fj.1365-246x.2004.02251.x

  29. Schuster et al. 2004 - Schuster, G.; Yu, J.; Sheng, J.; Rickett, J. (2004). "Interferometric/daylight seismic imaging". Geophysical Journal International. 157 (2): 838–852. Bibcode:2004GeoJI.157..838S. doi:10.1111/j.1365-246x.2004.02251.x. https://doi.org/10.1111%2Fj.1365-246x.2004.02251.x

  30. Yu, Followill & Schuster 2003 - Yu, J.; Followill, F.; Schuster, G. (2003). "Autocorrelogram migration of IVSPWD data: Field data test". Geophysics. 68 (1): 297–307. Bibcode:2003Geop...68..297Y. doi:10.1190/1.1543215. https://ui.adsabs.harvard.edu/abs/2003Geop...68..297Y

  31. Draganov, Wapenaar & Thorbecke 2006 - Draganov, D.; Wapenaar, K.; Thorbecke, J. (2006). "Seismic interferometry: Reconstructing the earth's reflection response". Geophysics. 71 (4): SI61 – SI70. Bibcode:2006Geop...71SI.61D. CiteSeerX 10.1.1.75.113. doi:10.1190/1.2209947. https://ui.adsabs.harvard.edu/abs/2006Geop...71SI.61D

  32. Curtis et al. 2006 - Curtis, A.; Gerstoft, P.; Sato, H.; Snieder, R.; Wapenaar, K. (2006). "Seismic interferometry turning noise into signal". The Leading Edge. 25 (9): 1082–1092. Bibcode:2006LeaEd..25.1082C. doi:10.1190/1.2349814. http://resolver.tudelft.nl/uuid:bc099609-47bd-4ad8-91ad-2a4a9a949c0f

  33. Snieder & Safak 2006 harvnb error: no target: CITEREFSniederSafak2006 (help)

  34. Chaput et al., 2012

  35. Brenguier et al. 2007

  36. Hornby & Yu 2007 - Hornby, B.; Yu, J. (2007). "Interferometric imaging of a salt flank using walkaway VSP data". The Leading Edge. 26 (6): 760–763. Bibcode:2007LeaEd..26..760H. doi:10.1190/1.2748493. https://ui.adsabs.harvard.edu/abs/2007LeaEd..26..760H

  37. Lu et al. 2006 - 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: 2191–2195. doi:10.1190/1.2369970. https://doi.org/10.1190/1.2369970

  38. Hornby & Yu 2007 - Hornby, B.; Yu, J. (2007). "Interferometric imaging of a salt flank using walkaway VSP data". The Leading Edge. 26 (6): 760–763. Bibcode:2007LeaEd..26..760H. doi:10.1190/1.2748493. https://ui.adsabs.harvard.edu/abs/2007LeaEd..26..760H

  39. Lu et al. 2006 - 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: 2191–2195. doi:10.1190/1.2369970. https://doi.org/10.1190/1.2369970

  40. Schuster et al. 2004 - Schuster, G.; Yu, J.; Sheng, J.; Rickett, J. (2004). "Interferometric/daylight seismic imaging". Geophysical Journal International. 157 (2): 838–852. Bibcode:2004GeoJI.157..838S. doi:10.1111/j.1365-246x.2004.02251.x. https://doi.org/10.1111%2Fj.1365-246x.2004.02251.x

  41. Curtis et al. 2006 - Curtis, A.; Gerstoft, P.; Sato, H.; Snieder, R.; Wapenaar, K. (2006). "Seismic interferometry turning noise into signal". The Leading Edge. 25 (9): 1082–1092. Bibcode:2006LeaEd..25.1082C. doi:10.1190/1.2349814. http://resolver.tudelft.nl/uuid:bc099609-47bd-4ad8-91ad-2a4a9a949c0f

  42. Curtis et al. 2006 - Curtis, A.; Gerstoft, P.; Sato, H.; Snieder, R.; Wapenaar, K. (2006). "Seismic interferometry turning noise into signal". The Leading Edge. 25 (9): 1082–1092. Bibcode:2006LeaEd..25.1082C. doi:10.1190/1.2349814. http://resolver.tudelft.nl/uuid:bc099609-47bd-4ad8-91ad-2a4a9a949c0f

  43. Schuster et al. 2004 - Schuster, G.; Yu, J.; Sheng, J.; Rickett, J. (2004). "Interferometric/daylight seismic imaging". Geophysical Journal International. 157 (2): 838–852. Bibcode:2004GeoJI.157..838S. doi:10.1111/j.1365-246x.2004.02251.x. https://doi.org/10.1111%2Fj.1365-246x.2004.02251.x

  44. Curtis et al. 2006 - Curtis, A.; Gerstoft, P.; Sato, H.; Snieder, R.; Wapenaar, K. (2006). "Seismic interferometry turning noise into signal". The Leading Edge. 25 (9): 1082–1092. Bibcode:2006LeaEd..25.1082C. doi:10.1190/1.2349814. http://resolver.tudelft.nl/uuid:bc099609-47bd-4ad8-91ad-2a4a9a949c0f