Detector coupling corrections for vector infidelity of multicomponent OBC data

Geophysics ◽  
2007 ◽  
Vol 72 (3) ◽  
pp. V67-V77 ◽  
Author(s):  
James E. Gaiser
Keyword(s):  
Vertical Component ◽  
Low Frequency ◽  
Radial Component ◽  
Horizontal Component ◽  
Ocean Bottom ◽  
S Wave ◽  
S Waves ◽  
Resonance Modes ◽  
The Time Domain ◽  
Phase Distortions

Differences in the frequency response of horizontal and vertical detectors (vector infidelity) in ocean bottom cable (OBC) surveys can cause problems for multicomponent processing, such as S-wave birefringence and amplitude variation with azimuth (AVA) analyses, and combining vertical and hydrophone data for water-born multiple suppression. One source of this problem is poor detector coupling with the seabed that produces resonances and phase distortions. Coupling and data quality are generally excellent for the inline component. However, the crossline component often exhibits low-frequency resonance compared to the inline. Also, OBCs are susceptible to rotational modes about the cable axis that produce spurious S-waves on the vertical component. I derive a method for correcting the crossline and vertical components based on a model of OBC detector coupling, and design vector operators in the frequency domain from shots over many offsets and azimuths. The crossline data are corrected,relative to the inline, assuming linear polarization of early, near-offset arrivals on the radial-horizontal component. Thus, the transverse-horizontal component provides a convenient error or objective function to be minimized for operator design. Using the corrected crossline, as a model of rotational modes, leads to an estimate of spurious S-waves on the vertical component, which are adaptively subtracted. Data examples from the Gulf of Mexico and offshore Nigeria are presented to illustrate improvements in crossline frequency content and match to inline data. Typically there is [Formula: see text] reduction in error using the rms ratio of transverse-to-radial component data computed in the time domain. Suppression of spurious S-waves from the vertical component without undesirable effects of low-cut or [Formula: see text] filters is shown for prestack and poststack data. Also, vector operators indicate they contain important information related to resonance modes of crossline coupling and rotational modes associated with seabed-deployed versus buried OBCs.

Robust Empirical Time–Frequency Relations for Seismic Spectral Amplitudes, Part 1: Application to Regional S Waves in Southeastern Iran

2020 ◽  
Author(s):  
Maryam Safarshahi ◽  
Igor B. Morozov
Keyword(s):  
Seismic Waves ◽  
Joint Inversion ◽  
Vertical Component ◽  
Transverse Component ◽  
Model Parameters ◽  
S Wave ◽  
Near Surface ◽  
S Waves ◽  
Time Frequency ◽  
Frequency Independent

ABSTRACT Empirical models of geometrical-, Q-, t-star, and kappa-type attenuation of seismic waves and ground-motion prediction equations (GMPEs) are viewed as cases of a common empirical standard model describing variation of wave amplitudes with time and frequency. Compared with existing parametric and nonparametric approaches, several new features are included in this model: (1) flexible empirical parameterization with possible nonmonotonous time or distance dependencies; (2) joint inversion for time or distance and frequency dependencies, source spectra, site responses, kappas, and Q; (3) additional constraints removing spurious correlations of model parameters and data residuals with source–receiver distances and frequencies; (4) possible kappa terms for sources as well as for receivers; (5) orientation-independent horizontal- and three-component amplitudes; and (6) adaptive filtering to reduce noise effects. The approach is applied to local and regional S-wave amplitudes in southeastern Iran. Comparisons with previous studies show that conventional attenuation models often contain method-specific biases caused by limited parameterizations of frequency-independent amplitude decays and assumptions about the models, such as smoothness of amplitude variations. Without such assumptions, the frequency-independent spreading of S waves is much faster than inferred by conventional modeling. For example, transverse-component amplitudes decrease with travel time t as about t−1.8 at distances closer than 90 km and as t−2.5 beyond 115 km. The rapid amplitude decay at larger distances could be caused by scattering within the near surface. From about 90 to 115 km distances, the amplitude increases by a factor of about 3, which could be due to reflections from the Moho and within the crust. With more accurate geometrical-spreading and kappa models, the Q factor for the study area is frequency independent and exceeds 2000. The frequency-independent and Q-type attenuation for vertical-component and multicomponent amplitudes is somewhat weaker than for the horizontal components. These observations appear to be general and likely apply to other areas.

