Elimination of the overburden response from multicomponent source and receiver seismic data, with source designature and decomposition into PP-, PS-, SP-, and SS-wave responses

Geophysics ◽  
2005 ◽  
Vol 70 (2) ◽  
pp. S43-S59 ◽  
Author(s):  
Egil Holvik ◽  
Lasse Amundsen

This paper shows that Betti's reciprocity theorem gives an integral equation procedure to eliminate from the physical multicomponent-source, multicomponent-receiver seismic measurements the effect of the physical source radiation pattern and the response of the physical overburden (that is, the medium above the receiver plane). The physical multicomponent sources are assumed to be orthogonally aligned anywhere above the multicomponent-receiver depth level. Other than the position of the sources, no source characteristics are required. The method, denoted the Betti designature/elastic demultiple, has the following additional characteristics: it preserves primary amplitudes while eliminating all waves scattered from the overburden; it requires no knowledge of the medium below the receiver level; it requires no knowledge of the medium above the receiver level; it requires information only of the local density and elastic wave propagation velocities at the receiver level to decompose the physical seismic measurements into upgoing and downgoing waves. Following the Betti designature/elastic demultiple step is an elastic wavefield decomposition step that decomposes the data into PP-, PS-, SP-, and SS-wave responses that would be recorded from idealized compressional-wave and shear-wave sources and receivers. The combined elastic wavefield decomposition on the source and receiver side gives data equivalent to data from a hypothetical survey with overburden absent, with single-component compressional and shear-wave sources, and single-component compressional and shear-wave receivers. When the medium is horizontally layered, the Betti designature/elastic demultiple scheme followed by the elastic source-receiver decomposition scheme greatly simplifies and is conveniently implemented as deterministic multidimensional deconvolution and elastic source-receiver wavefield decomposition of common-source gathers (or common-receiver gathers when source array variations are negligible). Betti designature/elastic demultiple followed by source-receiver wavefield decomposition applies to three different seismic experiments: a 9-component (9C) land seismic experiment, a 12-component (12C) ocean-bottom seismic experiment, and an 18-component (18C) borehole seismic experiment. For the land and ocean-bottom seismic experiments, an additional geophone should be deployed below the zero-offset geophone to predict the source-induced vertical traction vector at the source location. A numerical example for the 12C ocean-bottom seismic experiment over a horizontally layered medium validates the Betti designature/elastic demultiple scheme.

Geophysics ◽  
2003 ◽  
Vol 68 (3) ◽  
pp. 1091-1102 ◽  
Author(s):  
K. M. Schalkwijk ◽  
C. P. A. Wapenaar ◽  
D. J. Verschuur

With wavefield decomposition, the recorded wavefield at a certain depth level can be separated into upgoing and downgoing wavefields as well as into P‐ and S‐waves. The medium parameters at the considered depth level (e.g., just below the ocean‐bottom) need to be known in order to be able to do a decomposition. In general, these parameters are unknown and, in addition, measurement‐related issues, such as geophone coupling and crosstalk between the different components, need to be dealt with. In order to apply decomposition to field data, an adaptive five‐stage decomposition scheme was developed in which these issues are addressed. In this study, the adaptive decomposition scheme is tested on a data example with a relatively shallow water depth (∼120 m), consisting recordings from of a full line of ocean‐bottom receivers. Although some of the individual stages in the decomposition scheme are more difficult to apply because of stronger interference between events compared to data acquired over deeper water, the end result is satisfying. Also, a good decomposition result is obtained for the S‐waves. The extension of the decomposition scheme to a complete line of ocean‐bottom cable data consists of a repeated application of the procedure for each receiver. The resulting decomposed upgoing P‐ and S‐wavefields are processed, yielding poststack time migrated images of the subsurface. Comparison with the images obtained from the original (i.e., not decomposed) measurements shows that wavefield decomposition just below the ocean bottom leads to a strong attenuation of multiply reflected events at the sea surface and better event definition in both P‐ and S‐wave sections. Other decomposition effects like improved angle‐dependent amplitudes cannot be evaluated in this way.


Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. V281-V293 ◽  
Author(s):  
Qiang Zhao ◽  
Qizhen Du ◽  
Xufei Gong ◽  
Xiangyang Li ◽  
Liyun Fu ◽  
...  

