Multicomponent VSP imaging of tight-gas sands

Geophysics ◽  
2006 ◽  
Vol 71 (6) ◽  
pp. E83-E90 ◽  
Author(s):  
John O’Brien ◽  
Ron Harris
Keyword(s):  
Shear Wave ◽  
Negative Impact ◽  
Mode Conversion ◽  
Strong Support ◽  
Seismic Profile ◽  
Limited Range ◽  
P Wave ◽  
Wave Source ◽  
Cotton Valley ◽  
Vsp Data

Low-porosity Bossier and Cotton Valley sands of the East Texas Basin, U. S., have only a small acoustic impedance contrast with the encasing shales but a greater relative contrast in shear-wave impedance. Vertical seismic profile (VSP) data acquired with both a near-offset and far-offset P-wave source clearly demonstrate the P-P reflectivity and P-S mode conversions within the Bossier section. We designate conventional P-wave reflectivity as P-P, shear-wave reflectivity as S-S, and P-wave/shear-wave mode conversion data as P-S. While Bossier P-P reflectivity is low, it appears to be adequate for mapping thick sandbodies such as the York Sandstone, the main exploration target in this area. However, P-P reflectivity is even lower and is inadequate for imaging the overlying Cotton Valley Sands. In contrast, the far-offset VSP data acquired with a P-wave source demonstrate a high level of P-S-mode conversion, which is used to image this interval with definition that is not provided by P-P reflectivity. This provides strong support for the use of P-S-mode conversion imaging for seismic characterization of tight sand reservoirs. Near-offset shear-wave VSP data acquired with a shearwave source show low S/N ratio and limited bandwidth for the downgoing waveform because of the depth of the target; shear-wave energy appears to have a more limited range of propagation than P-waves. Such effects may also have a strong negative impact on multicomponent imaging of these sands using surface seismic techniques. Multicomponent 3D VSP imaging provides a superior solution by placing the geophones closer to the subsurface zone of interest.

Pure S-waves in land P-wave source VSP data

2007 ◽  
Vol 4 (3) ◽  
pp. 173-182 ◽  
Author(s):  
Liu Yang ◽  
Zhang Qinghong ◽  
Bao Leiying ◽  
Wei Xiucheng
Keyword(s):  
P Wave ◽  
Wave Source ◽  
S Waves ◽  
Vsp Data

Quantifying hydraulically induced fracture height and density from rapid time-lapse DAS VSP

Geophysics ◽  
2020 ◽  
pp. 1-26
Author(s):  
Xiaomin Zhao ◽  
Mark E. Willis ◽  
Tanya Inks ◽  
Glenn A. Wilson
Keyword(s):  
Rock Physics ◽  
Horizontal Wells ◽  
Time Lapse ◽  
Seismic Profile ◽  
P Wave ◽  
Fracture Properties ◽  
Vsp Data ◽  
Rock Physics Modeling ◽  
Physics Modeling ◽  
Fracture Height

Several recent studies have advanced the use of time-lapse distributed acoustic sensing (DAS) vertical seismic profile (VSP) data in horizontal wells for determining hydraulically stimulated fracture properties. Hydraulic fracturing in a horizontal well typically generates vertical fractures in the rock medium around each stage. We model the hydraulically stimulated formation with vertical fracture sets about the lateral wellbore as a horizontally transverse isotropic (HTI) medium. Rock physics modeling is used to relate the anisotropy parameters to fracture properties. This modeling was used to develop an inversion for P-wave time delay to fracture height and density of each stage. Field data from two horizontal wells were analyzed, and fracture height evaluated using this technique agreed with microseismic analysis.

Mode Conversion Technique Employed in Shear Wave Velocity Studies of Rock Samples Under Axial and Uniform Compression

1967 ◽  
Vol 7 (02) ◽  
pp. 136-148 ◽  
Author(s):  
A.R. Gregory
Keyword(s):  
Shear Wave ◽  
Wave Velocity ◽  
Shear Wave Velocity ◽  
Wave Energy ◽  
Pressure Range ◽  
Mode Conversion ◽  
P Wave ◽  
S Wave ◽  
Rock Samples ◽  
Mode Converters

