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

Labeling long‐period multiple reflections

Geophysics ◽  
1989 ◽  
Vol 54 (1) ◽  
pp. 122-126 ◽  
Author(s):  
R. J. J. Hardy ◽  
M. R. Warner ◽  
R. W. Hobbs
Keyword(s):  
Seismic Reflection ◽  
High Amplitude ◽  
Data Sets ◽  
Multiple Reflections ◽  
Multiple Series ◽  
Seismic Reflection Data ◽  
Long Period ◽  
Reflection Data ◽  
Zero Offset ◽  
Marine Seismic

The many techniques that have been developed to remove multiple reflections from seismic data all leave remnant energy which can cause ambiguity in interpretation. The removal methods are mostly based on periodicity (e.g., Sinton et al., 1978) or the moveout difference between primary and multiple events (e.g., Schneider et al., 1965). They work on synthetic and selected field data sets but are rather unsatisfactory when applied to high‐amplitude, long‐period multiples in marine seismic reflection data acquired in moderately deep (700 m to 3 km) water. Differential moveout is often better than periodicity at discriminating between types of events because, while a multiple series may look periodic to the eye, it is only exactly so on zero‐offset reflections from horizontal layers. The technique of seismic event labeling described below works by returning offset information from CDP gathers to a stacked section by color coding, thereby discriminating between seismic reflection events by differential normal moveout. Events appear as a superposition of colors; the direction of color fringes indicates whether an event has been overcorrected or undercorrected for its hyperbolic normal moveout.

EXPLORATION HORIZONS FROM NEW SEISMIC CONCEPTS OF CDP AND DIGITAL PROCESSING

Geophysics ◽  
1967 ◽  
Vol 32 (2) ◽  
pp. 207-224 ◽  
Author(s):  
John D. Marr ◽  
Edward F. Zagst
Keyword(s):  
Data Processing ◽  
Seismic Data ◽  
Seismic Reflection ◽  
Synthetic Data ◽  
Digital Processing ◽  
Seismic Data Processing ◽  
Multiple Reflections ◽  
Seismic Reflection Data ◽  
Reflection Data ◽  
Recent Developments

The more recent developments in common‐depth‐point techniques to attenuate multiple reflections have resulted in an exploration capability comparable to the development of the seismic reflection method. The combination of new concepts in digital seismic data processing with CDP techniques is creating unforeseen exploration horizons with vastly improved seismic data. Major improvements in multiple reflection and reverberation attenuation are now attainable with appropriate CDP geometry and special CDP stacking procedures. Further major improvements are clearly evident in the very near future with the use of multichannel digital filtering‐stacking techniques and the application of deconvolution as the first step in seismic data processing. CDP techniques are briefly reviewed and evaluated with real and experimental data. Synthetic data are used to illustrate that all seismic reflection data should be deconvolved as the first processing step.

Ground‐roll suppression from deep crustal seismic reflection data using a wavelet‐based approach: A case study from western Canada

Geophysics ◽  
2004 ◽  
Vol 69 (4) ◽  
pp. 877-884 ◽  
Author(s):  
J. Kim Welford ◽  
Rongfeng Zhang
Keyword(s):  
Seismic Reflection ◽  
Low Frequency ◽  
Frequency Separation ◽  
Frequency Noise ◽  
Data Sets ◽  
Data Set ◽  
Dispersive Waves ◽  
Seismic Reflection Data ◽  
Reflection Data ◽  
Ground Roll

Seismic reflection data sets recorded on land are often contaminated by coherent ground‐roll noise generated by the propagation of dispersive waves along the free surface. For crustal‐scale investigations, this ground‐roll contamination can be particularly harmful as the higher amplitude, low‐frequency noise overwhelms low‐frequency signals coming from deep reflectors. Consequently, conventional ground‐roll suppression techniques which rely on frequency separation of ground roll from signal become ineffective for crustal studies. This paper presents the successful use of a new 2D wavelet method based on frame theory (physical wavelet frame denoising) in removing ground roll from a deep 3D reflection data set intended for the study of upper crustal Precambrian mafic sills in southwestern Alberta, Canada.

