Sedimentary rocks especially shales show anisotropy if seismic waves pass through them at different directions. The presence of anisotropy has considerable influence on many aspects of seismic wave propagation and therefore has profound implications for both conventional processing schemes and interpretation. When the subsurface is anisotropic, the moveout of P-waves in horizontally layered media as observed in common midpoint gathers, deviates from being hyperbolic at far offsets. Thus the nonhyperbolic component of the moveout is routinely used for estimating relevant seismic parameters describing such anisotropy (Grechka and Tsvankin, 1998). For P-waves in Vertical Transverse Isotropic (VTI) media, these parameters are the anellipticity parameter (η) and the normal moveout velocity (Vnmo). Quantitative measurement of seismic anisotropy provides a valuable clue to lithology and degree of stratification, a tool for imaging and can as well provide important quantitative information about depth discrepancies observed from depths determined from seismic and well log data (Van der Baan and Kendall, 2002).
Anisotropic properties are estimated using a variety of ways. Laboratory tests
are the main techniques used in measuring anisotropic properties. Vertical Seismic
Profiling (VSP) and surface seismic data are the in situ techniques used
to determine anisotropic parameters (Li, 2006). Using surface seismic data,
the most popular approach is using modified Taylor series expressions (Alkhalifah,
1997) or tau-p (τ-ρ) based methods (Van der Baan and Kendall, 2002).
Recently, Douma and Calvert (2006) used a semblance based rational interpolation
method to analyse non-hyperbolic moveout in the time-offset (t-x) domain to
estimate anisotropy parameters and moveout velocity. These workers pointed out
that the picking of traveltimes in the t-x domain which is often laborious and
prone to error in the τ-ρ technique of Van der Baan and Kendall (2002)
is a practical disadvantage. To overcome this caveat, Douma and Van der Baan
(2006) combined the τ-ρ technique of Van der Baan and Kendall (2002)
and the rational interpolation method of Douma and Calvert (2006) to estimate
In this study, we use the approach of Douma and Van der Baan (2006) to determine the lateral variations of both effective and interval values of the anellipticity parameter (η), the NMO (vnmo) and horizontal (vhor) velocities at selected CDP locations in a marine 2D seismic dataset acquired in the Gulf of Mexico.
The τ-p transform relates traveltime t and offset x to intercept time τ and horizontal slowness p by:
Hake (1986) observed that the total intercept time τ for p-wave reflections in horizontally layered transverse isotropy with vertical axis of symmetry (VTI) media is a summation of the products of layer interval zero-offset time (Δt0,i), the vertical P-wave velocity (v0,i) and the vertical slowness (qi) expressed as:
Alkhalifah (1998) expressed the vertical slowness qi in terms of vnmo and the anellipticity parameter η by assuming that the interval shear-wave phase velocity Vso has negligible influence on the traveltimes of P-waves in VTI media. Alkhalifahs equation is expressed as:
||The horizontal component of the slowness vector
Grechka and Tsvankin (1998) however demonstrated that the horizontal velocity expressed as Vhor = Vnmo [1+2η]1/2 is better suited for semblance-based moveout analysis than η as used in Alkhalifahs equation. Thus, expressing Eq. 3 in terms of Vhor and Vnmo and substituting into 2, gives:
||Horizontal velocity in layer i
||NMO velocity in layer i
For semblance analysis to be carried out in the t-x domain requires the traveltime for each recorded offset. From Eq. 1,
Substituting Eq. 4 into 5 and simplifying,
The traveltime associated with this offset was obtained by combining Eq.
1, 2 and 6 giving:
NONHYPERBOLIC MOVEOUT ANALYSIS USING RATIONAL INTERPOLATION
A rational approximation to a function T(x) is generally written as Stoer and Bulirsch (1993)
where, NL(x) and DM(x) are polynomials of degree L and M. Such a rational approximation can be written as [L/M]. For a polynomial of degree 2, L = M = 2. Thus, Eq. 8 can be rewritten as:
For a [2/2] rational interpolation for a single horizontal VTI layer, Douma and Calvert (2006) derived the squared traveltimes T as a function of squared offset X as:
where, T0 = t20 is the squared zero-offset
two-way traveltime and n1, n1, d1 and d2
are the coefficients of the numerator and denominator of the rational approximant.
To determine the unknowns, four traveltimes (ti) and four associated
offsets (xi) are required. The traveltimes (ti) and the
associated offset (xi) are determined using Eq. 7
and 11, respectively.
||The offset-to-depth ratios.
The interpolation traveltimes (ti) and (xi) are used in determining the coefficients of Eq. 10. When once the coefficients of Eq. 10 are known, the traveltimes (t) for offset (x) can be evaluated. Douma and Van der Baan (2006) extended this approach to horizontally layered VTI media.
The data used in this study is a 2D marine seismic P-wave data set acquired
in the Gulf of Mexico. The data was processed by applying a bandpass filter
of frequency range 1-5-35-70 Hz to remove frequencies outside of the useful
signal bandwidth. Common Midpoint gathers were then collected into groups of
nine and stacked (vertically) to form supergathers. Frequency-wavenumber (F-K)
filtering was then applied to each supergather to suppress coherent noise and
semblance analysis performed in order to obtain stacking velocities for use
in normal moveout (NMO) correction. Common noise observed after F-K filtering
is random noise. To reduce the random noise further, frequency-offset (F-X)
deconvolution (Canales, 1984) was applied. Finally, a coherency filter was applied
to improve signal character and continuity within the frequency range of the
ANISOTROPY PARAMETER ESTIMATION
Ten supergathers (38559-38640) between CDPs 37999 and 38941 were chosen for parameter estimation. These ten supergathers hereinafter are referred to as ensembles 1-10. This zone was chosen because the stack section shows horizons that are relatively flat. Thus, it was assumed that the presence of any nonhyperbolic moveout will be due to anisotropy. In each supergather, three events were chosen for analysis, these events occur at 1.14, 1.44 and 2.4 sec (Fig. 1). The algorithm SUVEL2DF used in this study determines values of the effective and interval anellipticity parameter (η), horizontal velocity (vhmax) and NMO velocity (vmax) through semblance analysis. Semblance analysis using SUVEL2DF requires the computation of the traveltimes for each offset acquired in the field for a specific combination of zero-offset time (t0), NMO (vnmo) and horizontal (vhor) velocities.
The most important input parameters are the zero-offset time (t0)
of the layer of interest determined from the supergather gather and user specified
effective NMO (vnmo), horizontal (vhor) velocities and
four offset-to-depth ratios ranging from 1-4. For each particular vhor
and vnmo input into the algorithm, the four offset-to-depth ratios are
converted to target offsets (xj) using Eq. 11.
||Ensemble 38604 showing reflection events chosen for parameter
estimation in each supergather
||Semblance scan for layer 1 (ensemble 6) showing optimal semblance
The algorithm also computes the appropriate horizontal slowness (pj)
to target offsets (xj) using Eq. 6 and subsequently
determine the corresponding traveltimes (tj) by the use of Eq.
7. The four values of the target offsets (xj) and traveltimes
(tj) constitute the support points for the [2/2] rational interpolation.
The algorithm then determines the desired traveltimes t(x) at any offset x using
the rational approximant (Eq. 10). The resulting t-x curves
are then used for the semblance analysis.
We determined both effective and interval values of η, vhor and vnmo. In determining effective values, we treated the overburden of the reflection event of interest as a single horizontal homogeneous VTI layer. As result, the effective values are average values. Interval values were obtained using a layer-stripping approach. The accuracy of the inversion was confirmed by the high semblance maxima (Smax≈0.8-1.0) which coincided with the semblance contours nicely centred in the respective semblance scans (Fig. 2).
RESULTS AND DISCUSSION
Effective values of vnmo and vhor range between 1850-2250±111 and 1975-2475±122 m sec-1 while interval values of vnmo and vhor range between 1750-2650±120 and 2100-3100±126 m sec-1, respectively. Values of the anellipticity parameter (η) as a function of ensemble number are shown in Fig. 3 and 4, respectively. The effective values are mainly between 0.03-0.16 ±0.01 in all the three layers investigated. Interval values range mainly between 0.03-0.2±0.03. Variation in anellipticity is observed in both effective and interval values suggesting lateral variation in anisotropy along the seismic line presumably due to varying shale content since the presence of shale is widely accepted as the cause of anisotropy (Wang, 2002; Li, 2006). The recovered values of vhor and vnmo exhibit a similar behaviour. In layer 2 in Fig. 3 (ensemble 1) and layer 2 in Fig. 4 (ensembles 4 and 8), we observe an abrupt increase in anellipticity values where the anellipticity is as high as 0.3 and 0.4, respectively. These values of η are high compared to values typically reported from nonhyperbolic moveout analysis (Toldi et al., 1999).
Without detailed geologic information, we can only speculate as to the cause
of these localised extreme values. The relatively high η values might reflect
some internal property of the medium, presumably increased shale content or
that these are inaccurate values due to the inability to obtain optimal semblance
response occasioned by the use of inaccurate input model parameters. If the
first assumption was correct, the nonhyperbolic moveout in the input CDP gathers
in these events should have been substantial.
anellipticity parameter (η) against ensemble number for all three
layers showing lateral the lateral variation of anisotropy at ensemble
anellipticity parameter (η) against ensemble number. Layer 1 is not
shown because in layer 1, effective and interval values are the same
We observe that the nonhyperbolic moveout in these events was rather small
and therefore might not be able to cause such a high level of anellipticity,
since the strength of the anisotropy determines the amount of non-hyperbolic
moveout (Wookey et al., 2002). On the other hand, a careful look at the
semblance scans of these events show more than one semblance contour making
it difficult to obtain optimum semblance response. Since correct anellipticity
values are only obtained at optimum semblance response, we are convinced that
the anellipticity values of 0.3 observed in ensemble 1, layer 2 (Fig.
3) and ~0.4 in ensembles 4 and 8, layer 2 (Fig. 4) are
not correct and do not represent the actual subsurface anellipticity level at
these ensemble locations. The presence of more than one semblance contour might
be due to the presence of multiples.
Results from this study show that the three layers investigated have considerable variations in lithology. The effective and interval values of the anellipticity parameter vary mainly between 0.03-0.2±0.01. Effective and interval values of vnmo and vhor range between 1975-2825±111 and 1750-2650±126 m sec-1, respectively. We attribute the variation in anisotropy to varying shale content. These results can be employed as new constraints in processing algorithms that include anisotropy so as to improve the imaging of dipping reflectors, such as fault planes in the study area.