P-wave and S-wave azimuthal anisotropy at a naturally fractured gas reservoir, Bluebell‐Altamont Field, Utah

Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1312-1328 ◽  
Author(s):  
Heloise B. Lynn ◽  
Wallace E. Beckham ◽  
K. Michele Simon ◽  
C. Richard Bates ◽  
M. Layman ◽  
...  

Reflection P- and S-wave data were used in an investigation to determine the relative merits and strengths of these two data sets to characterize a naturally fractured gas reservoir in the Tertiary Upper Green River formation. The objective is to evaluate the viability of P-wave seismic to detect the presence of gas‐filled fractures, estimate fracture density and orientation, and compare the results with estimates obtained from the S-wave data. The P-wave response to vertical fractures must be evaluated at different source‐receiver azimuths (travelpaths) relative to fracture strike. Two perpendicular lines of multicomponent reflection data were acquired approximately parallel and normal to the dominant strike of Upper Green River fractures as obtained from outcrop, core analysis, and borehole image logs. The P-wave amplitude response is extracted from prestack amplitude variation with offset (AVO) analysis, which is compared to isotropic‐model AVO responses of gas sand versus brine sand in the Upper Green River. A nine‐component vertical seismic profile (VSP) was also obtained for calibration of S-wave reflections with P-wave reflections, and support of reflection S-wave results. The direction of the fast (S1) shear‐wave component from the reflection data and the VSP coincides with the northwest orientation of Upper Green River fractures, and the direction of maximum horizontal in‐situ stress as determined from borehole ellipticity logs. Significant differences were observed in the P-wave AVO gradient measured parallel and perpendicular to the orientation of Upper Green River fractures. Positive AVO gradients were associated with gas‐producing fractured intervals for propagation normal to fractures. AVO gradients measured normal to fractures at known waterwet zones were near zero or negative. A proportional relationship was observed between the azimuthal variation of the P-wave AVO gradient as measured at the tops of fractured intervals, and the fractional difference between the vertical traveltimes of split S-waves (the “S-wave anisotropy”) of the intervals.

2016 ◽  
Vol 4 (3) ◽  
pp. SH21-SH37 ◽  
Author(s):  
Mathieu J. Duchesne ◽  
André J.-M. Pugin ◽  
Gabriel Fabien-Ouellet ◽  
Mathieu Sauvageau

The combined use of P- and S-wave seismic reflection data is appealing for providing insights into active petroleum systems because P-waves are sensitive to fluids and S-waves are not. The method presented herein relies on the simultaneous acquisition of P- and S-wave data using a vibratory source operated in the inline horizontal mode. The combined analysis of P- and S-wave reflections is tested on two potential hydrocarbon seeps located in a prospective area of the St. Lawrence Lowlands in Eastern Canada. For both sites, P-wave data indicate local changes in the reflection amplitude and slow velocities, whereas S-wave data present an anomalous amplitude at one site. Differences between P- and S-wave reflection morphology and amplitude and the abrupt decrease in P-velocity are indirect lines of evidence for hydrocarbon migration toward the surface through unconsolidated sediments. Surface-gas analysis made on samples taken at one potential seeping site reveals the occurrence of thermogenic gas that presumably vents from the underlying fractured Utica Shale forming the top of the bedrock. The 3C shear data suggest that fluid migration locally disturbs the elastic properties of the matrix. The comparative analysis of P- and S-wave data along with 3C recordings makes this method not only attractive for the remote detection of shallow hydrocarbons but also for the exploration of how fluid migration impacts unconsolidated geologic media.


2018 ◽  
Vol 16 (2) ◽  
pp. 19 ◽  
Author(s):  
Natashia Christy Viony ◽  
Wahyu Triyoso

The application of converted-wave seismic method in hydrocarbon exploration has increased significantly. Since the conventional seismic ceases to provide an adequate result in complex geology area and it provides an ambiguous brightspot response. The main principle is that an incident P-wave produces reflected and converted P and SV wave when the downgoing P-wave impinges on an interface. Converted-wave seismic uses the multicomponent receiver that records both of vertical component and horizontal component. The vertical component is assumed to correspond to the compressional PP wave and the horizontal correspond to the PS converted-wave. In this research, a synthetic model with the shallow gas and the salt dome below are constructed. The purpose of this study is to analyze the brightspot due to the presence of shallow gas and its effect to the quality of PP and PS wave reflection below the gas zone. To achieve the goal, both vertical and horizontal seismic data processing are performed. In horizontal data processing, the best gamma function (Vp/Vs) value is estimated to produce the better and reliable image. The result shows that the brightspot response in conventional data doesn’t exist in converted-wave data and the imaging below the gas zone in converted-wave data is better than the conventional due to the attenuation and diffraction effect that caused by gas column. Processing is followed by AVO analysis to compare the AVO response of PP and PS data in characterizing gas reservoir. Both PP and PS AVO curve shows the consistency with synthetic AVO from well data. Gas reservoir is a class 1 AVO anomaly with positive intercept and negative gradient on PP data. However, PS AVO curve does not refer any anomaly. It is because S-wave is not sensitive to the existence of rock saturant.


Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1277-1292 ◽  
Author(s):  
Robert E. Grimm ◽  
Heloise B. Lynn ◽  
C. R. Bates ◽  
D. R. Phillips ◽  
K. M. Simon ◽  
...  

