Multiple water reflections recorded at the ocean bottom

1978 ◽  
Vol 15 (1) ◽  
pp. 78-85 ◽  
Author(s):  
George A. McMechan

The use of direct arrivals and multiple reflections that have travelled completely in water from source to receiver to determine epicentral distances is a standard technique in the analysis of marine seismic profiles. The configuration of a source at the air–water interface and a seismometer at the water–sediment interface is investigated in the ray parameter – distance plane and the travel time – distance plane. Vertical component synthetic seismograms are computed by the Cagniard – de Hoop algorithm and are compared with seismograms recorded at the ocean bottom. The results explain the prominent features of the observed wavetrains, including the asymptotic behaviour of arrivals, the location of caustics and the variable observability of arrivals as a fu nction of distance.

1976 ◽  
Vol 13 (11) ◽  
pp. 1545-1555 ◽  
Author(s):  
R. M. Clowes ◽  
S. J. Malecek

A marine seismic system for recording near-vertical incidence to wide-angle reflected waves and refracted waves with penetration from the ocean bottom to the upper mantle (deep seismic sounding or DSS) has been developed. Signals from six individual hydrophones suspended at 45 m depth from a 600 m cable trailed behind the receiving ship are recorded in digital form. Using charges ranging from 2.3 to 280 kg, two reversed DSS profiles were recorded in the region of Explorer ridge during 1974. A preliminary interpretation of the profiles based on first-arrival information in the range 4 to 80 km has been made.The reversed profile run across the ridge showed no anomalous effects as the ridge was crossed; the profile on Juan de Fuca plate, paralleling the ridge, exhibited traveltime branch offsets and delays. These have been interpreted as due to faulting with a vertical component of offset of about 4 km. The reversed upper mantle velocities are 7.85 and'7.30 km/s indirections perpendicular and parallel to the ridge. Anisotropy is proposed to explain these different velocities and gives a 7% anisotropic effect. The data require that 'layer 2' comprise at least two layers with velocities of 4.13 km/s and 5.25 km/s and individual depth extents ranging from 1 to 2 km. Compared with crustal sections from other ridge areas, the interpretation gives a thick 'layer 3' (up to 6 km) near the ridge crest. The sub-bottom thickness of the oceanic crust varies between 7 and 9 km, except in the faulted region, where the 7.30 km/s material is present less than 3 km from the bottom.


Geophysics ◽  
2001 ◽  
Vol 66 (1) ◽  
pp. 327-341 ◽  
Author(s):  
Lasse Amundsen

This paper presents a new, wave‐equation based method for eliminating the effect of the free surface from marine seismic data without destroying primary amplitudes and without any knowledge of the subsurface. Compared with previously published methods which require an estimate of the source wavelet, the present method has the following characteristics: it does not require any information about the marine source array and its signature, it does not rely on removal of the direct wave from the data, and it does not require any explicit deghosting. Moreover, the effect of the source signature is removed from the data in the multiple elimination process by deterministic signature deconvolution, replacing the original source signature radiated from the marine source array with any desired wavelet (within the data frequency‐band) radiated from a monopole point source. The fundamental constraint of the new method is that the vertical derivative of the pressure or the vertical component of the particle velocity is input to the free‐surface demultiple process along with pressure recordings. These additional data are routinely recorded in ocean‐bottom seismic surveys. The method can be applied to conventional towed streamer pressure data recorded in the water column at a depth which is greater than the depth of the source array only when the pressure derivative can be estimated, or even better, is measured. Since the direct wave and its source ghost is part of the free‐ surface demultiple, designature process, the direct arrival must be properly measured for the method to work successfully. In the case when the geology is close to horizontally layering, the free‐surface multiple elimination method greatly simplifies, reducing to a well‐known deterministic deconvolution process which can be applied to common shot gathers (or common receiver gathers or common midpoint gathers when source array variations are negligible) in the τ-p domain or frequency‐wavenumber domain.


2021 ◽  
Author(s):  
Pimpawee Sittipan ◽  
Pisanu Wongpornchai

Some of the important petroleum reservoirs accumulate beneath the seas and oceans. Marine seismic reflection method is the most efficient method and is widely used in the petroleum industry to map and interpret the potential of petroleum reservoirs. Multiple reflections are a particular problem in marine seismic reflection investigation, as they often obscure the target reflectors in seismic profiles. Multiple reflections can be categorized by considering the shallowest interface on which the bounces take place into two types: internal multiples and surface-related multiples. Besides, the multiples can be categorized on the interfaces where the bounces take place, a difference between long-period and short-period multiples can be considered. The long-period surface-related multiples on 2D marine seismic data of the East Coast of the United States-Southern Atlantic Margin were focused on this research. The seismic profile demonstrates the effectiveness of the results from predictive deconvolution and the combination of surface-related multiple eliminations (SRME) and parabolic Radon filtering. First, predictive deconvolution applied on conventional processing is the method of multiple suppression. The other, SRME is a model-based and data-driven surface-related multiple elimination method which does not need any assumptions. And the last, parabolic Radon filtering is a moveout-based method for residual multiple reflections based on velocity discrimination between primary and multiple reflections, thus velocity model and normal-moveout correction are required for this method. The predictive deconvolution is ineffective for long-period surface-related multiple removals. However, the combination of SRME and parabolic Radon filtering can attenuate almost long-period surface-related multiple reflections and provide a high-quality seismic images of marine seismic data.


2016 ◽  
Vol 187 (4-5) ◽  
pp. 187-215 ◽  
Author(s):  
François Fournier ◽  
Aurélie Tassy ◽  
Isabelle Thinon ◽  
Philippe Münch ◽  
Jean-Jacques Cornée ◽  
...  

