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.

Scalar reverse‐time depth migration of prestack elastic seismic data

Geophysics ◽  
2001 ◽  
Vol 66 (5) ◽  
pp. 1519-1527 ◽  
Author(s):  
Robert Sun ◽  
George A. McMechan
Keyword(s):  
Seismic Waves ◽  
Velocity Model ◽  
P Wave ◽  
Common Source ◽  
Reverse Time ◽  
S Wave ◽  
Wave Source ◽  
Reflected Waves ◽  
S Waves ◽  
Elastic Data

Reflected P‐to‐P and P‐to‐S converted seismic waves in a two‐component elastic common‐source gather generated with a P‐wave source in a two‐dimensional model can be imaged by two independent scalar reverse‐time depth migrations. The inputs to migration are pure P‐ and S‐waves that are extracted by divergence and curl calculations during (shallow) extrapolation of the elastic data recorded at the earth’s surface. For both P‐to‐P and P‐to‐S converted reflected waves, the imaging time at each point is the P‐wave traveltime from the source to that point. The extracted P‐wave is reverse‐time extrapolated and imaged with a P‐velocity model, using a finite difference solution of the scalar wave equation. The extracted S‐wave is reverse‐time extrapolated and imaged similarly, but with an S‐velocity model. Converted S‐wave data requires a polarity correction prior to migration to ensure constructive interference between data from adjacent sources. Synthetic examples show that the algorithm gives satisfactory results for laterally inhomogeneous models.

On-axis triplications in elastic orthorhombic media

2020 ◽  
Vol 224 (1) ◽  
pp. 449-467
Author(s):  
Shibo Xu ◽  
Alexey Stovas ◽  
Hitoshi Mikada ◽  
Junichi Takekawa
Keyword(s):  
Group Velocity ◽  
Transversely Isotropic ◽  
Model Complexity ◽  
Model Parameters ◽  
S Wave ◽  
Slowness Surface ◽  
S Waves ◽  
Redundant Signal ◽  
Velocity Surface ◽  
Offset Curve

SUMMARY Triplicated traveltime curve has three arrivals at a given distance with the bowtie shape in the traveltime-offset curve. The existence of the triplication can cause a lot of problems such as several arrivals for the same wave type, anomalous amplitudes near caustics, anomalous behaviour of rays near caustics, which leads to the structure imaging deviation and redundant signal in the inversion of the model parameters. Hence, triplication prediction becomes necessary when the medium is known. The research of the triplication in transversely isotropic medium with a vertical symmetry axis (VTI) has been well investigated and it has become clear that, apart from the point singularity case, the triplicated traveltime only occurs for S wave. On contrary to the VTI case, the triplication behaviour in the orthorhombic (ORT) medium has not been well focused due to the model complexity. In this paper, we derive the second-order coefficients of the slowness surface for two S waves in the vicinity of three symmetry axes and define the elliptic form function to examine the existence of the on-axis triplication in ORT model. The existence of the on-axis triplication is found by the sign of the defined curvature coefficients. Three ORT models are defined in the numerical examples to analyse the behaviour of the on-axis triplication. The plots of the group velocity surface in the vicinity of three symmetry axes are shown for different ORT models where different shapes: convex or the saddle-shaped (concave along one direction and convex along with another) indicates the existence of the on-axis triplication. We also show the traveltime plots (associated with the group velocity surface) to illustrate the effect of the on-axis triplication.

scholarly journals A new acoustic assumption for orthorhombic media

2020 ◽  
Vol 223 (2) ◽  
pp. 1118-1129
Author(s):  
Mohammad Mahdi Abedi ◽  
Alexey Stovas
Keyword(s):  
Wave Propagation ◽  
Conventional Approach ◽  
P Wave ◽  
Model Parameters ◽  
Governing Equations ◽  
S Wave ◽  
S Waves ◽  
Exploration Seismology ◽  
And Inversion ◽  
Approximate Equations

