Detection of fracture orientation using azimuthal variation of P-wave AVO responses

Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1253-1265 ◽  
Author(s):  
Maria A. Pérez ◽  
Richard L. Gibson ◽  
M. Nafi Toksöz
Keyword(s):  
Oil And Gas ◽  
Seismic Wave Propagation ◽  
P Wave ◽  
Fractured Reservoir ◽  
P Waves ◽  
Data Set ◽  
Amplitude Variation With Offset ◽  
Azimuthal Variation ◽  
Fracture Orientation ◽  
Converted Waves

Azimuthally dependent P-wave amplitude variation with offset (AVO) responses can be related theoretically to open fracture orientation and have been suggested as a geophysical tool to identify fracture orientation in fractured oil and gas reservoirs. A field experiment conducted recently over a fractured reservoir in the Barinas Basin, Venezuela provides data for an excellent test of this approach. Three lines of data were collected in three different azimuths, and three component receivers were used. The distribution of fractures in this reservoir was obtained previously using measurements of shear‐wave splitting from P-S converted waves from the same data set. In this work, we use P-wave data to see if the data can yield the same information using azimuthal variation of P-wave AVO responses. Results obtained from the azimuthal P-wave AVO analysis corroborate the fracture orientation obtained previously using P-S converted waves. This analysis with field data is an example of the high potential of P-waves to detect fracture effects on seismic wave propagation.

Basin-wide rock-physics analysis in Campeche Basin, Gulf of Mexico — Phase II: Reservoir rock and fluid properties

2021 ◽  
Vol 40 (10) ◽  
pp. 716-722
Author(s):  
Yangjun (Kevin) Liu ◽  
Michelle Ellis ◽  
Mohamed El-Toukhy ◽  
Jonathan Hernandez
Keyword(s):  
Oil And Gas ◽  
Rock Physics ◽  
Clay Content ◽  
American Petroleum Institute ◽  
P Wave ◽  
Reservoir Rock ◽  
Gas Oil ◽  
Amplitude Variation With Offset ◽  
Fluid Properties ◽  
The Impact

We present a basin-wide rock-physics analysis of reservoir rocks and fluid properties in Campeche Basin. Reservoir data from discovery wells are analyzed in terms of their relationship between P-wave velocity, density, porosity, clay content, Poisson's ratio (PR), and P-impedance (IP). The fluid properties are computed by using in-situ pressure, temperature, American Petroleum Institute gravity, gas-oil ratio, and volume of gas, oil, and water. Oil- and gas-saturated reservoir sands show strong PR anomalies compared to modeled water sand at equivalent depth. This suggests that PR anomalies can be used as a direct hydrocarbon indicator in the Tertiary sands in Campeche Basin. However, false PR anomalies due to residual gas or oil exist and compose about 30% of the total anomalies. The impact of fluid properties on IP and PR is calibrated using more than 30 discovery wells. These calibrated relationships between fluid properties and PR can be used to guide or constrain amplitude variation with offset inversion for better pore fluid discrimination.

Amplitude variation with offset inversion using the reflectivity method

Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. R185-R195 ◽  
Author(s):  
Hongxing Liu ◽  
Jingye Li ◽  
Xiaohong Chen ◽  
Bo Hou ◽  
Li Chen
Keyword(s):  
Wave Velocity ◽  
P Wave ◽  
Transmission Losses ◽  
Amplitude Variation ◽  
Data Set ◽  
Amplitude Variation With Offset ◽  
Avo Inversion ◽  
Propagation Effects ◽  
Reflectivity Method ◽  
Internal Multiples

Most existing amplitude variation with offset (AVO) inversion methods are based on the Zoeppritz’s equation or its approximations. These methods assume that the amplitude of seismic data depends only on the reflection coefficients, which means that the wave-propagation effects, such as geometric spreading, attenuation, transmission loss, and multiples, have been fully corrected or attenuated before inversion. However, these requirements are very strict and can hardly be satisfied. Under a 1D assumption, reflectivity-method-based inversions are able to handle transmission losses and internal multiples. Applications of these inversions, however, are still time-consuming and complex in computation of differential seismograms. We have evaluated an inversion methodology based on the vectorized reflectivity method, in which the differential seismograms can be calculated from analytical expressions. It is computationally efficient. A modification is implemented to transform the inversion from the intercept time and ray-parameter domain to the angle-gather domain. AVO inversion is always an ill-posed problem. Following a Bayesian approach, the inversion is stabilized by including the correlation of the P-wave velocity, S-wave velocity, and density. Comparing reflectivity-method-based inversion with Zoeppritz-based inversion on a synthetic data and a real data set, we have concluded that reflectivity-method-based inversion is more accurate when the propagation effects of transmission losses and internal multiples are not corrected. Model testing has revealed that the method is robust at high noise levels.

