scholarly journals Removal of surface‐related wave phenomena—The marine case

Geophysics ◽  
1996 ◽  
Vol 61 (1) ◽  
pp. 202-210 ◽  
Author(s):  
R. G. van Borselen ◽  
J. T. Fokkema ◽  
P. M. van den Berg
Keyword(s):  
Integral Equation ◽  
Free Surface ◽  
Synthetic Data ◽  
Computation Time ◽  
Neumann Series ◽  
Data Set ◽  
Seismic Reflection Data ◽  
Prestack Migration ◽  
Reflection Data ◽  
Wave Phenomena

Removal of the effects of the free surface from seismic reflection data is an essential preprocessing step before prestack migration. The problem can be formulated by means of Rayleigh’s reciprocity theorem which leads to an integral equation of the second kind for the desired pressure field that does not include these free‐surface effects. This integral equation can be solved numerically, both in the spatial domain and in the double Radon domain. Solving the integral equation in the double Radon domain has the advantage of reducing the computation time significantly since the kernel of the integral equation becomes dominant diagonally. Two methods are proposed to solve the integral equation: direct matrix inversion and a recursive subtraction of the free‐surface multiples using a Neumann series. Both methods have been developed and tested on a synthetic data set, which was computed with the help of an independent forward‐modeling scheme.

Broadband receiver response from dual‐streamer data and applications in deep reflection seismology

Geophysics ◽  
1996 ◽  
Vol 61 (1) ◽  
pp. 232-243 ◽  
Author(s):  
Satish C. Singh ◽  
R. W. Hobbs ◽  
D. B. Snyder
Keyword(s):  
Signal To Noise Ratio ◽  
Spectral Response ◽  
Synthetic Data ◽  
Real Data ◽  
Reflection Seismology ◽  
Data Set ◽  
Seismic Reflection Data ◽  
Reflection Data ◽  
Streamer Data ◽  
Lower Crustal

A method to process dual‐streamer data with under and over configuration is presented. The method combines the results of dephase‐sum and dephase‐subtraction methods. In the dephase methods, the response of one streamer is time shifted so that the primary arrivals on both streamers are aligned, and these responses are then summed or subtracted. The method provides a broad spectral response from dual‐streamer data and increases the signal‐to‐noise ratio by a factor of 1.5. Testing was done on synthetic data and then applied to a real data set collected by the British Institutions Reflection Profiling Syndicate (BIRPS). Its application to a deep seismic reflection data set from the British Isles shows that the reflections from the lower crust contain frequencies up to 80 Hz, suggesting that some of the lower crustal reflectors may have sharp boundaries and could be 20–30 m thick.

Wavelet-based detection of singularities in acoustic impedances from surface seismic reflection data

Geophysics ◽  
2008 ◽  
Vol 73 (1) ◽  
pp. V1-V9 ◽  
Author(s):  
Chun-Feng Li ◽  
Christopher Liner
Keyword(s):  
Seismic Data ◽  
Acoustic Impedance ◽  
Synthetic Data ◽  
Singularity Analysis ◽  
Real Surface ◽  
Multiple Reflections ◽  
Data Set ◽  
Seismic Reflection Data ◽  
Seismic Processing ◽  
Reflection Data

Although the passage of singularity information from acoustic impedance to seismic traces is now well understood, it remains unanswered how routine seismic processing, mode conversions, and multiple reflections can affect the singularity analysis of surface seismic data. We make theoretical investigations on the transition of singularity behaviors from acoustic impedances to surface seismic data. We also perform numerical, wavelet-based singularity analysis on an elastic synthetic data set that is processed through routine seismic processing steps (such as stacking and migration) and that contains mode conversions, multiple reflections, and other wave-equation effects. Theoretically, seismic traces can be approximated as proportional to a smoothed version of the [Formula: see text] derivative of acoustic impedance,where [Formula: see text] is the vanishing moment of the seismic wavelet. This theoretical approach forms the basis of linking singularity exponents (Hölder exponents) in acoustic impedance with those computable from seismic data. By using wavelet-based multiscale analysis with complex Morlet wavelets, we can estimate singularity strengths and localities in subsurface impedance directly from surface seismic data. Our results indicate that rich singularity information in acoustic impedance variations can be preserved by surface seismic data despite data-acquisition and processing activities. We also show that high-resolution detection of singularities from real surface seismic data can be achieved with a proper choice of the scale of the mother wavelet in the wavelet transform. Singularity detection from surface seismic data thus can play a key role in stratigraphic analysis and acoustic impedance inversion.

