Reducing source‐generated noise in shallow seismic data using linear and hyperbolic τ‐p transformations

Geophysics ◽  
2001 ◽  
Vol 66 (5) ◽  
pp. 1612-1621 ◽  
Author(s):  
Roman Spitzer ◽  
Frank O. Nitsche ◽  
Alan G. Green

High‐resolution seismic reflection data recorded at many locations on the earth are plagued by the overwhelming effects of direct, refracted, guided, and surface waves. These different components of source‐generated noise may completely mask reflections at traveltimes <∼50–100 ms. Conventional processing methods that include the time‐consuming application of mute functions may lead to the misprocessing of source‐generated noise (especially guided waves) as reflected events and/or the unintentional removal of important shallow reflections. We introduce a combined linear and hyperbolic τ‐p processing scheme that results in the effective separation of reflections from source‐generated noise. After applying linear moveout terms that adjust the direct, refracted, and guided arrivals to appear horizontal to subhorizontal, the reduced traveltime shot gathers are transformed into the linear τ‐p domain. It is then straightforward to design a single τ‐p filter that eliminates most of the source‐generated noise throughout the entire data set. Following inverse linear τ‐p transformation and removal of the linear moveout terms, the filtered shot gathers contain reflections and residual elements of the source‐generated noise. Because summing along hyperbolas favors reflections, transforming the filtered shot gathers into the hyperbolic τ‐p domain leads to significant enhancements in the S/N ratio. A simple rescaling of data values in the hyperbolic τ‐p domain, which results in the loss of true amplitude information, increases further the relative strength of the reflected signals. Finally, inverse hyperbolic transformation yields shot gathers dominated by reflections. In tests of the combined τ‐p processing scheme on a synthetic shot gather and on a complete shallow seismic reflection data set recorded in northern Switzerland, significant improvements in the quality of reflections in the prestacked data and on a fully processed section are readily apparent. According to the results of these tests, the new scheme works well for reflections originating from flat and dipping horizons.

Geophysics ◽  
1998 ◽  
Vol 63 (4) ◽  
pp. 1395-1407 ◽  
Author(s):  
Frank Büker ◽  
Alan G. Green ◽  
Heinrich Horstmeyer

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


Geophysics ◽  
2001 ◽  
Vol 66 (6) ◽  
pp. 1761-1773 ◽  
Author(s):  
Roman Spitzer ◽  
Alan G. Green ◽  
Frank O. Nitsche

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.


Geophysics ◽  
2012 ◽  
Vol 77 (5) ◽  
pp. B227-B236 ◽  
Author(s):  
Hassan Masoomzadeh ◽  
Satish C. Singh ◽  
Penny J. Barton

We developed a method of moveout correction in the [Formula: see text] domain to tackle some of the problems associated with processing wide-angle seismic reflection data, including residual moveout and normal-moveout stretching. We evaluated the concept of the shifted ellipse in the [Formula: see text] domain as an alternative to the well-known concept of the shifted hyperbola in the [Formula: see text] domain. We used this shifted-ellipse concept to address the problem of residual moveout caused by vertical heterogeneity in the subsurface. We also addressed the stretching problem associated with dynamic corrections by combining selected strips from a set of constant-moveout stacks generated using a shifted-ellipse equation. Application of this method to a wide-angle data set from the Faeroe-Shetland Basin provided an enhanced image of the subbasalt structure.


Geophysics ◽  
1998 ◽  
Vol 63 (4) ◽  
pp. 1434-1450 ◽  
Author(s):  
Frank Büker ◽  
Alan G. Green ◽  
Heinrich Horstmeyer

A comprehensive strategy of 3-D seismic reflection data acquisition and processing has been used in a study of glacial sediments deposited within a Swiss mountain valley. Seismic data generated by a downhole shotgun source were recorded with single 30-Hz geophones distributed at 3 m × 3 m intervals across a 357 m × 432 m area. For most common‐midpoint (CMP) bins, traces covering a full range of azimuths and source‐receiver distances of ∼2 to ∼125 m were recorded. A common processing scheme was applied to the entire data set and to various subsets designed to simulate data volumes collected with lower density source and receiver patterns. Comparisons of seismic sections extracted from the processed 3-D subsets demonstrated that high‐fold (>40) and densely spaced (CMP bin sizes ⩽ 3 m × 3 m) data with relatively large numbers (>6) of traces recorded at short (<20 m) source‐receiver offsets were essential for obtaining clear images of the shallowest (<100 ms) reflecting horizons. Reflections rich in frequencies >100 Hz at traveltimes of ∼20 to ∼170 ms provided a vertical resolution of 3 to 6 m over a depth range of ∼15 to ∼150 m. The shallowest prominent reflection at 20 to 35 ms (∼15 to 27 m depth) originated from the boundary between a near‐surface sequence of clays/silts and an underlying unit of heterogeneous sands/gravels.


Geophysics ◽  
1994 ◽  
Vol 59 (3) ◽  
pp. 402-410 ◽  
Author(s):  
Ross A. Black ◽  
Don W. Steeples ◽  
Richard D. Miller

We present an analysis of migration effects on seismic reflection images of very shallow targets such as those that are common objectives of engineering, groundwater, and environmental investigations. We use an example of seismic reflection data from depths of 5 to 15 m that show negligible effect from migration, despite the apparent steep dip on the seismic section. Our analysis of the question of when to migrate shallow reflection data indicates it is critical to take into account the highly variable near‐surface velocities and the vertical exaggeration on the seismic section. A simple set of calculations is developed as well as a flow chart based on the “migrator’s equation” that can predict whether migration of an arbitrary shallow seismic section is advisable. Because shallow reflection data are often processed on personal computers, unnecessary migration of a large data set can be prohibitively time‐consuming and wasteful.


Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. A25-A29
Author(s):  
Lele Zhang

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.


Sign in / Sign up

Export Citation Format

Share Document