Analysis of multiazimuthal VSP data for anisotropy and AVO

Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1172-1180 ◽  
Author(s):  
W. Scott Leaney ◽  
Colin M. Sayers ◽  
Douglas E. Miller
Keyword(s):  
Fractured Reservoirs ◽  
Vertical Axis ◽  
Horizontal Stress ◽  
Azimuthal Anisotropy ◽  
Carbonate Reservoir ◽  
P Wave ◽  
Fractured Reservoir ◽  
Data Set ◽  
Vsp Data ◽  
The Impact

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.

Joint inversion of P- and PS-waves in orthorhombic media: Theory and a physical modeling study

Geophysics ◽  
1999 ◽  
Vol 64 (1) ◽  
pp. 146-161 ◽  
Author(s):  
Vladimir Grechka ◽  
Stephen Theophanis ◽  
Ilya Tsvankin
Keyword(s):  
Shear Wave ◽  
Physical Modeling ◽  
Horizontal Plane ◽  
Joint Inversion ◽  
Symmetry Plane ◽  
P Wave ◽  
Simple Relationship ◽  
Converted Wave ◽  
Data Set ◽  
Converted Waves

Reflection traveltimes recorded over azimuthally anisotropic fractured media can provide valuable information for reservoir characterization. As recently shown by Grechka and Tsvankin, normal moveout (NMO) velocity of any pure (unconverted) mode depends on only three medium parameters and usually has an elliptical shape in the horizontal plane. Because of the limited information contained in the NMO ellipse of P-waves, it is advantageous to use moveout velocities of shear or converted modes in attempts to resolve the coefficients of realistic orthorhombic or lower‐symmetry fractured models. Joint inversion of P and PS traveltimes is especially attractive because it does not require shear‐wave excitation. Here, we show that for models composed of horizontal layers with a horizontal symmetry plane, the traveltime of converted waves is reciprocal with respect to the source and receiver positions (i.e., it remains the same if we interchange the source and receiver) and can be adequately described by NMO velocity on conventional‐length spreads. The azimuthal dependence of converted‐wave NMO velocity has the same form as for pure modes but requires the spatial derivatives of two-way traveltime for its determination. Using the generalized Dix equation of Grechka, Tsvankin, and Cohen, we derive a simple relationship between the NMO ellipses of pure and converted waves that provides a basis for obtaining shear‐wave information from P and PS data. For orthorhombic models, the combination of the reflection traveltimes of the P-wave and two split PS-waves makes it possible to reconstruct the azimuthally dependent NMO velocities of the pure shear modes and to find the anisotropic parameters that cannot be determined from P-wave data alone. The method is applied to a physical modeling data set acquired over a block of orthorhombic material—Phenolite XX-324. The inversion of conventional‐spread P and PS moveout data allowed us to obtain the orientation of the vertical symmetry planes and eight (out of nine) elastic parameters of the medium (the reflector depth was known). The remaining coefficient (c12 or δ(3) in Tsvankin’s notation) is found from the direct P-wave arrival in the horizontal plane. The inversion results accurately predict moveout curves of the pure S-waves and are in excellent agreement with direct measurements of the horizontal velocities.

AVO inversion using pseudoisotropic elastic properties

2021 ◽  
Vol 40 (1) ◽  
pp. 52-59
Author(s):  
Michinori Asaka
Keyword(s):  
Elastic Properties ◽  
Oil And Gas ◽  
Reservoir Characterization ◽  
Elastic Anisotropy ◽  
Data Set ◽  
Amplitude Variation With Offset ◽  
Avo Inversion ◽  
Characterization Study ◽  
Anisotropy Parameters ◽  
Using Data

