Land based S-wave reflection seismology with P sources — does it work?

Author(s):  
Bob A. Hardage
2021 ◽  
pp. 1-49
Author(s):  
Bob Hardage ◽  
Mike Graul ◽  
Tim Hall ◽  
Chris Hall ◽  
Mark Kelley ◽  
...  

We compared two methods for extracting the azimuth of maximum horizontal stress (SHmax) from 3D land-based seismic data generated by a P source and recorded with vertical geophones. In the first method, we used the direct-SV mode that is produced by all land-based P sources. P sources generate SV illumination that radiates in all azimuth directions from a source station and creates SV-P reflections that are recorded by vertical geophones. Unless stratigraphy has steep dip, SV-P raypaths recorded by vertical geophones are the reverse of P-SV raypaths recorded by horizontal geophones. Thus, SV-P data provide the same S-wave sensitivity to stress fields as popular P-SV data do. In the second method, we retrieved P-P reflections and then performed an amplitude-versus-incident-angle (AVA) analysis of the amplitude-gradient behavior of P-P reflection wavelets. We did this analysis in narrow azimuth corridors to determine the gradient of reflection-wavelet amplitudes as a function of azimuth. This P-P AVA amplitude-gradient method has been of great interest in the reflection seismology community since it was introduced in the late 1990s. Each of these methods, AVA analysis of the gradient of P-P reflection amplitudes and azimuth-dependent arrival times of SV-P reflections can be used to determine the azimuth of SHmax stress. We compare the results of the two methods with ground truth measurements of SHmax azimuth at a CO2 sequestration site in the Michigan Basin. SHmax azimuths were determined from P-P and SV-P data at three major boundaries at depths of approximately 3500 ft (1067 m), 5500 ft (1676 m), and 7500 ft (2286 m). Two estimates of SHmax azimuth (one using SV-P data and one using P-P data) were made at each stacking bin inside a 24 mi2 (62 km2) image space. The result was approximately 98,000 estimates of SHmax azimuth across each of these three boundaries for each of these two prediction strategies. Histogram displays of PP AVA gradient estimates had peaks at correct azimuths of SHmax at all three depths, but the spread of the distributions widened with depth and split into two peaks at the deepest boundary. In contrast, each histogram of SHmax azimuth predicted by azimuth-dependent SV-P traveltimes had a single, definitive peak that was positioned at the correct SHmax azimuth at all three boundary depths.


Geophysics ◽  
2009 ◽  
Vol 74 (2) ◽  
pp. Q1-Q16 ◽  
Author(s):  
Mark A. Meier ◽  
Paul J. Lee

Seismic-resolution theory used in survey design applies to P- and S-wave reflection seismology but is not readily applicable to converted waves. We rederived the theory with sufficient generality to include converted waves explicitly. Generalization of the theory requires that the inherent converted-wave properties of mode change and finite scattering angle be accommodated. The significance of amplitude fidelity in modern seismic applications also affects the resolution description. We considered resolution in the context of a single recorded trace (trace resolution) and an image derived from many traces (image resolution). Trace resolution is governed by wavelet width, which results in a minimum bandwidth requirement, depending on scattering angle and local propagation velocities of incident and scattered modes. Image resolution is governed additionally by the migration acceptance angle and results in minimum aperture and sampling requirements. The requirements for converted-wave surveys generally differ from those of P- and S-wave surveys. Our theory predicts converted-wave reflection seismology resolution and provides the minimum acquisition parameters required to achieve resolution objectives.


Geophysics ◽  
2002 ◽  
Vol 67 (3) ◽  
pp. 690-700 ◽  
Author(s):  
Josef Holzschuh

Compressional (P) wave and shear (S) wave seismic reflection techniques were used to delineate the sand and gravel aquifer within a highly saline clay‐filled paleochannel in the Eastern Goldfields of Western Australia. The seismic refraction and gravity methods were also used to investigate the paleochannel. The unsaturated loose fine‐grained sand up to 10 m in depth at the surface is a major factor in degrading subsurface imaging. The seismic processing needed to be precise, with accurate static corrections and normal moveout corrections. Deconvolution enhanced the aquifer and other paleochannel reflectors. P‐wave reflection and refraction layer depths had good correlation and showed a total of six boundaries: (1) water table, (2) change in velocity (compaction) in the paleochannel sediments, (3) sand and gravel aquifer, (4) red‐brown saprolite and green saprolite boundary, (5) weathered bedrock, and (6) unweathered bedrock. P‐wave explosive and hammer sources were found to have similar signal characteristics, and the aquifer and bedrock were both imaged using the hammer source. The deep shots below the water table have the most broadband frequency response for reflections, but stacking clear reflections was difficult. The S‐wave reflection results showed high lateral and vertical resolution of the basal saprolite clay, the sand and gravel aquifer, and very shallow clays above the aquifer. The S‐wave reflection stacking velocities were 10–20% of the P‐waves, increasing the resolution of the S‐wave section. The gravity data were modelled to fit the known drilling and P‐wave seismic reflection depths. The refraction results did not identify the top of bedrock, so refraction depths were not used for the gravity modeling in this highly weathered environment. The final gravity model mapped the bedrock topography beyond the lateral extent of the seismic and drilling data.


2016 ◽  
Vol 4 (4) ◽  
pp. T613-T625 ◽  
Author(s):  
Qizhen Du ◽  
Bo Zhang ◽  
Xianjun Meng ◽  
Chengfeng Guo ◽  
Gang Chen ◽  
...  

Three-term amplitude-variation with offset (AVO) inversion generally suffers from instability when there is limited prior geologic or petrophysical constraints. Two-term AVO inversion shows higher instability compared with three-term AVO inversion. However, density, which is important in the fluid-type estimation, cannot be recovered from two-term AVO inversion. To reliably predict the P- and S-waves and density, we have developed a robust two-step joint PP- and PS-wave three-term AVO-inversion method. Our inversion workflow consists of two steps. The first step is to estimate the P- and S-wave reflectivities using Stewart’s joint two-term PP- and PS-AVO inversion. The second step is to treat the P-wave reflectivity obtained from the first step as the prior constraint to remove the P-wave velocity related-term from the three-term Aki-Richards PP-wave approximated reflection coefficient equation, and then the reduced PP-wave reflection coefficient equation is combined with the PS-wave reflection coefficient equation to estimate the S-wave and density reflectivities. We determined the effectiveness of our method by first applying it to synthetic models and then to field data. We also analyzed the condition number of the coefficient matrix to illustrate the stability of the proposed method. The estimated results using proposed method are superior to those obtained from three-term AVO inversion.


2000 ◽  
Author(s):  
Deidra K. Begay ◽  
Richard D. Miller ◽  
W. Lynn Watney ◽  
Jianghai Xia

Sign in / Sign up

Export Citation Format

Share Document