A new method to estimate ocean-bottom-seismometer orientation using teleseismic receiver functions

2020 ◽  
Vol 221 (2) ◽  
pp. 893-904
Author(s):  
Hong Zheng ◽  
Jianke Fan ◽  
Dapeng Zhao ◽  
Cuilin Li ◽  
Dongdong Dong ◽  
...  
Keyword(s):  
Epicentral Distance ◽  
Synthetic Data ◽  
Receiver Functions ◽  
P Wave ◽  
New Method ◽  
Ocean Bottom ◽  
Ocean Bottom Seismometer ◽  
S Wave ◽  
P Waves ◽  
The Western Pacific

SUMMARY The orientation of an ocean-bottom-seismometer (OBS) is a critical parameter for analysing three-component seismograms, but it is difficult to estimate because of the uncontrollable OBS posture after its deployment. In this study, we develop a new and effective method to estimate the OBS orientation by fitting the amplitude of direct P wave of teleseismic receiver functions. The reliability of this method is verified using synthetic data and observed waveforms recorded at land seismic stations in Shandong Province, China. Our extensive synthetic tests show that our new method is little affected by a thin sedimentary layer that has a low S-wave velocity. The orientations of OBS stations that we deployed in the Yap subduction zone in the Western Pacific Ocean are estimated and corrected using our new method. After the correction, the direct P waves of teleseismic receiver functions show very good consistency. The effects of white and coloured noise in different levels, epicentral distance and backazimuth are also investigated, and the results show that these factors have small effects on the new method. We also examine the effect of sensor tilting on estimation of the OBS orientation, and find that a tilting correction should be made before the misorientation correction. We compare the OBS orientations determined with the new method and other methods and find that they are generally consistent with each other. We also discuss advantages and shortcomings of various methods, and think that our new method is more robust than the existing methods.

scholarly journals OBS Data Analysis to Quantify Gas Hydrate and Free Gas in the South Shetland Margin (Antarctica)

Energies ◽  
2018 ◽  
Vol 11 (12) ◽  
pp. 3290 ◽  
Author(s):  
Sha Song ◽  
Umberta Tinivella ◽  
Michela Giustiniani ◽  
Sunny Singhroha ◽  
Stefan Bünz ◽  
...  
Keyword(s):  
Velocity Field ◽  
Wave Velocity ◽  
Gas Hydrate ◽  
P Wave ◽  
Ocean Bottom ◽  
Ocean Bottom Seismometer ◽  
S Wave ◽  
Free Gas ◽  
Hydrate Reservoir ◽  
Gas Hydrate Reservoir

The presence of a gas hydrate reservoir and free gas layer along the South Shetland margin (offshore Antarctic Peninsula) has been well documented in recent years. In order to better characterize gas hydrate reservoirs, with a particular focus on the quantification of gas hydrate and free gas and the petrophysical properties of the subsurface, we performed travel time inversion of ocean-bottom seismometer data in order to obtain detailed P- and S-wave velocity estimates of the sediments. The P-wave velocity field is determined by the inversion of P-wave refractions and reflections, while the S-wave velocity field is obtained from converted-wave reflections received on the horizontal components of ocean-bottom seismometer data. The resulting velocity fields are used to estimate gas hydrate and free gas concentrations using a modified Biot‐Geertsma‐Smit theory. The results show that hydrate concentration ranges from 10% to 15% of total volume and free gas concentration is approximately 0.3% to 0.8% of total volume. The comparison of Poisson’s ratio with previous studies in this area indicates that the gas hydrate reservoir shows no significant regional variations.

scholarly journals Post-collisional mantle delamination in the Dinarides implied from staircases of Oligo-Miocene uplifted marine terraces

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Philipp Balling ◽  
Christoph Grützner ◽  
Bruno Tomljenović ◽  
Wim Spakman ◽  
Kamil Ustaszewski
Keyword(s):  
Balkan Peninsula ◽  
Receiver Functions ◽  
P Wave ◽  
S Wave ◽  
Slab Geometry ◽  
Surface Uplift ◽  
River Incision ◽  
Marine Terraces ◽  
Convergence System ◽  
Internal Dinarides