Flows in a rotating spherical shell: the equatorial case

1994 ◽  
Vol 276 ◽  
pp. 233-260 ◽  
Author(s):  
A. Colin de Verdière ◽  
R. Schopp
Keyword(s):  
Angular Momentum ◽  
Rotation Axis ◽  
Dynamic Pressure ◽  
Vertical Component ◽  
Low Frequency ◽  
Zonal Flow ◽  
Horizontal Component ◽  
Meridional Plane ◽  
Earth’S Rotation ◽  
Earth's Rotation

It is well known that the widely used powerful geostrophic equations that single out the vertical component of the Earth's rotation cease to be valid near the equator. Through a vorticity and an angular momentum analysis on the sphere, we show that if the flow varies on a horizontal scale L smaller than (Ha)1/2 (where H is a vertical scale of motion and a the Earth's radius), then equatorial dynamics must include the effect of the horizontal component of the Earth's rotation. In equatorial regions, where the horizontal plane aligns with the Earth's rotation axis, latitudinal variations of planetary angular momentum over such scales become small and approach the magnitude of its radial variations proscribing, therefore, vertical displacements to be freed from rotational constraints. When the zonal flow is strong compared to the meridional one, we show that the zonal component of the vorticity equation becomes (2Ω.Δ)u1 = g/ρ0)(∂ρ/a∂θ). This equation, where θ is latitude, expresses a balance between the buoyancy torque and the twisting of the full Earth's vorticity by the zonal flow u1. This generalization of the mid-latitude thermal wind relation to the equatorial case shows that u1 may be obtained up to a constant by integrating the ‘observed’ density field along the Earth's rotation axis and not along gravity as in common mid-latitude practice. The simplicity of this result valid in the finite-amplitude regime is not shared however by the other velocity components.Vorticity and momentum equations appropriate to low frequency and predominantly zonal flows are given on the equatorial β-plane. These equatorial results and the mid-latitude geostrophic approximation are shown to stem from an exact generalized relation that relates the variation of dynamic pressure along absolute vortex lines to the buoyancy field. The usual hydrostatic equation follows when the aspect ratio δ = H/L is such that tan θ/δ is much larger than one. Within a boundary-layer region of width (Ha)1/2 and centred at the equator, the analysis shows that the usually neglected Coriolis terms associated with the horizontal component of the Earth's rotation must be kept.Finally, some solutions of zonally homogeneous steady equatorial inertial jets are presented in which the Earth's vorticity is easily turned upside down by the shear flow and the correct angular momentum ‘Ωr2cos2(θ)+u1rCos(θ)’ contour lines close in the vertical–meridional plane.

scholarly journals Characteristics of the horizontal component of Rayleigh waves in multimode analysis of surface waves

Geophysics ◽  
2015 ◽  
Vol 80 (1) ◽  
pp. EN1-EN11 ◽  
Author(s):  
Tatsunori Ikeda ◽  
Toshifumi Matsuoka ◽  
Takeshi Tsuji ◽  
Toru Nakayama
Keyword(s):  
Surface Waves ◽  
Wave Velocity ◽  
Rayleigh Waves ◽  
Field Data ◽  
Vertical Component ◽  
Dispersion Curves ◽  
Horizontal Component ◽  
S Wave ◽  
Higher Modes ◽  
Explosive Source