SUMMARY In exploration seismology, the acquisition, processing and inversion of P-wave data is a routine. However, in orthorhombic anisotropic media, the governing equations that describe the P-wave propagation are coupled with two S waves that are considered as redundant noise. The main approach to free the P-wave signal from the S-wave noise is the acoustic assumption on the wave propagation. The conventional acoustic assumption for orthorhombic media zeros out the S-wave velocities along three orthogonal axes, but leaves significant S-wave artefacts in all other directions. The new acoustic assumption that we propose mitigates the S-wave artefacts by zeroing out their velocities along the three orthogonal symmetry planes of orthorhombic media. Similar to the conventional approach, our method reduces the number of required model parameters from nine to six. As numerical experiments on multiple orthorhombic models show, the accuracy of the new acoustic assumption also compares well to the conventional approach. On the other hand, while the conventional acoustic assumption simplifies the governing equations, the new acoustic assumption further complicates them—an issue that emphasizes the necessity of simple approximate equations. Accordingly, we also propose simpler rational approximate phase-velocity and eikonal equations for the new acoustic orthorhombic media. We show a simple ray tracing example and find out that the proposed approximate equations are still highly accurate.

Velocity anisotropy in shale determined from crosshole seismic and vertical seismic profile data

Geophysics ◽  
1990 ◽  
Vol 55 (4) ◽  
pp. 470-479 ◽  
Author(s):  
D. F. Winterstein ◽  
B. N. P. Paulsson
Keyword(s):  
Seismic Profile ◽  
P Wave ◽  
Isotropic Model ◽  
S Wave ◽  
Vertical Seismic Profile ◽  
Near Surface ◽  
S Waves ◽  
Velocity Anisotropy ◽  
Wave Velocities ◽  
Late Tertiary

Crosshole and vertical seismic profile (VST) data made possible accurate characterization of the elastic properties, including noticeable velocity anisotropy, of a near‐surface late Tertiary shale formation. Shear‐wave splitting was obvious in both crosshole and VSP data. In crosshole data, two orthologonally polarrized shear (S) waves arrived 19 ms in the uppermost 246 ft (75 m). Vertically traveling S waves of the VSP separated about 10 ms in the uppermost 300 ft (90 m) but remained at nearly constant separation below that level. A transversely isotropic model, which incorporates a rapid increase in S-wave velocities with depth but slow increase in P-wave velocities, closely fits the data over most of the measured interval. Elastic constants of the transvesely isotropic model show spherical P- and [Formula: see text]wave velocity surfaces but an ellipsoidal [Formula: see text]wave surface with a ratio of major to minor axes of 1.15. The magnitude of this S-wave anisotropy is consistent with and lends credence to S-wave anisotropy magnitudes deduced less directly from data of many sedimentary basins.

Crustal structures across the young and oblique North-eastern Gulf of Aden margin

2020 ◽  
Author(s):  
Louise Watremez ◽  
Sylvie Leroy ◽  
Elia d'Acremont ◽  
Stéphane Rouzo
Keyword(s):  
Seismic Waves ◽  
Seismic Refraction ◽  
P Wave ◽  
Gulf Of Aden ◽  
S Wave ◽  
Acoustic Basement ◽  
S Waves ◽  
P Wave Velocity ◽  
North Eastern ◽  
Lower Crustal

<p>The Gulf of Aden is a young and active oceanic basin, which separates the south-eastern margin of the Arabian Plate from the Somali Plate. The rifting leading to the formation of the north-eastern Gulf of Aden passive margin started ca. 34 Ma ago when the oceanic spreading in this area initiated at least 17.6 Ma ago. The opening direction (N26°E) is oblique to the mean orientation of the Gulf (N75°E), leading to a strong structural segmentation.</p><p>The Encens cruise (2006) allowed for the acquisition of a large seismic refraction dataset with profiles across (6 lines) and along (3 lines) the margin, between the Alula-Fartak and Socotra-Hadbeen fracture zones, which define a first order segment of the Gulf. P-wave velocity modelling already allowed us to image the crustal thinning and the structures, from continental to oceanic domains, along some of the profiles. A lower crustal intermediate body is observed in the Ashawq-Salalah segment, at the base of the transitional and oceanic crusts. The nature of this intermediate body is most probably mafic, linked to a post-rift thermal anomaly. The thin (1-2 km) sediment layer in the study area allows for a clear conversion of P-waves to S-waves at the top basement. Thus, most seismic refraction records show very clear S-wave arrivals.</p><p>In this study, we use both P-wave and S-wave arrivals to delineate the crustal structures and segmentation along and across the margin and add insight into the nature of the rocks below the acoustic basement. P-wave velocity modelling allows for the delineation of the structure variations across and along the margin. The velocity models are used as a base for the S-wave modelling, through the definition of Poisson’s ratios in the different areas of the models. Picking and modelling of S-wave arrivals allow us to identify two families of converted waves: (1) seismic waves converted at the basement interface on the way up, just before arriving to the OBS and (2) seismic waves converted at the basement on the way down, which travelled into the deep structures as S-waves. The first set of arrivals allows for the estimation the S-wave velocities (Poisson’s ratio) in the sediments, showing that the sediments in this area are unconsolidated and water saturated. The second set of arrivals gives us constraints on the S-wave velocities below the acoustic basement. This allows for an improved mapping of the transitional and oceanic domains and the confirmation of the mafic nature of the lower crustal intermediate body.</p>