AbstractThe Dinarides fold-thrust belt on the Balkan Peninsula resulted from convergence between the Adriatic and Eurasian plates since Mid-Jurassic times. Under the Dinarides, S-wave receiver functions, P-wave tomographic models, and shear-wave splitting data show anomalously thin lithosphere overlying a short down-flexed slab geometry. This geometry suggests a delamination of Adriatic lithosphere. Here, we link the evolution of this continental convergence system to hitherto unreported sets of extensively uplifted Oligocene–Miocene (28–17 Ma) marine terraces preserved at elevations of up to 600 m along the Dinaric coastal range. River incision on either side of the Mediterranean-Black Sea drainage divide is comparable to the amounts of terrace uplift. The preservation of the uplifted terraces implies that the most External Dinarides did not experience substantial deformation other than surface uplift in the Neogene. These observations and the contemporaneous emplacement of igneous rocks (33–22 Ma) in the internal Dinarides suggest that the Oligo-Miocene orogen-wide uplift was driven by post-break-off delamination of the Adriatic lithospheric mantle, this was followed by isostatic readjustment of the remaining crust. Our study details how lithospheric delamination exerts an important control on crustal deformation and that its crustal signature and geomorphic imprint can be preserved for millions of years.

Numerical simulation of azimuthal acoustic logging in a borehole penetrating a rock formation boundary

Geophysics ◽  
2016 ◽  
Vol 81 (3) ◽  
pp. D283-D291 ◽  
Author(s):  
Peng Liu ◽  
Wenxiao Qiao ◽  
Xiaohua Che ◽  
Xiaodong Ju ◽  
Junqiang Lu ◽  
...  
Keyword(s):  
Domain Model ◽  
P Wave ◽  
S Wave ◽  
Angle Variation ◽  
P Waves ◽  
Acoustic Logging ◽  
S Waves ◽  
Rock Formation ◽  
Logging Tool ◽  
Difference Time

We have developed a new 3D acoustic logging tool (3DAC). To examine the azimuthal resolution of 3DAC, we have evaluated a 3D finite-difference time-domain model to simulate a case in which the borehole penetrated a rock formation boundary when the tool worked at the azimuthal-transmitting-azimuthal-receiving mode. The results indicated that there were two types of P-waves with different slowness in waveforms: the P-wave of the harder rock (P1) and the P-wave of the softer rock (P2). The P1-wave can be observed in each azimuthal receiver, but the P2-wave appears only in the azimuthal receivers toward the softer rock. When these two types of rock are both fast formations, two types of S-waves also exist, and they have better azimuthal sensitivity compared with P-waves. The S-wave of the harder rock (S1) appears only in receivers toward the harder rock, and the S-wave of the softer rock (S2) appears only in receivers toward the softer rock. A model was simulated in which the boundary between shale and sand penetrated the borehole but not the borehole axis. The P-wave of shale and the S-wave of sand are azimuthally sensitive to the azimuth angle variation of two formations. In addition, waveforms obtained from 3DAC working at the monopole-transmitting-azimuthal-receiving mode indicate that the corresponding P-waves and S-waves are azimuthally sensitive, too. Finally, we have developed a field example of 3DAC to support our simulation results: The azimuthal variation of the P-wave slowness was observed and can thus be used to reflect the azimuthal heterogeneity of formations.

Multi-parameter model building for the Qiuyue structure using four-component ocean-bottom seismometer data

Geophysics ◽  
2021 ◽  
pp. 1-52
Author(s):  
Yuzhu Liu ◽  
Xinquan Huang ◽  
Jizhong Yang ◽  
Xueyi Liu ◽  
Bin Li ◽  
...  
Keyword(s):  
Wave Velocity ◽  
Waveform Inversion ◽  
Velocity Model ◽  
P Wave ◽  
Velocity Analysis ◽  
Full Waveform Inversion ◽  
Ocean Bottom ◽  
Ocean Bottom Seismometer ◽  
Full Waveform ◽  
P Wave Velocity

