From Anisotropic Elastic Tensors to Seismic Data and Back: Extraction of Azimuthal Anisotropy Parameters from a Field Sc

Author(s):  
P. Tillotson ◽  
D. Davies ◽  
L. Wang ◽  
M. Ball ◽  
A. Toomey ◽  
...  
Geophysics ◽  
2019 ◽  
Vol 84 (5) ◽  
pp. C229-C237 ◽  
Author(s):  
Shibo Xu ◽  
Alexey Stovas

The moveout approximations are commonly used in seismic data processing such as velocity analysis, modeling, and time migration. The anisotropic effect is very obvious for a converted wave when estimating the physical and processing parameters from the real data. To approximate the traveltime in an elastic orthorhombic (ORT) medium, we defined an explicit rational-form approximation for the traveltime of the converted [Formula: see text]-, [Formula: see text]-, and [Formula: see text]-waves. To obtain the expression of the coefficients, the Taylor-series approximation is applied in the corresponding vertical slowness for three pure-wave modes. By using the effective model parameters for [Formula: see text]-, [Formula: see text]-, and [Formula: see text]-waves, the coefficients in the converted-wave traveltime approximation can be represented by the anisotropy parameters defined in the elastic ORT model. The accuracy in the converted-wave traveltime for three ORT models is illustrated in numerical examples. One can see from the results that, for converted [Formula: see text]- and [Formula: see text]-waves, our rational-form approximation is very accurate regardless of the tested ORT model. For a converted [Formula: see text]-wave, due to the existence of cusps, triplications, and shear singularities, the error is relatively larger compared with PS-waves.


2019 ◽  
Vol 38 (2) ◽  
pp. 144-150 ◽  
Author(s):  
Marianne Rauch-Davies ◽  
David Langton ◽  
Michael Bradshaw ◽  
Allon Bartana ◽  
Dan Kosloff ◽  
...  

With readily available wide-azimuth, onshore, 3D seismic data, the search for attributes utilizing the azimuthal information is ongoing. Theoretically, in the presence of ordered fracturing, the seismic wavefront shape changes from spherical to nonspherical with the propagation velocity being faster parallel to the fracturing and slower perpendicular to the fracture direction. This concept has been adopted and is used to map fracture direction and density within unconventional reservoirs. More specifically, azimuthal variations in normal moveout velocity or migration velocity are often used to infer natural fracture orientation. Analyses of recent results have called into question whether azimuthal velocity linked to intrinsic azimuthal velocity variations can actually be detected from seismic data. By use of 3D orthorhombic anisotropic elastic simulation, we test whether fracture orientation and intensity can be detected from seismic data. We construct two subsurface models based on interpreted subsurface layer structure of the Anadarko Basin in Oklahoma. For the first model, the material parameters in the layers are constant vertically transverse isotropic (VTI) in all intervals. The second model was constructed the same way as the base model for all layers above the Woodford Shale Formation. For the shale layer, orthorhombic properties were introduced. In addition, a thicker wedge layer was added below the shale layer. Using the constructed model, synthetic seismic data were produced by means of 3D anisotropic elastic simulation resulting in two data sets: VTI and orthorhombic. The simulated data set was depth migrated using the VTI subsurface model. After migration, the residual moveouts on the migrated gathers were analyzed. The analysis of the depth-migrated model data indicates that for the typical layer thicknesses of the Woodford Shale layer in the Anadarko Basin, observed and modeled percentage of anisotropy and target depth, the effect of intrinsic anisotropy is too small to be detected in real seismic data.


Geophysics ◽  
1998 ◽  
Vol 63 (2) ◽  
pp. 692-706 ◽  
Author(s):  
Subhashis Mallick ◽  
Kenneth L. Craft ◽  
Laurent J. Meister ◽  
Ronald E. Chambers

In an azimuthally anisotropic medium, the principal directions of azimuthal anisotropy are the directions along which the quasi-P- and the quasi-S-waves propagate as pure P and S modes. When azimuthal anisotropy is induced by oriented vertical fractures imposed on an azimuthally isotropic background, two of these principal directions correspond to the directions parallel and perpendicular to the fractures. S-waves propagating through an azimuthally anisotropic medium are sensitive to the direction of their propagation with respect to the principal directions. As a result, primary or mode‐converted multicomponent S-wave data are used to obtain the principal directions. Apart from high acquisition cost, processing and interpretation of multicomponent data require a technology that the seismic industry has not fully developed. Anisotropy detection from conventional P-wave data, on the other hand, has been limited to a few qualitative studies of the amplitude variation with offset (AVO) for different azimuthal directions. To quantify the azimuthal AVO, we studied the amplitude variation with azimuth for P-wave data at fixed offsets. Our results show that such amplitude variation with azimuth is periodic in 2θ, θ being the orientation of the shooting direction with respect to one of the principal directions. For fracture‐induced anisotropy, this principal direction corresponds to the direction parallel or perpendicular to the fractures. We use this periodic azimuthal dependence of P-wave reflection amplitudes to identify two distinct cases of anisotropy detection. The first case is an exactly determined one, where we have observations from three azimuthal lines for every common‐midpoint (CMP) location. We derive equations to compute the orientation of the principal directions for such a case. The second case is an overdetermined one where we have observations from more than three azimuthal lines. Orientation of the principal direction from such an overdetermined case can be obtained from a least‐squares fit to the reflection amplitudes over all the azimuthal directions or by solving many exactly determined problems. In addition to the orientation angle, a qualitative measure of the degree of azimuthal anisotropy can also be obtained from either of the above two cases. When azimuthal anisotropy is induced by oriented vertical fractures, this qualitative measure of anisotropy is proportional to fracture density. Using synthetic seismograms, we demonstrate the robustness of our method in evaluating the principal directions from conventional P-wave seismic data. We also apply our technique to real P-wave data, collected over a wide source‐to‐receiver azimuth distribution. Computations using our method gave an orientation of the principal direction consistent with the general fracture orientation in the area as inferred from other geological and geophysical evidence.


Geophysics ◽  
2001 ◽  
Vol 66 (2) ◽  
pp. 582-597 ◽  
Author(s):  
Donald F. Winterstein ◽  
Gopa S. De ◽  
Mark A. Meadows

Since 1986, when industry scientists first publicly showed data supporting the presence of azimuthal anisotropy in sedimentary rock, we have studied vertical shear‐wave (S-wave) birefringence in 23 different wells in western North America. The data were from nine‐component vertical seismic profiles (VSPs) supplemented in recent years with data from wireline crossed‐dipole logs. This paper summarizes our results, including birefringence results in tabular form for 54 depth intervals in 19 of those 23 wells. In the Appendix we present our conclusions about how to record VSP data optimally for study of vertical birefringence. We arrived at four principal conclusions about vertical S-wave birefringence. First, birefringence was common but not universal. Second, birefringence ranged from 0–21%, but values larger than 4% occurred only in shallow formations (<1200 m) within 40 km of California’s San Andreas fault. Third, at large scales birefringence tended to be blocky. That is, both the birefringence magnitude and the S-wave polarization azimuth were often consistent over depth intervals of several tens to hundreds of meters but then changed abruptly, sometimes by large amounts. Birefringence in some instances diminished with depth and in others increased with depth, but in almost every case a layer near the surface was more birefringent than the layer immediately below it. Fourth, observed birefringence patterns generally do not encourage use of multicomponent surface reflection seismic data for finding fractured hydrocarbon reservoirs, but they do encourage use of crossed‐dipole logs to examine them. That is, most reservoirs were birefringent, but none we studied showed increased birefringence confined to the reservoir.


Sign in / Sign up

Export Citation Format

Share Document