Efficient acquisition, processing, and interpretation strategy for shallow 3D seismic surveying: A Case Study

Geophysics ◽  
2003 ◽  
Vol 68 (6) ◽  
pp. 1792-1806 ◽  
Author(s):  
Roman Spitzer ◽  
Frank O. Nitsche ◽  
Alan G. Green ◽  
Heinrich Horstmeyer

A new 3D seismic reflection data set has been used to map the shallow subsurface beneath a key region of the Swiss Rhine Valley. Seismic signals generated by a pipegun were recorded with single 30‐Hz geophones distributed across a 277.5 × 357.0‐m area. The dense distribution of sources and receivers resulted in a binning grid of 2.12 × 2.12 m and an average fold of ∼22. To improve the visibility and continuity of reflections, a novel processing strategy was designed and applied to the acquired data. A combination of regridding and sharing traces in the common midpoint (CMP) domain resulted in increased S/N ratios with only minor loss of resolution. This prestack interpolation method yielded composite CMPs distributed on a 1.5 × 1.5‐m binning grid and an increased average fold of ∼44. The composite CMPs were subjected to a combined linear and hyperbolic τ–p processing scheme that led to the effective separation of reflections from source‐generated noise. Finally, 3D depth migration of the stacked data produced high‐resolution images of the subsurface from ∼15 to ∼130 m depth. On the basis of characteristic seismic facies and information from nearby boreholes, four principal lithological units were identified. At increasing depths they were glaciofluvial sand and gravel, glaciolacustrine clay and silt, morainal deposits, and sandstone basement. These lithological units were separated by three principal reflecting boundaries that were mapped through the data volume using semiautomatic tracking procedures. The deepest boundary defined a trough‐shaped basement structure.

Geophysics ◽  
2007 ◽  
Vol 72 (6) ◽  
pp. B149-B160 ◽  
Author(s):  
Cedric Schmelzbach ◽  
Heinrich Horstmeyer ◽  
Christopher Juhlin

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.


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


2014 ◽  
Vol 2 (1) ◽  
pp. SA151-SA162 ◽  
Author(s):  
John H. McBride ◽  
R. William Keach ◽  
Eugene E. Wolfe ◽  
Hannes E. Leetaru ◽  
Clayton K. Chandler ◽  
...  

Because the confinement of [Formula: see text] in a storage reservoir depends on a stratigraphically continuous set of seals to isolate the fluid in the reservoir, the detection of structural anomalies is critical for guiding any assessment of a potential subsurface carbon storage site. Employing a suite of 3D seismic attribute analyses (as opposed to relying upon a single attribute) maximizes the chances of identifying geologic anomalies or discontinuities (e.g., faults) that may affect the integrity of a seal that will confine the stored [Formula: see text] in the reservoir. The Illinois Basin, a major area for potential carbon storage, presents challenges for target assessment because geologic anomalies can be ambiguous and easily misinterpreted when using 2D seismic reflection data, or even 3D data, if only conventional display techniques are used. We procured a small 3D seismic reflection data set in the central part of the basin (Stewardson oil field) to experiment with different strategies for enhancing the appearance of discontinuities by integrating 3D seismic attribute analyses with conventional visualizations. Focusing on zones above and below the target interval of the Cambrian Mt. Simon Sandstone, we computed attribute traveltime slices (combined with vertical views) based on discontinuity computations, crossline-directed amplitude change, azimuth of the dip, shaded relief, and fault likelihood attributes. The results provided instructive examples of how discontinuities (e.g., subseismic scale faults) may be almost “invisible” on conventional displays but become detectable and mappable using an appropriate integration of 3D attributes. Strong discontinuities in underlying Precambrian basement rocks do not necessarily propagate upward into the target carbon storage interval. The origin of these discontinuities is uncertain, but we explored a possible strike-slip role that also explains the localization of a structural embayment developed in Lower Paleozoic strata above the basement discontinuities.


Sign in / Sign up

Export Citation Format

Share Document