Shear Strain: Engineering and Tensor Definitions, Mohr Circle Analysis, and Reservoir Geomechanics in the Montney and Duvernay
Shear strain is the dimensionless measure of distortional deformation in a continuum body, defined as the tangent of the change in angle between two material lines that were initially perpendicular to each other; for small deformations relevant to almost all subsurface engineering applications, the tangent approximation reduces to the angle change itself measured in radians, and the symbol used is typically gamma for engineering shear strain or epsilon-xy for the tensorial shear strain, where engineering shear strain equals twice the tensor shear strain (gamma_xy equals 2 times epsilon_xy). In contrast to normal strain, which describes the fractional change in length along a single axis, shear strain captures the angular distortion experienced when two adjacent particles slide past each other while the material itself remains intact, producing the rhomboidal deformation of an initially rectangular element. The state of strain at any point in a rock body is fully described by the strain tensor with six independent components (three normal strains and three shear strains in three-dimensional space), and through eigenvalue decomposition the tensor yields three principal strains aligned with three mutually orthogonal principal strain axes along which shear strain is zero and only normal strain operates. Shear strain matters profoundly in oil and gas geomechanics because rock failure in compression, tension, and shear is controlled by deviatoric stress and strain rather than mean stress alone, and most reservoir, caprock, and wellbore failures from wellbore stability shear breakouts in the Duvernay Shale, to caprock fracture during SAGD steam injection in the McMurray Formation, to slip on bedding planes during Montney hydraulic fracturing are governed by shear-strain accumulation that drives the rock onto its Mohr-Coulomb or other failure envelope. WCSB operators use shear-strain analysis in 4D finite-element coupled reservoir-geomechanics simulators such as VISAGE, ABAQUS, and CMG GEM coupled with FLAC3D to predict caprock integrity for in-situ bitumen recovery under AER Directive 086, to design hydraulic fracture spacing in the Montney and Duvernay so that adjacent stages do not generate destructive interference, and to forecast compaction-induced subsidence and casing deformation in the Liard Basin gas pools and Lloydminster heavy oil thermal floods. Engineering shear strain is also the deformation quantity measured directly by downhole fiber-optic distributed acoustic and strain sensing (DAS and DSS) during fracture diagnostics, where strain rates on the order of microstrain per second (10 to the negative 6 per second) reveal hydraulic fracture geometry, propagation rate, and stress shadow interference with offset wells in real time.
Key Takeaways
- Engineering Versus Tensor Definition: Engineering shear strain (gamma) equals the tangent of the angle change between two initially perpendicular material lines, while tensor shear strain (epsilon_xy) equals half that quantity. The factor-of-two distinction is critical when applying Hooke's law, Mohr circle constructions, and finite-element formulations: shear stress equals shear modulus times engineering shear strain (tau equals G times gamma), but the strain tensor used in matrix operations uses the half-magnitude tensor form to maintain symmetry.
- Mohr Circle Representation: On a Mohr circle plot of normal stress versus shear stress, the diameter equals the principal stress difference (sigma_1 minus sigma_3), and the maximum shear stress equals half this difference and occurs on planes oriented at 45 degrees to the principal stress axes. The corresponding maximum shear strain in linear elasticity equals (sigma_1 minus sigma_3) divided by twice the shear modulus G, providing the working calculation for shear failure prediction in WCSB shales with G values of 8 to 20 GPa.
- Shear Modulus and Rock Stiffness: Shear strain relates to shear stress through the shear modulus G, with typical WCSB reservoir values ranging from 6 to 12 GPa for Montney siltstone, 8 to 14 GPa for Duvernay shale, 0.5 to 2.5 GPa for McMurray bitumen-saturated sand, and 15 to 25 GPa for Leduc reef dolomite. G correlates with Young's modulus E through Poisson's ratio nu as G equals E divided by 2 times (1 plus nu), with nu commonly 0.20 to 0.25 in tight shales and 0.30 to 0.40 in unconsolidated bitumen sands.
- Failure and Plasticity: Shear strain accumulation drives rock onto the Mohr-Coulomb failure envelope, defined by cohesion c and friction angle phi. Once failure initiates, plastic shear strain accumulates without further stress increase, producing the strain-softening behaviour seen in Montney and Duvernay cores at differential stresses above 35 to 55 MPa (5,000 to 8,000 psi). Localized shear bands form at angles of 45 degrees minus half the friction angle to the maximum principal stress direction, typically 25 to 35 degrees in WCSB reservoir rocks.
- Fiber Optic Strain Measurement: Distributed strain sensing along fiber-optic cables cemented behind casing measures axial strain at micrometre resolution along the entire wellbore length, capable of detecting fracture hits, stress shadows, and bedding slip during stimulation. WCSB operators including ARC Resources in the Montney and Chevron Canada in the Duvernay routinely deploy DAS and DSS to optimize stage spacing and identify shear strain anomalies indicating poor cluster efficiency.
Wellbore Stability and Shear Failure
Around a vertical or deviated wellbore, the in-situ stress field is locally redistributed, producing compressive tangential stress concentration at the borehole wall in the direction of the minimum horizontal stress and tensile stress in the direction of the maximum horizontal stress. When the shear strain at the wall exceeds the rock's shear capacity, breakouts form along the minimum stress direction. WCSB operators in the deep Montney run mud densities of 1,650 to 1,850 kg/m3 (13.8 to 15.4 ppg) to control shear strain at the wall, with mud weight design tied to elastic-plastic geomechanics models calibrated against image log breakouts and dipole sonic anisotropy measurements.
Hydraulic Fracture Geometry and Stress Shadow
Hydraulic fracturing produces pressurized planar tensile openings that perturb the local stress field of surrounding rock, raising the minimum horizontal stress by 0.5 to 3.5 MPa (75 to 510 psi) within one to two fracture half-lengths and creating shear strain on bedding planes and natural fractures up to 50 metres (165 ft) away. WCSB stage spacing in the Montney and Duvernay has tightened from 30 to 40 m (98 to 130 ft) in 2015 to 8 to 15 m (26 to 49 ft) by 2025 as operators learned to manage stress shadow interference through pump rate, proppant loading, and cluster efficiency monitoring with fiber optics.
Fast Facts
The 2014 Crooked Lake earthquake swarm near Fox Creek, Alberta, generated felt seismic events up to magnitude 4.4 on the Richter scale, linked by the AER and the Geological Survey of Canada to shear strain release on basement faults reactivated by Duvernay hydraulic fracturing in the Kaybob area. The events drove the AER to enact Subsurface Order 2 in 2015 requiring traffic-light response protocols, immediate shutdown for events at magnitude 4.0 or greater, and modified pumping schedules within 5 km (3.1 miles) of mapped basement structures, a regulatory innovation that has been adopted internationally as a model for induced-seismicity management.
Related Terms
Shear strain is intertwined with several rock mechanics concepts. Shear stress is the conjugate quantity to shear strain, linked through the shear modulus in elastic theory. Young's modulus describes uniaxial stiffness and connects to shear modulus via Poisson's ratio. Mohr's circle provides the graphical framework for transforming stress and strain between coordinate systems and identifying maximum shear conditions. Poisson's ratio describes the ratio of transverse to axial strain and controls how shear modulus scales with Young's modulus in linear elastic materials.
WCSB Application: Montney Stage Spacing Optimization
An operator in the British Columbia Montney completing a 3,200 m (10,500 ft) lateral with 60 stages of 80,000 lb (36,300 kg) proppant per stage observed early production decline 18 percent below type curve, triggering a fiber-optic DAS and DSS evaluation on the next pad at a cost of CAD 950,000 per well. The fiber sensing revealed that stages 12 through 24 generated shear strain perturbations exceeding 200 microstrain at offset wells located 150 m (490 ft) away, indicating that propped fractures were intersecting the heel of the offset lateral and short-circuiting drainage. Stage spacing was reduced from 50 m (164 ft) to 25 m (82 ft) and pump rate dropped from 16 m3/min (100 bbl/min) to 11 m3/min (69 bbl/min) on subsequent pads.
The redesigned completion increased EUR per well from 6.8 Bcf (193 e6m3) to 9.4 Bcf (266 e6m3), recovered the CAD 950,000 fiber cost within four months, and was applied across the operator's 240-well development plan for an aggregate present-value benefit of approximately CAD 380 million.