Receiver function analyses and joint inversion with gravity data: new constraints on the Ivrea geophysical body along a high-resolution profile

2020 ◽  
Author(s):  
Matteo Scarponi ◽  
György Hetényi ◽  
Jaroslava Plomerová ◽  
Stefano Solarino ◽  
Ludovic Baron
Keyword(s):  
Receiver Function ◽  
Joint Inversion ◽  
Gravity Data ◽  
Sampling Rate ◽  
Seismic Refraction ◽  
Western Alps ◽  
Transverse Component ◽  
Reference Points ◽  
Physical Parameters ◽  
S Wave

<p>We collected new seismological and gravity data in the Val Sesia and Lago Maggiore regions in NW Italy to constrain the geometry and properties of the Ivrea Geophysical Body. This piece of lower Adriatic lithosphere is known to be at anomalously shallow depth along the inner arc of the Western Alps, yet existing seismological constraints (vintage seismic refraction data, local earthquake tomography) are spatially sparse. With the aim to reach higher spatial resolution in imaging the structure of the IGB, we analyze the seismological data with various receiver function approaches to map the main velocity discontinuities, followed by joint inversion with gravity data to fill the bulk properties of bodies with densities.</p><p>The new data acquisition consisted of two type of campaigns. For seismology, we deployed 10 broadband seismic stations (MOBNET pool, IG CAS Prague) along a linear West-East profile at 5 km spacing along Val Sesia and across the Lago Maggiore. This network continuously recorded seismic data for 27 months at 100 Hz sampling rate. For gravimetry, we compiled existing datasets and then completed the spatial gaps by relative gravity surveys, tied to absolute reference points, to achieve 1 gravity point every 1-2 km along the profile.</p><p>The receiver function (RF) analyses aim at detecting velocity increases with depth: primarily the Moho and the shallow IGB interfaces and their crustal reverberations (multiples), together with their potential dip by analyzing the transverse component RFs. Furthermore, we aim at investigating the sharpness of the velocity gradient across the discontinuities by analyzing the frequency dependence of the corresponding RF peaks. We aim at reproducing the observations by simple synthetic models.</p><p>The 2D joint inversion combines S wave velocity V<sub>S</sub> and bulk density as physical parameters to match both the seismological and gravimetry data. The relationship between the two parameters is initially chosen from the literature, but depending on the first results the relation itself may be inverted for, considering the various high-grade metamorphic rocks observed at the surface in the area, whose properties may not align with classical V<sub>S</sub>–density equations. In conclusion, we propose new constraints on the IGB, demonstrating the advantage of using multi-disciplinary geophysical observations and improved data coverage across the study area.</p>

Robust Empirical Time–Frequency Relations for Seismic Spectral Amplitudes, Part 2: Model Uncertainty and Optimal Parameterization

2020 ◽  
Author(s):  
Maryam Safarshahi ◽  
Igor B. Morozov
Keyword(s):  
Measurement Errors ◽  
Nonparametric Model ◽  
Model Parameters ◽  
Geometrical Spreading ◽  
Model Quality ◽  
Frequency Dependent ◽  
Time Frequency ◽  
Trade Offs ◽  
Frequency Independent ◽  
Seismic Amplitudes

