A convenient approximation of transverse
isotropy for higherorder moveout,
prestack time migration,
and depth calibration
William S. Harlan
August 1998
Processors use different RMS velocity models for
three steps of time imaging: NMO, DMO, and poststack time migration.
To perform prestack time migration in a single step, we must
use a single velocity model. A single step avoids an extra
stationaryphase approximation and should produce more accurate
results. Nevertheless, results are usually worse with a single
velocity model, unless different velocities are used for flat and dipping
reflections [8]. Those velocities which best
fit prestack normal moveouts over offset (flat reflections) do not best
focus the tails of diffractions over midpoint (dipping reflections).
Conventional processing hides this difference with inconsistent
velocity models for prestack moveout analysis and poststack migration.
Occasionally, processors want to use a higherorder
normalmoveout equation to flatten prestack reflections with long offsets.
Conventional moveout analysis does a good job of fitting the
difference in traveltime between near and far offsets,
but a higherorder moveout can better flatten any residual
bulge in the middle. Our parameterization of this normal moveout
should be consistent with the model of anisotropy
used in full prestack time imaging.
It is recognized that the kinematics of surface reflection seismic data
are insensitive to component of transverse isotropy that
is essential for an accurate conversion of time to depth.
We can isolate the parameters needed to fit surface
reflection times and allow depths to be calibrated independently.
A single convenient analytic approximation of transverse isotropy
will allow us to generalize prestack moveout, time imaging, and
depth conversion. No more parameters will be introduced than
necessary. Parameters are decoupled so that each can be
estimated in turn, if the additional degrees of freedom are
required to fit the data.
Much of this material originally appeared in an appendix to
Harlan [7].
Parameters for approximate transverse isotropy
Assume that anisotropic velocities have a vertical axis of symmetry,
like the transversely isotropic (TI) media described in
Thomsen [11]. Although that paper is
titled ``Weak elastic anisotropy,'' the same parameterizations
can be applied to very strong anisotropy [12].
Three of Thomsen's parameters, , , and ,
are defined by the elastic constants of a general TI medium.
These constants can be used to specify three different
effective velocities at a single point in the model.
is the velocity of a wave traveling vertically along
the axis of symmetry. The velocity in any horizontal direction
is defined by
A ``normal moveout velocity'' (NMO) velocity is defined by
If these TI properties represent the equivalent medium of
many isotropic layers [3,10],
then we can expect
[4].
Using Backus averaging, Phil Anno of Conoco has also shown that
we can expect and , if
the ratio and
have a positive correlation. These inequalities imply
that
. Well calibrations have shown
that shales can produce . Such shales
possess ``intrinsic'' anisotropy that cannot be
modeled as the macroscopic equivalent of isotropic layers.
Researchers at the Colorado School of Mines
[12,2] have also defined a constant
For an equivalent medium of isotropic layers .
Many combinations of three of these parameters
can be used to describe
a TI medium for compressional P waves with a known axis of symmetry.
Such an approximation has already dropped a fourth constant (shear
wave velocity) to which compressional waves are insensitive.
The parameter is an excellent parameter choice for maximum
sensitivity to the kinematics of surface measurements only.
If TI velocities are parameterized by , , and ,
we find that surface measurements are very insensitive to .
We could not make the same claim if were used instead of
.
I prefer the three velocities and because
they share the same units. Surface measurements are very insensitive
to , given values for and .
This choice is not the most convenient for processing, however.
The exact equations for TI phase velocity as a function of
angle are rather clumsy, and no explicit form is available
for group velocity. Explicit approximate equations can
fit the same family of curves almost as well
as the original correct equations [9].
I use an approximate equation for group velocity
which appears to emulate closely the exact curves for
large ranges of positive and negative .
Estimated curves usually have larger statistical errors from noisy data
than introduced by these approximations.
Kirchhoff migrations can calculate traveltimes by integrating
the group velocity along a stationary Fermat raypath.
Fourierdomain implementations can use only
the equation for phase velocities. Some raytracing methods
use both, because phases must be matched across
discontinuous boundaries.
I choose approximate curves with the three velocities and .
Let be the group angle of a raypath from the vertical.
Then the group velocity can be expressed as
Compare Byun et al [5], who use
the same approximation with different parameters.
Greg Lazear of Conoco found that a symmetric equation
approximates the phase velocity
as a function of the phase angle ,
but with reciprocals of the same velocity parameters:



(15) 
Compare this phase equation closely to the group equation (11).
I know of no other approximation that allows such symmetry.
The normalmoveout (NMO) velocity has a physical
interpretation to justify its name.
Imagine an experiment on a homogeneous
and anisotropic medium, or imagine a small scale experiment
on a smooth model. Measure the traveltime between two
points placed on a vertical line, separated by a vertical
distance . Now displace the
upper point a distance along a horizontal line (a normalmoveout)
and measure the new traveltime .
Then according to equation (11) the traveltime
as a function of offset is exactly



(16) 
For small offsets , the value of in this
``moveout equation'' is
dominated by the NMO velocity rather than . For large
offsets , the raypath is almost horizontal
and dominates.
I find it convenient to define a stacking velocity
as a function of the offset for a fixed vertical distance
:
I use the term stacking velocity because I want to suggest
the bestfitting curve over a finite range of offsets, as
you would prefer for a stacking or semblance analysis.
Simplify the moveout equation (16) to fit a pseudohyperbola:



(20) 
The stacking velocity covers the range
for a Backus equivalent medium with
negative , increasing in value as increases.
When , the curve is exactly hyperbolic, and .
Notice that this stacking velocity can measure a local property
as well as an average to the surface.
To use twoway reflection times in (20)
we need only replace the half offset
by the full offset.
Three measurements of traveltimes at three
different offsets should uniquely determine the three
velocity constants , , . The
traveltimes are much more sensitive to , which determines
moveouts at small offsets, and to , which determines
moveout at larger offsets. The vertical velocity affects
only the rate at which the stacking velocity (18)
changes from one limit to the other. As long as has roughly the
correct magnitude, then we can fit all measured traveltimes
very well.
For imaging surface data in time, we acknowledge our insensitivity to
and can approximate it with another value.
We can approximate
and simplify
stacking velocity (19) even further, as in



(21) 
This new equation depends only on two parameters, and .
A better approximation might be
, or
equivalently
, which is probably
closer to commonly observed values.
NMO equation (16) now is equivalent to



(22) 
which is equivalent to equation (5) in
Alkhalifah [1] and
equation (7) in Grechka and Tsvankin
[6]. Both these publications derive from
Tsvankin and Thomsen [12], which
uses an asymptotic correction of a Taylor expansion to arrive
at this approximation.
The equivalent stacking velocity is then



(23) 
The differences between these two approximate
stacking velocities (23)
and (21) are negligible for numerical work.
I will use the simpler version (21).
Moveout analyses determine the stacking velocity for a specific aperture of
offsets. Define our best isotropic approximation of the velocity
to be the stacking velocity at the maximum offset
of the aperture:
Or we can assume that we know the aperture angle from the vertical,
so that
The best isotropic velocity is a simple function of
infinitesimaloffset NMO velocity and the anisotropic parameter
.
Similarly we can use the definition of in (8)
to rewrite the above equation (26) as



(27) 
This form will prove to be very useful when rewriting
our group velocity equations
(13) and (14).
A conventional velocity analysis produces
densely picked values for .
The anisotropic parameter adjusts the
moveout between near and far offsets.
For anisotropic moveout analysis, we could use the
following offsetdependent stacking velocity to scan for ,
holding constant:



(28) 
Conventional moveout analysis is not very sensitive
to the anisotropic parameter except
for unusually wideaperture data, with offsets greater than depth.
Much more anisotropic information is available by performing a full prestack
time migration.
For prestack time migration, we can expect that conventional
analysis for will best describe the moveouts of flat
reflections. Dipping reflections are difficult to pick in
prestack semblance analysis because they are sparser and
move across midpoints, from gather to gather, as migration
velocity changes.
Holding constant, we can perturb at a low
spatial resolution until the imaging of steep reflections improves.
As changes, there will be a small adjustment of the
bulge in flat reflections over offset, with little effect
on a fully timemigrated stack. The moveouts of dipping reflections,
on the other hand, will change drastically as changes, with swings
from positive to negative residual moveouts, and with lateral
movement over midpoint. As improves, you should see faultplane
reflections sharpen and focus in targeted prestack timemigrated images
Flat reflections should change little, and should require little
revision after updating . By contrast, an optimization of
and requires both to be adjusted simultaneously,
with equal resolution.
The group velocity equation (13) is useful
for Kirchhoff depth imaging.
We can focus images very well with good values for and ,
(or and ),
then adjust imaged depths to tie wells with
(or , or ) while
holding the other two parameters constant.
If we have used tomographic methods to estimate
isotropic interval velocities, then
we should attempt also to estimate the effective aperture angle
at each depth in the model.
To recapituate our algebra, we combine the definitions of
(9) and of stacking velocity
(19) to solve for the horizontal velocity:
At maximum offset, we recover our previous
relationship (27) between and .
Substitute this horizontal velocity into the group velocity
(13).
We can then adjust the isotropic velocities with and
according to



