Elastic (Rock Mechanics)

Key Takeaways

• Seismic waves propagate as elastic deformations and seismic velocities directly measure elastic moduli — seismic compressional (P) waves and shear (S) waves are stress waves that travel through rock by alternately compressing and extending (P-waves) or shearing (S-waves) the rock in a fully reversible, elastic manner; the velocities of these waves are directly related to the rock's elastic moduli and bulk density by: Vp = sqrt((K + 4G/3)/ρ) and Vs = sqrt(G/ρ); measuring Vp and Vs from sonic logs or seismic surveys and combining them with density allows calculation of K, G, E, and ν for the formation, converting the seismic measurement into the mechanical properties needed for geomechanical analysis; the relationship between seismic velocities and mechanical properties is one of the most powerful connections in applied geophysics, allowing wellbore-scale and seismic-scale measurements to be integrated in a common mechanical earth model that describes formation properties throughout the field volume rather than only at well locations.
• Static and dynamic elastic moduli differ systematically and the difference must be accounted for in geomechanical models — dynamic elastic moduli are measured at seismic or sonic frequencies (tens to thousands of Hertz) at very small strain amplitudes (microstrains), while static elastic moduli are measured in the laboratory on core samples deformed at low rates to large strains (hundreds of microstrains); dynamic moduli are consistently higher than static moduli for the same rock because: at high frequencies, the rock matrix and pore fluid are undrained (pore fluid does not have time to redistribute, making the rock appear stiffer than its drained state); at small strains, the rock's microcracks are closed by the confining stress and do not contribute to deformation (making the rock appear stiffer than it would under larger strains that open microcracks); converting dynamic moduli from logs to static moduli for use in wellbore stability and fracture design calculations requires empirically derived correlation equations specific to each formation type, typically calibrated using matched laboratory static and dynamic measurements on core samples from the formation of interest.
• Pore pressure and effective stress govern the elastic response of reservoir rocks during production — rock elastic moduli depend on the effective stress (the difference between total confining stress and pore pressure), with stiffness increasing as effective stress increases; as reservoir pressure depletes during production, effective stress increases on the reservoir rock framework, causing elastic compaction (volumetric reduction) of the reservoir that reduces pore space and drives some additional hydrocarbon production (compaction drive) — a significant production mechanism in chalk reservoirs (North Sea Ekofisk field, where chalk compaction has driven both production and seafloor subsidence of several meters) and unconsolidated sandstone reservoirs; in overpressured formations, the pore pressure is close to the total stress, effective stress is low, and the rock is relatively soft and deforms easily under applied loads; drilling into overpressured formations requires careful mud weight management because the low effective stress means the rock's elastic compressive strength threshold is lower than in normally pressured formations at the same depth.
• Fracture mechanics in hydraulic fracturing assumes elastic deformation of the host rock to predict fracture geometry — hydraulic fracture models (PKN, KGD, Radnitz-Sneddon, and pseudo-3D models) treat the reservoir rock as a linear elastic medium and use the elastic moduli (particularly Young's modulus and Poisson's ratio) to relate the net pressure in the fracture (injection pressure minus closure stress) to the fracture width and the stress intensity at the fracture tip that drives propagation; higher Young's modulus (stiffer rock) requires higher net pressure to achieve the same fracture width, which is relevant for design of fracturing treatments in stiff carbonates versus more compliant shales; Poisson's ratio controls the relationship between horizontal and vertical stresses in the formation, which determines the principal stress magnitudes and therefore the closure stress and fracture orientation; the elastic moduli from sonic logs (converted from dynamic to static values) are the inputs to fracture geometry models, and the accuracy of the fracture design depends directly on the accuracy of these converted moduli.
• Elastic anisotropy in layered formations creates directional dependence of seismic velocities and mechanical properties — most sedimentary rock is not isotropic (having the same properties in all directions) but rather transversely isotropic (TI), with different elastic properties in the vertical direction (perpendicular to bedding) than in the horizontal direction (parallel to bedding); shale is the most common example, with horizontal P-wave velocity typically 10-30% higher than vertical P-wave velocity due to preferred alignment of clay mineral platelets parallel to bedding; this seismic anisotropy affects the accuracy of seismic depth conversion (depth errors if isotropic velocities are assumed in anisotropic shale overburden) and the mechanical anisotropy affects wellbore stability (the failure stress depends on orientation relative to bedding) and hydraulic fracture height containment (fractures may be deflected or arrested at bedding-parallel planes of weakness); accounting for TI anisotropy in seismic processing and geomechanical modeling is increasingly standard in fields with significant shale sections, but requires cross-dipole acoustic data and careful integration of borehole and seismic measurements.