Reflection maxima for reflections from single interfaces

Geophysics ◽  
1988 ◽  
Vol 53 (2) ◽  
pp. 271-275 ◽  
Author(s):  
C. A. Rendleman ◽  
F. K. Levin

At a workshop on refraction and wide‐angle reflections, Hilterman (1985) pointed out that, in contrast to the plane‐wave case, when there is a point source, a P-wave reflected from a plane interface attains its maximum amplitude at an offset greater than that corresponding to the critical angle (Figure 1). The same conclusion had been drawn earlier by Červený (1967). However, neither Červený’s results, which were based on very complicated mathematical expressions derived by Brekhovskikh (1960), nor Hilterman’s computer‐generated data shed light on the physics implied by the shifted maximum.

Geophysics ◽  
1984 ◽  
Vol 49 (4) ◽  
pp. 398-410 ◽  
Author(s):  
D. P. Blair

Fourier‐Bessel theory is used to derive filters representing the influence of both empty and fluid‐filled cylindrical boreholes on particle motion induced in rock by a plane P-wave incident perpendicular to the borehole axis. For wavelengths greater than 10 times the borehole circumference, the effect of the borehole on particle motions is shown to be negligible; thus the results have little relevance for the long wavelengths commonly encountered in earthquake seismology. The results are, however, relevant to the study of stress wave propagation at ultrasonic frequencies in rock masses. For small wavelengths (αa > 3.0) the filter representing particle motion on the wave incident site of an empty borehole reduces to a linear phase filter which increases all amplitudes by a factor of 2 while the filter representing fluid stress at the center of a fluid‐filled borehole may be reduced to simple mathematical expressions. Experimental results were obtained for the interaction of a stress wave with either accelerometers mounted in an empty borehole or a hydrophone located centrally in a fluid‐filled borehole. Both theory and experiment show a similar distortion in the rise time of the pulse traveling past the borehole.


Geophysics ◽  
2005 ◽  
Vol 70 (4) ◽  
pp. R45-R56 ◽  
Author(s):  
Lars Nielsen ◽  
Hans Thybo ◽  
Martin Glendrup

Seismic wide-angle data were recorded to more than 300-km offset from powerful airgun sources during the MONA LISA experiments in 1993 and 1995 to determine the seismic-velocity structure of the crust and uppermost mantle along three lines in the southeastern North Sea with a total length of 850 km. We use the first arrivals observed out to an offset of 90 km to obtain high-resolution models of the velocity structure of the sedimentary layers and the upper part of the crystalline crust. Seismic tomographic traveltime inversion reveals 2–8-km-thick Paleozoic sedimentary sequences with P-wave velocities of 4.5–5.2 km/s. These sedimentary rocks are situated below a Mesozoic-Cenozoic sequence with variable thickness: ∼2–3 km on the basement highs, ∼2–4 km in the Horn Graben and the North German Basin, and ∼6–7 km in the Central Graben. The thicknesses of the Paleozoic sedimentary sequences are ∼3–5 km in the Central Graben, more than 4 km in the Horn Graben, up to ∼4 km on the basement highs, and up to 8 km in the North German Basin. The Paleozoic strata are clearly separated from the shallower and younger sequences with velocities of ∼1.8–3.8 km/s and the deeper crystalline crust with velocities of more than 5.8–6.0 km/s in the tomographic P-wave velocity model. Resolution tests show that the existence of the Paleozoic sediments is well constrained by the data. Hence, our wide-angle seismic models document the presence of Paleozoic sediments throughout the southeastern North Sea, both in the graben structures and in deep basins on the basement highs.


2021 ◽  
Author(s):  
Nicola Piana Agostinetti ◽  
Alberto Villa ◽  
Gilberto Saccorotti