Thin sand-mud-coal interbedded layers and multiples caused by shallow water pose great challenges to conventional 3D multi-channel seismic techniques used to detect the deeply buried reservoirs in the Qiuyue field. In 2017, a dense ocean-bottom seismometer (OBS) acquisition program acquired a four-component dataset in East China Sea. To delineate the deep reservoir structures in the Qiuyue field, we applied a full-waveform inversion (FWI) workflow to this dense four-component OBS dataset. After preprocessing, including receiver geometry correction, moveout correction, component rotation, and energy transformation from 3D to 2D, a preconditioned first-arrival traveltime tomography based on an improved scattering integral algorithm is applied to construct an initial P-wave velocity model. To eliminate the influence of the wavelet estimation process, a convolutional-wavefield-based objective function for the preprocessed hydrophone component is used during acoustic FWI. By inverting the waveforms associated with early arrivals, a relatively high-resolution underground P-wave velocity model is obtained, with updates at 2.0 km and 4.7 km depth. Initial S-wave velocity and density models are then constructed based on their prior relationships to the P-wave velocity, accompanied by a reciprocal source-independent elastic full-waveform inversion to refine both velocity models. Compared to a traditional workflow, guided by stacking velocity analysis or migration velocity analysis, and using only the pressure component or other single-component, the workflow presented in this study represents a good approach for inverting the four-component OBS dataset to characterize sub-seafloor velocity structures.

On the Green function of the Lord–Shulman thermoelasticity equations

2019 ◽  
Vol 220 (1) ◽  
pp. 393-403 ◽  
Author(s):  
Zhi-Wei Wang ◽  
Li-Yun Fu ◽  
Jia Wei ◽  
Wanting Hou ◽  
Jing Ba ◽  
...  
Keyword(s):  
Green’S Function ◽  
Green's Function ◽  
Parabolic Type ◽  
Second Order ◽  
P Wave ◽  
Thermal Mode ◽  
S Wave ◽  
Order Tensor ◽  
P Waves ◽  
Thermal Conductivities

SUMMARY Thermoelasticity extends the classical elastic theory by coupling the fields of particle displacement and temperature. The classical theory of thermoelasticity, based on a parabolic-type heat-conduction equation, is characteristic of an unphysical behaviour of thermoelastic waves with discontinuities and infinite velocities as a function of frequency. A better physical system of equations incorporates a relaxation term into the heat equation; the equations predict three propagation modes, namely, a fast P wave (E wave), a slow thermal P wave (T wave), and a shear wave (S wave). We formulate a second-order tensor Green's function based on the Fourier transform of the thermodynamic equations. It is the displacement–temperature solution to a point (elastic or heat) source. The snapshots, obtained with the derived second-order tensor Green's function, show that the elastic and thermal P modes are dispersive and lossy, which is confirmed by a plane-wave analysis. These modes have similar characteristics of the fast and slow P waves of poroelasticity. Particularly, the thermal mode is diffusive at low thermal conductivities and becomes wave-like for high thermal conductivities.

Near-surface scattering phenomena and implications for tunnel detection

2015 ◽  
Vol 3 (1) ◽  
pp. SF43-SF54 ◽  
Author(s):  
Shelby L. Peterie ◽  
Richard D. Miller
Keyword(s):  
Synthetic Data ◽  
P Wave ◽  
S Wave ◽  
Seismic Line ◽  
Data Set ◽  
Near Surface ◽  
Tunnel Detection ◽  
Impedance Contrast ◽  
Diffraction Imaging ◽  
Imaging Algorithms

Tunnel locations are accurately interpreted from diffraction sections of focused mode converted P- to S-wave diffractions from a perpendicular tunnel and P-wave diffractions from a nonperpendicular (oblique) tunnel. Near-surface tunnels are ideal candidates for diffraction imaging due to their small size relative to the seismic wavelength and large acoustic impedance contrast at the tunnel interface. Diffraction imaging algorithms generally assume that the velocities of the primary wave and the diffracted wave are approximately equal, and that the diffraction apex is recorded directly above the scatterpoint. Scattering phenomena from shallow tunnels with kinematic properties that violate these assumptions were observed in one field data set and one synthetic data set. We developed the traveltime equations for mode-converted and oblique diffractions and demonstrated a diffraction imaging algorithm designed for the roll-along style of acquisition. Potential processing and interpretation pitfalls specific to these diffraction types were identified. Based on our observations, recommendations were made to recognize and image mode-converted and oblique diffractions and accurately interpret tunnel depth, horizontal location, and azimuth with respect to the seismic line.

