Estimating shear‐wave velocities from P-wave and converted‐wave data

Geophysics ◽  
1999 ◽  
Vol 64 (2) ◽  
pp. 504-507 ◽  
Author(s):  
Franklyn K. Levin
Keyword(s):  
Wave Velocity ◽  
Seismic Data ◽  
P Wave ◽  
S Wave ◽  
Converted Wave ◽  
Wave Data ◽  
Wave Velocities ◽  
Shear Wave Velocities ◽  
P Wave Velocity ◽  
Growing Body

Tessmer and Behle (1988) show that S-wave velocity can be estimated from surface seismic data if both normal P-wave data and converted‐wave data (P-SV) are available. The relation of Tessmer and Behle is [Formula: see text] (1) where [Formula: see text] is the S-wave velocity, [Formula: see text] is the P-wave velocity, and [Formula: see text] is the converted‐wave velocity. The growing body of converted‐wave data suggest a brief examination of the validity of equation (1) for velocities that vary with depth.

Poststack, prestack, and joint inversion of P- and S-wave data for Morrow A sandstone characterization

2015 ◽  
Vol 3 (3) ◽  
pp. SZ59-SZ92 ◽  
Author(s):  
Paritosh Singh ◽  
Thomas L. Davis ◽  
Bryan DeVault
Keyword(s):  
Wave Velocity ◽  
Seismic Data ◽  
Joint Inversion ◽  
Wave Impedance ◽  
P Wave ◽  
S Wave ◽  
Wave Data ◽  
Wave Inversion ◽  
Density Maps ◽  
P Wave Velocity

Exploration for oil-bearing Morrow sandstones using conventional seismic data/methods has a startlingly low success rate of only 3%. The S-wave velocity contrast between the Morrow shale and A sandstone is strong compared with the P-wave velocity contrast, and, therefore, multicomponent seismic data could help to characterize these reservoirs. The SV and SH data used in this study are generated using S-wave data from horizontal source and horizontal receiver recording. Prestack P- and S-wave inversions, and joint P- and S-wave inversions, provide estimates of P- and S-wave impedances, and density for characterization of the Morrow A sandstone. Due to the weak P-wave amplitude-versus-angle response at the Morrow A sandstone top, the density and S-wave impedance estimated from joint P- and S-wave inversions were inferior to the prestack S-wave inversion. The inversion results were compared with the Morrow A sandstone thickness and density maps obtained from well logs to select the final impedance and density volume for interpretation. The P-wave impedance estimated from prestack P-wave data, as well as density and S-wave impedance estimated from prestack SV‐wave data were used to identify the distribution, thickness, quality, and porosity of the Morrow A sandstone. The stratal slicing method was used to get the P- and S-wave impedances and density maps. The S-wave impedance characterizes the Morrow A sandstone distribution better than the P-wave impedance throughout the study area. Density estimation from prestack inversion of SV data was able to distinguish between low- and high-quality reservoirs. The porosity volume was estimated from the density obtained from prestack SV-wave inversion. We found some possible well locations based on the interpretation.

Using Seismic Data to Predict Shale Pore Pressure and Overpressure Sweet Spots: A Case Study from the Lower Silurian Longmaxi Formation in Sichuan Basin, China

2021 ◽  
Author(s):  
Sheng Chen ◽  
Qingcai Zeng ◽  
Xiujiao Wang ◽  
Qing Yang ◽  
Chunmeng Dai ◽  
...  
Keyword(s):  
Prediction Model ◽  
Wave Velocity ◽  
Shale Gas ◽  
Seismic Data ◽  
P Wave ◽  
Formation Pressure ◽  
S Wave ◽  
P Wave Velocity ◽  
New Formation ◽  
Sweet Spots

Abstract Practices of marine shale gas exploration and development in south China have proved that formation overpressure is the main controlling factor of shale gas enrichment and an indicator of good preservation condition. Accurate prediction of formation pressure before drilling is necessary for drilling safety and important for sweet spots predicting and horizontal wells deploying. However, the existing prediction methods of formation pore pressures all have defects, the prediction accuracy unsatisfactory for shale gas development. By means of rock mechanics analysis and related formulas, we derived a formula for calculating formation pore pressures. Through regional rock physical analysis, we determined and optimized the relevant parameters in the formula, and established a new formation pressure prediction model considering P-wave velocity, S-wave velocity and density. Based on regional exploration wells and 3D seismic data, we carried out pre-stack seismic inversion to obtain high-precision P-wave velocity, S-wave velocity and density data volumes. We utilized the new formation pressure prediction model to predict the pressure and the spatial distribution of overpressure sweet spots. Then, we applied the measured pressure data of three new wells to verify the predicted formation pressure by seismic data. The result shows that the new method has a higher accuracy. This method is qualified for safe drilling and prediction of overpressure sweet spots for shale gas development, so it is worthy of promotion.