In surface-wave analysis, S-wave velocity estimations can be improved by the use of higher modes of the surface waves. The vertical component of P-SV waves is commonly used to estimate multimode Rayleigh waves, although Rayleigh waves are also included in horizontal components of P-SV waves. To demonstrate the advantages of using the horizontal components of multimode Rayleigh waves, we investigated the characteristics of the horizontal and vertical components of Rayleigh waves. We conducted numerical modeling and field data analyses rather than a theoretical study for both components of Rayleigh waves. As a result of a simulation study, we found that the estimated higher modes have larger relative amplitudes in the vertical and horizontal components as the source depth increases. In particular, higher-order modes were observed in the horizontal component data for an explosive source located at a greater depth. Similar phenomena were observed in the field data acquired by using a dynamite source at 15-m depth. Sensitivity analyses of dispersion curves to S-wave velocity changes revealed that dispersion curves additionally estimated from the horizontal components can potentially improve S-wave velocity estimations. These results revealed that when the explosive source was buried at a greater depth, the horizontal components can complement Rayleigh waves estimated from the vertical components. Therefore, the combined use of the horizontal component data with the vertical component data would contribute to improving S-wave velocity estimations, especially in the case of buried explosive source signal.

Realtime detection and location of very local seismicity using SeisComp3

2020 ◽  
Author(s):  
Marcelo Bianchi ◽  
Lucas Schirbel ◽  
Alexandre Ausgusto
Keyword(s):  
Real Time ◽  
Earthquake Swarm ◽  
Regional Scale ◽  
Vertical Component ◽  
Operation Time ◽  
Window Size ◽  
P Wave ◽  
Local Network ◽  
Ocean Bottom ◽  
S Waves

<p>We put SeisComp3 to test by using it to analyze a very dense (9 squared kilometers) local network of 712, four components sensors (stations). Each station had a 3-component accelerometer and a pressure sensor deployed at the ocean bottom, close to the Brazilian platform near an oil exploration field. Noise levels were extreme. During the two months of the operation time, the network recorded an earthquake swarm sequence, and later analysis indicated more than 1000 earthquakes detected in a one-hour interval employing a coherency stacking method. While still not a common practice, real-time earthquake detection and location in this situation would be beneficial since this could support decisions while drilling or oil recovering is in place. Traditional tools as SeisComp3 are routinely used and allows for real-time detection and location along with the rapid revision of regional and teleseismic events, but are not widely adapted to work in a very local environment. Our experience so far showed that SeisComp3 efficiently handled the data volume (4 components at 500 samples per second times 712 stations) with a modern average workstation. Traditional SEG-Y data can be routinely converted and fed in real-time to SeedLink FIFO using ObsPy. Still, data must be correctly rotated since SeisComp3 needs at least a vertical component. Processing workflow included parallel picking using scautopick with STA/LTA, nucleation of origins using scautoloc, and location using Locsat and Hypo71 tools. In this harsh environment, the optimal window size for STA is about the size of the P-wave (0.05-0.1 s) and, LTA is about 30-60 times the S-P times (60-120 s). Using those parameters, SeisComp3 managed to generate from 400-1200 readings per data channel. We fed all picks into scautoloc that handled origin nucleation and location. Despite parameters supplied to scautoloc, the tool has many limits and relations hardcoded that inhibit it from respecting maximum requested residuals. In other words, its nucleation algorithm is adapted to work on the teleseismic and regional scale. Actual results indicate that we were able to nucleate and locate only 10-20% of known origins. Due to the flexibility of the tool, we also developed a pipeline using S-waves only. S-waves had a higher SNR for the events of interest and, due to lower velocities, presents a larger moveout on the small array easing the location. Manually picked and relocated detections returned an RMS lower as 0.04 s. Additional tests performed using the Scanloc module (GEMPA closed source nucleator) showed a higher performance during the nucleation of new origins. In this case, Hypo71 was the used locator. We did not observe any clear difference between LocSat and Hypo71 performance once the earthquake source is nucleated, and a proper velocity model is supplied.</p>

Extraction of P‐ and S‐waves from the vertical component of the particle velocity at the sea floor

Geophysics ◽  
1995 ◽  
Vol 60 (1) ◽  
pp. 231-240 ◽  
Author(s):  
Lasse Amundsen ◽  
Arne Reitan
Keyword(s):  
Particle Velocity ◽  
Synthetic Data ◽  
Vertical Component ◽  
Transversely Isotropic ◽  
Vertical Axis ◽  
Elastic Stiffness ◽  
P Wave ◽  
S Wave ◽  
Sea Floor ◽  
S Waves

