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.

Separation and enhancement of split S‐waves on multicomponent shot records from the BIRPS WISPA experiment

Geophysics ◽  
1994 ◽  
Vol 59 (1) ◽  
pp. 131-139 ◽  
Author(s):  
M. Boulfoul ◽  
D. R. Watts
Keyword(s):  
Synthetic Data ◽  
P Wave ◽  
Moving Window ◽  
Wave Analysis ◽  
S Wave ◽  
Near Surface ◽  
S Waves ◽  
Wave Velocities ◽  
Wave Splitting ◽  
Explosive Source

Instantaneous rotations are combined with f-k filtering to extract coherent S‐wave events from multicomponent shot records recorded by British Institutions Reflection Profiling Syndicate (BIRPS) Weardale Integrated S‐wave and P‐wave analysis (WISPA) experiment. This experiment was an attempt to measure the Poisson’s ratio of the lower crest by measuring P‐wave and S‐wave velocities. The multihole explosive source technique did generate S‐waves although not of opposite polarization. Attempts to produce stacks of the S‐wave data are unsuccessful because S‐wave splitting in the near surface produced random polarizations from receiver group to receiver group. The delay between the split wavelets varies but is commonly between 20 to 40 ms for 10 Hz wavelets. Dix hyperbola are produced on shot records after instantaneous rotations are followed by f-k filtering. To extract the instantaneous polarization, the traces are shifted back by the length of a moving window over which the calculation is performed. The instantaneous polarization direction is computed from the shifted data using the maximum eigenvector of the covariance matrix over the computation window. Split S‐waves are separated by the instantaneous rotation of the unshifted traces to the directions of the maximum eigenvectors determined for each position of the moving window. F-K filtering is required because of the presence of mode converted S‐waves and S‐waves produced by the explosive source near the time of detonation. Examples from synthetic data show that the method of instantaneous rotations will completely separate split S‐waves if the length of the moving window over which the calculation is performed is the length of the combined split wavelets. Separation may be achieved on synthetic data for wavelet delays as small as two sample intervals.

Shallow velocity structure and poisson's ratio at the Tarzana, California, strong-motion accelerometer site

1996 ◽  
Vol 86 (6) ◽  
pp. 1704-1713 ◽  
Author(s):  
R. D. Catchings ◽  
W. H. K. Lee
Keyword(s):  
Urban Areas ◽  
Velocity Structure ◽  
Strong Motion ◽  
P Wave ◽  
S Wave ◽  
Near Surface ◽  
Velocity Models ◽  
Wave Velocities ◽  
Poisson's Ratios ◽  
Poisson’S Ratios

Abstract The 17 January 1994, Northridge, California, earthquake produced strong ground shaking at the Cedar Hills Nursery (referred to here as the Tarzana site) within the city of Tarzana, California, approximately 6 km from the epicenter of the mainshock. Although the Tarzana site is on a hill and is a rock site, accelerations of approximately 1.78 g horizontally and 1.2 g vertically at the Tarzana site are among the highest ever instrumentally recorded for an earthquake. To investigate possible site effects at the Tarzana site, we used explosive-source seismic refraction data to determine the shallow (<70 m) P-and S-wave velocity structure. Our seismic velocity models for the Tarzana site indicate that the local velocity structure may have contributed significantly to the observed shaking. P-wave velocities range from 0.9 to 1.65 km/sec, and S-wave velocities range from 0.20 and 0.6 km/sec for the upper 70 m. We also found evidence for a local S-wave low-velocity zone (LVZ) beneath the top of the hill. The LVZ underlies a CDMG strong-motion recording site at depths between 25 and 60 m below ground surface (BGS). Our velocity model is consistent with the near-surface (<30 m) P- and S-wave velocities and Poisson's ratios measured in a nearby (<30 m) borehole. High Poisson's ratios (0.477 to 0.494) and S-wave attenuation within the LVZ suggest that the LVZ may be composed of highly saturated shales of the Modelo Formation. Because the lateral dimensions of the LVZ approximately correspond to the areas of strongest shaking, we suggest that the highly saturated zone may have contributed to localized strong shaking. Rock sites are generally considered to be ideal locations for site response in urban areas; however, localized, highly saturated rock sites may be a hazard in urban areas that requires further investigation.

Comparison of P‐ and S‐wave velocities and Q’S from VSP and sonic log data

