"uuid","repository link","title","author","contributor","publication year","abstract","subject topic","language","publication type","publisher","isbn","issn","patent","patent status","bibliographic note","access restriction","embargo date","faculty","department","research group","programme","project","coordinates"
"uuid:c4772e2c-caae-4123-b4af-c462f23f3489","http://resolver.tudelft.nl/uuid:c4772e2c-caae-4123-b4af-c462f23f3489","Seismic interferometry using multidimensional deconvolution and crosscorrelation for crosswell seismic reflection data without borehole sources","Minato, S.; Matsuoka, T.; Tsuji, T.; Draganov, D.S.; Hunziker, J.W.; Wapenaar, C.P.A.","","2011","Crosswell reflection method is a high-resolution seismic imaging method that uses recordings between boreholes. The need for downhole sources is a restrictive factor in its application, for example, to time-lapse surveys. An alternative is to use surface sources in combination with seismic interferometry. Seismic interferometry (SI) could retrieve the reflection response at one of the boreholes as if from a source inside the other borehole. We investigate the applicability of SI for the retrieval of the reflection response between two boreholes using numerically modeled field data. We compare two SI approaches — crosscorrelation (CC) and multidimensional deconvolution (MDD). SI by MDD is less sensitive to underillumination from the source distribution, but requires inversion of the recordings at one of the receiver arrays from all the available sources. We find that the inversion problem is ill-posed, and propose to stabilize it using singular-value decomposition. The results show that the reflections from deep boundaries are retrieved very well using both the CC and MDD methods. Furthermore, the MDD results exhibit more realistic amplitudes than those from the CC method for downgoing reflections from shallow boundaries. We find that the results retrieved from the application of both methods to field data agree well with crosswell seismic-reflection data using borehole sources and with the logged P-wave velocity.","geophysical techniques; interferometry; seismic waves; seismology","en","journal article","Society of Exploration Geophysicists","","","","","","","","Civil Engineering and Geosciences","Geoscience & Engineering","","","",""
"uuid:bffc466e-3073-47c0-a35d-32c6366ae2bf","http://resolver.tudelft.nl/uuid:bffc466e-3073-47c0-a35d-32c6366ae2bf","High-resolution lithospheric imaging with seismic interferometry","Ruigrok, E.N.; Campman, X.; Draganov, D.S.; Wapenaar, K.","","2010","In recent years, there has been an increase in the deployment of relatively dense arrays of seismic stations. The availability of spatially densely sampled global and regional seismic data has stimulated the adoption of industry-style imaging algorithms applied to converted- and scattered-wave energy from distant earthquakes, leading to relatively high-resolution images of the lower crust and upper mantle.We use seismic interferometry to extract reflection responses from the coda of transmitted energy from distant earthquakes. In theory, higher resolution images can be obtained when migrating reflections obtained with seismic interferometry rather than with conversions, traditionally used in lithospheric imaging methods. Moreover, reflection data allow the straightforward application of algorithms previously developed in exploration seismology. In particular, the availability of reflection data allows us to extract from it a velocity model using standard multichannel data-processing methods. However, the success of our approach relies mainly on a favourable distribution of earthquakes. In this paper, we investigate how the quality of the reflection response obtained with interferometry is influenced by the distribution of earthquakes and the complexity of the transmitted wavefields. Our analysis shows that a reasonable reflection response could be extracted if (1) the array is approximately aligned with an active zone of earthquakes, (2) different phase responses are used to gather adequate angular illumination of the array and (3) the illumination directions are properly accounted for during processing. We illustrate our analysis using a synthetic data set with similar illumination and source-side reverberation characteristics as field data recorded during the 2000–2001 Laramie broad-band experiment. Finally, we apply our method to the Laramie data, retrieving reflection data. We extract a 2-D velocity model from the reflections and use this model to migrate the data. On the final reflectivity image, we observe a discontinuity in the reflections. We interpret this discontinuity as the Cheyenne Belt, a suture zone between Archean and Proterozoic terranes.","seismology; interferometry; body waves; crustal structure","en","journal article","Wiley-Blackwell","","","","","","","","Civil Engineering and Geosciences","Geotechnology","","","",""
"uuid:a7510463-5446-4a49-a8fc-81f44db1d984","http://resolver.tudelft.nl/uuid:a7510463-5446-4a49-a8fc-81f44db1d984","Tutorial on seismic interferometry: Part 1 — Basic principles and applications","Wapenaar, C.P.A.; Draganov, D.S.; Snieder, R.; Campman, X.; Verdel, A.","","2010","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.","geophysical techniques; Green's function methods; interferometry; seismic waves; seismology","en","journal article","Society of Exploration Geophysicists","","","","","","","","Civil Engineering and Geosciences","Geoscience & Engineering","","","",""
"uuid:07504f32-d9fb-46b9-8095-dcfa5b3e817b","http://resolver.tudelft.nl/uuid:07504f32-d9fb-46b9-8095-dcfa5b3e817b","Tutorial on seismic interferometry: Part 2 — Underlying theory and new advances","Wapenaar, C.P.A.; Slob, E.C.; Snieder, R.; Curtis, A.","","2010","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.","deconvolution; geophysical techniques; Green's function methods; interferometry; seismic waves; seismology","en","journal article","Society of Exploration Geophysicists","","","","","","","","Civil Engineering and Geosciences","Geoscience & Engineering","","","",""
"uuid:9bf7fc0a-55ba-4cd0-91dd-b51b2064b638","http://resolver.tudelft.nl/uuid:9bf7fc0a-55ba-4cd0-91dd-b51b2064b638","On seismic interferometry, the generalized optical theorem, and the scattering matrix of a point scatterer","Wapenaar, C.P.A.; Slob, E.C.; Snieder, R.","","2010","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.","geophysical techniques; Green's function methods; interferometry; seismic waves; seismology","en","journal article","Society of Exploration Geophysicists","","","","","","","","Civil Engineering and Geosciences","Geoscience & Engineering","","","",""
"uuid:bc099609-47bd-4ad8-91ad-2a4a9a949c0f","http://resolver.tudelft.nl/uuid:bc099609-47bd-4ad8-91ad-2a4a9a949c0f","Seismic interferometry-turning noise into signal","Curtis, A.; Gerstoft, P.; Sato, H.; Snieder, R.; Wapenaar, C.P.A.","","2006","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.","geophysical techniques; seismology; structural engineering; earthquakes; interferometry","en","journal article","Society of Exploration Geophysicists","","","","","","","","Civil Engineering and Geosciences","","","","",""
"uuid:fc9a5a03-cbfa-40ca-8f4d-9652ecd325f5","http://resolver.tudelft.nl/uuid:fc9a5a03-cbfa-40ca-8f4d-9652ecd325f5","Seismic interferometry: Reconstructing the earth's reflection response","Draganov, D.S.; Wapenaar, C.P.A.; Thorbecke, J.W.","","2006","In 1968, Jon Claerbout showed that the reflection response of a 1D acoustic medium can be reconstructed by autocorrelating the transmission response. Since then, several authors have derived relationships for reconstructing Green's functions at the surface, using crosscorrelations of (noise) recordings that were taken at the surface and that derived from subsurface sources.For acoustic media, we review relations between the reflection response and the transmission response in 3D inhomogeneous lossless media. These relations are derived from a one-way wavefield reciprocity theorem. We use modeling results to show how to reconstruct the reflection response in the presence of transient subsurface sources with distinct excitation times, as well as in the presence of simultaneously acting noise sources in the subsurface. We show that the quality of reconstructed reflections depends on the distribution of the subsurface sources. For a situation with enough subsurface sources — that is, for a distribution that illuminates the subsurface area of interest from nearly alldirections — the reconstructed reflection responses and the migrated depth image exhibit all the reflection events and the subsurface structures of interest, respectively. With only a few subsurface sources, that is, with insufficient illumination, the reconstructed reflection responses are noisy and can even become kinematically incorrect. At the same time, however, the depth image, which was obtained from their migration, still shows clearly all the illuminated subsurface structures at their correct positions.For the elastic case, we review a relationship between the reflection Green's functions and the transmission Green's functions derived from a two-way wavefield reciprocity theorem. Using modeling examples, we show how to reconstruct the different components of the particle velocity observed at the surface and resulting from a surface traction source. This reconstruciton is achieved using crosscorrelations of particle velocity components measured at the surface and resulting from separate P- and S-wave sources in the subsurface.","seismology; interferometry; seismic waves; Green's function methods","en","journal article","Society of Exploration Geophysicists","","","","","","","","Civil Engineering and Geosciences","Geotechnology","","","",""
"uuid:68e13eb2-52e7-499c-8a04-ddee6fa0d6dd","http://resolver.tudelft.nl/uuid:68e13eb2-52e7-499c-8a04-ddee6fa0d6dd","Green's function representations for seismic interferometry","Wapenaar, C.P.A.; Fokkema, J.T.","","2006","The term seismic interferometry refers to the principle of generating new seismic responses by crosscorrelating seismic observations at different receiver locations. The first version of this principle was derived by Claerbout (1968), who showed that the reflection response of a horizontally layered medium can be synthesized from the autocorrelation of its transmission response. For an arbitrary 3D inhomogeneous lossless medium it follows from Rayleigh's reciprocity theorem and the principle of time-reversal invariance that the acoustic Green's function between any two points in the medium can be represented by an integral of crosscorrelations of wavefield observations at those two points. The integral is along sources on an arbitrarily shaped surface enclosing these points. No assumptions are made with respect to the diffusivity of the wavefield. The Rayleigh-Betti reciprocity theorem leads to a similar representation of the elastodynamic Green's function. When a part of the enclosing surface is the earth's free surface, the integral needs only to be evaluated over the remaining part of the closed surface. In practice, not all sources are equally important: The main contributions to the reconstructed Green's function come from sources at stationary points. When the sources emit transient signals, a shaping filter can be applied to correct for the differences in source wavelets. When the sources are uncorrelated noise sources, the representation simplifies to a direct crosscorrelation of wavefield observations at two points, similar as in methods that retrieve Green's functions from diffuse wavefields in disordered media or in finite media with an irregular bounding surface.","seismology; interferometry; seismic waves; Green's function methods","en","journal article","Society of Exploration Geophysicists","","","","","","","","Civil Engineering and Geosciences","Geotechnology","","","",""
"uuid:c0512fe5-692e-47e3-98a2-a24cf29c0d09","http://resolver.tudelft.nl/uuid:c0512fe5-692e-47e3-98a2-a24cf29c0d09","Spurious multiples in seismic interferometry of primaries","Snieder, R.; Wapenaar, C.P.A.; Larner, K.","","2006","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; seismic waves; seismology","en","journal article","Society of Exploration Geophysicists","","","","","","","","Civil Engineering and Geosciences","Geotechnology","","","",""