On seismic interferometry, the generalized optical theorem, and the scattering matrix of a point scatterer
TL;DR: In this article, a far-field approximation of the Green's function representation for seismic interferometry is analyzed, and the generalized optical theorem is derived from the nonlinear scattering matrix of a point scatterer.
Abstract: 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.
Citations
More filters
TL;DR: In this paper, the authors present a new MATLAB code based on the Foldy method that models theoretically exact, direct, and scattered parts of a wavefield, while all multiple scattering interactions are taken into account,unlikefinite-difference or finite-element methods.
Abstract: Improvements in industrial seismic, seismological, acoustic, or interferometric theory and applications often result in quite subtle changes in sound quality, seismic images, or information which are nevertheless crucial for improved interpretation or experience. When evaluating new theories and algorithms using synthetic data, an important aspect of related researchisthereforethat numerical errors due towavefield modeling are reduced to a minimum. We present a new MATLAB code based on the Foldy method that models theoretically exact, direct, and scattered parts of a wavefield. Its main advantage lies in the fact that while all multiple scattering interactions aretakenintoaccount,unlikefinite-difference or finite-element methods, numerical dispersion errors are avoided. The method is therefore ideal for testing new theory in industrial seismics, seismology, acoustics, and in wavefield interferometry in particular because the latter is particularly sensitive to the dynamics of scattering interactions. We present the theory behind the Foldy acoustic modeling method and provide examples of its implementation. We also benchmark the code against a good finite-difference code. Because our Foldy code was written and optimized to test new theory in seismic interferometry, examples of its application to seismicinterferometryarealsopresented, showing its validity and importance when exact modeling results are needed.
19 citations
Cites background from "On seismic interferometry, the gene..."
...Each of these arrivals is the result of crosscorrelation of the scattered monopole wavefields to A and B, and its amplitude is proportional to the energy scattered by the diffractor that produced it (Snieder et al., 2008; Wapenaar et al., 2010b)....
[...]
TL;DR: It is shown that the reciprocity theorem for the scattered field is responsible for the cancellation of non-physical (spurious) arrivals in seismic interferometry, and as such provides the mathematical description of such arrivals.
Abstract: By analyzing correlation-type reciprocity theorems for wavefields in perturbed media, it is shown that the correlation-type reciprocity theorem for the scattered field is the progenitor of the generalized optical theorem. This reciprocity theorem, in contrast to the generalized optical theorem, allows for inhomogeneous background properties and does not make use of a far-field condition. This theorem specializes to the generalized optical theorem when considering a finite-size scatterer embedded in a homogeneous background medium and when utilizing the far-field condition. Moreover, it is shown that the reciprocity theorem for the scattered field is responsible for the cancellation of non-physical (spurious) arrivals in seismic interferometry, and as such provides the mathematical description of such arrivals. Even though here only acoustic waves are treated, the presented treatment is not limited to such wavefields and can be generalized to general wavefields. Therefore, this work provides the framework for deriving equivalents of the generalized optical theorem for general wavefields.
15 citations
TL;DR: Two-sided non-linear imaging provides the best elastic puremode (PP and SS) and converted-mode (PS) images, having higher resolution and more uniform illumination than those obtained from both one-sided linear imaging and from other intermediate steps of imaging as discussed by the authors.
Abstract: SUMMARY Source–receiver interferometric imaging can be used to synthesize a subsurface acoustic or elastic image, consisting of a zero-time, zero-offset response (or Green’s function) between a colocated pseudo-source and pseudo-receiver placed at each point in the subsurface image. However, if the imaging process does not properly account for multiple reflections, and enclosing boundaries of sources and receivers are not available, the image shows artefacts, poorly illuminated areas and distorted image amplitudes. Here we demonstrate with numerical examples that two-sided non-linear imaging provides the best elastic pure-mode (PP and SS) and converted-mode (PS) images, having higher resolution and more uniform illumination than those obtained from both one-sided linear imaging and from other intermediate steps of imaging (e.g. non-linear one-sided, linear two-sided). We also propose practical approaches to construct the additional fields required by two-sided non-linear imaging without the need for a detailed velocity model and receivers (and/or sources) in the subsurface. Moreover, when conversions are used for imaging, ‘true-amplitude’ images (here true-amplitude means properlyretrievingamplitudesthatrepresentthezero-time,zero-offsetelasticresponse)should theoretically vanish because neither P-to-S or S-to-P conversions arise at zero-time and zerooffset. Applying a correction procedure that accounts for the polarity reversal in PS (or SP) single-shot images helps with their structural interpretation but results in an unphysical estimate of the subsurface response and uninterpretable amplitudes. This suggests that there are advantages in exploiting pure-mode SS reflections/transmissions, in addition to converted waves only, because they require no polarity correction and the resulting image contains meaningful amplitudes that are proportional to the local shear-wave properties of the medium.
15 citations
TL;DR: In this paper, the authors show how to fingerprint individual diffractors inside an acoustic medium using interrogative wave energy from arrays of sources and receivers, and decompose complex multiply diffracted wavefields into constituent, single-diffraction interactions.
Abstract: SUMMARY We show how to ‘fingerprint’ individual diffractors inside an acoustic medium using interrogative wave energy from arrays of sources and receivers. For any recorded multiply diffracted wave observed between any source and any receiver, the set of such fingerprints is sufficient information to identify all diffractors involved in the corresponding diffraction path, and the sequential order in which diffractors are encountered. The method herein thus decomposes complex, multiply diffracted wavefields into constituent, single-diffraction interactions.
13 citations
Cites background from "On seismic interferometry, the gene..."
...Diffractors are generally characterized by a variety of properties, for example, spatial position and scattering amplitudes at U niversity of E http://gji.oxfordjournals.org/ D ow nloaded from or cross-sections (Groeneboom & Snieder 1995; Wapenaar et al. 2010)....
[...]
DOI•
07 Jun 2016
TL;DR: In this paper, the authors present a matrix-vector representation of the wave equation and derive new, alternative power-flux normalized eigenvector sets that prove to be numerically more stable and accurate.
Abstract: In this thesis, I study coupled poroelastic waves and electromagnetic fields in layered media. The focus is two-fold: 1. Increase the theoretical and physical understanding of the seismo-electromagnetic phenomenon by analytically-based numerical modeling. 2. Investigate the potential of seismo-electromagnetic interferometry. After presenting the governing equations that form the basis of the theoretical framework, I capture this system into a matrix-vector representation of the wave equation. I first use literature eigenvector sets, which I normalize with respect to power-flux. I then derive new, alternative power-flux normalized eigenvector sets that I prove to be numerically more stable and accurate. The eigenvector sets form the basis of the analytically-based numerical modeling code `ESSEMOD' that I developed to model seismo-electromagnetic wave/field propagation/diffusion in layered-Earth media. The alternative eigenvector set models scenarios with no seismo-electromagnetic coupling correctly, where the literature eigenvector sets fail. In addition, the alternative set properly deals with scenarios where both small amplitude signals and large amplitude signals occur in the record, whereas the literature eigenvector sets result in noise levels masking the small events. The same holds for scenarios with a small seismo-electromagnetic coupling coefficient. I design an effective global reflection scheme that properly describes the primary and multiple reflections in the models. I implement the correct boundary conditions to account for scenarios with a free-surface, and also for scenarios containing fluid/porous medium/fluid transitions. To transform all the seismo-electromagnetic source-receiver combinations in a numerically effective way back from the horizontal wavenumber-frequency domain to the space-frequency domain, I derive and implement explicit Fourier-Bessel transformations. I then validate the developed modeling code in numerous ways. First of all, I compare the results of seismo-EM layer-code modeling in a homogeneous medium with explicit homogeneous space Green's function expressions. This comparison provides a clear validation that the layer-code models the dynamic responses in homogeneous scenarios correctly. Next, I check numerical consistency by carrying out reciprocity checks. I study homogeneous space models, models containing a free-surface and models with interfaces. As a next step, I validate the modeling results of seismo-EM layer-code modeling for typical seismo-electromagnetic laboratory configurations, i.e. models containing fluid/porous medium/fluid transitions. I first compare the purely electromagnetic part of the seismo-EM layer-code with an independently developed purely electromagnetic layered-Earth code. The results match perfectly in both phase and amplitude for full transmission and pure reflection experiments, as well as for a combination of both. I then carry out a seismo-electromagnetic reciprocity test for a fluid halfspace overlying a porous medium halfspace, proving that the coupled poroelastic and electromagnetic fields are modeled consistently and yield the expected results. As a final validation step, I compare ESSEMOD with an independently developed seismo-electromagnetic layered-Earth modeling code. The results display an almost perfect match in both phase and relative amplitudes, and a constant amplitude correction factor of 4 needs to be applied to let the absolute amplitudes match. I then carry out a small feasibility test to study the potential of the seismo-electromagnetic effect for exploration purposes. I investigate different source-receiver combinations for the same model, and focus on the signal strength recorded at different distances from the target depth level. I conclude that for the source-receiver combinations studied, the electric field due to a volume injection monopole source, as well as the magnetic field due to a seismic bulk force source, yield the strongest converted signals. The receiver-distance from the target of interest plays an important role in the signal measurability. The closer the receivers to the target, the higher the signal strengths. However, when the receivers are located too close to the target, the coseismic reflected fields can mask the interface response fields that we are mainly interested in. Next, I study if nature itself can help us to overcome the very low signal-to-noise ratio of seismo-electromagnetic converted fields, by investigating the effects of thin-bed geological structures on the seismo-electromagnetic signal. To investigate the effects of bed-thinning on the seismo-electromagnetic interference patterns, I numerically simulate seismo-electromagnetic wave propagation through horizontally layered media with different amounts and thicknesses of thin-beds. I demonstrate seismo-electromagnetic sensitivity to changes in medium parameters on a spatial scale much smaller than the seismic resolution. By simulating moving oil/water contacts during production, where the oil layer is gradually being thinned, seismo-electromagnetic signals are proven very sensitive to oil/water contacts. I now explore the application of interferometric techniques to the seismo-electromagnetic system, which might eventually lead to an improved signal-to-noise ratio of the weak converted fields. I derive the theory for interferometric retrieval of 2D SH-TE seismo-electromagnetic Green's functions. Using both a circular source configuration and a line source configuration, I show that it is possible to correctly retrieve the dynamic seismo-electromagnetic 2D SH-TE response in a homogeneous medium, using seismic boundary sources only. Using seismo-EM layer-code data, I then show that it is also possible to correctly retrieve the direct shear wave-related causal coseismic field in a homogeneous medium, in both phase and amplitude. To obtain a perfect match in absolute amplitudes, I apply a single linear scaling factor. I finally carry out interferometric experiments in a model containing a single interface at 800 m depth, proving that it is possible to correctly retrieve all 2D SH-TE causal seismic-related direct and reflected coseismic fields, as well as interface response fields, by cross-correlation interferometry, using seismic boundary sources only. These results are promising for the application of 3D seismo-electromagnetic interferometry using seismo-EM layer-code modeling, and later on, in the field. Next, I present an alternative way to effectively decompose fields into their up- and downgoing components and different field types, using recordings at multiple depth levels. I present the theory of this MDL decomposition scheme, followed by successful decomposition of synthetic elastodynamic data sets. I additionally study the implications of laterally-varying media on the horizontal wavenumber-frequency domain MDL decomposition scheme. I demonstrate successful decomposition, using an acoustic approximation and applying a combined multi-component / MDL decomposition approach, of a field data set recorded in Annerveen, in the North of the Netherlands. I address how to effectively use the MDL decomposition scheme in a unified fashion, applied to all wave phenomena including seismo-electromagnetic phenomena. I then make a step towards seismo-electromagnetic inversion, presenting an effective way to carry out a seismo-electromagnetic sensitivity analysis using resolution functions. I start by explaining the theory of resolution functions using a seismo-electromagnetic example. I define the seismo-electromagnetic resolution function for inversion for a bulk density perturbation. I demonstrate the effectiveness of this method by first carrying out a purely electromagnetic sensitivity analysis for a point perturbation in conductivity, located in an isotropic homogeneous half-space. These results are compared with literature results based on analytical homogeneous space Green's function expressions. The result using the seismo-EM layer-code is nearly identical to the literature result. The position of the scatterer is correctly resolved. At the end of this section, I present the results of the fully-coupled seismo-electromagnetic senstivity analysis for a bulk density contrast for a specific source-receiver combination, using single-frequency multi-component line data. I show that the coupled seismo-electromagnetic system is sensitive to a perturbation in bulk density and that the position of the perturbation can be correctly recovered. I finalize this thesis by discussing potential seismo-electromagnetic applications, as well as by providing a brief outlook for future research.
12 citations
Cites background from "On seismic interferometry, the gene..."
...inverse problem, the contrast cannot be too large with respect to the background model, since otherwise the higher order scattering terms cannot be neglected (Born & Wolf, 1980; Habashy et al., 1993; Wapenaar et al., 2010)....
[...]
References
More filters
Book•
[...]
01 Jan 1959
TL;DR: In this paper, the authors discuss various topics about optics, such as geometrical theories, image forming instruments, and optics of metals and crystals, including interference, interferometers, and diffraction.
Abstract: The book is comprised of 15 chapters that discuss various topics about optics, such as geometrical theories, image forming instruments, and optics of metals and crystals. The text covers the elements of the theories of interference, interferometers, and diffraction. The book tackles several behaviors of light, including its diffraction when exposed to ultrasonic waves.
19,815 citations
[...]
01 Oct 1999
TL;DR: In this article, the authors discuss various topics about optics, such as geometrical theories, image forming instruments, and optics of metals and crystals, including interference, interferometers, and diffraction.
Abstract: The book is comprised of 15 chapters that discuss various topics about optics, such as geometrical theories, image forming instruments, and optics of metals and crystals. The text covers the elements of the theories of interference, interferometers, and diffraction. The book tackles several behaviors of light, including its diffraction when exposed to ultrasonic waves.
19,503 citations
9,185 citations
TL;DR: In this article, it was shown 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.
Abstract: 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.
700 citations