Elastic reflectivity vectors and the geometry of intercept-gradient crossplots

2019 ◽  
Vol 38 (10) ◽  
pp. 762-769
Author(s):  
Patrick Connolly
Keyword(s):  
Elastic Property ◽  
Elastic Properties ◽  
Wave Velocity ◽  
Seismic Data ◽  
Measurement Errors ◽  
Rock Physics ◽  
P Wave ◽  
S Wave ◽  
P Wave Velocity ◽  
Well Log Analysis

Reflectivities of elastic properties can be expressed as a sum of the reflectivities of P-wave velocity, S-wave velocity, and density, as can the amplitude-variation-with-offset (AVO) parameters, intercept, gradient, and curvature. This common format allows elastic property reflectivities to be expressed as a sum of AVO parameters. Most AVO studies are conducted using a two-term approximation, so it is helpful to reduce the three-term expressions for elastic reflectivities to two by assuming a relationship between P-wave velocity and density. Reduced to two AVO components, elastic property reflectivities can be represented as vectors on intercept-gradient crossplots. Normalizing the lengths of the vectors allows them to serve as basis vectors such that the position of any point in intercept-gradient space can be inferred directly from changes in elastic properties. This provides a direct link between properties commonly used in rock physics and attributes that can be measured from seismic data. The theory is best exploited by constructing new seismic data sets from combinations of intercept and gradient data at various projection angles. Elastic property reflectivity theory can be transferred to the impedance domain to aid in the analysis of well data to help inform the choice of projection angles. Because of the effects of gradient measurement errors, seismic projection angles are unlikely to be the same as theoretical angles or angles derived from well-log analysis, so seismic data will need to be scanned through a range of angles to find the optimum.

Wavefield simulation and data-acquisition-scheme analysis for LWD acoustic tools in very slow formations

Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. E59-E68 ◽  
Author(s):  
Hua Wang ◽  
Guo Tao
Keyword(s):  
Wave Velocity ◽  
Cutoff Frequency ◽  
P Wave ◽  
Flexural Wave ◽  
S Wave ◽  
Low Frequencies ◽  
Wave Velocities ◽  
Dipole System ◽  
P Wave Velocity ◽  
Wavefield Simulation

Propagating wavefields from monopole, dipole, and quadrupole acoustic logging-while-drilling (LWD) tools in very slow formations have been studied using the discrete wavenumber integration method. These studies examine the responses of monopole and dipole systems at different source frequencies in a very slow surrounding formation, and the responses of a quadrupole system operating at a low source frequency in a slow formation with different S-wave velocities. Analyses are conducted of coherence-velocity/slowness relationships (semblance spectra) in the time domain and of the dispersion characteristics of these waveform signals from acoustic LWD array receivers. These analyses demonstrate that, if the acoustic LWD tool is centralized properly and is operating at low frequencies (below 3 kHz), a monopole system can measure P-wave velocity by means of a “leaky” P-wave for very slow formations. Also, for very slow formations a dipole system can measure the P-wave velocity via a leaky P-wave and can measure the S-wave velocity from a formation flexural wave. With a quadrupole system, however, the lower frequency limit (cutoff frequency) of the drill-collar interference wave would decrease to 5 kHz and might no longer be neglected if the surrounding formation becomes a very slow formation, with S-wave velocities at approximately 500 m/s.

Acoustic measurements in unconsolidated sands at low effective pressure and overpressure detection

Geophysics ◽  
2002 ◽  
Vol 67 (2) ◽  
pp. 405-412 ◽  
Author(s):  
Manika Prasad
Keyword(s):  
Wave Velocity ◽  
P Wave ◽  
Acoustic Measurements ◽  
S Wave ◽  
Deepwater Drilling ◽  
Unconsolidated Sands ◽  
Wave Velocities ◽  
P Wave Velocity ◽  
Shallow Water Flows ◽  
Low Pressures