Migration‐based traveltime waveform inversion of 2-D simple structures: A synthetic example

Geophysics ◽  
2001 ◽  
Vol 66 (3) ◽  
pp. 845-860 ◽  
Author(s):  
François Clément ◽  
Guy Chavent ◽  
Susana Gómez
Keyword(s):  
Least Squares ◽  
Waveform Inversion ◽  
Domain Of Attraction ◽  
Optimization Methods ◽  
Data Sets ◽  
Seismic Reflection Data ◽  
Prestack Migration ◽  
Reflection Data ◽  
Full Wave ◽  
The Time Domain

Migration‐based traveltime (MBTT) formulation provides algorithms for automatically determining background velocities from full‐waveform surface seismic reflection data using local optimization methods. In particular, it addresses the difficulty of the nonconvexity of the least‐squares data misfit function. The method consists of parameterizing the reflectivity in the time domain through a migration step and providing a multiscale representation for the smooth background velocity. We present an implementation of the MBTT approach for a 2-D finite‐difference (FD) full‐wave acoustic model. Numerical analysis on a 2-D synthetic example shows the ability of the method to find much more reliable estimates of both long and short wavelengths of the velocity than the classical least‐squares approach, even when starting from very poor initial guesses. This enlargement of the domain of attraction for the global minima of the least‐squares misfit has a price: each evaluation of the new objective function requires, besides the usual FD full‐wave forward modeling, an additional full‐wave prestack migration. Hence, the FD implementation of the MBTT approach presented in this paper is expected to provide a useful tool for the inversion of data sets of moderate size.

Shallow seismic reflection study of a glaciated valley

Geophysics ◽  
1998 ◽  
Vol 63 (4) ◽  
pp. 1395-1407 ◽  
Author(s):  
Frank Büker ◽  
Alan G. Green ◽  
Heinrich Horstmeyer
Keyword(s):  
Seismic Reflection ◽  
High Amplitude ◽  
Data Sets ◽  
Seismic Section ◽  
Data Set ◽  
Seismic Reflection Data ◽  
Northern Switzerland ◽  
Reflection Data ◽  
Shallow Seismic ◽  
Glaciated Valley

Shallow seismic reflection data were recorded along two long (>1.6 km) intersecting profiles in the glaciated Suhre Valley of northern Switzerland. Appropriate choice of source and receiver parameters resulted in a high‐fold (36–48) data set with common midpoints every 1.25 m. As for many shallow seismic reflection data sets, upper portions of the shot gathers were contaminated with high‐amplitude, source‐generated noise (e.g., direct, refracted, guided, surface, and airwaves). Spectral balancing was effective in significantly increasing the strength of the reflected signals relative to the source‐generated noise, and application of carefully selected top mutes ensured guided phases were not misprocessed and misinterpreted as reflections. Resultant processed sections were characterized by distributions of distinct seismic reflection patterns or facies that were bounded by quasi‐continuous reflection zones. The uppermost reflection zone at 20 to 50 ms (∼15 to ∼40 m depth) originated from a boundary between glaciolacustrine clays/silts and underlying glacial sands/gravels (till) deposits. Of particular importance was the discovery that the deepest part of the valley floor appeared on the seismic section at traveltimes >180 ms (∼200 m), approximately twice as deep as expected. Constrained by information from boreholes adjacent to the profiles, the various seismic units were interpreted in terms of unconsolidated glacial, glaciofluvial, and glaciolacustrine sediments deposited during two principal phases of glaciation (Riss at >100 000 and Würm at ∼18 000 years before present).

Migration with reduced artifacts from internal multiple reflections

Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. A25-A29
Author(s):  
Lele Zhang
Keyword(s):  
Seismic Reflection ◽  
Computational Cost ◽  
Data Sets ◽  
Square Root ◽  
Multiple Reflections ◽  
Data Set ◽  
Seismic Reflection Data ◽  
Reflection Data ◽  
Recent Developments ◽  
Multiple Elimination