Geophysics ◽  
1994 ◽  
Vol 59 (10) ◽  
pp. 1512-1529 ◽  
Author(s):  
Gopa S. De ◽  
Donald F. Winterstein ◽  
Mark A. Meadows
Keyword(s):  
Velocity Dispersion ◽  
P Wave ◽  
S Wave ◽  
S Waves ◽  
Wave Velocities ◽  
Sonic Log ◽  
Transit Times ◽  
Wave Slope ◽  
Northern Alberta ◽  
Q Values

We compared P‐ and S‐wave velocities and quality factors (Q’S) from vertical seismic profiling (VSP) and sonic log measurements in five wells, three from the southwest San Joaquin Basin of California, one from near Laredo, Texas, and one from northern Alberta. Our purpose was to investigate the bias between sonic log and VSP velocities and to examine to what degree this bias might be a consequence of dispersion. VSPs and sonic logs were recorded in the same well in every case. Subsurface formations were predominantly clastic. The bias found was that VSP transit times were greater than sonic log times, consistent with normal dispersion. For the San Joaquin wells, differences in S‐wave transit times averaged 1–2 percent, while differences in P‐wave transit times averaged 6–7 percent. For the Alberta well, the situation was reversed, with differences in S‐wave transit times being about 6 percent, while those for P‐waves were 2.5 percent. For the Texas well, the differences averaged about 4 percent for both P‐ and S‐waves. Drift‐curve slopes for S‐waves tended to be low where the P‐wave slopes were high and vice versa. S‐wave drift‐curve slopes in the shallow California wells were 5–10 μs/ft (16–33 μs/m) and the P‐wave slopes were 15–30 μs/ft (49–98 μs/m). The S‐wave slope in sandstones in the northern Alberta well was up to 50 μs/ft (164 μs/m), while the P‐wave slope was about 5 μs/ft (16 μs/m). In the northern Alberta well the slopes for both P‐ and S‐waves flattened in the carbonate. In the Texas well, both P‐ and S‐wave drifts were comparable. We calculated (Q’s) from a velocity dispersion formula and from spectral ratios. When the two Q’s agreed, we concluded that velocity dispersion resulted solely from absorption. These Q estimation methods were reliable only for Q values smaller than 20. We found that, even with data of generally outstanding quality, Q values determined by standard methods can have large uncertainties, and negative Q’s may be common.

Practical concept of traveltime inversion of simulated P-wave vertical seismic profile data in weak to moderate arbitrary anisotropy

Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. C107-C123
Author(s):  
Ivan Pšenčík ◽  
Bohuslav Růžek ◽  
Petr Jílek
Keyword(s):  
Waveform Inversion ◽  
Anisotropic Media ◽  
Compressional Wave ◽  
Seismic Profile ◽  
P Wave ◽  
Source Distribution ◽  
S Wave ◽  
Vertical Seismic Profile ◽  
Weak Anisotropy ◽  
Traveltime Inversion

We have developed a practical concept of compressional wave (P-wave) traveltime inversion in weakly to moderately anisotropic media of arbitrary symmetry and orientation. The concept provides sufficient freedom to explain and reproduce observed anisotropic seismic signatures to a high degree of accuracy. The key to this concept is the proposed P-wave anisotropy parameterization (A-parameters) that, together with the use of the weak-anisotropy approximation, leads to a significantly simplified theory. Here, as an example, we use a simple and transparent formula relating P-wave traveltimes to 15 P-wave A-parameters describing anisotropy of arbitrary symmetry. The formula is used in the inversion scheme, which does not require any a priori information about anisotropy symmetry and its orientation, and it is applicable to weak and moderate anisotropy. As the first step, we test applicability of the proposed scheme on a blind inversion of synthetic P-wave traveltimes generated in vertical seismic profile experiments in homogeneous models. Three models of varying anisotropy are used: tilted orthorhombic and triclinic models of moderate anisotropy (approximately 10%) and an orthorhombic model of strong anisotropy (>25%) with a horizontal plane of symmetry. In all cases, the inversion yields the complete set of 15 P-wave A-parameters, which make reconstruction of corresponding phase-velocity surfaces possible with high accuracy. The inversion scheme is robust with respect to noise and the source distribution pattern. Its quality depends on the angular illumination of the medium; we determine how the absence of nearly horizontal propagation directions affects inversion accuracy. The results of the inversion are applicable, for example, in migration or as a starting model for inversion methods, such as full-waveform inversion, if a model refinement is desired. A similar procedure could be designed for the inversion of S-wave traveltimes in anisotropic media of arbitrary symmetry.

