; Oh+'0HP
$TU Delft Repository search results0TU Delft Repository search results (max. 1000)TU Delft LibraryTU Delft Library@[^\@[^\՜.+,0HPX`hp
x
WorksheetFeuilles de calcul
B=%r8X"1Calibri1Calibri1Calibri1
Calibri 83ffff̙̙3f3fff3f3f33333f33333.JTU Delft Repositoryg fuuidrepository linktitleauthorcontributorpublication yearabstract
subject topiclanguagepublication type publisherisbnissnpatent
patent statusbibliographic noteaccess restrictionembargo datefaculty
departmentresearch group programmeprojectcoordinates)uuid:b959853fcc3443e7add28e2c68620e16Dhttp://resolver.tudelft.nl/uuid:b959853fcc3443e7add28e2c68620e16(Adaptive tracking in seismic explorationGeerlings, A.C.Berkhout, A.J. (promotor)1seismology; heuristic search; pattern recognitionendoctoral thesisApplied Sciences)uuid:09df29e294874fe687e0c2a9470933d7Dhttp://resolver.tudelft.nl/uuid:09df29e294874fe687e0c2a9470933d7;Surfacerelated multiple elimination, an inversion approachVerschuur, D.J.
seismology)uuid:1b7ad6beb31e4b0888cd201049bbf8feDhttp://resolver.tudelft.nl/uuid:1b7ad6beb31e4b0888cd201049bbf8feLinear avo inversion by prestack depth migration: Imaging angle dependent reflectivity as a tool for lithostratigraphic inversionDe Bruin, C.G.M.seismology; lithology)uuid:10ed7f4d81574cdba97621dc65175ee5Dhttp://resolver.tudelft.nl/uuid:10ed7f4d81574cdba97621dc65175ee55Controlled Illumination in Prestack Seismic MigrationRietveld, W.E.A.4seismology; wave field synthesis; prestack migration)uuid:15232dc10d984ba2bdaef59f520134b4Dhttp://resolver.tudelft.nl/uuid:15232dc10d984ba2bdaef59f520134b4*Crosswell Seismic Tomography and MigrationGiling, E.J.M.(seismology; tomography; wave propagation)uuid:c3fd77aa1a0f4a7192d286a722ed1366Dhttp://resolver.tudelft.nl/uuid:c3fd77aa1a0f4a7192d286a722ed1366fRetrieving the Green s function in an open system by cross correlation: A comparison of approaches (L)*Wapenaar, C.P.A.; Fokkema, J.; Snieder, R.&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 timereversal 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 timereversal approach can be obtained as an approximation of the result of the reciprocity approach.Green's function methods; acoustic wave propagation; acoustic wave scattering; vibrations; structural acoustics; acoustic signal processing; seismologyjournal articleAcoustical Society of America!Civil Engineering and Geosciences
Geotechnology)uuid:2933280d74b44bbbb46d77b5b26a59aaDhttp://resolver.tudelft.nl/uuid:2933280d74b44bbbb46d77b5b26a59aa@A new elastic model for ground coupling of geophones with spikes=Drijkoningen, G.G.; Rademakers, F.; Slob, E.C.; Fokkema, J.T.Ground coupling are terms that describe the transfer from seismic ground motion to the motion of a geophone. In previous models, ground coupling was mainly considered as a disk lying on top of a halfspace, not considering the fact that in current practice geophones are spiked and are buried for optimal response. In this paper we introduce a new model that captures the spike added to the geophone and models the effect of geophone burial. The geophone is modeled as a rigid, movable cylinder embedded in a halfspace near or at the surface. The coupling problem is then tackled by a scattering approach using the elastic form of reciprocity; we consider the vertical component only. The main feature in the coupling function is a resonance whose location and shape depend on the different parameters of the geophone and the soil. In accordance with previous models, adding mass reduces the frequency of resonance. However, we show that pure mass loading assumption is too restrictive for standard geophones. Our new model shows that increasing the spike radius and length decreases the frequency of resonance and the resonance is more peaked. Furthermore, burying the geophone decreases the frequency of resonance, but when one takes into account that the soil at depth is more compact, then the behavior i< s as observed in practice namely, an increase in frequency of resonance. As for the properties of the soil, the shearwave velocity has the largest effect; when increased, it shifts the frequency of resonance to the highfrequency end as desired.8seismic waves; soil; seismometers; seismology; resonance$Society of Exploration Geophysicists)uuid:eefac5acf44046808978a26848f2230fDhttp://resolver.tudelft.nl/uuid:eefac5acf44046808978a26848f2230fYDiscrimination between phase and amplitude attributes in timelapse seismic streamer dataSpetzler, J.; Kvam, O.Timelapse seismic experiments aim to obtain information about productionrelated effects in hydrocarbon reservoirs to increase the recovery percentage. However, nonrepeatability problems such as acquisition differences, overburden effects, and noise are often significantly stronger than the imprint of production changes in timelapse seismic data sets. Consequently, it is very difficult to appraise the changes in petrophysical reservoir parameters over time. We introduce a 4D monitoring approach based on the spectral ratio method. This method produces two timelapse attributes: the relative change in reflection coefficient and the traveltime shift at reflecting interfaces. These attributes can be used for appraising productionrelated changes in the subsurface. The approach corrects for timeinvariant nonrepeatability effects in the overburden and sourcereceiver coupling problems in timelapse surveys. The validity of the method is limited to structurally simple overburden and reservoirs with weak lateral variations. First, we validate the methodology using a synthetic timelapse seismic experiment. Next, we apply the method to a real timelapse data set from the Troll West gas province in the North Sea. In the real example, we could not detect movement in the fluid contact of 5 15 m. The expected change in amplitude is less than 10%, which is probably below the background noise level for this data set.rocks; seismology)uuid:85b144578770463aaf0cff2243ac1b01Dhttp://resolver.tudelft.nl/uuid:85b144578770463aaf0cff2243ac1b01,Seismic processing in the inverse data spaceBerkhout, A.J.
Until now, seismic processing has been carried out by applying inverse filters in the forward data space. Because the acquired data of a seismic survey is always discrete, seismic measurements in the forward data space can be arranged conveniently in a data matrix (P). Each column in the data matrix represents one shot record. If we represent seismic data in the temporal frequency domain, then each matrix element consists of a complexvalued number. Considering the dominant role of multiple scattering in seismic data, it is proposed to replace data matrix P by its inverse P 1 before starting seismic processing. Making use of the feedback model for seismic data, multiple scattered energy is mapped onto the zero time axis of the inverse data space. The practical consequence of this remarkable property may be significant: multiple elimination in the inverse data space simplifies to removing data at zero time only. Moving to the inverse data space may cause a fundamental change in the way we preprocess and image seismic data.eseismology; inverse problems; geophysical techniques; geophysical signal processing; matrix inversion)uuid:c0512fe5692e47e398a2a24cf29c0d09Dhttp://resolver.tudelft.nl/uuid:c0512fe5692e47e398a2a24cf29c0d099Spurious multiples in seismic interferometry of primaries)Snieder, R.; Wapenaar, C.P.A.; Larner, K.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 c< orrectly 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)uuid:fc9a5a03cbfa40ca8f4d9652ecd325f5Dhttp://resolver.tudelft.nl/uuid:fc9a5a03cbfa40ca8f4d9652ecd325f5FSeismic interferometry: Reconstructing the earth's reflection response1Draganov, D.S.; Wapenaar, C.P.A.; Thorbecke, J.W.PIn 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 oneway 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 twoway 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 Swave sources in the subsurface.Cseismology; interferometry; seismic waves; Green's function methods)uuid:68e13eb252e7499c8a04ddee6fa0d6ddDhttp://resolver.tudelft.nl/uuid:68e13eb252e7499c8a04ddee6fa0d6dd;Green's function representations for seismic interferometryWapenaar, C.P.A.; Fokkema, J.T.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 timereversal 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 RayleighBetti reciprocity theorem leads to a similar representation of the elastodynamic Green's function. Whe< n 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.)uuid:4e602a55c4b647c78d81bce52c48df6fDhttp://resolver.tudelft.nl/uuid:4e602a55c4b647c78d81bce52c48df6fImaging of multiple reflectionsBerkhout, A.J.; Verschuur, D.J.Current multipleremoval algorithms in seismic processing use either differential moveout or predictability. If the differential moveout between primaries and multiples is small, prediction is the only option available. In the last decade, multidimensional predictionerror filtering by weighted convolution, such as surfacerelated multiple elimination (SRME), have proved to be very successful in practice. So far, multiples have been considered as noise and have been discarded after the removal process. In this paper, we argue that multiple reflections contain a wealth of information that can be used in seismic processing to improve the resolution of reservoir images beyond current capability. In the near future, one may expect that the socalled weightedcrosscorrelation (WCC) concept may offer an attractive alternative in approaching the multiple problem. WCC creates an option to avoid the adaptive subtraction process as applied in predictionerror algorithms. Moreover, it allows the transformation of multiples into primaries. The latter means that seismic imaging with primaries and multiples (nonlinear process) can be implemented by a sequence of linear processes, including the transformation of multiples into primaries and the imaging of primaries.4seismology; imaging; seismic waves; Radon transforms)uuid:bc09960947bd4ad891ad2a4a9a949c0fDhttp://resolver.tudelft.nl/uuid:bc09960947bd4ad891ad2a4a9a949c0f0Seismic interferometryturning noise into signalACurtis, A.; Gerstoft, P.; Sato, H.; Snieder, R.; Wapenaar, C.P.A.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 complexlooking 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.Wgeophysical techniques; seismology; structural engineering; earthquakes; interferometry)uuid:a3b3a766c9214ac9826276ce9b68cf6fDhttp://resolver.tudelft.nl/uuid:a3b3a766c9214ac9826276ce9b68cf6fQFocal transformation, an imaging concept for signal restoration and noise removalInterpolation of data beyond aliasing limits and removal of noise that occurs within the seismic bandwidth are still important problems in seismic processing. The focal transform is introduced as a promising tool in data interpolation and noise removal, allowing the incorporation of macroinformation about the involved wavefields. From a physical point of view, the principal action of the forward focal ope< rator is removing the spatial phase of the signal content from the input data, and the inverse focal operator restores what the forward operator has removed. The strength of the method is that in the transformed domain, the focused signals at the focal area can be separated from the dispersed noise away from the focal area. Applications of particular interest in preprocessing are interpolation of missing offsets and reconstruction of signal beyond aliasing. The latter can be seen as the removal of aliasing noise.wgeophysical signal processing; signal reconstruction; signal restoration; imaging; seismology; interference suppression)uuid:97664d6d4fb1494aa99f16f75e5f9e9dDhttp://resolver.tudelft.nl/uuid:97664d6d4fb1494aa99f16f75e5f9e9dJFourier reconstruction of marinestreamer data in four spatial coordinatesZwartjes, P.M.; Gisolf, A.Many methods exist for interpolation of seismic data in one and two spatial dimensions, but few can interpolate properly in three or four spatial dimensions. Marine multistreamer data typically are sampled relatively well in the midpoint and absolute offset coordinates but not in the azimuth because the crossline shot coordinate is significantly under sampled. We approach the problem of interpolation of marinestreamer data in four spatial dimensions by splitting the problem into a 1D interpolation along the densely sampled streamers and a 3D Fourier reconstruction for the remaining spatial coordinates. In Fourier reconstruction, the Fourier coefficients that synthesize the nonuniformly sampled seismic data are estimated in a leastsquares inversion. The method is computationally efficient, requires no subsurface information, and can handle uniform grids with missing data as well as nonuniform grids or random sampling.The output grid of the 1D interpolation in the first step is arbitrary. When the output grid has uniform inline midpoints spacing, the 3D Fourier reconstruction in the second step is performed in the crossline midpoint, absolute offset, and azimuth coordinates. When the first step outputs to uniform absolute offset, the 3D Fourier reconstruction handles the crossline/inline midpoint and the azimuth coordinates. In both cases, the main innovation is the inclusion of the azimuthal coordinate in the Fourier reconstruction. The azimuth multiplicity must be increased for the method to be successful, which means that overlap shooting is required. We have tested the algorithm on synthetic streamer data for which the proposed method outperforms an approach where the azimuthal coordinate is ignored. Potential applications are interpolation of marine streamer data to decrease the crossline source sampling for the benefit of 3D multiple prediction and regularization to reduce samplingrelated differences in processing of timelapse data.seismology; Fourier transforms)uuid:6d29c07992eb43e0b1c478b461598ebfDhttp://resolver.tudelft.nl/uuid:6d29c07992eb43e0b1c478b461598ebfPlanewave depth migrationFStoffa, P.L.; Sen, M.K.; Seifoullaev, R.K.; Pestana, R.; Fokkema, J.T.]We present fast and efficient planewave migration methods for densely sampled seismic data in both the source and receiver domains. The methods are based on slant stacking over both shot and receiver positions (or offsets) for all the recorded data. If the dataacquisition geometry permits, both inline and crossline source and receiver positions can be incorporated into a multidimensional phasevelocity space, which is regular even for randomly positioned input data. By noting the maximum time dips present in the shot and receiver gathers and constantoffset sections, the number of plane waves required can be estimated, and this generally results in a reduction of the data volume used for migration. The required traveltime computations for depth imaging are independent for each particular planewave component. It thus can be used for either the source or the receiver plane waves during extrapolation in phase space, reducing considerably the computational burden. Since only vertical delay times are required, many travelti< me techniques can be employed, and the problems with multipathing and first arrivals are either reduced or eliminated. Further, the planewave integrals can be pruned to concentrate the image on selected targets. In this way, the computation time can be further reduced, and the technique lends itself naturally to a velocitymodeling scheme where, for example, horizontal and then steeply dipping events are gradually introduced into the velocity analysis. The migration method also lends itself to imaging in anisotropic media because phase space is the natural domain for such an analysis.seismic waves; seismology)uuid:78eea899c7f443839eb89eac4bf7ec0eDhttp://resolver.tudelft.nl/uuid:78eea899c7f443839eb89eac4bf7ec0edValidation of firstorder diffraction theory for the traveltimes and amplitudes of propagating waves9Jocker, J.; Spetzler, J.; Smeulders, D.M.J.; Trampert, J.Ultrasonic measurements of acoustic wavefields scattered by single spheres placed in a homogenous background medium (water) are presented. The dimensions of the spheres are comparable to the wavelength and the wavelength and represent both positive (rubber) and negative (teflon) velocity anomalies with respect to the background medium. The sensitivity of the recorded wavefield to scattering in terms of traveltime delay and amplitude variation is investigated. The results validate a linear (firstorder) diffraction theory for wavefields propagating in heterogeneous media with anomaly contrasts on the order of 15%. The diffraction theory is compared further with the exact results known from literature for scattering from an elastic sphere, formulated in terms of Legendre polynomials. To investigate the 2D case, the firstorder scattering theory is tested against 2D elastic finitedifference calculations. As the presented theory involves a volume integral, it is applicable to any geometric shape, and the scattering object does not need to be spherical or any other specific symmetrical shape. Furthermore, it can be implemented easily in seismic data inversion schemes, which is illustrated with examples from seismic crosswell tomography and a reflection experiment.Lseismic waves; seismology; acoustic waves; rubber; finite difference methods)uuid:4c09ff8864be40529645bd08a02b1424Dhttp://resolver.tudelft.nl/uuid:4c09ff8864be40529645bd08a02b14244Recursive prestack depth migration using CFP gathersThorbecke, J.; Berkhout, A.J.The commonfocuspoint technology (CFP) describes prestack migration by focusing in two steps: emission and detection. The output of the first focusing step represents a CFP gather. This gather defines a shot record that represents the subsurface response resulting from a focused source wavefield. We propose applying the recursive shotrecord, depthmigration algorithm to the CFP gathers of a seismic data volume and refer to this process as CFPgather migration. In the situation of complex geology and/or low signaltonoise ratio, CFPbased image gathers are easier to interpret for nonalignment than the conventional image gathers. This makes the CFPbased image gathers better suited for velocity analysis. This important property is illustrated by examples on the Marmousi model.seismology; seismic waves)uuid:389fad2146c34368bf75ea6f7b1f8586Dhttp://resolver.tudelft.nl/uuid:389fad2146c34368bf75ea6f7b1f8586The dynamics of statics<Turhan Taner, M.; Berkhout, A.J.; Treitel, S.; Kelamis, P.G.The statics problem, whether short wavelength, long wavelength, residual, or trim, has always been one of the more timeconsuming and problematic steps in seismic data processing. We routinely struggle with issues such as poor signaltonoise (S/N) ratio, cycle skipping, truncated refractors, wavelets with ambiguous first arrival times, etc. Elevation variations create their own problems and impact the choice of datum floating, phantom or recourse to a zerovelocity layer. Even if we can overcome some of these problems, we still have a catch 22 situation in which accurate velocity estimation requires good statics, while good statics estima< tion requires accurate velocities. To characterize these ambiguities, we have come up the oxymoron timevarying statics. Ygeophysical techniques; seismology; seismic waves; geophysical signal processing; statics)uuid:d318d89877994c6ab3393fde71306c9eDhttp://resolver.tudelft.nl/uuid:d318d89877994c6ab3393fde71306c9ezApplication of a linear finitefrequency theory to timelapse crosswell tomography in ultrasonic and numerical experiments)Spetzler, J.; Sijacic, D.; Wolf, K.H.A.A.MTimelapse seismic monitoring is the geophysical discipline whereby multiple data sets recorded at the same location but at different times are used to locate and quantify temporal changes in the elastic parameters of the subsurface. We validate a timelapse monitoring method by crosswell tomography using two types of wavefieldmodeling experiments: (1) a 3D real ultrasonic waveform experiment and (2) 2D synthetic finitedifference wavefield simulations. For both wavefield experiments, a timelapse structure simulating a fluid sweep in a reservoir layer is applied. The timelapse tomographic monitoring approach is based on the standard ray theory and a finitefrequency wavefield theory, where the latter takes into account the finitefrequency properties of recorded wavefields. The inverted timelapse models compiled with either the ray theory or the finitefrequency wavefield theory locate and correctly quantify the flooding zone in the simulated fluid sweep model. Both wavefield theories provide an adequate result because the flooding zone is comparable in size to the Fresnel volume.dfinite difference methods; geophysical techniques; hydrocarbon reservoirs; seismic waves; seismology)uuid:e599c89e0186472394b5364e85ddb0afDhttp://resolver.tudelft.nl/uuid:e599c89e0186472394b5364e85ddb0afTOn the relation between seismic interferometry and the migration resolution function!Thorbecke, J.W.; Wapenaar, C.P.A.Seismic interferometry refers to the process of retrieving new seismic responses by crosscorrelating seismic observations at different receiver locations. Seismic migration is the process of forming an image of the subsurface by wavefield extrapolation. Comparing the expressions for backward propagation known from migration literature with the Green's function representations for seismic interferometry reveals that these seemingly distinct concepts are mathematically equivalent. The frequencydomain representation for the resolution function of migration is identical to that for the Green's function retrieved by seismic interferometry (or its square, in the case of double focusing). In practice, they differ because the involved Green's functions in seismic interferometry are all defined in the actual medium, whereas in migration one of the Green's functions is defined in a background medium.<geophysical techniques; Green's function methods; seismology)uuid:b68278f691fb471b9d49decc2927f766Dhttp://resolver.tudelft.nl/uuid:b68278f691fb471b9d49decc2927f766KSimulating migrated and inverted seismic data by filtering a geologic modelXToxopeus, G.; Thorbecke, J.W.; Wapenaar, C.P.A.; Petersen, S.; Slob, E.C.; Fokkema, J.T.%The simulation of migrated and inverted data is hampered by the high computational cost of generating 3D synthetic data, followed by processes of migration and inversion. For example, simulating the migrated seismic signature of subtle stratigraphic traps demands the expensive exercise of 3D forward modeling, followed by 3D migration of the synthetic seismograms. This computational cost can be overcome using a strategy for simulating migrated and inverted data by filtering a geologic model with 3D spatialresolution and angle filters, respectively. A key property of the approach is this: The geologic model that describes a target zone is decoupled from the macrovelocity model used to compute the filters. The process enables a targetorientedapproach, by which a geologically detailed earth model describing a reservoir is adjusted without having to recalculate the filters. Because a spatialresolution filter combines the resul< ts of the modeling and migration operators, the simulated images can be compared directly to a real migration image. We decompose the spatialresolution filter into two parts and show that applying one of those parts produces output directly comparable to 1D inverted real data. Twodimensional synthetic examples that include seismic uncertainties demonstrate the usefulness of the approach. Results from a real data example show that horizontal smearing, which is not simulated by the 1D convolution model result, is essential to understand the seismic expression of the deformation related to sulfate dissolution and karst collapse.Bdeformation; geochemistry; seismic waves; seismology; stratigraphy)uuid:7d656a13c0f1425a8835616a1ccfdb0aDhttp://resolver.tudelft.nl/uuid:7d656a13c0f1425a8835616a1ccfdb0aBGlobalscale seismic interferometry: Theory and numerical examples+Ruigrok, E.N.; Draganov, D.S.; Wapenaar, K.Progress in the imaging of the mantle and core is partially limited by the sparse distribution of natural sources; the earthquake hypocenters are mainly along the active lithospheric plate boundaries. This problem can be approached with seismic interferometry. In recent years, there has been considerable progress in the development of seismic interferometric techniques. The term seismic interferometry refers to the principle of generating new seismic responses by crosscorrelating seismic observations at different receiver locations. The application of interferometric techniques on a global scale could create sources at locations where no earthquakes occur. In this way, yet unknown responses would become available for the application of traveltime tomography and surfacewave dispersion studies. The retrieval of a denseenough sampling of source gathers would largely benefit the application of reflection imaging. We derive new elastodynamic representation integrals for globalscale seismic interferometry. The relations are different from other seismic interferometry relations for transient sources, in the sense that they are suited for a rotating closed system like the Earth. We use a correlation of an observed response with a response to which freesurface multiple elimination has been applied to account for the closed system. Despite the fact that the rotation of the Earth breaks sourcereceiver reciprocity, the seismic interferometry relations are shown to be valid. The Coriolis force is included without the need to evaluate an extra term. We synthesize globalscale earthquake responses and use them to illustrate the acoustic versions of the new interferometric relations. When the sampling of real source locations is dense enough, then both the responses with and without freesurface multiples are retrieved. When we do not take into account the responses from the sources in the direct neighborhood of the seismic interferometryconstructed source location, the response with freesurface multiples can still be retrieved. Even when only responses from sources at a certain range of epicentral distances are available, some events in the Green s function between two receiver locations can still be retrieved. The retrieved responses are not perfect, but the artefacts can largely be ascribed to numerical errors. The reconstruction of internal events the response as if there was a source and a receiver on (major) contrasts within the model could possibly be of use for imaging. With modelling it is possible to discover in which region of the correlation panel stationary phases occur that contribute to the retrieval of events. This knowledge opens up a new way of filtering out undesired events and of discovering whether specific events could be retrieved with a given sourcereceiver configuration..seismology; body waves; seismic interferometryWileyBlackwell)uuid:cacc1e07acd743a68d39b16ad3a3abf0Dhttp://resolver.tudelft.nl/uuid:cacc1e07acd743a68d39b16ad3a3abf04Adaptive curveletdomain primarymultiple separation)Herrmann, F.J.; Wang, D.; Verschuur, D.J.In many exploration areas, successful separation of primaries and m< ultiples greatly determines the quality of seismic imaging. Despite major advances made by surfacerelated multiple elimination (SRME), amplitude errors in the predicted multiples remain a problem. When these errors vary for each type of multiple in different ways (as a function of offset, time, and dip), they pose a serious challenge for conventional leastsquares matching and for the recently introduced separation by curveletdomain thresholding. We propose a dataadaptive method that corrects amplitude errors, which vary smoothly as a function of location, scale (frequency band), and angle. With this method, the amplitudes can be corrected by an elementwise curveletdomain scaling of the predicted multiples. We show that this scaling leads to successful estimation of primaries, despite amplitude, sign, timing, and phase errors in the predicted multiples. Our results on synthetic and real data show distinct improvements over conventional leastsquares matching in terms of better suppression of multiple energy and highfrequency clutter and better recovery of estimated primaries.@geophysical techniques; least squares approximations; seismologyImaging Science and Technology)uuid:b2b336646da1418f9af45ab98fdd439dDhttp://resolver.tudelft.nl/uuid:b2b336646da1418f9af45ab98fdd439d1Acquisition geometry analysis in complex 3D media4Van Veldhuizen, E.J.; Blacquiere, G.; Berkhout, A.J.WIncreasingly, we must deal with complex subsurface structures in seismic exploration, often resulting in poor illumination and, therefore, poor image quality. Consequently, it is desirable to take into consideration the effects of wave propagation in the subsurface structure when designing an acquisition geometry. We developed a new, modelbased implementation of the previously introduced focalbeam analysis method. The method's objective is to provide quantitative insight into the combined influence of acquisition geometry, overburden structure, and migration operators on image resolution and angledependent amplitude accuracy. This is achieved by simulation of migrated gridpoint responses using focal beams. Note that the seismic response of any subsurface can be composed of a linear sum of gridpoint responses. The focal beams have been chosen because any migration process represents double focusing. In addition, the focal source beam and focal detector beam relate migration quality to illumination properties of the source geometry and sensing properties of the detector geometry, respectively. Waveequation modeling ensures that frequencydependent effects in the seismicfrequency range are incorporated. We tested our method by application to a 3D salt model in the Gulf of Mexico. Investigation of wellsampled, allazimuth, longoffset acquisition geometries revealed fundamental illumination and sensing limitations. Further results exposed the shortcomings of narrowazimuth data acquisition. The method also demonstrates how acquisitionrelated amplitude errors affect seismic inversion results.Tdata acquisition; geophysical prospecting; geophysical signal processing; seismology)uuid:4ce33f95c28047f2b443d39eb24f7eeaDhttp://resolver.tudelft.nl/uuid:4ce33f95c28047f2b443d39eb24f7eeakThe spatial dataadaptive minimumvariance distortionlessresponse beamformer on seismic singlesensor dataPanea, I.; Drijkoningen, G.G.Coherent noise generated by surface waves or ground roll within a heterogeneous near surface is a major problem in land seismic data. Array forming based on singlesensor recordings might reduce such noise more robustly than conventional hardwired arrays. We use the minimumvariance distortionlessresponse (MVDR) beamformer to remove (aliased) surfacewave energy from singlesensor data. This beamformer is data adaptive and robust when the presumed and actual desired signals are mismatched. We compute the intertrace covariance for the desired signal, and then for the total signal (desired signal+noise) to obtain optimal weights. We use the raw data of only one array for the covariance of the total signal, and the wavenumberfilte< red version of a full seismic singlesensor record for the covariance of the desired signal. In the determination of optimal weights, a parameter that controls the robustness of the beamformer against an arbitrary desired signal mismatch has to be chosen so that the results are optimal. This is similar to stabilization in deconvolution problems. This parameter needs to be smaller than the largest eigenvalue provided by the singular value decomposition of the presumed desired signal covariance. We compare results of MVDR beamforming with standard array forming on singlesensor synthetic and field seismic data. We apply 2D and 3D beamforming and show prestack and poststack results. MVDR beamformers are superior to conventional hardwired arrays for all examples.array signal processing; covariance analysis; geophysical prospecting; geophysical signal processing; seismology; signal denoising; singularalue decomposition)uuid:1873cc6f61604a1ba484b56a3b99ece9Dhttp://resolver.tudelft.nl/uuid:1873cc6f61604a1ba484b56a3b99ece9VAssessing estimated velocitydepth models: Finding error bars in tomographic inversion*Chitu, D.A.; AlAli, M.N.; Verschuur, D.J.xIn conventional migration velocity analysis methods, a velocity model is estimated that results in flattened events in commonimage gathers. However, after this process, no information is available on the accuracy of this velocity model. A statistical analysis of velocitymodel parameters is very difficult because of the integrated nature of the process. In commonfocuspoint technology, velocity estimation is split into two processes: a first step to estimate oneway focusing operators from the seismic data and a second step to translate these oneway propagation operators into a velocitydepth model. Because the second step does not involve seismic data and uses a handsoff model parameterization, a statistical analysis of the inversion result becomes rather straightforward. We developed a methodology for obtaining a suite of possible solutions, from which statistical measures can be extracted. By varying initial settings, the inversion of oneway traveltimes provides a space of solutions. Rather than having a single estimated model, we can obtain an ensemble of models. By performing statistical analysis of this ensemble, the error bars of the estimated velocity model can be retrieved. The procedure was tested for a 2D synthetic and field data set, for which the latter compares favorably to a conventional twoway traveltime tomography approach. The information provided by such an analysis is important because it shows the reliability of the final estimated model and could provide feedback for acquisition geometry design. More or better data might be needed to obtain a model to which a smaller degree of ambiguity is associated.mgeophysical techniques; inverse problems; seismic waves; seismology; statistical analysis; topography (Earth)8Electrical Engineering, Mathematics and Computer Science)uuid:1b484198cf8543508c306f271e6d9680Dhttp://resolver.tudelft.nl/uuid:1b484198cf8543508c306f271e6d9680CNearsurface attenuation estimation using wavepropagation modeling$El Yadari, N.; Ernst, F.; Mulder, W.GThe effect of the near surface on seismic land data can be so severe that static corrections are insufficient. Fullwaveform inversion followed by redatuming may be an alternative, but inversion will work only if the starting model is sufficiently close to the true model. As a first step toward determining a viscoelastic nearsurface model, we assume that existing methods can provide a horizontally layered velocity and density model. Because nearsurface attenuation is strongest, we propose a method to estimate the Pwave attenuation based on viscoacoustic finitedifference modeling. We compare energy decay along traveltime curves of reflection and refraction events in the modeled and observed seismic data for a range of attenuation parameters. The best match provides an estimate of the attenuation. First, we estimate only the attenuation of the top layer and study the sensitivity t< o depth and velocity perturbations. Then, we consider multiple layers. We apply the method to synthetic and real data and investigate the effect of source wavelet and topography. The method is robust against depth and velocity perturbations smaller than 10%. The results are sensitive to the source wavelet. Incorporating the surface topography in the computed traveltimes reduces the uncertainty of the attenuation estimates, especially for deeper layers.Zfinite difference methods; inverse problems; seismic waves; seismology; wavelet transforms)uuid:015bf386811f4587b50b9bac69f19699Dhttp://resolver.tudelft.nl/uuid:015bf386811f4587b50b9bac69f19699@Passive seismic interferometry by multidimensional deconvolution3Wapenaar, C.P.A.; Van der Neut, J.R.; Ruigrok, E.N.We introduce seismic interferometry of passive data by multidimensional deconvolution (MDD) as an alternative to the crosscorrelation method. Interferometry by MDD has the potential to correct for the effects of source irregularity, assuming the first arrival can be separated from the full response. MDD applications can range from reservoir imaging using microseismicity to crustal imaging with teleseismic data.Udeconvolution; geophysical techniques; multidimensional signal processing; seismology)uuid:14eeb991c4fc4959b8d9e15371864dc6Dhttp://resolver.tudelft.nl/uuid:14eeb991c4fc4959b8d9e15371864dc6~Stochastic joint inversion of 2D seismic and seismoelectric signals in linear poroelastic materials: A numerical investigation/Jardani, A.; Revil, A.; Slob, E.C.; Sllner, W.The interpretation of seismoelectrical signals is a difficult task because coseismic and seismoelectric converted signals are recorded simultaneously and the seismoelectric conversions are typically several orders of magnitude smaller than the coseismic electrical signals. The seismic and seismoelectric signals are modeled using a finiteelement code with perfectly matched layer boundary conditions assuming a linear poroelastic body. We present a stochastic joint inversion of the seismic and seismoelectrical data based on the adaptive Metropolis algorithm, to obtain the posterior probability density functions of the material properties of each geologic unit. This includes the permeability, porosity, electrical conductivity, bulk modulus of the dry porous frame, bulk modulus of the fluid, bulk modulus of the solid phase, and shear modulus of the formations. A test of this approach is performed with a synthetic model comprising two horizontal layers and a reservoir partially saturated with oil, which is embedded in the second layer. The result of the joint inversion shows that we can invert the permeability of the reservoir and its mechanical properties.elastic moduli; finite element analysis; geophysical prospecting; geophysical signal processing; hydrocarbon reservoirs; permeability; porosity; seismology; terrestrial electricity)uuid:b71c6656d7614fb4b0795c8f363383e9Dhttp://resolver.tudelft.nl/uuid:b71c6656d7614fb4b0795c8f363383e9_Raybased stochastic inversion of prestack seismic data for improved reservoir characterization.Van der Burg, D.; Verdel, A.; Wapenaar, C.P.A.Trace inversion for reservoir parameters is affected by angle averaging of seismic data and wavelet distortion on the migration image. In an alternative approach to stochastic trace inversion, the data are inverted prestack before migration using 3D dynamic ray tracing. This choice makes it possible to interweave trace inversion with Kirchhoff migration. The new method, called raybased stochastic inversion, is a generalization of current amplitude versus offset/amplitude versus angle (AVO/AVA) inversion techniques. The new method outperforms standard stochastic inversion techniques in cases of reservoir parameter estimation in a structurally complex subsurface with substantial lateral velocity variations and significant reflector dips. A simplification of the method inverts the normalincidence response from reservoirs with approximately planar layering at the subsurface target locations selected for inversion. It oper< ates along raypaths perpendicular to the reflectors, the direction that offers optimal resolution to discern layering in a reservoir. In a test on field data from the Gulf of Mexico, reservoir parameter estimates obtained with the simplified method, the estimates found by conventional stochastic inversion, and the actual values at a well drilled after the inversion are compared. Although the new method uses only 2% of the prestack data, the result indicates it improves accuracy on the dipping part of the reservoir, where conventional stochastic inversion suffers from wavelet stretch caused by migration._geophysical techniques; hydrocarbon reservoirs; seismic waves; seismology; stochastic processes)uuid:84703eaca0504e85bf3268bbae218732Dhttp://resolver.tudelft.nl/uuid:84703eaca0504e85bf3268bbae218732,Reflection images from ambient seismic noiseJDraganov, D.S.; Campman, X.; Thorbecke, J.W.; Verdel, A.; Wapenaar, C.P.A.`One application of seismic interferometry is to retrieve the impulse response (Green's function) from crosscorrelation of ambient seismic noise. Various researchers show results for retrieving the surfacewave part of the Green's function. However, reflection retrieval has proven more challenging. We crosscorrelate ambient seismic noise, recorded along eight parallel lines in the Sirte basin east of Ajdabeya, Libya, to obtain shot gathers that contain reflections. We take advantage of geophone groups to suppress part of the undesired surfacewave noise and apply frequencywavenumber filtering before crosscorrelation to suppress surface waves further. After comparing the retrieved results with data from an active seismic exploration survey along the same lines, we use the retrieved reflection data to obtain a migrated reflection image of the subsurface.dgeophysical signal processing; interference suppression; seismic waves; seismology; signal denoising)uuid:a2747c0b4ff741b291b19e4445c21b11Dhttp://resolver.tudelft.nl/uuid:a2747c0b4ff741b291b19e4445c21b11dEstimation of primaries and nearoffset reconstruction by sparse inversion: Marine data applications(Van Groenestijn, G.J.A.; Verschuur, D.J.ost waveequationbased multiple removal algorithms are based on prediction and subtraction of multiples. Especially for shallow water, the prediction strongly relies on a correct interpolation of the missing near offsets. The subtraction of predicted multiples from the data can easily lead to the distortion of primaries if primaries and multiples overlap. Recently, a new approach for surfacerelated multiple removal was proposed: the estimation of primaries by sparse inversion (EPSI), which is based on a full waveform inversion approach. EPSI is based on the same primarymultiple model as surfacerelated multiple elimination (SRME) and does not require a subsurface model. In contrast to SRME, EPSI estimates the primaries as unknowns in a multidimensional inversion process rather than a subtraction process.The multidimensional primary impulse responses are parameterized by bandlimited spikes, which are estimated such that they, along with their corresponding multiples, match the input data. An interesting aspect of the EPSI method is that it produces a residual, which is the part of the input data not explained by primaries and multiples. This residual can be analyzed and may provide useful information on the primary estimation process. Furthermore, it has been demonstrated that EPSI is also capable of reconstructing the missing near offsets from the multiples. The proposed method is applied to a field data set with moderate water depth, where it is demonstrated that the results are comparable with SRME. This data set is used to illustrate the residual. For a shallowwater field data set, it is shown that EPSI gives a better result than the standard SRME result caused by EPSI's capability to reconstruct the missing near offsets.1geophysical techniques; seismic waves; seismology)uuid:828be3c982fb4ec3b66385d09d0a130fDhttp://resolver.tudelft.nl/uuid:828be3c982fb4ec3b66385d09d0a130ffApplication of the finite< difference contrastsource inversion algorithm to seismic fullwaveform data7Abubakar, A.; Hu, W.; Habashy, T.M.; Van den Berg, P.M.We have applied the finitedifference contrastsource inversion (FDCSI) method to seismic fullwaveform inversion problems. The FDCSI method is an iterative nonlinear inversion algorithm. However, unlike the nonlinear conjugate gradient method and the GaussNewton method, FDCSI does not solve any full forward problem explicitly in each iterative step of the inversion process. This feature makes the method very efficient in solving largescale computational problems. It is shown that FDCSI, with a significant lower computation cost, can produce inversion results comparable in quality to those produced by the GaussNewton method and better than those produced by the nonlinear conjugate gradient method. Another attractive feature of the FDCSI method is that it is capable of employing an inhomogeneous background medium without any extra or special effort. This feature is useful when dealing with timelapse inversion problems where the objective is to reconstruct changes between the baseline and the monitor model. By using the baseline model as the background medium in crosswell seismic monitoring problems, high quality timelapse inversion results are obtained.Ofinite difference methods; geophysical techniques; inverse problems; seismology"IST/Imaging Science and Technology)uuid:9bf7fc0a55ba4cd091ddb51b2064b638Dhttp://resolver.tudelft.nl/uuid:9bf7fc0a55ba4cd091ddb51b2064b638jOn seismic interferometry, the generalized optical theorem, and the scattering matrix of a point scatterer)Wapenaar, C.P.A.; Slob, E.C.; Snieder, R.We have analyzed the farfield 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 threeterm 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 fourterm 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 fourterm 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 wellestablished 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; seismologyGeoscience & Engineering)uuid:869588aeb5e84f41b70af5b362e02d3eDhttp://resolver.tudelft.nl/uuid:869588aeb5e84f41b70af5b362e02d3eEEstimation of primaries by sparse inversion from passive seismic data"For passive seismic data, surface multiples are used to obtain an estimate of the subsurface responses, usually by a crosscorrelation process. This crosscorrelation process relies on the assumption that the surface has been uniformly illuminated by subsurface sources in terms of incident angles and strengths. If this is not the case, the crosscorrelation process cannot give a true amplitude estimation of the subsurface response. Furthermore, cross terms in the crosscorrelation result are not related to actual subsurface inhomogeneities. We have developed a meth< od that can obtain true amplitude subsurface responses without a uniform surfaceillumination assumption. Our methodology goes beyond the crosscorrelation process and estimates primaries only from the surfacerelated multiples in the available signal. We use the recently introduced estimation of primaries by sparse inversion (EPSI) methodology, in which the primary impulse responses are considered to be the unknowns in a largescale inversion process. With some modifications, the EPSI method can be used for passive seismic data. The output of this process is primary impulse responses with point sources and receivers at the surface, which can be used directly in traditional imaging schemes. The methodology was tested on 2D synthetic data."geophysical techniques; seismology)uuid:f8d8b93ed90c4d44addf1873edf600ffDhttp://resolver.tudelft.nl/uuid:f8d8b93ed90c4d44addf1873edf600ff8A perspective on 3D surfacerelated multiple elimination4Dragoset, B.; Verschuur, D.J.; Moore, I.; Bisley, R.lSurfacerelated multiple elimination (SRME) is an algorithm that predicts all surface multiples by a convolutional process applied to seismic field data. Only minimal preprocessing is required. Once predicted, the multiples are removed from the data by adaptive subtraction. Unlike other methods of multiple attenuation, SRME does not rely on assumptions or knowledge about the subsurface, nor does it use event properties to discriminate between multiples and primaries. In exchange for this freedom from the subsurface, SRME requires knowledge of the acquisition wavelet and a dense spatial distribution of sources and receivers. Although a 2D version of SRME sometimes suffices, most field data sets require 3D SRME for accurate multiple prediction. All implementations of 3D SRME face a serious challenge: The sparse spatial distribution of sources and receivers available in typical seismic field data sets does not conform to the algorithmic requirements. There are several approaches to implementing 3D SRME that address the data sparseness problem. Among those approaches are preSRME data interpolation, onthefly data interpolation, zeroazimuth SRME, and trueazimuth SRME. Field data examples confirm that (1) multiples predicted using trueazimuth 3D SRME are more accurate than those using zeroazimuth 3D SRME and (2) onthefly interpolation produces excellent results.)geophysical signal processing; seismology)uuid:a751046354464a49a8fc81f44db1d984Dhttp://resolver.tudelft.nl/uuid:a751046354464a49a8fc81f44db1d984NTutorial on seismic interferometry: Part 1 Basic principles and applicationsFWapenaar, C.P.A.; Draganov, D.S.; Snieder, R.; Campman, X.; Verdel, A.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 xaxis, the crosscorrelation of the responses at two receivers along the xaxis 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. Directwave 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 directwave interferometry is theretrieval of seismic surfacewave responses from ambient noise and the subsequent tomographic determination of the surfacewave 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 planewave reflection response of that medium. This is essentially 1D< reflectedwave 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 reflectedwave interferometry is retrieving the seismic reflection response from ambient noise and imaging of the reflectors in the subsurface. A common aspect of direct and reflectedwave 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.)uuid:07504f32d9fb46b98095dcfa5b3e817bDhttp://resolver.tudelft.nl/uuid:07504f32d9fb46b98095dcfa5b3e817bOTutorial on seismic interferometry: Part 2 Underlying theory and new advances5Wapenaar, C.P.A.; Slob, E.C.; Snieder, R.; Curtis, A.`In the 1990s, the method of timereversed 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 timereversed responses are not physically sent into the earth, but they are convolved with other measured responses. The effect is essentially the same: The timereversed 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 different 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 tracebytrace deconvolution process compensates for complex source functions and the attenuation of the medium. Interferometry by multidimensional deconvolution also compensates for the effects of onesided and/or irregular illumination.jdeconvolution; geophysical techniques; Green's function methods; interferometry; seismic waves; seismology)uuid:bffc466e307347c0a35d32c6366ae2bfDhttp://resolver.tudelft.nl/uuid:bffc466e307347c0a35d32c6366ae2bf@Highresolution lithospheric imaging with seismic interferometry8Ruigrok, E.N.; Campman, X.; Draganov, D.S.; Wapenaar, K.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 industrystyle imaging algorithms applied to converted and scatteredwave energy from distant earthquakes, leading to relatively highresolution 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 dataprocessing me< thods. 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 sourceside reverberation characteristics as field data recorded during the 2000 2001 Laramie broadband experiment. Finally, we apply our method to the Laramie data, retrieving reflection data. We extract a 2D 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.9seismology; interferometry; body waves; crustal structure)uuid:c43b79fe17134fdca310ee70325cfaffDhttp://resolver.tudelft.nl/uuid:c43b79fe17134fdca310ee70325cfaffHighresolution reservoir characterization by an acoustic impedance inversion of a Tertiary deltaic clinoform system in the North Sea>Tetyukhina, D.; Van Vliet, L.J.; Luthi, S.M.; Wapenaar, C.P.A.IFluviodeltaic sedimentary systems are of great interest for explorationists because they can form prolific hydrocarbon plays. However, they are also among the most complex and heterogeneous ones encountered in the subsurface, and potential reservoir units are often close to or below seismic resolution. For seismic inversion, it is therefore important to integrate the seismic data with higher resolution constraints obtained from well logs, whereby not only the acoustic properties are used but also the detailed layering characteristics. We have applied two inversion approaches for poststack, timemigrated seismic data to a clinoform sequence in the North Sea. Both methods are recursive tracebased techniques that use well data as a priori constraints but differ in the way they incorporate structural information. One method uses a discrete layer model from the well that is propagated laterally along the clinoform layers, which are modeled as sigmoids. The second method uses a constant sampling rate from the well data and uses horizontal and vertical regularization parameters for lateral propagation. The first method has a low level of parameterization embedded in a geologic framework and is computationally fast. The second method has a much higher degree of parameterization but is flexible enough to detect deviations in the geologic settings of the reservoir; however, there is no explicit geologic significance and the method is computationally much less efficient. Forward seismic modeling of the two inversion results indicates a good match of both methods with the actual seismic data.Ngeology; geophysical techniques; hydrocarbon reservoirs; sediments; seismology)uuid:fdd34e7403a7458fbb8940640a312e74Dhttp://resolver.tudelft.nl/uuid:fdd34e7403a7458fbb8940640a312e74tMethodology for dense spatial sampling of multicomponent recording of converted waves in shallow marine environments7El Allouche, N.; Drijkoningen, G.G.; Van der Neut, J.R.NA widespread use of converted waves for shallow marine applications is hampered by spatial aliasing and field efficiency. Their short wavelengths require dense spatial sampling which often needs to be achieved by receivers deployed on the seabed. We adopted a new methodology where the dense spatial sampling is achieved in the commonreceiver domain by reducing the shot spacing. This is done by shooting one track multiple times and merging the shot lines in an effective manner in a separate processing step. This processing step is essential because positioning errors< introduced during the field measurement can become significant in the combined line, particularly when they exceed the distance between two adjacent shot positions. For this processing step, a particular shot line is used as a reference line and relative variations in source and receiver positions in the other shot lines are corrected for using crosscorrelation. The combined shot line can subsequently be regularized for further processing. The methodology is adopted in a field experiment conducted in the Danube River in Hungary. The aim of the seismic experiment was to acquire properly sampled convertedwave data using a multicomponent receiver array. The dense spatial sampling was achieved by sailing one track 14 times. After correcting for the underwater receiver positions using the direct arrival, the crosscorrelation step was applied to merge the different shot lines. The successfully combined result is regularized into a densely sampled data set free of visible spatial aliasing and suitable for convertedwave processing.9geophysical techniques; rivers; seismic waves; seismology)uuid:fd60045ded3b487fa251085255996931Dhttp://resolver.tudelft.nl/uuid:fd60045ded3b487fa251085255996931RConverted waves in a shallow marine environment: Experimental and modeling studies<El Allouche, N.; Drijkoningen, G.G.; Versteeg, W.; Ghose, R.Seismic waves converted from compressional to shear mode in the shallow subsurface can be useful not only for obtaining shearwave velocity information but also for improved processing of deeper reflection data. These waves generated at deep seas have been used successfully in hydrocarbon exploration; however, acquisition of goodquality convertedwave data in shallow marine environments remains challenging. We have looked into this problem through field experiments and synthetic modeling. A highresolution seismic survey was conducted in a shallowwater canal using different types of seismic sources; data were recorded with a fourcomponent waterbottom cable. Observed events in the field data were validated through modeling studies. Compressional waves converted to shear waves at the water bottom and at shallow reflectors were identified. The shear waves showed distinct linear polarization in the horizontal plane and low velocities in the marine sediments. Modeling results indicated the presence of a nongeometric shearwave arrival excited only when the dominant wavelength exceeded the height of the source with respect to the water/sediment interface, as observed in airgun data. This type of shear wave has a traveltime that corresponds to the raypath originating not at the source but at the interface directly below the source. Thus, these shear waves, excited by the source/waterbottom coupled system, kinematically behave as if they were generated by an Swave source placed at the water bottom. In a shallowwater environment, the condition appears to be favorable for exciting such shear waves with nongeometric arrivals. These waves can provide useful information of shearwave velocity in the sediments.gdata acquisition; geophysical prospecting; hydrocarbon reservoirs; sediments; seismic waves; seismology)uuid:ee28ded7276045178a38b19a52d6b918Dhttp://resolver.tudelft.nl/uuid:ee28ded7276045178a38b19a52d6b918RSeismic migration of blended shot records with surfacerelated multiple scatteringVerschuur, D.J.; Berkhout, A.J.This paper focuses on the concept of using blended data and multiple scattering directly in the migration process, meaning that the blended input data for the proposed migration algorithm includes blended surfacerelated multiples. It also means that both primary and multiple scattering contribute to the seismic image of the subsurface. Essential in our approach is that multiples are not included in the Green's functions but are part of the incident wavefields, utilizing the socalled double illumination property. We find that complex incident wavefields, such as blended primaries and/or blended multiples, require a reformulation of the imaging principle in order to provide < broadband angledependent reflection properties.Delft University of Technology)uuid:c4772e2ccaae4123b4afc462f23f3489Dhttp://resolver.tudelft.nl/uuid:c4772e2ccaae4123b4afc462f23f3489Seismic interferometry using multidimensional deconvolution and crosscorrelation for crosswell seismic reflection data without borehole sourcesUMinato, S.; Matsuoka, T.; Tsuji, T.; Draganov, D.S.; Hunziker, J.W.; Wapenaar, C.P.A.Crosswell reflection method is a highresolution seismic imaging method that uses recordings between boreholes. The need for downhole sources is a restrictive factor in its application, for example, to timelapse 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 illposed, and propose to stabilize it using singularvalue 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 seismicreflection data using borehole sources and with the logged Pwave velocity.Ageophysical techniques; interferometry; seismic waves; seismology)uuid:853667221d954c9e9ff5de071270e49fDhttp://resolver.tudelft.nl/uuid:853667221d954c9e9ff5de071270e49f\Estimation of changes in saturation and pressure from 4D seismic AVO and timeshift analysis2Trani, M.; Arts, R.; Leeuwenburgh, O.; Brouwer, J.XA reliable estimate of reservoir pressure and fluid saturation changes from timelapse seismic data is difficult to obtain. Existing methods generally suffer from leakage between the estimated parameters. We propose a new method using different combinations of timelapse seismic attributes based on four equations: two expressing changes in prestack AVO attributes (zerooffset and gradient reflectivities), and two expressing poststack timeshifts of compressional and shear waves as functions of productioninduced changes in fluid properties. The effect of using different approximations of these equations was tested on a realistic, synthetic reservoir, where seismic data have been simulated during the 30year lifetime of a waterflooded oil reservoir. Results found the importance of the porosity in the inversion with a clear attenuation of the porosity imprint on the final estimates in case the porosity field or the vertically averaged porosity field is known a priori. The use of a firstorder approximation of the gradient reflectivity equation leads to severely biased estimates of changes in saturation and leakage between the two different parameters. Both the bias and the leakage can be reduced, if not eliminated, by including higherorder terms in the description of the gradient, or by replacing the gradient equation with P and/or Swave timeshift data. The final estimates are relatively robust to random noise, as they present fairly high accuracy in the presence of white noise with a standard deviation of 15%. The introduction of systematic noise decreases the inversion accuracy more severely.Vgeophysical techniques; hydrocarbon reservoirs; seismic waves; seismology; white noise)uuid:dc1eb3721686416c89df29bd007da4dfDhttp://resolver.tudelft.nl/uuid:dc1eb3721686416c89df29bd007da4dfpSource depopulation potential and surfacewave tomography using a crosscorrelation method in a scatte< ring mediumPGouedard, P.; Roux, P.; Campillo, M.; Verdel, A.R.; Yao, H.; Van der Hilst, R.D.kWe use seismic prospecting data on a 40 40 regular grid of sources and receivers deployed on a 1 km 1 km area to assess the feasibility and advantages of velocity analysis of the shallow subsurface by means of surfacewave tomography with Green's functions estimated from crosscorrelation. In a first application we measure Rayleighwave dispersion curves in a 1D equivalent medium. The assumption that the medium is laterally homogeneous allows using a simple projection scheme and averaging of crosscorrelation functions over the whole network. Because averaging suppresses noise, this method yields better signaltonoise ratio than traditional activesource approaches, and the improvement can be estimated a priori from acquisition parameters. We find that highquality dispersion curves can be obtained even when we reduce the number of active sources used as input for the correlations. Such source depopulation can achieve significant reduction in the cost of active source acquisition. In a second application we compare Rayleighwave group velocity tomography from raw and reconstructed data. We can demonstrate that the crosscorrelation approach yields group velocity maps that are similar to active source maps. Scattering has an importance here as it may enhance the crosscorrelation performance. We quantify the scattering properties of the medium using mean free path measurements from coherent and incoherent parts of the signal. We conclude that for firstorder velocity analysis of the shallow subsurface, the use of crosscorrelation offers a costeffective alternative to methods that rely exclusively on active sources.correlation methods; geophysical prospecting; geophysical signal processing; Green's function methods; Rayleigh waves; seismology; tomography)uuid:ffc0289e5fdc408b9c8d8cc8e14ecdbdDhttp://resolver.tudelft.nl/uuid:ffc0289e5fdc408b9c8d8cc8e14ecdbdoLaboratory measurements and theoretical modeling of seismoelectric interface response and coseismic wave fields<Schakel, M.D.; Smeulders, D.M.J.; Slob, E.C.; Heller, H.K.J.wA fullwaveform seismoelectric numerical model incorporating the directivity pattern of a pressure source is developed. This model provides predictions of coseismic electric fields and the electromagnetic waves that originate from a fluid/porousmedium interface. An experimental setup in which coseismic electric fields and interface responses are measured is constructed. The seismoelectric origin of the signals is confirmed. The numerically predicted polarity reversal of the interfacial signal and seismoelectric effects due to multiple scattering are detected in the measurements. Both the simulated coseismic electric fields and the electromagnetic waves originating from interfaces agree with the measurements in terms of travel times, waveform, polarity, amplitude, and spatial amplitude decay, demonstrating that seismoelectric effects are comprehensively described by theory.Jgeophysical techniques; seismic waves; seismology; terrestrial electricityAmerican Institute of Physics)uuid:1d6d4d70a4584e3b954c969cef7dc2e8Dhttp://resolver.tudelft.nl/uuid:1d6d4d70a4584e3b954c969cef7dc2e8aSeparation of blended data by iterative estimation and subtraction of blending interference noise)Mahdad, A.; Doulgeris, P.; Blacquiere, G.vSeismic acquisition is a tradeoff between economy and quality. In conventional acquisition the time intervals between successive records are large enough to avoid interference in time. To obtain an efficient survey, the spatial source sampling is therefore often (too) large. However, in blending, or simultaneous acquisition, temporal overlap between shot records is allowed. This additional degree of freedom in survey design significantly improves the quality or the economics or both. Deblending is the procedure of recovering the data as if they were acquired in the conventional, unblended way. A simple leastsquares procedure, however, does not remove the interference due to other source< s, or blending noise. Fortunately, the character of this noise is different in different domains, e.g., it is coherent in the common source domain, but incoherent in the common receiver domain. This property is used to obtain a considerable improvement. We propose to estimate the blending noise and subtract it from the blended data. The estimate does not need to be perfect because our procedure is iterative. Starting with the leastsquares deblended data, the estimate of the blending noise is obtained via the following steps: sort the data to a domain where the blending noise is incoherent; apply a noise suppression filter; apply a threshold to remove the remaining noise, ending up with (part of) the signal; compute an estimate of the blending noise from this signal. At each iteration, the threshold can be lowered and more of the signal is recovered. Promising results were obtained with a simple implementation of this method for both impulsive and vibratory sources. Undoubtedly, in the future algorithms will be developed for the direct processing of blended data. However, currently a highquality deblending procedure is an important step allowing the application of contemporary processing flows~data acquisition; geophysical signal processing; iterative methods; least squares approximations; seismology; signal denoising)uuid:55fabcd8043548baaba0a0bad1e05033Dhttp://resolver.tudelft.nl/uuid:55fabcd8043548baaba0a0bad1e05033tControlledsource interferometric redatuming by crosscorrelation and multidimensional deconvolution in elastic mediaLVan der Neut, J.R.; Thorbecke, J.W.; Mehta, K.; Slob, E.C.; Wapenaar, C.P.A.{Various researchers have shown that accurate redatuming of controlled seismic sources to downhole receiver locations can be achieved without requiring a velocity model. By placing receivers in a horizontal or deviated well and turning them into virtual sources, accurate images can be obtained even below a complex nearsubsurface. Examples include controlledsource interferometry and the virtualsource method, both based on crosscorrelated signals at two downhole receiver locations, stacked over source locations at the surface. Because the required redatuming operators are taken directly from the data, even multiple scattered waveforms can be focused at the virtualsource location, and accurate redatuming can be achieved. To reach such precision in a solid earth, representations for elastic wave propagation that require multicomponent sources and receivers must be implemented. Wavefield decomposition prior to crosscorrelation allows us to enforce virtual sources to radiate only downward or only upward. Virtualsource focusing and undesired multiples from the overburden can be diagnosed with the interferometric pointspread function (PSF), which can be obtained directly from the data if an array of subsurface receivers is deployed. The quality of retrieved responses can be improved by filtering with the inverse of the PSF, a methodology referred to as multidimensional deconvolution.acoustic wave interferometry; correlation methods; deconvolution; filtering theory; geophysical signal processing; geophysical techniques; seismic waves; seismology)uuid:00af288facaf48ecac90812f1c8c4988Dhttp://resolver.tudelft.nl/uuid:00af288facaf48ecac90812f1c8c4988ISeismoelectric interface response: Experimental results and forward modelUnderstanding the seismoelectric interface response is important for developing seismoelectric field methods for oil exploration and environmental/engineering geophysics. The existing seismoelectric theory has never been validated systematically by controlled experiments. We have designed and developed an experimental setup in which acoustictoelectromagnetic wave conversions at interfaces are measured. An acoustic source emits a pressure wave that impinges upon a porous sample. The reflected electricwave potential is recorded by a wire electrode. We have also developed a fullwaveform electrokinetic theoretical model based on the Sommerfeld approach and have compared it with measurements at pos< itions perpendicular and parallel to the fluid/porousmedium interface. We performed experiments at several salinities. For 103 and 102 M sodium chloride (NaCl) solutions, both waveforms and amplitudes agree. For 104 M NaCl, however, amplitude deviations occur. We found that a single amplitude field scaling factor describes these discrepancies. We also checked the repeatability of experiments. The amplitudes are constant for the duration of an experiment (1 4 hours) but decrease on longer time scales (~24 hours). However, the waveforms and spatial amplitude pattern of the electric wavefield are preserved over time. Our results validate electrokinetic theory for the seismictoelectromagneticwave conversion at interfaces for subsurface exploration purposes.jacoustoelectric effects; geophysical prospecting; geophysical signal processing; seismic waves; seismology)uuid:a16f2acac0454894bdcd2ec6dfddc75cDhttp://resolver.tudelft.nl/uuid:a16f2acac0454894bdcd2ec6dfddc75cSBodywave seismic interferometry applied to earthquake and storminduced wavefield
Ruigrok, E.N.Wapenaar, C.P.A. (promotor)
&Seismology is the study of the vibration of the Earth. Seismologists pay much attention to the main source of Earth vibration: earthquakes. But also other seismic sources, like mining blasts, ocean storms and windmills, are studied. All these sources induce seismic waves, which can eventually be recorded as ground vibrations. These seismic records contain not only information about the sources, but also about the part of the Earth through which the waves have propagated. This thesis focuses on a main subclass of seismology: seismic imaging. With seismic imaging, seismic records are studied with the aim to unravel the structure and composition of the Earth. Seismic imaging finds its main use in extracting information about the solid Earth that is within our reach. With the current state of the art, we can only mine minerals from the upper 0.2 percent. Imaging the deep (i.e., unreachable) Earth finds its main use in hazard assessment. Through imaging the deep Earth, we can better understand the dynamics of our planet, manifested in, amongst others, earthquakes and volcanism. By forecasting how the Earth will reshape, we can sensibly adapt. Seismic characterization of the subsurface also improves our ability to assess the direct impact of earthquakes. A large part of the current deepEarth images are obtained through tomographic inversion. These images are successful for identifying largescale structural anomalies, like (remnants of) subducting slabs. However, they lack the resolution to accurately image discontinuities in space. Sharp discontinuities, like the interface between sedimentary and crystalline rock, lead to conversions and reflections of seismic waves. It is those reflections that are most suitable for imaging the discontinuities. The technique developed to do this is called seismic reflection imaging (SRI). Though SRI is considered to be an advanced technique for imaging the subsurface, the dense and wide source (and receiver) distribution requirement thus far limited the use for deepEarth imaging. DeepEarth imaging is done primarily with natural sources. Their limited and irregular distribution and their unpredictable occurrence time is inconvenient for SRI. In this thesis we work on a methodology to alleviate the stringent requirement on the source distribution. This methodology is called seismic interferometry (SI). SI is a remapping operation. With SI, the responses from many sources are combined to create the response as if there were a source at a receiver position. The response from this "virtual" source could contain both surface waves and body waves. We apply SI such that especially (reflected) body waves are retrieved, because it are these waves that are used for SRI. When we apply SI to a regularlyspaced array of receivers, we could turn uncontrolled natural sources into a wellorganized succession of virtual sources, which is amenable for SRI. The aim of this work is to improve the imaging of the interior of the < Earth by the application of bodywave SI to naturally induced seismic data. A prior study focused on the explorationscale (i.e., the part of the Earth that is being mined). In this thesis, we focus on lowfrequency waves (< 2 Hz), which contain especially information about the deep Earth. In Part I we consider applications to earthquake recordings. In Part II we evaluate the retrieval of reflections using microseisms. Part I starts by looking into the globalscale configuration. For this configuration, two seismic stations (or two arrays of stations) may be located anywhere on the globe and reflectors may be located anywhere in depth. We derive and numerically test relations to retrieve the complete response between two stations, using worldwide seismicity. The SI relations are valid for a closed entity for seismic waves, like the Earth is by approximation. The relations do not account for inelastic losses, by which eventually all seismic energy is transformed to heat. Though normally of negligible influence, the SI relations do account for the rotation of the Earth. In practice, the distribution of the larger earthquakes is not wide enough to retrieve complete responses between any two points on the Earth. However, still relevant primary reflections can be retrieved with a limited dense distribution of seismicity, like from the Tohoku (2011) aftershock sequence. Part I proceeds by zooming in to the lithosphericscale configuration. The lithosphere is the part of the solid Earth where most (seismic) activity takes place. It makes up the crust and the upper part of the mantle. It is this part of the Earth that, broken into plates and fueled by convection in the mantle, undergoes birth (at midoceanic ridges), collisions (resulting in stunning topography) and death (subduction back into the mantle). Inspecting the Earth in depth unveils a part of the history of the lithosphere and gives an indication of the processes still ongoing. Information about the structure of the lithosphere is hidden in the coda of wavefields that arrive due to distant seismicity. We test a few different SI approaches to extract this information. We use synthetic data with similar characteristics as field data recorded during the Laramie broadband experiment (20002001). This was an array of seismic stations in Wyoming, USA, to study an suture zone. We evaluate the requirements for obtaining multioffset reflection responses. These multioffset reflections are important for obtaining a velocity model. We estimate this model for the subsurface below the Laramie array and we use it to map the extracted reflections to a reflectivity image of the lithosphere. In the remainder of part I, we focus again on the lithospheric scale. This time, we assume that a velocity model is already available through other means. Consequently, only singleoffset reflection responses are required to make an image. We use global phases to obtain zerooffset reflection responses. We show the robustness of the method with data from the HiCLIMB experiment (20022005). This was a large and wellsampled seismic array, passing the Himalayas and a significant part of the Tibetan Plateau. The successive application of SRI leads to an image of the IndianEurasian Plate collision. The rediscovery of body waves in lowfrequency noise (<1 Hz) opened up the research for Part II of this thesis. The noise are in fact microseisms: Earth vibrations that are indirectly caused by ocean gravity waves. Highamplitude ocean waves are caused when storms, leading to persistent wind fields, cross an ocean. We retrieve reflection responses from bodywave noise using SI approaches similar to the ones used for earthquake recordings. The main difference between applications to noise and (large) earthquake responses is that the origin of large earthquakes is known, whereas the origin of noise is generally not known. Therefore, the main challenge for noise applications is to unravel the noise illumination. For this purpose, we use a wellsampled areal array in Egypt to study the noise illumination in different frequ<jency bands. We only select the noise records with a favorable bodywave content and process them into separate reflection responses, of both the lithosphere and the upper crust. We further evaluate what basinscale information can be extracted from the retrieved reflections. The wavelengths of the reflections extracted from lowfrequency noise are too long to image a sedimentary basin in detail. However, the sharp boundary between the sediments and the underlying crystalline rocks can be delineated. For the same dataset from Egypt, we compare the information that is extracted with SI with the information that is obtained when using two other passive seismic techniques, namely horizontaltovertical spectral ratio and receiver function. In the appendices we discuss a few spinoff developments from the main research. In Appendix A we work out an alternative SI relation that integrates over midpoints between receiver positions, instead of over sources. Using this relation, reflections can already be retrieved using only a single source. However, a wellsampled array of receivers needs to be available and an additional evaluation needs to be performed. In Appendix B we develop a technique to estimate loss factors from planewave transmission responses. The losses are estimated from amplitude ratios, obtained before and after applying autocorrelation. In Appendix C we present a method that was inspired by SI. Crosscorrelations of diffractions from the coremantle boundary are used to swiftly scan anomalies in the crust and mantle. Finally, in Appendix D we discuss an alternative implementation of SI: SI by multidimensional deconvolution. With this technique the outcome of SI may be improved. However, wellsampled arrays of receivers are required and additional processing needs to be implemented. In this thesis, we show that the application of SI to earthquake recordings is sufficiently mature to yield highresolution images of the lithosphere. Relations for imaging deeper structures are worked out and modus operandi are thought through. However, the reality check is still to be encountered with fielddata applications. Furthermore, we show that the applications of SI to microseisms can be used to yield information about a basin depth. Moreover, we show the great promise of microseism applications for unveiling a lithosphere in depth, using only a few days of noise recordings. The more repetitive character of microseism sources with respect to earthquakes, in fact makes them more amenable for monitoring applications, if their yearly variations and radiation characteristics are well understood._seismic interferometry; seismic imaging; seismology; microseism; lithosphere; sedimentary basinBOXPress
20121001)uuid:50b487928ebf456380d80b70b151318eDhttp://resolver.tudelft.nl/uuid:50b487928ebf456380d80b70b151318evEvanescent wave coupling in a geophysical system: Airborne acoustic signals from the Mw 8.1 Macquarie Ridge earthquakeNEvers, L.G.; Brown, D.; Heaney, K.D.; Assink, J.D.; Smets, P.S.M.; Snellen, M.bAtmospheric lowfrequency sound, i.e., infrasound, from underwater events has not been considered thus far, due to the high impedance contrast of the waterair interface making it almost fully reflective. Here we report for the first time on atmospheric infrasound from a large underwater earthquake (Mw 8.1) near the Macquarie Ridge, which was recorded at 1325 km from the epicenter. Seismic waves coupled to hydroacoustic waves at the ocean floor, after which the energy entered the Sound Fixing and Ranging channel and was detected on a hydrophone array. The energy was diffracted by a seamount and an oceanic ridge, which acted as a secondary source, into the water column followed by coupling into the atmosphere. The latter results from evanescent wave coupling and the attendant anomalous transparency of the sea surface for very low frequency acoustic waves.2earthquake; seismology; infrasound; hydroacousticsAmerican Geophysical Union
20140910
*+&ffffff?'ffffff?(?)?"dXX333333?333333?U}}}}}}}}}} }
}}}
}}}}}}}}}}}}
@
!
"
#
$
@
%
&
'
(
)
@
*
+
,