Abstract A shear wave velocity laboratory apparatus and techniques for testing rock samples under simulated subsurface conditions have been developed. In the apparatus, two electromechanical transducers operating in the frequency range 0.5 to 5.0 megahertz (MHz: megacycles per second) are mounted in contact with each end of the sample. Liquid-solid interfaces of Drakeol-aluminum are used as mode converters. In the generator transducer, there is total mode conversion from P-wave energy to plain S-wave energy, S-wave energy is converted back to P-wave energy in the motor transducer. Similar transducers without mode converters are used to measure P-wave velocities. The apparatus is designed for testing rock samples under axial or uniform loading in the pressure range 0 to 12,000 psi. The transducers have certain advantages over those used by King,1 and the measurement techniques are influenced less by subjective elements than other methods previously reported. An electronic counter-timer having a resolution of 10 nanoseconds measures the transit time of ultrasonic pulses through the sample; elastic wave velocities of most homogeneous materials can be measured with errors of less than 1 percent. S- and P-wave velocity measurements on Bandera sandstone and Solenhofen limestone are reported for the axial pressure range 0 to 6,000 psi and for the uniform pressure range 0 to 10,000 psi. The influence of liquid pore saturants on P- and S-wave velocity is investigated and found to be in broad agreement with Biot's theory. In specific areas, the measurements do not conform to theory. Velocities of samples measured under axial and uniform loading are compared and, in general, velocities measured under uniform stress are higher than those measured under axial stress. Liquid pore fluids cause increases in Poisson's ratio and the bulk modulus but reduce the rigidity modulus, Young's modulus and the bulk compressibility. INTRODUCTION Ultrasonic pulse methods for measuring the shear wave velocity of rock samples in the laboratory have been gradually improved during the last few years. Early experimental pulse techniques reported by Hughes et al.2, and by Gregory3 were beset by uncertainties in determining the first arrival of the shear wave (S-wave) energy. Much of this ambiguity was caused by the multiple modes propagated by piezoelectric crystals and by boundary conversions in the rock specimens. Shear wave velocity data obtained from the critical angle method, described by Schneider and Burton4 and used later by King and Fatt5 and by Gregory,3,6 are of limited accuracy, and interpreting results is too complicated for routine laboratory work. The mode conversion method described by Jamieson and Hoskins7 was recently used by King1 for measuring the S-wave velocities of dry and liquid-saturated rock samples. Glass-air interfaces acted as mode converters in the apparatus, and much of the compressional (P-wave) energy apparently was eliminated from the desired pure shear mode. A more detailed discussion of the current status of laboratory pulse methods applied to geological specimens is given in a review by Simmons.8

Near-surface velocities and attenuation at two boreholes near Anza, California, from logging data

1990 ◽  
Vol 80 (4) ◽  
pp. 807-831 ◽  
Author(s):  
Jon B. Fletcher ◽  
Tom Fumal ◽  
Hsi-Ping Liu ◽  
Linda C. Carroll
Keyword(s):  
Shear Wave ◽  
Wave Velocity ◽  
Spectral Ratio ◽  
Body Waves ◽  
P Wave ◽  
Wave Source ◽  
Near Surface ◽  
Wave Velocities ◽  
Rock Structure ◽  
Over The Top