Simultaneous source acquisition has attracted more and more attention from geophysicists because of its cost savings, whereas it also brings some challenges that have never been addressed before. Deblending of simultaneous source data is usually considered as an underdetermined inverse problem, which can be effectively solved with a least-squares (LS) iterative procedure between data consistency ([Formula: see text]-norm) and regularization ([Formula: see text]-norm or [Formula: see text]-norm). However, when it comes to abnormal noise that follows non-Gaussian distribution and possesses high-amplitude features (e.g., erratic noise, swell noise, and power line noise), the [Formula: see text]-norm is a nonrobust statistic that can easily lead to suboptimal deblended results. Although abnormal noise can be attenuated in the common source domain at first, it is still challenging to apply a coherency-based filter due to the sparse receiver or crossline sampling, e.g., that commonly found in ocean bottom node (OBN) acquisition. To address this problem, we have developed a normalized shaping regularization to make the inversion-based deblending approach robust for the separation of blended data when abnormal noise exists. Its robustness comes from the normalized shaping operator defined by the confidence interval of normal distribution, which minimizes the abnormal risk to a normal level to satisfy the assumption of LS shaping regularization. In special cases, the proposed approach will revert to the classic LS shaping regularization once the normalized coefficient is large enough. Experimental results on synthetic and field data indicate that the proposed method can effectively restore the separated records from blended data at essentially the same convergence rate as the LS shaping regularization for the abnormal noise-free scenario, but it can obtain better deblending performance and less energy leakage when abnormal noise exists.


Solid Earth ◽  
2015 ◽  
Vol 6 (1) ◽  
pp. 33-47 ◽  
Author(s):  
T. Burschil ◽  
T. Beilecke ◽  
C. M. Krawczyk

Abstract. High-resolution reflection seismic methods are an established non-destructive tool for engineering tasks. In the near surface, shear-wave reflection seismic measurements usually offer a higher spatial resolution in the same effective signal frequency spectrum than P-wave data, but data quality varies more strongly. To discuss the causes of these differences, we investigated a P-wave and a SH-wave seismic reflection profile measured at the same location on the island of Föhr, Germany and applied seismic reflection processing to the field data as well as finite-difference modelling of the seismic wave field. The simulations calculated were adapted to the acquisition field geometry, comprising 2 m receiver distance (1 m for SH wave) and 4 m shot distance along the 1.5 km long P-wave and 800 m long SH-wave profiles. A Ricker wavelet and the use of absorbing frames were first-order model parameters. The petrophysical parameters to populate the structural models down to 400 m depth were taken from borehole data, VSP (vertical seismic profile) measurements and cross-plot relations. The simulation of the P-wave wave-field was based on interpretation of the P-wave depth section that included a priori information from boreholes and airborne electromagnetics. Velocities for 14 layers in the model were derived from the analysis of five nearby VSPs (vP =1600–2300 m s-1). Synthetic shot data were compared with the field data and seismic sections were created. Major features like direct wave and reflections are imaged. We reproduce the mayor reflectors in the depth section of the field data, e.g. a prominent till layer and several deep reflectors. The SH-wave model was adapted accordingly but only led to minor correlation with the field data and produced a higher signal-to-noise ratio. Therefore, we suggest to consider for future simulations additional features like intrinsic damping, thin layering, or a near-surface weathering layer. These may lead to a better understanding of key parameters determining the data quality of near-surface shear-wave seismic measurements.


Geophysics ◽  
1985 ◽  
Vol 50 (1) ◽  
pp. 37-48 ◽  
Author(s):  
Ross Alan Ensley

Shear waves differ from compressional waves in that their velocity is not significantly affected by changes in the fluid content of a rock. Because of this relationship, a gas‐related compressional‐wave “bright spot” or direct hydrocarbon indicator will have no comparable shear‐wave anomaly. In contrast, a lithology‐related compressional‐wave anomaly will have a corresponding shear‐wave anomaly. Thus, it is possible to use shear‐wave seismic data to evaluate compressional‐wave direct hydrocarbon indicators. This case study presents data from Myrnam, Alberta which exhibit the relationship between compressional‐ and shear‐wave seismic data over a gas reservoir and a low‐velocity coal.


Sign in / Sign up

Export Citation Format

Share Document