Apparent Anisotropy: Definition, NMO Velocity, and VTI Correction

Apparent anisotropy is a seismic velocity measurement that quantifies the ratio of the normal-moveout (NMO) stacking velocity to the true vertical interval velocity derived from well logs or a vertical seismic profile (VSP). When the two velocity measurements diverge, the subsurface appears anisotropic to the seismic acquisition system even if portions of the rock are individually isotropic. Geophysicists and petrophysicists rely on this parameter to diagnose whether velocity differences stem from genuine rock-fabric anisotropy, from fine-scale layering that the seismic wavelet cannot resolve, or from processing artifacts that have crept into the velocity field. Correctly identifying the source of apparent anisotropy is fundamental to accurate depth conversion, reliable pore-pressure prediction, and meaningful amplitude interpretation.

Key Takeaways

• Apparent anisotropy equals the ratio V_NMO / V_interval, where V_NMO is measured from surface seismic moveout and V_interval is measured vertically by a sonic log or zero-offset VSP.
• The primary physical cause in shale-rich basins is intrinsic VTI (vertical transverse isotropy): horizontal velocities exceed vertical velocities, so NMO sampling biases toward faster horizontal travel paths.
• Fine layering below seismic resolution produces an equivalent anisotropic response through Backus averaging, even when every individual layer is itself isotropic.
• The Alkhalifah-Tsvankin eta (η) parameter quantifies the degree of anelliptic moveout and is the standard correction factor applied during seismic processing to flatten gathers and improve depth conversion.
• Calibrating anisotropy corrections against check-shot or VSP interval velocities reduces depth-conversion errors in structurally complex areas by 5 to 15 percent in typical shale-carbonate sequences.

How Apparent Anisotropy Arises

When a compressional seismic wave travels from a surface source to a reflector and back, it samples a range of propagation angles. Near-offset traces sample nearly vertical paths; far-offset traces travel at angles that can exceed 40 degrees from vertical. Normal-moveout analysis fits a hyperbolic (or higher-order) curve to the arrival times across this offset range, producing a stacking velocity. In an isotropic medium, this stacking velocity equals the root-mean-square (RMS) velocity, which can be converted to interval velocity through the Dix equation. In an anisotropic medium, however, the horizontal velocity component that dominates the far-offset arrivals is faster than the vertical velocity sampled by a sonic log. The stacking velocity therefore overestimates the vertical velocity, and the ratio V_app / V_interval rises above unity. Values of this ratio typically range from 1.02 to 1.15 in organic-rich shales and laminated sequences, although extreme cases in the Devonian shales of the Appalachian Basin and the Vaca Muerta of Argentina have yielded ratios approaching 1.25.

Processing artifacts introduce a second class of apparent anisotropy that is unrelated to rock physics. Incorrect velocity picks on dipping reflectors cause NMO velocities to deviate from their true RMS values. Residual dip moveout, cycle-skipping during semblance analysis, and anisotropy mis-handling in pre-stack depth migration all contaminate the velocity field. Distinguishing artifact-driven apparent anisotropy from rock-physics-driven apparent anisotropy requires calibration wells with VSP surveys or at minimum sonic log ties combined with synthetic seismograms. In exploration settings without nearby well control, regional trends derived from analogue formations can be used as a guide, though uncertainty bounds must be widened accordingly.

A third mechanism involves structural complexity: horizontal or gently dipping layers produce standard hyperbolic moveout, but folded or faulted reflectors generate non-hyperbolic arrivals that a simple NMO velocity model misrepresents. This effect compounds genuine anisotropy in thrust-belt plays such as the Foothills of the Canadian Rockies, the Zagros Mountains of Iran and Iraq, and the Fold Belt of Papua New Guinea, where both rock-fabric anisotropy and structural dip must be addressed simultaneously before depth conversion can be trusted.

Intrinsic VTI Anisotropy in Shales and Laminated Sands

Most sedimentary basins contain thick shale intervals that display transverse isotropy about a vertical symmetry axis, commonly abbreviated VTI. In VTI media, elastic properties are identical in any horizontal direction but differ along the vertical axis. The horizontal P-wave velocity (V_P,h) exceeds the vertical P-wave velocity (V_P,v) because clay platelets align sub-horizontally during compaction, and the bond stiffness along the bedding plane is greater than stiffness across bedding. Sonic tools in boreholes measure vertical velocity; surface seismic gathers measure an offset-dependent mix dominated by sub-horizontal propagation at far offsets. The gap between the two measurements defines the apparent anisotropy attributable to intrinsic fabric.

Alkhalifah and Tsvankin (1995) introduced a convenient two-parameter VTI parameterization that separates the anelliptic character of the moveout from the NMO velocity itself. Their eta parameter is defined as:

η = (V_NMO² − V_interval²) / (2 × V_interval²)

A value of η = 0 implies no anellipticity (the medium is either isotropic or elliptically anisotropic). Values of η between 0.05 and 0.15 are common in shale-dominated basins worldwide. In the Montney Formation of northeastern British Columbia, η values averaging 0.08 to 0.12 have been documented from multi-offset VSPs, consistent with the strong horizontal layering fabric of that siltstone-shale alternation. In the Haynesville Shale of Louisiana and East Texas, η values in the 0.10 to 0.18 range reflect the highly organic and finely laminated character of the formation. Applying the η correction during pre-stack time migration flattens far-offset residuals on common-image gathers, improving stacking coherence and signal-to-noise ratio on final migrated volumes.

