Extensive Dilatancy Anisotropy: Crampin's Crack Model, Fast Shear Polarization, and Fracture Characterization
Extensive dilatancy anisotropy, abbreviated EDA, is a rock-physics model and a form of azimuthal seismic anisotropy in which a rock mass is pervaded by stress-aligned, fluid-filled microcracks and fractures oriented other than horizontally, causing seismic velocity to vary with the direction of propagation through the rock. The model holds that throughout much of the upper crust the in-situ stress field opens and aligns a population of thin, near-vertical, parallel cracks whose planes strike parallel to the maximum horizontal compressive stress and lie perpendicular to the minimum horizontal stress. A seismic wave travelling parallel to those aligned cracks moves faster than a wave travelling perpendicular to them, because crossing a fluid-filled crack is mechanically more compliant than moving along the intact rock framework between cracks. The phrase itself encodes the physics: dilatancy is the stress-driven opening of cracks that slightly increases rock volume, extensive signals that the cracks are not isolated flaws but a pervasive distributed fabric, and anisotropy names the resulting directional dependence of wave speed. Stuart Crampin introduced the concept in 1978 and developed it through the 1980s to explain why crustal shear waves observed in refraction and earthquake studies almost everywhere arrive split into two polarizations, a behaviour difficult to account for with isotropic or simple layered models. The hallmark observable is shear-wave splitting: a shear wave entering an EDA medium separates into a fast wave polarized parallel to the dominant crack and fracture strike and a slow wave polarized perpendicular to it, the delay between them growing with the density and thickness of the cracked zone. Because the fast polarization azimuth tracks the open-crack strike and therefore the maximum horizontal stress, and because the delay time measures crack intensity, extensive dilatancy anisotropy converts a seismic measurement into two reservoir-critical quantities, fracture orientation and fracture density. In the Western Canadian Sedimentary Basin this is precisely the information needed to plan horizontal wells and multi-stage hydraulic fracture treatments in the Montney, the Duvernay, and naturally fractured Devonian carbonates such as the Slave Point and Nisku. The model also connects to anisotropic poroelasticity, which describes how the crack population reorganizes geometry as pore pressure rises and falls, opening or closing fluid-filled microcracks and changing the splitting over time. Crampin further proposed that monitoring such temporal changes in shear-wave splitting might act as a stress gauge with earthquake-precursor potential, an ambitious and debated claim that nevertheless drove the development of multicomponent seismology. For petroleum geophysicists the practical legacy is concrete: extensive dilatancy anisotropy is the theoretical foundation under azimuthal fracture detection, the reason wide-azimuth and multicomponent surveys can map stress and natural fractures before drilling.
Key Takeaways
- Pervasive Stress-Aligned Cracks: Extensive dilatancy anisotropy attributes directional velocity to a distributed fabric of near-vertical, fluid-filled microcracks opened and aligned by the in-situ stress, striking parallel to maximum horizontal stress. Extensive marks the fabric as pervasive rather than localized, distinguishing the model from a single fracture set and explaining anisotropy seen across whole crustal volumes, not just discrete fault zones.
- Crampin's 1978 Origin: Stuart Crampin introduced EDA in 1978 to explain the near-universal splitting of crustal shear waves that isotropic models could not reproduce. The model unified scattered observations of birefringent shear arrivals under a single mechanism, fluid-filled stress-aligned cracks, and became the conceptual backbone of multicomponent anisotropy studies in both earthquake and exploration seismology.
- Fast And Slow Shear Polarizations: A shear wave in an EDA medium splits into a fast wave polarized along crack strike and a slow wave polarized across it. The fast polarization azimuth reveals open-fracture orientation and maximum horizontal stress direction, while the time delay between the two scales with crack density and the thickness of the anisotropic interval, giving both orientation and intensity.
- Anisotropic Poroelasticity Link: EDA connects to the anisotropic poroelasticity model in which microcracks reorganize as pore pressure changes, opening and closing fluid-filled cracks and altering the splitting over time. This time-dependence is why Crampin argued shear-wave splitting could track evolving stress, and why repeat surveys can in principle monitor pressure and fracture changes during WCSB production.
- Foundation Of Azimuthal Fracture Detection: The practical payoff is that wide-azimuth and multicomponent seismic surveys map natural fractures and stress by measuring the velocity and amplitude variations EDA predicts. Amplitude-versus-offset-and-azimuth workflows built on this physics let operators target fractured sweet spots and align horizontal wells and fracture stages with the stress field before drilling.
Dilatancy, Cracks, And The Meaning Of The Name
Each word in extensive dilatancy anisotropy carries physical weight. Dilatancy is the well-known tendency of stressed rock to increase in volume as microcracks open, a behaviour measured in rock-mechanics laboratories long before it was invoked seismically. Extensive signals that, in Crampin's model, these opened cracks are not rare isolated flaws but a distributed, pervasive population filling most of the rock below a few hundred metres. Anisotropy is the consequence: with cracks preferentially aligned by stress, the rock transmits seismic energy faster along the crack planes than across them. Put together, the name compactly states the hypothesis that pervasive stress-opened cracks make ordinary crustal rock directionally variable to seismic waves, which is exactly what shear-wave splitting reveals.
Measuring Fracture Density And Orientation For Completions
Turning EDA into a completion input means quantifying two attributes from azimuthal data. Fracture orientation comes from the fast shear-wave polarization azimuth or, in P-wave AVAZ, from the azimuth of maximum velocity and amplitude, both tracking the open-crack strike and the maximum horizontal stress. Fracture intensity comes from the magnitude of the azimuthal variation or the fast-slow delay time, a proxy for crack density times interval thickness. In a Montney or Duvernay program these attributes guide two decisions: drilling the lateral perpendicular to maximum stress so transverse fractures open efficiently, and concentrating stages where natural fracture intensity is highest, improving stimulated reservoir volume per dollar of completion spend.
Fast Facts
Crampin's EDA carried a provocative corollary that still divides geophysicists: he argued the upper crust is so pervasively crack-pervaded that shear-wave splitting is essentially ubiquitous below 500 to 1,000 m, and that measured changes in splitting delay could serve as a stress monitor capable of forecasting earthquakes. The earthquake-precursor claim was never widely accepted, yet the underlying physics proved commercially durable: the same stress-aligned cracks that may or may not warn of seismicity reliably reveal fracture orientation in tight gas and shale reservoirs, giving an esoteric crustal model a profitable second life in WCSB exploration.
Related Terms
Extensive dilatancy anisotropy is the full name behind the abbreviation EDA, and the two glossary entries treat the same crack-controlled phenomenon. It is observed through the splitting of S-wave energy on multicomponent data and is a specific cause of azimuthal seismic anisotropy, the broader directional dependence of wave velocity. The fracture orientation and stress azimuth it yields directly shape hydraulic fracturing design and horizontal well planning in fractured WCSB reservoirs.
Real-World WCSB Scenario: Pre-Drill Fracture Characterization In The Montney
A Montney operator near Gold Creek, Alberta, integrates an EDA-based azimuthal anisotropy study into pre-drill planning for a six-well pad, spending about CAD 280,000 on wide-azimuth reprocessing and AVAZ inversion. The fast shear and P-wave azimuthal attributes confirm a maximum horizontal stress trending roughly northeast-southwest and reveal that fracture intensity is markedly higher in the eastern half of the block, where a deeper structural fabric appears to enhance the stress-aligned crack density predicted by the model.
Acting on the maps, the operator shifts three of the six laterals eastward and orients all of them northwest-southeast, perpendicular to the interpreted stress. The completion design loads more stages into the high-intensity corridor, and the CAD 33 million pad delivers stronger early liquids than the type curve, a result the asset team credits in part to letting extensive dilatancy anisotropy guide both well placement and fracture staging.