Shallow near-surface effects

Geophysics ◽  
2016 ◽  
Vol 81 (5) ◽  
pp. T221-T231 ◽  
Author(s):  
Christine E. Krohn ◽  
Thomas J. Murray
Keyword(s):  
Time Delay ◽  
Velocity Gradient ◽  
P Wave ◽  
Unconsolidated Sediments ◽  
P Waves ◽  
Large Computer ◽  
Near Surface ◽  
S Waves ◽  
High Attenuation ◽  
Pulse Broadening

The top 6 m of the near surface has a surprisingly large effect on the behavior of P- and S-waves. For unconsolidated sediments, the P-wave velocity gradient and attenuation can be quite large. Computer modeling should include these properties to accurately reproduce seismic effects of the near surface. We have used reverse VSP data and computer simulations to demonstrate the following effects for upgoing P-waves. Near the surface, we have observed a large time delay, indicating low velocity ([Formula: see text]), and considerable pulse broadening, indicating high attenuation ([Formula: see text]). Consequently, shallowly buried geophones have greater high-frequency bandwidth compared with surface geophones. In addition, there is a large velocity gradient in the shallow near surface (factor of 10 in 5 m), resulting in the rotation of P-waves to the vertical with progressively smaller amplitudes recorded on horizontal phones. Finally, we have found little indication of a reflection or ghost from the surface, although downgoing reflections have been observed from interfaces within the near surface. In comparison, the following have been observed for upgoing S-waves: There is a small increase in the time delay or pulse broadening near the surface, indicating a smaller velocity gradient and less change in attenuation. In addition, the surface reflection coefficient is nearly one with a prominent surface ghost.

Seismic vertical transversely isotropic parameter inversion from P- and S-wave cross-borehole measurements in an aquifer environment

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. D101-D116
Author(s):  
Julius K. von Ketelhodt ◽  
Musa S. D. Manzi ◽  
Raymond J. Durrheim ◽  
Thomas Fechner
Keyword(s):  
Transversely Isotropic ◽  
P Wave ◽  
Perturbation Equation ◽  
S Wave ◽  
P Waves ◽  
Data Set ◽  
Near Surface ◽  
S Waves ◽  
Wave Measurements ◽  
Wave Splitting

Joint P- and S-wave measurements for tomographic cross-borehole analysis can offer more reliable interpretational insight concerning lithologic and geotechnical parameter variations compared with P-wave measurements on their own. However, anisotropy can have a large influence on S-wave measurements, with the S-wave splitting into two modes. We have developed an inversion for parameters of transversely isotropic with a vertical symmetry axis (VTI) media. Our inversion is based on the traveltime perturbation equation, using cross-gradient constraints to ensure structural similarity for the resulting VTI parameters. We first determine the inversion on a synthetic data set consisting of P-waves and vertically and horizontally polarized S-waves. Subsequently, we evaluate inversion results for a data set comprising jointly measured P-waves and vertically and horizontally polarized S-waves that were acquired in a near-surface ([Formula: see text]) aquifer environment (the Safira research site, Germany). The inverted models indicate that the anisotropy parameters [Formula: see text] and [Formula: see text] are close to zero, with no P-wave anisotropy present. A high [Formula: see text] ratio of up to nine causes considerable SV-wave anisotropy despite the low magnitudes for [Formula: see text] and [Formula: see text]. The SH-wave anisotropy parameter [Formula: see text] is estimated to be between 0.05 and 0.15 in the clay and lignite seams. The S-wave splitting is confirmed by polarization analysis prior to the inversion. The results suggest that S-wave anisotropy may be more severe than P-wave anisotropy in near-surface environments and should be taken into account when interpreting cross-borehole S-wave data.

Simulation of 3D elastic wave propagation in the Salt Lake Basin

1995 ◽  
Vol 85 (6) ◽  
pp. 1688-1710 ◽  
Author(s):  
Kim B. Olsen ◽  
James C. Pechmann ◽  
Gerard T. Schuster
Keyword(s):  
Salt Lake ◽  
Incidence Angle ◽  
P Wave ◽  
Surface Velocity ◽  
Unconsolidated Sediments ◽  
Lake Basin ◽  
P Waves ◽  
Near Surface ◽  
Rock Site ◽  
Particle Velocities