ABSTRACT In a companion article, Safarshahi and Morozov (2020) argued that construction of distance- and frequency-dependent models for seismic-wave amplitudes should include four general elements: (1) a sufficiently detailed (parametric or nonparametric) model of frequency-independent spreading, capturing all essential features of observations; (2) model parameters with well-defined and nonoverlapping physical meanings; (3) joint inversion for multiple parameters, including the geometrical spreading, Q, κ, and source and receiver couplings; and (4) the use of additional dataset-specific criteria of model quality, while fitting the logarithms of seismic amplitudes. Some of these elements are present in existing models, but, taken together, they are poorly understood and require an integrated approach. Such an approach was illustrated by detailed analysis of an S-wave amplitude dataset from southern Iran. The resulting model is based on a frequency-independent Q, and matches the data closer than conventional models and across the entire epicentral-distance range. Here, we complete the analysis of this model by evaluating the uncertainties and trade-offs of its parameters. Two types of trade-offs are differentiated: one caused by a (possibly) limited model parameterization and the second due to statistical data errors. Data bootstrapping shows that with adequate parameterization, attenuation properties Q, κ, and geometrical spreading parameters are resolved well and show moderate trade-offs due to measurement errors. Using the principal component analysis of these trade-offs, an optimal (trade-off free) parameterization of seismic amplitudes is obtained. By contrast, when assuming theoretical values for certain model parameters and using multistep inversion procedures (as commonly done), parameter trade-offs increase dramatically and become difficult to assess. In particular, the frequency-dependent Q correlates with the distribution of the source and receiver-site factors, and also with biases in the resulting median data residuals. In the new model, these trade-offs are removed using an improved parameterization of geometrical spreading, constant Q, and model quality constraints.

Thin-bed prestack spectral inversion

Geophysics ◽  
2009 ◽  
Vol 74 (4) ◽  
pp. R49-R57 ◽  
Author(s):  
J. Germán Rubino ◽  
Danilo Velis
Keyword(s):  
Time Window ◽  
A Priori ◽  
Amplitude Spectrum ◽  
Wavelet Spectrum ◽  
Model Parameters ◽  
Data Sets ◽  
S Wave ◽  
P Waves ◽  
S Waves ◽  
Zero Offset

Prestack seismic data has been used in a new method to fully determine thin-bed properties, including the estimation of its thickness, P- and S-wave velocities, and density. The approach requires neither phase information nor normal-moveout (NMO) corrections, and assumes that the prestack seismic response of the thin layer can be isolated using an offset-dependent time window. We obtained the amplitude-versus-angle (AVA) response of the thin bed considering converted P-waves, S-waves, and all the associated multiples. We carried out the estimation of the thin-bed parameters in the frequency (amplitude spectrum) domain using simulated annealing. In contrast to using zero-offset data, the use of AVA data contributes to increase the robustness of this inverse problem under noisy conditions, as well as to significantly reduce its inherent nonuniqueness. To further reduce the nonuniqueness, and as a means to incorporate a priori geologic or geophysical information (e.g., well-log data), we imposed appropriate bounding constraints to the parameters of the media lying above and below the thin bed, which need not be known accurately. We tested the method by inverting noisy synthetic gathers corresponding to simple wedge models. In addition, we stochastically estimated the uncertainty of the solutions by inverting different data sets that share the same model parameters but are contaminated with different noise realizations. The results suggest that thin beds can be characterized fully with a moderate to high degree of confidence below tuning, even when using an approximate wavelet spectrum.

scholarly journals Converted‐wave seismic exploration: Applications

Geophysics ◽  
2003 ◽  
Vol 68 (1) ◽  
pp. 40-57 ◽  
Author(s):  
Robert R. Stewart ◽  
James E. Gaiser ◽  
R. James Brown ◽  
Don C. Lawton
Keyword(s):  
Seismic Waves ◽  
Seismic Exploration ◽  
Fracture Density ◽  
Converted Wave ◽  
P Waves ◽  
Near Surface ◽  
S Waves ◽  
Subsurface Fluid ◽  
Processing Techniques ◽  
Gas Bearing

Converted seismic waves (specifically, downgoing P‐waves that convert on reflection to upcoming S‐waves are increasingly being used to explore for subsurface targets. Rapid advancements in both land and marine multicomponent acquisition and processing techniques have led to numerous applications for P‐S surveys. Uses that have arisen include structural imaging (e.g., “seeing” through gas‐bearing sediments, improved fault definition, enhanced near‐surface resolution), lithologic estimation (e.g., sand versus shale content, porosity), anisotropy analysis (e.g., fracture density and orientation), subsurface fluid description, and reservoir monitoring. Further applications of P‐S data and analysis of other more complicated converted modes are developing.

Sign in / Sign up

Export Citation Format

Share Document