Abstract To investigate near-surface site effects in granite rock, we drilled 300-m-deep boreholes at two sites which are collocated with stations from the digital array at Anza, California. The first borehole was sited at station KNW (Keenwild fire station), which is located along a ridge line about 8.7 km east of the San Jacinto Fault zone. Station PFO (Piñon Flat Observatory), chosen for the second site, is another 6 km further to the east of station KNW and is located on a gently sloping hillside. We logged each borehole for P- and S-wave velocities, as well as for crack density and orientation. P waves were generated by striking a plate with a hammer at the surface. A tool consisting of weighted anvils driven by compressed air against end plates along a 3.5-m beam was used to generate shear waves. Signals were recorded downhole with a three-component sensor package at 2.5-m intervals from the surface to 50 m depth, and at 5-m intervals from 50 m depth to the bottom of the hole. Velocities were determined by differencing the measured arrival times of first arrivals or peaks over each interval in depth. Travel times were computed for the first breaks at shallow depths, however, below about 100 m depth, times were computed for the first peaks rather than for first breaks since the first arrival was no longer clearly distinguishable. The KNW site yielded a shear velocity of 1.9 km/sec by only 30 m in depth and reached close to 2.6 km/sec at the bottom of the hole. P-wave velocities at KNW were also high at 5.4 km/sec starting at 120 m depth. The PFO site had similar but slightly higher shear-wave velocities. The bottom-hole shear-wave velocity reached 3.0 km/sec, and the P-wave velocity was 5.4 km/sec. Shear-wave attenuation was computed using both the pulse rise time and spectral ratio methods. At station KNW, attenuation was significant only in an interval between 17.5 and approximately 40 m in depth. Over the top 50 m, attenuation corresponding to a Q of about 8 was obtained. A total T* of 0.004 sec was measured for this interval. Pulse rise times also increased rapidly in this zone. The spectral ratio data for station PFO yields two peaks in attenuation above 50 m. Similar to the attenuation found for station KNW, the peak in attenuation corresponds to a Q of about 11, averaged over the top 50 m. Spectra of the seismic pulses produced by the hammer give good signal between 20 to 80 Hz. Significant motion perpendicular to the polarizations of the first shear-wave arrival was recorded within a few meters of the surface. Apparently, the rock structure is sufficiently complicated that body waves are being converted (SH to SV at oblique incidence) very close to the surface. The presence of these elliptical particle motions within a mere few m of the pure shear-wave source suggests that the detection of polarizations perpendicular to the main shear arrival at a single location at the surface is not, by itself, a good method for detecting shearwave splitting within the upper few tens of kilometers of the earth's crust. Crack densities and orientations were determined from televiewer records. These records showed cracks with a preferred direction at station KNW and of a greater density than at station PFO. At station PFO, crack densities were smaller and more diffuse in orientation.

SV-P extraction and imaging for far-offset vertical seismic profile data

2015 ◽  
Vol 3 (3) ◽  
pp. SW27-SW35 ◽  
Author(s):  
Yandong Li ◽  
Bob A. Hardage
Keyword(s):  
Ray Tracing ◽  
Marcellus Shale ◽  
Seismic Profile ◽  
P Wave ◽  
Time Varying ◽  
Profile Data ◽  
Vertical Seismic Profile ◽  
Optimal Receiver ◽  
Common Depth Point ◽  
Vsp Data

We have analyzed vertical seismic profile (VSP) data acquired across a Marcellus Shale prospect and found that SV-P reflections could be extracted from far-offset VSP data generated by a vertical-vibrator source using time-variant receiver rotations. Optimal receiver rotation angles were determined by a dynamic steering of geophones to the time-varying approach directions of upgoing SV-P reflections. These SV-P reflections were then imaged using a VSP common-depth-point transformation based on ray tracing. Comparisons of our SV-P image with P-P and P-SV images derived from the same offset VSP data found that for deep targets, SV-P data created an image that extended farther from the receiver well than P-P and P-SV images and that spanned a wider offset range than P-P and P-SV images do. A comparison of our VSP SV-P image with a surface-based P-SV profile that traversed the VSP well demonstrated that SV-P data were equivalent to P-SV data for characterizing geology and that a VSP-derived SV-P image could be used to calibrate surface-recorded SV-P data that were generated by P-wave sources.

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.

Two methods for determining geophone orientations from VSP data