Amplitude balancing in separating P- and S-waves in 2D and 3D elastic seismic data

Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. S103-S113 ◽  
Author(s):  
Robert Sun ◽  
George A. McMechan ◽  
Han-Hsiang Chuang
Keyword(s):  
Wave Amplitude ◽  
Amplitude Ratio ◽  
Velocity Ratio ◽  
Displacement Component ◽  
P Wave ◽  
Elastic Displacement ◽  
Reverse Time ◽  
S Wave ◽  
Near Surface ◽  
S Waves

The reflected P- and S-waves in elastic displacement component data recorded at the earth’s surface are separated by reverse-time (downward) extrapolation of the data in an elastic computational model, followed by calculations to give divergence (dilatation) and curl (rotation) at a selected reference depth. The surface data are then reconstructed by separate forward-time (upward) scalar extrapolations, from the reference depth, of the magnitude of the divergence and curl wavefields, and extraction of the separated P- and S-waves, respectively, at the top of the models. A P-wave amplitude will change by a factor that is inversely proportional to the P-velocity when it is transformed from displacement to divergence, and an S-wave amplitude will change by a factor that is inversely proportional to the S-velocity when it is transformed from displacement to curl. Consequently, the ratio of the P- to the S-wave amplitude (the P-S amplitude ratio) in the form of divergence and curl (postseparation) is different from that in the (preseparation) displacement form. This distortion can be eliminated by multiplying the separated S-wave (curl) by a relative balancing factor (which is the S- to P-velocity ratio); thus, the postseparation P-S amplitude ratio can be returned to that in the preseparation data. The absolute P- and S-wave amplitudes are also recoverable by multiplying them by a factor that depends on frequency, on the P-velocity α, and on the unit of α and is location-dependent if the near-surface P-velocity is not constant.

EXPERIMENTAL EVALUATION OF ULTRASONIC VELOCITIES AND ANISOTROPY IN THE TUSCALOOSA MARINE SHALE FORMATION

2020 ◽  
pp. 1-62 ◽  
Author(s):  
Jamal Ahmadov ◽  
Mehdi Mokhtari
Keyword(s):  
Wave Velocity ◽  
Mineralogical Composition ◽  
Organic Content ◽  
P Wave ◽  
S Wave ◽  
Velocity Anisotropy ◽  
Wave Velocities ◽  
Marine Shale ◽  
P Wave Velocity ◽  
Degree Of Anisotropy

Tuscaloosa Marine Shale (TMS) formation is a clay- and organic-rich emerging shale play with a considerable amount of hydrocarbon resources. Despite the substantial potential, there have been only a few wells drilled and produced in the formation over the recent years. The analyzed TMS samples contain an average of 50 wt% total clay, 27 wt% quartz and 14 wt% calcite and the mineralogy varies considerably over the small intervals. The high amount of clay leads to pronounced anisotropy and the frequent changes in mineralogy result in the heterogeneity of the formation. We studied the compressional (VP) and shear-wave (VS) velocities to evaluate the degree of anisotropy and heterogeneity, which impact hydraulic fracture growth, borehole instabilities, and subsurface imaging. The ultrasonic measurements of P- and S-wave velocities from five TMS wells are the best fit to the linear relationship with R2 = 0.84 in the least-squares criteria. We observed that TMS S-wave velocities are relatively lower when compared to the established velocity relationships. Most of the velocity data in bedding-normal direction lie outside constant VP/VS lines of 1.6–1.8, a region typical of most organic-rich shale plays. For all of the studied TMS samples, the S-wave velocity anisotropy exhibits higher values than P-wave velocity anisotropy. In the samples in which the composition is dominated by either calcite or quartz minerals, mineralogy controls the velocities and VP/VS ratios to a great extent. Additionally, the organic content and maturity account for the velocity behavior in the samples in which the mineralogical composition fails to do so. The results provide further insights into TMS Formation evaluation and contribute to a better understanding of the heterogeneity and anisotropy of the play.

Shear-wave anisotropy in the Parkfield Varian well VSP

1990 ◽  
Vol 80 (4) ◽  
pp. 857-869 ◽  
Author(s):  
T. M. Daley ◽  
T. V. McEvilly
Keyword(s):  
Shear Wave ◽  
San Andreas Fault ◽  
Monitoring Program ◽  
Seismic Profile ◽  
S Wave ◽  
S Waves ◽  
Velocity Anisotropy ◽  
San Andreas ◽  
The Difference ◽  
Fault Trace