Ground Motion From A Magnitude 3.5 Earthquake Near Massena, New York: Evidence For Poor Resolution of Corner Frequency For Small Events

1989 ◽  
Vol 60 (3) ◽  
pp. 95-100 ◽  
Author(s):  
S.E. Hough ◽  
K. Jacob ◽  
R. Busby ◽  
P.A. Friberg
Keyword(s):  
Earthquake Engineering ◽  
Epicentral Distance ◽  
Wave Spectrum ◽  
P Wave ◽  
Corner Frequency ◽  
Source Spectrum ◽  
Wave Spectra ◽  
S Wave ◽  
High Quality ◽  
Wide Range

Abstract We present analysis of a magnitude 3.5 event which occurred at 9 km epicentral distance from a digital strong motion instrument operated by the National Center for Earthquake Engineering Research. Although the size of this isolated event is such that it can scarcely be considered to be a significant earthquake, a careful analysis of this high quality recording does yield several interesting results: 1) the S-wave spectra can be interpreted in terms of a simple omega-squared source spectrum and frequency-independent attenuation, 2) there is the suggestion of a poorly-resolved resonance in the P-wave spectrum, and perhaps most importantly, 3) the apparently simple S-wave spectra can be fit almost equally well with a surprisingly wide range of seismic corner frequencies, from roughly 5 to 25 Hz. This uncertainty in corner frequency translates into uncertainties in inferred Q values of almost an order of magnitude, and into uncertainties in stress drop of two orders of magnitude. Given the high quality of the data and the short epicentral distance to the station, we consider it likely that resolution of spectral decay and corner frequency will be at least as poor for any other recording of earthquakes with comparable or smaller magnitudes.

Seismic reflectivity of hydraulic fractures approximated by thin fluid layers

Geophysics ◽  
2013 ◽  
Vol 78 (4) ◽  
pp. T79-T87 ◽  
Author(s):  
A. Oelke ◽  
D. Alexandrov ◽  
I. Abakumov ◽  
S. Glubokovskikh ◽  
R. Shigapov ◽  
...  
Keyword(s):  
Analytical Solutions ◽  
Fluid Layer ◽  
Synthetic Data ◽  
P Wave ◽  
Hydraulic Fractures ◽  
Reflection Coefficients ◽  
S Wave ◽  
Reservoir Properties ◽  
Thin Fluid ◽  
Theoretical Results

We have analyzed the angle-dependent reflectivity of microseismic wavefields at a hydraulic fracture, which we modeled as an ideal thin fluid layer embedded in an elastic, isotropic solid rock. We derived full analytical solutions for the reflections of an incident P-wave, the P-P and P-S reflection coefficients, as well as for an incident S-wave, and the S-S and S-P reflection coefficients. The rather complex analytical solutions were then approximated and we found that these zero-thickness limit approximations are in good agreement with the linear slip model, representing a fracture at slip contact. We compared the analytical solutions for the P-P reflections with synthetic data that were derived using finite-difference modeling and found that the modeling confirmed our theoretical results. For typical parameters of microseismic monitoring by hydraulic fracturing, e.g., a layer thickness of [Formula: see text] and frequencies of [Formula: see text], the reflection coefficients depend on the Poisson’s ratio. Furthermore, the reflection coefficients of an incident S-wave are remarkably high. Theoretical results suggested that it is feasible to image hydraulic fractures using microseismic events as a source and to solve the inverse problem, that is, to interpret reflection coefficients extracted from microseismic data in terms of reservoir properties.

Sign in / Sign up

Export Citation Format

Share Document