Geophysics ◽  
2006 ◽  
Vol 71 (4) ◽  
pp. V87-V97 ◽  
Author(s):  
Xiaoxian Zeng ◽  
George A. McMechan
Keyword(s):  
A Priori ◽  
Polarization Plane ◽  
Seismic Profile ◽  
P Wave ◽  
Earth Model ◽  
S Wave ◽  
Coherent Signals ◽  
Plane Method ◽  
Wave Separation ◽  
Vsp Data

Vertical seismic profile (VSP) data are usually acquired with three-component geophones of unknown azimuthal orientation. The geophone orientation must be estimated from the recorded data as a prerequisite to processing such as P- and S-wave separation, calculation of wave-incident directions, and 3D migration. We compare and combine two methods for estimating azimuthal orientation by least-squares fitting over a large number of shots. Combining the two methods can be done in an automated manner, which provides more accurate estimates of the geophone orientations than previous methods. In the polarization-plane method, we calculate the polarization plane of the first P-wave arrival. Then we subtract the source azimuth to determine the geophone orientation, independently for each geophone, with an angular uncertainty of [Formula: see text], and with no accumulated errors. In the relative-angle method, we obtain relative angles between adjacent geophone pairs using trace crosscorrelations, and operate on all coherent signals (even noise). Swapped geophone components can be detected automatically using the polarization-plane method. The main limitation of these (and all other known) methods is that uncertainties associated with path refraction are not estimated, unless some geophones have a priori known orientations, or we have a known earth model to correct for refraction.

Interpretation of total wave‐field data over Lost Hills field, Kern County, California

Geophysics ◽  
1989 ◽  
Vol 54 (11) ◽  
pp. 1420-1429 ◽  
Author(s):  
S. G. Squires ◽  
C. D. Y. Kim ◽  
D. Y. Kim
Keyword(s):  
Wave Field ◽  
Shear Wave ◽  
Delay Time ◽  
Field Data ◽  
P Wave ◽  
Kern County ◽  
Fracture Orientation ◽  
Total Wave ◽  
Fracture Intensity ◽  
Vsp Data

Approximately 5 miles (8 km) of total wave‐field data were acquired by Production Geophysical Services (then Kim Tech., Inc.), using the OMNIPULSE® Multi‐mode Shear‐wave Generator source over the southern end of Lost Hills field, Kern County, California. The quality of the shear‐wave sections was excellent. They represent a significant improvement over conventional P‐wave sections from the area in that they provide better reflection continuity and imaging of the Lost Hills anticline. A multicomponent VSP, which was acquired close to the line, provided crucial P‐wave to S‐wave correlation, as well as fracture information. [Formula: see text] ratios computed from interval times ranged from 2.79 to 1.63. An anomalously low [Formula: see text] ratio of 1.65 in the zone of interest (Lower Reef Ridge to McDonald shale), confirmed by multicomponent VSP data, corresponds to the producing interval. Evidence of shear‐wave splitting due to azimuthal anisotropy was observed, so the SV‐wave and SH‐wave data sets were rotated into principal‐component axes of N45E for S1 and N45W for S2. The predominant fracture orientation changes from N45E at depth to N45W near the surface. This change in fracture orientation with depth was confirmed by multicomponent VSP data. Delay‐time ratios (used as a measure of fracture intensity) ranged from a maximum of 11.71 percent to a minimum of −5.48 percent across the structure. These ratios are interpreted to show changes in fracture intensity and orientation across the structure. Delay‐time ratios in the zone of interest were anomalously high (1.55–6.53 percent). Comparison of fracture intensity on the flanks of the structure with that on the crest indicates that the flanks have undergone greater deformation than the crest. The total wave‐field data set and associated analyses have provided significant structural and stratigraphic information on the Miocene Monterey formation over the Lost Hills field, highlighting the productive interval.

Shear‐wave studies in glacial till

