Coarse‐layer stripping of vertically variable azimuthal anisotropy from shear‐wave data

Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1126-1138 ◽  
Author(s):  
Leon Thomsen ◽  
Ilya Tsvankin ◽  
Michael C. Mueller

Alford rotation analysis of 2C × 2C shear‐wave data (two source components, two receiver components) for azimuthal anisotropy is valid only when the orientation of that azimuthal anisotropy is invariant with depth. The Winterstein and Meadows method of layer stripping vertical seismic profiling (VSP) data relaxes this restriction for coarse‐layer variation of the orientation of the anisotropy. Here we present a tensor generalization of the conventional convolutional model of scalar wave propagation and use it to derive generalizations of Winterstein and Meadows layer stripping, valid for 2C × 2C data and for the restricted 2C-only case, in the VSP and reflection contexts. In the 2C × 2C VSP application, the result reduces to that of Winterstein and Meadows in the case where both fast and slow shear modes have the same attenuation and dispersion; otherwise, a balancing of mode spectra and amplitudes is required. The 2C × 2C reflection result differs from the 2C × 2C VSP result since two applications of the mode‐balancing and mode‐advance operations are required (since the waves travel up as well as down). Application to a synthetic data set confirms these results. The 2C × 2C reflection algorithm enables the exploration for sweet spots of high fracture intensity ahead of the bit without the restrictive assumption that the anisotropy orientation is depth invariant.

Geophysics ◽  
2008 ◽  
Vol 73 (5) ◽  
pp. A39-A43 ◽  
Author(s):  
Andrey Bakulin ◽  
Albena Mateeva

Measuring shear-wave splitting from vertical seismic profiling (VSP) data can benefit fracture and stress characterization as well as seismic processing and interpretation. The classic approach to measuring azimuthal anisotropy at depth involves layer stripping. Its inherent weakness is the need to measure and undo overburden effects before arriving at an anisotropy estimate at depth. That task is challenging when the overburden is complex and varies quickly with depth. Moreover, VSP receivers are rarely present all the way from the surface to the target. That necessitates the use of simplistic assumptions about the uninstrumented part of the overburden that limit the quality of the result. We propose a new technique for measuring shear-wave splitting at depth that does not require any knowledge of the overburden. It is based on a multicomponent version of the virtual source method in which each two-component (2-C) VSP receiver is turned into a 2-C shear source and recorded at deeper geophones. The resulting virtual data set is affected only by the properties of the medium between the receivers. A simple Alford rotation transforms the data set into fast and slow shear virtual check shots from which shear-wave splitting can be measured easily and accurately under arbitrarily complex overburden.


Geophysics ◽  
2011 ◽  
Vol 76 (3) ◽  
pp. B79-B87 ◽  
Author(s):  
Christian Stotter ◽  
Erika Angerer

A 2D vibroseis line was acquired in the Vienna Basin (Austria) for the purpose of comparing the data of digital multicomponent single sensors based on microelectromechanical system (MEMS) sensors alongside conventional vertical-component geophone arrays. For efficient removal of coherent noise during processing, all source points were recorded in single-sweep mode, i.e., no vertical stacking was performed in the field. On this densely sampled data set, several noise-reduction techniques, such as digital array forming, frequency-wavenumber (f-k) filtering in shot and receiver domains, and polarization filters, proved to be valuable in reducing source-generated noise. The results showed that, with the use of single-sweep recording and polarization filter techniques, it is possible to produce seismic sections for a single-receiver three-component (3C) MEMS line that are comparable to a conventional geophone array line in signal-to-noise ratio. However, the higher number of single geophones and hence the stronger attenuation of random noise in the conventional array resulted in an advantage for the analog geophone data set. The second goal for this survey was to evaluate additional information contained in the horizontal components of the MEMS data. The multicomponent data allowed for the processing of mode-converted shear-wave data, performed for the first time in the Vienna Basin. Azimuthal anisotropy related to horizontal stresses was observed in the Neogene section of the shear-wave data set. A PP-PS event correlation allowed the identification of major shallow horizons. Interpretation of the final sections confirmed that the PS data are useful to distinguish between gas reservoirs and high-porosity water sands, which can cause similar P-wave amplitude variation with offset (AVO) effects.


Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1172-1180 ◽  
Author(s):  
W. Scott Leaney ◽  
Colin M. Sayers ◽  
Douglas E. Miller