Migration of seismic reflection data leads to artifacts due to the presence of internal multiple reflections. Recent developments have shown that these artifacts can be avoided using Marchenko redatuming or Marchenko multiple elimination. These are powerful concepts, but their implementation comes at a considerable computational cost. We have derived a scheme to image the subsurface of the medium with significantly reduced computational cost and artifacts. This scheme is based on the projected Marchenko equations. The measured reflection response is required as input, and a data set with primary reflections and nonphysical primary reflections is created. Original and retrieved data sets are migrated, and the migration images are multiplied with each other, after which the square root is taken to give the artifact-reduced image. We showed the underlying theory and introduced the effectiveness of this scheme with a 2D numerical example.

Shallow 3D seismic-reflection imaging of fracture zones in crystalline rock

Geophysics ◽  
2007 ◽  
Vol 72 (6) ◽  
pp. B149-B160 ◽  
Author(s):  
Cedric Schmelzbach ◽  
Heinrich Horstmeyer ◽  
Christopher Juhlin
Keyword(s):  
Seismic Reflection ◽  
Crystalline Rock ◽  
Three Dimensions ◽  
Disposal Site ◽  
3D Seismic ◽  
Fracture Zones ◽  
Data Set ◽  
Seismic Reflection Data ◽  
Reflection Data ◽  
Receiver Array

A limited 3D seismic-reflection data set was used to map fracture zones in crystalline rock for a nuclear waste disposal site study. Seismic-reflection data simultaneously recorded along two roughly perpendicular profiles (1850 and [Formula: see text] long) and with a [Formula: see text] receiver array centered at the intersection of the lines sampled a [Formula: see text] area in three dimensions. High levels of source-generated noise required a processing sequence involving surface-consistent deconvolution, which effectively increased the strength of reflected signals, and a linear [Formula: see text] filtering scheme to suppress any remaining direct [Formula: see text]-wave energy. A flexible-binning scheme significantly balanced and increased the CMP fold, but the offset and azimuth distributions remain irregular; a wide azimuth range and offsets [Formula: see text] are concentrated in the center of the survey area although long offsets [Formula: see text] are only found at the edges of the site. Three-dimensional dip moveout and 3D poststack migration were necessary to image events with conflicting dips up to about 40°. Despite the irregular acquisition geometry and the high level of source-generated noise, we obtained images rich in structural detail. Seven continuous to semicontinuous reflection events were traced through the final data volume to a maximum depth of around [Formula: see text]. Previous 2D seismic-reflection studies and borehole data indicate that fracture zones are the most likely cause of the reflections.

Feasibility of approximate broadband estimation of acoustic impedance profile from first principles and band-limited reflection data

Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. R57-R74 ◽  
Author(s):  
Santi Kumar Ghosh ◽  
Animesh Mandal
Keyword(s):  
First Principles ◽  
Acoustic Impedance ◽  
Linear Equations ◽  
Synthetic Data ◽  
Reflection Coefficients ◽  
Limited Data ◽  
Seismic Reflection Data ◽  
Reflection Data ◽  
Impedance Profile ◽  
Band Limited

Because seismic reflection data are band limited, acoustic impedance profiles derived from them are nonunique. The conventional inversion methods counter the nonuniqueness either by stabilizing the answer with respect to an initial model or by imposing mathematical constraints such as sparsity of the reflection coefficients. By making a nominal assumption of an earth model locally consisting of a stack of homogeneous and horizontal layers, we have formulated a set of linear equations in which the reflection coefficients are the unknowns and the recursively integrated seismic trace constitute the data. Drawing only on first principles, the Zoeppritz equation in this case, the approach makes a frontal assault on the problem of reconstructing reflection coefficients from band-limited data. The local layer-cake assumption and the strategy of seeking a singular value decomposition solution of the linear equations counter the nonuniqueness, provided that the objective is to reconstruct a smooth version of the impedance profile that includes only its crude structures. Tests on synthetic data generated from elementary models and from measured logs of acoustic impedance demonstrated the efficacy of the method, even when a significant amount of noise was added to the data. The emergence of consistent estimates of impedance, approximating the original impedance, from synthetic data generated for several frequency bands has inspired our confidence in the method. The other attractive outputs of the method are as follows: (1) an accurate estimate of the impedance mean, (2) an accurate reconstruction of the direct-current (DC) frequency of the reflectivity, and (3) an acceptable reconstruction of the broad outline of the original impedance profile. These outputs can serve as constraints for either more refined inversions or geologic interpretations. Beginning from the restriction of band-limited data, we have devised a method that neither requires a starting input model nor imposes mathematical constraints on the earth reflectivity and still yielded significant and relevant geologic information.