Amplitude variation with offset (AVO) inversion of an anisotropic data set is a challenging task. Nonnegligible differences in the anisotropy parameters between the various lithologies make the seismic data AVO response completely different from the isotropic synthetic seismogram. In this case, it is difficult to invert for VP/VS and density consistent with well-log data. AVO inversion using pseudoisotropic elastic properties is a practical solution to this problem. Verification of this method was performed using data from an offshore Western Australia field. It was found that wavelet extraction and density inversion are improved significantly by replacing the isotropic elastic properties with the pseudoisotropic properties. Inverted density shows reasonable quality and therefore can be included in the reservoir characterization study. Postinversion analyses can be performed effectively on the pseudoisotropic elastic properties because crossplot analysis shows the increased separation of different lithofacies due to contrasts in anisotropy parameters. This result could have significant implications for other fields, as shale constitutes most of the overburden in conventional oil and gas fields and often shows strong elastic anisotropy.

AVO and AVA inversion for fractured reservoir characterization

Geophysics ◽  
2002 ◽  
Vol 67 (1) ◽  
pp. 300-306 ◽  
Author(s):  
Matteo Mario Beretta ◽  
Giancarlo Bernasconi ◽  
Giuseppe Drufuca
Keyword(s):  
Reservoir Characterization ◽  
Spectral Component ◽  
Plane Waves ◽  
Synthetic Data ◽  
Transversely Isotropic ◽  
Fractured Reservoir ◽  
Fracture Density ◽  
P Waves ◽  
Data Set ◽  
Inversion Procedure

Seismic wave reflection amplitudes are used to detect fluids and fracture properties in reservoirs. This paper studies the characterization of a vertically fractured fluid‐filled reservoir by analyzing the reflection amplitudes of P‐waves with varying incident and azimuthal angles. The reservoir is modeled as a horizontal transversely isotropic medium embedded in an isotropic background, and the linearized P‐waves reflection coefficient are considered. The conditioning of the inverse problem is analyzed, and fracture density is found to be the best conditioned parameter. Using diffraction tomography under the Born approximation, an inversion procedure is proposed in the transformed k–ω domain to detect fracture density variations within the reservoir. Seismic data are rearranged in pairs of incident and reflected plane waves, enlightening only one spectral component of the fracture density field at a time. Only the observable spectral components are inverted. Moreover, working in the transformed domain, picking reflection amplitudes is not required. An example of the inversion applied to a synthetic data set is presented. The limitation of source and receiver numbers and the finite bandwidth of the wavelet produce a loss of resolution, but the overall fracture density variations are recovered correctly.

Velocity model building for a single-offset 3C VSP data via deformable-layer tomography: a Texas salt dome example

Geophysics ◽  
2021 ◽  
pp. 1-44
Author(s):  
Yukai Wo ◽  
Jingjing Zong ◽  
Hao Hu ◽  
Hua-Wei Zhou ◽  
Robert R. Stewart
Keyword(s):  
Velocity Model ◽  
Salt Dome ◽  
P Wave ◽  
Velocity Variation ◽  
Shear Mode ◽  
P Waves ◽  
Data Set ◽  
Image Domain ◽  
Vsp Data ◽  
Geometric Variation

We have applied multiscale deformable-layer tomography (DLT) to build a laterally varying velocity model, using a single-offset vertical seismic profile (VSP) data set acquired for a salt proximity survey in southern Texas. The purpose of the VSP survey is to delineate the 2D salt flank using the P-wave reflections. Previous study has identified an anhydrate layer as the cap rock of the salt dome. The large impedance contrasts of this anhydrite layer generate strong downgoing P (sediment)-S (anhydrite)-P (salt) waves recorded by downhole geophones. Incidentally, the P-S-P-waves have similar traveltimes as those of the P-wave salt flank reflections, thus contaminating the imaging of the salt flank. Identifying shear-mode contamination requires an accurate velocity model of anhydrite. However, the extremely poor coverage of the single-offset VSP greatly challenges tomographic techniques to determine the lateral velocity variation. We tackle this problem using multiscale DLT, which characterizes the velocity field by a set of deformable layers. We constrain the layer velocities using the check-shot data and invert for the geometric variation. The inverted model indicates that the anhydrite layer has a “thick-thin-thick” lateral variation with offset, and the S-wave in the anhydrite layer helps in imaging the P-S-P-waves along the well track. The estimated anhydrite layer geometry is validated by the kinematic accuracies of P-waves in the data domain and P-S-P-waves in the image domain. Some in-salt dipping structures are determined by multiscale DLT as well. This field data example indicates that multiscale DLT is feasible for estimating velocities using VSP data of the single-offset situation. An accurate velocity model is the key for modeling and adaptive subtraction of the shear-mode contamination related to the salt geometry.