Multiazimuth binning of 3-D P-wave reflection data is a relatively simple but robust way of characterizing the spatial distribution of gas‐producing natural fractures. In our survey, data were divided into two volumes by ray azimuth (approximately perpendicular and parallel (±45° to the dominant fracture strike) and separately processed. Azimuthal differences or ratios of attributes provided a rough measure of anisotropy. Improved imaging was also attained in the more coherent fracture‐parallel volume. A neural network using azimuthally dependent velocity, reflectivity, and frequency attributes identified commercial gas wells with greater than 85% success. Furthermore, we were able to interpret the physical mechanisms of most of these correlations and so better generalize the approach. The apparent velocity anisotropy was compared to that derived from other P- and S-wave methods in an inset three‐component survey. Prestack determination of the azimuthal moveout ellipse will best quantify velocity anisotropy, but simple two‐ or four‐azimuth poststack analysis can adequately identify regions of high fracture density and gas yield.


Geophysics ◽  
2015 ◽  
Vol 80 (1) ◽  
pp. C21-C35 ◽  
Author(s):  
Faranak Mahmoudian ◽  
Gary F. Margrave ◽  
Joe Wong ◽  
David C. Henley

We evaluated a quantitative amplitude analysis of 3D physical model reflection data acquired over an experimental phenolic layer that modeled a fractured medium with one set of vertical fractures. The phenolic layer was overlain by two isotropic layers, the uppermost being water, and the data acquisition was designed to avoid the interference of the primary and ghost events. The elastic stiffness coefficients and hence the anisotropy of the phenolic layer were known in advance from a previous traveltime analysis. The reflection amplitudes from the top of the phenolic layer required corrections to make them suitable for an amplitude study. In addition to the usual amplitude corrections applied to seismic field data, a directivity correction specific to the physical model transducers was applied. The corrected amplitudes along different azimuths showed a clear azimuthal variation caused by the phenolic layer and agreed with amplitudes predicted theoretically. An amplitude variation with angle and azimuth inversion was performed for horizontal transverse isotropy (HTI) parameters of the phenolic layer. We determined from the inversion results that from the azimuthally varying P-wave reflectivity response, it was possible to estimate HTI parameters that compared favorably to those obtained previously by a traveltime analysis. This result made it possible to compute the S-wave splitting parameter [Formula: see text] (historically determined from S-wave data and directly related to fracture density) from a quantitative analysis of the PP data.


2022 ◽  
Vol 41 (1) ◽  
pp. 47-53
Author(s):  
Zhiwen Deng ◽  
Rui Zhang ◽  
Liang Gou ◽  
Shaohua Zhang ◽  
Yuanyuan Yue ◽  
...  

The formation containing shallow gas clouds poses a major challenge for conventional P-wave seismic surveys in the Sanhu area, Qaidam Basin, west China, as it dramatically attenuates seismic P-waves, resulting in high uncertainty in the subsurface structure and complexity in reservoir characterization. To address this issue, we proposed a workflow of direct shear-wave seismic (S-S) surveys. This is because the shear wave is not significantly affected by the pore fluid. Our workflow includes acquisition, processing, and interpretation in calibration with conventional P-wave seismic data to obtain improved subsurface structure images and reservoir characterization. To procure a good S-wave seismic image, several key techniques were applied: (1) a newly developed S-wave vibrator, one of the most powerful such vibrators in the world, was used to send a strong S-wave into the subsurface; (2) the acquired 9C S-S data sets initially were rotated into SH-SH and SV-SV components and subsequently were rotated into fast and slow S-wave components; and (3) a surface-wave inversion technique was applied to obtain the near-surface shear-wave velocity, used for static correction. As expected, the S-wave data were not affected by the gas clouds. This allowed us to map the subsurface structures with stronger confidence than with the P-wave data. Such S-wave data materialize into similar frequency spectra as P-wave data with a better signal-to-noise ratio. Seismic attributes were also applied to the S-wave data sets. This resulted in clearly visible geologic features that were invisible in the P-wave data.


Geophysics ◽  
2006 ◽  
Vol 71 (5) ◽  
pp. D161-D170 ◽  
Author(s):  
Xiaoxia Xu ◽  
Ilya Tsvankin

Compensation for geometrical spreading along a raypath is one of the key steps in AVO (amplitude-variation-with-offset) analysis, in particular, for wide-azimuth surveys. Here, we propose an efficient methodology to correct long-spread, wide-azimuth reflection data for geometrical spreading in stratified azimuthally anisotropic media. The P-wave geometrical-spreading factor is expressed through the reflection traveltime described by a nonhyperbolic moveout equation that has the same form as in VTI (transversely isotropic with a vertical symmetry axis) media. The adapted VTI equation is parameterized by the normal-moveout (NMO) ellipse and the azimuthally varying anellipticity parameter [Formula: see text]. To estimate the moveout parameters, we apply a 3D nonhyperbolic semblance algorithm of Vasconcelos and Tsvankin that operates simultaneously with traces at all offsets andazimuths. The estimated moveout parameters are used as the input in our geometrical-spreading computation. Numerical tests for models composed of orthorhombic layers with strong, depth-varying velocity anisotropy confirm the high accuracy of our travetime-fitting procedure and, therefore, of the geometrical-spreading correction. Because our algorithm is based entirely on the kinematics of reflection arrivals, it can be incorporated readily into the processing flow of azimuthal AVO analysis. In combination with the nonhyperbolic moveout inversion, we apply our method to wide-azimuth P-wave data collected at the Weyburn field in Canada. The geometrical-spreading factor for the reflection from the top of the fractured reservoir is clearly influenced by azimuthal anisotropy in the overburden, which should cause distortions in the azimuthal AVO attributes. This case study confirms that the azimuthal variation of the geometrical-spreading factor often is comparable to or exceeds that of the reflection coefficient.


2021 ◽  
Author(s):  
Janneke de Laat ◽  
Sergei Lebedev ◽  
Bruna Chagas de Melo ◽  
Nicolas Celli ◽  
Raffaele Bonadio

<p>We present a new S-wave velocity tomographic model of the Australian Plate, Aus21.  It is constrained by waveforms of 0.9 million seismograms with both the corresponding sources and stations located within the half-hemisphere centred at the Australian continent. Waveform inversion extracts structural information from surface, S- and multiple S-waves on the seismograms in the form of a set of linear equations. These equations are then combined in a large linear system and inverted jointly to obtain a tomographic model of S- and P-wave speeds and S-wave azimuthal anisotropy of the crust and upper mantle. The model has been validated by resolution tests and, for particular locations in Australia with notable differences with previous models, by independent inter-station measurements of surface-wave phase velocities, which we performed using available array data. </p><p> </p><p>Aus21 offers new insights into the structure and evolution of the Australian Plate and its boundaries. The Australian cratonic lithosphere occupies nearly all of the western and central Australia but shows substantial lateral heterogeneity. It extends up to the northern edge of the plate, where it is colliding with island arcs, without subducting. The rugged eastern boundary of the cratonic lithosphere provides a lithospheric definition of the Tasman Line. The thin, warm lithosphere below the eastern part of the continent, east of the Tasman Line, underlies the Cenozoic volcanism locations in the area. The lithosphere is also thin and warm below much of the Tasman Sea, underlying the Lord Howe hotspot and the submerged part of western Zealandia. A low velocity anomaly that may indicate the single source of the Lord Howe and Tasmanid hotspots is observed in the transition zone offshore the Australian continent, possibly also sourcing the East Australia hotspot. Another potential hotspot source is identified below the Kermadec Trench, causing an apparent slab gap in the overlying slab and possibly related to the Samoa Hotspot to the north. Below a portion of the South East Indian Ridge (the southern boundary of the Australian Plate) a pronounced high velocity anomaly is present in the 200-400 km depth range just east of the Australian-Antarctic Discordance (AAD), probably linked to the evolution of this chaotic ridge system.</p>


2002 ◽  
Vol 42 (1) ◽  
pp. 627
Author(s):  
R.G. Williams ◽  
G. Roberts ◽  
K. Hawkins

Seismic energy that has been mode converted from pwave to s-wave in the sub-surface may be recorded by multi-component surveys to obtain information about the elastic properties of the earth. Since the energy converted to s-wave is missing from the p-wave an alternative to recording OBC multi-component data is to examine p-wave data for the missing energy. Since pwave velocities are generally faster than s-wave velocities, then for a given reflection point the converted s-wave signal reaches the surface at a shorter offset than the equivalent p-wave information. Thus, it is necessary to record longer offsets for p-wave data than for multicomponent data in order to measure the same information.A non-linear, wide-angle (including post critical) AVO inversion has been developed that allows relative changes in p-wave velocities, s-wave velocities and density to be extracted from long offset p-wave data. To extract amplitudes at long offsets for this inversion it is necessary to image the data correctly, including correcting for higher order moveout and possibly anisotropy if it is present.The higher order moveout may itself be inverted to yield additional information about the anisotropy of the sub-surface.


Sign in / Sign up

Export Citation Format

Share Document