.
,@
/
0
1
2
3
,@
4
5
6
7
8T@
9
:
;
<
=
>
?
@
A
BX@
C
D
;
E
=
>
F
G
H
IX@
J
K
;
E
=
>
L
M
N
O X@
P
Q
;
E
=
>
R
S
T
U
X@
V
W
;
E
=
>
X
Y
Z
[X@
\
]
;
E
=
>
^
_
`
aX@
b
]
;
E
=
>
c
d
e
f
X@
g
h
;
E
=
>
i
j
k
lX@
m
n
;
E
=
o
p
q
fX@
r
s
;
E
=
>
t
u
v
wX@
x
y
;
E
z
{

}X@
~
;
E
=
X@
;
E
=
>
X@
;
E
=
>
\@
;
E
=
>
\@
;
E
=
>
\@
;
E
=
>
`@
;
E
=
>
`@
;
=
>
`@
;
E
=
`@
;
E
=
>
`@
;
E
=
>
`@
;
E
`@
;
E
=
>
`@
;
E
=
>
d@
;
E
=
>
d@
;
E
=
>
!
!
!
!!d@
!
!
!
! ;
!
E
!=
!>
"
"
"
""d@
"
"
"
" ;
"
E
"
"
#
#
#
##d@
#
#
#
# ;
#
E
#
#
$
$
$
$$h@
$
$
$
$ ;
$
E
$=
$
%
%
%
%%h@
%
%
%
% ;
%
E
%
%
&
&
&
&&h@
&
&
&
& ;
&
E
&
&
'
'
'
''h@
'
'
'
' ;
'
E
'=
'
(
(
(
(
(h@
(
(
(
( ;
(
E
(=
(
)
)
)
))h@
)
)
)
) ;
)
)=
)>
*
*
*
**h@
*
*
*
* ;
*
E
*=
*
+
+
+
++h@
+
+
+
+ ;
+
E
+=
+
,
,
,!
,",l@
,#
,$
,
, ;
,
E
,=
,
%
&
'
(l@
)


 ;

E
*
.+
.,
.
...l@
./
.0
.
. ;
.
E
.=
.
/1
/2
/3
/4/l@
/5
/6
/
/ ;
/
E
/=
/
07
08
09
0:0l@
0;
0<
0
0 ;
0
E
0=
0
1=
1>
1?
1@1l@
1A
1B
1
1 ;
1
C
1=
1>
2D
2E
2F
2G2l@
2H
2I
2
2 ;
2
E
2=
2
3J
3K
3L
3M3l@
3N
3O
3
3 ;
3
E
3=
3
4P
4Q
4R
4@4l@
4S
4T
4
4 ;
4
E
4=
4
5U
5V
5W
5X
5Y5p@
5Z
5[
5
5
5
\
5]
5=
5
6^
6_
6`
6a6x@
6b
6c
6
6 ;
6
d
6e
6=
6789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~>@dd666yKyKhttp://resolver.tudelft.nl/uuid:b959853fcc3443e7add28e2c68620e16yKyKhttp://resolver.tudelft.nl/uuid:09df29e294874fe687e0c2a9470933d7yKyKhttp://resolver.tudelft.nl/uuid:1b7ad6beb31e4b0888cd201049bbf8feyKyKhttp://resolver.tudelft.nl/uuid:10ed7f4d81574cdba97621dc65175ee5yKyKhttp://resolver.tudelft.nl/uuid:15232dc10d984ba2bdaef59f520134b4yKyKhttp://resolver.tudelft.nl/uuid:c3fd77aa1a0f4a7192d286a722ed1366yKyKhttp://resolver.tudelft.nl/uuid:2933280d74b44bbbb46d77b5b26a59aayKyKhttp://resolver.tudelft.nl/uuid:eefac5acf44046808978a26848f2230f yKyKhttp://resolver.tudelft.nl/uuid:85b144578770463aaf0cff2243ac1b01
yKyKhttp://resolver.tudelft.nl/uuid:c0512fe5692e47e398a2a24cf29c0d09yKyKhttp://resolver.tudelft.nl/uuid:fc9a5a03cbfa40ca8f4d9652ecd325f5yKyKhttp://resolver.tudelft.nl/uuid:68e13eb252e7499c8a04ddee6fa0d6dd
yKyKhttp://resolver.tudelft.nl/uuid:4e602a55c4b647c78d81bce52c48df6fyKyKhttp://resolver.tudelft.nl/uuid:bc09960947bd4ad891ad2a4a9a949c0fyKyKhttp://resolver.tudelft.nl/uuid:a3b3a766c9214ac9826276ce9b68cf6fyKyKhttp://resolver.tudelft.nl/uuid:97664d6d4fb1494aa99f16f75e5f9e9dyKyKhttp://resolver.tudelft.nl/uuid:6d29c07992eb43e0b1c478b461598ebfyKyKhttp://resolver.tudelft.nl/uuid:78eea899c7f443839eb89eac4bf7ec0eyKyKhttp://resolver.tudelft.nl/uuid:4c09ff8864be40529645bd08a02b1424yKyKhttp://resolver.tudelft.nl/uuid:389fad2146c34368bf75ea6f7b1f8586yKyKhttp://resolver.tudelft.nl/uuid:d318d89877994c6ab3393fde71306c9eyKyKhttp://resolver.tudelft.nl/uuid:e599c89e0186472394b5364e85ddb0afyKyKhttp://resolver.tudelft.nl/uuid:b68278f691fb471b9d49decc2927f766yKyKhttp://resolver.tudelft.nl/uuid:7d656a13c0f1425a8835616a1ccfdb0ayKyKhttp://resolver.tudelft.nl/uuid:cacc1e07acd743a68d39b16ad3a3abf0yKyKhttp://resolver.tudelft.nl/uuid:b2b336646da1418f9af45ab98fdd439dyKyKhttp://resolver.tudelft.nl/uuid:4ce33f95c28047f2b443d39eb24f7eeayKyKhttp://resolver.tudelft.nl/uuid:1873cc6f61604a1ba484b56a3b99ece9yKyKhttp://resolver.tudelft.nl/uuid:1b484198cf8543508c306f271e6d9680yKyKhttp://resolver.tudelft.nl/uuid:015bf386811f4587b50b9bac69f19699yKyKhttp://resolver.tudelft.nl/uuid:14eeb991c4fc4959b8d9e15371864dc6 yKyKhttp://resolver.tudelft.nl/uuid:b71c6656d7614fb4b0795c8f363383e9!!yKyKhttp://resolver.tudelft.nl/uuid:84703eaca0504e85bf3268bbae218732""yKyKhttp://resolver.tudelft.nl/uuid:a2747c0b4ff741b291b19e4445c21b11##yKyKhttp://resolver.tudelft.nl/uuid:828be3c982fb4ec3b66385d09d0a130f$$yKyKhttp://resolver.tudelft.nl/uuid:9bf7fc0a55ba4cd091ddb51b2064b638%%yKyKhttp://resolver.tudelft.nl/uuid:869588aeb5e84f41b70af5b362e02d3e&&yKyKhttp://resolver.tudelft.nl/uuid:f8d8b93ed90c4d44addf1873edf600ff''yKyKhttp://resolver.tudelft.nl/uuid:a751046354464a49a8fc81f44db1d984((yKyKhttp://resolver.tudelft.nl/uuid:07504f32d9fb46b98095dcfa5b3e817b))yKyKhttp://resolver.tudelft.nl/uuid:bffc466e307347c0a35d32c6366ae2bf**yKyKhttp://resolver.tudelft.nl/uuid:c43b79fe17134fdca310ee70325cfaff++yKyKhttp://resolver.tudelft.nl/uuid:fdd34e7403a7458fbb8940640a312e74,,yKyKhttp://resolver.tudelft.nl/uuid:fd60045ded3b487fa251085255996931yKyKhttp://resolver.tudelft.nl/uuid:ee28ded7276045178a38b19a52d6b918..yKyKhttp://resolver.tudelft.nl/uuid:c4772e2ccaae4123b4afc462f23f3489//yKyKhttp://resolver.tudelft.nl/uuid:853667221d954c9e9ff5de071270e49f00yKyKhttp://resolver.tudelft.nl/uuid:dc1eb3721686416c89df29bd007da4df11yKyKhttp://resolver.tudelft.nl/uuid:ffc0289e5fdc408b9c8d8cc8e14ecdbd22yKyKhttp://resolver.tudelft.nl/uuid:1d6d4d70a4584e3b954c969cef7dc2e833yKyKhttp://resolver.tudelft.nl/uuid:55fabcd8043548baaba0a0bad1e0503344yKyKhttp://resolver.tudelft.nl/uuid:00af288facaf48ecac90812f1c8c498855yKyKhttp://resolver.tudelft.nl/uuid:a16f2acac0454894bdcd2ec6dfddc75c66yKyKhttp://resolver.tudelft.nl/uuid:50b487928ebf456380d80b70b151318egg
Root Entry F[^\[^\@SummaryInformation( F<Workbook FDocumentSummaryInformation8 F
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~
!"#$%&'()*+,./0123456789:;<=>?@ABCDEFGHIJKLMNOPQRSTUVWXYZ[\]^_`abcdefghijklmnopqrstuvwxyz{}~