Abstract. We use PoroTOMO experimental data to compare the performance of Distributed Acoustic Sensing (DAS) and geophone data in executing standard exploration and monitoring activities. The PoroTOMO experiment consists of two "seismic systems": (a) a 8.6 km long optical fibre cable deployed across the Brady geothermal field and covering an area of 1.5 x 0.5 km with 100 m long segments, and (b) an array of 238 co-located geophones with an average spacing of 60 m. The PoroTOMO experiment recorded continuous seismic data between March 10th and March 25th 2016. During such period, a ML 4.3 regional event occurred in the southwest, about 150 km away from the geothermal field, together with several microseismic local events related to the geothermal activity. The seismic waves generated from such seismic events have been used as input data in this study. For the exploration tasks, we compare the propagation of the ML 4.3 event across the geothermal field in both seismic systems in term of relative time-delay, for a number of configurations and segments. Defined the propagation, we analyse and compare the amplitude and the signal-to-noise ratio (SNR) of the P-wave in the two systems at high resolution. For testing the potential in monitoring local seismicity, we first perform an analysis of the geophone data for locating a microseismic event, based on expert opinion. Then, we a adopt different workflow for the automatic location of the same microseismic event using DAS data. For testing the potential in monitoring distant event, data from the regional earthquake are used for retrieving both the propagation direction and apparent velocity of the wavefield, using a standard plane-wave-fitting approach. Our results indicate that: (1) at a local scale, the seismic P-waves propagation and their characteristics (i.e. SNR and amplitude) along a single cable segment are robustly consistent with recordings from co-located geophones (delay-times δt ∼ 0.3 over 400 m for both seismic systems) ; (2) the interpretation of seismic wave propagation across multiple separated segments is less clear, due to the heavy contamination of scattering sources and local velocity heterogeneities; nonetheless, results from the plane-wave fitting still indicate the possibility for a consistent detection and location of the event; (3) at high-resolution (10 m), large amplitude variations along the fibre cable seem to robustly correlate with near surface geology; (4) automatic monitoring of microseismicity can be performed with DAS recordings with results comparable to manual analysis of geophone recordings (i.e. maximum horizontal error on event location around 70 m for both geophones and DAS data) ; and (5) DAS data pre-conditioning (e.g., temporal sub-sampling and channel-stacking) and dedicated processing techniques are strictly necessary for making any real-time monitoring procedure feasible and trustable.


2021 ◽  
Author(s):  
Irene DeFelipe ◽  
Puy Ayarza ◽  
Imma Palomeras ◽  
Juvenal Andrés ◽  
Mario Ruiz ◽  
...  

<p>The Iberian Central System represents an outstanding topographic feature in the central Iberian Peninsula. It is an intraplate mountain range formed by igneous and metasedimentary rocks of the Variscan Iberian Massif that has been exhumed since the Eocene in the context of the Alpine orogeny. The Iberian Central System has been conventionally interpreted as a thick-skinned pop-up mountain range thrust over the Duero and Tajo foreland basins. However, its lithospheric structure and the P-wave velocity distribution are not yet fully resolved. In order to place geophysical constraints on this relevant topographic feature, to identify lithospheric discontinuities, and to unravel the crustal deformation mechanisms, a wide-angle seismic reflection and refraction experiment, CIMDEF (Central Iberian Mechanism of DEFormation), was acquired in 2017 and 2019. It is a NNW-SSE oriented 360-km long profile that runs through the Duero basin, the Iberian Central System and the Tajo basin. First results based on forward modeling by raytracing show an irregularly layered lithosphere and allow to infer the depth extent of the northern Iberian Central System batholith. The crust is ~ 31 km thick under the Duero and Tajo basins and thickens to ~ 39 km under the Iberian Central System. A conspicuous thinning of the lower crust towards the south of the Iberian Central System is also modeled. Along this transect, a continuous and high amplitude upper mantle feature is observed and modeled as the reflection of an interface dipping from 58 to 62 km depth featuring a P-wave velocity contrast of 8.2 to 8.3 km/s. Our preliminary results complement previous models based on global-phase seismic and noise interferometry and gravity data, provide new constraints to validate the accuracy of passive seismic methods at lithospheric scale, and contribute with a resolute P-wave velocity model of the study area to unravel the effect of the Alpine reactivation on the central Iberian Massif.<br>This project has been funded by the EIT-RawMaterials 17024 (SIT4ME) and the MINECO projects: CGL2016-81964-REDE, CGL2014-56548-P.</p>


1974 ◽  
Vol 64 (6) ◽  
pp. 1887-1899
Author(s):  
George A. McMechan ◽  
Warren G. Workman

abstract The observed behavior of P-wave relative amplitudes, as a function of epicentral distance, between 10° and 35°, is controlled primarily by the velocity-depth structure of the upper mantle. P-wave synthetic seismograms calculated by the new quantized ray theory technique are used to determine theoretical log (A/T) versus log Δ curves from a number of upper mantle models. Maximum amplitude arrivals show less model dependence than the first arrivals in the same wave trains, and hence are more consistent magnitude indicators for regions where the upper mantle structure is poorly known. Log (A/T) versus log Δ curves vary considerably, but predictably, from model to model. This model-dependent variation can account for a major part of the large standard deviations usually associated with the calculation of magnitudes from body waves.


