EDA: Stress-Aligned Microcracks, Shear-Wave Splitting, and Azimuthal Anisotropy

EDA, extensive-dilatancy anisotropy, is a form of azimuthal seismic anisotropy that arises when a rock volume is pervaded by stress-aligned, fluid-filled microcracks and fractures that are not horizontal, so that the velocity of a seismic wave depends on the compass direction in which it travels. The acronym is the working shorthand geophysicists use for the rock-physics model in which a population of near-vertical, parallel cracks is held open by the in-situ stress field, with their planes oriented perpendicular to the minimum horizontal stress and parallel to the maximum horizontal stress. Waves polarized or propagating parallel to those crack faces travel faster than waves crossing them, because crossing a fluid-filled crack is mechanically softer than travelling along intact rock between cracks. The diagnostic signature of EDA is shear-wave splitting, also called birefringence: a single shear wave entering the cracked medium divides into a fast component polarized parallel to the dominant crack and fracture strike and a slow component polarized perpendicular to it, and the two arrive separated by a measurable delay. That delay time is a direct measure of the intensity and thickness of the cracked interval, while the polarization azimuth of the fast shear wave reveals the orientation of the open fractures and, by extension, the direction of maximum horizontal stress. This makes EDA enormously useful in the Western Canadian Sedimentary Basin, where open natural fractures and the present-day stress azimuth govern how a Montney or Duvernay horizontal well will be landed and how its multi-stage hydraulic fracture will propagate. Multicomponent and wide-azimuth seismic surveys are designed specifically to capture azimuthal velocity and amplitude variations diagnostic of EDA, feeding amplitude-versus-offset-and-azimuth, or AVAZ, analysis that maps fracture density and orientation across a lease before a single well is spudded. The concept was introduced by Stuart Crampin in the late 1970s and refined through the 1980s as a unifying explanation for shear-wave splitting observed almost universally in crustal rocks below a few hundred metres, and it later connected to the anisotropic poroelasticity framework describing how cracks reorganize as pore pressure changes. For the working interpreter, EDA is the bridge between a seismic observable, split shear waves, and a reservoir property that matters for drilling and completion: the orientation, intensity, and fluid state of the fracture fabric. Recognizing and quantifying it turns azimuthal seismic data into a predictive tool for natural fracture sweet spots and for aligning well azimuth and frac design with the stress field.

Key Takeaways

  • Stress-Aligned Crack Model: EDA describes anisotropy produced by near-vertical, fluid-filled microcracks held open parallel to the maximum horizontal stress and perpendicular to the minimum stress. Because the cracks are not horizontal, velocity varies with azimuth: travel parallel to crack strike is fast, travel perpendicular is slow, which is the defining contrast that distinguishes azimuthal anisotropy from simple flat-layer velocity layering.
  • Shear-Wave Splitting Signature: The diagnostic of EDA is birefringence, where one shear wave splits into a fast component polarized along fracture strike and a slow component polarized across it. The two waves arrive separated by a delay, and detecting that split in multicomponent WCSB data is the primary field evidence that a stress-aligned crack fabric is present in the interval.
  • Delay Time Maps Intensity: The time separation between fast and slow shear waves scales with both the crack density and the thickness of the anisotropic interval. A larger delay over a Montney or Duvernay section implies a more intensely or more thickly fractured rock volume, giving interpreters a quantitative proxy for fracture density rather than a yes-or-no fracture flag.
  • Fast Polarization Gives Stress Azimuth: The polarization direction of the fast shear wave aligns with the open fracture strike and therefore the maximum horizontal stress direction. In the WCSB that azimuth, often northeast to southwest in much of Alberta, dictates the optimal direction to drill a horizontal well perpendicular to it so that transverse hydraulic fractures open against the minimum stress.
  • Feeds AVAZ Fracture Mapping: Wide-azimuth and multicomponent seismic surveys exploit EDA through amplitude-versus-offset-and-azimuth analysis, converting azimuthal velocity and amplitude variation into maps of fracture density and orientation. These maps let operators high-grade natural-fracture sweet spots and tune completion design across a lease before committing capital to a horizontal program.

Why The Cracks Must Be Non-Horizontal

The azimuthal character of EDA depends entirely on the cracks departing from horizontal. A medium of perfectly flat, parallel layers or horizontal cracks produces vertical transverse isotropy, where velocity varies with the angle from vertical but is the same in every compass direction, so a vertically travelling wave sees no azimuthal change. EDA instead invokes near-vertical, stress-aligned cracks, which break that rotational symmetry: a wave now cares not just about its inclination but about its azimuth relative to the crack strike. This is why a vertically propagating shear wave, which would be blind to flat layering, splits under EDA and carries fracture-orientation information, making it the multicomponent geophysicist's most direct probe of the in-situ stress and fracture system.

From Crustal Observation To Reservoir Tool

EDA began as an explanation for a puzzle in earthquake seismology, the near-universal observation that shear waves travelling through crustal rock below a few hundred metres arrive split. Crampin's model attributed this to a pervasive fabric of stress-aligned fluid-filled microcracks, and proposed that temporal changes in splitting might even precede earthquakes as cracks reorganize under changing stress. The petroleum industry adopted the same physics for a commercial purpose: if split shear waves reveal crack orientation and density, then azimuthal seismic can map natural fractures in tight reservoirs. That translation turned a crustal-stress concept into a routine WCSB exploration input for the Montney, Duvernay, and fractured carbonate plays.

Fast Facts

One of the more striking claims attached to EDA is its near-universality: Crampin argued that essentially all sedimentary, igneous, and metamorphic rocks below roughly 500 to 1,000 m show some degree of stress-aligned shear-wave splitting, implying that the upper crust is everywhere a weakly anisotropic, crack-pervaded medium rather than an isotropic solid. The same framework led to the bold proposal that monitoring changes in shear-wave delay could serve as an earthquake stress gauge, a controversial idea that nonetheless pushed multicomponent seismology forward and seeded the fracture-detection methods now used commercially.

EDA is the abbreviated name for extensive dilatancy anisotropy, and the two entries describe the same crack-controlled phenomenon from complementary angles. It is detected through shear-wave splitting of S-wave energy recorded on multicomponent surveys, and it is the azimuthal cousin of the broader idea of seismic anisotropy, where wave speed depends on direction. The fracture and stress information EDA yields directly informs hydraulic hydraulic fracturing design and horizontal well azimuth in WCSB tight reservoirs.

Real-World WCSB Scenario: Azimuthal Fracture Mapping On A Duvernay Block

Ahead of a Duvernay development near Willesden Green, Alberta, an operator commissions a wide-azimuth 3D reprocessing and AVAZ study at roughly CAD 250,000 to extract EDA-driven fracture attributes across the 90 square kilometre block. The fast shear-wave polarization azimuth maps a consistent northeast-southwest maximum horizontal stress, and the fast-slow delay time highlights a corridor of elevated fracture intensity along a subtle structural flexure. That corridor becomes the high-grade target for the first pad.

The operator orients the horizontal laterals northwest-southeast, perpendicular to the mapped stress azimuth, so that transverse hydraulic fractures open cleanly against the minimum stress. The CAD 11 million two-well pad outperforms an offset pad drilled before the EDA work by a meaningful margin on first-year liquids, validating azimuthal anisotropy as a pre-drill predictive tool.