Shallow water flows and over‐pressured zones are a major hazard in deepwater drilling projects. Their detection prior to drilling would save millions of dollars in lost drilling costs. I have investigated the sensitivity of seismic methods for this purpose. Using P‐wave information alone can be ambiguous, because a drop in P‐wave velocity (Vp) can be caused both by overpressure and by presence of gas. The ratio of P‐wave velocity to S‐wave velocity (Vp/Vs), which increases with overpressure and decreases with gas saturation, can help differentiate between the two cases. Since P‐wave velocity in a suspension is slightly below that of the suspending fluid and Vs=0, Vp/Vs and Poisson's ratio must increase exponentially as a load‐bearing sediment approaches a state of suspension. On the other hand, presence of gas will also decrease Vp but Vs will remain unaffected and Vp/Vs will decrease. Analyses of ultrasonic P‐ and S‐wave velocities in sands show that the Vp/Vs ratio, especially at low effective pressures, decreases rapidly with pressure. At very low pressures, Vp/Vs values can be as large as 100 and higher. Above pressures greater than 2 MPa, it plateaus and does not change much with pressure. There is significant change in signal amplitudes and frequency of shear waves below 1 MPa. The current ultrasonic data shows that Vp/Vs values can be invaluable indicators of low differential pressures.

Converted-wave model building and imaging based on common-focus-point methodology

Geophysics ◽  
2020 ◽  
Vol 85 (6) ◽  
pp. U139-U149
Author(s):  
Hongwei Liu ◽  
Mustafa Naser Al-Ali ◽  
Yi Luo
Keyword(s):  
Elastic Wave ◽  
Wave Velocity ◽  
Seismic Data ◽  
Model Building ◽  
P Wave ◽  
Rock Properties ◽  
S Wave ◽  
Converted Wave ◽  
Seismic Images ◽  
Focus Point

Seismic images can be viewed as photographs for underground rocks. These images can be generated from different reflections of elastic waves with different rock properties. Although the dominant seismic data processing is still based on the acoustic wave assumption, elastic wave processing and imaging have become increasingly popular in recent years. A major challenge in elastic wave processing is shear-wave (S-wave) velocity model building. For this reason, we have developed a sequence of procedures for estimating seismic S-wave velocities and the subsequent generation of seismic images using converted waves. We have two main essential new supporting techniques. The first technique is the decoupling of the S-wave information by generating common-focus-point gathers via application of the compressional-wave (P-wave) velocity on the converted seismic data. The second technique is to assume one common VP/ VS ratio to approximate two types of ratios, namely, the ratio of the average earth layer velocity and the ratio of the stacking velocity. The benefit is that we reduce two unknown ratios into one, so it can be easily scanned and picked in practice. The PS-wave images produced by this technology could be aligned with the PP-wave images such that both can be produced in the same coordinate system. The registration between the PP and PS images provides cross-validation of the migrated structures and a better estimation of underground rock and fluid properties. The S-wave velocity, computed from the picked optimal ratio, can be used not only for generating the PS-wave images, but also to ensure well registration between the converted-wave and P-wave images.

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.

Wavelet distortion correction due to domain conversion

Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. V77-V87 ◽  
Author(s):  
Rishi Bansal ◽  
Mike Matheney
Keyword(s):  
Wave Velocity ◽  
Seismic Data ◽  
Seismic Velocity ◽  
Distortion Correction ◽  
P Wave ◽  
Converted Wave ◽  
Shaping Filter ◽  
Seismic Wavelet ◽  
Attribute Analysis ◽  
P Wave Velocity

Converted-wave (PS) data, when converted to PP time, develop time- and location-varying compression of the seismic wavelet due to a variable subsurface [Formula: see text] [Formula: see text]. The time-dependent compression distorts the wavelet in a seismic trace. The lack of a consistent seismic wavelet in a domain-converted PS volume can eventually lead to an erroneous joint PP/PS inversion result. Depth-converted seismic data also have wavelet distortion due to velocity-dependent wavelet stretch. A high value of seismic velocity produces more stretch in a seismic wavelet than a low value. Variable wavelet stretch renders the depth data unsuitable for attribute analysis. A filtering scheme is proposed that corrects for distortion in seismic wavelets due to domain conversions (PS to PP time and time-to-depth) of seismic data in an amplitude-preserving manner. The method uses a Fourier scaling theorem to predict the seismic wavelet in the converted domain and calculates a shaping filter for each time/depth sample that corrects for the distortion in the wavelet. The filter is applied to the domain-converted data using the method of nonstationary filtering. We provide analytical expressions for the squeeze factor [Formula: see text] that is used to predict the wavelet in the converted domain. The squeeze factor [Formula: see text] for PS to PP time conversion is a function of the subsurface [Formula: see text] whereas for PP time-to-depth conversion [Formula: see text] is dependent on subsurface P-wave velocity. After filtering, the squeezed wavelets in domain-converted PS data appear to have resulted from a constant subsurface [Formula: see text], which we denote as [Formula: see text]. Similarly, the filtered depth-converted data appear to have resulted from a constant subsurface P-wave velocity [Formula: see text].

