 |
 |
 |
 | A fast algorithm for 3D azimuthally anisotropic velocity scan |  |
![[pdf]](icons/pdf.png) |
Next: Basic formulation
Up: Hu et al.: Fast
Previous: Introduction
As explained by Grechka and Tsvankin (1998), a pure-mode (P or S) reflection event in an effectively azimuthally anisotropic medium can be described by
 |
(1) |
where
is two-way CMP traveltime,
is two-way zero-offset traveltime,
is the full source-receiver offset in surface survey coordinates, and
 |
(2) |
is the slowness matrix. Equation 1 follows from a truncated 2D Taylor expansion. Geometrically, it represents a curved surface that is hyperbolic in cross section and elliptic in map view.
Ideally, one can perform a semblance scan (Taner and Koehler, 1969) over the three parameters
,
, and
simultaneously to estimate the velocity and perform NMO correction. However, this approach, if not impossible, is extremely expensive for large-size seismic data. Furthermore, since these parameters are not orthogonal, the semblance plots might appear to be extended and ambiguous, hence presenting difficulties for picking (Fowler et al., 2006).
Davidson et al. (2011) proposed a stable way of detecting azimuthal anisotropy using an orthogonal parametrization of the moveout function, which is based on an equivalent reformulation of equation 1,
 |
(3) |
The cosine and sine dependent slownesses
and
are usually much smaller than the averaged slowness
. Therefore, a possible workflow for anisotropic velocity analysis and NMO correction can proceed in three steps:
- Perform an isotropic velocity scan to estimate
and flatten seismic events.
- Perform a residual anisotropic moveout to account for
and
dependent terms simultaneously.
- Convert orthogonal parameters to more intuitive anisot-ropy parameters. For instance, the normal moveout velocity at azimuth
can be recovered by
 |
(4) |
In this procedure, the first two steps require a velocity-scan process. In fact, because
and
are symmetric in
, the single-parameter isotropic scan involved in the first step can be handled efficiently by a 2D butterfly algorithm, as discussed in our previous work (Hu et al., 2013,2012). Our goal in this paper is to speed up the more expensive, two-parameter velocity scan in the second step.
To be specific, what we need for residual moveout is to compute a semblance as follows (assuming that the
part has been moved out from the previous step):
 |
(5) |
where
is a 3D CMP dataset after isotropic moveout and
 |
(6) |
Subsections
 |
 |
 |
 | A fast algorithm for 3D azimuthally anisotropic velocity scan |  |
![[pdf]](icons/pdf.png) |
Next: Basic formulation
Up: Hu et al.: Fast
Previous: Introduction
2015-03-27