Multioffset vertical seismic profile (VSP) experiments, commonly referred to as walkaways, enable anisotropy to be measured reliably in the field. The results can be fed into modeling programs to study the impact of anisotropy on velocity analysis, migration, and amplitude versus offset (AVO). Properly designed multioffset VSPs can also provide the target AVO response measured under optimum conditions, since the wavelet is recorded just above the reflectors of interest with minimal reflection point dispersal. In this paper, the multioffset VSP technique is extended to include multioffset azimuths, and a multiazimuthal multiple VSP data set acquired over a carbonate reservoir is analyzed for P-wave anisotropy and AVO. Direct arrival times down to the overlying shale and reflection times and amplitudes from the carbonate are analyzed. Data analysis involves a three‐term fit to account for nonhyperbolic moveout, dip, and azimuthal anisotropy. Results indicate that the overlying shale is transversely isotropic with a vertical axis of symmetry (VTI), while the carbonate shows 4–5% azimuthal anisotropy in traveltimes. The fast direction is consistent with the maximum horizontal stress orientation determined from break‐out logs and is also consistent with the strike of major faults. AVO analysis of the reflection from the top of the carbonate layer shows a critical angle reduction in the fast direction and maximum gradient in the slow direction. This agrees with modeling and indicates a greater amplitude sensitivity in the slow direction—the direction perpendicular to fracture strike. In principle, 3-D surveys should have wide azimuthal coverage to characterize fractured reservoirs. If this is not possible, it is important to have azimuthal line coverage in the minimum horizontal stress direction to optimize the use of AVO for fractured reservoir characterization. This direction can be obtained from multiazimuthal walkaways using the azimuthal P-wave analysis techniques presented.


1995 ◽  
Author(s):  
Leon Thomsen ◽  
Ilya Tsvankin ◽  
Michael C. Mueller

Geophysics ◽  
1990 ◽  
Vol 55 (7) ◽  
pp. 948-948

Figure 12 in our paper is not complete. Unfortunately, the top part was cut off and all the words do not appear. The entire illustration is shown here with caption.


Geophysics ◽  
2020 ◽  
Vol 85 (4) ◽  
pp. D121-D132
Author(s):  
Yang-Hu Li ◽  
Song Xu ◽  
Can Jiang ◽  
Yuan-Da Su ◽  
Xiao-Ming Tang

Seismic-wave anisotropy has long been an important topic in the exploration and development of unconventional reservoirs, especially in shales, which are commonly characterized as transversely isotropic ([TI] or vertical TI [VTI]) media. At present, the shear-wave (S-wave) TI properties have mainly been determined from monopole Stoneley- or dipole flexural-wave measurements in wireline acoustic logging, but the feasibility of those obtained from logging-while-drilling (LWD) acoustic data needs to be established. We have developed a joint inversion method for simultaneously determining formation S-wave transverse isotropy and vertical velocity from LWD multipole acoustic data. Our theoretical analysis shows that the presence of anisotropy strongly influences LWD Stoneley- and quadrupole-wave dispersion characteristics. Although the monopole Stoneley and quadrupole waves are sensitive to the formation S-wave TI parameters, they suffer from the typical nonuniqueness problem when using the individual-wave data to invert parameters alone. Thus, the respective dispersion data can be jointly used to estimate the formation S-wave TI properties. By the joint inversion, the nonuniqueness problem in the parameter inversion can also be effectively alleviated. The feasibility of the method has been verified by the processing results of theoretical synthetic data and field LWD acoustic-wave data. Therefore, the result offers an effective method for evaluating VTI formation anisotropy from acoustic LWD data.


Author(s):  
Babak Sayad Noghretab ◽  
Mohammad Kamal Ghassem-Alaskari

The purpose of this article was to generate and compare seismic modeling results with real vertical seismic profiling data (VSP data) based on Dipole Shear Imager (DSI) data in the reservoir zone (Kangan and upper Dalan Formations) of a well in South Pars gas field. Estimation of shear wave velocity (Vs) and density for layers above the reservoir zone, for which; DSI data did not exist, was also done by the applied modeling method to estimate elastic parameters of the layers. In this method, modeling for X-component of the VSP survey was run by utilizing the DSI data set of reservoir zone and the VSP survey report of the studied well with high precision. Computed results for the proposed modeling method led to achieving highly accurate, close to the reality of VSP model around the studied well. According to compression wave velocity (VP) attained from VSP survey reports of the well and Vp/Vs ratio obtained from Dipole Shear Imager (DSI), modeling was done. Afterward, shear wave velocity (Vs) for upper layers of reservoir zone estimated with high precision, then density and elastic moduli for the above layers and the reservoir zone were calculated.


Sign in / Sign up

Export Citation Format

Share Document