Backus Averaging and Sub-Resolution Fine Layering

Even when each individual bed in a thinly interbedded sequence is isotropic, the sequence as a whole can appear anisotropic to a long-wavelength seismic wave. This is the Backus averaging effect, formalized by George Backus in 1962. When layer thicknesses are small compared to the dominant seismic wavelength (typically less than one-tenth of the wavelength, or roughly 5 to 15 m for shallow targets and 15 to 40 m for deeper targets at 50 Hz dominant frequency), the seismic wave responds to an effective medium whose elastic constants are the thickness-weighted averages of the constituent layer stiffnesses. The effective medium is VTI even though no single layer has intrinsic anisotropy, and the horizontal velocity of the effective medium exceeds its vertical velocity whenever the sequence contains alternating fast and slow layers.

The Backus averaging equations for a two-component stack of layers with P-wave moduli M₁ and M₂ (where M = ρV_P²) and volume fractions f₁ and f₂ = 1 − f₁ yield a vertical effective modulus:

1 / C_33,eff = f₁ / M₁ + f₂ / M₂

and a horizontal effective modulus that is always greater than or equal to C_33,eff. The difference between horizontal and vertical effective moduli grows as the contrast between fast and slow layer velocities increases and as the layering becomes more regular. In cyclic turbidite sequences, regularly alternating sandstone and shale beds produce Backus anisotropy that rivals the intrinsic shale anisotropy in magnitude. This has direct consequences for the Deepwater Gulf of Mexico, where turbidite reservoirs in the Wilcox and Miocene sections display apparent anisotropy from both sources simultaneously, complicating both velocity model building and amplitude-variation-with-offset (AVO) analysis.

International Jurisdictions and Regional Context

Canada (Western Canada Sedimentary Basin): The Alberta Deep Basin and its extension into northeastern British Columbia host thick Cretaceous and Triassic shale-siltstone sequences including the Montney, Doig, Duvernay, and Nordegg formations. These units display apparent anisotropy values that consistently affect seismic depth conversion in the Foothills thrust belt. The Alberta Energy Regulator (AER) requires depth uncertainty documentation in well license applications for complex structural areas, and operators routinely employ anisotropic depth migration workflows to meet this requirement. The Canadian Society of Exploration Geophysicists (CSEG) has published several studies benchmarking anisotropy corrections for Horn River and Montney shale plays, with η values in the range 0.06 to 0.14 reported across the Dawson Creek and Fort St. John corridors.

United States (Appalachian Basin and Gulf of Mexico): The Marcellus and Utica shales of the northeastern US exhibit some of the strongest apparent anisotropy values documented in North America. VSP surveys in West Virginia and Pennsylvania have recorded η values above 0.20 in the most organic-rich Marcellus intervals. The U.S. Geological Survey (USGS) and industry consortium studies have shown that ignoring anisotropy in Marcellus depth conversion shifts predicted reservoir depths by up to 50 m (165 ft) at moderate offsets, a discrepancy that translates directly into wellbore positioning errors in pad-drilled horizontal wells. In the deepwater Gulf of Mexico, Backus anisotropy in Wilcox turbidites requires careful calibration before seismic reservoir characterization workflows can be trusted at 1:1 well-to-seismic ties.

Australia (Northwest Shelf and Cooper-Eromanga Basin): The Browse Basin and Carnarvon Basin on the Northwest Shelf contain thick Jurassic and Triassic shale sequences overlying major carbonate reservoirs. Operators developing LNG targets in the Ichthys, Browse, and Pluto fields have documented apparent anisotropy in the overburden shales that must be accounted for in the velocity model to correctly image the deeper carbonate reservoirs. The Cooper-Eromanga Basin in South Australia and Queensland features tightly interbedded Permian shales and coals above gas reservoirs, with Backus anisotropy in the coal measures producing η values of 0.04 to 0.10 that affect time-to-depth conversion of the Patchawarra and Tirrawarra sandstone targets.

Middle East (Arabian Platform): The giant carbonate fields of Saudi Arabia, Kuwait, UAE, and Iraq are covered by thick Aruma, Rus, and Damman shale-evaporite sequences. While carbonates themselves show minimal intrinsic anisotropy, the interbedded anhydrite and shale intervals in the overburden produce Backus anisotropy that biases depth conversion of the Arab Zone and Khuff carbonate reservoirs. Saudi Aramco and Abu Dhabi National Energy Company (TAQA) have published internal benchmarking studies noting depth conversion errors of 10 to 30 m in the Arab-D reservoir when overburden anisotropy is ignored, a significant concern given the lateral structural closure tolerances in these fields.

Norway and North Sea: The Norwegian Continental Shelf hosts thick Cretaceous chalk and shale overburdens above Jurassic Brent and Statfjord sandstone reservoirs. The Kimmeridge Clay Formation, a major source rock and seal, exhibits η values of 0.07 to 0.13 as documented in multi-offset VSP surveys on the Varg, Sleipner, and Oseberg fields. Equinor, TotalEnergies, and Shell have incorporated VTI anisotropy into their pre-stack depth migration workflows for all major Norwegian shelf projects since the mid-2000s. The Norwegian Petroleum Directorate (NPD, now renamed the Norwegian Offshore Directorate) includes velocity model documentation requirements in its exploration well reporting guidelines, which in practice necessitates explicit anisotropy parameterization for fields with significant shale overburden.