Abstract We have used a 3D finite-difference method to model 0.2 to 1.2 Hz elastodynamic site amplification in the Salt Lake Valley, Utah. The valley is underlain by a sedimentary basin, which in our model has dimensions of 48 by 25 by 1.3 km. Simulations are carried out for a P wave propagating vertically from below and for P waves propagating horizontally to the north, south, east, and west in a two-layer model consisting of semi-consolidated sediments surrounded by bedrock. Results show that in general, sites with the largest particle velocities, cumulative kinetic energies, duration times of motion, and spectral magnitudes overlie the deepest parts of the basin. The maximum values of these parameters are generally found above steeply dipping parts of the basin walls. The largest vector particle velocities are associated with P or SV waves that come from within 10° of the source azimuth. Low-energy S and surface waves follow the strongest arrivals. The largest peak particle velocities, cumulative kinetic energies, signal durations, and spectral magnitudes in the simulations are, respectively, 2.9, 15.9, 40.0, and 3.5 times greater than the values at a rock site measured on the component parallel to the propagation direction of the incident P wave. Scattering and/or mode conversions at the basin boundaries contribute significantly to the signal duration times. As a check on the validity of our simulations, we compared our 3D synthetic seismograms for the vertically incident plane P wave to seismograms of nearly vertically incident teleseismic P waves recorded at an alluvium site in the valley and at a nearby rock site. The 3D synthetics for the alluvium site overestimate the relatively small amplification of the initial P wave and underestimate the large amplification of the coda. Using 2D simulations, we find that most of the discrepancies between the 3D synthetic and observed records can be explained by an apparently incorrect total sediment thickness, omission from the model of the near-surface low-velocity unconsolidated sediments and of attenuation, and the inexact modeling of the incidence angle of the teleseism. The records from a 2D simulation in which these deficiencies are remedied (with Q = 65), and which also includes topography and a near-surface velocity gradient in the bedrock, provide a better match to the teleseismic data than the records from the simple two-layer 3D simulation. Our results suggest that for steeply incident P waves, the impedance decrease and resonance effects associated with the deeper basin structure control the amplification of the initial P-wave arrival, whereas reverberations in the near-surface unconsolidated sediments generate the large-amplitude coda. These reverberations are caused mainly by P-to-S converted waves, and their strength is therefore highly sensitive to the incidence angle of the source.

A simple method for resolving large converted‐wave (P-SV) statics

Geophysics ◽  
1993 ◽  
Vol 58 (3) ◽  
pp. 429-433 ◽  
Author(s):  
Peter W. Cary ◽  
David W. S. Eaton
Keyword(s):  
P Wave ◽  
S Wave ◽  
Converted Wave ◽  
P Waves ◽  
Simple Method ◽  
Near Surface ◽  
S Waves ◽  
Static Solutions ◽  
Surface Fluctuations

The processing of converted‐wave (P-SV) seismic data requires certain special considerations, such as commonconversion‐point (CCP) binning techniques (Tessmer and Behle, 1988) and a modified normal moveout formula (Slotboom, 1990), that makes it different for processing conventional P-P data. However, from the processor’s perspective, the most problematic step is often the determination of residual S‐wave statics, which are commonly two to ten times greater than the P‐wave statics for the same location (Tatham and McCormack, 1991). Conventional residualstatics algorithms often produce numerous cycle skips when attempting to resolve very large statics. Unlike P‐waves, the velocity of S‐waves is virtually unaffected by near‐surface fluctuations in the water table (Figure 1). Hence, the P‐wave and S‐wave static solutions are largely unrelated to each other, so it is generally not feasible to approximate the S‐wave statics by simply scaling the known P‐wave static values (Anno, 1986).

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.

Crosswell seismic imaging in a contaminated basalt aquifer

Geophysics ◽  
2004 ◽  
Vol 69 (1) ◽  
pp. 16-24 ◽  
Author(s):  
Thomas M. Daley ◽  
Ernest L. Majer ◽  
John E. Peterson
Keyword(s):  
Seismic Source ◽  
P Wave ◽  
Low Velocity ◽  
S Waves ◽  
High Attenuation ◽  
Well Spacing ◽  
Simultaneous Acquisition ◽  
Environmental Laboratory ◽  
New Type ◽  
Borehole Seismic