AbstractThe seaward extension of onshore formations and structures were previously almost unknown in Provence. The interpretation of 2D high-resolution marine seismic profiles together with the integration of sea-bottom rock samples provides new insights into the stratigraphic, structural and paleogeographic framework of pre-Messinian Salinity Crisis (MSC) deposits of the Provence continental shelf. Seven post-Jurassic seismic units have been identified on seismic profiles, mapped throughout the offshore Provence area and correlated with the onshore series. The studied marine surface and sub-surface database provided new insights into the mid and late Cretaceous paleogeography and structural framework as well as into the syn- and post-rift deformation in Provence. Thick (up to 2000 m) Aptian-Albian series whose deposition is controlled by E-W-trending faults are evidenced offshore. The occurrence and location of the Upper Cretaceous South-Provence basin is confirmed by the thick (up to 1500 m) basinal series downlaping the Aptian-Albian unit. This basin was fed in terrigenous sediments by a southern massif (“Massif Méridional”) whose present-day relict is the Paleozoic basement and its sedimentary cover from the Sicié imbricate. In the bay of Marseille, thick syn-rift (Rupelian to Aquitanian) deposition occurred (>1000 m). During the rifting phase, syn-sedimentary deformations consist of dominant N040 to N060 sub-vertical faults with a normal component and N050 drag-synclines and anticlines. The syn-rift and early post-rift units (Rupelian to early Burdigalian) are deformed and form a set of E-W-trending en echelon folds that may result from sinistral strike-slip reactivation of N040 to N060 normal faults during a N-S compressive phase of early-to-mid Burdigalian age (18–20 Ma). Finally, minor fault reactivation and local folding affect post-rift deposits within a N160-trending corridor localized south of La Couronne, and could result from a later, post-Burdigalian and pre-Pliocene compressive phase.


Geophysics ◽  
1988 ◽  
Vol 53 (4) ◽  
pp. 453-465 ◽  
Author(s):  
Michel Dietrich

In many cases, comparison of real data with synthetic seismograms provides additional constraints on the velocity‐depth profiles obtained with simple inversion techniques. Obtaining a satisfactory match between the real and computed data usually requires several trials with different models but can be performed rapidly if the theoretical seismograms are themselves easily interpretable, i.e., if the major contributions which make up the synthetic traces can be identified and separated. In horizontally stratified media, modeling is further simplified and is faster if the simulation techniques are implemented in the domain of the intercept time τ and the ray parameter p. The generalized reflection and transmission matrix method is well suited for these purposes and may be used to generate synthetic seismograms in both the t-x and τ-p planes. Computation of plane‐wave seismograms is straightforward and merely corresponds to an inverse Fourier transform of the overall reflectivity matrix for each ray parameter. The construction of point‐source seismograms can be carried out in several ways. In this paper, I combine the generalized reflection and transmission matrix method with a discrete wavenumber integration of the reflectivity function, extending previous work to include fluid‐solid interfaces. The introduction of scaling parameters also simplifies the reflectivity matrices. Simple numerical experiments demonstrate that the relations between the earth model and the corresponding seismic response are simpler in the τ-p domain than in the t-x domain. In particular, calculation of the PP, SP, and SS contributions to the complete seismic response shows that shear and converted waves may have a clearer expression in the τ-p plane than in the t-x plane and can in some cases provide discrimination between several earth models. Finally, the main features of a real marine seismic profile are well reproduced by synthetic sections in both the t-x and τ-p domains.


Geophysics ◽  
1989 ◽  
Vol 54 (1) ◽  
pp. 122-126 ◽  
Author(s):  
R. J. J. Hardy ◽  
M. R. Warner ◽  
R. W. Hobbs

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.


Geophysics ◽  
1956 ◽  
Vol 21 (2) ◽  
pp. 320-336 ◽  
Author(s):  
George P. Sarrafian

A model for the study of marine seismic phenomena is described. Study of multiple‐reflection phenomena forms the basis for the course of experiments. It is shown that the multiple‐reflection phenomenon of a disturbance with slowly decaying amplitude may be duplicated in the model. Multiple‐reflection problems are studied in which the bottom of the water layer is tilted or thin. A mass of air bubbles is shown to be of use in attenuating multiple reflections. The possible application of the marine model in a search for information about certain problems in field prospecting is suggested.


Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. WB87-WB94 ◽  
Author(s):  
Lasse Amundsen ◽  
Harald Westerdahl ◽  
Mark Thompson ◽  
Jon Andre Haugen ◽  
Arne Reitan ◽  
...  

In ocean-bottom seismic and vertical-cable surveying, receiver stations are stationary on the sea floor while a source vessel shoots on a predetermined [Formula: see text] grid on the sea surface. To reduce exploration cost, the shot point interval often is so coarse that the data recorded at a given receiver station are undersampled and thus irrecoverably aliased. However, when the pressure field and its [Formula: see text]- and [Formula: see text]-derivatives are measured in the water column, the nonaliased pressure field can be reconstructed by interpolation. Likewise, if the vertical component of the particle velocity (or acceleration) and its [Formula: see text]- and [Formula: see text]-derivatives are measured, then this component can also be reconstructed by interpolation. The interpolation scheme can be any scheme that reconstructs the field from its sampled values and sampled derivatives. In the case that the two fields’ first-order derivatives are recorded, the total number of components is six. When also their second-order derivatives are measured, the number of components is 10. The properly interpolated measurements of pressure and vertical component of particle velocity from the multicomponent measurements allow true 3D up/down wavefield decomposition (deghosting) and wave-equation demultiple before wave-equation migration.


Sign in / Sign up

Export Citation Format

Share Document