At the boundary between two solid media in welded contact, all three components of particle velocity and vertical traction are continuous through the boundary. Across the boundary between a fluid and a solid, however, only the vertical component of particle velocity is continuous while the horizontal components can be discontinuous. Furthermore, the pressure in the fluid is the negative of the vertical component of traction in the solid, while the horizontal components of traction vanish at the interface. Taking advantage of this latter fact, we show that total P‐ and S‐waves can be computed from the vertical component of the particle velocity recorded by single component geophones planted on the sea floor. In the case when the sea floor is transversely isotropic with a vertical axis of symmetry, the computation requires the five independent elastic stiffness components and the density. However, when the sea floor material is fully isotropic, the only material parameter needed is the local shear wave velocity. The analysis of the extraction problem is done in the slowness domain. We show, however, that the S‐wave section can be obtained by a filtering operation in the space‐frequency domain. The P‐wave section is then the difference between the vertical component of the particle velocity and the S‐wave component. A synthetic data example demonstrates the performance of the algorithm.

Calibrating the 2016 IRIS Wavefields Experiment Nodal Sensors for Amplitude Statics and Orientation Errors

2021 ◽  
Author(s):  
Oluwaseyi J. Bolarinwa ◽  
Charles A. Langston
Keyword(s):  
Signal To Noise Ratio ◽  
Vertical Component ◽  
Calibration Method ◽  
Horizontal Component ◽  
High Signal ◽  
Correction Factors ◽  
Small Aperture ◽  
S Waves ◽  
Amplitude Correction ◽  
Array Calibration

ABSTRACT We used teleseismic P and S waves recorded in the course of the 2016 Incorporated Research Institutions for Seismology (IRIS) community-planned experiment in northern Oklahoma, to estimate amplitude correction factors (ACFs) and orientation correction factors (OCFs) for the gradiometer’s three-component Fairfield nodal sensors and two other gradiometer-styled subarray nodal sensors. These subarrays were embedded in the 13 km aperture nodal array that was also fielded during the 2016 IRIS experiment. The array calibration method we used in this study is based on the premise that a common wavefield should be recorded over a small-aperture array using teleseismic observation. In situ estimates of ACF for the gradiometer vary by 2.3% (standard deviation) for the vertical components and, typically, variability is less than 4.3% for the horizontal components; associated OCFs generally dispersed by 3°. For the two subarrays, the vertical-component ACF usually vary up to 2.4%; their horizontal-component ACFs largely spread up to 3.6%. OCFs for the subarrays generally disperse by 6.5°. ACF and OCF estimates for the gradiometer are seen to be stable across frequency bands having high signal coherence and/or signal-to-noise ratio. Gradiometry analyses of calibrated and uncalibrated gradiometer records from a local event revealed notable improvements in accuracy of attributes obtained from analyzing the calibrated horizontal-component waveforms in the light of catalog epicenter-derived azimuth. The improved waveform relative amplitudes after calibration, coupled with the enhanced wave attribute accuracy, suggests that instrument calibration for amplitude statics and orientation errors should be encouraged prior to doing gradiometry analysis in future studies.

Adaptive decomposition of multicomponent ocean‐bottom seismic data into downgoing and upgoing P‐ and S‐waves

Geophysics ◽  
2003 ◽  
Vol 68 (3) ◽  
pp. 1091-1102 ◽  
Author(s):  
K. M. Schalkwijk ◽  
C. P. A. Wapenaar ◽  
D. J. Verschuur
Keyword(s):  
Ocean Bottom ◽  
S Wave ◽  
Depth Level ◽  
S Waves ◽  
Decomposition Scheme ◽  
Adaptive Decomposition ◽  
The Individual ◽  
Event Definition ◽  
Wavefield Decomposition ◽  
Shallow Water Depth

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.

scholarly journals Amplitude of PcP, PcS, ScS, and ScP in deep-focus earthquakes*

1953 ◽  
Vol 43 (1) ◽  
pp. 63-83
Author(s):  
Kazim Ergin
Keyword(s):  
Systematic Study ◽  
Transverse Wave ◽  
Vertical Component ◽  
Horizontal Component ◽  
P Waves ◽  
S Waves ◽  
Core Boundary ◽  
Deep Focus ◽  
The Waves ◽  
Good Agreement