Geophysics ◽  
1998 ◽  
Vol 63 (4) ◽  
pp. 1273-1284 ◽  
Author(s):  
Bradley J. Carr ◽  
Zoltan Hajnal ◽  
Arnfinn Prugger
Keyword(s):  
Shear Wave ◽  
P Wave ◽  
Vertical Resolution ◽  
S Wave ◽  
Glacial Deposits ◽  
Seismic Line ◽  
S Waves ◽  
Additional Information ◽  
Vsp Data ◽  
Clay Percentage

Within a high‐resolution shallow reflection survey program in Saskatchewan, Canada, S-waves were produced using a single seismo‐electric blasting cap and were found to be distinguishable from surface wave phases. The local glacial deposits have average velocities of 450 m/s. [Formula: see text] ratios average 3.6 in these sequences, but they vary laterally, according to the velocity analyses done in two boreholes drilled along the seismic line. Vertical resolution for S-wave reflections are 0.75 m [in the vertical seismic profiling (VSP) data] and 1.5 m (in the CDP data). Yet, the S-wave CDP results are still better than corresponding P-wave data, which had a vertical resolution of 2.6 m. S-wave anisotropy is inferred in the glacial deposits on the basis of particle motion analysis and interpretations of S-wave splitting. However, the amount of observed splitting is small (∼2–6 ms over 5–10 m) and could go undetected for seismic surveys with larger sampling intervals. VSPs indicate that S-wave reflectivity is caused by both distinct and subtle lithologic changes (e.g., clay/sand contacts or changes in clay percentage within a particular till unit) and changes in bulk porosity. Migrated S-wave sections from line 1 and line 2 image reflections from sand layers within the tills as well as the first “bedrock” sequence (known as the Judith River Formation). Shear wave images are not only feasible in unconsolidated materials, but provide additional information about structural relationships within these till units.

Seismic shear‐wave observations in a physical model experiment

Geophysics ◽  
1983 ◽  
Vol 48 (6) ◽  
pp. 688-701 ◽  
Author(s):  
Robert H. Tatham ◽  
Donald V. Goolsbee ◽  
Wulf F. Massell ◽  
H. Roice Nelson
Keyword(s):  
Reflection Coefficient ◽  
Physical Model ◽  
Shear Wave ◽  
Model Experiment ◽  
Mode Conversion ◽  
Shear Waves ◽  
P Wave ◽  
Wave Reflections ◽  
Sea Floor ◽  
Physical Model Experiment

The observation and common‐depth‐point (CDP) processing of mode‐converted shear waves is demonstrated for real data collected in a physical model experiment. The model, submerged in water, represented water depth scaled to 250 ft, the first subsea reflector at 4000 ft, and the last reflector at 7000 ft below the sea floor with a structural wedge at the center. Very efficient mode conversion, from P to SV and back to P, is anticipated for angles of incidence at the liquid‐solid interface (sea floor) between 35 and 80 degrees. The model, constructed of Plexiglas and 3180 resin, will support elastic shear‐wave propagation. One anticipated problem, internal reflections from the sides of the model, was solved by tapering the sides of the model to 45 degrees off vertical. The P wave reflection coefficient at an interface between Plexiglas and water is 35 percent for vertical incidence, but it diminishes to very nearly zero between 43 and 75 degrees. Thus, by tapering the sides of the model, any undesired internal P wave reflections had to undergo at least two reflections at angles of incidence in the low reflection coefficient range for P waves. Data were collected in both an end‐on CDP mode, with offsets from 1000 ft to 20,000 ft, and a variety of walkaway experiments with scaled ranges from 1000 ft to 31,000 ft. Processing and analysis of the data confirm the existence of mode‐converted shear‐wave reflections in a modeled marine environment. In particular, the S wave reflections from all interfaces are identified on both the 100 percent gathered records and the final stacked records. These SV wave reflections were isolated for stacking by considering those portions of the gathered records, both offset and arrival time, that correspond to optimum angles of incidence. In addition, τ-p processing isolated particular angles of incidence, further confirming the incidence angle‐range criterion. Thus, the desired events are unambiguously identified as mode‐converted shear waves.

Sign in / Sign up

Export Citation Format

Share Document