scholarly journals Improved Interpretation of Deep Seismic Reflection Data in Areas of Complex Geology Through Integration With Passive Seismic Data Sets

2018 ◽  
Vol 123 (12) ◽  
pp. 10,810-10,830
Author(s):  
Michael Dentith ◽  
Huaiyu Yuan ◽  
Ruth Elaine Murdie ◽  
Perla Pina-Varas ◽  
Simon P. Johnson ◽  
...  
Keyword(s):  
Seismic Data ◽  
Seismic Reflection ◽  
Data Sets ◽  
Seismic Reflection Data ◽  
Deep Seismic Reflection ◽  
Reflection Data ◽  
Passive Seismic ◽  
Complex Geology

scholarly journals Data-driven internal multiple elimination and its consequences for imaging: A comparison of strategies

Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. S365-S372 ◽  
Author(s):  
Lele Zhang ◽  
Jan Thorbecke ◽  
Kees Wapenaar ◽  
Evert Slob
Keyword(s):  
Inverse Scattering ◽  
Computational Cost ◽  
Multiple Reflection ◽  
Data Driven ◽  
Transmission Losses ◽  
Multiple Reflections ◽  
Data Set ◽  
Numerical Examples ◽  
Inverse Scattering Series ◽  
Multiple Elimination

We have compared three data-driven internal multiple reflection elimination schemes derived from the Marchenko equations and inverse scattering series (ISS). The two schemes derived from Marchenko equations are similar but use different truncation operators. The first scheme creates a new data set without internal multiple reflections. The second scheme does the same and compensates for transmission losses in the primary reflections. The scheme derived from ISS is equal to the result after the first iteration of the first Marchenko-based scheme. It can attenuate internal multiple reflections with residuals. We evaluate the success of these schemes with 2D numerical examples. It is shown that Marchenko-based data-driven schemes are relatively more robust for internal multiple reflection elimination at a higher computational cost.

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.

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 Target-enclosed seismic imaging

Geophysics ◽  
2017 ◽  
Vol 82 (6) ◽  
pp. Q53-Q66 ◽  
Author(s):  
Joost van der Neut ◽  
Matteo Ravasi ◽  
Yi Liu ◽  
Ivan Vasconcelos
Keyword(s):  
Seismic Reflection ◽  
Target Volume ◽  
Multiple Reflections ◽  
Reflection And Transmission ◽  
Seismic Reflection Data ◽  
Computational Burden ◽  
Reflection Data ◽  
Seismic Resolution ◽  
Marchenko Equation ◽  
Deconvolution Problem

Seismic reflection data can be redatumed to a specified boundary in the subsurface by solving an inverse (or multidimensional deconvolution) problem. The redatumed data can be interpreted as an extended image of the subsurface at the redatuming boundary, depending on the subsurface offset and time. We retrieve target-enclosed extended images by using two redatuming boundaries, which are selected above and below a specified target volume. As input, we require the upgoing and downgoing wavefields at both redatuming boundaries due to impulsive sources at the earth’s surface. These wavefields can be obtained from actual measurements in the subsurface, they can be numerically modeled, or they can be retrieved by solving a multidimensional Marchenko equation. As output, we retrieved virtual reflection and transmission responses as if sources and receivers were located at the two target-enclosing boundaries. These data contain all orders of reflections inside the target volume but exclude all interactions with the part of the medium outside this volume. The retrieved reflection responses can be used to image the target volume from above or from below. We found that the images from above and below are similar (given that the Marchenko equation is used for wavefield retrieval). If a model with sharp boundaries in the target volume is available, the redatumed data can also be used for two-sided imaging, where the retrieved reflection and transmission responses are exploited. Because multiple reflections are used by this strategy, seismic resolution can be improved significantly. Because target-enclosed extended images are independent on the part of the medium outside the target volume, our methodology is also beneficial to reduce the computational burden of localized inversion, which can now be applied inside the target volume only, without suffering from interactions with other parts of the medium.

Sign in / Sign up

Export Citation Format

Share Document