Multiple seismic crosswell surveys have been acquired and analyzed in a fractured basalt aquifer at Idaho National Engineering and Environmental Laboratory. Most of these surveys used a high‐frequency (1000–10,000 Hz) piezoelectric seismic source to obtain P‐wave velocity tomograms. The P‐wave velocities range from less than 3200 m/s to more than 5000 m/s. Additionally, a new type of borehole seismic source was deployed as part of the subsurface characterization program at this contaminated groundwater site. This source, known as an orbital vibrator, allows simultaneous acquisition of P‐ and S‐waves at frequencies of 100 to 400 Hz, and acquisition over larger distances. The velocity tomograms show a relationship to contaminant transport in the groundwater; zones of high contaminant concentration are coincident with zones of low velocity and high attenuation and are interpreted to be fracture zones at the boundaries between basalt flows. The orbital vibrator data show high Vp/Vs values, from 1.8 to 2.8. In spite of the lower resolution of orbital vibrator data, these data were sufficient for constraining hydrologic models at this site while achieving imaging over large interwell distances. The combination of piezoelectric data for closer well spacing and orbital vibrator data for larger well spacings has provided optimal imaging capability and has been instrumental in our understanding of the site aquifer's hydrologic properties and its scale of heterogeneity.

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.

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.

scholarly journals MORE ABOUT THE WILSONIAN ANALYSIS ON THE PIONLESS NEFT

2009 ◽  
Vol 24 (16n17) ◽  
pp. 3191-3225 ◽  
Author(s):  
KOJI HARADA ◽  
HIROFUMI KUBO ◽  
ATSUSHI NINOMIYA
Keyword(s):  
Fixed Points ◽  
Scattering Amplitude ◽  
Flow Equation ◽  
Effective Field ◽  
P Wave ◽  
P Waves ◽  
S Waves ◽  
Partial Waves ◽  
The One ◽  
Power Divergence

We extend our Wilsonian renormalization group (RG) analysis on the pionless nuclear effective field theory in the two-nucleon sector in two ways; on the one hand, (1) we enlarge the space of operators up to including those of [Formula: see text] in the S waves, and, on the other hand, (2) we consider the RG flows in higher partial waves (P and D waves). In the larger space calculations, we find, in addition to nontrivial fixed points, two "fixed lines" and a "fixed surface" which are related to marginal operators. In the higher partial wave calculations, we find similar phase structures to that of the S waves, but there are two relevant directions in the P waves at the nontrivial fixed points and three in the D waves. We explain the physical meaning of the P-wave phase structure by explicitly calculating the low-energy scattering amplitude. We also discuss the relation between the Legendre flow equation which we employ and the RG equation by Birse, McGovern and Richardson, and possible implementation of power divergence subtraction in higher partial waves.

Imaging beneath a high‐velocity layer using converted waves

Geophysics ◽  
1992 ◽  
Vol 57 (11) ◽  
pp. 1444-1452 ◽  
Author(s):  
Guy W. Purnell
Keyword(s):  
High Velocity ◽  
Energy Partitioning ◽  
P Wave ◽  
Converted Wave ◽  
P Waves ◽  
Lower Velocity ◽  
S Waves ◽  
High Velocity Layer ◽  
The Waves ◽  
Velocity Layer

High‐velocity layers (HVLs) often hinder seismic imaging of deeper reflectors using conventional techniques. A major factor is often the unusual energy partitioning of waves incident at an HVL boundary from lower‐velocity material. Using elastic physical modeling, I demonstrate that one effect of this factor is to limit the range of dips beneath an HVL that can be imaged using unconverted P‐wave arrivals. At the same time, however, partitioning may also result in P‐waves outside the HVL coupling efficiently with S‐waves inside. By exploiting some of the waves that convert upon transmission into and/or out of the physical‐model HVL, I am able to image a much broader range of underlying dips. This is accomplished by acoustic migration tailored (via the migration velocities used) for selected families of converted‐wave arrivals.

Rapid map and inversion of P-SV waves

Geophysics ◽  
1991 ◽  
Vol 56 (6) ◽  
pp. 859-862 ◽  
Author(s):  
Robert R. Stewart
Keyword(s):  
Spatial Resolution ◽  
Seismic Data ◽  
P Wave ◽  
Rock Properties ◽  
Low Velocity ◽  
P Waves ◽  
Near Surface ◽  
Seismic Profiling ◽  
Vertical Seismic Profiling ◽  
And Inversion

Multicomponent seismic recordings are currently being analyzed in an attempt to improve conventional P‐wave sections and to find and use rock properties associated with shear waves (e.g. Dohr, 1985; Danbom and Dominico, 1986). Mode‐converted (P-SV) waves hold a special interest for several reasons: They are generated by conventional P‐wave sources and have only a one‐way travel path as a shear wave through the typically low velocity and attenuative near surface. For a given frequency, they will have a shorter wavelength than the original P wave, and thus offer higher spatial resolution; this has been observed in several vertical seismic profiling (VSP) cases (e.g., Geis et al., 1990). However, for surface seismic data, converted waves are often found to be of lower frequency than P-P waves (e.g., Eaton et al., 1991).

Sign in / Sign up

Export Citation Format

Share Document