2010 ◽  
pp. 29-59
Author(s):  
Philippe Gatignol ◽  
Catherine Potel ◽  
Jean Franois de Belleval
Keyword(s):  

2020 ◽  
Vol 221 (3) ◽  
pp. 1765-1776 ◽  
Author(s):  
Jia Wei ◽  
Li-Yun Fu ◽  
Zhi-Wei Wang ◽  
Jing Ba ◽  
José M Carcione

SUMMARY The Lord–Shulman thermoelasticity theory combined with Biot equations of poroelasticity, describes wave dissipation due to fluid and heat flow. This theory avoids an unphysical behaviour of the thermoelastic waves present in the classical theory based on a parabolic heat equation, that is infinite velocity. A plane-wave analysis predicts four propagation modes: the classical P and S waves and two slow waves, namely, the Biot and thermal modes. We obtain the frequency-domain Green's function in homogeneous media as the displacements-temperature solution of the thermo-poroelasticity equations. The numerical examples validate the presence of the wave modes predicted by the plane-wave analysis. The S wave is not affected by heat diffusion, whereas the P wave shows an anelastic behaviour, and the slow modes present a diffusive behaviour depending on the viscosity, frequency and thermoelasticity properties. In heterogeneous media, the P wave undergoes mesoscopic attenuation through energy conversion to the slow modes. The Green's function is useful to study the physics in thermoelastic media and test numerical algorithms.


Geophysics ◽  
1991 ◽  
Vol 56 (5) ◽  
pp. 664-674 ◽  
Author(s):  
F. Kormendi ◽  
M. Dietrich

We present a method for determining the elastic parameters of a horizontally stratified medium from its plane‐wave reflectivity. The nonlinear inverse problem is iteratively solved by using a generalized least‐squares formalism. The proposed method uses the (relatively) fast convergence properties of the conjugate gradient algorithm and achieves computational efficiency through analytical solutions for calculating the reference and perturbational wavefields. The solution method is implemented in the frequency‐wave slowness domain and can be readily adapted to various source‐receiver configurations. The behavior of the algorithm conforms to the predictions of generalized least‐squares inverse theory: the inversion scheme yields satisfactory results as long as the correct velocity trends are introduced in the starting model. In practice, the inversion algorithm should be applied first in the precritical region because of the strong nonlinear behavior of postcritical data with respect to velocity perturbations. The suggested inversion strategy consists of first inverting for the density and P‐wave velocity (or P‐wave impedance) by considering plane waves in the low slowness region (near‐normal angles of incidence), then in optimizing for the S‐wave velocity by progressively including contributions from the high slowness region (steep angles of incidence). Numerical experiments performed with noise‐free synthetic data prove that the proposed inversion method satisfactorialy reconstructs the elastic properties of a stratified medium from a limited set of plane‐wave components, at a reasonable computing cost.


Geophysics ◽  
1989 ◽  
Vol 54 (10) ◽  
pp. 1339-1343 ◽  
Author(s):  
S. C. Singh ◽  
G. F. West ◽  
C. H. Chapman

The delay‐time (τ‐p) parameterization, which is also known as the plane‐wave decomposition (PWD) of seismic data, has several advantages over the more traditional time‐distance (t‐x) representation (Schultz and Claerbout, 1978). Plane‐wave seismograms in the (τ, p) domain can be used for obtaining subsurface elastic properties (P‐wave and S‐wave velocities and density as functions of depth) from inversion of the observed oblique‐incidence seismic data (e.g., Yagle and Levy, 1985; Carazzone, 1986; Carrion, 1986; Singh et al., 1989). Treitel et al. (1982) performed time migration of plane‐wave seismograms. Diebold and Stoffa (1981) used plane‐wave seismograms to derive a velocity‐depth function. Decomposing seismic data also allows more rapid modeling, since it is faster to compute synthetic seismograms in the (τ, p) than in the (t, x) domain. Unfortunately, the transformation of seismic data from the (t, x) to the (τ, p) domain may produce artifacts, such as those caused by discrete sampling, of the data in space.


Sign in / Sign up

Export Citation Format

Share Document