(31) 
The isotropic velocity best explains the moveouts and traveltimes
of relatively flat reflections over the finite aperture.
The parameter modifies these velocities at high dips, to image steep
reflections better without degrading the imaging of lowdip reflections.
The third parameter has little effect on measured
surface traveltimes at any dip (holding and constant),
but can be adjusted as necessary to tie wells.
We can also use layered medium theory to
predict this from estimated
and . Or if shale dominates, with strong
intrinsic anisotropy and , then correlations
can be calibrated for a given area.
At worst, we know , so we can assume a default
value of
, as appropriate
for a given area. Such a default value is
still better than a default value of .
These three parameters , , and
do not completely decouple the steps of anisotropic
velocity analysis, but they should minimize the number
of iterations necessary for revisions.
If you prefer to use instead of
as the third degree of freedom, then simply use
the definition of in (7) for



(32) 
Although I greatly prefer the approach in the preceding
section, many prefer to treat their estimated isotropic
velocities as equivalent to NMO velocities .
Such an assumption is not a bad one if angles are limited
during interval velocity estimation. Dix inversion
of stacking velocities may be closer to if stacking
velocities were consciously optimized for inner offsets.
The Common Reflecting Surface tomography of Karlsruhe University
inverts only the curvature of reflection traveltimes
around zerooffsets.
For such approaches, one might prefer a different triplet
of velocity parameters: , , and .
Group velocity can be described by replacing in
the approximation (14) with
and as in equation (10).
Additionally we can replace with ,
using the definition of in (7):
Again and should be sufficient to model
all surface reflection traveltimes, for all dips.
Holding these two constant, we can adjust either
or to tie known well depths.
Many thanks to Greg Lazear and Phil Anno, who spent
months with me at Conoco trying to reconcile and simplify
different parameterizations of transverse isotropy.
Thanks to Andreas Rüger for a very careful review of this paper.
 1

Tariq Alkhalifah.
Velocity analysis using nonhyperbolic moveout in transversely
isotropic media.
Geophysics, 62(6):18391854, 1997.
 2

Tariq Alkhalifah and Ilya Tsvankin.
Velocity analysis for transversely isotropic media.
Geophysics, 60(5):15501566, 1995.
 3

George E. Backus.
Longwave elastic anisotropy produced by horizontal layering.
Journal of Geophysical Research, 67(11):44274440, 1962.
 4

J.G. Berryman.
Longwave elastic anisotropy in transversely isotropic media.
Geophysics, 44:896917, 1979.
 5

B. S. Byun, D. Corrigan, and J. E. Gaiser.
Anisotropic velocity analysis for lithology discrimination.
Geophysics, 54(12):15641574, 1989.
 6

V. Grechka and I. Tsvankin.
Feasibility of nonhyperbolic moveout inversion in transversely
isotropic media.
Geophysics, 63(3):957969, 1998.
 7

William S. Harlan.
Flexible seismic traveltime tomography applied to diving waves.
Stanford Exploration Project Report,
http://sepwww.stanford.edu/research/reports/,
89, 1995.
 8

Walt Lynn, A. González, and S. Mackay.
Where are the faultplane reflections?
In 61st Annual International Meeting, SEG, Expanded Abstracts,
volume 91, pages 11511154. Soc. Expl. Geophys., 1991.
 9

R. J. Michelena, F. Muir, and J. M. Harris.
Anisotropic traveltime tomography.
Geophys. Prosp., 41(4):381412, 1993.
 10

M. Schoenberg and F. Muir.
A calculus for finely layered anisotropic media.
Geophysics, 54(5):581589, 1989.
Discussion and reply by authors in GEO564572576.
 11

L. Thomsen.
Weak elastic anisotropy.
Geophysics, 51(10):19541966, 1986.
Discussion in GEO530405580560 with reply by author.
 12

I. Tsvankin and L. Thomsen.
Nonhyperbolic reflection moveout in anisotropic media.
Geophysics, 59(8):12901304, 1994.
harlan
20180808