Minimizing field operations in shallow 3‐D seismic reflection surveying

Geophysics ◽  
2001 ◽  
Vol 66 (6) ◽  
pp. 1761-1773 ◽  
Author(s):  
Roman Spitzer ◽  
Alan G. Green ◽  
Frank O. Nitsche
Keyword(s):  
Seismic Reflection ◽  
Cost Effective ◽  
Unconsolidated Sediments ◽  
Data Set ◽  
Seismic Reflection Data ◽  
Northern Switzerland ◽  
Shallow Subsurface ◽  
Reflection Data ◽  
Effective Manner ◽  
Original Survey

By appropriately decimating a comprehensive shallow 3‐D seismic reflection data set recorded across unconsolidated sediments in northern Switzerland, we have investigated the potential and limitations of four different source‐receiver acquisition patterns. For the original survey, more than 12 000 shots and 18 000 receivers deployed on a [Formula: see text] grid resulted in common midpoint (CMP) data with an average fold of ∼40 across a [Formula: see text] area. A principal goal of our investigation was to determine an acquisition strategy capable of producing reliable subsurface images in a more efficient and cost‐effective manner. Field efforts for the four tested acquisition strategies were approximately 50%, 50%, 25%, and 20% of the original effort. All four data subsets were subjected to a common processing sequence. Static corrections, top‐mute functions, and stacking velocities were estimated individually for each subset. Because shallow reflections were difficult to discern on shot and CMP gathers generated with the lowest density acquisition pattern (20% field effort) such that dependable top‐mute functions could not be estimated, data resulting from this acquisition pattern were not processed to completion. Of the three fully processed data subsets, two (50% field effort and 25% field effort) yielded 3‐D migrated images comparable to that derived from the entire data set, whereas the third (50% field effort) resulted in good‐quality images only in the shallow subsurface because of a lack of far‐offset data. On the basis of these results, we concluded that all geological objectives associated with our particular study site, which included mapping complex lithological units and their intervening shallow dipping boundaries, would have been achieved by conducting a 3‐D seismic reflection survey that was 75% less expensive than the original one.

scholarly journals Reconstructing the primary reflections in seismic data by Marchenko redatuming and convolutional interferometry

Geophysics ◽  
2016 ◽  
Vol 81 (2) ◽  
pp. Q15-Q26 ◽  
Author(s):  
Giovanni Angelo Meles ◽  
Kees Wapenaar ◽  
Andrew Curtis
Keyword(s):  
Velocity Model ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Seismic Reflection Data ◽  
Surface Reflection ◽  
Reflection Data ◽  
Reference Velocity ◽  
Time Migration ◽  
Recorded Data

State-of-the-art methods to image the earth’s subsurface using active-source seismic reflection data involve reverse time migration. This and other standard seismic processing methods such as velocity analysis provide best results only when all waves in the data set are primaries (waves reflected only once). A variety of methods are therefore deployed as processing to predict and remove multiples (waves reflected several times); however, accurate removal of those predicted multiples from the recorded data using adaptive subtraction techniques proves challenging, even in cases in which they can be predicted with reasonable accuracy. We present a new, alternative strategy to construct a parallel data set consisting only of primaries, which is calculated directly from recorded data. This obviates the need for multiple prediction and removal methods. Primaries are constructed by using convolutional interferometry to combine the first-arriving events of upgoing and direct-wave downgoing Green’s functions to virtual receivers in the subsurface. The required upgoing wavefields to virtual receivers are constructed by Marchenko redatuming. Crucially, this is possible without detailed models of the earth’s subsurface reflectivity structure: Similar to the most migration techniques, the method only requires surface reflection data and estimates of direct (nonreflected) arrivals between the virtual subsurface sources and the acquisition surface. We evaluate the method on a stratified synclinal model. It is shown to be particularly robust against errors in the reference velocity model used and to improve the migrated images substantially.

Sign in / Sign up

Export Citation Format

Share Document