3-D moveout inversion in azimuthally anisotropic media with lateral velocity variation: Theory and a case study

Geophysics ◽  
1999 ◽  
Vol 64 (4) ◽  
pp. 1202-1218 ◽  
Author(s):  
Vladimir Grechka ◽  
Ilya Tsvankin
Keyword(s):  
Correction Term ◽  
Anisotropic Media ◽  
Azimuthal Anisotropy ◽  
P Wave ◽  
Velocity Variation ◽  
Variation Theory ◽  
Fractured Reservoir ◽  
Lateral Velocity ◽  
Data Set ◽  
Zero Offset

Reflection moveout recorded over an azimuthally anisotropic medium (e.g., caused by vertical or dipping fractures) varies with the azimuth of the source‐receiver line. Normal‐moveout (NMO) velocity, responsible for the reflection traveltimes on conventional‐length spreads, forms an elliptical curve in the horizontal plane. While this result remains valid in the presence of arbitrary anisotropy and heterogeneity, the inversion of the NMO ellipse for the medium parameters has been discussed so far only for horizontally homogeneous models above a horizontal or dipping reflector. Here, we develop an analytic moveout correction for weak lateral velocity variation in horizontally layered azimuthally anisotropic media. The correction term is proportional to the curvature of the zero‐offset traveltime surface at the common midpoint and, therefore, can be estimated from surface seismic data. After the influence of lateral velocity variation on the effective NMO ellipses has been stripped, the generalized Dix equation can be used to compute the interval ellipses and evaluate the magnitude of azimuthal anisotropy (measured by P-wave NMO velocity) within the layer of interest. This methodology was applied to a 3-D “wide‐azimuth” data set acquired over a fractured reservoir in the Powder River Basin, Wyoming. The processing sequence included 3-D semblance analysis (based on the elliptical NMO equation) for a grid of common‐midpoint “supergathers,” spatial smoothing of the effective NMO ellipses and zero‐offset traveltimes, correction for lateral velocity variation, and generalized Dix differentiation. Our estimates of depth‐varying fracture trends in the survey area, based on the interval P-wave NMO ellipses, are in good agreement with the results of outcrop and borehole measurements and the rotational analysis of four‐ component S-wave data.

Seismic vertical transversely isotropic parameter inversion from P- and S-wave cross-borehole measurements in an aquifer environment

Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. D101-D116
Author(s):  
Julius K. von Ketelhodt ◽  
Musa S. D. Manzi ◽  
Raymond J. Durrheim ◽  
Thomas Fechner
Keyword(s):  
Transversely Isotropic ◽  
P Wave ◽  
Perturbation Equation ◽  
S Wave ◽  
P Waves ◽  
Data Set ◽  
Near Surface ◽  
S Waves ◽  
Wave Measurements ◽  
Wave Splitting

Joint P- and S-wave measurements for tomographic cross-borehole analysis can offer more reliable interpretational insight concerning lithologic and geotechnical parameter variations compared with P-wave measurements on their own. However, anisotropy can have a large influence on S-wave measurements, with the S-wave splitting into two modes. We have developed an inversion for parameters of transversely isotropic with a vertical symmetry axis (VTI) media. Our inversion is based on the traveltime perturbation equation, using cross-gradient constraints to ensure structural similarity for the resulting VTI parameters. We first determine the inversion on a synthetic data set consisting of P-waves and vertically and horizontally polarized S-waves. Subsequently, we evaluate inversion results for a data set comprising jointly measured P-waves and vertically and horizontally polarized S-waves that were acquired in a near-surface ([Formula: see text]) aquifer environment (the Safira research site, Germany). The inverted models indicate that the anisotropy parameters [Formula: see text] and [Formula: see text] are close to zero, with no P-wave anisotropy present. A high [Formula: see text] ratio of up to nine causes considerable SV-wave anisotropy despite the low magnitudes for [Formula: see text] and [Formula: see text]. The SH-wave anisotropy parameter [Formula: see text] is estimated to be between 0.05 and 0.15 in the clay and lignite seams. The S-wave splitting is confirmed by polarization analysis prior to the inversion. The results suggest that S-wave anisotropy may be more severe than P-wave anisotropy in near-surface environments and should be taken into account when interpreting cross-borehole S-wave data.

Sign in / Sign up

Export Citation Format

Share Document