abstract A systematic study has been made of the ratios of (displacementperiod) of PcP, PcS, ScS, and ScP to that of the corresponding incident wave {e.g.,(displacementperiod)PcP/(displacementperiod)P}, using intermediate and deep-focus earthquake seismograms. The results indicate that the observed ratios of the horizontal components of the waves that are reflected as P waves (i.e., PcP/P and ScP/S) and that of the vertical component of the waves that are reflected as S waves (i.e., ScS/S and PcS/P) at the mantle-core boundary are considerably larger than the theoretical ones, whereas the observed ratios of the vertical component of the first group and that of the horizontal component of the second group are in fairly good agreement with the theoretical values. Theoretical computations were based on the assumption that in the case of a longitudinal wave the vibration is in the direction of propagation and in the case of a transverse wave the vibration is perpendicular to the direction of propagation. It is further found that the behavior of the direct P and S waves is in accord with the theory, but the vibration of the ground is not in the direction of propagation for PcP and ScP and is not perpendicular to the direction of propagation for PcS and ScS.

Theoretical seismic wave radiation from a fluid‐filled borehole

Geophysics ◽  
1982 ◽  
Vol 47 (9) ◽  
pp. 1308-1314 ◽  
Author(s):  
Myung W. Lee ◽  
A. H. Balch
Keyword(s):  
Radiation Pattern ◽  
Wave Velocity ◽  
Surrounding Medium ◽  
Low Frequency ◽  
P Wave ◽  
Wave Radiation ◽  
S Wave ◽  
S Waves ◽  
Field Radiation ◽  
Tube Wave

This paper concerns far‐field radiation of compressional P and shear S waves into the surrounding medium from a fluid‐filled borehole in an infinite medium and tube waves propagating along a borehole, using a low‐frequency approximation. Two kinds of sources are considered: (1) a volume displacement point source acting on the axis of a borehole, and (2) a uniform radial stress source acting on the wall of a borehole. When the tube‐wave velocity is close to the shear‐wave velocity, the effect of the borehole fluid on the P‐wave radiation pattern and on the S‐wave radiation pattern is substantial.

scholarly journals Seismic Ambient Noise Imaging of a Quasi-Amagmatic Ultra-Slow Spreading Ridge

2021 ◽  
Vol 13 (14) ◽  
pp. 2811
Author(s):  
Mohamadhasan Mohamadian Sarvandani ◽  
Emanuel Kästle ◽  
Lapo Boschi ◽  
Sylvie Leroy ◽  
Mathilde Cannat
Keyword(s):  
Wave Velocity ◽  
Vertical Component ◽  
Velocity Model ◽  
Dispersion Curves ◽  
Southwest Indian Ridge ◽  
Spreading Ridge ◽  
Ocean Bottom ◽  
S Wave ◽  
Crust And Upper Mantle ◽  
Phase Velocity Dispersion

Passive seismic interferometry has become very popular in recent years in exploration geophysics. However, it has not been widely applied in marine exploration. The purpose of this study is to investigate the internal structure of a quasi-amagmatic portion of the Southwest Indian Ridge by interferometry and to examine the performance and reliability of interferometry in marine explorations. To reach this goal, continuous vertical component recordings from 43 ocean bottom seismometers were analyzed. The recorded signals from 200 station pairs were cross-correlated in the frequency domain. The Bessel function method was applied to extract phase–velocity dispersion curves from the zero crossings of the cross-correlations. An average of all the dispersion curves was estimated in a period band 1–10 s and inverted through a conditional neighborhood algorithm which led to the final 1D S-wave velocity model of the crust and upper mantle. The obtained S-wave velocity model is in good agreement with previous geological and geophysical studies in the region and also in similar areas. We find an average crustal thickness of 7 km with a shallow layer of low shear velocities and high Vp/Vs ratio. We infer that the uppermost 2 km are highly porous and may be strongly serpentinized.

Sign in / Sign up

Export Citation Format

Share Document