Abstract A vertical seismic profile (VSP) survey was run to 1334 m depth in the instrumented Varian well, 1.4 km from the San Andreas fault trace at Parkfield, California, to test the sensor string shortly after its permanent installation. The cable subsequently failed near the 1000 m level, so the test survey represents the deepest data acquired in the study. A shear-wave vibrator source was used at three ofsets and two orthogonal orientations, and the data have been processed for P- and S-wave velocities and for S-wave velocity anisotropy. Velocities are well-determined (3.3 and 1.9 km/sec, respectively, at the deeper levels), and the S waves are seen clearly to be split by anisotropy below about 400 m. Some 8 per cent velocity difference is apparent between polarizations parallel to and perpendicular to the San Andreas fault (faster and slower, respectively), and the difference seems to decrease with distance from the fault, suggesting that the cause may be the fabric of the fault zone. Repeated surveys at the 1000 m depth are being conducted as part of the Parkfield monitoring program.

P- and S-wave characterization of near‐surface reflectivity from glacial tills using vertical seismic profiles

Geophysics ◽  
1999 ◽  
Vol 64 (3) ◽  
pp. 970-980 ◽  
Author(s):  
Bradley J. Carr ◽  
Zoltan Hajnal
Keyword(s):  
Wave Energy ◽  
Acoustic Property ◽  
Seismic Profile ◽  
P Wave ◽  
S Wave ◽  
Glacial Deposits ◽  
Near Surface ◽  
Borehole Geophysics ◽  
Vsp Data

Fundamental reflectivity properties are established within the glacial deposits of central Saskatchewan, Canada. Multicomponent vertical seismic profile (VSP) data collected in three shallow boreholes are used to obtain detailed acoustic property information within the first 80 m of the near‐surface strata. The integration of both P- and S-wave VSP data, in conjunction with other borehole geophysics, provided a unique opportunity to obtain in‐situ seismic reflection response properties in layered clay and sand tills. P- and S-wave interval velocity profiles, in conjunction with P- and S-wave VSP reflectivities are analyzed to provide insight into seismic wavefield behavior within ∼80 m of the surface. In general, shear wave energy identifies more reflective intervals than the P-wave energy because of better vertical resolution for S-wave energy (0.75 m) compared to P-wave energy (2.3 m) based on quarter wavelength criterion. For these saturated, unconsolidated glacial deposits, more details about the lithologic constituents and in‐situ porosity are detectable from the S-wave reflectivity, but P-wave reflections provide a good technique for mapping the bulk changes. The principal cause of seismic reflectivity is the presence and/or amount of sand, and the degree of fluid‐filled porosity within the investigated formations.

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.

Seismological studies in a rock body at Chalk River, Ontario and their relation to fractures

1982 ◽  
Vol 19 (8) ◽  
pp. 1535-1547 ◽  
Author(s):  
C. Wright
Keyword(s):  
Wave Velocity ◽  
Test Site ◽  
P Wave ◽  
Seismic Velocities ◽  
S Wave ◽  
P Waves ◽  
Rock Body ◽  
Near Surface ◽  
Wave Velocities ◽  
P Wave Velocity

Seismological experiments have been undertaken at a test site near Chalk River, Ontario that consists of crystalline rocks covered by glacial sediments. Near-surface P and S wave velocity and amplitude variations have been measured along profiles less than 2 km in length. The P and S wave velocities were generally in the range 4.5–5.6 and 2.9–3.2 km/s, respectively. These results are consistent with propagation through fractured gneiss and monzonite, which form the bulk of the rock body. The P wave velocity falls below 5.0 km/s in a region where there is a major fault and in an area of high electrical conductivity; such velocity minima are therefore associated with fracture systems. For some paths, the P and 5 wave velocities were in the ranges 6.2–6.6 and 3.7–4.1 km/s, respectively, showing the presence of thin sheets of gabbro. Temporal changes in P travel times of up to 1.4% over a 12 h period were observed where the sediment cover was thickest. The cause may be changes in the water table. The absence of polarized SH arrivals from specially designed shear wave sources indicates the inhomogeneity of the test site. A Q value of 243 ± 53 for P waves was derived over one relatively homogeneous profile of about 600 m length. P wave velocity minima measured between depths of 25 and 250 m in a borehole correlate well with the distribution of fractures inferred from optical examination of borehole cores, laboratory measurements of seismic velocities, and tube wave studies.

Sign in / Sign up

Export Citation Format

Share Document