Determining P- and S-wave velocities and Q-values from single ultrasound transmission measurements performed on cylindrical rock samples: it’s possible, when…

2020 ◽  
Author(s):  
Marc S. Boxberg ◽  
Mandy Duda ◽  
Katrin Löer ◽  
Wolfgang Friederich ◽  
Jörg Renner
Keyword(s):  
Wave Velocity ◽  
P Wave ◽  
S Wave ◽  
Standard Material ◽  
Rock Samples ◽  
Intrinsic Attenuation ◽  
Finite Element Software ◽  
Wave Velocities ◽  
P Wave Velocity

<p>Determining elastic wave velocities and intrinsic attenuation of cylindrical rock samples by transmission of ultrasound signals appears to be a simple experimental task, which is performed routinely in a range of geoscientific and engineering applications requiring characterization of rocks in field and laboratory. P- and S-wave velocities are generally determined from first arrivals of signals excited by specifically designed transducers. A couple of methods exist for determining the intrinsic attenuation, most of them relying either on a comparison between the sample under investigation and a standard material or on investigating the same material for various geometries.</p><p>Of the three properties of interest, P-wave velocity is certainly the least challenging one to determine, but dispersion phenomena lead to complications with the consistent identification of frequency-dependent first breaks. The determination of S-wave velocities is even more hampered by converted waves interfering with the S-wave arrival. Attenuation estimates are generally subject to higher uncertainties than velocity measurements due to the high sensitivity of amplitudes to experimental procedures. The achievable accuracy of determining S-wave velocity and intrinsic attenuation using standard procedures thus appears to be limited.</p><p>We pursue the determination of velocity and attenuation of rock samples based on full waveform modeling and inversion. Assuming the rock sample to be homogeneous - an assumption also underlying standard analyses - we quantify P-wave velocity, S-wave velocity and intrinsic P- and S-wave attenuation from matching a single ultrasound trace with a synthetic one numerically modelled using the spectral finite-element software packages SPECFEM2D and SPECFEM3D. We find that enough information on both velocities is contained in the recognizable reflected and converted phases even when nominal P-wave sensors are used. Attenuation characteristics are also inherently contained in the relative amplitudes of these phases due to their different travel paths. We present recommendations for and results from laboratory measurements on cylindrical samples of aluminum and rocks with different geometries that we also compare with various standard analysis methods. The effort put into processing for our approach is particularly justified when accurate values and/or small variations, for example in response to changing P-T-conditions, are of interest or when the amount of sample material is limited.</p>

Predicting reservoir quality in the Bakken Formation, North Dakota, using petrophysics and 3C seismic data

2020 ◽  
Vol 8 (4) ◽  
pp. T851-T868
Author(s):  
Andrea G. Paris ◽  
Robert R. Stewart
Keyword(s):  
North Dakota ◽  
Seismic Data ◽  
Rock Physics ◽  
P Wave ◽  
Well Logs ◽  
Reservoir Quality ◽  
Bakken Formation ◽  
S Wave ◽  
Converted Wave ◽  
Wave Velocities

Combining rock-property analysis with multicomponent seismic imaging can be an effective approach for reservoir quality prediction in the Bakken Formation, North Dakota. The hydrocarbon potential of shale is indicated on well logs by low density, high gamma-ray response, low compressional-wave (P-wave) and shear-wave (S-wave) velocities, and high neutron porosity. We have recognized the shale intervals by cross plotting sonic velocities versus density. Intervals with total organic carbon (TOC) content higher than 10 wt% deviate from lower TOC regions in the density domain and exhibit slightly lower velocities and densities (<2.30 g/cm3). We consider TOC to be the principal factor affecting changes in the density and P- and S-wave velocities in the Bakken shales, where VP/ VS ranges between 1.65 and 1.75. We generate the synthetic seismic data using an anisotropic version of the Zoeppritz equations, including estimated Thomsen’s parameters. For the tops of the Upper and Lower Bakken, the amplitude shows a negative intercept and a positive gradient, which corresponds to an amplitude variation with offset of class IV. The P-impedance error decreases by 14% when incorporating the converted-wave information in the inversion process. A statistical approach using multiattribute analysis and neural networks delimits the zones of interest in terms of P-impedance, density, TOC content, and brittleness. The inverted and predicted results show reasonable correlations with the original well logs. The integration of well log analysis, rock physics, seismic modeling, constrained inversions, and statistical predictions contributes to identifying the areas of highest reservoir quality within the Bakken Formation.

Sign in / Sign up

Export Citation Format

Share Document